VERITAS-Observatory · GernotMaier · Jan 4, 2026 · Dec 12, 2025 · Dec 12, 2025 · Dec 12, 2025
diff --git a/CITATION.cff b/CITATION.cff
@@ -5,7 +5,7 @@ message: If you use this software, please cite it using these metadata.
 title: "Eventdisplay Analysis Files for VERITAS"
 doi:  10.5281/zenodo.10616288
 version: 4.92.0
-date-released: 2025-12-31
+date-released: 2026-01-31
 keywords:
 - "gamma-ray astronomy"
 - "astronomy software"

diff --git a/DispXGB/.gitignore b/DispXGB/.gitignore
@@ -0,0 +1,3 @@
+V6*
+V4*
+V5*
diff --git a/DispXGB/AP/.gitkeep b/DispXGB/AP/.gitkeep
diff --git a/DispXGB/copy_dispXGB.sh b/DispXGB/copy_dispXGB.sh
@@ -0,0 +1,76 @@
+#!/bin/bash
+# copy dispXGB model files from IRF production site
+# Stereo and gamma/hadron analysis
+#
+# hardwired
+# - SIMTYPE (e.g., CARE_June2020)
+#
+
+IRFVERSION=$(cat ../IRFVERSION)
+ANALYSISTYPE="${VERITAS_ANALYSIS_TYPE:0:2}"
+SIMTYPE="CARE_UV_2212"
+SIMTYPE="GRISU"
+SIMTYPE="CARE_RedHV"
+SIMTYPE="CARE_June2020"
+SIMTYPE="CARE_24_20"
+SIMTYPE="CARE_RedHV_Feb2024"
+SIMTYPE="CARE_202404"
+
+echo "COPY dispXGB  for ${IRFVERSION}, analysis type ${ANALYSISTYPE}, and simulation type ${SIMTYPE}"
+
+for Z in XZE LZE MZE SZE
+do
+    if [[ $Z == "XZE" ]]; then
+        ZE="60deg"
+    elif [[ $Z == "LZE" ]]; then
+        ZE="55deg"
+    elif [[ $Z == "MZE" ]]; then
+        ZE="45deg"
+    elif [[ $Z == "SZE" ]]; then
+        ZE="20deg"
+    fi
+    echo "Zenith bin $Z $ZE"
+    for A in ATM61 ATM62
+    do
+        if [[ ${SIMTYPE} == "GRISU" ]]; then
+            A=${A/6/2}
+            EPOCHS="V4 V5"
+        elif [[ ${SIMTYPE} == *"UV"* ]]; then
+            EPOCHS=$(cat ../IRF_EPOCHS_obsfilter.dat | sort -u)
+            if [[ ${A} == "ATM62" ]]; then
+                continue
+            fi
+        else
+            if [[ ${A} == "ATM62" ]]; then
+                EPOCHS=$(cat ../IRF_EPOCHS_SUMMER.dat | sort -u)
+            else
+                EPOCHS=$(cat ../IRF_EPOCHS_WINTER.dat | sort -u)
+            fi
+            EPOCHS="V6_2016_2017"
+        fi
+        for E in $EPOCHS
+        do
+            echo "EPOCH ${E} ATMO ${A} (**NOTE EPOCHS ARE IGNORED**)"
+            if [[ ${SIMTYPE} == *"RedHV"* ]]; then
+                ODIR="${VERITAS_ANALYSIS_TYPE:0:2}/${E}_${A}_redHV"
+            elif [[ ${SIMTYPE} == *"UV"* ]]; then
+                ODIR="${VERITAS_ANALYSIS_TYPE:0:2}/${E}_${A}_UV"
+            else
+                ODIR="${VERITAS_ANALYSIS_TYPE:0:2}/${E}_${A}"
+            fi
+            IDIR="${VERITAS_IRFPRODUCTION_DIR}/${IRFVERSION}/${ANALYSISTYPE}/${SIMTYPE}"
+
+            # Stereo analysis
+            ZDIR="${ODIR}/${Z}"
+            mkdir -p ${ZDIR}
+            SDIR="${IDIR}/${E}_${A}_gamma/TrainXGBStereoAnalysis/ze${ZE}"
+            cp -v -i ${SDIR}/* ${ZDIR}/
+
+            # Gamma/hadron BDTs (zenith angle independent)
+            if [[ $Z == "SZE" ]]; then
+                GDIR="${IDIR}/${E}_${A}_gamma/TrainXGBGammaHadron"
+                cp -v -i ${GDIR}/*.{joblib,log} ${ODIR}/
+            fi
+        done
+   done
+done
diff --git a/GammaHadronCutFiles/ANASUM.GammaHadron-Cut-NTel2-Extended040-Hard-TMVA-BDT.dat b/GammaHadronCutFiles/ANASUM.GammaHadron-Cut-NTel2-Extended040-Hard-TMVA-BDT.dat
@@ -91,7 +91,6 @@ v) TMVA Cuts
 * TMVAPARAMETER V5 BDT 0 3 0 2 /GammaHadronBDTs/AP/V5/NTel2-Hard/  BDT
 * TMVAPARAMETER V6 BDT 0 8 0 2 /GammaHadronBDTs/AP/V6/NTel2-Hard/  BDT
 
- TMVACUTS <required significance> <required number of events> <observing time> <particle rate file> <fixed signal efficiency> <minimum source strength> <conversion factor>
 
 * TMVA_MVACut V4_ATM21 00 0.142 01 0.142 02 0.142 10 0.139 11 0.133 12 0.156 20 0.12 21 0.138 22 0.134 30 0.0951 31 0.129 32 0.147
 * TMVA_MVACut V4_ATM22 00 0.148 01 0.148 02 0.148 10 0.144 11 0.136 12 0.174 20 0.125 21 0.148 22 0.155 30 0.0983 31 0.129 32 0.16

diff --git a/GammaHadronCutFiles/ANASUM.GammaHadron-Cut-NTel2-Extended040-Moderate-TMVA-BDT.dat b/GammaHadronCutFiles/ANASUM.GammaHadron-Cut-NTel2-Extended040-Moderate-TMVA-BDT.dat
@@ -91,7 +91,6 @@ v) TMVA Cuts
 * TMVAPARAMETER V5 BDT 0 3 0 2 /GammaHadronBDTs/AP/V5/NTel2-Moderate/   BDT
 * TMVAPARAMETER V6 BDT 0 8 0 2 /GammaHadronBDTs/AP/V6/NTel2-Moderate/   BDT
 
- TMVACUTS <required significance> <required number of events> <observing time> <particle rate file> <fixed signal efficiency> <minimum source strength> <conversion factor>
 
 * TMVA_MVACut V4_ATM21 00 0.0719 01 0.0719 02 0.0719 10 0.135 11 0.142 12 0.141 20 0.111 21 0.125 22 0.138 30 0.0628 31 0.113 32 0.124
 * TMVA_MVACut V4_ATM22 00 0.0719 01 0.0719 02 0.0719 10 0.133 11 0.146 12 0.148 20 0.106 21 0.128 22 0.145 30 0.0633 31 0.113 32 0.129

diff --git a/GammaHadronCutFiles/ANASUM.GammaHadron-Cut-NTel2-PointSource-Hard-TMVA-BDT.dat b/GammaHadronCutFiles/ANASUM.GammaHadron-Cut-NTel2-PointSource-Hard-TMVA-BDT.dat
@@ -91,7 +91,6 @@ v) TMVA Cuts
 * TMVAPARAMETER V5 BDT 0 3 0 2 /GammaHadronBDTs/AP/V5/NTel2-Hard/  BDT
 * TMVAPARAMETER V6 BDT 0 8 0 2 /GammaHadronBDTs/AP/V6/NTel2-Hard/  BDT
 
- TMVACUTS <required significance> <required number of events> <observing time> <particle rate file> <fixed signal efficiency> <minimum source strength> <conversion factor>
 
 * TMVA_MVACut V6_2016_2017_ATM61 10 0.102 11 0.0611 20 0.106 21 0.0827 30 0.113 31 0.0975 32 0.078 40 0.126 41 0.108 42 0.0678 50 0.127 51 0.119 52 0.0914 60 0.102 61 0.123 62 0.113 70 0.0938 71 0.121 72 0.121 80 0.0916 81 0.12 82 0.141
 * TMVA_MVACut V6_2024_2024s_ATM62 10 0.104 11 -0.0155 20 0.106 21 0.0847 30 0.108 31 0.0975 32 0.119 40 0.113 41 0.104 42 0.0624 50 0.106 51 0.116 52 0.0855 60 0.0896 61 0.116 62 0.105 70 0.0917 71 0.112 72 0.113 80 0.0875 81 0.111 82 0.136
diff --git a/GammaHadronCutFiles/ANASUM.GammaHadron-Cut-NTel2-PointSource-Moderate-TMVA-BDT.dat b/GammaHadronCutFiles/ANASUM.GammaHadron-Cut-NTel2-PointSource-Moderate-TMVA-BDT.dat
@@ -91,7 +91,6 @@ v) TMVA Cuts
 * TMVAPARAMETER V5 BDT 0 3 0 2 /GammaHadronBDTs/AP/V5/NTel2-Moderate/   BDT
 * TMVAPARAMETER V6 BDT 0 8 0 2 /GammaHadronBDTs/AP/V6/NTel2-Moderate/   BDT
 
- TMVACUTS <required significance> <required number of events> <observing time> <particle rate file> <fixed signal efficiency> <minimum source strength> <conversion factor>
 
 * TMVA_MVACut V6_2016_2017_ATM61 00 0.118 01 0.099 10 0.131 11 0.113 12 0.0393 20 0.142 21 0.129 22 0.0442 30 0.151 31 0.139 32 0.0956 40 0.148 41 0.152 42 0.131 50 0.136 51 0.153 52 0.145 60 0.114 61 0.143 62 0.154 70 0.1 71 0.135 72 0.158 80 0.0823 81 0.126 82 0.165
 * TMVA_MVACut V6_2024_2024s_ATM62 00 0.122 01 0.0971 10 0.127 11 0.123 12 0.13 20 0.14 21 0.133 22 0.0464 30 0.141 31 0.141 32 0.0959 40 0.133 41 0.147 42 0.127 50 0.116 51 0.142 52 0.139 60 0.104 61 0.135 62 0.139 70 0.0961 71 0.128 72 0.144 80 0.0842 81 0.121 82 0.16
diff --git a/GammaHadronCutFiles/ANASUM.GammaHadron-Cut-NTel2-PointSource-Soft-TMVA-BDT.dat b/GammaHadronCutFiles/ANASUM.GammaHadron-Cut-NTel2-PointSource-Soft-TMVA-BDT.dat
@@ -91,7 +91,6 @@ v) TMVA Cuts
 * TMVAPARAMETER V5 BDT 0 3 0 2 /GammaHadronBDTs/AP/V5/NTel2-Soft/   BDT
 * TMVAPARAMETER V6 BDT 0 8 0 2 /GammaHadronBDTs/AP/V6/NTel2-Soft/   BDT
 
- TMVACUTS <required significance> <required number of events> <observing time> <particle rate file> <fixed signal efficiency> <minimum source strength> <conversion factor>
 cut values for different epochs and its corresponding directory where to find .root and weight files
 
 * TMVA_MVACut V6_2016_2017_ATM61 00 0.0963 01 0.0795 02 0.229 10 0.113 11 0.0914 12 0.0223 20 0.125 21 0.11 22 0.0163 30 0.128 31 0.118 32 0.066 40 0.116 41 0.127 42 0.0901 50 0.0992 51 0.127 52 0.102 60 0.0793 61 0.112 62 0.113 70 0.0685 71 0.104 72 0.116 80 0.0498 81 0.0829 82 0.117

diff --git a/GammaHadronCutFiles/ANASUM.GammaHadron-Cut-NTel2-PointSource-SuperSoft-NN-TMVA-BDT.dat b/GammaHadronCutFiles/ANASUM.GammaHadron-Cut-NTel2-PointSource-SuperSoft-NN-TMVA-BDT.dat
@@ -91,7 +91,6 @@ v) TMVA Cuts
 * TMVAPARAMETER V5 BDT 0 3 0 2 /GammaHadronBDTs/NN/V5/NTel2-SuperSoft/   BDT
 * TMVAPARAMETER V6 BDT 0 8 0 2 /GammaHadronBDTs/NN/V6/NTel2-SuperSoft/   BDT
 
- TMVACUTS <required significance> <required number of events> <observing time> <particle rate file> <fixed signal efficiency> <minimum source strength> <conversion factor>
 cut values for different epochs and its corresponding directory where to find .root and weight files
 
 * TMVA_MVACut V4_ATM21 00 0.104 01 0.104 02 0.104 10 0.13 11 0.135 12 0.145 20 0.0936 21 0.0998 22 0.112 30 0.0328 31 0.0584 32 0.0631

diff --git a/GammaHadronCutFiles/ANASUM.GammaHadron-Cut-NTel3-Extended040-Hard-TMVA-BDT.dat b/GammaHadronCutFiles/ANASUM.GammaHadron-Cut-NTel3-Extended040-Hard-TMVA-BDT.dat
@@ -91,8 +91,6 @@ v) TMVA Cuts
 * TMVAPARAMETER V5 BDT 0 3 0 2 /GammaHadronBDTs/AP/V5/NTel3-Hard/  BDT
 * TMVAPARAMETER V6 BDT 0 8 0 2 /GammaHadronBDTs/AP/V6/NTel3-Hard/  BDT
 
- TMVACUTS <required significance> <required number of events> <observing time> <particle rate file> <fixed signal efficiency> <minimum source strength> <conversion factor>
-
 * TMVA_MVACut V4_ATM21 00 0.154 01 0.154 02 0.154 10 0.154 11 0.155 12 0.147 20 0.125 21 0.158 22 0.0956 30 0.124 31 0.159 32 0.138
 * TMVA_MVACut V4_ATM22 00 0.158 01 0.158 02 0.158 10 0.158 11 0.157 12 0.166 20 0.134 21 0.158 22 0.108 30 0.126 31 0.164 32 0.14
 

diff --git a/GammaHadronCutFiles/ANASUM.GammaHadron-Cut-NTel3-PointSource-Hard-TMVA-BDT.dat b/GammaHadronCutFiles/ANASUM.GammaHadron-Cut-NTel3-PointSource-Hard-TMVA-BDT.dat
@@ -91,7 +91,6 @@ v) TMVA Cuts
 * TMVAPARAMETER V5 BDT 0 3 0 2 /GammaHadronBDTs/AP/V5/NTel3-Hard/  BDT
 * TMVAPARAMETER V6 BDT 0 8 0 2 /GammaHadronBDTs/AP/V6/NTel3-Hard/  BDT
 
- TMVACUTS <required significance> <required number of events> <observing time> <particle rate file> <fixed signal efficiency> <minimum source strength> <conversion factor>
 
  * TMVA_MVACut V6_2016_2017_ATM61 10 0.0932 11 0.0455 20 0.0988 21 0.0774 30 0.105 31 0.0931 32 0.093 40 0.116 41 0.104 42 0.0711 50 0.118 51 0.114 52 0.0905 60 0.0952 61 0.116 62 0.109 70 0.0905 71 0.114 72 0.112 80 0.0867 81 0.112 82 0.135
  * TMVA_MVACut V6_2024_2024s_ATM62 10 0.0967 20 0.095 21 0.0757 30 0.103 31 0.0944 32 0.128 40 0.107 41 0.102 42 0.0597 50 0.0996 51 0.111 52 0.0827 60 0.0863 61 0.111 62 0.103 70 0.0855 71 0.108 72 0.109 80 0.0791 81 0.109 82 0.131
diff --git a/LICENSE b/LICENSE
@@ -1,6 +1,6 @@
 BSD 3-Clause License
 
-Copyright (c) 2022-2025 Eventdisplay Developers
+Copyright (c) 2022-2026 Eventdisplay Developers
 All rights reserved.
 
 Redistribution and use in source and binary forms, with or without

diff --git a/ParameterFiles/ANASUM.runparameter b/ParameterFiles/ANASUM.runparameter
@@ -104,3 +104,5 @@ Spectral reconstruction
 Energy spectra binning (logarithmic axis)
 * ENERGYBINSIZE 0.05
 * DEADTIMECALCULATIONMETHOD 0
+
+XGBFILESUFFIX xgb_stereo
-Original file line number
+Diff line change
@@ -0,0 +1,3 @@
+    V6*
+    V4*
+    V5*