25
25
26
26
class Run3ScoutingPFJetIDSelectionFunctor {
27
27
public: // interface
28
- enum Version_t {
29
- RUN3Scouting,
30
- N_VERSIONS
31
- };
28
+ enum Version_t { RUN3Scouting, N_VERSIONS };
32
29
enum Quality_t { TIGHT, TIGHTLEPVETO, N_QUALITY };
33
30
34
31
Run3ScoutingPFJetIDSelectionFunctor () {}
@@ -45,17 +42,19 @@ class Run3ScoutingPFJetIDSelectionFunctor {
45
42
if (versionStr == " RUN3Scouting" )
46
43
version_ = RUN3Scouting;
47
44
else {
48
- std::cout << " JetID version not specified -- setting default to RUN3Scouting" << std::endl;
49
- version_ = RUN3Scouting; // set RUN3Scouting as default //this is extremely unsafe --> similarly it's done in PFJetIDSelectionFunctor.h
45
+ edm::LogWarning (" BadJetIDQuality" ) << " JetID quality not specified -- setting default to RUN3Scouting" ;
46
+ version_ =
47
+ RUN3Scouting; // set RUN3Scouting as default //this is extremely unsafe --> similarly it's done in PFJetIDSelectionFunctor.h
50
48
}
51
-
49
+
52
50
if (qualityStr == " TIGHT" )
53
51
quality_ = TIGHT;
54
52
else if (qualityStr == " TIGHTLEPVETO" )
55
53
quality_ = TIGHTLEPVETO;
56
54
else {
57
- edm::LogWarning (" BadJetIDQuality" ) << " JetID quality not specified -- setting default to RUN3Scouting" ;
58
- quality_ = TIGHT; // set TIGHT as default //this is extremely unsafe --> similarly it's done in PFJetIDSelectionFunctor.h
55
+ edm::LogWarning (" BadJetIDQuality" ) << " JetID quality not specified -- setting default to TIGHT" ;
56
+ quality_ =
57
+ TIGHT; // set TIGHT as default //this is extremely unsafe --> similarly it's done in PFJetIDSelectionFunctor.h
59
58
}
60
59
61
60
initCuts ();
@@ -72,18 +71,17 @@ class Run3ScoutingPFJetIDSelectionFunctor {
72
71
edm::ParameterSetDescription desc;
73
72
74
73
desc.ifValue (
75
- edm::ParameterDescription<std::string>(" version" , " RUN3Scouting" , true , edm::Comment (" " )), // default "version"
76
- edm::allowedValues<std::string>(" RUN3Scouting" )); // more options about "version"
77
- desc.ifValue (edm::ParameterDescription<std::string>(" quality" , " TIGHT" , true , edm::Comment (" " )), // default "quality"
78
- edm::allowedValues<std::string>(" TIGHT" , " TIGHTLEPVETO" )); // more options about "quality"
74
+ edm::ParameterDescription<std::string>(" version" , " RUN3Scouting" , true , edm::Comment (" " )), // default "version"
75
+ edm::allowedValues<std::string>(" RUN3Scouting" )); // more options about "version"
76
+ desc.ifValue (
77
+ edm::ParameterDescription<std::string>(" quality" , " TIGHT" , true , edm::Comment (" " )), // default "quality"
78
+ edm::allowedValues<std::string>(" TIGHT" , " TIGHTLEPVETO" )); // more options about "quality"
79
79
desc.addOptional <std::vector<std::string>>(" cutsToIgnore" )->setComment (" " );
80
80
81
81
return desc;
82
82
}
83
83
84
-
85
84
bool operator ()(const Run3ScoutingPFJet &jet) {
86
-
87
85
// cache some variables
88
86
// float pt = 0.0;
89
87
float eta = 0.0 ;
@@ -107,18 +105,20 @@ class Run3ScoutingPFJetIDSelectionFunctor {
107
105
eta = scoutingpfjet->eta ();
108
106
// phi = scoutingpfjet->phi();
109
107
double jetEnergyUncorrected = scoutingpfjet->chargedHadronEnergy () + scoutingpfjet->neutralHadronEnergy () +
110
- scoutingpfjet->photonEnergy () + scoutingpfjet->electronEnergy () + scoutingpfjet-> muonEnergy () +
111
- scoutingpfjet->HFEMEnergy ();
108
+ scoutingpfjet->photonEnergy () + scoutingpfjet->electronEnergy () +
109
+ scoutingpfjet->muonEnergy () + scoutingpfjet-> HFEMEnergy ();
112
110
if (jetEnergyUncorrected > 0 .) {
113
111
chf = scoutingpfjet->chargedHadronEnergy () / jetEnergyUncorrected;
114
112
nhf = scoutingpfjet->neutralHadronEnergy () / jetEnergyUncorrected;
115
113
// cef= scoutingpfjet->electronEnergy() / jetEnergyUncorrected; // for now: electron energy is 0, since HLT scouting jets by construction don't contain electrons
116
114
nef = (scoutingpfjet->photonEnergy () + scoutingpfjet->HFEMEnergy ()) / jetEnergyUncorrected;
117
115
muf = scoutingpfjet->muonEnergy () / jetEnergyUncorrected;
118
116
}
119
-
117
+
120
118
nch = scoutingpfjet->chargedHadronMultiplicity () + scoutingpfjet->electronMultiplicity ();
121
- nconstituents = scoutingpfjet->chargedHadronMultiplicity () + scoutingpfjet->electronMultiplicity () + scoutingpfjet->neutralHadronMultiplicity () + scoutingpfjet->photonMultiplicity () + scoutingpfjet->HFEMMultiplicity ();
119
+ nconstituents = scoutingpfjet->chargedHadronMultiplicity () + scoutingpfjet->electronMultiplicity () +
120
+ scoutingpfjet->neutralHadronMultiplicity () + scoutingpfjet->photonMultiplicity () +
121
+ scoutingpfjet->HFEMMultiplicity ();
122
122
// nneutrals = scoutingpfjet->neutralHadronMultiplicity() + scoutingpfjet->photonMultiplicity() + scoutingpfjet->HFEMMultiplicity();
123
123
}
124
124
@@ -159,7 +159,6 @@ class Run3ScoutingPFJetIDSelectionFunctor {
159
159
}
160
160
161
161
private: // member variables
162
-
163
162
int nConstituents_ = 0 ;
164
163
double CHF_ = 0.0 ;
165
164
double NHF_ = 0.0 ;
@@ -170,7 +169,7 @@ class Run3ScoutingPFJetIDSelectionFunctor {
170
169
double MUF_TR_ = 0.0 ;
171
170
double NEF_EC_ = 0.0 ;
172
171
double NEF_FW_ = 0.0 ;
173
-
172
+
174
173
void initCuts () {
175
174
if (quality_ == TIGHT) {
176
175
if (version_ == RUN3Scouting) {
@@ -200,10 +199,8 @@ class Run3ScoutingPFJetIDSelectionFunctor {
200
199
}
201
200
}
202
201
203
-
204
202
Version_t version_;
205
203
Quality_t quality_;
206
-
207
204
};
208
205
209
206
#endif
0 commit comments