Skip to content

Commit

Permalink
Update Ecal 2nd correction for 8 GeV beam (#1509)
Browse files Browse the repository at this point in the history
  • Loading branch information
tvami authored Jan 13, 2025
1 parent c9b2459 commit 5d7ad33
Showing 1 changed file with 5 additions and 4 deletions.
9 changes: 5 additions & 4 deletions Ecal/python/digi.py
Original file line number Diff line number Diff line change
Expand Up @@ -198,13 +198,14 @@ def v12(self) :
def v14(self) :
"""Generated for the v14 geometry
The secondOrderEnergyCorrection is deteremined by generating 1M single 4GeV
The secondOrderEnergyCorrection is deteremined by generating 1M single 4GeV or 8GeV
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 = 4000. / 3940.5;
#self.secondOrderEnergyCorrection = 4000. / 3940.5;
self.secondOrderEnergyCorrection = 8000. / 7998.3;
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,
Expand All @@ -216,10 +217,10 @@ def reduced_v2(self) :
"""Generated for the reduced v2 geometry
TODO: The secondOrderEnergyCorrection for this geometry has yet to be calculated,
so the v14 secondOrderEnergyCorrection is being used as a placeholder.
so unity is being used as a placeholder.
"""

self.secondOrderEnergyCorrection = 4000. / 3940.5;
self.secondOrderEnergyCorrection = 1.;
self.layerWeights = [
2.312, 5.417, 9.837, 11.910, 11.910, 11.910
]

0 comments on commit 5d7ad33

Please sign in to comment.