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
|
#ifndef L1GCTINTERNEMCAND_H
#define L1GCTINTERNEMCAND_H
#include <ostream>
#include <string>
#include "DataFormats/L1CaloTrigger/interface/L1CaloEmCand.h"
/*! \class L1GctEmCand
* \brief Level-1 Trigger EM candidate within GCT processing
*
*/
/*! \author Jim Brooke
* \date June 2006
*/
class L1GctInternEmCand { //: public L1GctCand {
public:
/// default constructor (for vector initialisation etc.)
L1GctInternEmCand();
/// construct from raw data
L1GctInternEmCand(uint16_t data, bool iso, unsigned block, unsigned index, int16_t bx);
// construct from eta, phi etc
L1GctInternEmCand(
unsigned rank, unsigned eta, unsigned etaSgn, unsigned phi, bool iso, unsigned block, unsigned index, int16_t bx);
/// destructor (virtual to prevent compiler warnings)
virtual ~L1GctInternEmCand();
/// region associated with the candidate
L1CaloRegionDetId regionId() const;
/// name of object
std::string name() const;
/// was an object really found?
bool empty() const;
/// get the raw data
uint16_t raw() const { return m_data; }
/// get rank bits
unsigned rank() const { return m_data & 0x3f; }
/// get eta index - Barrel 1:4, Endcap 5:7, HF 8:11
unsigned etaIndex() const { return (m_data >> 6) & 0x1f; }
/// get eta sign (1 for -ve Z, 0 for +ve Z)
unsigned etaSign() const { return (m_data >> 10) & 0x1; }
/// get phi index (0-17)
unsigned phiIndex() const { return (m_data >> 11) & 0x1f; }
/// which stream did this come from
bool isolated() const { return m_iso; }
/// which capture block did this come from
unsigned capBlock() const { return m_captureBlock; }
/// what index within capture block
unsigned capIndex() const { return m_captureIndex; }
/// get bunch-crossing index
int16_t bx() const { return m_bx; }
/// equality operator
int operator==(const L1GctInternEmCand& c) const {
return ((m_data == c.raw() && m_iso == c.isolated()) || (this->empty() && c.empty()));
}
/// inequality operator
int operator!=(const L1GctInternEmCand& c) const {
return ((m_data != c.raw() || m_iso != c.isolated()) && (!this->empty() || !c.empty()));
}
private:
// set internal data from rank and region ieta, iphi
void construct(unsigned rank, unsigned eta, unsigned etaSgn, unsigned phi);
private:
uint16_t m_data;
uint16_t m_captureBlock;
uint8_t m_captureIndex;
bool m_iso;
int16_t m_bx;
};
std::ostream& operator<<(std::ostream& s, const L1GctInternEmCand& cand);
#endif
|