update second order energy correction given 1M electron sample

LDMX-Software · Dec 5, 2022 · bb668db · bb668db
1 parent e7e24cd
commit bb668db
Showing 1 changed file with 5 additions and 3 deletions.
diff --git a/python/digi.py b/python/digi.py
@@ -198,11 +198,13 @@ def v12(self) :
     def v14(self) :
         """Generated for the v14 geometry
 
-        The secondOrderEnergyCorrection has not been calculated yet so it is set to the
-        suspicious value of 1.
+        The secondOrderEnergyCorrection is deteremined by generating 1M single 4GeV
+        electron events shot directly into the front of the ECal from immediately upstream.
+        The mean of the resulting total recon energy is found by fitting a two-sided normal
+        distribution (one mean, a low and high deviation) to the histogram.
         """
 
-        self.secondOrderEnergyCorrection = 1.;
+        self.secondOrderEnergyCorrection = 4000. / 3940.5;
         self.layerWeights = [ 
                 2.312, 4.312, 6.522, 7.490, 8.595, 10.253, 10.915, 10.915, 10.915, 10.915, 10.915,
                 10.915, 10.915, 10.915, 10.915, 10.915, 10.915, 10.915, 10.915, 10.915, 10.915,