Merge pull request #64 from VERITAS-Observatory/v490.7-preprocessing-dev2

v490.7 preprocessing dev2

Author: GernotMaier

.pre-commit-config.yaml

@@ -1,24 +1,24 @@
 repos:
   # https://pycqa.github.io/isort/docs/configuration/black_compatibility.html#integration-with-pre-commit
   - repo: https://github.com/pycqa/isort
-    rev: 5.12.0
+    rev: 5.13.2
     hooks:
       - id: isort
         args: ["--profile", "black", "--filter-files"]
   - repo: https://github.com/psf/black
-    rev: 23.7.0
+    rev: 24.3.0
     hooks:
       - id: black
         args: ["--line-length=100"]
   # https://black.readthedocs.io/en/stable/guides/using_black_with_other_tools.html?highlight=other%20tools#flake8
   - repo: https://github.com/PyCQA/flake8
-    rev: 6.1.0
+    rev: 7.0.0
     hooks:
       - id: flake8
         args: ["--max-line-length=100", "--extend-ignore=E203,E712"]
   # https://github.com/pre-commit/pre-commit-hooks
   - repo: https://github.com/pre-commit/pre-commit-hooks
-    rev: v4.4.0
+    rev: v4.5.0
     hooks:
       - id: trailing-whitespace
       - id: end-of-file-fixer

CITATION.cff

@@ -4,8 +4,8 @@ cff-version: 1.2.0
 message: If you use this software, please cite it using these metadata.
 title: "Eventdisplay Analysis Files for VERITAS"
 doi: 10.5281/zenodo.10616289
-version: 4.9
-date-released: 2024-02-05
+version: 4.9.7
+date-released: 2024-03-22
 keywords:
 - "gamma-ray astronomy"
 - "astronomy software"

GammaHadronCutFiles/ANASUM.GammaHadron-Cut-NTel2-PointSource-Hard-TMVA-BDT.dat

@@ -7,13 +7,13 @@ Hard cuts (minimum 2 tel)
 All lines without a '*' in the first column are ignored!
 --------------------------------------------------------
 i) stereo quality cuts
-successfull core and direction reconstruction (>=0.)
+successful core and direction reconstruction (>=0.)
 * arraychi2 0. 1.e30
 distance of core position to telescopes [m]
 * arraycore -1. 400.
 
 size of the camera [deg], all events with direction
-outside of that intervall are ignored
+outside of that interval are ignored
 * cameraedge -1. 2.00
 
 second largest image size
@@ -73,7 +73,7 @@ ii) direction cut [deg2]
 * arraytheta2 -1. 0.008
 
 iii) gamma hadron cuts:
-b) mean scaled width cut (use cut selector 0 in runlist)
+b) mean scaled width cut (use cut selector 0 in run list)
 * arraymscw -2. 2.
 * arraymscl -2. 5.
 

GammaHadronCutFiles/ANASUM.GammaHadron-Cut-NTel2-PointSource-Hard-TMVA-Preselection.dat

@@ -7,13 +7,13 @@ BDT preselection cuts for TMVA BDT cuts optimization  (hard cuts)
 All lines without a '*' in the first column are ignored!
 --------------------------------------------------------
 i) stereo quality cuts
-successfull core and direction reconstruction (>=0.)
+successful core and direction reconstruction (>=0.)
 * arraychi2 0. 1.e30
 distance of core position to telescopes [m]
 * arraycore -1. 400.
 
 size of the camera [deg], all events with direction
-outside of that intervall are ignored
+outside of that interval are ignored
 * cameraedge -1. 2.00
 
 second largest image size
@@ -73,7 +73,7 @@ ii) direction cut [deg2]
 * arraytheta2 -1. 0.008
 
 iii) gamma hadron cuts:
-b) mean scaled width cut (use cut selector 0 in runlist)
+b) mean scaled width cut (use cut selector 0 in run list)
 * arraymscw -2. 2.
 * arraymscl -2. 5.
 

GammaHadronCutFiles/ANASUM.GammaHadron-Cut-NTel2-PointSource-Hard.dat

