|
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
143
144
145
146
147
148
149
150
151
152
153
154
|
#include "DataFormats/L1GlobalCaloTrigger/interface/L1GctInternHtMiss.h"
// PUBLIC METHODS
// Default ctor
L1GctInternHtMiss::L1GctInternHtMiss() : type_(nulltype), capBlock_(0), capIndex_(0), bx_(0), data_(0) {}
// Destructor
L1GctInternHtMiss::~L1GctInternHtMiss() {}
// Named ctor for making missing Ht x-component object from unpacker raw data.
L1GctInternHtMiss L1GctInternHtMiss::unpackerMissHtx(const uint16_t capBlock,
const uint16_t capIndex,
const int16_t bx,
const uint32_t data) {
return L1GctInternHtMiss(miss_htx, capBlock, capIndex, bx, data & kSingleComponentRawMask);
}
// Named ctor for making missing Ht y-component object from unpacker raw data.
L1GctInternHtMiss L1GctInternHtMiss::unpackerMissHty(const uint16_t capBlock,
const uint16_t capIndex,
const int16_t bx,
const uint32_t data) {
return L1GctInternHtMiss(miss_hty, capBlock, capIndex, bx, data & kSingleComponentRawMask);
}
// Named ctor for making missing Ht x & y components object from unpacker raw data.
L1GctInternHtMiss L1GctInternHtMiss::unpackerMissHtxHty(const uint16_t capBlock,
const uint16_t capIndex,
const int16_t bx,
const uint32_t data) {
return L1GctInternHtMiss(miss_htx_and_hty, capBlock, capIndex, bx, data & kDoubleComponentRawMask);
}
// Named ctor for making missing Ht x & y components object from unpacker raw data.
L1GctInternHtMiss L1GctInternHtMiss::emulatorJetMissHt(const int htx,
const int hty,
const bool overFlow,
const int16_t bx) {
int32_t xdata = (htx & kJetFinderComponentHtMask);
int32_t ydata = (hty & kJetFinderComponentHtMask) << kDoubleComponentHtyShift;
int32_t odata = 0;
if (overFlow || (htx >= kJetFinderComponentHtMask / 2) || (htx < -kJetFinderComponentHtMask / 2) ||
(hty >= kJetFinderComponentHtMask / 2) || (hty < -kJetFinderComponentHtMask / 2))
odata = kDoubleComponentOflowMask;
return L1GctInternHtMiss(jf_miss_htx_and_hty, 0, 0, bx, xdata | ydata | odata);
}
/// Named ctor for making missing Ht x & y components object from emulator (wheel input).
L1GctInternHtMiss L1GctInternHtMiss::emulatorMissHtxHty(const int htx,
const int hty,
const bool overFlow,
const int16_t bx) {
int32_t xdata = (htx & kDoubleComponentHtMask);
int32_t ydata = (hty & kDoubleComponentHtMask) << kDoubleComponentHtyShift;
int32_t odata = 0;
if (overFlow || (htx >= kDoubleComponentHtMask / 2) || (htx < -kDoubleComponentHtMask / 2) ||
(hty >= kDoubleComponentHtMask / 2) || (hty < -kDoubleComponentHtMask / 2))
odata = kDoubleComponentOflowMask;
return L1GctInternHtMiss(miss_htx_and_hty, 0, 0, bx, xdata | ydata | odata);
}
/// Named ctor for making missing Ht x component object from emulator
L1GctInternHtMiss L1GctInternHtMiss::emulatorMissHtx(const int htx, const bool overFlow, const int16_t bx) {
int32_t xdata = (htx & kSingleComponentHtMask);
int32_t odata = 0;
if (overFlow || (htx >= kSingleComponentHtMask / 2) || (htx < -kSingleComponentHtMask / 2))
odata = kSingleComponentOflowMask;
return L1GctInternHtMiss(miss_htx, 0, 0, bx, xdata | odata);
}
/// Named ctor for making missing Ht y component object from emulator
L1GctInternHtMiss L1GctInternHtMiss::emulatorMissHty(const int hty, const bool overFlow, const int16_t bx) {
int32_t ydata = (hty & kSingleComponentHtMask);
int32_t odata = 0;
if (overFlow || (hty >= kSingleComponentHtMask / 2) || (hty < -kSingleComponentHtMask / 2))
odata = kSingleComponentOflowMask;
return L1GctInternHtMiss(miss_hty, 0, 0, bx, ydata | odata);
}
// Get Ht x-component
int16_t L1GctInternHtMiss::htx() const {
if (type() == miss_htx) {
return static_cast<int16_t>(raw() & kSingleComponentHtMask);
}
if (type() == miss_htx_and_hty) {
return convert14BitTwosCompTo16Bit(raw() & kDoubleComponentHtMask);
}
return 0;
}
// Get Ht y-component
int16_t L1GctInternHtMiss::hty() const {
if (type() == miss_hty) {
return static_cast<int16_t>(raw() & kSingleComponentHtMask);
}
if (type() == miss_htx_and_hty) {
return convert14BitTwosCompTo16Bit((raw() >> kDoubleComponentHtyShift) & kDoubleComponentHtMask);
}
return 0;
}
// Get overflow
bool L1GctInternHtMiss::overflow() const {
if (type() == miss_htx || type() == miss_hty) {
return (raw() & kSingleComponentOflowMask) != 0;
}
if (type() == miss_htx_and_hty) {
return (raw() & kDoubleComponentOflowMask) != 0;
}
return false;
}
// PRIVATE METHODS
L1GctInternHtMiss::L1GctInternHtMiss(const L1GctInternHtMissType type,
const uint16_t capBlock,
const uint16_t capIndex,
const int16_t bx,
const uint32_t data)
: type_(type), capBlock_(capBlock), capIndex_(capIndex), bx_(bx), data_(data) {}
int16_t L1GctInternHtMiss::convert14BitTwosCompTo16Bit(const uint16_t data) const {
// If bit 13 is high, set bits 13, 14, 15 high.
if ((data & 0x2000) != 0) {
return static_cast<int16_t>(data | 0xe000);
}
// Else, bit 13 must be low, so set bits 13, 14, 15 low.
return static_cast<int16_t>(data & 0x1fff);
}
// PRETTY PRINTOUT OPERATOR
std::ostream& operator<<(std::ostream& os, const L1GctInternHtMiss& rhs) {
os << " L1GctInternHtMiss: htx=";
if (rhs.isThereHtx()) {
os << rhs.htx();
} else {
os << "n/a";
}
os << ", hty=";
if (rhs.isThereHty()) {
os << rhs.hty();
} else {
os << "n/a";
}
if (rhs.overflow()) {
os << "; overflow set";
}
os << "; cap block=0x" << std::hex << rhs.capBlock() << std::dec << ", index=" << rhs.capIndex()
<< ", BX=" << rhs.bx();
return os;
}
|