ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/MitPhysics/Utils/src/JetTools.cc
(Generate patch)

Comparing UserCode/MitPhysics/Utils/src/JetTools.cc (file contents):
Revision 1.24 by ceballos, Thu Apr 26 11:14:16 2012 UTC vs.
Revision 1.28 by pharris, Tue May 1 16:53:38 2012 UTC

# Line 517 | Line 517 | Double_t JetTools::frac(const PFJet *iJe
517      Double_t pDR = MathUtils::DeltaR(iJet->Mom(),pCand->Mom());
518      if(pDR > iDRMax) continue;
519      if(pDR < iDRMax-0.1) continue;
520 <    lFrac += pCand->Pt()/iJet->Pt();
520 >    lFrac += pCand->Pt()/iJet->RawMom().Pt();
521    }
522    return lFrac;
523   }
# Line 544 | Line 544 | Double_t JetTools::betaStar(const PFJet
544    return lPileup/(lTotal);
545   }
546   Bool_t  JetTools::passPFLooseId(const PFJet *iJet) {
547 <  if(iJet->E()                              == 0)       return false;
548 <  if(iJet->NeutralHadronEnergy()/iJet->E()  >  0.99)    return false;
549 <  if(iJet->NeutralEmEnergy()/iJet->E()      >  0.99)    return false;
550 <  if(iJet->NConstituents()                  <  2)       return false;
551 <  if(iJet->ChargedHadronEnergy()/iJet->E()  <= 0     && fabs(iJet->Eta()) < 2.4 ) return false;
552 <  if(iJet->ChargedEmEnergy()/iJet->E()      >  0.99  && fabs(iJet->Eta()) < 2.4 ) return false;
553 <  if(iJet->ChargedMultiplicity()            < 1      && fabs(iJet->Eta()) < 2.4 ) return false;
547 >  if(iJet->RawMom().E()                              == 0)       return false;
548 >  if(iJet->NeutralHadronEnergy()/iJet->RawMom().E()  >  0.99)    return false;
549 >  if(iJet->NeutralEmEnergy()/iJet->RawMom().E()      >  0.99)    return false;
550 >  if(iJet->NConstituents()                           <  2)       return false;
551 >  if(iJet->ChargedHadronEnergy()/iJet->RawMom().E()  <= 0     && fabs(iJet->Eta()) < 2.4 ) return false;
552 >  if(iJet->ChargedEmEnergy()/iJet->RawMom().E()      >  0.99  && fabs(iJet->Eta()) < 2.4 ) return false;
553 >  if(iJet->ChargedMultiplicity()                     < 1      && fabs(iJet->Eta()) < 2.4 ) return false;
554 >  if(fabs(iJet->Eta())                               > 4.99) return false;
555    return true;
556   }
557  

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines