ATLAS Offline Software
Loading...
Searching...
No Matches
CaloSwEtaoff_g3.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/CaloSwEtaoff_g3.py
5# Created: Nov 2006, sss
6# Purpose: Eta offset (s-shape) corrections, original G3 version.
7#
8# This is the original G3-based eta offset correction, translated from the
9# original fortran version.
10#
11
12
13from CaloClusterCorrection.constants import \
14 CALOCORR_EMB1, CALOCORR_EMB2, CALOCORR_EME1, CALOCORR_EME2, \
15 EMB1, EME1, EMB2, EME2
16
17
18
22 region = CALOCORR_EMB1
23 correction_coef = 0.001
24 correction_degree = 2
25 interp_degree = 1
26
27 correction_bins = [
28 0.14915, 0.35091, 0.54866, 0.70023, 0.95074, 1.15057
29 ]
30
31 correction = [
32 [0.00060625, 0.46496, 0.43951, 0.44668, 0.46457, 0.44269, 0.45247],
33 [0.00081875, 0.48602, 0.476, 0.4554, 0.4935, 0.47243, 0.4684 ],
34 [0.00103125, 0.47323, 0.43763, 0.45227, 0.45387, 0.4561, 0.42372],
35 [0.00124375, 0.35678, 0.33161, 0.3366, 0.3769, 0.36993, 0.36012],
36 [0.00145625, 0.15322, 0.17447, 0.20115, 0.19289, 0.26117, 0.25023],
37 [0.00166875, -0.1021, -0.0798, -0.00301, 0.02587, 0.08334, 0.10294],
38 [0.00188125, -0.32048, -0.26631, -0.22049, -0.18293, -0.1057, -0.06259],
39 [0.00209375, -0.4314, -0.39734, -0.3599, -0.32175, -0.25388, -0.2258 ],
40 [0.00230625, -0.48008, -0.46247, -0.43731, -0.39469, -0.30764, -0.32463],
41 [0.00251875, -0.44999, -0.45039, -0.39664, -0.40431, -0.32203, -0.34574]
42 ]
43
44
45
49 region = CALOCORR_EME1
50 correction_coef = 0.001
51 correction_degree = 3
52
53 correction = [
54 [ [0.00015, 0.03531],
55 [0.00045, 0.17828],
56 [0.00075, 0.33171],
57 [0.00105, 0.55264],
58 [0.00135, 0.56879],
59 [0.00165, 0.65316],
60 [0.00195, 0.59085],
61 [0.00225, 0.33315],
62 [0.00255, 0.09746],
63 [0.00285, -0.16709],
64 [0.00315, -0.51071],
65 [0.00345, -0.79086],
66 [0.00375, -1.0527],
67 [0.00405, -1.1633],
68 [0.00435, -1.2246],
69 [0.00465, -1.1661],
70 [0.00495, -1.0835],
71 [0.00525, -0.93026],
72 [0.00555, -0.63431],
73 [0.00585, -0.42498] ],
74
75 [ [0.0001, 0.0],
76 [0.0003, 0.25306],
77 [0.0005, 0.34187],
78 [0.0007, 0.50228],
79 [0.0009, 0.57427],
80 [0.0011, 0.55935],
81 [0.0013, 0.52874],
82 [0.0015, 0.39925],
83 [0.0017, 0.20678],
84 [0.0019, -0.04625],
85 [0.0021, -0.33231],
86 [0.0023, -0.62167],
87 [0.0025, -0.78044],
88 [0.0027, -0.85615],
89 [0.0029, -0.9053],
90 [0.0031, -0.84229],
91 [0.0033, -0.77421],
92 [0.0035, -0.66805],
93 [0.0037, -0.53387],
94 [0.0039, 0.0] ],
95
96 [ [7.5e-05, 0.26476],
97 [0.000225, 0.35028],
98 [0.000375, 0.46987],
99 [0.000525, 0.51902],
100 [0.000675, 0.56768],
101 [0.000825, 0.55523],
102 [0.000975, 0.53262],
103 [0.001125, 0.43148],
104 [0.001275, 0.27557],
105 [0.001425, 0.12852],
106 [0.001575, -0.11432],
107 [0.001725, -0.24506],
108 [0.001875, -0.39161],
109 [0.002025, -0.44415],
110 [0.002175, -0.48602],
111 [0.002325, -0.50372],
112 [0.002475, -0.50988],
113 [0.002625, -0.32567],
114 [0.002775, -0.20379],
115 [0.002925, 0.0] ]
116 ]
117
118
119
123 region = CALOCORR_EMB2
124 correction_coef = 0.01
125 residual_coef = 0.001
126 correction_degree = 3
127 residual_eval_degree = 1
128 residual_degree = 2
129
130 correction = [
131 [0.0003125, 0.00637],
132 [0.0009375, 0.04211],
133 [0.0015625, 0.08828],
134 [0.0021875, 0.11382],
135 [0.0028125, 0.16287],
136 [0.0034375, 0.20353],
137 [0.0040625, 0.23561],
138 [0.0046875, 0.27112],
139 [0.0053125, 0.30883],
140 [0.0059375, 0.34944],
141 [0.0065625, 0.36506],
142 [0.0071875, 0.38478],
143 [0.0078125, 0.39651],
144 [0.0084375, 0.39869],
145 [0.0090625, 0.38781],
146 [0.0096875, 0.35214],
147 [0.0103125, 0.30744],
148 [0.0109375, 0.23783],
149 [0.0115625, 0.14325],
150 [0.0121875, 0.03548],
151 [0.0128125, -0.08761],
152 [0.0134375, -0.19445],
153 [0.0140625, -0.28279],
154 [0.0146875, -0.34365],
155 [0.0153125, -0.38791],
156 [0.0159375, -0.41733],
157 [0.0165625, -0.43066],
158 [0.0171875, -0.42128],
159 [0.0178125, -0.41288],
160 [0.0184375, -0.39362],
161 [0.0190625, -0.37207],
162 [0.0196875, -0.34845],
163 [0.0203125, -0.31187],
164 [0.0209375, -0.274],
165 [0.0215625, -0.23781],
166 [0.0221875, -0.195],
167 [0.0228125, -0.15279],
168 [0.0234375, -0.10913],
169 [0.0240625, -0.08176],
170 [0.0246875, -0.04005] ]
171
172 residuals = [
173 [0.00125, 0.27067, 0.12679, 0.0027, -0.09542, -0.06591, -0.34557],
174 [0.00375, 0.27145, 0.16919, 0.07353, 0.00906, -0.1, -0.29619],
175 [0.00625, 0.23725, 0.15188, 0.04369, -0.09484, -0.15978, -0.31441],
176 [0.00875, 0.28479, 0.20336, 0.08277, -0.00911, -0.05585, -0.35899],
177 [0.01125, 0.29884, 0.21377, 0.14735, 0.02146, -0.14185, -0.41464],
178 [0.01375, 0.09564, 0.07504, -0.05496, -0.15623, -0.07018, -0.13584],
179 [0.01625, 0.09158, -0.00519, -0.04834, -0.07372, -0.07021, -0.06172],
180 [0.01875, 0.12776, 0.10608, -0.01719, -0.10394, -0.10419, -0.21907],
181 [0.02125, 0.13389, 0.10165, 0.03954, -0.08366, -0.06337, -0.30239],
182 [0.02375, 0.11144, 0.04308, 0.03373, 0.07283, -0.17986, -0.3195] ]
183
184 residual_bins = [
185 0.14915, 0.35091, 0.54866, 0.70023, 0.95074, 1.15057
186 ]
187
188
189
190
194 region = CALOCORR_EME2
195 correction_degree = 3
196
197 correction = [
198 [0.0003125, 0.00046],
199 [0.0009375, 0.0009],
200 [0.0015625, 0.00153],
201 [0.0021875, 0.00213],
202 [0.0028125, 0.00253],
203 [0.0034375, 0.00294],
204 [0.0040625, 0.00322],
205 [0.0046875, 0.0035],
206 [0.0053125, 0.0037],
207 [0.0059375, 0.0039],
208 [0.0065625, 0.00404],
209 [0.0071875, 0.00414],
210 [0.0078125, 0.00411],
211 [0.0084375, 0.00401],
212 [0.0090625, 0.0038],
213 [0.0096875, 0.00342],
214 [0.0103125, 0.00294],
215 [0.0109375, 0.00219],
216 [0.0115625, 0.00141],
217 [0.0121875, 0.00052],
218 [0.0128125, -0.00038],
219 [0.0134375, -0.00131],
220 [0.0140625, -0.00198],
221 [0.0146875, -0.00265],
222 [0.0153125, -0.00316],
223 [0.0159375, -0.00344],
224 [0.0165625, -0.00366],
225 [0.0171875, -0.00383],
226 [0.0178125, -0.00369],
227 [0.0184375, -0.00363],
228 [0.0190625, -0.00352],
229 [0.0196875, -0.00337],
230 [0.0203125, -0.00311],
231 [0.0209375, -0.00282],
232 [0.0215625, -0.00242],
233 [0.0221875, -0.00212],
234 [0.0228125, -0.0017],
235 [0.0234375, -0.00139],
236 [0.0240625, -0.00095],
237 [0.0246875, 7e-05] ]
238
239
240CaloSwEtaoff_g3_parms = { EMB1 : CaloSwEtaoff_g3_b1_parms,
241 EMB2 : CaloSwEtaoff_g3_b2_parms,
242 EME1 : CaloSwEtaoff_g3_e1_parms,
243 EME2 : CaloSwEtaoff_g3_e2_parms,
244 }
245
246
barrel eta position (S-shape) correction for sampling 1