|
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
|
#ifndef ECALDETID_ESDETID_H
#define ECALDETID_ESDETID_H
#include <iosfwd>
#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/EcalDetId/interface/EcalSubdetector.h"
#include "FWCore/Utilities/interface/Exception.h"
/** \class ESDetId
Det id for a preshower (endcap) strip
*/
class ESDetId : public DetId {
public:
enum { Subdet = EcalPreshower };
/** Constructor of a null id */
ESDetId() {}
/** Constructor from a raw value */
ESDetId(uint32_t rawid) : DetId(rawid) {}
/// constructor from strip, ix, iy, plane, and iz
ESDetId(int strip, int ixs, int iys, int plane, int iz, bool doverify = false) : DetId(Ecal, EcalPreshower) {
id_ |= (strip & 0x3F) | ((ixs & 0x3F) << 6) | ((iys & 0x3F) << 12) | (((plane - 1) & 0x1) << 18) |
((iz > 0) ? (1 << 19) : (0));
if (doverify)
verify(strip, ixs, iys, plane, iz);
}
/** constructor from a generic DetId */
ESDetId(const DetId& id);
/** assignment from a generic DetId */
ESDetId& operator=(const DetId& id);
/// get the subdetector
EcalSubdetector subdet() const { return EcalSubdetector(subdetId()); }
/** get the zside */
int zside() const { return (id_ & 0x80000) ? (1) : (-1); }
/** get the plane */
int plane() const { return ((id_ >> 18) & 0x1) + 1; }
/** get the sensor ix */
int six() const { return (id_ >> 6) & 0x3F; }
/** get the sensor iy */
int siy() const { return (id_ >> 12) & 0x3F; }
/** get the strip */
int strip() const { return (id_ & 0x3F); }
/// get a compact index for arrays [TODO: NEEDS WORK]
int hashedIndex() const;
uint32_t denseIndex() const { return hashedIndex(); }
static bool validDenseIndex(uint32_t din) { return validHashIndex(din); }
static ESDetId detIdFromDenseIndex(uint32_t din) { return unhashIndex(din); }
/// get a DetId from a compact index for arrays
static ESDetId unhashIndex(int hi);
static bool validHashIndex(int hi) { return (hi < kSizeForDenseIndexing); }
/// check if a valid index combination
static bool validDetId(int istrip, int ixs, int iys, int iplane, int iz);
static void verify(int istrip, int ixs, int iys, int iplane, int iz);
static const int IX_MIN = 1;
static const int IY_MIN = 1;
static const int IX_MAX = 40;
static const int IY_MAX = 40;
static const int ISTRIP_MIN = 1;
static const int ISTRIP_MAX = 32;
static const int PLANE_MIN = 1;
static const int PLANE_MAX = 2;
static const int IZ_NUM = 2;
private:
enum {
kXYMAX = 1072,
kXYMIN = 1,
kXMAX = 40,
kYMAX = 40,
kXMIN = 1,
kYMIN = 1,
// now normalize to A-D notation for ease of use
kNa = IZ_NUM,
kNb = PLANE_MAX - PLANE_MIN + 1,
kNc = kXYMAX - kXYMIN + 1,
kNd = ISTRIP_MAX - ISTRIP_MIN + 1,
kLd = kNd,
kLc = kLd * kNc,
kLb = kLc * kNb,
kLa = kLb * kNa
};
static const unsigned short hxy1[kXMAX][kYMAX];
static const unsigned short hx1[kXYMAX];
static const unsigned short hy1[kXYMAX];
static const unsigned short hxy2[kXMAX][kYMAX];
static const unsigned short hx2[kXYMAX];
static const unsigned short hy2[kXYMAX];
public:
static constexpr int kSizeForDenseIndexing = kLa;
};
std::ostream& operator<<(std::ostream&, const ESDetId& id);
inline ESDetId::ESDetId(const DetId& gen) : DetId(gen) {
#ifdef EDM_ML_DEBUG
if (!gen.null() && (gen.det() != Ecal || gen.subdetId() != EcalPreshower)) {
throw cms::Exception("InvalidDetId");
}
#endif
}
inline ESDetId& ESDetId::operator=(const DetId& gen) {
#ifdef EDM_ML_DEBUG
if (!gen.null() && (gen.det() != Ecal || gen.subdetId() != EcalPreshower)) {
throw cms::Exception("InvalidDetId");
}
#endif
id_ = gen.rawId();
return *this;
}
#endif
|