24 |
|
kEleTrigger_Ele_CaloIdT_TrkIdVL_CaloIsoVL_TrkIsoVL = 0x000004 |
25 |
|
}; |
26 |
|
|
27 |
+ |
//******************************************* |
28 |
+ |
//=== tree versions === |
29 |
+ |
//******************************************* |
30 |
+ |
enum ElectronTreeVersion { kTreeV1, kTreeV2 }; |
31 |
+ |
|
32 |
|
/// variables |
33 |
|
Float_t fWeight; |
34 |
|
UInt_t fRunNumber; |
45 |
|
Bool_t fEleIsEcalDriven; |
46 |
|
UInt_t fEleTriggerBit; |
47 |
|
Float_t fRho; |
48 |
< |
UInt_t fNVertices; |
48 |
> |
Float_t fNVertices; |
49 |
|
|
50 |
|
// Conversion and IP |
51 |
|
Float_t fEleD0; |
55 |
|
Bool_t fEleMatchedConversion; |
56 |
|
Float_t fEleConvDCot; |
57 |
|
Float_t fEleConvDist; |
58 |
< |
UInt_t fEleNMissHits; |
58 |
> |
Float_t fEleNMissHits; |
59 |
|
|
60 |
|
// E/P variables |
61 |
|
Float_t fEleNBrem; |
85 |
|
//track quality |
86 |
|
Float_t fEleGsfTrackChi2OverNdof; |
87 |
|
Float_t fEleKFTrackChi2OverNDoF; |
88 |
< |
Int_t fEleKFTrackNHits; |
89 |
< |
Int_t fEleKFTrackNLayersWithMeasurement; |
88 |
> |
Float_t fEleKFTrackNHits; |
89 |
> |
Float_t fEleKFTrackNLayersWithMeasurement; |
90 |
|
Float_t fEleOneMinusSeedE1x5OverE5x5; |
91 |
|
|
92 |
|
//Isolation Variables |
116 |
|
|
117 |
|
Bool_t fElePassTriggerDenominator; |
118 |
|
|
119 |
+ |
//Regression Variables |
120 |
+ |
Bool_t fIsEB; |
121 |
+ |
Bool_t fIsEE; |
122 |
+ |
Float_t fSCRawEnergy; |
123 |
+ |
Float_t fNClusters; |
124 |
+ |
Float_t fEtaSeed; |
125 |
+ |
Float_t fPhiSeed; |
126 |
+ |
Float_t fESeed; |
127 |
+ |
Float_t fE3x3Seed; |
128 |
+ |
Float_t fE5x5Seed; |
129 |
+ |
Float_t fEMaxSeed; |
130 |
+ |
Float_t fE2ndSeed; |
131 |
+ |
Float_t fETopSeed; |
132 |
+ |
Float_t fEBottomSeed; |
133 |
+ |
Float_t fELeftSeed; |
134 |
+ |
Float_t fERightSeed; |
135 |
+ |
Float_t fE2x5MaxSeed; |
136 |
+ |
Float_t fE2x5TopSeed; |
137 |
+ |
Float_t fE2x5BottomSeed; |
138 |
+ |
Float_t fE2x5LeftSeed; |
139 |
+ |
Float_t fE2x5RightSeed; |
140 |
+ |
Float_t fIEtaSeed; |
141 |
+ |
Float_t fIPhiSeed; |
142 |
+ |
Float_t fEtaCrySeed; |
143 |
+ |
Float_t fPhiCrySeed; |
144 |
+ |
Float_t fEcalEnergyError; |
145 |
+ |
Float_t fGsfTrackPIn; |
146 |
+ |
Float_t fTrackMomentumError; |
147 |
+ |
Float_t fpmeangsf; |
148 |
+ |
Float_t fpmeankf; |
149 |
+ |
Float_t fCharge; |
150 |
+ |
Float_t fGeneratedEnergy; |
151 |
+ |
Float_t fGeneratedEnergyStatus1; |
152 |
+ |
Float_t fGeneratedEnergyStatus3; |
153 |
+ |
|
154 |
|
public: |
155 |
|
/// this is the main element |
156 |
|
TTree *tree_; |
174 |
|
fEventNumber = 0.0; |
175 |
|
fEleEventNumberParity = 0.0; |
176 |
|
fElePt = 0.0; |
177 |
< |
fEleEta = 0.0; |
178 |
< |
fElePhi = 0.0; |
179 |
< |
fEleSCEt = 0.0; |
177 |
> |
fEleEta = 0.0; |
178 |
> |
fElePhi = 0.0; |
179 |
> |
fEleSCEt = 0.0; |
180 |
|
fEleSCEta = 0.0; |
181 |
|
fEleSCPhi = 0.0; |
182 |
|
fEleEcalEnergy = 0.0; |
183 |
< |
fEleIsEcalDriven = 0.0; |
183 |
> |
fEleIsEcalDriven = 0.0; |
184 |
|
fEleTriggerBit = 0.0; |
185 |
|
fRho = 0.0; |
186 |
|
fNVertices = 0.0; |
187 |
|
fEleD0 = 0.0; |
188 |
|
fEleDZ = 0.0; |
189 |
< |
fEleIP3d = 0.0; |
189 |
> |
fEleIP3d = 0.0; |
190 |
|
fEleIP3dSig = 0.0; |
191 |
|
fEleMatchedConversion = 0.0; |
192 |
|
fEleConvDCot = 0.0; |
195 |
|
fEleNBrem = 0.0; |
196 |
|
fEleFBrem = 0.0; |
197 |
|
fEleEOverP = 0.0; |
198 |
< |
fEleESeedClusterOverPIn = 0.0; |
199 |
< |
fEleESeedClusterOverPout = 0.0; |
200 |
< |
fEleEEleClusterOverPout = 0.0; |
201 |
< |
fEleOneOverEMinusOneOverP = 0.0; |
198 |
> |
fEleESeedClusterOverPIn = 0.0; |
199 |
> |
fEleESeedClusterOverPout = 0.0; |
200 |
> |
fEleEEleClusterOverPout = 0.0; |
201 |
> |
fEleOneOverEMinusOneOverP = 0.0; |
202 |
|
fEleDEtaIn = 0.0; |
203 |
|
fEleDPhiIn = 0.0; |
204 |
|
fEledEtaCalo = 0.0; |
205 |
|
fEledPhiCalo = 0.0; |
206 |
|
fEleSigmaIEtaIEta = 0.0; |
207 |
|
fEleSigmaIPhiIPhi = 0.0; |
208 |
< |
fEleSigmaIEtaIPhi = 0.0; |
208 |
> |
fEleSigmaIEtaIPhi = 0.0; |
209 |
|
fEleSCEtaWidth = 0.0; |
210 |
|
fEleSCPhiWidth = 0.0; |
211 |
|
fEleR9 = 0.0; |
212 |
|
fElePreShowerOverRaw = 0.0; |
213 |
|
fEleHoverE = 0.0; |
214 |
< |
fEleGsfTrackChi2OverNdof = 0.0; |
214 |
> |
fEleGsfTrackChi2OverNdof = 0.0; |
215 |
|
fEleKFTrackChi2OverNDoF = 0.0; |
216 |
< |
fEleKFTrackNHits = 0.0; |
216 |
> |
fEleKFTrackNHits = 0.0; |
217 |
|
fEleKFTrackNLayersWithMeasurement = 0.0; |
218 |
< |
fEleOneMinusSeedE1x5OverE5x5 = 0.0; |
219 |
< |
fElePFMVA = 0.0; |
218 |
> |
fEleOneMinusSeedE1x5OverE5x5 = 0.0; |
219 |
> |
fElePFMVA = 0.0; |
220 |
|
fEleTrkIso03 = 0.0; |
221 |
|
fEleEMIso03 = 0.0; |
222 |
|
fEleHadIso03 = 0.0; |
223 |
|
fEleTrkIso04 = 0.0; |
224 |
|
fEleEMIso04 = 0.0; |
225 |
< |
fEleHadIso04 = 0.0; |
225 |
> |
fEleHadIso04 = 0.0; |
226 |
|
fElePFIso04 = 0.0; |
227 |
|
fChargedIso_DR0p0To0p1 = 0.0; |
228 |
|
fChargedIso_DR0p1To0p2 = 0.0; |
234 |
|
fGammaIso_DR0p2To0p3 = 0.0; |
235 |
|
fGammaIso_DR0p3To0p4 = 0.0; |
236 |
|
fGammaIso_DR0p4To0p5 = 0.0; |
237 |
< |
fNeutralHadronIso_DR0p0To0p1 = 0.0; |
238 |
< |
fNeutralHadronIso_DR0p1To0p2 = 0.0; |
239 |
< |
fNeutralHadronIso_DR0p2To0p3 = 0.0; |
240 |
< |
fNeutralHadronIso_DR0p3To0p4 = 0.0; |
241 |
< |
fNeutralHadronIso_DR0p4To0p5 = 0.0; |
242 |
< |
fElePassTriggerDenominator = 0.0; |
237 |
> |
fNeutralHadronIso_DR0p0To0p1 = 0.0; |
238 |
> |
fNeutralHadronIso_DR0p1To0p2 = 0.0; |
239 |
> |
fNeutralHadronIso_DR0p2To0p3 = 0.0; |
240 |
> |
fNeutralHadronIso_DR0p3To0p4 = 0.0; |
241 |
> |
fNeutralHadronIso_DR0p4To0p5 = 0.0; |
242 |
> |
fElePassTriggerDenominator = 0.0; |
243 |
> |
fIsEB = 0.0; |
244 |
> |
fIsEE = 0.0; |
245 |
> |
fSCRawEnergy = 0.0; |
246 |
> |
fNClusters = 0.0; |
247 |
> |
fEtaSeed = 0.0; |
248 |
> |
fPhiSeed = 0.0; |
249 |
> |
fESeed = 0.0; |
250 |
> |
fE3x3Seed = 0.0; |
251 |
> |
fE5x5Seed = 0.0; |
252 |
> |
fEMaxSeed = 0.0; |
253 |
> |
fE2ndSeed = 0.0; |
254 |
> |
fETopSeed = 0.0; |
255 |
> |
fEBottomSeed = 0.0; |
256 |
> |
fELeftSeed = 0.0; |
257 |
> |
fERightSeed = 0.0; |
258 |
> |
fE2x5MaxSeed = 0.0; |
259 |
> |
fE2x5TopSeed = 0.0; |
260 |
> |
fE2x5BottomSeed = 0.0; |
261 |
> |
fE2x5LeftSeed = 0.0; |
262 |
> |
fE2x5RightSeed = 0.0; |
263 |
> |
fIEtaSeed = 0.0; |
264 |
> |
fIPhiSeed = 0.0; |
265 |
> |
fEtaCrySeed = 0.0; |
266 |
> |
fPhiCrySeed = 0.0; |
267 |
> |
fEcalEnergyError = 0.0; |
268 |
> |
fGsfTrackPIn = 0.0; |
269 |
> |
fTrackMomentumError = 0.0; |
270 |
> |
fpmeangsf = 0.0; |
271 |
> |
fpmeankf = 0.0; |
272 |
> |
fCharge = 0.0; |
273 |
> |
fGeneratedEnergy = 0.0; |
274 |
> |
fGeneratedEnergyStatus1 = 0.0; |
275 |
> |
fGeneratedEnergyStatus3 = 0.0; |
276 |
> |
|
277 |
> |
|
278 |
|
} |
279 |
|
|
280 |
|
/// load a ElectronTree |
281 |
< |
void LoadTree(const char* file){ |
281 |
> |
void LoadTree(const char* file, int version=kTreeV1){ |
282 |
|
f_ = TFile::Open(file); |
283 |
|
assert(f_); |
284 |
< |
tree_ = dynamic_cast<TTree*>(f_->Get("Electrons")); |
284 |
> |
if(version==kTreeV1) tree_ = dynamic_cast<TTree*>(f_->Get("Electrons")); |
285 |
> |
else if(version==kTreeV2) tree_ = dynamic_cast<TTree*>(f_->Get("eleIDdir/T1")); |
286 |
|
assert(tree_); |
287 |
|
} |
288 |
|
|
307 |
|
tree_->Branch("ecaldriven",&fEleIsEcalDriven,"ecaldriven/O"); |
308 |
|
tree_->Branch("triggerBit",&fEleTriggerBit,"triggerBit/i"); |
309 |
|
tree_->Branch("rho",&fRho,"rho/F"); |
310 |
< |
tree_->Branch("vertices",&fNVertices,"vertices/i"); |
310 |
> |
tree_->Branch("vertices",&fNVertices,"vertices/F"); |
311 |
|
tree_->Branch("d0",&fEleD0,"d0/F"); |
312 |
|
tree_->Branch("dz",&fEleDZ,"dz/F"); |
313 |
|
tree_->Branch("ip3d",&fEleIP3d,"ip3d/F"); |
315 |
|
tree_->Branch("matchConv",&fEleMatchedConversion,"matchConv/O"); |
316 |
|
tree_->Branch("dcot",&fEleConvDCot,"dcot/F"); |
317 |
|
tree_->Branch("dist",&fEleConvDist,"dist/F"); |
318 |
< |
tree_->Branch("missHits",&fEleNMissHits,"NMissHits/i"); |
318 |
> |
tree_->Branch("missHits",&fEleNMissHits,"NMissHits/F"); |
319 |
|
tree_->Branch("nbrems",&fEleNBrem,"nbrems/F"); |
320 |
|
tree_->Branch("fbrem",&fEleFBrem,"fbrem/F"); |
321 |
|
tree_->Branch("EoP",&fEleEOverP,"EoP/F"); |
322 |
|
tree_->Branch("EoPin",&fEleESeedClusterOverPIn,"EoPin/F"); |
323 |
|
tree_->Branch("ESeedoPout",&fEleESeedClusterOverPout,"ESeedoPout/F"); |
324 |
|
tree_->Branch("EEleoPout",&fEleEEleClusterOverPout,"EEleoPout/F"); |
325 |
< |
// tree_->Branch("eleEoPout",&fEleEEleClusterOverPout,"eleEoPout/F"); |
326 |
< |
tree_->Branch("IoEmIoP",&fEleOneOverEMinusOneOverP,"IoEmIoP/F"); |
325 |
> |
// tree_->Branch("eleEoPout",&fEleEEleClusterOverPout,"eleEoPout/F"); //old name of the branch |
326 |
> |
tree_->Branch("IoEmIoP",&fEleOneOverEMinusOneOverP,"IoEmIoP/F"); |
327 |
|
tree_->Branch("deta",&fEleDEtaIn,"deta/F"); |
328 |
|
tree_->Branch("dphi",&fEleDPhiIn,"dphi/F"); |
329 |
|
tree_->Branch("detacalo",&fEledEtaCalo,"detacalo/F"); |
338 |
|
tree_->Branch("HoE",&fEleHoverE,"HoE/F"); |
339 |
|
tree_->Branch("gsfchi2",&fEleGsfTrackChi2OverNdof,"gsfchi2/F"); |
340 |
|
tree_->Branch("kfchi2",&fEleKFTrackChi2OverNDoF,"kfchi2/F"); |
341 |
< |
tree_->Branch("kfhits",&fEleKFTrackNHits,"kfhits/I"); |
342 |
< |
tree_->Branch("kflayers",&fEleKFTrackNLayersWithMeasurement,"kflayers/I"); |
341 |
> |
tree_->Branch("kfhits",&fEleKFTrackNHits,"kfhits/F"); |
342 |
> |
tree_->Branch("kflayers",&fEleKFTrackNLayersWithMeasurement,"kflayers/F"); |
343 |
|
tree_->Branch("OneMinusSeedE1x5OverE5x5",&fEleOneMinusSeedE1x5OverE5x5,"OneMinusSeedE1x5OverE5x5/F"); |
344 |
|
tree_->Branch("PFMVA",&fElePFMVA,"PFMVA/F"); |
345 |
|
tree_->Branch("trkIso03",&fEleTrkIso03,"trkIso03/F"); |
365 |
|
tree_->Branch("NeutralHadronIso_DR0p3To0p4",&fNeutralHadronIso_DR0p3To0p4,"NeutralHadronIso_DR0p3To0p4/F"); |
366 |
|
tree_->Branch("NeutralHadronIso_DR0p4To0p5",&fNeutralHadronIso_DR0p4To0p5,"NeutralHadronIso_DR0p4To0p5/F"); |
367 |
|
tree_->Branch("PassTriggerDenominator",&fElePassTriggerDenominator,"PassTriggerDenominator/O"); |
368 |
+ |
tree_->Branch("IsEB",&fIsEB,"IsEB/O"); |
369 |
+ |
tree_->Branch("IsEE",&fIsEE,"IsEE/O"); |
370 |
+ |
tree_->Branch("SCRawEnergy",&fSCRawEnergy,"SCRawEnergy/F"); |
371 |
+ |
tree_->Branch("NClusters",&fNClusters,"NClusters/F"); |
372 |
+ |
tree_->Branch("EtaSeed",&fEtaSeed,"EtaSeed/F"); |
373 |
+ |
tree_->Branch("PhiSeed",&fPhiSeed,"PhiSeed/F"); |
374 |
+ |
tree_->Branch("ESeed",&fESeed,"ESeed/F"); |
375 |
+ |
tree_->Branch("E3x3Seed",&fE3x3Seed,"E3x3Seed/F"); |
376 |
+ |
tree_->Branch("E5x5Seed",&fE5x5Seed,"E5x5Seed/F"); |
377 |
+ |
tree_->Branch("EMaxSeed",&fEMaxSeed,"EMaxSeed/F"); |
378 |
+ |
tree_->Branch("E2ndSeed",&fE2ndSeed,"E2ndSeed/F"); |
379 |
+ |
tree_->Branch("ETopSeed",&fETopSeed,"ETopSeed/F"); |
380 |
+ |
tree_->Branch("EBottomSeed",&fEBottomSeed,"EBottomSeed/F"); |
381 |
+ |
tree_->Branch("ELeftSeed",&fELeftSeed,"ELeftSeed/F"); |
382 |
+ |
tree_->Branch("ERightSeed",&fERightSeed,"ERightSeed/F"); |
383 |
+ |
tree_->Branch("E2x5MaxSeed",&fE2x5MaxSeed,"E2x5MaxSeed/F"); |
384 |
+ |
tree_->Branch("E2x5TopSeed",&fE2x5TopSeed,"E2x5TopSeed/F"); |
385 |
+ |
tree_->Branch("E2x5BottomSeed",&fE2x5BottomSeed,"E2x5BottomSeed/F"); |
386 |
+ |
tree_->Branch("E2x5LeftSeed",&fE2x5LeftSeed,"E2x5LeftSeed/F"); |
387 |
+ |
tree_->Branch("E2x5RightSeed",&fE2x5RightSeed,"E2x5RightSeed/F"); |
388 |
+ |
tree_->Branch("IEtaSeed",&fIEtaSeed,"IEtaSeed/F"); |
389 |
+ |
tree_->Branch("IPhiSeed",&fIPhiSeed,"IPhiSeed/F"); |
390 |
+ |
tree_->Branch("EtaCrySeed",&fEtaCrySeed,"EtaCrySeed/F"); |
391 |
+ |
tree_->Branch("PhiCrySeed",&fPhiCrySeed,"PhiCrySeed/F"); |
392 |
+ |
tree_->Branch("ecalenergyerror",&fEcalEnergyError,"ecalenergyerror/F"); |
393 |
+ |
tree_->Branch("pmodegsf",&fGsfTrackPIn,"pmodegsf/F"); |
394 |
+ |
tree_->Branch("perror",&fTrackMomentumError,"perror/F"); |
395 |
+ |
tree_->Branch("pmeangsf",&fpmeangsf,"pmodegsf/F"); |
396 |
+ |
tree_->Branch("pmeankf",&fpmeankf,"pmeankf/F"); |
397 |
+ |
tree_->Branch("Charge",&fCharge,"Charge/F"); |
398 |
+ |
tree_->Branch("GeneratedEnergy",&fGeneratedEnergy,"GeneratedEnergy/F"); |
399 |
+ |
tree_->Branch("GeneratedEnergyStatus1",&fGeneratedEnergyStatus1,"GeneratedEnergyStatus1/F"); |
400 |
+ |
tree_->Branch("GeneratedEnergyStatus3",&fGeneratedEnergyStatus3,"GeneratedEnergyStatus3/F"); |
401 |
|
|
402 |
|
|
403 |
|
} |
483 |
|
tree_->SetBranchAddress("NeutralHadronIso_DR0p3To0p4",&fNeutralHadronIso_DR0p3To0p4); |
484 |
|
tree_->SetBranchAddress("NeutralHadronIso_DR0p4To0p5",&fNeutralHadronIso_DR0p4To0p5); |
485 |
|
tree_->SetBranchAddress("PassTriggerDenominator",&fElePassTriggerDenominator); |
486 |
+ |
tree_->SetBranchAddress("IsEB",&fIsEB); |
487 |
+ |
tree_->SetBranchAddress("IsEE",&fIsEE); |
488 |
+ |
tree_->SetBranchAddress("SCRawEnergy",&fSCRawEnergy); |
489 |
+ |
tree_->SetBranchAddress("NClusters",&fNClusters); |
490 |
+ |
tree_->SetBranchAddress("EtaSeed",&fEtaSeed); |
491 |
+ |
tree_->SetBranchAddress("PhiSeed",&fPhiSeed); |
492 |
+ |
tree_->SetBranchAddress("ESeed",&fESeed); |
493 |
+ |
tree_->SetBranchAddress("E3x3Seed",&fE3x3Seed); |
494 |
+ |
tree_->SetBranchAddress("E5x5Seed",&fE5x5Seed); |
495 |
+ |
tree_->SetBranchAddress("EMaxSeed",&fEMaxSeed); |
496 |
+ |
tree_->SetBranchAddress("E2ndSeed",&fE2ndSeed); |
497 |
+ |
tree_->SetBranchAddress("ETopSeed",&fETopSeed); |
498 |
+ |
tree_->SetBranchAddress("EBottomSeed",&fEBottomSeed); |
499 |
+ |
tree_->SetBranchAddress("ELeftSeed",&fELeftSeed); |
500 |
+ |
tree_->SetBranchAddress("ERightSeed",&fERightSeed); |
501 |
+ |
tree_->SetBranchAddress("E2x5MaxSeed",&fE2x5MaxSeed); |
502 |
+ |
tree_->SetBranchAddress("E2x5TopSeed",&fE2x5TopSeed); |
503 |
+ |
tree_->SetBranchAddress("E2x5BottomSeed",&fE2x5BottomSeed); |
504 |
+ |
tree_->SetBranchAddress("E2x5LeftSeed",&fE2x5LeftSeed); |
505 |
+ |
tree_->SetBranchAddress("E2x5RightSeed",&fE2x5RightSeed); |
506 |
+ |
tree_->SetBranchAddress("IEtaSeed",&fIEtaSeed); |
507 |
+ |
tree_->SetBranchAddress("IPhiSeed",&fIPhiSeed); |
508 |
+ |
tree_->SetBranchAddress("EtaCrySeed",&fEtaCrySeed); |
509 |
+ |
tree_->SetBranchAddress("PhiCrySeed",&fPhiCrySeed); |
510 |
+ |
tree_->SetBranchAddress("ecalenergyerror",&fEcalEnergyError); |
511 |
+ |
tree_->SetBranchAddress("pmodegsf",&fGsfTrackPIn); |
512 |
+ |
tree_->SetBranchAddress("perror",&fTrackMomentumError); |
513 |
+ |
tree_->SetBranchAddress("pmeangsf",&fpmeankf); |
514 |
+ |
tree_->SetBranchAddress("pmeankf",&fpmeankf); |
515 |
+ |
tree_->SetBranchAddress("Charge",&fCharge); |
516 |
+ |
tree_->SetBranchAddress("GeneratedEnergy",&fGeneratedEnergy); |
517 |
+ |
tree_->SetBranchAddress("GeneratedEnergyStatus1",&fGeneratedEnergyStatus1); |
518 |
+ |
tree_->SetBranchAddress("GeneratedEnergyStatus3",&fGeneratedEnergyStatus3); |
519 |
|
|
520 |
|
gErrorIgnoreLevel = currentState; |
521 |
|
} |