Fork me on GitHub

Changes in / [0203656:e4cc3f7] in git


Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed
  • modules/StatusPidFilter.cc

    r0203656 re4cc3f7  
    104104      return false;
    105105  }
    106 
    107   bool isTauDaughter(int pdgCode, int M1, const TObjArray *fInputArray){
    108     //not needed, just to speed up the code - can be further refined but gives only negligible improvement:
    109     if ( pdgCode==15 || pdgCode<11 || (pdgCode > 22 && pdgCode < 100) || pdgCode>1000 )
    110       return false;
    111 
    112     if ( M1 < 0 )
    113       return false;
    114 
    115     Candidate *mother;
    116     mother = static_cast<Candidate*>(fInputArray->At( M1 ));
    117     if ( TMath::Abs(mother->PID) == 15 )
    118       return true;
    119 
    120     return false;
    121   }
    122 
    123106}
    124107
     
    194177    bool is_b_quark  = (pdgCode == 5);
    195178
    196     bool is_tau_daughter = isTauDaughter(pdgCode, candidate->M1, fInputArray);
    197 
    198179    if (is_b_hadron)
    199180      pass = kTRUE;
    200181
    201     if (is_tau_daughter)
    202       pass = kTRUE;
    203 
    204182    // fPTMin not applied to b_hadrons / b_quarks to allow for b-enriched sample stitching
    205     // fPTMin not applied to tau decay products to allow visible-tau four momentum determination
    206     if(!pass || (candidate->Momentum.Pt() < fPTMin && !(is_b_hadron || is_b_quark || is_tau_daughter)) ) continue;
     183    if(!pass || (candidate->Momentum.Pt() < fPTMin && !(is_b_hadron || is_b_quark)) ) continue;
    207184
    208185    fOutputArray->Add(candidate);
Note: See TracChangeset for help on using the changeset viewer.