| 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; |
| 218 |
|
vector<vector<map<string, TH2D*>>> twoDHists_; |
| 219 |
|
vector<TTree*> BNTrees_; // one tree per channel |
| 220 |
|
map<string, vector<float>> BNTreeBranchVals_; // data structure to hold the values of the branches to be stored in the BNTrees_; the string is the name of a variable |
| 221 |
+ |
long BNTreeBranchVals_evtLong_; // event number |
| 222 |
+ |
int BNTreeBranchVals_runInt_; // run number |
| 223 |
+ |
int BNTreeBranchVals_lumiInt_; // lumi number |
| 224 |
|
|
| 225 |
|
edm::Service<TFileService> fs_; |
| 226 |
|
|
| 269 |
|
vector<string> objectsToGet; |
| 270 |
|
vector<string> objectsToCut; |
| 271 |
|
vector<string> objectsToPlot; |
| 272 |
+ |
vector<string> objectsToFlag; |
| 273 |
+ |
|
| 274 |
|
vector<channel> channels; |
| 275 |
|
vector<histogram> histograms; |
| 276 |
|
vector<BranchSpecs> treeBranches_; |
| 286 |
|
StopCTauWeight *stopCTauWeight_; |
| 287 |
|
double stopCTauScaleFactor_; |
| 288 |
|
|
| 289 |
+ |
ofstream* findEventsLog; |
| 290 |
+ |
bool isFirstEvent_; |
| 291 |
+ |
|
| 292 |
|
template <class InputCollection> void setObjectFlags(cut &, uint, flagMap &, flagMap &, InputCollection, string); |
| 293 |
< |
template <class InputCollection1, class InputCollection2> void setObjectFlags(cut &, uint, flagMap &, flagMap &, InputCollection1, InputCollection2, flagPair, flagPair, string); |
| 293 |
> |
template <class InputCollection1, class InputCollection2> void setObjectFlags(cut &, uint, flagMap &, flagMap &, InputCollection1, InputCollection2, string); |
| 294 |
|
|
| 295 |
|
template <class InputCollection> void assignTreeBranch(BranchSpecs brSpecs, InputCollection inputCollection, flagPair flags); |
| 296 |
|
template <class InputCollection> void fill1DHistogram(TH1*, histogram, InputCollection, flagPair, double); |
| 297 |
< |
template <class InputCollection1, class InputCollection2> void fill1DHistogram(TH1*, histogram, InputCollection1, InputCollection2, flagPair, flagPair, flagPair, double); |
| 297 |
> |
template <class InputCollection1, class InputCollection2> void fill1DHistogram(TH1*, histogram, InputCollection1, InputCollection2, flagPair, double); |
| 298 |
|
template <class InputCollection> void fill2DHistogram(TH2*, histogram, InputCollection, flagPair, double); |
| 299 |
< |
template <class InputCollection1, class InputCollection2> void fill2DHistogram(TH2*, histogram, InputCollection1, InputCollection2, flagPair, flagPair, flagPair, double); |
| 299 |
> |
template <class InputCollection1, class InputCollection2> void fill2DHistogram(TH2*, histogram, InputCollection1, InputCollection2, flagPair, double); |
| 300 |
|
bool getPreviousCumulativeFlags(uint currentCutIndex, flagMap &individualFlags, string obj1Type, uint object1, string flagType); |
| 301 |
|
|
| 302 |
|
|
| 312 |
|
const BNmet *chosenMET (); |
| 313 |
|
const BNelectron *chosenElectron (); |
| 314 |
|
const BNmuon *chosenMuon (); |
| 315 |
+ |
double chosenHT (); |
| 316 |
+ |
pair<const BNmuon *, const BNmuon*> leadMuonPair (); |
| 317 |
+ |
pair<const BNelectron *, const BNelectron*> leadElectronPair (); |
| 318 |
|
|
| 319 |
|
}; |
| 320 |
|
|
