|
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
|
#ifndef DataFormats_L1TCalorimeterPhase2_CaloCrystalsCluster_h
#define DataFormats_L1TCalorimeterPhase2_CaloCrystalsCluster_h
#include <ap_int.h>
#include <vector>
#include <map>
#include <string>
#include <algorithm>
#include "DataFormats/L1Trigger/interface/L1Candidate.h"
#include "DataFormats/DetId/interface/DetId.h"
namespace l1tp2 {
class CaloCrystalCluster : public l1t::L1Candidate {
public:
CaloCrystalCluster()
: l1t::L1Candidate(),
calibratedPt_(0.),
hovere_(0.),
iso_(0.),
puCorrPt_(0.),
bremStrength_(0.),
e2x2_(0.),
e2x5_(0.),
e3x5_(0.),
e5x5_(0.),
standaloneWP_(0.),
electronWP98_(0.),
photonWP80_(0.),
electronWP90_(0.),
looseL1TkMatchWP_(0.),
stage2effMatch_(0.) {}
CaloCrystalCluster(const PolarLorentzVector &p4,
float calibratedPt,
float hovere,
float iso,
DetId seedCrystal,
float puCorrPt = 0.,
float bremStrength = 0.,
float e2x2 = 0.,
float e2x5 = 0.,
float e3x5 = 0.,
float e5x5 = 0.,
bool standaloneWP = false,
bool electronWP98 = false,
bool photonWP80 = false,
bool electronWP90 = false,
bool looseL1TkMatchWP = false,
bool stage2effMatch = false)
: l1t::L1Candidate(p4),
calibratedPt_(calibratedPt),
hovere_(hovere),
iso_(iso),
seedCrystal_(seedCrystal),
puCorrPt_(puCorrPt),
bremStrength_(bremStrength),
e2x2_(e2x2),
e2x5_(e2x5),
e3x5_(e3x5),
e5x5_(e5x5),
standaloneWP_(standaloneWP),
electronWP98_(electronWP98),
photonWP80_(photonWP80),
electronWP90_(electronWP90),
looseL1TkMatchWP_(looseL1TkMatchWP),
stage2effMatch_(stage2effMatch) {}
~CaloCrystalCluster() override {}
inline float calibratedPt() const { return calibratedPt_; };
inline float hovere() const { return hovere_; };
inline float isolation() const { return iso_; };
inline float puCorrPt() const { return puCorrPt_; };
inline float bremStrength() const { return bremStrength_; };
inline DetId seedCrystal() const { return seedCrystal_; };
void setCrystalPtInfo(std::vector<float> info) {
std::sort(info.begin(), info.end());
std::reverse(info.begin(), info.end());
crystalPt_ = std::move(info);
};
void setExperimentalParams(const std::map<std::string, float> ¶ms) { experimentalParams_ = params; };
const std::map<std::string, float> &getExperimentalParams() const { return experimentalParams_; };
inline float experimentalParam(const std::string &name) const {
auto iter = experimentalParams_.find(name);
if (iter != experimentalParams_.end()) {
return iter->second;
} else {
warningNoMapping(name);
return -99.;
}
};
inline float e2x2() const { return e2x2_; };
inline float e2x5() const { return e2x5_; };
inline float e3x5() const { return e3x5_; };
inline float e5x5() const { return e5x5_; };
inline float standaloneWP() const { return standaloneWP_; };
inline float electronWP98() const { return electronWP98_; };
inline float photonWP80() const { return photonWP80_; };
inline float electronWP90() const { return electronWP90_; };
inline float looseL1TkMatchWP() const { return looseL1TkMatchWP_; };
inline float stage2effMatch() const { return stage2effMatch_; };
// The index range depends on the algorithm eta,phi window, currently 3x5
// The pt should always be ordered.
inline float crystalPt(unsigned int index) const { return (index < crystalPt_.size()) ? crystalPt_[index] : 0.; };
private:
static void warningNoMapping(const std::string &name);
// pT calibrated to Stage-2 (Phase-I) L1EG Objects. NOTE
// all working points are defined with respect to cluster.pt(),
// not cluster.calibratedPt()
float calibratedPt_;
// HCal energy in region behind cluster (for size, look in producer) / ECal energy in cluster
float hovere_;
// ECal isolation (for outer window size, again look in producer)
float iso_;
// DetId of seed crystal used to make cluster (could be EBDetId or EEDetId)
DetId seedCrystal_;
// Pileup-corrected energy deposit, not studied carefully yet, don't use
float puCorrPt_;
// Bremstrahlung strength, should be proportional to the likelihood of a brem.
float bremStrength_;
// Shower shape variable - max 2x2 energy containing seed crystal
float e2x2_;
// Shower shape variable - max 2x5 energy containing seed crystal, phi centered
float e2x5_;
// Shower shape variable - 3x5 energy containing seed crystal, phi centered
float e3x5_;
// Shower shape variable - 5x5 energy containing centered on seed crystal
float e5x5_;
// Standalone L1EG WP
bool standaloneWP_;
// 98% efficient electron WP, for electrons above 35 GeV
bool electronWP98_;
// 80% efficient photon WP, for photons above 35 GeV
bool photonWP80_;
// 90% efficient electron based WP, early rise to efficiency plateau
bool electronWP90_;
// loose isolation and shower shape requirements to be used in conjunction with L1Trk matching
bool looseL1TkMatchWP_;
// Stage-2 L1EG efficiency matched WP, for rate comparisons
bool stage2effMatch_;
// Crystal pt (in order of strength) for all crystals in the cluster
std::vector<float> crystalPt_;
// For investigating novel algorithm parameters
std::map<std::string, float> experimentalParams_;
};
// Concrete collection of output objects (with extra tuning information)
typedef std::vector<l1tp2::CaloCrystalCluster> CaloCrystalClusterCollection;
} // namespace l1tp2
#endif
|