ATLAS Offline Software
CaloSwEtamod_v4.py
Go to the documentation of this file.
1 # Copyright (C) 2002-2020 CERN for the benefit of the ATLAS collaboration
2 
3 #
4 # File: CaloClusterCorrection/python/CaloSwEtamod_v4.py
5 # Created: Nov 2006, sss
6 # Purpose: Eta modulation corrections, v4.
7 #
8 # These corrections were derived using 12.0.3 simulation and reconstruction,
9 # for cluster sizes of 5x5, 3x5, and 3x7, for both electrons
10 # and photons.
11 #
12 # This was added in CaloClusterCorrection-00-02-64, in 12.0.4.
13 #
14 
15 
16 from CaloClusterCorrection.constants import CALOCORR_COMBINED2
17 from AthenaCommon.SystemOfUnits import GeV
18 
19 
20 # Energies at which corrections were tabulated.
21 etamod_ele_sample_energies = [10*GeV, 25*GeV, 50*GeV, 75*GeV,
22  100*GeV, 200*GeV, 500*GeV, 1000*GeV]
23 etamod_gam_sample_energies = [5*GeV, 10*GeV, 25*GeV, 50*GeV, 75*GeV,
24  100*GeV, 200*GeV, 500*GeV, 1000*GeV]
25 
26 
27 etamod_v4_ele5_55 = [0.983806, 3.750856, -154.635676]
28 etamod_v4_ele5_35 = [0.999497, 0.188522, -8.791095]
29 etamod_v4_ele5_37 = [1.001259, -0.474412, 9.681194]
30 etamod_v4_ele10_55 = [0.993712, 1.209621, -49.526417]
31 etamod_v4_ele10_35 = [0.991334, 1.695083, -60.905850]
32 etamod_v4_ele10_37 = [0.989896, 2.063848, -73.906376]
33 etamod_v4_ele25_55 = [0.995942, 0.555437, -15.500041]
34 etamod_v4_ele25_35 = [0.992321, 1.538730, -52.050354]
35 etamod_v4_ele25_37 = [0.994053, 1.392923, -54.094068]
36 etamod_v4_ele50_55 = [0.996982, 0.528204, -17.990771]
37 etamod_v4_ele50_35 = [0.993117, 1.512039, -56.901355]
38 etamod_v4_ele50_37 = [0.993183, 1.473668, -54.724100]
39 etamod_v4_ele75_55 = [0.997754, 0.377698, -16.414560]
40 etamod_v4_ele75_35 = [0.993308, 1.318025, -46.984286]
41 etamod_v4_ele75_37 = [0.992406, 1.573704, -57.762698]
42 etamod_v4_ele100_55 = [0.998076, 0.308318, -9.139902]
43 etamod_v4_ele100_35 = [0.992809, 1.425346, -51.864830]
44 etamod_v4_ele100_37 = [0.992503, 1.586192, -58.698449]
45 etamod_v4_ele200_55 = [0.997283, 0.478227, -16.154217]
46 etamod_v4_ele200_35 = [0.993058, 1.496153, -56.248150]
47 etamod_v4_ele200_37 = [0.992762, 1.558639, -58.741619]
48 etamod_v4_ele500_55 = [0.997707, 0.421512, -14.367443]
49 etamod_v4_ele500_35 = [0.992741, 1.499762, -55.405371]
50 etamod_v4_ele500_37 = [0.992597, 1.570775, -58.663627]
51 etamod_v4_ele1000_55 = [0.997539, 0.412089, -12.911872]
52 etamod_v4_ele1000_35 = [0.992670, 1.557133, -57.989423]
53 etamod_v4_ele1000_37 = [0.992992, 1.484024, -55.313061]
54 
55 etamod_v4_gam5_55 = [0.990666, 1.975982, -93.145708]
56 etamod_v4_gam5_35 = [0.979576, 3.897743, -146.429330]
57 etamod_v4_gam5_37 = [0.990993, 1.612493, -68.488132]
58 etamod_v4_gam10_55 = [0.990695, 1.738726, -61.737552]
59 etamod_v4_gam10_35 = [0.997539, 1.247956, -56.066169]
60 etamod_v4_gam10_37 = [0.991926, 2.244305, -91.768508]
61 etamod_v4_gam25_55 = [0.996626, 0.767360, -29.799398]
62 etamod_v4_gam25_35 = [0.989307, 2.261203, -85.400993]
63 etamod_v4_gam25_37 = [0.993203, 1.603192, -64.949502]
64 etamod_v4_gam50_55 = [0.997944, 0.265441, -5.366243]
65 etamod_v4_gam50_35 = [0.991916, 1.677999, -60.542668]
66 etamod_v4_gam50_37 = [0.991638, 1.689945, -61.928597]
67 etamod_v4_gam75_55 = [0.997735, 0.386117, -13.184488]
68 etamod_v4_gam75_35 = [0.992081, 1.694850, -63.205380]
69 etamod_v4_gam75_37 = [0.992619, 1.604587, -60.797382]
70 etamod_v4_gam100_55 = [0.996571, 0.666099, -24.579692]
71 etamod_v4_gam100_35 = [0.991244, 1.791867, -67.731569]
72 etamod_v4_gam100_37 = [0.991843, 1.740992, -64.891654]
73 etamod_v4_gam200_55 = [0.996892, 0.582237, -19.943738]
74 etamod_v4_gam200_55 = [0.996909, 0.573404, -19.243274]
75 etamod_v4_gam200_35 = [0.992165, 1.597763, -58.265427]
76 etamod_v4_gam200_37 = [0.992049, 1.642700, -60.023122]
77 etamod_v4_gam500_55 = [0.997435, 0.456722, -15.115561]
78 etamod_v4_gam500_35 = [0.992351, 1.594982, -58.915334]
79 etamod_v4_gam500_37 = [0.992352, 1.592001, -59.154386]
80 etamod_v4_gam1000_55 = [0.997663, 0.399569, -12.693380]
81 etamod_v4_gam1000_35 = [0.993040, 1.504923, -56.415771]
82 etamod_v4_gam1000_37 = [0.993451, 1.568644, -58.555394]
83 
84 
85 etamod_v4_ele55_correction = [
86  #etamod_v4_ele5_55,
87  etamod_v4_ele10_55,
88  etamod_v4_ele25_55,
89  etamod_v4_ele50_55,
90  etamod_v4_ele75_55,
91  etamod_v4_ele100_55,
92  etamod_v4_ele200_55,
93  etamod_v4_ele500_55,
94  etamod_v4_ele1000_55,
95  ]
96 
97 
98 etamod_v4_ele35_correction = [
99  #etamod_v4_ele5_35,
100  etamod_v4_ele10_35,
101  etamod_v4_ele25_35,
102  etamod_v4_ele50_35,
103  etamod_v4_ele75_35,
104  etamod_v4_ele100_35,
105  etamod_v4_ele200_35,
106  etamod_v4_ele500_35,
107  etamod_v4_ele1000_35,
108  ]
109 
110 
111 etamod_v4_ele37_correction = [
112  #etamod_v4_ele5_37,
113  etamod_v4_ele10_37,
114  etamod_v4_ele25_37,
115  etamod_v4_ele50_37,
116  etamod_v4_ele75_37,
117  etamod_v4_ele100_37,
118  etamod_v4_ele200_37,
119  etamod_v4_ele500_37,
120  etamod_v4_ele1000_37,
121  ]
122 
123 
124 etamod_v4_gam55_correction = [
125  etamod_v4_gam5_55,
126  etamod_v4_gam10_55,
127  etamod_v4_gam25_55,
128  etamod_v4_gam50_55,
129  etamod_v4_gam75_55,
130  etamod_v4_gam100_55,
131  etamod_v4_gam200_55,
132  etamod_v4_gam500_55,
133  etamod_v4_gam1000_55,
134  ]
135 
136 
137 etamod_v4_gam35_correction = [
138  etamod_v4_gam5_35,
139  etamod_v4_gam10_35,
140  etamod_v4_gam25_35,
141  etamod_v4_gam50_35,
142  etamod_v4_gam75_35,
143  etamod_v4_gam100_35,
144  etamod_v4_gam200_35,
145  etamod_v4_gam500_35,
146  etamod_v4_gam1000_35,
147  ]
148 
149 
150 etamod_v4_gam37_correction = [
151  etamod_v4_gam5_37,
152  etamod_v4_gam10_37,
153  etamod_v4_gam25_37,
154  etamod_v4_gam50_37,
155  etamod_v4_gam75_37,
156  etamod_v4_gam100_37,
157  etamod_v4_gam200_37,
158  etamod_v4_gam500_37,
159  etamod_v4_gam1000_37,
160  ]
161 
162 
163 class CaloSwEtamod_v4_parms:
164  region = CALOCORR_COMBINED2
165  energy_degree = 2
166  rfac = [ [0, 1], [3, 1] ]
167  rfac_degree = 1
168 
169  energies = {'ele55' : etamod_ele_sample_energies,
170  'ele35' : etamod_ele_sample_energies,
171  'ele37' : etamod_ele_sample_energies,
172  'gam55' : etamod_gam_sample_energies,
173  'gam35' : etamod_gam_sample_energies,
174  'gam37' : etamod_gam_sample_energies,
175  }
176 
177  correction = {'ele55' : etamod_v4_ele55_correction,
178  'ele35' : etamod_v4_ele35_correction,
179  'ele37' : etamod_v4_ele37_correction,
180  'gam55' : etamod_v4_gam55_correction,
181  'gam35' : etamod_v4_gam35_correction,
182  'gam37' : etamod_v4_gam37_correction,
183  }
184 
185 
186 # Rescale all the fits so that the integral over a cell is exactly 1.
187 def _etamod_rescale (pars, cellwidth = 0.025):
188  scale = pars[0] + pars[1]*cellwidth/2 + pars[2]*cellwidth**2/3
189  pars[0] = pars[0] /scale
190  pars[1] = pars[1] /scale
191  pars[2] = pars[2] /scale
192  return
193 def _etamod_rescale_all(corr):
194  for v in corr.values():
195  for p in v:
196  _etamod_rescale (p)
197  return
198 _etamod_rescale_all (CaloSwEtamod_v4_parms.correction)
SystemOfUnits
CaloSwEtamod_v4._etamod_rescale
def _etamod_rescale(pars, cellwidth=0.025)
Definition: CaloSwEtamod_v4.py:187
CaloSwEtamod_v4.CaloSwEtamod_v4_parms
Definition: CaloSwEtamod_v4.py:163
CaloSwEtamod_v4._etamod_rescale_all
def _etamod_rescale_all(corr)
Definition: CaloSwEtamod_v4.py:193
Generated on Fri Jun 28 2024 21:12:00 for ATLAS Offline Software by   doxygen 1.8.18