ATLAS Offline Software
Loading...
Searching...
No Matches
MuidCaloEnergyParam.cxx
Go to the documentation of this file.
1/*
2 Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
3*/
4
6// MuidCaloEnergyParam
7// AlgTool returning the parametrized value for the calorimeter energy
8// deposited by a traversing muon. The energy is parametrized according
9// to trackParameters given at the calo mid-surface.
10//
12
13//<<<<<< INCLUDES >>>>>>
14
15#include "MuidCaloEnergyParam.h"
16
17#include <cmath>
18
19#include "AthenaKernel/Units.h"
20#include "muonEvent/CaloEnergy.h"
21
22namespace Units = Athena::Units;
23namespace {
24 constexpr double FifteenGeV = 15 * Units::GeV;
25 constexpr double PerGeV = 1. / Units::GeV;
26} // namespace
27namespace Rec {
28
29 //<<<<<< CLASS STRUCTURE INITIALIZATION >>>>>>
30
31 MuidCaloEnergyParam::MuidCaloEnergyParam(const std::string& type, const std::string& name, const IInterface* parent) :
32 AthAlgTool(type, name, parent),
33 m_cosmics(false),
34 m_smoothingFraction(0.5),
35 m_binWidth(0.),
36 m_etaOffset(0.),
37 m_inverseWidth(0.),
38 m_etaGranularity(0.) {
39 declareInterface<IMuidCaloEnergyParam>(this);
40 declareProperty("Cosmics", m_cosmics);
41 declareProperty("SmoothingFraction", m_smoothingFraction);
42 }
43
44 MuidCaloEnergyParam::~MuidCaloEnergyParam() = default;
45
46 //<<<<<< PUBLIC MEMBER FUNCTION DEFINITIONS >>>>>>
47
48 StatusCode MuidCaloEnergyParam::initialize() {
49 ATH_MSG_DEBUG(" Initializing ");
50
51 // set up eta bin width and precomputed smoothing constants
52 m_binWidth = 0.1;
53 if (m_smoothingFraction < 0.) m_smoothingFraction = 0.;
54 if (m_smoothingFraction > 1.) m_smoothingFraction = 1.;
55 m_etaOffset = 0.5 * m_binWidth * m_smoothingFraction;
56 m_inverseWidth = 1. / m_binWidth;
57
58 // fix to constant term for ATLAS-GEO-10-00-00 --> derived using a pt=5GeV
59 // single muon sample comparing the actual eloss and the expected according
60 // to the parametrisation after comparison with data, this is only to be
61 // applied to barrel region
62 constexpr double fix[26] = {0.21892, 0.13412, 0.064942, 0.154388, 0.16903, 0.212424, 0.246268, 0.234471, 0.23724,
63 0.209337, 0.186304, 0.100489, 0.0740633, 0.0319171, 0.0352365, 0.0218442, 0.19691, 0.124486,
64 0.0882254, 0.132367, 0.181742, 0.213001, 0.192657, 0.219348, 0.261353, 0.142727};
65
66 // Initialization of the mean energy deposition arrays
67 //(function [0]+[1]*Log(0.0067*p*p)+[2]*p with GeV units)
68 constexpr double meanEnergyLoss_p0[26] = {2.7909, 2.82994, 2.91726, 3.00969, 3.09838, 3.25242, 3.32701, 2.99668, 3.05648,
69 3.50282, 4.41603, 4.64789, 3.90948, 3.29712, 2.86013, 2.8295, 3.13659, 3.6972,
70 3.69033, 3.76311, 3.62784, 3.60054, 3.67615, 3.57151, 3.59287, 3.40605};
71 constexpr double meanEnergyLoss_p1[26] = {0.0497463, 0.0440959, 0.0331098, 0.0565649, 0.0721456, 0.0948259, 0.0825271,
72 0.0543386, 0.0790654, 0.0494549, 0.118818, 0.139855, 0.00587402, 0.0589792,
73 0.0461476, 0.0285649, 0.0460483, 0.130774, 0.0999097, 0.126468, 0.0302406,
74 0.0645617, 0.0416555, 0.0503394, -0.00448376, 0.135223};
75 constexpr double meanEnergyLoss_p2[26] = {0.0102245, 0.0113924, 0.0090614, 0.0105087, 0.0104329, 0.010759, 0.0103308,
76 0.011115, 0.00944668, 0.0126312, 0.0127012, 0.0162508, 0.0131209, 0.00842482,
77 0.00889747, 0.00755592, 0.00944529, 0.0106429, 0.0111037, 0.0113795, 0.0129818,
78 0.0128469, 0.0138801, 0.015034, 0.0135201, 0.0109176};
79 // tabulated coefficients for left assymetric error: pol1 with GeV units
80 constexpr double meanEnergyLoss_left_p0[26] = {0.290254, 0.253085, 0.252006, 0.279803, 0.284166, 0.269267, 0.370476,
81 0.309723, 0.339191, 0.392122, 0.504243, 0.452033, 0.389579, 0.325196,
82 0.241338, 0.230355, 0.319177, 0.341208, 0.262456, 0.300411, 0.255884,
83 0.257605, 0.272939, 0.247617, 0.267175, 0.377071};
84 constexpr double meanEnergyLoss_left_p1[26] = {0.00570519, 0.00594466, 0.00519917, 0.0054666, 0.00552764, 0.00631414, 0.00539921,
85 0.00573473, 0.00508457, 0.00624818, 0.00691296, 0.00864038, 0.00640888, 0.00488251,
86 0.00507001, 0.00462523, 0.00529166, 0.00594752, 0.00640289, 0.00655527, 0.00631339,
87 0.00643479, 0.00749501, 0.00766192, 0.00637578, 0.00567618};
88 // tabulated coefficients for right assymetric error: pol1 with GeV units
89 constexpr double meanEnergyLoss_right_p0[26] = {0.478493, 0.587466, 0.678961, 0.723157, 0.536312, 0.60896, 0.570827,
90 0.641256, 0.713216, 0.648768, 0.684193, 0.472149, 0.697594, 0.648713,
91 0.577021, 0.52603, 0.536587, 0.573851, 0.603372, 0.817339, 0.631103,
92 0.491342, 0.546488, 0.46241, 0.46241, 0.951949};
93 constexpr double meanEnergyLoss_right_p1[26] = {0.0250521, 0.0284547, 0.0202018, 0.0202702, 0.0257285, 0.0230366, 0.0240065,
94 0.0209904, 0.0220515, 0.0296286, 0.0314476, 0.0452839, 0.0343686, 0.0197284,
95 0.0217681, 0.018329, 0.0193914, 0.0305364, 0.0258112, 0.0272972, 0.0279697,
96 0.0228007, 0.0268285, 0.0312118, 0.0312118, 0.0205223};
97
98 for (int i = 0; i < 26; ++i) {
99 if (i < 13) {
100 m_meanEnergyLossP0[i] = meanEnergyLoss_p0[i] + fix[i];
101 } else {
102 m_meanEnergyLossP0[i] = meanEnergyLoss_p0[i];
103 }
104 m_meanEnergyLossP1[i] = meanEnergyLoss_p1[i];
105 m_meanEnergyLossP2[i] = meanEnergyLoss_p2[i];
106 m_meanEnergyLossErrorLeftP0[i] = meanEnergyLoss_left_p0[i];
107 m_meanEnergyLossErrorLeftP1[i] = meanEnergyLoss_left_p1[i];
108 m_meanEnergyLossErrorRightP0[i] = meanEnergyLoss_right_p0[i];
109 m_meanEnergyLossErrorRightP1[i] = meanEnergyLoss_right_p1[i];
110 }
111
112 // Initialization of the mop energy deposition arrays
113 // tabulated coefficients for MOP energy loss and for its error
114 // (function [0]+[1]*Log(0.0067*p*p)+[2]*p with GeV units)
115 constexpr double mopEnergyLoss_p0[26] = {2.57291, 2.57581, 2.65592, 2.73852, 2.81242, 2.94167, 3.01323, 2.72711, 2.75402,
116 3.17494, 4.07206, 4.28525, 3.58247, 2.99343, 2.59807, 2.56842, 2.83854, 3.4031,
117 3.40046, 3.3927, 3.38412, 3.38123, 3.4004, 3.37481, 3.3171, 3.28507};
118 constexpr double mopEnergyLoss_p1[26] = {0.0464874, 0.0213596, 0.0113974, 0.0314924, 0.0283828, 0.0363767, 0.0405601,
119 0.0211039, 0.0273863, 0.0351794, 0.0381777, 0.079885, 0.0176414, 0.0280359,
120 0.0173828, 0.0280728, 0.0366678, 0.0416954, 0.029448, 0.0300134, 0.0262876,
121 0.0229897, 0.0386583, 0.0305293, 0.0214848, 0.0316924};
122 constexpr double mopEnergyLoss_p2[26] = {0.00172448, 0.00225414, 0.00225085, 0.00196091, 0.00229175, 0.00232876, 0.00248703,
123 0.00239371, 0.00204022, 0.00238259, 0.00387238, 0.00398731, 0.0030019, 0.00220271,
124 0.00173978, 0.00153614, 0.00224971, 0.00291137, 0.00297603, 0.00313257, 0.00330651,
125 0.00319543, 0.00287297, 0.00333396, 0.00301831, 0.0027259};
126
127 constexpr double mopEnergyLoss_sigma_p0[26] = {0.129825, 0.120123, 0.110096, 0.129288, 0.13066, 0.128621, 0.139855,
128 0.136144, 0.143521, 0.192355, 0.22055, 0.24224, 0.195584, 0.147344,
129 0.121476, 0.112492, 0.175842, 0.127176, 0.12163, 0.123569, 0.116129,
130 0.118597, 0.126694, 0.129753, 0.128578, 0.109232};
131 constexpr double mopEnergyLoss_sigma_p1[26] = {0.00117338, 0.00114104, 0.00116575, 0.00106617, 0.00121612, 0.0013242, 0.00133937,
132 0.00125934, 0.00114077, 0.00122792, 0.00190011, 0.00203161, 0.00127366, 0.00110695,
133 0.000916009, 0.000920279, 0.00131875, 0.00147309, 0.0015842, 0.0016292, 0.00166247,
134 0.00151376, 0.00149124, 0.00163482, 0.00156624, 0.00148398};
135
136 for (int i = 0; i < 26; ++i) {
137 if (i < 13) {
138 m_mopEnergyLossP0[i] = mopEnergyLoss_p0[i] + fix[i];
139 } else {
140 m_mopEnergyLossP0[i] = mopEnergyLoss_p0[i];
141 }
142 m_mopEnergyLossP1[i] = mopEnergyLoss_p1[i];
143 m_mopEnergyLossP2[i] = mopEnergyLoss_p2[i];
144 m_mopEnergyLossErrorP0[i] = mopEnergyLoss_sigma_p0[i];
145 m_mopEnergyLossErrorP1[i] = mopEnergyLoss_sigma_p1[i];
146 }
147
148 // not sure the mop symmetric is still needed...
149 constexpr double mopSymmetricEnergyLoss_p0[26] = {2.70345, 2.68656, 2.75213, 2.85555, 2.95905, 3.09844, 3.19676, 2.87172, 2.92657,
150 3.38863, 4.28484, 4.54497, 3.77225, 3.18316, 2.72541, 2.68793, 3.08594, 3.52879,
151 3.5142, 3.48087, 3.50184, 3.4893, 3.52484, 3.48467, 3.44116, 3.43226};
152 constexpr double mopSymmetricEnergyLoss_p1[26] = {
153 0.0388166, 0.00777747, 0.0261084, 0.0406787, 0.0322767, 0.038223, 0.0409763, 0.00096891, 0.0439047,
154 0.0430484, 0.0419927, 0.0766375, -0.000214255, 0.0420977, -0.0143302, -0.0022226, 0.0574891, 0.0260815,
155 -0.00223455, 0.0016692, -0.00404322, -0.00620859, 0.00603695, -0.0271089, 0.0386722, 0.0273018};
156 constexpr double mopSymmetricEnergyLoss_p2[26] = {
157 0.00263823, 0.00360428, 0.00293494, 0.00235128, 0.00245869, 0.00304957, 0.00317534, 0.00428433, 0.00231015,
158 0.0029778, 0.00518089, 0.00522931, 0.00446165, 0.0023943, 0.00365962, 0.00310547, 0.00236659, 0.00405503,
159 0.00488483, 0.00499606, 0.0050601, 0.0045985, 0.00515479, 0.00607963, 0.00405829, 0.00440702};
160 for (int i = 0; i < 26; ++i) {
161 m_mopSymmetricEnergyLossP0[i] = mopSymmetricEnergyLoss_p0[i];
162 m_mopSymmetricEnergyLossP1[i] = mopSymmetricEnergyLoss_p1[i];
163 m_mopSymmetricEnergyLossP2[i] = mopSymmetricEnergyLoss_p2[i];
164 }
165
166 // Initialize the em energy deposit fraction
167 constexpr double emMopFraction[26] = {0.100105, 0.101756, 0.1033, 0.103466, 0.10543, 0.103732, 0.106252, 0.120923, 0.116363,
168 0.116984, 0.118591, 0.104436, 0.100136, 0.0934386, 0.104249, 0.122974, 0.122026, 0.113485,
169 0.0978895, 0.0967849, 0.100603, 0.101762, 0.104887, 0.106119, 0.104617, 0.100805};
170 for (int i = 0; i < 26; ++i) { m_emMopFraction[i] = emMopFraction[i]; }
171
172 // Initialization of the dead material map arrays
173 constexpr double dead_eta[277] = {
174 0.005, 0.015, 0.025, 0.035, 0.045, 0.055, 0.065, 0.075, 0.085, 0.095, 0.105, 0.115, 0.125, 0.135, 0.145, 0.155, 0.165, 0.175,
175 0.185, 0.195, 0.205, 0.215, 0.225, 0.235, 0.245, 0.255, 0.265, 0.275, 0.285, 0.295, 0.305, 0.315, 0.325, 0.335, 0.345, 0.355,
176 0.365, 0.375, 0.385, 0.395, 0.405, 0.415, 0.425, 0.435, 0.445, 0.455, 0.465, 0.475, 0.485, 0.495, 0.505, 0.515, 0.525, 0.535,
177 0.545, 0.555, 0.565, 0.575, 0.585, 0.595, 0.605, 0.615, 0.625, 0.635, 0.645, 0.655, 0.665, 0.675, 0.685, 0.695, 0.705, 0.715,
178 0.725, 0.735, 0.745, 0.755, 0.765, 0.775, 0.785, 0.795, 0.805, 0.815, 0.825, 0.835, 0.845, 0.855, 0.865, 0.875, 0.885, 0.895,
179 0.905, 0.915, 0.925, 0.935, 0.945, 0.955, 0.965, 0.975, 0.985, 0.995, 1.005, 1.015, 1.025, 1.035, 1.045, 1.055, 1.065, 1.075,
180 1.085, 1.095, 1.105, 1.115, 1.125, 1.135, 1.145, 1.155, 1.165, 1.175, 1.185, 1.195, 1.205, 1.215, 1.225, 1.235, 1.245, 1.255,
181 1.265, 1.275, 1.285, 1.295, 1.305, 1.315, 1.325, 1.335, 1.345, 1.355, 1.365, 1.375, 1.385, 1.395, 1.405, 1.415, 1.425, 1.435,
182 1.445, 1.455, 1.465, 1.475, 1.485, 1.495, 1.505, 1.515, 1.525, 1.535, 1.545, 1.555, 1.565, 1.575, 1.585, 1.595, 1.605, 1.615,
183 1.625, 1.635, 1.645, 1.655, 1.665, 1.675, 1.685, 1.695, 1.705, 1.715, 1.725, 1.735, 1.745, 1.755, 1.765, 1.775, 1.785, 1.795,
184 1.805, 1.815, 1.825, 1.835, 1.845, 1.855, 1.865, 1.875, 1.885, 1.895, 1.905, 1.915, 1.925, 1.935, 1.945, 1.955, 1.965, 1.975,
185 1.985, 1.995, 2.005, 2.015, 2.025, 2.035, 2.045, 2.055, 2.065, 2.075, 2.085, 2.095, 2.105, 2.115, 2.125, 2.135, 2.145, 2.155,
186 2.165, 2.175, 2.185, 2.195, 2.205, 2.215, 2.225, 2.235, 2.245, 2.255, 2.265, 2.275, 2.285, 2.295, 2.305, 2.315, 2.325, 2.335,
187 2.345, 2.355, 2.365, 2.375, 2.385, 2.395, 2.405, 2.415, 2.425, 2.435, 2.445, 2.455, 2.465, 2.475, 2.485, 2.495, 2.505, 2.515,
188 2.525, 2.535, 2.545, 2.555, 2.565, 2.575, 2.585, 2.595, 2.605, 2.615, 2.625, 2.635, 2.645, 2.655, 2.665, 2.675, 2.685, 2.695,
189 2.705, 2.715, 2.725, 2.735, 2.745, 2.755, 2.765};
190 constexpr double dead_x01[277] = {
191 1.8469, 1.7213, 1.7217, 1.7222, 1.7366, 1.7402, 1.7667, 1.7847, 1.7927, 1.7968, 1.8054, 1.7959, 1.7737, 1.7626, 1.7475, 1.7259,
192 1.7211, 1.7032, 1.6776, 1.6731, 1.6688, 1.6516, 1.6474, 1.6497, 1.6394, 1.6384, 1.6322, 1.6288, 1.6187, 1.6251, 1.6474, 1.6547,
193 1.6712, 1.6758, 1.6901, 1.6935, 1.6883, 1.6966, 1.7005, 1.6978, 1.7044, 1.7184, 1.7257, 1.7322, 1.7327, 1.7402, 1.7483, 1.7332,
194 1.7477, 1.7599, 1.783, 1.7907, 1.8299, 1.8343, 1.8427, 1.8513, 1.8271, 1.8276, 1.8345, 1.8702, 1.8821, 1.8889, 1.9048, 1.9538,
195 2.0558, 2.1443, 2.2201, 2.2642, 2.2656, 2.2472, 2.2776, 2.2798, 2.3388, 2.3809, 2.3763, 2.3869, 2.4254, 2.4758, 2.5282, 2.5364,
196 2.5637, 2.5917, 2.5999, 2.6376, 2.6386, 2.6423, 2.6479, 2.674, 2.676, 2.6806, 2.7002, 2.6985, 2.704, 2.6898, 2.7206, 2.7414,
197 2.7647, 2.7856, 2.7867, 2.8078, 2.8187, 2.8461, 2.8572, 2.8544, 2.8487, 2.8796, 2.9303, 2.9401, 2.928, 2.9376, 2.9437, 3.0083,
198 2.9823, 2.9945, 2.9809, 3.0485, 3.0494, 3.0595, 3.0642, 3.1098, 3.1232, 3.1306, 3.0995, 3.1497, 3.1399, 3.1285, 3.1489, 3.1748,
199 3.2728, 3.2219, 3.2252, 3.2826, 3.3505, 3.4002, 3.4492, 3.4564, 3.4622, 3.4631, 3.4628, 3.5027, 3.5506, 3.6616, 3.8859, 4.0115,
200 4.1715, 4.3591, 4.4486, 4.4909, 4.5445, 4.5871, 4.6396, 4.398, 3.6563, 4.672, 3.9715, 3.0968, 2.7105, 2.7465, 2.7719, 2.7361,
201 2.8049, 2.7947, 2.7961, 2.8328, 2.8847, 2.8935, 3.1969, 3.5616, 3.9473, 3.3614, 2.5151, 2.3433, 2.3557, 2.0129, 2.0082, 2.0186,
202 1.9549, 1.9559, 1.873, 1.9516, 1.9924, 2.0921, 2.1637, 2.0858, 2.0419, 1.9849, 2.0455, 1.9273, 1.8368, 1.7655, 1.779, 1.726,
203 1.724, 1.7487, 1.7251, 1.6957, 1.7127, 1.6828, 1.6412, 1.6795, 1.6421, 1.631, 1.6366, 1.6835, 1.6884, 1.6811, 1.6724, 1.6534,
204 1.6404, 1.6498, 1.6268, 1.6205, 1.6453, 1.657, 1.6372, 1.6172, 1.5865, 1.561, 1.5822, 1.6654, 1.6161, 1.57, 1.5173, 1.4629,
205 1.5068, 1.4886, 1.496, 1.5121, 1.4521, 1.4178, 1.4192, 1.4095, 1.42, 1.4151, 1.4129, 1.3897, 1.4006, 1.3767, 1.3668, 1.352,
206 1.3715, 1.3606, 1.3526, 1.3903, 1.4008, 1.3913, 1.363, 1.3655, 1.3634, 1.366, 1.3687, 1.3626, 1.3925, 1.3764, 1.3442, 1.3366,
207 1.3301, 1.3277, 1.3036, 1.304, 1.3068, 1.3114, 1.303, 1.345, 1.3567, 1.3424, 1.3092, 1.3122, 1.3137, 1.3231, 1.3236, 1.2951,
208 1.2928, 1.3456, 1.8221, 2.0029, 2.0075};
209 constexpr double dead_x02[277] = {
210 27.3438, 27.2208, 27.2263, 27.243, 31.8853, 28.1255, 27.5756, 27.0132, 27.0413, 27.0637, 27.1041, 27.1223, 27.1304, 27.1563,
211 27.1702, 27.1865, 27.2222, 27.8551, 28.488, 30.206, 28.5837, 28.0049, 27.4459, 27.5156, 27.5572, 27.6216, 27.6856, 27.7451,
212 27.8147, 27.8949, 27.9914, 28.078, 28.1787, 28.2678, 28.3671, 28.4603, 29.1869, 29.9369, 31.3911, 29.6731, 29.6085, 29.0849,
213 29.2069, 29.3229, 29.4405, 29.5618, 29.6891, 29.7915, 29.9375, 30.0806, 30.2322, 30.3764, 30.5561, 30.6942, 30.8476, 31.6978,
214 31.8183, 33.6375, 33.149, 32.3431, 32.5023, 31.9797, 32.1619, 32.3803, 32.6547, 32.9188, 33.1731, 33.3989, 33.5851, 33.7547,
215 33.9763, 34.1728, 34.4293, 34.676, 35.3402, 35.8687, 36.8943, 37.3228, 37.0672, 35.8788, 30.0848, 30.2997, 30.4957, 30.7282,
216 30.917, 31.1227, 31.3252, 31.5502, 31.7627, 31.9752, 32.2057, 32.4178, 32.6401, 32.8458, 33.9831, 34.4918, 34.9247, 35.4172,
217 35.0033, 35.2645, 34.6033, 34.8723, 35.1359, 35.8106, 37.0775, 37.4231, 37.5118, 37.6158, 38.0132, 38.5355, 38.9488, 40.1903,
218 40.5554, 39.3216, 40.47, 40.5203, 40.3146, 40.0119, 40.0371, 40.0744, 40.0102, 40.0563, 39.9966, 40.0213, 40.0384, 39.9564,
219 38.3485, 37.0049, 36.011, 34.3675, 32.746, 31.6373, 30.2222, 28.8545, 31.2748, 30.5887, 29.1005, 28.1896, 28.9692, 29.0944,
220 29.5132, 29.8182, 30.2405, 30.4836, 30.6649, 31.4598, 31.8968, 30.9419, 32.2847, 32.2785, 29.4552, 28.8184, 27.8418, 26.9469,
221 26.5213, 26.0314, 26.0441, 26.4897, 26.9313, 27.3156, 27.7043, 27.9047, 28.1472, 28.4043, 28.6613, 28.8634, 29.3494, 29.8874,
222 30.4377, 30.0112, 29.3163, 29.2913, 29.4448, 29.2385, 29.3655, 29.5045, 29.5655, 29.6873, 29.7209, 29.9143, 30.0671, 30.2626,
223 30.4579, 30.4513, 30.5094, 30.5531, 30.7121, 30.6904, 30.6949, 30.7176, 30.8232, 30.9326, 31.0953, 31.376, 31.6177, 31.8488,
224 32.1461, 32.4569, 32.7735, 33.1196, 33.4516, 33.7533, 34.0885, 34.4624, 34.7467, 34.9918, 35.2164, 35.4157, 35.6083, 35.812,
225 35.9719, 36.1537, 36.3483, 36.5131, 36.65, 36.7811, 36.8951, 37.0119, 37.1696, 37.3852, 37.4656, 37.5443, 37.6143, 37.6783,
226 37.839, 37.9339, 38.0519, 38.1968, 38.2431, 38.3125, 38.4156, 38.5057, 38.6145, 38.7049, 38.797, 38.8656, 38.9668, 39.031,
227 39.108, 39.177, 39.2796, 39.2613, 39.3312, 39.4452, 39.6109, 39.6561, 39.7265, 39.7582, 39.7801, 39.8325, 40.0373, 30.445,
228 30.4259, 30.3236, 30.335, 30.4099, 30.4834, 30.4998, 30.549, 30.6477, 30.6725, 30.7221, 30.7792, 30.8878, 30.9677, 31.025,
229 31.0644, 31.1425, 31.2212, 31.3097, 31.392, 31.4477, 31.5309, 31.6727, 32.2409, 32.5151, 32.6164};
230 constexpr double dead_x0tot[277] = {
231 108.064, 114.015, 102.328, 103.979, 108.575, 99.904, 107.609, 114.499, 110.58, 105.503, 111.6, 109.618, 107.369, 108.663,
232 108.656, 108.724, 107.947, 111.123, 113.801, 114.927, 110.352, 113.048, 109.105, 112.754, 110.352, 112.808, 113.401, 111.948,
233 113.811, 111.447, 114.069, 115.163, 115.728, 113.215, 116.055, 115.433, 116.233, 118.157, 119.787, 118.035, 117.5, 117.333,
234 116.777, 119.61, 119.342, 121.258, 119.331, 121.887, 122.346, 122.314, 122.692, 122.296, 123.452, 124.217, 125.029, 126.766,
235 127.559, 129.726, 129.524, 130.828, 128.679, 129.477, 129.923, 131.236, 132.793, 132.212, 132.597, 133.999, 131.29, 129.22,
236 126.212, 126.427, 127.45, 124.279, 119.29, 115.547, 114.276, 116.323, 117.272, 117.058, 112.01, 112.727, 114.061, 115.384,
237 115.753, 114.081, 112.184, 111.034, 109.526, 109.854, 111.763, 113.335, 114.989, 116.799, 119.382, 121.044, 122.683, 123.843,
238 123.798, 124.567, 122.571, 122.951, 124.179, 124.977, 128.058, 131.607, 135.404, 139.132, 142.167, 146.221, 150.063, 154.225,
239 158.159, 160.614, 164.357, 168.218, 171.728, 174.947, 177.505, 181.286, 183.258, 186.901, 192.562, 189.001, 189.092, 185.551,
240 181.202, 177.128, 168.268, 159.719, 156.125, 151.878, 149.137, 144.408, 141.419, 138.169, 133.679, 130.259, 128.061, 126.093,
241 123.544, 121.135, 119.087, 116.403, 113.911, 112.066, 109.821, 105.606, 104.274, 101.77, 100.493, 99.9533, 98.9401, 98.3901,
242 97.6717, 97.3916, 97.613, 98.0975, 98.6056, 99.3877, 99.9002, 100.231, 100.728, 101.391, 102.93, 105.539, 106.374, 107.412,
243 108.678, 108.733, 111.575, 110.964, 113.989, 117.078, 120.737, 124.857, 129.56, 134.323, 138.298, 142.47, 146.55, 149.97,
244 153.422, 153.182, 152.819, 152.886, 152.96, 152.822, 152.713, 152.624, 152.621, 152.623, 152.681, 152.86, 153.001, 153.134,
245 153.335, 153.551, 153.775, 154.031, 154.275, 154.49, 154.74, 155.03, 155.233, 155.398, 155.545, 155.667, 155.785, 155.915,
246 156.003, 156.114, 156.24, 156.337, 156.407, 156.474, 156.524, 162.991, 163.283, 163.632, 163.845, 164.055, 164.256, 164.851,
247 165.542, 165.713, 165.796, 165.9, 166.052, 166.246, 166.473, 166.686, 166.917, 167.129, 167.341, 167.53, 167.75, 167.932,
248 168.201, 168.312, 168.131, 168.228, 168.412, 168.64, 168.918, 169.075, 169.257, 169.399, 169.53, 169.691, 170.004, 160.519,
249 160.606, 160.61, 160.726, 160.905, 161.082, 161.201, 161.352, 161.552, 161.677, 161.826, 161.982, 162.189, 162.366, 162.52,
250 162.655, 162.829, 163.002, 163.184, 163.36, 163.508, 163.683, 163.916, 164.574, 164.938, 165.128};
251 constexpr double dead_cryo0[277] = {
252 0.6732, 0.6734, 0.6735, 0.6737, 0.6739, 0.6743, 0.6746, 0.6751, 0.6757, 0.6763, 0.677, 0.6777, 0.6785, 0.6794, 0.6804, 0.6813,
253 0.6825, 0.6836, 0.6848, 0.6861, 0.6874, 0.6889, 0.6904, 0.6919, 0.6936, 0.6952, 0.697, 0.6988, 0.7007, 0.7028, 0.7048, 0.707,
254 0.7091, 0.7114, 0.7137, 0.7162, 0.7186, 0.7212, 0.7237, 0.7265, 0.7293, 0.7321, 0.735, 0.7379, 0.741, 0.7442, 0.7474, 0.7507,
255 0.754, 0.7575, 0.7609, 0.7646, 0.7682, 0.7719, 0.7758, 0.7796, 0.7836, 0.7876, 0.7918, 0.796, 0.8002, 0.8046, 0.8091, 0.8136,
256 0.8182, 0.8229, 0.8277, 0.8326, 0.8375, 0.8425, 0.8476, 0.8528, 0.8581, 0.8634, 0.869, 0.8744, 0.88, 0.8857, 0.8916, 0.8975,
257 0.9034, 0.9095, 0.9157, 0.9219, 0.9283, 0.9347, 0.9412, 0.9479, 0.9546, 0.9614, 0.9683, 0.9753, 0.9824, 0.9896, 0.997, 1.0043,
258 1.0119, 1.0194, 1.0271, 1.0349, 1.0429, 1.0509, 1.059, 1.0672, 1.0756, 1.084, 1.0925, 1.1012, 1.1099, 1.1189, 1.1278, 1.137,
259 1.1461, 1.1556, 1.1649, 1.1746, 1.1842, 1.1941, 1.2039, 1.214, 1.2242, 1.2344, 1.2449, 1.2553, 1.266, 1.2768, 1.2877, 1.2988,
260 1.31, 1.3212, 1.3327, 1.3443, 1.3559, 1.3678, 1.3798, 1.3919, 1.4041, 1.4165, 1.429, 1.4418, 1.4545, 1.4675, 1.4807, 1.4938,
261 1.5073, 1.5208, 1.5346, 1.5485, 1.5625, 1.5766, 1.5601, 1.2354, 0.4486, 0.5128, 0.3507, 0.1844, 0.1841, 0.1837, 0.1835, 0.1832,
262 0.1828, 0.1826, 0.1823, 0.182, 0.1817, 0.1814, 0.1812, 0.181, 0.1807, 0.1804, 0.1802, 0.1799, 0.1797, 0.1795, 0.1793, 0.1791,
263 0.1789, 0.1787, 0.1784, 0.1782, 0.178, 0.1778, 0.1777, 0.1775, 0.1773, 0.1771, 0.177, 0.1768, 0.1766, 0.1764, 0.1763, 0.1762,
264 0.176, 0.1759, 0.1758, 0.1756, 0.1754, 0.1753, 0.1752, 0.1751, 0.1749, 0.1748, 0.1746, 0.1745, 0.1745, 0.1743, 0.1742, 0.174,
265 0.174, 0.1739, 0.1737, 0.1736, 0.1735, 0.1735, 0.1734, 0.1732, 0.1731, 0.1731, 0.173, 0.1729, 0.1728, 0.1728, 0.1727, 0.1726,
266 0.1725, 0.1725, 0.1724, 0.1722, 0.1722, 0.1721, 0.172, 0.172, 0.1719, 0.1719, 0.1718, 0.1717, 0.1716, 0.1716, 0.1716, 0.1715,
267 0.1714, 0.1714, 0.1713, 0.1713, 0.1712, 0.1711, 0.1711, 0.1711, 0.1711, 0.171, 0.1709, 0.1709, 0.1709, 0.1708, 0.1708, 0.1708,
268 0.1707, 0.1706, 0.1706, 0.1706, 0.1706, 0.1705, 0.1704, 0.1704, 0.1704, 0.1703, 0.1703, 0.1703, 0.1702, 0.1702, 0.1702, 0.1702,
269 0.1701, 0.1701, 0.1701, 0.17, 0.17};
270 constexpr double dead_cryo1[277] = {
271 1.6192, 1.6194, 1.6197, 1.6203, 6.3509, 2.5841, 1.9874, 1.3905, 1.3916, 1.3929, 1.3943, 1.3957, 1.3975, 1.3992, 1.4013, 1.4033,
272 1.4055, 2.0138, 2.6243, 4.3032, 2.6344, 2.0294, 1.4221, 1.4251, 1.4284, 1.432, 1.4357, 1.4394, 1.4433, 1.4475, 1.4516, 1.4561,
273 1.4605, 1.4652, 1.4702, 1.475, 2.117, 2.7636, 4.1287, 2.3159, 2.1483, 1.5079, 1.5139, 1.5199, 1.5261, 1.5326, 1.5393, 1.546,
274 1.553, 1.5599, 1.5673, 1.5747, 1.5821, 1.5899, 1.5977, 2.297, 2.3084, 3.9766, 3.3404, 2.3449, 2.345, 1.6573, 1.6664, 1.6757,
275 1.6853, 1.6949, 1.7047, 1.7148, 1.7249, 1.7352, 1.7457, 1.7564, 1.7673, 1.7784, 2.2539, 2.576, 3.3648, 3.5384, 3.0253, 2.6439,
276 1.8606, 1.8733, 1.8859, 1.8989, 1.9119, 1.9251, 1.9387, 1.9521, 1.966, 1.98, 1.9943, 2.0088, 2.0235, 2.0383, 2.937, 3.214,
277 3.4082, 3.6655, 3.026, 3.0489, 2.148, 2.1643, 2.1812, 2.198, 3.6572, 3.7318, 3.5266, 3.3742, 3.534, 3.7906, 3.7431, 4.3293,
278 4.4603, 2.9506, 3.8466, 3.5582, 3.078, 2.4882, 2.2279, 1.7277, 1.2731, 1.0272, 0.7105, 0.3933, 0.2912, 0.2902, 0.2893, 0.2883,
279 0.2875, 0.2866, 0.2857, 0.2849, 0.284, 0.2834, 2.6585, 2.6785, 2.6989, 2.7197, 2.7405, 2.7614, 2.7828, 2.8045, 2.8265, 2.8488,
280 2.8712, 2.8939, 2.9114, 2.7992, 2.9324, 3.4047, 0.414, 0.252, 0.1168, 0.9172, 0.9156, 0.914, 0.9123, 0.9108, 0.9091, 0.9076,
281 0.9062, 0.9048, 0.9033, 0.9018, 0.9006, 0.8993, 0.8979, 0.8966, 0.8955, 0.8942, 0.893, 0.8918, 0.8907, 0.8896, 0.8886, 0.8874,
282 0.8865, 0.8853, 0.8843, 0.8833, 0.8823, 0.9288, 0.9279, 0.927, 0.926, 0.9251, 0.9242, 0.9233, 0.9224, 0.9217, 0.9208, 0.9199,
283 0.9192, 0.9184, 0.9176, 0.9169, 0.9162, 0.9154, 0.9148, 0.9141, 0.9135, 0.9128, 0.9121, 0.9115, 0.9109, 0.9103, 0.9097, 0.9091,
284 0.9086, 0.9079, 0.9074, 0.9069, 0.9064, 0.9059, 0.9053, 0.9048, 0.9043, 0.9038, 0.9034, 0.903, 0.9025, 0.9022, 0.9017, 0.9013,
285 0.9009, 0.9004, 0.9002, 0.8996, 0.8994, 0.899, 0.8986, 0.8982, 0.8979, 0.8976, 0.8971, 0.897, 0.8966, 0.8963, 0.896, 0.8957,
286 0.8954, 0.8951, 0.8948, 0.8945, 0.8942, 0.8941, 0.8937, 0.8934, 0.8932, 0.8929, 0.8927, 0.8924, 0.8922, 0.892, 0.8917, 0.8915,
287 0.8914, 0.8912, 0.891, 0.8907, 0.8905, 0.8903, 0.8902, 0.89, 0.8898, 0.8897, 0.8894, 0.8893, 0.889, 0.8889, 0.8888, 0.8886,
288 0.8886, 0.8883, 0.8882, 0.888, 0.8879};
289 constexpr double dead_cryo2[277] = {
290 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
291 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
292 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
293 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
294 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
295 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.2895, 0.9101, 1.5272, 1.8249, 1.8279, 1.8312,
296 1.8348, 1.8385, 1.8424, 1.8466, 1.851, 1.8555, 1.8604, 1.5309, 0, 0, 0, 0, 0, 0, 0, 0,
297 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0.0763, 0.6379, 0.8382, 1.0003, 1.2294,
298 1.5831, 2.0439, 2.3796, 2.6456, 2.9827, 3.3447, 1.141, 1.138, 1.1349, 1.132, 1.1291, 1.1264, 1.1236, 1.1208, 1.1183, 1.1156,
299 1.1132, 1.1108, 1.1084, 1.106, 1.1037, 1.1014, 4.514, 4.4934, 4.5226, 3.5091, 3.5412, 3.5737, 3.6066, 3.6398, 3.6734, 3.7074,
300 3.7418, 3.7764, 3.8115, 4.423, 5.2744, 6.1316, 6.6266, 6.6847, 6.7435, 6.803, 9.6759, 9.0086, 8.1971, 7.7761, 7.5583, 7.7543,
301 8.5776, 9.3644, 9.4149, 9.466, 9.4281, 8.6493, 7.8651, 7.5418, 7.2376, 7.3746, 7.4009, 7.3938, 7.3869, 7.3801, 7.3735, 7.3671,
302 7.3608, 7.3546, 7.3485, 7.3426, 7.3367, 7.331, 7.3254, 7.3199, 7.3146, 7.3093, 7.3042, 7.2991, 7.2943, 7.2894, 7.2846, 7.28,
303 7.2754, 7.271, 7.2666, 7.2624, 7.2582, 7.2541, 7.2501, 7.2461, 7.2423, 7.2386, 7.2348, 7.2312, 7.2276, 7.2242, 7.2207, 7.2174,
304 7.2141, 7.211, 7.2078, 7.2047, 7.2017, 7.1987, 7.1958, 7.1931, 7.1902, 7.1875, 7.1848, 7.1822, 7.1796, 7.1772, 7.1747, 7.1723,
305 7.1699, 7.1676, 7.1653, 7.1632, 7.161, 7.1588, 7.1567, 7.1547, 7.1527, 7.1507, 7.1487, 7.1469, 7.145, 7.1433, 7.1415, 7.1397,
306 7.138, 7.1364, 7.1347, 7.1331, 7.1315, 7.13, 7.1284, 7.1269, 7.1255, 7.1241, 7.1227, 7.1213, 7.12, 7.1187, 7.1175, 7.1162,
307 7.1149, 7.1137, 7.1126, 7.1114, 7.1102};
308 constexpr double dead_girder[277] = {
309 11.9514, 11.9525, 11.9549, 11.9584, 11.9632, 11.9693, 11.9764, 11.9848, 11.9943, 12.0051, 12.0171, 12.0303, 12.0446, 12.0602,
310 12.077, 12.095, 12.1142, 12.1346, 12.1563, 12.1791, 12.2032, 12.2285, 12.255, 12.2827, 12.3116, 12.3418, 12.3732, 12.4059,
311 12.4398, 12.475, 12.5113, 12.549, 12.5879, 12.6281, 12.6695, 12.7122, 12.7561, 12.8014, 12.8479, 12.8957, 12.9448, 12.9952,
312 13.0469, 13.0998, 13.1541, 13.2098, 13.2667, 13.325, 13.3845, 13.4455, 13.5078, 13.5714, 13.6364, 13.7028, 13.7705, 13.8396,
313 13.91, 13.9818, 14.0552, 14.1298, 14.206, 14.2835, 12.3788, 7.5701, 7.4957, 6.6037, 3.8717, 0, 0, 0,
314 0, 3.5585, 7.2483, 7.505, 8.4513, 10.9561, 15.539, 15.7235, 15.8265, 15.9311, 16.0372, 16.1449, 16.2543, 16.3653,
315 16.4779, 16.5921, 16.7081, 16.8256, 16.9449, 17.0659, 17.1886, 17.313, 17.4392, 17.567, 17.6966, 17.8281, 17.9612, 18.0962,
316 18.233, 18.3716, 18.5121, 18.6543, 18.7985, 18.9446, 19.0925, 19.2424, 19.3942, 19.5479, 19.7036, 19.8612, 20.0208, 20.1824,
317 20.3461, 20.5118, 20.6795, 20.8493, 21.0212, 16.6712, 11.697, 11.2949, 10.114, 9.906, 8.5687, 3.7101, 0, 0,
318 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
319 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
320 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
321 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
322 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
323 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
324 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
325 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
326 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
327 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
328 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0};
329 for (int i = 0; i < 277; i++) {
330 m_dead_eta[i] = dead_eta[i];
331 m_dead_x01[i] = dead_x01[i];
332 m_dead_x02[i] = dead_x02[i];
333 m_dead_x0tot[i] = dead_x0tot[i];
334 m_dead_cryo0[i] = dead_cryo0[i];
335 m_dead_cryo1[i] = dead_cryo1[i];
336 m_dead_cryo2[i] = dead_cryo2[i];
337 m_dead_girder[i] = dead_girder[i];
338 }
339 m_etaGranularity = 0.01;
340 //---- End of material map arrays initialization -----//
341
342 return StatusCode::SUCCESS;
343 }
344 CaloEnergy* MuidCaloEnergyParam::meanParametrizedEnergy(double trackMomentum, double eta, double phi) const {
345 CaloEnergy* caloEnergy = nullptr;
346 double meanEloss = meanEnergyLoss(eta, trackMomentum);
347 std::pair<double, double> parametrizedErrors = meanEnergyLossError(eta, trackMomentum);
348
349 ATH_MSG_DEBUG(" Mean_deposited_energy = " << meanEloss / Units::GeV);
350
351 // cosmics: energy sign flip in upper hemisphere
352 if (m_cosmics && phi > 0.) {
353 caloEnergy = new CaloEnergy(-meanEloss, 0.5 * (parametrizedErrors.first + parametrizedErrors.second), parametrizedErrors.second,
354 parametrizedErrors.first, CaloEnergy::Parametrized);
355 } else {
356 caloEnergy = new CaloEnergy(meanEloss, 0.5 * (parametrizedErrors.first + parametrizedErrors.second), parametrizedErrors.first,
357 parametrizedErrors.second, CaloEnergy::Parametrized);
358 }
359 // set parametrized energy loss in em
360 double emLoss = caloEnergy->deltaE() * emMopFraction(eta);
361 caloEnergy->set_fsrCandidateEnergy(emLoss);
362 return caloEnergy;
363 }
364
365 CaloEnergy* MuidCaloEnergyParam::mopParametrizedEnergy(double trackMomentum, double eta, double phi) const {
366 CaloEnergy* caloEnergy = nullptr;
367
368 // mop and mean energy deposition
369 double mopEloss = mopEnergyLoss(eta, trackMomentum);
370 double meanEloss = meanEnergyLoss(eta, trackMomentum);
371
372 ATH_MSG_DEBUG(" Most probable (mean) deposited energy " << mopEloss / Units::GeV << " (" << meanEloss / Units::GeV << ")");
373
374 // mop energy deposition uncertainty (Gaussian fit around mop value +- 2 sigma)
375 double mopError = mopEnergyLossError(eta, trackMomentum);
376
377 // scale low transverse momentum to allow for fluctuations
378 double pt = trackMomentum / std::cosh(eta);
379 double scale = 1.;
380 if (pt < FifteenGeV) scale = (2.5 - pt * (1.5 / (FifteenGeV)));
381
382 // asymmetric error to allow for Landau
383
384 constexpr double alpha = 0.2;
385 constexpr double beta = 300. * Units::MeV;
386 double mopErrorPlus = mopError + alpha * (meanEloss - mopEloss) + beta;
387
388 // cosmics: energy sign flip in upper hemisphere
389 if (m_cosmics && phi > 0.) {
390 caloEnergy = new CaloEnergy(-mopEloss, scale * mopError, scale * mopErrorPlus, scale * mopError, CaloEnergy::Parametrized);
391 } else {
392 caloEnergy = new CaloEnergy(mopEloss, scale * mopError, scale * mopError, scale * mopErrorPlus, CaloEnergy::Parametrized);
393 }
394 // set parametrized energy loss in em
395 double emLoss = caloEnergy->deltaE() * emMopFraction(eta);
396 caloEnergy->set_fsrCandidateEnergy(emLoss);
397 return caloEnergy;
398 }
399
400 CaloEnergy* MuidCaloEnergyParam::mopPeakEnergy(double trackMomentum, double eta, double phi) const {
401 CaloEnergy* caloEnergy = nullptr;
402
403 // mop and mean energy deposition
404 double mopEloss = mopEnergyLoss(eta, trackMomentum);
405 double meanEloss = meanEnergyLoss(eta, trackMomentum);
406
407 ATH_MSG_DEBUG(" Most probable (mean) deposited energy " << mopEloss / Units::GeV << " (" << meanEloss / Units::GeV << ")");
408
409 // mop energy deposition uncertainty (Gaussian fit around mop value +- 2 sigma)
410 double mopError = mopEnergyLossError(eta, trackMomentum);
411
412 // scale low transverse momentum to allow for fluctuations
413 double pt = trackMomentum / std::cosh(eta);
414 double scale = 1.;
415 if (pt < FifteenGeV) scale = (2.5 - pt * (1.5 / (FifteenGeV)));
416
417 // asymmetric error to allow for Landau
418 constexpr double alpha = 0.1;
419 constexpr double beta = 150. * Units::MeV;
420 double mopErrorPlus = mopError + alpha * (meanEloss - mopEloss) + beta;
421
422 // cosmics: energy sign flip in upper hemisphere
423 if (m_cosmics && phi > 0.) {
424 caloEnergy = new CaloEnergy(-mopEloss, scale * mopError, scale * mopErrorPlus, scale * mopError, CaloEnergy::MOP);
425 } else {
426 caloEnergy = new CaloEnergy(mopEloss, scale * mopError, scale * mopError, scale * mopErrorPlus, CaloEnergy::MOP);
427 }
428 // set parametrized energy loss in em
429 double emLoss = caloEnergy->deltaE() * emMopFraction(eta);
430 caloEnergy->set_fsrCandidateEnergy(emLoss);
431 return caloEnergy;
432 }
433
434 CaloEnergy* MuidCaloEnergyParam::mopDepositedEnergy(double trackMomentum, double eta, double phi) const {
435 CaloEnergy* caloEnergy = nullptr;
436
437 std::pair<double, double> eta_phi;
438 eta_phi.first = eta;
439 eta_phi.second = phi;
440
441 double mop_eloss = mopEnergyLoss(eta_phi.first, trackMomentum);
442
443 // double mop_eloss = symmetricMopEnergyLoss(eta_phi.first,trackMomentum);
444 double mop_eloss_error = mopEnergyLossError(eta_phi.first, trackMomentum);
445
446 ATH_MSG_DEBUG(" Most probable deposited energy... " << mop_eloss / Units::GeV);
447
448 // cosmics: energy sign flip in upper hemisphere
449 if (m_cosmics && phi > 0.) {
450 caloEnergy = new CaloEnergy(-mop_eloss, mop_eloss_error, mop_eloss_error, mop_eloss_error, CaloEnergy::Parametrized);
451 } else {
452 caloEnergy = new CaloEnergy(mop_eloss, mop_eloss_error, mop_eloss_error, mop_eloss_error, CaloEnergy::Parametrized);
453 }
454 // set parametrized energy loss in em
455 double emLoss = caloEnergy->deltaE() * emMopFraction(eta);
456 caloEnergy->set_fsrCandidateEnergy(emLoss);
457 return caloEnergy;
458 }
459
460 MuidCaloEnergyParam::BinsWeights MuidCaloEnergyParam::etaFixedBin(double eta) const {
461 MuidCaloEnergyParam::BinsWeights BW{};
462 eta = std::abs(eta);
463 double offsetEta1 = eta - m_etaOffset;
464 double offsetEta2 = eta + m_etaOffset;
465 BW.etaWeight2 = 0.5;
466 if (offsetEta2 < m_binWidth) // first (i.e. central) bin #0 - no interpolation
467 {
468 BW.etaBin1 = 0;
469 BW.etaBin2 = 0;
470 } else if (offsetEta1 >= 2.50) // last bin #25 (eta > 2.5) - no interpolation
471 {
472 BW.etaBin1 = 25;
473 BW.etaBin2 = 25;
474 } else if (m_smoothingFraction > 0.) {
475 double eta1 = floor(m_inverseWidth * offsetEta1);
476 BW.etaBin1 = static_cast<int>(eta1);
477 double eta2 = floor(m_inverseWidth * offsetEta2);
478 BW.etaBin2 = static_cast<int>(eta2);
479 if (BW.etaBin2 > BW.etaBin1) BW.etaWeight2 = (m_inverseWidth * offsetEta2 - eta2) / m_smoothingFraction;
480 } else {
481 double eta1 = floor(m_inverseWidth * eta);
482 BW.etaBin1 = static_cast<int>(eta1);
483 BW.etaBin2 = BW.etaBin1;
484 }
485
486 BW.etaWeight1 = 1. - BW.etaWeight2;
487 return BW;
488 }
489
490 MuidCaloEnergyParam::BinsWeights MuidCaloEnergyParam::etaVariableBin(double eta) const {
491 eta = std::abs(eta);
492 if (eta < 1.05) // barrel: contains the variable bin region
493 {
494 if (eta < 0.75) // identical to fixed binning
495 {
496 return etaFixedBin(eta);
497 } else if (eta < 0.80) // wide bin#7
498 {
499 return etaFixedBin(0.75);
500 } else if (eta < 1.00) // similar to fixed binning (with half bin offset)
501 {
502 return etaFixedBin(eta - 0.5 * m_binWidth);
503 } else // wide bin#9
504 {
505 return etaFixedBin(0.95);
506 }
507 } else // endcap: fixed binning (with full bin offset)
508 {
509 return etaFixedBin(eta - m_binWidth);
510 }
511 }
512
513 double MuidCaloEnergyParam::meanEnergyLoss(double eta, double momentum) const {
514 // Mean Energy Loss parametrization with asymmetric gaussian fit
515 // tabulated coefficients for function [0]+[1]*Log(0.0067*p*p)+[2]*p with GeV units
516 MuidCaloEnergyParam::BinsWeights BW = etaVariableBin(eta);
517 double p0 = BW.etaWeight1 * m_meanEnergyLossP0[BW.etaBin1] + BW.etaWeight2 * m_meanEnergyLossP0[BW.etaBin2];
518 double p1 = BW.etaWeight1 * m_meanEnergyLossP1[BW.etaBin1] + BW.etaWeight2 * m_meanEnergyLossP1[BW.etaBin2];
519 double p2 = BW.etaWeight1 * m_meanEnergyLossP2[BW.etaBin1] + BW.etaWeight2 * m_meanEnergyLossP2[BW.etaBin2];
520 double parametrisedDeposit = p0 * Units::GeV + p1 * Units::GeV * std::log(0.0067 * std::pow(momentum * PerGeV, 2)) + p2 * momentum;
521
522 // additional offset from high-statistics Z->mumu MC (measured by Peter K 30/11/2011, 03/2012)
523 // NOTE: constant bin width used here
524 BW = etaFixedBin(eta);
525 // mc11 tune
526 // double fixFromPeter[26] = { -0.203781 , -0.146775 , -0.0952301 , -0.186141 , -0.180687 ,
527 // -0.225095 , -0.195912 , -0.221158 , -0.198556 , -0.271551 ,
528 // -0.200737 , -0.151854 , -0.0799826 , -0.237688 , -0.255126 ,
529 // -0.313466 , -0.35815 , -0.334661 , -0.458948 , -0.469967 ,
530 // -0.516059 , -0.522978 , -0.53896 , -0.550683 , -0.51722 ,
531 // -0.543858 };
532 // mc12 tune
533 constexpr double fixFromPeter[26] = {-0.159195, -0.039180, 0.071644, -0.118753, -0.144227, -0.182778, -0.158818,
534 -0.188853, -0.180997, -0.269571, -0.227767, -0.191740, -0.121575, -0.304938,
535 -0.296636, -0.312171, -0.352982, -0.303020, -0.357250, -0.396249, -0.433535,
536 -0.417214, -0.424093, -0.415557, -0.325372, -0.553009};
537 double fix = BW.etaWeight1 * fixFromPeter[BW.etaBin1] + BW.etaWeight2 * fixFromPeter[BW.etaBin2];
538 parametrisedDeposit += fix * Units::GeV;
539
540 return parametrisedDeposit;
541 }
542
543 std::pair<double, double> MuidCaloEnergyParam::meanEnergyLossError(double eta, double momentum) const {
544 MuidCaloEnergyParam::BinsWeights BW = etaVariableBin(eta);
545 double left_p0 = BW.etaWeight1 * m_meanEnergyLossErrorLeftP0[BW.etaBin1] + BW.etaWeight2 * m_meanEnergyLossErrorLeftP0[BW.etaBin2];
546 double left_p1 = BW.etaWeight1 * m_meanEnergyLossErrorLeftP1[BW.etaBin1] + BW.etaWeight2 * m_meanEnergyLossErrorLeftP1[BW.etaBin2];
547 double right_p0 =
548 BW.etaWeight1 * m_meanEnergyLossErrorRightP0[BW.etaBin1] + BW.etaWeight2 * m_meanEnergyLossErrorRightP0[BW.etaBin2];
549 double right_p1 =
550 BW.etaWeight1 * m_meanEnergyLossErrorRightP1[BW.etaBin1] + BW.etaWeight2 * m_meanEnergyLossErrorRightP1[BW.etaBin2];
551 double sigma_left = left_p0 * Units::GeV + left_p1 * momentum;
552 double sigma_right = right_p0 * Units::GeV + right_p1 * momentum;
553
554 return std::make_pair(sigma_left, sigma_right);
555 }
556
557 double MuidCaloEnergyParam::mopEnergyLoss(double eta, double momentum) const {
558 // Mop Energy Loss parametrization with asymmetric gaussian fit in the mop region
559 // tabulated coefficients for function [0]+[1]*Log(0.0067*p*p)+[2]*p with GeV units
560 MuidCaloEnergyParam::BinsWeights BW = etaVariableBin(eta);
561 double p0 = BW.etaWeight1 * m_mopEnergyLossP0[BW.etaBin1] + BW.etaWeight2 * m_mopEnergyLossP0[BW.etaBin2];
562 double p1 = BW.etaWeight1 * m_mopEnergyLossP1[BW.etaBin1] + BW.etaWeight2 * m_mopEnergyLossP1[BW.etaBin2];
563 double p2 = BW.etaWeight1 * m_mopEnergyLossP2[BW.etaBin1] + BW.etaWeight2 * m_mopEnergyLossP2[BW.etaBin2];
564 double parametrisedDeposit = p0 * Units::GeV + p1 * Units::GeV * std::log(0.0067 * std::pow(momentum * PerGeV, 2)) + p2 * momentum;
565
566 // additional offset from high-statistics Z->mumu MC (measured by Peter K 30/11/2011)
567 // NOTE: constant bin width used here
568 BW = etaFixedBin(eta);
569 // mc11 tune
570 // double fixFromPeter[26] = { -0.203781 , -0.146775 , -0.0952301 , -0.186141 , -0.180687 ,
571 // -0.225095 , -0.195912 , -0.221158 , -0.198556 , -0.271551 ,
572 // -0.200737 , -0.151854 , -0.0799826 , -0.237688 , -0.255126 ,
573 // -0.313466 , -0.35815 , -0.334661 , -0.458948 , -0.469967 ,
574 // -0.516059 , -0.522978 , -0.53896 , -0.550683 , -0.51722 ,
575 // -0.543858 };
576 // mc12 tune
577 constexpr double fixFromPeter[26] = {-0.159195, -0.039180, 0.071644, -0.118753, -0.144227, -0.182778, -0.158818,
578 -0.188853, -0.180997, -0.269571, -0.227767, -0.191740, -0.121575, -0.304938,
579 -0.296636, -0.312171, -0.352982, -0.303020, -0.357250, -0.396249, -0.433535,
580 -0.417214, -0.424093, -0.415557, -0.325372, -0.553009};
581
582 double fix = BW.etaWeight1 * fixFromPeter[BW.etaBin1] + BW.etaWeight2 * fixFromPeter[BW.etaBin2];
583 parametrisedDeposit += fix * Units::GeV;
584
585 return parametrisedDeposit;
586 }
587
588 double MuidCaloEnergyParam::mopEnergyLossError(double eta, double momentum) const {
589 eta = std::abs(eta);
590 MuidCaloEnergyParam::BinsWeights BW = etaVariableBin(eta);
591 double sigma_p0 = BW.etaWeight1 * m_mopEnergyLossErrorP0[BW.etaBin1] + BW.etaWeight2 * m_mopEnergyLossErrorP0[BW.etaBin2];
592 double sigma_p1 = BW.etaWeight1 * m_mopEnergyLossErrorP1[BW.etaBin1] + BW.etaWeight2 * m_mopEnergyLossErrorP1[BW.etaBin2];
593 double sigma = sigma_p0 * Units::GeV + sigma_p1 * momentum;
594
595 // fix for mc11 as underestimate in barrel, overestimate in endcap
596 if (eta < 1.7) {
597 sigma *= 1.1;
598 } else {
599 sigma *= 0.9;
600 }
601
602 return sigma;
603 }
604
605 double MuidCaloEnergyParam::symmetricMopEnergyLoss(double eta, double momentum) const {
606 // Mop Energy Loss parametrization with symmetric gaussian fit in the mop region
607 // new coefficients with symmetric fits
608 MuidCaloEnergyParam::BinsWeights BW = etaVariableBin(eta);
609 double p0 = BW.etaWeight1 * m_mopSymmetricEnergyLossP0[BW.etaBin1] + BW.etaWeight2 * m_mopSymmetricEnergyLossP0[BW.etaBin2];
610 double p1 = BW.etaWeight1 * m_mopSymmetricEnergyLossP1[BW.etaBin1] + BW.etaWeight2 * m_mopSymmetricEnergyLossP1[BW.etaBin2];
611 double p2 = BW.etaWeight1 * m_mopSymmetricEnergyLossP2[BW.etaBin1] + BW.etaWeight2 * m_mopSymmetricEnergyLossP2[BW.etaBin2];
612 double parametrisedDeposit = p0 * Units::GeV + p1 * Units::GeV * log(0.0067 * std::pow(momentum * PerGeV, 1)) + p2 * momentum;
613 return parametrisedDeposit;
614 }
615
616 double MuidCaloEnergyParam::x0mapInertMaterial(double eta) const {
617 int index = static_cast<int>(floor(std::abs(eta) / m_etaGranularity));
618 if (index >= 277) index = 276;
619
620 double x01 = m_dead_x01[index];
621 // double x02 = m_dead_x02[index];
622 double x0tot = m_dead_x0tot[index];
623 double x0cryo1 = m_dead_cryo1[index];
624 double x0cryo2 = m_dead_cryo2[index];
625 double x0girder = m_dead_girder[index];
626 double sum = x01 + x0cryo1 + x0cryo2 + x0girder;
627 return sum / x0tot;
628 }
629
630 double MuidCaloEnergyParam::x0mapEmMaterial(double eta) const {
631 int index = static_cast<int>(floor(std::abs(eta) / m_etaGranularity));
632 if (index >= 277) index = 276;
633
634 double x01 = m_dead_x01[index];
635 double x02 = m_dead_x02[index];
636 double x0tot = m_dead_x0tot[index];
637 double x0cryo1 = m_dead_cryo1[index];
638 // double x0cryo2 = m_dead_cryo2[index];
639 // double x0girder = m_dead_girder[index];
640 double sum = x02 - x01 - x0cryo1;
641 return sum / x0tot;
642 }
643
644 double MuidCaloEnergyParam::x0mapHecMaterial(double eta) const {
645 int index = static_cast<int>(floor(std::abs(eta) / m_etaGranularity));
646 if (index >= 277) index = 276;
647
648 // double x01 = m_dead_x01[index];
649 double x02 = m_dead_x02[index];
650 double x0tot = m_dead_x0tot[index];
651 // double x0cryo1 = m_dead_cryo1[index];
652 double x0cryo2 = m_dead_cryo2[index];
653 double x0girder = m_dead_girder[index];
654 double sum = x0tot - x02 - x0cryo2 - x0girder;
655 return sum / x0tot;
656 }
657
658 double MuidCaloEnergyParam::caloCompartmentDepth(int icalo, int icomp) const {
659 /* ------------------
660 icalo calorimeter
661 0 Tile
662 1 LArHEC
663 2 LArEM
664 ------------------ */
665 if (icalo == 0)
666 return caloCompartmentDepthTile(icomp);
667 else if (icalo == 1)
668 return caloCompartmentDepthLArHEC(icomp);
669 else if (icalo == 2)
670 return caloCompartmentDepthLArEM(icomp);
671 else
672 ATH_MSG_WARNING(" Incorrect selection of calorimeter... returning 0.");
673 return 0.;
674 }
675
676 double MuidCaloEnergyParam::caloCompartmentDepthTile(int icomp) {
677 // depth in absorption lengths (from Tile TDR p 22) eta = 0
678 constexpr double comp[4] = {1.5, 4.2, 1.9, 0.}; // assign 0. to ITC for the moment
679 double total = 0.;
680 for (int i = 0; i < 4; i++) total += comp[i];
681 if (icomp >= 0 && icomp < 4) { return comp[icomp] / total; }
682 return 0.;
683 }
684
685 double MuidCaloEnergyParam::caloCompartmentDepthLArHEC(int icomp) {
686 // depth in absorption lengths (from Liquid Argon TDR Table 1-3 p 11)
687 constexpr double comp[4] = {1.4, 2.9, 5.7 / 2., 5.7 / 2.}; // divide 5.7 rad lengths equally to 3rd and 4th sampling
688 double total = 0.;
689 for (int i = 0; i < 4; i++) total += comp[i];
690 const double inv_total = 1. / total;
691 if (icomp >= 0 && icomp < 4) { return comp[icomp] * inv_total; }
692 return 0.;
693 }
694
695 double MuidCaloEnergyParam::caloCompartmentDepthLArEM(int icomp) {
696 // depth in radiation lengths (from Liquid Argon TDR p 5) eta = 0
697 constexpr double comp[4] = {0.2, 4.3, 16., 2.};
698 double total = 0.;
699 for (int i = 0; i < 4; i++) total += comp[i];
700 if (icomp >= 0 && icomp < 4) { return comp[icomp] / total; }
701 return 0.;
702 }
703
704 double MuidCaloEnergyParam::emMopFraction(double eta) const {
705 int ieta = static_cast<int>(std::abs(eta) / 0.10);
706 if (ieta > 25) ieta = 25;
707
708 return m_emMopFraction[ieta];
709 }
710
711 double MuidCaloEnergyParam::em1MopFraction(double eta) const {
712 // FIXME!!
713 double em1Fraction = 0.1 * emMopFraction(eta);
714 return em1Fraction;
715 }
716
717} // namespace Rec
eta
Scalar eta() const
pseudorapidity method
Definition AmgMatrixBasePlugin.h:83
phi
Scalar phi() const
phi method
Definition AmgMatrixBasePlugin.h:67
ATH_MSG_WARNING
#define ATH_MSG_WARNING(x)
Definition AthMsgStreamMacros.h:32
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition AthMsgStreamMacros.h:29
Units.h
Wrapper to avoid constant divisions when using units.
AthAlgTool::AthAlgTool
AthAlgTool(const std::string &type, const std::string &name, const IInterface *parent)
Constructor with parameters:
Definition AthAlgTool.cxx:16
AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T, V, H > &t)
Definition AthCommonDataStore.h:145
CaloEnergy
class extending the basic Trk::EnergyLoss to describe the measured or parameterised muon energy loss ...
Definition CaloEnergy.h:28
CaloEnergy::set_fsrCandidateEnergy
void set_fsrCandidateEnergy(const float fs)
FSR Candidate Energy.
Definition CaloEnergy.h:118
Rec::MuidCaloEnergyParam::caloCompartmentDepthTile
static double caloCompartmentDepthTile(int)
Definition MuidCaloEnergyParam.cxx:676
Rec::MuidCaloEnergyParam::etaFixedBin
BinsWeights etaFixedBin(double) const
Definition MuidCaloEnergyParam.cxx:460
Rec::MuidCaloEnergyParam::mopParametrizedEnergy
CaloEnergy * mopParametrizedEnergy(double trackMomentum, double eta, double phi) const
IMuidCaloEnergyParam interface: to get the most probable (peak value) energy loss in the calorimeter ...
Definition MuidCaloEnergyParam.cxx:365
Rec::MuidCaloEnergyParam::caloCompartmentDepth
double caloCompartmentDepth(int, int) const
IMuidCaloEnergyParam interface: depth of calo compartment.
Definition MuidCaloEnergyParam.cxx:658
Rec::MuidCaloEnergyParam::mopEnergyLossError
double mopEnergyLossError(double, double) const
Definition MuidCaloEnergyParam.cxx:588
Rec::MuidCaloEnergyParam::symmetricMopEnergyLoss
double symmetricMopEnergyLoss(double, double) const
Definition MuidCaloEnergyParam.cxx:605
Rec::MuidCaloEnergyParam::MuidCaloEnergyParam
MuidCaloEnergyParam(const std::string &type, const std::string &name, const IInterface *parent)
Definition MuidCaloEnergyParam.cxx:31
Rec::MuidCaloEnergyParam::em1MopFraction
double em1MopFraction(double eta) const
IMuidCaloEnergyParam interface: fraction of the full mop deposition corresponding to the em presample...
Definition MuidCaloEnergyParam.cxx:711
Rec::MuidCaloEnergyParam::meanParametrizedEnergy
CaloEnergy * meanParametrizedEnergy(double trackMomentum, double eta, double phi) const
IMuidCaloEnergyParam interface: to get the mean muon energy loss in the calorimeter (with asymmetric ...
Definition MuidCaloEnergyParam.cxx:344
Rec::MuidCaloEnergyParam::x0mapHecMaterial
double x0mapHecMaterial(double) const
IMuidCaloEnergyParam interface: returns the fraction of hec calorimeter inert material between InDet ...
Definition MuidCaloEnergyParam.cxx:644
Rec::MuidCaloEnergyParam::etaVariableBin
BinsWeights etaVariableBin(double) const
Definition MuidCaloEnergyParam.cxx:490
Rec::MuidCaloEnergyParam::emMopFraction
double emMopFraction(double eta) const
IMuidCaloEnergyParam interface: fraction of the full mop deposition corresponding to the em.
Definition MuidCaloEnergyParam.cxx:704
Rec::MuidCaloEnergyParam::meanEnergyLoss
double meanEnergyLoss(double, double) const
Definition MuidCaloEnergyParam.cxx:513
Rec::MuidCaloEnergyParam::mopDepositedEnergy
CaloEnergy * mopDepositedEnergy(double trackMomentum, double eta, double phi) const
IMuidCaloEnergyParam interface: to get the most probable (peak value) energy loss and error in the ca...
Definition MuidCaloEnergyParam.cxx:434
Rec::MuidCaloEnergyParam::x0mapEmMaterial
double x0mapEmMaterial(double) const
IMuidCaloEnergyParam interface: returns the fraction of em calorimeter inert material between InDet a...
Definition MuidCaloEnergyParam.cxx:630
Rec::MuidCaloEnergyParam::caloCompartmentDepthLArEM
static double caloCompartmentDepthLArEM(int)
Definition MuidCaloEnergyParam.cxx:695
Rec::MuidCaloEnergyParam::mopPeakEnergy
CaloEnergy * mopPeakEnergy(double trackMomentum, double eta, double phi) const
IMuidCaloEnergyParam interface: to get the most probable (peak value) energy loss in the calorimeter ...
Definition MuidCaloEnergyParam.cxx:400
Rec::MuidCaloEnergyParam::caloCompartmentDepthLArHEC
static double caloCompartmentDepthLArHEC(int)
Definition MuidCaloEnergyParam.cxx:685
Rec::MuidCaloEnergyParam::x0mapInertMaterial
double x0mapInertMaterial(double) const
IMuidCaloEnergyParam interface: returns the fraction of inert material between InDet and MS.
Definition MuidCaloEnergyParam.cxx:616
Rec::MuidCaloEnergyParam::mopEnergyLoss
double mopEnergyLoss(double, double) const
Definition MuidCaloEnergyParam.cxx:557
Rec::MuidCaloEnergyParam::meanEnergyLossError
std::pair< double, double > meanEnergyLossError(double, double) const
Definition MuidCaloEnergyParam.cxx:543
Trk::EnergyLoss::deltaE
double deltaE() const
returns the
Rec
Gaudi Tools.
Definition FakeTrackBuilder.h:10
index
Definition index.py:1
type
Generated on for ATLAS Offline Software by doxygen 1.14.0