Belle II Software development
ECLCrystalData.cc
1/**************************************************************************
2 * basf2 (Belle II Analysis Software Framework) *
3 * Author: The Belle II Collaboration *
4 * *
5 * See git log for contributors and copyright holders. *
6 * This file is licensed under LGPL-3.0, see LICENSE.md. *
7 **************************************************************************/
8
9/* Own header. */
10#include <ecl/modules/eclBackgroundStudy/ECLCrystalData.h>
11
12/* Basf2 headers. */
13#include <framework/logging/Logger.h>
14
15/* ROOT headers. */
16#include <TMath.h>
17
18/* C++ headers. */
19#include <cmath>
20
21using namespace Belle2;
22
23//ECL has 16-fold symmetry in phi. Multiply this by 16 to get the total number of crystals at each theta location.
24int const ECLCrystalData::Ring[69] = {
25 3, 3, 4, 4, 4, 6, 6, 6, 6, 6, 6, 9, 9, //forward calorimeter
26 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, //barrel
27 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, //barrel
28 9, 9, 9, 9, 9, 9, //barrel
29 9, 9, 6, 6, 6, 6, 6, 4, 4, 4 //backward calorimeter
30};
31
32int const ECLCrystalData::sumPrevious[24] = {0, 3, 6, 10, 14, 18, 24, 30, 36, 42, 48, 54, 63, 72, 81, 90, 96, 102, 108, 114, 120, 124, 128, 132};
33
34//endcap crystal masses
35float const ECLCrystalData::CrystalMassEndcap[132] = {4.99, 5.01, 4.99, 5.82, 5.94, 5.83, 4.75, 4.76, 4.77, 4.74, 5.14, 5.16, 5.16,
36 5.14, 5.76, 5.8, 5.79, 5.78, 4.14, 4.28, 4.15, 4.16, 4.28, 4.14, 4.32, 4.34,
37 4.34, 4.34, 4.34, 4.32, 4.85, 4.87, 4.87, 4.87, 4.87, 4.85, 4.99, 5.02, 5.02,
38 5.02, 5.02, 4.99, 5.37, 5.39, 5.39, 5.39, 5.39, 5.37, 5.89, 5.92, 5.92, 5.92,
39 5.92, 5.9, 4.03, 4.12, 4.05, 4.05, 4.12, 4.05, 4.05, 4.12, 4.03, 4.16, 4.18,
40 4.18, 4.18, 4.18, 4.18, 4.18, 4.18, 4.16, 4.47, 4.49, 4.49, 4.49, 4.49, 4.49,
41 4.49, 4.49, 4.47, 4.26, 4.34, 4.28, 4.28, 4.34, 4.28, 4.28, 4.34, 4.26, 5.88,
42 5.91, 5.91, 5.91, 5.91, 5.88, 5.46, 5.48, 5.48, 5.48, 5.48, 5.46, 5.02, 5.04,
43 5.04, 5.04, 5.04, 5.02, 4.63, 4.65, 4.65, 4.65, 4.65, 4.63, 4.4, 4.53, 4.42,
44 4.42, 4.53, 4.4, 5.74, 5.77, 5.76, 5.74, 5.09, 5.11, 5.11, 5.09, 4.56, 4.58,
45 4.58, 4.56
46 };
47
48//Barrel crystal masses
49float const ECLCrystalData::CrystalMassBarrel[46] = {4.599, 4.619, 4.639, 4.659, 4.680, 4.701, 4.722, 4.744, 4.765,
50 4.786, 4.807, 4.828, 4.848, 4.868, 4.887, 4.906, 4.923, 4.940,
51 4.955, 4.969, 4.981, 4.992, 5.001, 5.008, 5.013, 5.016, 5.017,
52 5.016, 4.455, 4.455, 5.016, 5.017, 5.016, 5.013, 5.008, 5.001,
53 4.992, 4.981, 4.969, 4.955, 4.940, 4.923, 4.906, 4.887, 4.868,
54 4.848
55 };
56
57
58//Theta ID to theta(deg)
59double const ECLCrystalData::theta[69] = {13.9486666667, 15.4862, 17.110575, 18.757625, 20.255425, 21.8637333333,
60 23.4259666667, 24.96125, 26.4777833333, 27.9407, 29.36275, 30.8382555556,
61 32.2618444444, 33.667039, 35.062886, 36.509337, 38.007345, 39.557769,
62 41.161373, 42.818633, 44.530027, 46.295838, 48.115876, 49.990235, 51.91882,
63 53.900865, 55.935415, 58.021324, 60.157158, 62.3412, 64.571442, 66.8455,
64 69.160698, 71.513983, 73.902011, 76.321253, 78.767888, 81.237718, 83.726351,
65 86.229301, 88.741891, 90, 90, 91.258109, 93.770699, 96.273649, 98.762282,
66 101.232112, 103.678747, 106.097989, 108.486017, 110.839302, 113.154501,
67 115.428558, 117.658801, 119.842842, 121.978676, 124.064585, 126.099135,
68 128.2902222222, 130.4424444444, 132.6918333333, 134.9551666667, 137.3811666667,
69 139.9081666667, 142.5951666667, 145.374, 148.0185, 150.8355
70 };
71
73double const ECLCrystalData::EndcapRadius[132] = {
74// 3,3,4,4,4,6,6,6,6,6,6,9,9, //forward calorimeter
75 217.25, 217.25, 217.25,
76 218.72, 218.72, 218.72,
77 220.29, 220.29, 220.29, 220.29,
78 222.00, 222.00, 222.00, 222.00,
79 223.96, 223.96, 223.96, 223.96,
80 226.01, 226.01, 226.01, 226.01, 226.01, 226.01,
81 228.21, 228.21, 228.21, 228.21, 228.21, 228.21,
82 230.67, 230.67, 230.67, 230.67, 230.67, 230.67,
83 233.20, 233.20, 233.20, 233.20, 233.20, 233.20,
84 235.89, 235.89, 235.89, 235.89, 235.89, 235.89,
85 238.89, 238.89, 238.89, 238.89, 238.89, 238.89,
86 241.94, 241.94, 241.94, 241.94, 241.94, 241.94, 241.94, 241.94, 241.94,
87 245.16, 245.16, 245.16, 245.16, 245.16, 245.16, 245.16, 245.16, 245.16,
88// 9,9,6,6,6,6,6,4,4,4 //backward calorimeter
89 172.06, 172.06, 172.06, 172.06, 172.06, 172.06, 172.06, 172.06, 172.06,
90 165.69, 165.69, 165.69, 165.69, 165.69, 165.69, 165.69, 165.69, 165.69,
91 159.78, 159.78, 159.78, 159.78, 159.78, 159.78,
92 154.36, 154.36, 154.36, 154.36, 154.36, 154.36,
93 149.37, 149.37, 149.37, 149.37, 149.37, 149.37,
94 144.82, 144.82, 144.82, 144.82, 144.82, 144.82,
95 140.67, 140.67, 140.67, 140.67, 140.67, 140.67,
96 136.89, 136.89, 136.89, 136.89,
97 133.49, 133.49, 133.49, 133.49,
98 133.49, 133.49, 133.49, 133.49
99};
100
101
103double const ECLCrystalData::BarrelZ[46] = {210.15, 200.14, 190.44, 181.05, 171.95, 163.03, 154.45, 146.10, 137.98, 129.99, 122.27, 114.72, 107.34, 100.05, 92.96, 86.00, 79.15, 72.36, 65.72, 59.16, 52.69, 46.22, 39.87, 33.56, 27.29, 21.01, 14.81, 8.62, 2.78, -2.78, -14.81, -21.01, -27.29, -33.56, -39.87, -46.22, -52.69, -59.16, -65.72, -72.36, -79.15, -86.00, -92.95, -100.05, -107 - 28};
104
105
106
107
109ECLCrystalData::ECLCrystalData()
110{
111 m_cell_ID = -1;
112 m_phi_ID = -1;
113 m_theta_ID = -1;
114 m_phi_idx = -1;
115 m_theta_idx = -1;
116
117 m_mass = -1;
118 m_volume = -1;
119 m_PosR = -1;
120 m_PosP = -1;
121 m_PosX = -1;
122 m_PosY = -1;
123 m_PosZ = -1;
124}
125
126ECLCrystalData::ECLCrystalData(int cid)
127{
128 m_cell_ID = cid;
129 Mapping(cid);
130
131 Eval();
132
133}
134
135ECLCrystalData::ECLCrystalData(int tid, int pid)
136{
137 m_phi_idx = -1;
138 m_theta_idx = -1;
139 m_theta_ID = tid;
140 m_phi_ID = pid;
141 m_cell_ID = GetCellID(tid, pid);
142
143 Eval();
144}
145
146ECLCrystalData::~ECLCrystalData()
147{
148}
149
150
152void ECLCrystalData::Eval()
153{
154 //Crystal properties
155 m_mass = EvalMass();
156 m_volume = m_mass / 4.51;
157
158 //Crystal Location
159 m_PosP = EvalPhi();
160 m_PosR = EvalR();
161
162 m_PosZ = EvalZ();
163 m_PosX = EvalX();
164 m_PosY = EvalY();
165
166
167}
168
169double ECLCrystalData::EvalX()
170{
171 return m_PosR * cos(TMath::DegToRad() * m_PosP);
172}
173
174double ECLCrystalData::EvalY()
175{
176 return m_PosR * sin(TMath::DegToRad() * m_PosP);
177}
178
179
180//Determine mass of crystal
181float ECLCrystalData::EvalMass()
182{
183 if (IsEndCap()) {
184 return CrystalMassEndcap[GetCrystalIndex()];
185
186 } else {
187 return CrystalMassBarrel[GetCrystalIndex()];
188 }
189}
190
191int ECLCrystalData::GetCrystalIndex()
192{
193 if (m_theta_ID <= 12) { //forward ECL crystal
194 return m_phi_ID % Ring[m_theta_ID] + sumPrevious[m_theta_ID];
195
196 } else if (m_theta_ID >= 59) { //backward ECL crystal
197 return m_phi_ID % Ring[m_theta_ID] + sumPrevious[m_theta_ID - 46];
198
199 } else if (12 < m_theta_ID && m_theta_ID < 59) { //barrel ECL crystal
200 return m_theta_ID - 13;
201 }
202
203 B2ERROR("ERROR: invalid theta_ID=" << m_theta_ID);
204 return -1;
205}
206
207inline bool ECLCrystalData::IsEndCap()
208{
209 return m_theta_ID < 13 || m_theta_ID > 58 ;
210}
211
212inline bool ECLCrystalData::IsBarrel()
213{
214 return !IsEndCap();
215}
216
217
218inline int ECLCrystalData::GetCellID()
219{
220 return m_cell_ID;
221}
222
223double ECLCrystalData::EvalPhi()
224{
225 return 22.5 * m_phi_ID / Ring[m_theta_ID];
226}
227
228double ECLCrystalData::EvalZ()
229{
230 if (IsEndCap()) {
231 return EndcapRadius[GetCrystalIndex()] * cos(TMath::DegToRad() * theta[m_theta_ID]);
232 }
233
234 return BarrelZ[GetCrystalIndex()];
235}
236
237double ECLCrystalData::EvalR()
238{
239 if (IsEndCap()) {
240 return EndcapRadius[GetCrystalIndex()] * sin(TMath::DegToRad() * theta[m_theta_ID]);
241 }
242
243 return 147.7;
244}
245
246
247//Get CellID from theta and phi IDs
248int ECLCrystalData::GetCellID(int ThetaId, int PhiId)
249{
253 if (ThetaId < 13) {
254 const int forwRing[13] = {0, 3, 6, 10, 14, 18, 24, 30, 36, 42, 48, 54, 63 };
255 return forwRing[ThetaId] * 16 + PhiId;
256
257 } else if (ThetaId > 58) {
258 const int backRing[10] = {0, 9, 18, 24, 30, 36, 42, 48, 52, 56} ;
259 return ECLElementNumbers::c_NCrystalsForwardBarrel + backRing[ThetaId - 59] * 16 + PhiId;
260
261 }
262
263 return ECLElementNumbers::c_NCrystalsForward + 144 * (ThetaId - 13) + PhiId;
264}
265
266void ECLCrystalData::Mapping(int cid)
267{
268 if (cid < 0) {
269 B2ERROR("ECL ECLGeometryPar Mapping " << cid << ". Out of range.");
270
271 } else if (cid < 3 * 16) { //Forkward start
272 m_theta_ID = 0;
273 m_phi_ID = cid;
274 m_phi_idx = m_phi_ID % 3;
275 m_theta_idx = m_phi_ID / 3;
276 } else if (cid < 6 * 16) {
277 m_theta_ID = 1;
278 m_phi_ID = cid - 3 * 16;
279 m_phi_idx = m_phi_ID % 3 + 3;
280 m_theta_idx = m_phi_ID / 3;
281 } else if (cid < 10 * 16) {
282 m_theta_ID = 2;
283 m_phi_ID = cid - 6 * 16 ;
284 m_phi_idx = m_phi_ID % 4 + 6;
285 m_theta_idx = m_phi_ID / 4;
286 } else if (cid < 14 * 16) {
287 m_theta_ID = 3;
288 m_phi_ID = cid - 10 * 16 ;
289 m_phi_idx = m_phi_ID % 4 + 10;
290 m_theta_idx = m_phi_ID / 4;
291 } else if (cid < 18 * 16) {
292 m_theta_ID = 4;
293 m_phi_ID = cid - 14 * 16 ;
294 m_phi_idx = m_phi_ID % 4 + 14;
295 m_theta_idx = m_phi_ID / 4;
296 } else if (cid < 24 * 16) {
297 m_theta_ID = 5;
298 m_phi_ID = cid - 18 * 16 ;
299 m_phi_idx = m_phi_ID % 6 + 18;
300 m_theta_idx = m_phi_ID / 6;
301 } else if (cid < 30 * 16) {
302 m_theta_ID = 6;
303 m_phi_ID = cid - 24 * 16 ;
304 m_phi_idx = m_phi_ID % 6 + 24;
305 m_theta_idx = m_phi_ID / 6;
306 } else if (cid < 36 * 16) {
307 m_theta_ID = 7;
308 m_phi_ID = cid - 30 * 16 ;
309 m_phi_idx = m_phi_ID % 6 + 30;
310 m_theta_idx = m_phi_ID / 6;
311 } else if (cid < 42 * 16) {
312 m_theta_ID = 8;
313 m_phi_ID = cid - 36 * 16 ;
314 m_phi_idx = m_phi_ID % 6 + 36;
315 m_theta_idx = m_phi_ID / 6;
316 } else if (cid < 48 * 16) {
317 m_theta_ID = 9;
318 m_phi_ID = cid - 42 * 16 ;
319 m_phi_idx = m_phi_ID % 6 + 42;
320 m_theta_idx = m_phi_ID / 6;
321 } else if (cid < 54 * 16) {
322 m_theta_ID = 10;
323 m_phi_ID = cid - 48 * 16 ;
324 m_phi_idx = m_phi_ID % 6 + 48;
325 m_theta_idx = m_phi_ID / 6;
326 } else if (cid < 63 * 16) {
327 m_theta_ID = 11;
328 m_phi_ID = cid - 54 * 16 ;
329 m_phi_idx = m_phi_ID % 9 + 54;
330 m_theta_idx = m_phi_ID / 9;
331 } else if (cid < 72 * 16) {
332 m_theta_ID = 12;
333 m_phi_ID = cid - 63 * 16 ;
334 m_phi_idx = m_phi_ID % 9 + 63;
335 m_theta_idx = m_phi_ID / 9;
336 } else if (cid < ECLElementNumbers::c_NCrystalsForwardBarrel) {//Barrel start
337 m_phi_ID = (cid - ECLElementNumbers::c_NCrystalsForward) % 144;
338 m_theta_ID = (cid - ECLElementNumbers::c_NCrystalsForward) / 144 + 13;
339 m_theta_idx = (cid - ECLElementNumbers::c_NCrystalsForward) / 144;
340 m_phi_idx = m_phi_ID;
341 } else if (cid < ECLElementNumbers::c_NCrystalsForwardBarrel + 9 * 16) { //Backward start
342 m_theta_ID = 59;
343 m_phi_ID = cid - ECLElementNumbers::c_NCrystalsForwardBarrel ;
344 m_phi_idx = m_phi_ID % 9 + 72;
345 m_theta_idx = m_phi_ID / 9;
346 } else if (cid < ECLElementNumbers::c_NCrystalsForwardBarrel + 18 * 16) {
347 m_theta_ID = 60;
348 m_phi_ID = cid - ECLElementNumbers::c_NCrystalsForwardBarrel - 9 * 16 ;
349 m_phi_idx = m_phi_ID % 9 + 81;
350 m_theta_idx = m_phi_ID / 9;
351 } else if (cid < ECLElementNumbers::c_NCrystalsForwardBarrel + 24 * 16) {
352 m_theta_ID = 61;
353 m_phi_ID = cid - ECLElementNumbers::c_NCrystalsForwardBarrel - 18 * 16 ;
354 m_phi_idx = m_phi_ID % 6 + 90;
355 m_theta_idx = m_phi_ID / 6;
356 } else if (cid < ECLElementNumbers::c_NCrystalsForwardBarrel + 30 * 16) {
357 m_theta_ID = 62;
358 m_phi_ID = cid - ECLElementNumbers::c_NCrystalsForwardBarrel - 24 * 16 ;
359 m_phi_idx = m_phi_ID % 6 + 96;
360 m_theta_idx = m_phi_ID / 6;
361 } else if (cid < ECLElementNumbers::c_NCrystalsForwardBarrel + 36 * 16) {
362 m_theta_ID = 63;
363 m_phi_ID = cid - ECLElementNumbers::c_NCrystalsForwardBarrel - 30 * 16 ;
364 m_phi_idx = m_phi_ID % 6 + 102;
365 m_theta_idx = m_phi_ID / 6;
366 } else if (cid < ECLElementNumbers::c_NCrystalsForwardBarrel + 42 * 16) {
367 m_theta_ID = 64;
368 m_phi_ID = cid - ECLElementNumbers::c_NCrystalsForwardBarrel - 36 * 16 ;
369 m_phi_idx = m_phi_ID % 6 + 108;
370 m_theta_idx = m_phi_ID / 6;
371 } else if (cid < ECLElementNumbers::c_NCrystalsForwardBarrel + 48 * 16) {
372 m_theta_ID = 65;
373 m_phi_ID = cid - ECLElementNumbers::c_NCrystalsForwardBarrel - 42 * 16 ;
374 m_phi_idx = m_phi_ID % 6 + 114;
375 m_theta_idx = m_phi_ID / 6;
376 } else if (cid < ECLElementNumbers::c_NCrystalsForwardBarrel + 52 * 16) {
377 m_theta_ID = 66;
378 m_phi_ID = cid - ECLElementNumbers::c_NCrystalsForwardBarrel - 48 * 16 ;
379 m_phi_idx = m_phi_ID % 4 + 120;
380 m_theta_idx = m_phi_ID / 4;
381 } else if (cid < ECLElementNumbers::c_NCrystalsForwardBarrel + 56 * 16) {
382 m_theta_ID = 67;
383 m_phi_ID = cid - ECLElementNumbers::c_NCrystalsForwardBarrel - 52 * 16 ;
384 m_phi_idx = m_phi_ID % 4 + 124;
385 m_theta_idx = m_phi_ID / 4;
386 } else if (cid < ECLElementNumbers::c_NCrystalsForwardBarrel + 60 * 16) {
387 m_theta_ID = 68;
388 m_phi_ID = cid - ECLElementNumbers::c_NCrystalsForwardBarrel - 56 * 16 ;
389 m_phi_idx = m_phi_ID % 4 + 128;
390 m_theta_idx = m_phi_ID / 4;
391 } else {
392 B2ERROR("ECL ECLCrystalData Mapping " << cid << ". Out of range.");
393 }
394}
Belle2::ECLCrystalData::Eval
void Eval()
Set values for the crystal.
Definition: ECLCrystalData.cc:152
Belle2::ECLCrystalData::theta
static const double theta[69]
Theta ID to theta(deg)
Definition: ECLCrystalData.h:90
Belle2::ECLCrystalData::m_PosY
float m_PosY
Crystal Position - Y.
Definition: ECLCrystalData.h:57
Belle2::ECLCrystalData::m_mass
float m_mass
Crystal Mass.
Definition: ECLCrystalData.h:43
Belle2::ECLCrystalData::m_PosZ
float m_PosZ
Crystal Position - Z.
Definition: ECLCrystalData.h:59
Belle2::ECLCrystalData::EvalPhi
double EvalPhi()
Set phi value of crystal.
Definition: ECLCrystalData.cc:223
Belle2::ECLCrystalData::EvalR
double EvalR()
Set R value of crystal.
Definition: ECLCrystalData.cc:237
Belle2::ECLCrystalData::EvalMass
float EvalMass()
Evaluate mass of crystal.
Definition: ECLCrystalData.cc:181
Belle2::ECLCrystalData::EvalZ
double EvalZ()
Set Z value of crystal.
Definition: ECLCrystalData.cc:228
Belle2::ECLCrystalData::GetCrystalIndex
int GetCrystalIndex()
return crystal index
Definition: ECLCrystalData.cc:191
Belle2::ECLCrystalData::m_volume
float m_volume
Crystal Volume.
Definition: ECLCrystalData.h:45
Belle2::ECLCrystalData::EvalY
double EvalY()
Set Y value of crystal.
Definition: ECLCrystalData.cc:174
Belle2::ECLCrystalData::CrystalMassEndcap
static const float CrystalMassEndcap[132]
Barrel crystal mass.
Definition: ECLCrystalData.h:84
Belle2::ECLCrystalData::m_cell_ID
int m_cell_ID
Cell ID of crystal.
Definition: ECLCrystalData.h:29
Belle2::ECLCrystalData::sumPrevious
static const int sumPrevious[24]
Number of crystals in preceding theta rings (endcaps only)
Definition: ECLCrystalData.h:81
Belle2::ECLCrystalData::IsEndCap
bool IsEndCap()
True if crystal is in endcaps.
Definition: ECLCrystalData.cc:207
Belle2::ECLCrystalData::Ring
static const int Ring[69]
ECL has 16-fold symmetry in phi.
Definition: ECLCrystalData.h:79
Belle2::ECLCrystalData::m_PosR
float m_PosR
Crystal Position - R.
Definition: ECLCrystalData.h:51
Belle2::ECLCrystalData::m_PosP
float m_PosP
Crystal Position - P.
Definition: ECLCrystalData.h:53
Belle2::ECLCrystalData::EndcapRadius
static const double EndcapRadius[132]
EndCapCrystalID to radius (spherical)
Definition: ECLCrystalData.h:93
Belle2::ECLCrystalData::Mapping
void Mapping(int cid)
set theta and phi value of crystal
Definition: ECLCrystalData.cc:266
Belle2::ECLCrystalData::ECLCrystalData
ECLCrystalData()
Blank Constructor.
Definition: ECLCrystalData.cc:109
Belle2::ECLCrystalData::m_phi_ID
int m_phi_ID
phi ID of crystal
Definition: ECLCrystalData.h:31
Belle2::ECLCrystalData::CrystalMassBarrel
static const float CrystalMassBarrel[46]
Endcap crystal mass.
Definition: ECLCrystalData.h:86
Belle2::ECLCrystalData::m_theta_ID
int m_theta_ID
theta ID of crystal
Definition: ECLCrystalData.h:33
Belle2::ECLCrystalData::IsBarrel
bool IsBarrel()
True if crystal is in barrel.
Definition: ECLCrystalData.cc:212
Belle2::ECLCrystalData::EvalX
double EvalX()
Set X value of crystal.
Definition: ECLCrystalData.cc:169
Belle2::ECLCrystalData::m_PosX
float m_PosX
Crystal Position - X.
Definition: ECLCrystalData.h:55
Belle2::ECLCrystalData::GetCellID
int GetCellID()
return cell ID
Definition: ECLCrystalData.cc:218
Belle2::ECLCrystalData::BarrelZ
static const double BarrelZ[46]
BarrelCrystalID to z-coordinate.
Definition: ECLCrystalData.h:96
Belle2::ECLElementNumbers::c_NCrystalsForwardBarrel
const int c_NCrystalsForwardBarrel
Number of crystals in the forward and barrel ECL.
Definition: ECLElementNumbers.h:35
Belle2::ECLElementNumbers::c_NCrystalsForward
const int c_NCrystalsForward
Number of crystals in the forward ECL.
Definition: ECLElementNumbers.h:26
Belle2
Abstract base class for different kinds of events.
Definition: MillepedeAlgorithm.h:17