@@ -0,0 +1,90 @@
+cut values for quality and gamma hadron cuts
+========================================================
+
+Hard point source cuts (box cuts)
+
+========================================================
+All lines without a '*' in the first column are ignored!
+--------------------------------------------------------
+i) stereo quality cuts
+successful core and direction reconstruction (>=0.)
+* arraychi2 0. 1.e30
+distance of core position to telescopes [m]
+* arraycore -1. 400.
+
+size of the camera [deg], all events with direction
+outside of that  are ignored
+* cameraedge -1. 2.00
+
+second largest image size
+(this depends on the state of the instrument)
+* sizesecondmax 900. 1.e30 V4
+* sizesecondmax 900. 1.e30 V5
+* sizesecondmax 1400. 1.e30 V6
+
+cut on mean emission height
+* emissionheight 6. 1.e10
+
+cut on direction reconstruction quality
+(applies to Disp method only)
+Third parameter: require sucessful reconstruction for both methods
+* arraydispdiff -1. 0.5 1
+
+Telescope multiplicity
+arrayltrig <bit coded telescope combination> <0=off, 1=on> <epoch> <telescope combinations for which this selection is valid>
+* arrayltrig 0 0
+T1
+* arrayltrig 1 0
+T2
+* arrayltrig 2 0
+T1/T2
+* arrayltrig 3 1
+T3
+* arrayltrig 4 0
+T1/T3 V4: unlikely to get a 2-tel gamma-ray event with T1 and T3
+* arrayltrig 5 0 V4 1234
+* arrayltrig 5 1 V5
+* arrayltrig 5 1 V6
+T2/T3
+* arrayltrig 6 1 V4
+T2/T3 and V5/V6: unlikely to get a 2-tel gamma-ray event with T2 and T3
+* arrayltrig 6 0 V5 1234
+* arrayltrig 6 0 V6 1234
+T1/T2/T3
+* arrayltrig 7 1
+T4
+* arrayltrig 8 0
+T1/T4 and V4: telescopes too close to each other
+* arrayltrig 9 0 V4
+T1/T4 and V5/V6: unlikely to get a 2-tel gamma-ray event with T1 and T4
+* arrayltrig 9 0 V5 1234
+* arrayltrig 9 0 V6 1234
+T2/T4 and V4: unlikely to get a 2-tel gamma-ray event with T2 and T4
+* arrayltrig 10 0 V4 1234
+* arrayltrig 10 1 V5
+* arrayltrig 10 1 V6
+T1/T2/T4
+* arrayltrig 11 1
+T3/T4
+* arrayltrig 12 1
+T1/T3/T4
+* arrayltrig 13 1
+T2/T3/T4
+* arrayltrig 14 1
+T1/T2/T3/T4
+* arrayltrig 15 1
+
+ii) direction cut [deg2]
+* arraytheta2 -1. 0.008
+
+iii) gamma hadron cuts:
+b) mean scaled width cut (use cut selector 0 in run list)
+* arraymscw -1.2 0.3
+* arraymscl -1.2 0.5
+
+iv) energy reconstruction cuts
+echi2 is a measure of the difference between the energies
+estimated per individual telescope
+* arrayechi2 0. 99999.
+cut on energy of reconstruction shower (on linear energy scale, in [TeV])
+* arrayerec 0. 1.e10

GammaHadronCutFiles/ANASUM.GammaHadron-Cut-NTel2-PointSource-Moderate-TMVA-BDT.dat

@@ -7,13 +7,13 @@ Moderate cuts
 All lines without a '*' in the first column are ignored!
 --------------------------------------------------------
 i) stereo quality cuts
-successfull core and direction reconstruction (>=0.)
+successful core and direction reconstruction (>=0.)
 * arraychi2 0. 1.e30
 distance of core position to telescopes [m]
 * arraycore -1. 400.
 
 size of the camera [deg], all events with direction
-outside of that intervall are ignored
+outside of that  are ignored
 * cameraedge -1. 2.00
 
 second largest image size
@@ -73,7 +73,7 @@ ii) direction cut [deg2]
 * arraytheta2 -1. 0.008
 
 iii) gamma hadron cuts:
-b) mean scaled width cut (use cut selector 0 in runlist)
+b) mean scaled width cut (use cut selector 0 in run list)
 * arraymscw -2. 2.
 * arraymscl -2. 5.
 

