32 |
|
UInt_t fRunNumber; |
33 |
|
UInt_t fLumiSectionNumber; |
34 |
|
UInt_t fEventNumber; |
35 |
+ |
UInt_t fNPU; |
36 |
|
Float_t fRho; |
37 |
|
UInt_t fNVertices; |
38 |
|
UInt_t fEventTriggerBits; |
51 |
|
Float_t fEle1EnergyCorrAndSmeared; |
52 |
|
Float_t fEle1Energy; |
53 |
|
Float_t fEle1EnergyRegression; |
54 |
< |
Float_t fEle1EnergyRegressionWithTrkVarTwoPtBins; |
55 |
< |
Float_t fEle1EnergyRegressionWithTrkVar; |
56 |
< |
Float_t fEle1EnergyRegressionNoTrkVarTwoPtBins; |
56 |
< |
Float_t fEle1EnergyRegressionNoTrkVar; |
54 |
> |
Float_t fEle1EnergyRegressionV0; |
55 |
> |
Float_t fEle1EnergyRegressionV1; |
56 |
> |
Float_t fEle1EnergyRegressionV2; |
57 |
|
Int_t fEle1Charge; |
58 |
|
Float_t fEle1HZZICHEP2012IDMVA; |
59 |
|
Float_t fEle1PFIso04; |
74 |
|
Float_t fEle2EnergyCorrAndSmeared; |
75 |
|
Float_t fEle2Energy; |
76 |
|
Float_t fEle2EnergyRegression; |
77 |
< |
Float_t fEle2EnergyRegressionWithTrkVarTwoPtBins; |
78 |
< |
Float_t fEle2EnergyRegressionWithTrkVar; |
79 |
< |
Float_t fEle2EnergyRegressionNoTrkVarTwoPtBins; |
80 |
< |
Float_t fEle2EnergyRegressionNoTrkVar; |
77 |
> |
Float_t fEle2EnergyRegressionV0; |
78 |
> |
Float_t fEle2EnergyRegressionV1; |
79 |
> |
Float_t fEle2EnergyRegressionV2; |
80 |
|
Int_t fEle2Charge; |
81 |
|
Float_t fEle2HZZICHEP2012IDMVA; |
82 |
|
Float_t fEle2PFIso04; |
125 |
|
fRunNumber = 0.0; |
126 |
|
fLumiSectionNumber = 0.0; |
127 |
|
fEventNumber = 0.0; |
128 |
+ |
fNPU = 0.0; |
129 |
|
fRho = 0.0; |
130 |
|
fNVertices = 0.0; |
131 |
|
fEventTriggerBits = 0; |
140 |
|
fEle1EnergyCorrAndSmeared = 0.0; |
141 |
|
fEle1Energy = 0.0; |
142 |
|
fEle1EnergyRegression = 0.0; |
143 |
< |
fEle1EnergyRegressionWithTrkVarTwoPtBins = 0.0; |
144 |
< |
fEle1EnergyRegressionWithTrkVar = 0.0; |
145 |
< |
fEle1EnergyRegressionNoTrkVarTwoPtBins = 0.0; |
146 |
< |
fEle1EnergyRegressionNoTrkVar = 0.0; |
143 |
> |
fEle1EnergyRegressionV0 = 0.0; |
144 |
> |
fEle1EnergyRegressionV1 = 0.0; |
145 |
> |
fEle1EnergyRegressionV2 = 0.0; |
146 |
|
fEle1Charge = 0; |
147 |
|
fEle1HZZICHEP2012IDMVA = 0.0; |
148 |
|
fEle1PFIso04 = 0.0; |
161 |
|
fEle2EnergyCorrAndSmeared = 0.0; |
162 |
|
fEle2Energy = 0.0; |
163 |
|
fEle2EnergyRegression = 0.0; |
164 |
< |
fEle2EnergyRegressionWithTrkVarTwoPtBins = 0.0; |
165 |
< |
fEle2EnergyRegressionWithTrkVar = 0.0; |
166 |
< |
fEle2EnergyRegressionNoTrkVarTwoPtBins = 0.0; |
168 |
< |
fEle2EnergyRegressionNoTrkVar = 0.0; |
164 |
> |
fEle2EnergyRegressionV0 = 0.0; |
165 |
> |
fEle2EnergyRegressionV1 = 0.0; |
166 |
> |
fEle2EnergyRegressionV2 = 0.0; |
167 |
|
fEle2Charge = 0; |
168 |
|
fEle2HZZICHEP2012IDMVA = 0.0; |
169 |
|
fEle2PFIso04 = 0.0; |
202 |
|
} else { |
203 |
|
cout << "Warning: Type " << Type << " is not supported \n"; |
204 |
|
} |
205 |
+ |
InitTree(Type); |
206 |
|
assert(tree_); |
207 |
|
} |
208 |
|
|
216 |
|
tree_->Branch("run",&fRunNumber,"run/i"); |
217 |
|
tree_->Branch("lumi",&fLumiSectionNumber,"lumi/i"); |
218 |
|
tree_->Branch("event",&fEventNumber,"event/i"); |
219 |
+ |
tree_->Branch("npu",&fNPU,"npu/i"); |
220 |
|
tree_->Branch("rho",&fRho,"rho/F"); |
221 |
|
tree_->Branch("vertices",&fNVertices,"vertices/i"); |
222 |
|
tree_->Branch("triggerbits",&fEventTriggerBits,"triggerbits/i"); |
232 |
|
tree_->Branch("Ele1EnergyCorrAndSmeared",&fEle1EnergyCorrAndSmeared,"Ele1EnergyCorrAndSmeared/F"); |
233 |
|
tree_->Branch("Ele1Energy",&fEle1Energy,"Ele1Energy/F"); |
234 |
|
tree_->Branch("Ele1EnergyRegression",&fEle1EnergyRegression,"Ele1EnergyRegression/F"); |
235 |
< |
tree_->Branch("Ele1EnergyRegressionWithTrkVarTwoPtBins",&fEle1EnergyRegressionWithTrkVarTwoPtBins,"Ele1EnergyRegressionWithTrkVarTwoPtBins/F"); |
236 |
< |
tree_->Branch("Ele1EnergyRegressionWithTrkVar",&fEle1EnergyRegressionWithTrkVar,"Ele1EnergyRegressionWithTrkVar/F"); |
237 |
< |
tree_->Branch("Ele1EnergyRegressionNoTrkVarTwoPtBins",&fEle1EnergyRegressionNoTrkVarTwoPtBins,"Ele1EnergyRegressionNoTrkVarTwoPtBins/F"); |
238 |
< |
tree_->Branch("Ele1EnergyRegressionNoTrkVar",&fEle1EnergyRegressionNoTrkVar,"Ele1EnergyRegressionNoTrkVar/F"); |
235 |
> |
tree_->Branch("Ele1EnergyRegressionV0",&fEle1EnergyRegressionV0,"Ele1EnergyRegressionV0/F"); |
236 |
> |
tree_->Branch("Ele1EnergyRegressionV1",&fEle1EnergyRegressionV1,"Ele1EnergyRegressionV1/F"); |
237 |
> |
tree_->Branch("Ele1EnergyRegressionV2",&fEle1EnergyRegressionV2,"Ele1EnergyRegressionV2/F"); |
238 |
|
tree_->Branch("Ele1Charge",&fEle1Charge,"Ele1Charge/I"); |
239 |
|
tree_->Branch("Ele1HZZICHEP2012IDMVA",&fEle1HZZICHEP2012IDMVA,"Ele1HZZICHEP2012IDMVA/F"); |
240 |
|
tree_->Branch("Ele1PFIso04",&fEle1PFIso04,"Ele1PFIso04/F"); |
253 |
|
tree_->Branch("Ele2EnergyCorrAndSmeared",&fEle2EnergyCorrAndSmeared,"Ele2EnergyCorrAndSmeared/F"); |
254 |
|
tree_->Branch("Ele2Energy",&fEle2Energy,"Ele2Energy/F"); |
255 |
|
tree_->Branch("Ele2EnergyRegression",&fEle2EnergyRegression,"Ele2EnergyRegression/F"); |
256 |
< |
tree_->Branch("Ele2EnergyRegressionWithTrkVarTwoPtBins",&fEle2EnergyRegressionWithTrkVarTwoPtBins,"Ele2EnergyRegressionWithTrkVarTwoPtBins/F"); |
257 |
< |
tree_->Branch("Ele2EnergyRegressionWithTrkVar",&fEle2EnergyRegressionWithTrkVar,"Ele2EnergyRegressionWithTrkVar/F"); |
258 |
< |
tree_->Branch("Ele2EnergyRegressionNoTrkVarTwoPtBins",&fEle2EnergyRegressionNoTrkVarTwoPtBins,"Ele2EnergyRegressionNoTrkVarTwoPtBins/F"); |
260 |
< |
tree_->Branch("Ele2EnergyRegressionNoTrkVar",&fEle2EnergyRegressionNoTrkVar,"Ele2EnergyRegressionNoTrkVar/F"); |
256 |
> |
tree_->Branch("Ele2EnergyRegressionV0",&fEle2EnergyRegressionV0,"Ele2EnergyRegressionV0/F"); |
257 |
> |
tree_->Branch("Ele2EnergyRegressionV1",&fEle2EnergyRegressionV1,"Ele2EnergyRegressionV1/F"); |
258 |
> |
tree_->Branch("Ele2EnergyRegressionV2",&fEle2EnergyRegressionV2,"Ele2EnergyRegressionV2/F"); |
259 |
|
tree_->Branch("Ele2Charge",&fEle2Charge,"Ele2Charge/I"); |
260 |
|
tree_->Branch("Ele2HZZICHEP2012IDMVA",&fEle2HZZICHEP2012IDMVA,"Ele2HZZICHEP2012IDMVA/F"); |
261 |
|
tree_->Branch("Ele2PFIso04",&fEle2PFIso04,"Ele2PFIso04/F"); |
277 |
|
Int_t currentState = gErrorIgnoreLevel; |
278 |
|
|
279 |
|
if (Type == kCITZeeEvent) { |
280 |
+ |
cout << "CITZeeEvent\n"; |
281 |
|
tree_->SetBranchAddress("weight",&fWeight); |
282 |
|
tree_->SetBranchAddress("run",&fRunNumber); |
283 |
|
tree_->SetBranchAddress("lumi",&fLumiSectionNumber); |
284 |
|
tree_->SetBranchAddress("event",&fEventNumber); |
285 |
+ |
tree_->SetBranchAddress("npu",&fNPU); |
286 |
|
tree_->SetBranchAddress("rho",&fRho); |
287 |
|
tree_->SetBranchAddress("vertices",&fNVertices); |
288 |
|
tree_->SetBranchAddress("triggerbits",&fEventTriggerBits); |
297 |
|
tree_->SetBranchAddress("Ele1EnergyCorrAndSmeared",&fEle1EnergyCorrAndSmeared); |
298 |
|
tree_->SetBranchAddress("Ele1Energy",&fEle1Energy); |
299 |
|
tree_->SetBranchAddress("Ele1EnergyRegression",&fEle1EnergyRegression); |
300 |
< |
tree_->SetBranchAddress("Ele1EnergyRegressionWithTrkVarTwoPtBins",&fEle1EnergyRegressionWithTrkVarTwoPtBins); |
301 |
< |
tree_->SetBranchAddress("Ele1EnergyRegressionWithTrkVar",&fEle1EnergyRegressionWithTrkVar); |
302 |
< |
tree_->SetBranchAddress("Ele1EnergyRegressionNoTrkVarTwoPtBins",&fEle1EnergyRegressionNoTrkVarTwoPtBins); |
303 |
< |
tree_->SetBranchAddress("Ele1EnergyRegressionNoTrkVar",&fEle1EnergyRegressionNoTrkVar); |
300 |
> |
tree_->SetBranchAddress("Ele1EnergyRegressionV0",&fEle1EnergyRegressionV0); |
301 |
> |
tree_->SetBranchAddress("Ele1EnergyRegressionV1",&fEle1EnergyRegressionV1); |
302 |
> |
tree_->SetBranchAddress("Ele1EnergyRegressionV2",&fEle1EnergyRegressionV2); |
303 |
|
tree_->SetBranchAddress("Ele1Charge",&fEle1Charge); |
304 |
|
tree_->SetBranchAddress("Ele1HZZICHEP2012IDMVA",&fEle1HZZICHEP2012IDMVA); |
305 |
|
tree_->SetBranchAddress("Ele1PFIso04",&fEle1PFIso04); |
318 |
|
tree_->SetBranchAddress("Ele2EnergyCorrAndSmeared",&fEle2EnergyCorrAndSmeared); |
319 |
|
tree_->SetBranchAddress("Ele2Energy",&fEle2Energy); |
320 |
|
tree_->SetBranchAddress("Ele2EnergyRegression",&fEle2EnergyRegression); |
321 |
< |
tree_->SetBranchAddress("Ele2EnergyRegressionWithTrkVarTwoPtBins",&fEle2EnergyRegressionWithTrkVarTwoPtBins); |
322 |
< |
tree_->SetBranchAddress("Ele2EnergyRegressionWithTrkVar",&fEle2EnergyRegressionWithTrkVar); |
323 |
< |
tree_->SetBranchAddress("Ele2EnergyRegressionNoTrkVarTwoPtBins",&fEle2EnergyRegressionNoTrkVarTwoPtBins); |
325 |
< |
tree_->SetBranchAddress("Ele2EnergyRegressionNoTrkVar",&fEle2EnergyRegressionNoTrkVar); |
321 |
> |
tree_->SetBranchAddress("Ele2EnergyRegressionV0",&fEle2EnergyRegressionV0); |
322 |
> |
tree_->SetBranchAddress("Ele2EnergyRegressionV1",&fEle2EnergyRegressionV1); |
323 |
> |
tree_->SetBranchAddress("Ele2EnergyRegressionV2",&fEle2EnergyRegressionV2); |
324 |
|
tree_->SetBranchAddress("Ele2Charge",&fEle2Charge); |
325 |
|
tree_->SetBranchAddress("Ele2HZZICHEP2012IDMVA",&fEle2HZZICHEP2012IDMVA); |
326 |
|
tree_->SetBranchAddress("Ele2PFIso04",&fEle2PFIso04); |
330 |
|
tree_->SetBranchAddress("Ele2PassTightSimpleCuts",&fEle2PassTightSimpleCuts); |
331 |
|
tree_->SetBranchAddress("Ele2PassHZZICHEP2012",&fEle2PassHZZICHEP2012); |
332 |
|
} else if (Type == kHggFutyanZeeDataEvent || Type == kHggFutyanZeeMCEvent ) { |
333 |
< |
|
333 |
> |
cout << "other type\n"; |
334 |
|
//************************************************************ |
335 |
|
//The commented out branches don't exist in this version |
336 |
|
//************************************************************ |