Fork me on GitHub

Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed
  • modules/Isolation.cc

    r004aa60 rfdfac34  
    206206      {
    207207        sumAllParticles += isolationMomentum.Pt();
    208         if(isolation->Charge !=0)
    209         {
    210           sumCharged += isolationMomentum.Pt();
    211           if(isolation->IsRecoPU != 0) sumChargedPU += isolationMomentum.Pt();
     208        if(isolation->Charge != 0)
     209        {
     210          if(isolation->IsRecoPU)
     211          {
     212            sumChargedPU += isolationMomentum.Pt();
     213          }
     214          else
     215          {
     216            sumCharged += isolationMomentum.Pt();
     217          }
    212218        }
    213219        else
     
    234240
    235241    // correct sum for pile-up contamination
    236     sumDBeta = sumCharged + TMath::Max(sumNeutral-0.5*sumChargedPU, 0.0);
    237     sumRhoCorr = sumCharged + TMath::Max(sumNeutral-TMath::Max(rho, 0.0)*fDeltaRMax*fDeltaRMax*TMath::Pi(), 0.0);
     242    sumDBeta = sumCharged + TMath::Max(sumNeutral - 0.5*sumChargedPU, 0.0);
     243    sumRhoCorr = sumCharged + TMath::Max(sumNeutral - TMath::Max(rho, 0.0)*fDeltaRMax*fDeltaRMax*TMath::Pi(), 0.0);
    238244    ratioDBeta = sumDBeta/candidateMomentum.Pt();
    239245    ratioRhoCorr = sumRhoCorr/candidateMomentum.Pt();
Note: See TracChangeset for help on using the changeset viewer.