15 |
|
#include "TH1D.h" |
16 |
|
#include "TH2D.h" |
17 |
|
#include "TLorentzVector.h" |
18 |
+ |
#include "TVector3.h" |
19 |
+ |
#include "TVector2.h" |
20 |
|
#include "TTree.h" |
21 |
|
|
22 |
|
#include "FWCore/Framework/interface/EDProducer.h" |
77 |
|
double valueLookup (const BNelectron* object1, const BNjet* object2, string variable, string function, string &stringValue); |
78 |
|
double valueLookup (const BNelectron* object1, const BNphoton* object2, string variable, string function, string &stringValue); |
79 |
|
double valueLookup (const BNmuon* object1, const BNphoton* object2, string variable, string function, string &stringValue); |
80 |
+ |
double valueLookup (const BNmuon* object1, const BNevent* object2, string variable, string function, string &stringValue); |
81 |
|
double valueLookup (const BNphoton* object1, const BNjet* object2, string variable, string function, string &stringValue); |
82 |
|
double valueLookup (const BNelectron* object1, const BNtrack* object2, string variable, string function, string &stringValue); |
83 |
|
double valueLookup (const BNmuon* object1, const BNjet* object2, string variable, string function, string &stringValue); |
106 |
|
|
107 |
|
int getTrkIsIso (const BNtrack* track1, const BNtrackCollection* trackColl); |
108 |
|
double getTrkPtTrue (const BNtrack* track1, const BNmcparticleCollection* genPartColl); |
109 |
+ |
double getHt (const BNjetCollection* jetColl); |
110 |
|
double getTrkPtRes (const BNtrack* track1); |
111 |
|
double getTrkIsIso (const BNtrack* track1); |
112 |
|
double getTrkCaloTotRhoCorr(const BNtrack* track); |
113 |
+ |
double getTrkDepTrkRp5RhoCorr(const BNtrack* track); |
114 |
+ |
double getTrkDepTrkRp3RhoCorr(const BNtrack* track); |
115 |
|
void WriteDeadEcal (); |
116 |
|
int getTrkIsMatchedDeadEcal (const BNtrack* track1); |
117 |
|
|
131 |
|
//pair(bool) first bool counts towards the event passing, second bool determines whether to plot the object |
132 |
|
typedef map<string, vector<vector<pair<bool,bool> > > > flagMap; |
133 |
|
|
134 |
+ |
//flagPair (corresponds to flagMap): |
135 |
+ |
//vector corresponds to each object in input collection |
136 |
+ |
//pair(bool) first bool counts towards the event passing, second bool determines whether to plot the object |
137 |
|
typedef vector<pair<bool,bool> > flagPair; |
138 |
|
|
139 |
|
//counterMap: |
176 |
|
bool applyLeptonSF_; |
177 |
|
bool applyBtagSF_; |
178 |
|
int minBtag_; |
170 |
– |
int maxBtag_; |
179 |
|
bool printEventInfo_; |
180 |
|
bool printAllTriggers_; |
181 |
|
bool useTrackCaloRhoCorr_; // to use the calo-based rho correction for the by-hand calculation of the track isolation energy |
182 |
|
vector<double> stopCTau_; |
183 |
|
bool GetPlotsAfterEachCut_; |
184 |
+ |
int verbose_; |
185 |
|
|
186 |
|
struct DeadEcal { |
187 |
|
double etaEcal; |
266 |
|
vector<string> objectsToGet; |
267 |
|
vector<string> objectsToCut; |
268 |
|
vector<string> objectsToPlot; |
269 |
+ |
vector<string> objectsToFlag; |
270 |
+ |
|
271 |
|
vector<channel> channels; |
272 |
|
vector<histogram> histograms; |
273 |
|
vector<BranchSpecs> treeBranches_; |
283 |
|
StopCTauWeight *stopCTauWeight_; |
284 |
|
double stopCTauScaleFactor_; |
285 |
|
|
286 |
+ |
ofstream* findEventsLog; |
287 |
+ |
bool isFirstEvent_; |
288 |
+ |
|
289 |
|
template <class InputCollection> void setObjectFlags(cut &, uint, flagMap &, flagMap &, InputCollection, string); |
290 |
< |
template <class InputCollection1, class InputCollection2> void setObjectFlags(cut &, uint, flagMap &, flagMap &, InputCollection1, InputCollection2, flagPair, flagPair, string); |
290 |
> |
template <class InputCollection1, class InputCollection2> void setObjectFlags(cut &, uint, flagMap &, flagMap &, InputCollection1, InputCollection2, string); |
291 |
|
|
292 |
|
template <class InputCollection> void assignTreeBranch(BranchSpecs brSpecs, InputCollection inputCollection, flagPair flags); |
293 |
|
template <class InputCollection> void fill1DHistogram(TH1*, histogram, InputCollection, flagPair, double); |
294 |
< |
template <class InputCollection1, class InputCollection2> void fill1DHistogram(TH1*, histogram, InputCollection1, InputCollection2, flagPair, flagPair, flagPair, double); |
294 |
> |
template <class InputCollection1, class InputCollection2> void fill1DHistogram(TH1*, histogram, InputCollection1, InputCollection2, flagPair, double); |
295 |
|
template <class InputCollection> void fill2DHistogram(TH2*, histogram, InputCollection, flagPair, double); |
296 |
< |
template <class InputCollection1, class InputCollection2> void fill2DHistogram(TH2*, histogram, InputCollection1, InputCollection2, flagPair, flagPair, flagPair, double); |
296 |
> |
template <class InputCollection1, class InputCollection2> void fill2DHistogram(TH2*, histogram, InputCollection1, InputCollection2, flagPair, double); |
297 |
|
bool getPreviousCumulativeFlags(uint currentCutIndex, flagMap &individualFlags, string obj1Type, uint object1, string flagType); |
298 |
|
|
299 |
|
|
309 |
|
const BNmet *chosenMET (); |
310 |
|
const BNelectron *chosenElectron (); |
311 |
|
const BNmuon *chosenMuon (); |
312 |
+ |
double chosenHT (); |
313 |
+ |
pair<const BNmuon *, const BNmuon*> leadMuonPair (); |
314 |
+ |
pair<const BNelectron *, const BNelectron*> leadElectronPair (); |
315 |
|
|
316 |
|
}; |
317 |
|
|