Fork me on GitHub

Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed
  • modules/Isolation.cc

    rfdfac34 r004aa60  
    206206      {
    207207        sumAllParticles += 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           }
     208        if(isolation->Charge !=0)
     209        {
     210          sumCharged += isolationMomentum.Pt();
     211          if(isolation->IsRecoPU != 0) sumChargedPU += isolationMomentum.Pt();
    218212        }
    219213        else
     
    240234
    241235    // correct sum for pile-up contamination
    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);
     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);
    244238    ratioDBeta = sumDBeta/candidateMomentum.Pt();
    245239    ratioRhoCorr = sumRhoCorr/candidateMomentum.Pt();
Note: See TracChangeset for help on using the changeset viewer.