|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
|
#include "DataFormats/EcalDetId/interface/EcalTrigTowerDetId.h"
#include "FWCore/Utilities/interface/Exception.h"
#include <cassert>
EcalTrigTowerDetId::EcalTrigTowerDetId() {}
EcalTrigTowerDetId::EcalTrigTowerDetId(uint32_t rawid) : DetId(rawid) {}
EcalTrigTowerDetId::EcalTrigTowerDetId(int zside, EcalSubdetector subDet, int i, int j, int mode)
: DetId(Ecal, EcalTriggerTower) {
int tower_i = 0;
int tower_j = 0;
if (mode == SUBDETIJMODE) {
tower_i = i;
tower_j = j;
} else if (mode == SUBDETDCCTTMODE) {
throw cms::Exception("InvalidDetId")
<< "EcalTriggerTowerDetId: Cannot create object. SUBDETDCCTTMODE not yet implemented.";
} else
throw cms::Exception("InvalidDetId") << "EcalTriggerTowerDetId: Cannot create object. Unknown mode for (int, "
"EcalSubdetector, int, int) constructor.";
if (tower_i > MAX_I || tower_i < MIN_I || tower_j > MAX_J || tower_j < MIN_J)
throw cms::Exception("InvalidDetId") << "EcalTriggerTowerDetId: Cannot create object. Indexes out of bounds.";
id_ |= ((zside > 0) ? (0x8000) : (0x0)) | ((subDet == EcalBarrel) ? (0x4000) : (0x0)) | (tower_i << 7) |
(tower_j & 0x7F);
}
EcalTrigTowerDetId::EcalTrigTowerDetId(const DetId& gen) {
if (!gen.null() && (gen.det() != Ecal || gen.subdetId() != EcalTriggerTower)) {
throw cms::Exception("InvalidDetId");
}
id_ = gen.rawId();
}
EcalTrigTowerDetId& EcalTrigTowerDetId::operator=(const DetId& gen) {
if (!gen.null() && (gen.det() != Ecal || gen.subdetId() != EcalTriggerTower)) {
throw cms::Exception("InvalidDetId");
}
id_ = gen.rawId();
return *this;
}
//New SM numbering scheme. Avoids discontinuity in phi crossing \eta=0
int EcalTrigTowerDetId::iDCC() const {
if (subDet() == EcalBarrel) {
//Correction since iphi is uniformized with HB convention
int iphi_simple = iphi() + 2;
if (iphi_simple > 72)
iphi_simple = iphi_simple % 72;
int id = (iphi_simple - 1) / kEBTowersInPhi + 1;
if (zside() < 0)
id += 18;
return id;
} else
throw cms::Exception("MethodNotImplemented") << "EcalTriggerTowerDetId: iDCC not yet implemented";
}
int EcalTrigTowerDetId::iTT() const {
if (subDet() == EcalBarrel) {
int ie = ietaAbs() - 1;
int ip;
int iphi_simple = iphi() + 2;
if (iphi_simple > 72)
iphi_simple = iphi_simple % 72;
if (zside() < 0) {
ip = ((iphi_simple - 1) % kEBTowersInPhi) + 1;
} else {
ip = kEBTowersInPhi - ((iphi_simple - 1) % kEBTowersInPhi);
}
return (ie * kEBTowersInPhi) + ip;
} else
throw cms::Exception("MethodNotImplemented") << "EcalTriggerTowerDetId: iTT not yet implemented";
}
int EcalTrigTowerDetId::iquadrant() const {
if (subDet() == EcalEndcap)
return int((iphi() - 1) / kEETowersInPhiPerQuadrant) + 1;
else
throw cms::Exception("MethodNotApplicable") << "EcalTriggerTowerDetId: iquadrant not applicable";
}
bool EcalTrigTowerDetId::validDetId(int iz, EcalSubdetector sd, int i, int j) {
return (
1 == abs(iz) && 0 < i && 0 < j && kEETowersInPhiPerEndcap >= j &&
((EcalBarrel == sd && kEBTowersInEta >= i) || (EcalEndcap == sd && kEEOuterEta <= i && kEEInnerEta >= i &&
(27 > i || ((0 > iz && 0 == j % 2) || (0 < iz && 1 == j % 2))))));
}
int EcalTrigTowerDetId::hashedIndex() const {
const unsigned int iea(ietaAbs());
const unsigned int iph(iphi());
return (subDet() == EcalBarrel ? (iDCC() - 1) * kEBTowersPerSM + iTT() - 1
: kEBTotalTowers + ((zside() + 1) / 2) * kEETowersPerEndcap +
((iea < 27 ? iea : 27) - kEEOuterEta) * kEETowersInPhiPerEndcap +
(iea < 27 ? iph : // for iphi=27,28 only half TT present, odd for EE-, even EE+
(iea - 27) * kEETowersInPhiPerEndcap / 2 + (iph + 1) / 2) -
1);
}
EcalTrigTowerDetId EcalTrigTowerDetId::detIdFromDenseIndex(uint32_t di) {
const EcalSubdetector sd(di < kEBTotalTowers ? EcalBarrel : EcalEndcap);
const int iz(di < kEBTotalTowers ? (di < kEBHalfTowers ? 1 : -1)
: (di - kEBTotalTowers < kEETowersPerEndcap ? -1 : 1));
int i;
int j;
if (di < kEBTotalTowers) // barrel
{
const unsigned int itt(di % kEBTowersPerSM);
const unsigned int idc(di / kEBTowersPerSM);
j = (idc % 18) * kEBTowersInPhi + ((1 + iz) / 2) * kEBTowersInPhi - iz * (itt % kEBTowersInPhi) + 1 - (1 + iz) / 2 -
2;
if (j < 1)
j += 72;
i = 1 + itt / kEBTowersInPhi;
} else {
const int eonly((di - kEBTotalTowers) % kEETowersPerEndcap);
i = kEEOuterEta + eonly / kEETowersInPhiPerEndcap;
j = 1 + eonly % kEETowersInPhiPerEndcap;
if (27 == i) // last two rings have half of normal phi elementes
{
if (j > kEETowersInPhiPerEndcap / 2) {
++i;
j -= kEETowersInPhiPerEndcap / 2;
}
j = 2 * j;
if (0 < iz)
--j;
}
}
assert(validDetId(iz, sd, i, j));
return EcalTrigTowerDetId(iz, sd, i, j);
}
#include <ostream>
std::ostream& operator<<(std::ostream& s, const EcalTrigTowerDetId& id) {
return s << "(EcalTT subDet " << ((id.subDet() == EcalBarrel) ? ("Barrel") : ("Endcap")) << " iz "
<< ((id.zside() > 0) ? ("+ ") : ("- ")) << " ieta " << id.ietaAbs() << " iphi " << id.iphi() << ')';
}
|