Fork me on GitHub

Changeset a3a556c in git for examples


Ignore:
Timestamp:
Jun 28, 2016, 2:03:11 PM (8 years ago)
Author:
Alexandre Mertens <alexandre.mertens@…>
Branches:
ImprovedOutputFile, Timing, dual_readout, llp, master
Children:
f01dee6
Parents:
fe3908f
Message:

fit like CMS

File:
1 edited

Legend:

Unmodified
Added
Removed
  • examples/Validation.C

    rfe3908f ra3a556c  
    101101    obj = object;
    102102
    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);
    104104    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);
    105105
     
    365365      }
    366366
    367       if(deltaR < 0.3)
     367      if(deltaR < 0.25)
    368368      {
    369369        pt  = genJetMomentum.Pt();
     
    372372        for (bin = 0; bin < Nbins; bin++)
    373373        {
    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)
    375375            {
    376                 histos->at(bin).cenResolHist->Fill((bestGenJetMomentum.E()-jetMomentum.E())/bestGenJetMomentum.E());
     376                histos->at(bin).cenResolHist->Fill(jetMomentum.Pt()/bestGenJetMomentum.Pt());
    377377            }
    378378        }
Note: See TracChangeset for help on using the changeset viewer.