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 |
|
} |
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 |
|
|