|
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
|
#include "DataFormats/L1CaloTrigger/interface/L1CaloRegion.h"
using std::dec;
using std::endl;
using std::hex;
using std::ostream;
// default constructor
L1CaloRegion::L1CaloRegion() : m_id(), m_data(0), m_bx(0) {}
// constructor for RCT emulator (HB/HE regions)
L1CaloRegion::L1CaloRegion(
unsigned et, bool overFlow, bool tauVeto, bool mip, bool quiet, unsigned crate, unsigned card, unsigned rgn)
: m_id(crate, card, rgn),
m_data(0), // over-ridden below
m_bx(0) {
pack(et, overFlow, tauVeto, mip, quiet);
}
// constructor for RCT emulator (HF regions)
L1CaloRegion::L1CaloRegion(unsigned et, bool fineGrain, unsigned crate, unsigned rgn)
: m_id(crate, 999, rgn),
m_data(0), // over-ridden below
m_bx(0) {
pack((et & 0xff), (et >= 0xff), fineGrain, false, false);
}
// construct from global eta, phi indices
L1CaloRegion::L1CaloRegion(
unsigned et, bool overFlow, bool fineGrain, bool mip, bool quiet, unsigned ieta, unsigned iphi)
: m_id(ieta, iphi),
m_data(0), // over-ridden below
m_bx(0) {
pack(et, overFlow, fineGrain, mip, quiet);
}
//constructor for unpacking
L1CaloRegion::L1CaloRegion(uint16_t data, unsigned ieta, unsigned iphi, int16_t bx)
: m_id(ieta, iphi), m_data(data), m_bx(bx) {}
// destructor
L1CaloRegion::~L1CaloRegion() {}
// named ctors
// for HB/HE from RCT indices
L1CaloRegion L1CaloRegion::makeHBHERegion(
unsigned et, bool overFlow, bool tauVeto, bool mip, bool quiet, unsigned crate, unsigned card, unsigned rgn) {
L1CaloRegion r;
r.setRegionId(L1CaloRegionDetId(crate, card, rgn));
r.setBx(0);
r.pack(et, overFlow, tauVeto, mip, quiet);
return r;
}
// for HF from RCT indices
L1CaloRegion L1CaloRegion::makeHFRegion(unsigned et, bool fineGrain, unsigned crate, unsigned rgn) {
L1CaloRegion r;
r.setRegionId(L1CaloRegionDetId(crate, 999, rgn));
r.setBx(0);
r.pack((et & 0xff), (et >= 0xff), fineGrain, false, false);
return r;
}
// HB/HE/HF from GCT indices
L1CaloRegion L1CaloRegion::makeRegionFromGctIndices(
unsigned et, bool overFlow, bool fineGrain, bool mip, bool quiet, unsigned ieta, unsigned iphi) {
L1CaloRegion r;
r.setRegionId(L1CaloRegionDetId(ieta, iphi));
r.setBx(0);
r.pack(et, overFlow, fineGrain, mip, quiet);
return r;
}
//constructor for unpacking
L1CaloRegion L1CaloRegion::makeRegionFromUnpacker(
uint16_t data, unsigned ieta, unsigned iphi, uint16_t block, uint16_t index, int16_t bx) {
L1CaloRegion r;
r.setRegionId(L1CaloRegionDetId(ieta, iphi));
r.setRawData(data);
r.setCaptureBlock(block);
r.setCaptureIndex(index);
r.setBx(bx);
return r;
}
L1CaloRegion L1CaloRegion::makeGctJetRegion(const unsigned et,
const bool overFlow,
const bool fineGrain,
const unsigned ieta,
const unsigned iphi,
const int16_t bx) {
L1CaloRegion r;
r.setRegionId(L1CaloRegionDetId(ieta, iphi));
r.setBx(bx);
r.pack12BitsEt(et, overFlow, fineGrain, false, false);
return r;
}
// set BX
void L1CaloRegion::setBx(int16_t bx) { m_bx = bx; }
// set mip bit
void L1CaloRegion::setMip(bool mip) {
if (mip) {
m_data |= 0x1000;
} else {
m_data &= 0xefff;
}
}
// set quiet bit
void L1CaloRegion::setQuiet(bool quiet) {
if (quiet) {
m_data |= 0x2000;
} else {
m_data &= 0xdfff;
}
}
void L1CaloRegion::pack(unsigned et, bool overFlow, bool fineGrain, bool mip, bool quiet) {
bool checkOvF = overFlow || (et >= 0x400);
m_data = (et & 0x3ff) | ((checkOvF) ? 0x400 : 0x0) | ((fineGrain) ? 0x800 : 0x0) | ((mip) ? 0x1000 : 0x0) |
((quiet) ? 0x2000 : 0x0);
}
void L1CaloRegion::pack12BitsEt(unsigned et, bool overFlow, bool fineGrain, bool mip, bool quiet) {
bool checkOvF = overFlow || (et >= 0x400);
m_data = (et & 0xfff) | ((checkOvF) ? 0x400 : 0x0) | ((fineGrain) ? 0x800 : 0x0) | ((mip) ? 0x1000 : 0x0) |
((quiet) ? 0x2000 : 0x0);
}
// print to stream
ostream& operator<<(ostream& os, const L1CaloRegion& reg) {
os << "L1CaloRegion:";
os << " Et=" << reg.et();
os << " o/f=" << reg.overFlow();
os << " f/g=" << reg.fineGrain();
os << " tau=" << reg.tauVeto() << endl;
os << " RCT crate=" << reg.rctCrate();
os << " RCT card=" << reg.rctCard();
os << " RCT rgn=" << reg.rctRegionIndex();
os << " RCT eta=" << reg.rctEta();
os << " RCT phi=" << reg.rctPhi() << endl;
os << " GCT eta=" << reg.gctEta();
os << " GCT phi=" << reg.gctPhi() << endl;
os << hex << " cap block=" << reg.capBlock() << dec << ", index=" << reg.capIndex() << ", BX=" << reg.bx();
return os;
}
|