6 |
|
#include <string> |
7 |
|
#include <vector> |
8 |
|
|
9 |
+ |
#include "TH1.h" |
10 |
+ |
#include "TH2.h" |
11 |
|
#include "TH1D.h" |
12 |
|
#include "TH2D.h" |
13 |
|
#include "TLorentzVector.h" |
20 |
|
#include "DataFormats/Common/interface/Handle.h" |
21 |
|
#include "CommonTools/UtilAlgos/interface/TFileService.h" |
22 |
|
#include "FWCore/Framework/interface/MakerMacros.h" |
23 |
+ |
#include "DataFormats/Math/interface/deltaPhi.h" |
24 |
+ |
#include "DataFormats/Math/interface/deltaR.h" |
25 |
|
|
26 |
|
#include "ProductArea/BNcollections/interface/BNbxlumi.h" |
27 |
|
#include "ProductArea/BNcollections/interface/BNelectron.h" |
41 |
|
#include "ProductArea/BNcollections/interface/BNgenjet.h" |
42 |
|
|
43 |
|
#include "OSUT3Analysis/AnaTools/interface/CutFlow.h" |
44 |
+ |
#include "OSUT3Analysis/AnaTools/interface/PUWeight.h" |
45 |
|
|
46 |
|
|
47 |
|
|
59 |
|
bool evaluateTriggers (vector<string>,const BNtriggerCollection*); |
60 |
|
double applyFunction(string, double); |
61 |
|
|
62 |
< |
double valueLookup (const BNjet*, string, string); |
63 |
< |
double valueLookup (const BNmuon*, string, string); |
64 |
< |
double valueLookup (const BNelectron*, string, string); |
65 |
< |
double valueLookup (const BNevent*, string, string); |
66 |
< |
double valueLookup (const BNtau*, string, string); |
67 |
< |
double valueLookup (const BNmet*, string, string); |
68 |
< |
double valueLookup (const BNtrack*, string, string); |
69 |
< |
double valueLookup (const BNgenjet*, string, string); |
70 |
< |
double valueLookup (const BNmcparticle*, string, string); |
71 |
< |
double valueLookup (const BNprimaryvertex*, string, string); |
72 |
< |
double valueLookup (const BNbxlumi*, string, string); |
73 |
< |
double valueLookup (const BNphoton*, string, string); |
74 |
< |
double valueLookup (const BNsupercluster*, string, string); |
62 |
> |
double valueLookup (const BNjet* object, string variable, string function = ""); |
63 |
> |
double valueLookup (const BNmuon* object, string variable, string function = ""); |
64 |
> |
double valueLookup (const BNmuon* object1, const BNmuon* object2, string variable, string function = ""); |
65 |
> |
double valueLookup (const BNelectron* object, string variable, string function = ""); |
66 |
> |
double valueLookup (const BNelectron* object1, const BNelectron* object2, string variable, string function = ""); |
67 |
> |
double valueLookup (const BNelectron* object1, const BNmuon* object2, string variable, string function = ""); |
68 |
> |
double valueLookup (const BNevent* object, string variable, string function = ""); |
69 |
> |
double valueLookup (const BNtau* object, string variable, string function = ""); |
70 |
> |
double valueLookup (const BNmet* object, string variable, string function = ""); |
71 |
> |
double valueLookup (const BNtrack* object, string variable, string function = ""); |
72 |
> |
double valueLookup (const BNgenjet* object, string variable, string function = ""); |
73 |
> |
double valueLookup (const BNmcparticle* object, string variable, string function = ""); |
74 |
> |
double valueLookup (const BNprimaryvertex* object, string variable, string function = ""); |
75 |
> |
double valueLookup (const BNbxlumi* object, string variable, string function = ""); |
76 |
> |
double valueLookup (const BNphoton* object, string variable, string function = ""); |
77 |
> |
double valueLookup (const BNsupercluster* object, string variable, string function = ""); |
78 |
|
//BNskimbits |
71 |
– |
//BNtrigger |
79 |
|
//BNtrigobj |
80 |
|
|
81 |
|
vector<string> splitString (string); |
111 |
|
edm::InputTag photons_; |
112 |
|
edm::InputTag superclusters_; |
113 |
|
edm::InputTag triggers_; |
114 |
< |
|
115 |
< |
|
114 |
> |
std::string puFile_; |
115 |
> |
std::string dataPU_; |
116 |
> |
std::string dataset_; |
117 |
> |
std::string datasetType_; |
118 |
|
vector<edm::ParameterSet> channels_; |
119 |
+ |
vector<edm::ParameterSet> histogramSets_; |
120 |
+ |
bool plotAllObjectsInPassingEvents_; |
121 |
+ |
|
122 |
+ |
BNprimaryvertex *primaryVertex_; |
123 |
|
|
124 |
|
vector<map<string, TH1D*> > oneDHists_; |
125 |
+ |
vector<map<string, TH2D*> > twoDHists_; |
126 |
|
|
127 |
|
edm::Service<TFileService> fs_; |
128 |
|
|
129 |
|
CutFlow *masterCutFlow_; |
130 |
|
vector<CutFlow *> cutFlows_; |
131 |
|
|
132 |
+ |
typedef struct { |
133 |
+ |
string name; |
134 |
+ |
string title; |
135 |
+ |
vector<double> bins; |
136 |
+ |
string inputCollection; |
137 |
+ |
vector<string> inputVariables; |
138 |
+ |
} histogram; |
139 |
|
|
140 |
|
struct cut { |
141 |
|
string inputCollection; |
142 |
< |
string variable; |
143 |
< |
string comparativeOperator; |
144 |
< |
string function; |
145 |
< |
double cutValue; |
142 |
> |
vector<string> functions; |
143 |
> |
vector<string> variables; |
144 |
> |
vector<string> comparativeOperators;// >, <, =, etc. |
145 |
> |
vector<double> cutValues; |
146 |
> |
vector<string> logicalOperators;//and, or |
147 |
> |
int numSubcuts; |
148 |
|
int numberRequired; |
149 |
|
string eventComparativeOperator; |
150 |
|
string name; |
154 |
|
string name; |
155 |
|
vector<string> triggers; |
156 |
|
vector<cut> cuts; |
134 |
– |
vector<string> inputCollections; |
157 |
|
}; |
158 |
|
|
159 |
+ |
vector<string> objectsToPlot; |
160 |
|
vector<string> allNecessaryObjects; |
161 |
|
vector<channel> channels; |
162 |
+ |
vector<histogram> histograms; |
163 |
|
|
164 |
< |
template <class InputCollection> void setObjectFlags(cut &, uint, flagMap &, flagMap &, InputCollection, string); |
164 |
> |
PUWeight *puWeight_; |
165 |
|
|
166 |
+ |
template <class InputCollection> void setObjectFlags(cut &, uint, flagMap &, flagMap &, InputCollection, string); |
167 |
+ |
template <class InputCollection1, class InputCollection2> void setObjectFlags(cut &, uint, flagMap &, flagMap &, InputCollection1, InputCollection2, string); |
168 |
+ |
template <class InputCollection> void fill1DHistogram(TH1*, histogram, InputCollection, vector<bool>, double); |
169 |
+ |
template <class InputCollection1, class InputCollection2> void fill1DHistogram(TH1*, histogram, InputCollection1, InputCollection2, vector<bool>, vector<bool>, vector<bool>, double); |
170 |
+ |
template <class InputCollection> void fill2DHistogram(TH2*, histogram, InputCollection, vector<bool>, double); |
171 |
+ |
template <class InputCollection1, class InputCollection2> void fill2DHistogram(TH2*, histogram, InputCollection1, InputCollection2, vector<bool>, vector<bool>, vector<bool>, double); |
172 |
|
|
173 |
|
|
174 |
|
}; |