|
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
|
#include "DataFormats/ForwardDetId/interface/HGCHEDetId.h"
#include "FWCore/Utilities/interface/Exception.h"
#include <ostream>
const HGCHEDetId HGCHEDetId::Undefined(ForwardEmpty, 0, 0, 0, 0, 0);
HGCHEDetId::HGCHEDetId() : DetId() {}
HGCHEDetId::HGCHEDetId(uint32_t rawid) : DetId(rawid) {}
HGCHEDetId::HGCHEDetId(ForwardSubdetector subdet, int zp, int lay, int sec, int subsec, int cell)
: DetId(Forward, subdet) {
id_ |= ((cell & kHGCHECellMask) << kHGCHECellOffset);
id_ |= ((sec & kHGCHESectorMask) << kHGCHESectorOffset);
if (subsec < 0)
subsec = 0;
id_ |= ((subsec & kHGCHESubSectorMask) << kHGCHESubSectorOffset);
id_ |= ((lay & kHGCHELayerMask) << kHGCHELayerOffset);
if (zp > 0)
id_ |= ((zp & kHGCHEZsideMask) << kHGCHEZsideOffset);
}
HGCHEDetId::HGCHEDetId(const DetId& gen) {
if (!gen.null()) {
ForwardSubdetector subdet = (ForwardSubdetector(gen.subdetId()));
if ((gen.det() != Forward) || (subdet != HGCHEF && subdet != HGCHEB)) {
throw cms::Exception("Invalid DetId") << "Cannot initialize HGCHEDetId from " << std::hex << gen.rawId()
<< std::dec << " Det|SubDet " << gen.det() << "|" << subdet;
}
}
id_ = gen.rawId();
}
HGCHEDetId& HGCHEDetId::operator=(const DetId& gen) {
if (!gen.null()) {
ForwardSubdetector subdet = (ForwardSubdetector(gen.subdetId()));
if ((gen.det() != Forward) || (subdet != HGCHEF && subdet != HGCHEB)) {
throw cms::Exception("Invalid DetId") << "Cannot assign HGCHEDetId from " << std::hex << gen.rawId() << std::dec
<< " Det|SubDet " << gen.det() << "|" << subdet;
}
}
id_ = gen.rawId();
return (*this);
}
HGCHEDetId HGCHEDetId::geometryCell() const {
int sub = ((subdet() == HGCHEF) ? 0 : ((id_ >> kHGCHESubSectorOffset) & kHGCHESubSectorMask));
return HGCHEDetId(subdet(), zside(), layer(), sector(), sub, 0);
}
std::ostream& operator<<(std::ostream& s, const HGCHEDetId& id) {
if (id.subdet() == HGCHEF || id.subdet() == HGCHEB) {
return s << "isHE=" << id.isHE() << " zpos=" << id.zside() << " layer=" << id.layer()
<< " phi subSector=" << id.subsector() << " sector=" << id.sector() << " cell=" << id.cell();
} else {
return s << std::hex << id.rawId() << std::dec;
}
}
|