|
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
|
#include "DataFormats/L1GlobalCaloTrigger/interface/L1GctEmCand.h"
#include <iostream>
using std::dec;
using std::hex;
using std::ostream;
using std::string;
// default constructor
L1GctEmCand::L1GctEmCand() : m_data(0), m_iso(false), m_captureBlock(0), m_captureIndex(0), m_bx(0) {}
// construct from raw data, no source (i.e. no capBlock/capIndex); used in GT
L1GctEmCand::L1GctEmCand(uint16_t rawData, bool iso)
: m_data(rawData & 0x7fff), // 0x7fff is to mask off bit 15, which is not data that needs to be stored
m_iso(iso),
m_captureBlock(0),
m_captureIndex(0),
m_bx(0) {}
// construct from raw data with source - used in GCT unpacker
L1GctEmCand::L1GctEmCand(uint16_t rawData, bool iso, uint16_t block, uint16_t index, int16_t bx)
: m_data(rawData & 0x7fff), // 0x7fff is to mask off bit 15, which is not data that needs to be stored
m_iso(iso),
m_captureBlock(block & 0xfff),
m_captureIndex(index & 0xff),
m_bx(bx) {}
// construct from content - used in GCT emulator
// eta = -6 to -0, +0 to +6. Sign is bit 3, 1 means -ve Z, 0 means +ve Z
L1GctEmCand::L1GctEmCand(unsigned rank, unsigned phi, unsigned eta, bool iso)
: m_data(0), // override below
m_iso(iso),
m_captureBlock(0),
m_captureIndex(0),
m_bx(0)
{
construct(rank, eta, phi);
}
// construct from content, with source (i.e. capBlock/capIndex); will be used in GCT emulator one day?
// eta = -6 to -0, +0 to +6. Sign is bit 3, 1 means -ve Z, 0 means +ve Z
L1GctEmCand::L1GctEmCand(unsigned rank, unsigned phi, unsigned eta, bool iso, uint16_t block, uint16_t index, int16_t bx)
: m_data(0), // override below
m_iso(iso),
m_captureBlock(block & 0xfff),
m_captureIndex(index & 0xff),
m_bx(bx) {
construct(rank, eta, phi);
}
// construct from RCT output candidate
L1GctEmCand::L1GctEmCand(L1CaloEmCand& c)
: m_data(0), // override below
m_iso(c.isolated()),
m_captureBlock(0),
m_captureIndex(0),
m_bx(c.bx()) {
unsigned eta = ((c.regionId().rctEta() & 0x7) | (c.regionId().ieta() < 11 ? 0x8 : 0x0));
construct(c.rank(), eta, c.regionId().iphi());
}
// destructor
L1GctEmCand::~L1GctEmCand() {}
// name of candidate type
string L1GctEmCand::name() const { return (isolated() ? "iso EM" : "non iso EM"); }
// return region object
L1CaloRegionDetId L1GctEmCand::regionId() const {
// get global eta
unsigned eta = (etaSign() == 1 ? 10 - (etaIndex() & 0x7) : 11 + (etaIndex() & 0x7));
return L1CaloRegionDetId(eta, phiIndex());
}
// construct from rank, eta, phi
void L1GctEmCand::construct(unsigned rank, unsigned eta, unsigned phi) {
if (rank > 0) {
m_data = (rank & 0x3f) + ((eta & 0xf) << 6) + ((phi & 0x1f) << 10);
} else {
// Default values for zero rank electrons,
// different in hardware for positive and negative eta
if ((eta & 0x8) == 0) {
m_data = 0x7000;
} else {
m_data = 0x7400;
}
}
}
// pretty print
ostream& operator<<(ostream& s, const L1GctEmCand& cand) {
s << "L1GctEmCand : ";
s << "rank=" << cand.rank();
s << ", etaSign=" << cand.etaSign() << ", eta=" << (cand.etaIndex() & 0x7) << ", phi=" << cand.phiIndex();
s << ", iso=" << cand.isolated();
s << hex << " cap block=" << cand.capBlock() << dec << ", index=" << cand.capIndex() << ", BX=" << cand.bx();
return s;
}
|