45 |
|
|
46 |
|
#include "OSUT3Analysis/AnaTools/interface/CutFlow.h" |
47 |
|
#include "OSUT3Analysis/AnaTools/interface/PUWeight.h" |
48 |
+ |
#include "OSUT3Analysis/AnaTools/interface/SFWeight.h" |
49 |
|
|
50 |
|
|
51 |
|
|
69 |
|
double valueLookup (const BNelectron* object, string variable, string function = ""); |
70 |
|
double valueLookup (const BNelectron* object1, const BNelectron* object2, string variable, string function = ""); |
71 |
|
double valueLookup (const BNelectron* object1, const BNmuon* object2, string variable, string function = ""); |
72 |
+ |
double valueLookup (const BNelectron* object1, const BNtrack* object2, string variable, string function = ""); |
73 |
+ |
double valueLookup (const BNmuon* object1, const BNtrack* object2, string variable, string function = ""); |
74 |
+ |
double valueLookup (const BNmuon* object1, const BNtau* object2, string variable, string function = ""); |
75 |
+ |
double valueLookup (const BNtau* object1, const BNtau* object2, string variable, string function = ""); |
76 |
|
double valueLookup (const BNevent* object, string variable, string function = ""); |
77 |
|
double valueLookup (const BNtau* object, string variable, string function = ""); |
78 |
|
double valueLookup (const BNmet* object, string variable, string function = ""); |
105 |
|
//inner vector corresponds to each object in input collection |
106 |
|
//bool tells if object passes or fails cuts |
107 |
|
typedef map<string, vector < vector<bool> > > flagMap; |
108 |
< |
|
108 |
> |
|
109 |
|
//counterMap: |
110 |
|
//string holds input collection type |
111 |
|
//vector corresponds to each cut |
129 |
|
edm::InputTag triggers_; |
130 |
|
std::string puFile_; |
131 |
|
std::string deadEcalFile_; |
132 |
+ |
std::string muonSFFile_; |
133 |
|
std::string dataPU_; |
134 |
+ |
std::string electronSFID_; |
135 |
+ |
std::string muonSF_; |
136 |
|
std::string dataset_; |
137 |
|
std::string datasetType_; |
138 |
|
vector<edm::ParameterSet> channels_; |
139 |
|
vector<edm::ParameterSet> histogramSets_; |
140 |
|
bool plotAllObjectsInPassingEvents_; |
141 |
+ |
bool doPileupReweighting_; |
142 |
+ |
bool printEventInfo_; |
143 |
|
|
144 |
|
struct DeadEcal { |
145 |
< |
double etaEcal; |
146 |
< |
double phiEcal; |
145 |
> |
double etaEcal; |
146 |
> |
double phiEcal; |
147 |
|
}; |
148 |
< |
|
148 |
> |
|
149 |
|
std::vector<DeadEcal> DeadEcalVec; |
150 |
|
|
151 |
|
|
165 |
|
edm::Handle<BNphotonCollection> photons; |
166 |
|
edm::Handle<BNsuperclusterCollection> superclusters; |
167 |
|
|
168 |
< |
const BNprimaryvertex *chosenPrimaryVertex; |
159 |
< |
|
168 |
> |
flagMap cumulativeFlags; |
169 |
|
|
170 |
|
vector<map<string, TH1D*> > oneDHists_; |
171 |
|
vector<map<string, TH2D*> > twoDHists_; |
176 |
|
vector<CutFlow *> cutFlows_; |
177 |
|
|
178 |
|
typedef struct { |
179 |
< |
string name; |
180 |
< |
string title; |
181 |
< |
vector<double> bins; |
182 |
< |
string inputCollection; |
183 |
< |
vector<string> inputVariables; |
179 |
> |
string name; |
180 |
> |
string title; |
181 |
> |
vector<double> bins; |
182 |
> |
string inputCollection; |
183 |
> |
vector<string> inputVariables; |
184 |
|
} histogram; |
185 |
|
|
186 |
|
struct cut { |
187 |
< |
string inputCollection; |
188 |
< |
vector<string> functions; |
189 |
< |
vector<string> variables; |
190 |
< |
vector<string> comparativeOperators;// >, <, =, etc. |
191 |
< |
vector<double> cutValues; |
192 |
< |
vector<string> logicalOperators;//and, or |
193 |
< |
int numSubcuts; |
194 |
< |
int numberRequired; |
195 |
< |
string eventComparativeOperator; |
196 |
< |
string name; |
187 |
> |
string inputCollection; |
188 |
> |
vector<string> functions; |
189 |
> |
vector<string> variables; |
190 |
> |
vector<string> comparativeOperators;// >, <, =, etc. |
191 |
> |
vector<double> cutValues; |
192 |
> |
vector<string> logicalOperators;//and, or |
193 |
> |
int numSubcuts; |
194 |
> |
int numberRequired; |
195 |
> |
string eventComparativeOperator; |
196 |
> |
string name; |
197 |
|
}; |
198 |
|
|
199 |
|
struct channel { |
200 |
< |
string name; |
201 |
< |
vector<string> triggers; |
202 |
< |
vector<cut> cuts; |
200 |
> |
string name; |
201 |
> |
vector<string> triggers; |
202 |
> |
vector<cut> cuts; |
203 |
|
}; |
204 |
|
|
205 |
|
vector<string> objectsToGet; |
209 |
|
vector<histogram> histograms; |
210 |
|
|
211 |
|
PUWeight *puWeight_; |
212 |
+ |
MuonSFWeight *muonSFWeight_; |
213 |
+ |
ElectronSFWeight *electronSFWeight_; |
214 |
|
|
215 |
|
template <class InputCollection> void setObjectFlags(cut &, uint, flagMap &, flagMap &, InputCollection, string); |
216 |
|
template <class InputCollection1, class InputCollection2> void setObjectFlags(cut &, uint, flagMap &, flagMap &, InputCollection1, InputCollection2, vector<bool>, vector<bool>, string); |
224 |
|
int getPdgIdBinValue(int); |
225 |
|
int findTauMotherIndex(const BNmcparticle*); |
226 |
|
|
227 |
+ |
template <class InputObject> double getGenDeltaRLowest(InputObject); |
228 |
+ |
|
229 |
+ |
|
230 |
+ |
const BNprimaryvertex *chosenVertex (); |
231 |
+ |
const BNmet *chosenMET (); |
232 |
+ |
const BNelectron *chosenElectron (); |
233 |
+ |
const BNmuon *chosenMuon (); |
234 |
|
|
235 |
|
}; |
236 |
|
|