Fork me on GitHub

Changeset 06a83b3 in git for examples


Ignore:
Timestamp:
Oct 7, 2016, 10:54:51 AM (8 years ago)
Author:
Michele Selvaggi <michele.selvaggi@…>
Branches:
ImprovedOutputFile, Timing, dual_readout, llp, master
Children:
386f526, cb159ed
Parents:
f3be533
Message:

fixed very high energy b and tau tagging, and light mistag collection

File:
1 edited

Legend:

Unmodified
Added
Removed
  • examples/Validation.cpp

    rf3be533 r06a83b3  
    263263        //if(Momentum.Px() == 0 && genMomentum.Py() == 0) continue;
    264264
     265
     266
    265267          // take the closest parton candidate
    266268          if(TMath::Abs(pdgID) == 5)
     
    277279            if( !(jet->BTag & (1 << 0)) ) continue;
    278280          }
    279 
    280            if(TMath::Abs(pdgID) == 1)
     281       
     282          if(TMath::Abs(pdgID) == 1)
    281283          {
    282284            Jet *jet = (Jet *)recoObj;
    283             if( !(jet->BTag & (1 << 0)) ) continue;
     285            if( !(jet->BTag & (1 << 0)) ) continue;
    284286          }
    285287
     
    13491351  TClonesArray *branchPFJet = treeReaderJet->UseBranch("PFJet");
    13501352  TClonesArray *branchCaloJet = treeReaderJet->UseBranch("CaloJet");
    1351 
     1353  TClonesArray *branchJet = treeReaderJet->UseBranch("Jet");
     1354 
    13521355  TClonesArray *branchParticleBJet = treeReaderBJet->UseBranch("Particle");
    13531356  TClonesArray *branchPFBJet = treeReaderBJet->UseBranch("Jet");
     
    25572560    for (k = 0; k < etaVals.size()-1; k++)
    25582561    {
    2559 
    2560        h_recbjet_lmis_pt = GetEffPt<Jet>(branchPFJet, branchParticleJet, "Jet", 1, ptMin, ptMax, etaVals.at(k), etaVals.at(k+1), treeReaderJet);
     2562       
     2563       h_recbjet_lmis_pt = GetEffPt<Jet>(branchJet, branchParticleJet, "Jet", 1, ptMin, ptMax, etaVals.at(k), etaVals.at(k+1), treeReaderJet);
    25612564       gr_recbjet_lmis_pt[k] = TGraphErrors(h_recbjet_lmis_pt);
    25622565
     
    25742577    for (k = 0; k < ptVals.size(); k++)
    25752578    {
    2576        h_recbjet_lmis_eta = GetEffEta<Jet>(branchPFJet, branchParticleJet, "Jet", 1, 0.5*ptVals.at(k), 2.0*ptVals.at(k) ,etaMin, etaMax , treeReaderJet);
     2579       h_recbjet_lmis_eta = GetEffEta<Jet>(branchJet, branchParticleJet, "Jet", 1, 0.5*ptVals.at(k), 2.0*ptVals.at(k) ,etaMin, etaMax , treeReaderJet);
    25772580       gr_recbjet_lmis_eta[k] = TGraphErrors(h_recbjet_lmis_eta);
    25782581
     
    26872690    {
    26882691
    2689        h_rectaujet_mis_pt = GetTauEffPt<Jet>(branchPFJet, branchParticleJet, "TauJet", 1, ptMin, ptMax, etaVals.at(k), etaVals.at(k+1), treeReaderJet);
     2692       h_rectaujet_mis_pt = GetTauEffPt<Jet>(branchJet, branchParticleJet, "TauJet", 1, ptMin, ptMax, etaVals.at(k), etaVals.at(k+1), treeReaderJet);
    26902693       gr_rectaujet_mis_pt[k] = TGraphErrors(h_rectaujet_mis_pt);
    26912694
     
    27032706    for (k = 0; k < ptVals.size(); k++)
    27042707    {
    2705        h_rectaujet_mis_eta = GetTauEffEta<Jet>(branchPFJet, branchParticleJet, "TauJet", 1, 0.5*ptVals.at(k), 2.0*ptVals.at(k) ,etaMin, etaMax , treeReaderJet);
     2708       h_rectaujet_mis_eta = GetTauEffEta<Jet>(branchJet, branchParticleJet, "TauJet", 1, 0.5*ptVals.at(k), 2.0*ptVals.at(k) ,etaMin, etaMax , treeReaderJet);
    27062709       gr_rectaujet_mis_eta[k] = TGraphErrors(h_rectaujet_mis_eta);
    27072710
Note: See TracChangeset for help on using the changeset viewer.