ATLAS Offline Software
Loading...
Searching...
No Matches
GeoCaloCalibHit.icc
Go to the documentation of this file.
1/*
2 Copyright (C) 2002-2021 CERN for the benefit of the ATLAS collaboration
3*/
4
5//----------------------------------------------------------//
6// //
7// Adaptor for CaloCalibHits. //
8// //
9// Joe Boudreau Apr 14, 2005 //
10// Mikhail Leltchouk Apr 30, 2005 //
11// //
12// //
13//----------------------------------------------------------//
14
15#include "CaloIdentifier/CaloIdManager.h"
16#include "CaloDetDescr/CaloDetDescrElement.h"
17#include "CaloSimEvent/CaloCalibrationHit.h"
18#include "CaloIdentifier/CaloDM_ID.h"
19#include "CaloDetDescr/CaloDetDescrManager.h"
20
21#include "CLHEP/Units/SystemOfUnits.h"
22
23#undef DEBUG_GEOHIT
24
25inline GeoCaloCalibHit::GeoCaloCalibHit(const CaloCalibrationHit & hit
26 , const std::string & collectionName
27 , const CaloDetDescrManager* caloMgr)
28 : m_hit(&hit)
29 , m_ddManager(caloMgr)
30{
31 init(collectionName);
32}
33
34inline GeoCaloCalibHit::operator bool () const {
35 return !m_failed;
36}
37
38inline double GeoCaloCalibHit::energyEM() const
39{
40 return m_hit->energyEM();
41}
42
43inline double GeoCaloCalibHit::energyNonEM() const
44{
45 return m_hit->energyNonEM();
46}
47
48inline double GeoCaloCalibHit::energyInvisible() const
49{
50 return m_hit->energyInvisible();
51}
52
53inline double GeoCaloCalibHit::energyEscaped() const
54{
55 return m_hit->energyEscaped();
56}
57
58
59inline double GeoCaloCalibHit::energyTotal() const
60{
61 return m_hit->energyTotal();
62}
63
64
65inline double GeoCaloCalibHit::phiMin() const {
66 if (m_ddElement) { // Active and Inactive hits
67 double phi = m_ddElement->phi();
68 double dphi = m_ddElement->dphi();
69 return phi-dphi/2.0;
70 }
71 const CaloDM_ID* dmid = m_ddManager->getCalo_Mgr()->getDM_ID();
72 int phiIndex=dmid->phi(m_hit->cellID()); // Dead material hits
73 int n_phi=32;
74 if(abs(m_subdet) == 4 && m_type==2 && (m_region == 2 || m_region == 3) ) {
75 n_phi=16;
76 return (phiIndex-0.5)*M_PI/n_phi;
77 }
78 return phiIndex*M_PI/n_phi;
79}
80
81inline double GeoCaloCalibHit::phiMax() const {
82 if (m_ddElement) { // Active and Inactive hits
83 double phi = m_ddElement->phi();
84 double dphi = m_ddElement->dphi();
85 return phi+dphi/2.0;
86 }
87 const CaloDM_ID* dmid = m_ddManager->getCalo_Mgr()->getDM_ID();
88 int phiIndex=dmid->phi(m_hit->cellID()); // Dead material hits
89 int n_phi=32;
90 if(abs(m_subdet) == 4 && m_type==2 && (m_region == 2 || m_region == 3) ) {
91 n_phi=16;
92 return (phiIndex+0.5)*M_PI/n_phi;
93 }
94 return (phiIndex+1.)*M_PI/n_phi;
95}
96
97inline double GeoCaloCalibHit::phi() const {
98 if (m_ddElement) { // Active and Inactive hits
99 double phi = m_ddElement->phi();
100 return phi;
101 }
102 const CaloDM_ID* dmid = m_ddManager->getCalo_Mgr()->getDM_ID();
103 int phiIndex=dmid->phi(m_hit->cellID()); // Dead material hits
104 int n_phi=32;
105 if(abs(m_subdet) == 4 && m_type==2 && (m_region == 2 || m_region == 3) ) {
106 n_phi=16;
107 return (phiIndex)*M_PI/n_phi; // phiBin=0 means between hec module 0 and hec module 1, i.e. 11.25deg
108 }
109 return (phiIndex+0.5)*M_PI/n_phi;
110}
111
112
113inline double GeoCaloCalibHit::etaMin() const {
114 if (m_ddElement) { // Active and Inactive hits
115 double eta = m_ddElement->eta();
116 double deta = m_ddElement->deta();
117 return eta-deta/2.0;
118 }
119 double halfEta = 0.05; // Dead material hits
120 if ((m_type==1 && m_sampling==3 && m_region==1) ||
121 (m_type==2 && m_region==3)) {
122 halfEta = 0.1;
123 }
124 return eta()-halfEta;
125}
126
127inline double GeoCaloCalibHit::etaMax() const {
128 if (m_ddElement) { // Active and Inactive hits
129 double eta = m_ddElement->eta();
130 double deta = m_ddElement->deta();
131 return eta+deta/2.0;
132 }
133 double halfEta = 0.05; // Dead material hits
134 if ((m_type==1 && m_sampling==3 && m_region==1) ||
135 (m_type==2 && m_region==3)) {
136 halfEta = 0.1;
137 }
138 return eta()+halfEta;
139}
140
141inline double GeoCaloCalibHit::eta() const {
142 if (m_ddElement) { // Active and Inactive hits
143 double eta = m_ddElement->eta();
144 return eta;
145 }
146 const CaloDM_ID* dmid = m_ddManager->getCalo_Mgr()->getDM_ID();
147 int etaIndex=dmid->eta(m_hit->cellID()); // Dead material hits
148 double dEta = 0.1;
149 if ((m_type==1 && m_sampling==3 && m_region==1) ||
150 (m_type==2 && m_region==3)) {
151 dEta = 0.2;
152 }
153 double etaCenter = m_etaOffset + (etaIndex+0.5)*dEta;
154
155 bool isNegative = dmid->pos_neg_z(m_hit->cellID())<0;
156 if (isNegative) {
157 etaCenter = -etaCenter;
158 }
159 return etaCenter;
160}
161
162inline std::string GeoCaloCalibHit::hitType() const {
163 if (m_hitType==ACTIVE) return "LArCalibrationHitActive";
164 if (m_hitType==INACTIVE) return "LArCalibrationHitInactive";
165 if (m_hitType==DEAD) return "LArCalibrationHitDeadMaterial";
166 if (m_hitType==TILEAI) return "TileCalibrationCellHitCnt"; // Old format
167 if (m_hitType==TILEDM) return "TileCalibrationDMHitCnt";// Old format
168 if (m_hitType==TILEACTIVE) return "TileCalibHitActiveCell"; // New format
169 if (m_hitType==TILEINACTIVE) return "TileCalibHitInactiveCell";// New format
170 if (m_hitType==TILEDEAD) return "TileCalibHitDeadMaterial";// New format
171 if (m_hitType==UNRECOGNIZED) return "UNRECOGNIZED hit collection";
172 return "UNRECOGNIZED hit collection";
173}
174
175inline void GeoCaloCalibHit::init(const std::string & collectionName) {
176
177 if (collectionName=="LArCalibrationHitActive") {
178 m_hitType=ACTIVE;
179 }
180 else if (collectionName=="LArCalibrationHitInactive") {
181 m_hitType=INACTIVE;
182 }
183 else if (collectionName=="LArCalibrationHitDeadMaterial") {
184 m_hitType=DEAD;
185 }
186 else if (collectionName=="TileCalibrationCellHitCnt") {
187 m_hitType=TILEAI;
188 }
189 else if (collectionName=="TileCalibrationDMHitCnt") {
190 m_hitType=TILEDM;
191 }
192 else if (collectionName=="TileCalibHitDeadMaterial") {
193 m_hitType=TILEDEAD;
194 }
195 else if (collectionName=="TileCalibHitActiveCell") {
196 m_hitType=TILEACTIVE;
197 }
198 else if (collectionName=="TileCalibHitInactiveCell") {
199 m_hitType=TILEINACTIVE;
200 }
201 else {
202 m_hitType=UNRECOGNIZED;
203 m_failed=true;
204 return;
205 }
206
207 m_ddElement = getDetDescrElement();
208
209 //bool inactive = collectionName=="LArCalibrationHitInactive";
210
211 // Default values of m_etaOffset, m_distance, m_isBarrel(hit orientation)
212 m_etaOffset = 0.;
213 m_distance=230*CLHEP::mm;
214 m_isBarrel=true;
215
216 // if (abs(m_subdet)>= 4) { // Dead hits: LAr: m_subdet= +-4, Tile: m_subdet= +-5)
217 if (m_hitType==DEAD || m_hitType==TILEDM || m_hitType==TILEDEAD) { // Dead hits: LAr: m_subdet= +-4, Tile: m_subdet= +-5)
218 // from Calorimeter/ CaloIdentifier/ CaloIdentifier/ CaloDM_ID.h
219
220 const CaloDM_ID* dmid = m_ddManager->getCalo_Mgr()->getDM_ID();
221 m_subdet = dmid->pos_neg_z(m_hit->cellID());
222 m_type = dmid->dmat(m_hit->cellID());
223 m_sampling = dmid->sampling(m_hit->cellID());
224 m_region = dmid->region(m_hit->cellID());
225 m_etaBin = dmid->eta(m_hit->cellID());
226 m_phiBin = dmid->phi(m_hit->cellID());
227
228 } else { // Active and Inactive hit ID from Calorimeter/ CaloIdentifier/ CaloIdentifier/ CaloCell_ID.h
229
230 const CaloCell_ID* cellid = m_ddManager->getCaloCell_ID();
231 m_subdet = 1 + cellid->sub_calo(m_hit->cellID());
232 m_type = cellid->pos_neg(m_hit->cellID());
233 m_sampling = cellid->calo_sample(m_hit->cellID());
234 m_region = cellid->region(m_hit->cellID());
235 m_etaBin = cellid->eta(m_hit->cellID());
236 m_phiBin = cellid->phi(m_hit->cellID());
237
238 }
239
240 // Values of m_etaOffset, m_distance, m_isBarrel(hit orientation) which are different from default
241 // To help navigate through statements we will use in comments below
242 // labeling key=abs(subdet)*1000+type*100+sampling*10+region;
243
244 // for LAr Dead hits:
245 if (abs(m_subdet) == 4) { // subdet = 4 LAr.
246 switch (m_type) {
247 case 1: // type = 1
248 switch (m_sampling) {
249 case 0: // sampling = 0 Inner Detector
250 switch (m_region) {
251 case 0: // 4100
252 m_distance= 250*CLHEP::mm;
253 return;
254 case 1: // 4101
255 m_distance= 270*CLHEP::mm; // largest radius for this zone
256 return;
257 case 2: // 4102
258 m_distance= 405*CLHEP::mm;
259 return;
260 case 3: // 4103
261 m_distance= 595*CLHEP::mm;
262 return;
263 case 4: // 4104
264 m_distance= 815*CLHEP::mm;
265 return;
266 case 5: // 4105
267 m_distance=1065*CLHEP::mm;
268 return;
269 default:
270#ifdef DEBUG_GEOHIT
271 std::cout <<"Warning 0, GeoCaloCalibHit with (subdet,type,sampling,region,eta,phi) = "
272 <<m_subdet<<", "<<m_type<<", "<<m_sampling<<", "<<m_region
273 <<m_etaBin<<", "<<m_phiBin
274 << " not decoded" <<std::endl;
275#endif
276 m_failed=true;
277 return;
278 }
279 break;
280 case 1: // sampling = 1 - in front of and in active LAr calorimeters
281 switch (m_region) {
282 case 0: // 4110
283 m_distance=1250*CLHEP::mm;
284 return;
285 case 1: // 4111
286 m_distance=1375*CLHEP::mm;
287 return;
288 case 2: // 4112
289 m_distance=1410*CLHEP::mm;
290 return;
291 case 3: // 4113
292 m_distance=1490*CLHEP::mm;
293 return;
294 case 4: // 4114
295 m_etaOffset = 1.5;
296 m_distance=3490*CLHEP::mm;
297 m_isBarrel=false;
298 return;
299 case 5: // 4115
300 m_etaOffset = 1.5;
301 m_distance=3511*CLHEP::mm;
302 m_isBarrel=false;
303 return;
304 case 6: // 4116
305 m_etaOffset = 1.3;
306 m_distance=3640*CLHEP::mm;
307 m_isBarrel=false;
308 return;
309 case 7: // 4117
310 m_etaOffset = 3.2;
311 m_distance=3700*CLHEP::mm;
312 m_isBarrel=false;
313 return;
314 default:
315#ifdef DEBUG_GEOHIT
316 std::cout <<"Warning 1, GeoCaloCalibHit with (subdet,type,sampling,region,eta,phi) = "
317 <<m_subdet<<", "<<m_type<<", "<<m_sampling<<", "<<m_region
318 <<m_etaBin<<", "<<m_phiBin
319 << " not decoded" <<std::endl;
320#endif
321 m_failed=true;
322 return;
323 }
324 break;
325 case 2: // sampling = 2 - dead materials between active calorimeters
326 switch (m_region) {
327 case 0: // 4120
328 m_distance=2170*CLHEP::mm;
329 return;
330 case 1: // 4121
331 m_distance=2990*CLHEP::mm; // radius of lower part - feedthrough
332 return;
333 case 2: // 4122
334 m_etaOffset = 1.0;
335 m_distance=3300*CLHEP::mm; // z
336 m_isBarrel=false; // partly around EMB but in front of Extended Barrel of Tile
337 return;
338 case 3: // 4123
339 m_etaOffset = 1.5;
340 m_distance=4250*CLHEP::mm;
341 m_isBarrel=false;
342 return;
343 case 4: // 4124
344 m_etaOffset = 1.5;
345 m_distance=5110*CLHEP::mm;
346 m_isBarrel=false;
347 return;
348 case 5: // 4125
349 m_etaOffset = 1.0;
350 m_distance=2100*CLHEP::mm; // R
351 return;
352 default:
353#ifdef DEBUG_GEOHIT
354 std::cout <<"Warning 2, GeoCaloCalibHit with (subdet,type,sampling,region,eta,phi) = "
355 <<m_subdet<<", "<<m_type<<", "<<m_sampling<<", "<<m_region
356 <<m_etaBin<<", "<<m_phiBin
357 << " not decoded" <<std::endl;
358#endif
359 m_failed=true;
360 return;
361 }
362 break;
363 case 3: // sampling = 3
364 // Eta granularity in region 1 sampling 3 is twice that of previous cases.
365 switch(m_region) {
366 case 0: // 4130
367 m_etaOffset = 1.7;
368 m_distance=6200*CLHEP::mm;
369 m_isBarrel=false;
370 return;
371 case 1: // 4131
372 m_etaOffset = 5.0;
373 m_distance=6300*CLHEP::mm;
374 m_isBarrel=false;
375 return;
376 case 2: // 4132
377 m_etaOffset = 8.0;
378 m_distance=3500*CLHEP::mm; //or can be attached to end of InDet cavity
379 m_isBarrel=false;
380 return;
381 default:
382#ifdef DEBUG_GEOHIT
383 std::cout <<"Warning 4, GeoCaloCalibHit with (subdet,type,sampling,region,eta,phi) = "
384 <<m_subdet<<", "<<m_type<<", "<<m_sampling<<", "<<m_region
385 <<m_etaBin<<", "<<m_phiBin
386 << " not decoded" <<std::endl;
387#endif
388 m_failed=true;
389 return;
390 }
391 default:
392#ifdef DEBUG_GEOHIT
393 std::cout <<"Warning 5, GeoCaloCalibHit with (subdet,type,sampling,region,eta,phi) = "
394 <<m_subdet<<", "<<m_type<<", "<<m_sampling<<", "<<m_region
395 <<m_etaBin<<", "<<m_phiBin
396 << " not decoded" <<std::endl;
397#endif
398 m_failed=true;
399 return;
400 }
401 case 2: // type = 2
402 switch (m_region) {
403 case 0: // region = 0 crack between EMB halfs
404 switch (m_sampling) {
405 case 0: // 4200
406 // distances should be checked with experts
407 m_distance= 1550*CLHEP::mm;
408 return;
409 case 1: // 4210
410 m_distance= 1600*CLHEP::mm;
411 return;
412 case 2: // 4220
413 m_distance= 1750*CLHEP::mm;
414 return;
415 case 3: // 4230
416 m_distance= 1900*CLHEP::mm;
417 return;
418 default:
419#ifdef DEBUG_GEOHIT
420 std::cout <<"Warning 6, GeoCaloCalibHit with (subdet,type,sampling,region,eta,phi) = "
421 <<m_subdet<<", "<<m_type<<", "<<m_sampling<<", "<<m_region
422 <<m_etaBin<<", "<<m_phiBin
423 << " not decoded" <<std::endl;
424#endif
425 m_failed=true;
426 return;
427 }
428 case 1: // region = 1 dead materials between inner and outer EMEC wheels
429 m_isBarrel=false;
430 m_etaOffset = 2.45;
431 switch (m_sampling) {
432 case 0: // 4201
433 // distances should be checked with experts
434 m_distance=3700*CLHEP::mm;
435 return;
436 case 1: // 4211
437 m_distance=3800*CLHEP::mm;
438 return;
439 case 2: // 4221
440 m_distance=3950*CLHEP::mm;
441 return;
442 case 3: // 4231
443 m_distance=4150*CLHEP::mm;
444 return;
445 default:
446#ifdef DEBUG_GEOHIT
447 std::cout <<"Warning 7, GeoCaloCalibHit with (subdet,type,sampling,region,eta,phi) = "
448 <<m_subdet<<", "<<m_type<<", "<<m_sampling<<", "<<m_region
449 <<m_etaBin<<", "<<m_phiBin
450 << " not decoded" <<std::endl;
451#endif
452 m_failed=true;
453 return;
454 }
455 case 2: // region = 2 - azimuthal cracks between HEC modules (wedges) - outer
456 m_isBarrel=false;
457 m_etaOffset = 1.5;
458 switch (m_sampling) {
459 case 0: // 4202
460 // distances should be checked with experts
461 m_distance=4417*CLHEP::mm;
462 return;
463 case 1: // 4212
464 m_distance=4825*CLHEP::mm;
465 return;
466 case 2: // 4222
467 m_distance=5393*CLHEP::mm;
468 return;
469 case 3: // 4232
470 m_distance=5861*CLHEP::mm;
471 return;
472 default:
473#ifdef DEBUG_GEOHIT
474 std::cout <<"Warning 8, GeoCaloCalibHit with (subdet,type,sampling,region,eta,phi) = "
475 <<m_subdet<<", "<<m_type<<", "<<m_sampling<<", "<<m_region
476 <<m_etaBin<<", "<<m_phiBin
477 << " not decoded" <<std::endl;
478#endif
479 m_failed=true;
480 return;
481 }
482 case 3: // region = 3 - azimuthal cracks between HEC modules (wedges) - inner
483 m_isBarrel=false;
484 m_etaOffset = 2.5;
485 switch (m_sampling) {
486 case 0: // 4203
487 // distances should be checked with experts
488 m_distance=4417*CLHEP::mm;
489 return;
490 case 1: // 4213
491 m_distance=4825*CLHEP::mm;
492 return;
493 case 2: // 4223
494 m_distance=5393*CLHEP::mm;
495 return;
496 case 3: // 4233
497 m_distance=5861*CLHEP::mm;
498 return;
499 default:
500#ifdef DEBUG_GEOHIT
501 std::cout <<"Warning 9, GeoCaloCalibHit with (subdet,type,sampling,region,eta,phi) = "
502 <<m_subdet<<", "<<m_type<<", "<<m_sampling<<", "<<m_region
503 <<m_etaBin<<", "<<m_phiBin
504 << " not decoded" <<std::endl;
505#endif
506 m_failed=true;
507 return;
508 }
509 case 4: // region = 4 - dead materials between HEC and FCal
510 m_isBarrel=false;
511 m_etaOffset = 3.0;
512 switch (m_sampling) {
513 case 0: // 4204
514 // distances should be checked with experts
515 m_distance=4600*CLHEP::mm;
516 return;
517 case 1: // 4214
518 m_distance=4920*CLHEP::mm;
519 return;
520 case 2: // 4224
521 m_distance=5400*CLHEP::mm;
522 return;
523 case 3: // 4234
524 m_distance=5820*CLHEP::mm;
525 return;
526 default:
527#ifdef DEBUG_GEOHIT
528 std::cout <<"Warning 10, GeoCaloCalibHit with (subdet,type,sampling,region,eta,phi) = "
529 <<m_subdet<<", "<<m_type<<", "<<m_sampling<<", "<<m_region
530 <<m_etaBin<<", "<<m_phiBin
531 << " not decoded" <<std::endl;
532#endif
533 m_failed=true;
534 return;
535 }
536 case 5: // region = 4 - dead materials between HEC and FCal
537 m_isBarrel=false;
538 m_etaOffset = 3.0;
539 switch (m_sampling) {
540 case 0: // 4205
541 // distances should be checked with experts
542 m_distance=4340*CLHEP::mm;
543 return;
544 case 1: // 4215
545 m_distance=4668*CLHEP::mm;
546 return;
547 case 2: // 4225
548 m_distance=5138*CLHEP::mm;
549 return;
550 case 3: // 4235
551 m_distance=5602*CLHEP::mm;
552 return;
553 default:
554#ifdef DEBUG_GEOHIT
555 std::cout <<"Warning 11, GeoCaloCalibHit with (subdet,type,sampling,region,eta,phi) = "
556 <<m_subdet<<", "<<m_type<<", "<<m_sampling<<", "<<m_region
557 <<m_etaBin<<", "<<m_phiBin
558 << " not decoded" <<std::endl;
559#endif
560 m_failed=true;
561 return;
562 }
563
564 }
565 }
566 }
567 // --------------------------- Tile Dead hits --------------------------------------
568 // Values of m_etaOffset, m_distance, m_isBarrel(hit orientation) which are different from default
569 else if (abs(m_subdet) == 5) { //for Tile Dead hits
570
571 switch (m_type) {
572 case 1: // type = 1
573 switch (m_sampling) {
574 case 0: // sampling = 0, front plate of the module (1 CLHEP::cm Fe along R)
575 switch (m_region) {
576 case 0: // 5100, region = 0 - barrel 0.0 < eta < 1.0
577 m_distance= 2295*CLHEP::mm;
578 return;
579 case 1: // 5101, region = 1 - ext.barrel 1.1 < eta < 1.6
580 m_distance= 2295*CLHEP::mm;
581 m_etaOffset = 1.1;
582 return;
583 default:
584#ifdef DEBUG_GEOHIT
585 std::cout <<"Warning 12, GeoCaloCalibHit with (subdet,type,sampling,region,eta,phi) = "
586 <<m_subdet<<", "<<m_type<<", "<<m_sampling<<", "<<m_region
587 <<m_etaBin<<", "<<m_phiBin
588 << " not decoded" <<std::endl;
589#endif
590 m_failed=true;
591 return;
592 }
593 break;
594 case 1: // sampling = 1, end plates of the module ( Fe at fixed Z )
595 m_isBarrel=false;
596 switch (m_region) {
597 case 0: // 5110 - region = 0 - end of barrel 0.7 < eta < 1.1
598 m_distance= 2810*CLHEP::mm; // z, thickness = 20 mm
599 m_etaOffset = 0.7;
600 return;
601 case 1: // 5111 - region = 1 - smaller Z end of ext barrel 0.9 < eta < 1.2
602 m_distance= 3527*CLHEP::mm; // z, thickness = 15 mm
603 m_etaOffset = 0.9;
604 return;
605 case 2: // 5112 - region = 2 - higher Z end of ext barrel 1.3 < eta < 1.7
606 m_distance= 6100*CLHEP::mm; // z, thickness = 20 mm
607 m_etaOffset = 1.3;
608 return;
609 default:
610#ifdef DEBUG_GEOHIT
611 std::cout <<"Warning 13, GeoCaloCalibHit with (subdet,type,sampling,region,eta,phi) = "
612 <<m_subdet<<", "<<m_type<<", "<<m_sampling<<", "<<m_region
613 <<m_etaBin<<", "<<m_phiBin
614 << " not decoded" <<std::endl;
615#endif
616 m_failed=true;
617 return;
618 }
619 break;
620 case 2: // sampling = 2, so-called "girder" - iron at higher R
621 m_distance= 3900*CLHEP::mm; // R, should be checked with experts
622 switch (m_region) {
623 case 0: // 5120 - region = 0 - barrel 0. < eta < 0.7
624 // nothing to do
625 return;
626 case 1: // 5121 - region = 1 - ext barrel 0.9 < eta < 1.3
627 m_etaOffset = 0.9;
628 return;
629 case 2: // 5122 - region = 2 - gap between barrel and ext barrel 0.7 < eta < 0.9
630 m_etaOffset = 0.7;
631 return;
632 default:
633#ifdef DEBUG_GEOHIT
634 std::cout <<"Warning 14, GeoCaloCalibHit with (subdet,type,sampling,region,eta,phi) = "
635 <<m_subdet<<", "<<m_type<<", "<<m_sampling<<", "<<m_region
636 <<m_etaBin<<", "<<m_phiBin
637 << " not decoded" <<std::endl;
638#endif
639 m_failed=true;
640 return;
641 }
642 break;
643 case 3: // 5130 - sampling = 3, leakage outside Tile calorimeters, 0.0 < eta < 1.7; usually is in LAr DM hit container
644 m_distance = 3980*CLHEP::mm;
645 return;
646 default:
647 m_failed = true;
648 return;
649 } // over samplings
650 case 2: // type == 2, 5200
651 if (m_region !=0 || m_sampling !=0) {
652 m_failed = true;
653 return;
654 }
655 m_etaOffset = 0.0;
656 m_distance = 4000*CLHEP::mm;
657 m_isBarrel = true;
658 return;
659 default:
660#ifdef DEBUG_GEOHIT
661 std::cout <<"Warning 15, GeoCaloCalibHit with (subdet,type,sampling,region,eta,phi) = "
662 <<m_subdet<<", "<<m_type<<", "<<m_sampling<<", "<<m_region
663 <<m_etaBin<<", "<<m_phiBin
664 << " not decoded" <<std::endl;
665#endif
666 m_failed=true;
667 return;
668 } // over types for Tile DM hits
669#ifdef DEBUG_GEOHIT
670 std::cout <<"Warning 16, GeoCaloCalibHit with (subdet,type,sampling,region,eta,phi) = "
671 <<m_subdet<<", "<<m_type<<", "<<m_sampling<<", "<<m_region
672 <<m_etaBin<<", "<<m_phiBin
673 << " not decoded" <<std::endl;
674#endif
675
676 }
677
678 // ------------------------ Active and Inactive hits --------------------------------
679 // Values of m_etaOffset, m_distance are calculated in eta(), distance()
680 // Value of m_isBarrel(hit orientation) which is different from default is defined here
681 else if (m_subdet == 1) { // LAr EM hits
682 if (abs(m_type)>1){ // EMEC
683 m_isBarrel=false;
684 }
685 if (abs(m_type) == 2){ // EMEC Outer Wheel
686 m_sampling = m_sampling - 4; // adjustment to ATL-SOFT-2001-004 definition
687 }
688 if (abs(m_type) == 3){ // EMEC Outer Wheel
689 m_sampling = m_sampling - 6; // adjustment to ATL-SOFT-2001-004 definition
690 }
691 }
692 else if (m_subdet == 2) { // LAr HEC hits
693 m_isBarrel=false;
694 m_sampling = m_sampling - 8; // adjustment to ATL-SOFT-2001-004 definition
695 }
696 else if (m_subdet == 3) { // LAr FCal hits
697 m_isBarrel=false;
698 }
699 else {
700 }
701
702}
703
704inline double GeoCaloCalibHit::distance() const {
705
706 if (m_ddElement) { // Active and Inactive hits
707 return m_isBarrel ? m_ddElement->r() : fabs(m_ddElement->z());
708 }
709 return m_distance; // Dead hits
710}
711
712inline bool GeoCaloCalibHit::isBarrel() const {
713 return m_isBarrel;
714}
715
716inline bool GeoCaloCalibHit::isEndcap() const {
717 return !isBarrel();
718}
719
720inline int GeoCaloCalibHit::subdet() const {
721 return m_subdet;
722}
723
724inline int GeoCaloCalibHit::type() const {
725 return m_type;
726}
727
728inline int GeoCaloCalibHit::sampling() const {
729 return m_sampling;
730}
731
732inline int GeoCaloCalibHit::region() const {
733 return m_region;
734}
735
736inline const CaloDetDescrElement* GeoCaloCalibHit::getDetDescrElement() const
737{
738 return (m_hitType!=DEAD && m_hitType!=TILEDM && m_hitType!=TILEDEAD)
739 ? m_ddManager->get_element(m_hit->cellID())
740 : nullptr;
741}
Generated on for ATLAS Offline Software by doxygen 1.14.0