|
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
|
#include "DataFormats/L1TCalorimeter/interface/CaloCluster.h"
l1t::CaloCluster::CaloCluster(const LorentzVector p4, int pt, int eta, int phi)
: L1Candidate(p4, pt, eta, phi),
m_clusterFlags(0x7FF) // first 11 flags at 1
{}
l1t::CaloCluster::~CaloCluster() {}
void l1t::CaloCluster::setClusterFlag(ClusterFlag flag, bool val) {
if (val) {
m_clusterFlags |= (0x1 << flag);
} else {
m_clusterFlags &= ~(0x1 << flag);
}
};
void l1t::CaloCluster::setHwPtEm(int pt) { m_hwPtEm = pt; }
void l1t::CaloCluster::setHwPtHad(int pt) { m_hwPtHad = pt; }
void l1t::CaloCluster::setHwSeedPt(int pt) { m_hwSeedPt = pt; }
void l1t::CaloCluster::setFgEta(int fgEta) { m_fgEta = fgEta; }
void l1t::CaloCluster::setFgPhi(int fgPhi) { m_fgPhi = fgPhi; }
void l1t::CaloCluster::setHOverE(int hOverE) { m_hOverE = hOverE; }
void l1t::CaloCluster::setFgECAL(int fgECAL) { m_fgECAL = fgECAL; }
bool l1t::CaloCluster::checkClusterFlag(ClusterFlag flag) const { return (m_clusterFlags & (0x1 << flag)); };
bool l1t::CaloCluster::isValid() const { return (checkClusterFlag(INCLUDE_SEED)); }
int l1t::CaloCluster::hwPtEm() const { return m_hwPtEm; }
int l1t::CaloCluster::hwPtHad() const { return m_hwPtHad; }
int l1t::CaloCluster::hwSeedPt() const { return m_hwSeedPt; }
int l1t::CaloCluster::fgEta() const { return m_fgEta; }
int l1t::CaloCluster::fgPhi() const { return m_fgPhi; }
int l1t::CaloCluster::hOverE() const { return m_hOverE; }
int l1t::CaloCluster::fgECAL() const { return m_fgECAL; }
bool l1t::CaloCluster::operator<(const CaloCluster& cl) const {
bool res = false;
// Favour high pT
if (hwPt() < cl.hwPt())
res = true;
else if (hwPt() == cl.hwPt()) {
// Favour central clusters
if (abs(hwEta()) > abs(cl.hwEta()))
res = true;
else if (abs(hwEta()) == abs(cl.hwEta())) {
// Favour small phi (arbitrary)
if (hwPhi() > cl.hwPhi())
res = true;
}
}
return res;
}
|