- Timestamp:
- Jun 28, 2016, 2:03:11 PM (8 years ago)
- Branches:
- ImprovedOutputFile, Timing, dual_readout, llp, master
- Children:
- f01dee6
- Parents:
- fe3908f
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
examples/Validation.C
rfe3908f ra3a556c 101 101 obj = object; 102 102 103 cenResolHist = new TH1D(obj+"_delta_pt_"+Form("%4.2f",ptmin)+"_"+Form("%4.2f",ptmax)+"_cen", obj+"_delta_pt_"+Form("%4.2f",ptmin)+"_"+Form("%4.2f",ptmax)+"_cen", 500, -1, 1);103 cenResolHist = new TH1D(obj+"_delta_pt_"+Form("%4.2f",ptmin)+"_"+Form("%4.2f",ptmax)+"_cen", obj+"_delta_pt_"+Form("%4.2f",ptmin)+"_"+Form("%4.2f",ptmax)+"_cen", 500, 0, 2); 104 104 fwdResolHist = new TH1D(obj+"_delta_pt_"+Form("%4.2f",ptmin)+"_"+Form("%4.2f",ptmax)+"_fwd", obj+"_delta_pt_"+Form("%4.2f",ptmin)+"_"+Form("%4.2f",ptmax)+"_fwd", 500, 0.4, 0.4); 105 105 … … 365 365 } 366 366 367 if(deltaR < 0. 3)367 if(deltaR < 0.25) 368 368 { 369 369 pt = genJetMomentum.Pt(); … … 372 372 for (bin = 0; bin < Nbins; bin++) 373 373 { 374 if(pt > histos->at(bin).ptmin && pt < histos->at(bin).ptmax && eta > 0.0 && eta < 2.5)374 if(pt > histos->at(bin).ptmin && pt < histos->at(bin).ptmax && eta > 0.0 && eta < 1.5) 375 375 { 376 histos->at(bin).cenResolHist->Fill( (bestGenJetMomentum.E()-jetMomentum.E())/bestGenJetMomentum.E());376 histos->at(bin).cenResolHist->Fill(jetMomentum.Pt()/bestGenJetMomentum.Pt()); 377 377 } 378 378 }
Note:
See TracChangeset
for help on using the changeset viewer.