GammaHadronCutFiles/ANASUM.GammaHadron-Cut-NTel2-PointSource-Moderate-TMVA-Preselection.dat

@@ -7,13 +7,13 @@ BDT preselection cuts for TMVA BDT cuts optimization  (moderate cuts)
 All lines without a '*' in the first column are ignored!
 --------------------------------------------------------
 i) stereo quality cuts
-successfull core and direction reconstruction (>=0.)
+successful core and direction reconstruction (>=0.)
 * arraychi2 0. 1.e30
 distance of core position to telescopes [m]
 * arraycore -1. 400.
 
 size of the camera [deg], all events with direction
-outside of that intervall are ignored
+outside of that  are ignored
 * cameraedge -1. 2.00
 
 second largest image size
@@ -73,7 +73,7 @@ ii) direction cut [deg2]
 * arraytheta2 -1. 0.008
 
 iii) gamma hadron cuts:
-b) mean scaled width cut (use cut selector 0 in runlist)
+b) mean scaled width cut (use cut selector 0 in run list)
 * arraymscw -2. 2.
 * arraymscl -2. 5.
 

GammaHadronCutFiles/ANASUM.GammaHadron-Cut-NTel2-PointSource-Moderate.dat

@@ -7,13 +7,13 @@ Moderate point source cuts (box cuts)
 All lines without a '*' in the first column are ignored!
 --------------------------------------------------------
 i) stereo quality cuts
-successfull core and direction reconstruction (>=0.)
+successful core and direction reconstruction (>=0.)
 * arraychi2 0. 1.e30
 distance of core position to telescopes [m]
 * arraycore -1. 400.
 
 size of the camera [deg], all events with direction
-outside of that intervall are ignored
+outside of that  are ignored
 * cameraedge -1. 2.00
 
 second largest image size
@@ -78,7 +78,7 @@ ii) direction cut [deg2]
 * arraytheta2 -1. 0.008
 
 iii) gamma hadron cuts:
-b) mean scaled width cut (use cut selector 0 in runlist)
+b) mean scaled width cut (use cut selector 0 in run list)
 * arraymscw -1.2 0.3
 * arraymscl -1.2 0.5
 

GammaHadronCutFiles/ANASUM.GammaHadron-Cut-NTel2-PointSource-Soft-TMVA-BDT.dat

@@ -7,13 +7,13 @@ Soft cuts
 All lines without a '*' in the first column are ignored!
 --------------------------------------------------------
 i) stereo quality cuts
-successfull core and direction reconstruction (>=0.)
+successful core and direction reconstruction (>=0.)
 * arraychi2 0. 1.e30
 distance of core position to telescopes [m]
 * arraycore -1. 400.
 
 size of the camera [deg], all events with direction
-outside of that intervall are ignored
+outside of that  are ignored
 * cameraedge -1. 2.00
 
 second largest image size
@@ -73,7 +73,7 @@ ii) direction cut [deg2]
 * arraytheta2 -1. 0.008
 
 iii) gamma hadron cuts:
-b) mean scaled width cut (use cut selector 0 in runlist)
+b) mean scaled width cut (use cut selector 0 in run list)
 * arraymscw -2. 2.
 * arraymscl -2. 5.
 

GammaHadronCutFiles/ANASUM.GammaHadron-Cut-NTel2-PointSource-Soft-TMVA-Preselection.dat

@@ -7,13 +7,13 @@ BDT preselection cuts for TMVA BDT cuts optimization  (soft cuts)
 All lines without a '*' in the first column are ignored!
 --------------------------------------------------------
 i) stereo quality cuts
-successfull core and direction reconstruction (>=0.)
+successful core and direction reconstruction (>=0.)
 * arraychi2 0. 1.e30
 distance of core position to telescopes [m]
 * arraycore -1. 400.
 
 size of the camera [deg], all events with direction
-outside of that intervall are ignored
+outside of that  are ignored
 * cameraedge -1. 2.00
 
 second largest image size
@@ -73,7 +73,7 @@ ii) direction cut [deg2]
 * arraytheta2 -1. 0.008
 
 iii) gamma hadron cuts:
-b) mean scaled width cut (use cut selector 0 in runlist)
+b) mean scaled width cut (use cut selector 0 in run list)
 * arraymscw -2. 2.
 * arraymscl -2. 5.