- Timestamp:
- Nov 20, 2015, 1:43:00 PM (9 years ago)
- Branches:
- ImprovedOutputFile, Timing, dual_readout, llp, master
- Children:
- 715ab7c
- Parents:
- 004aa60
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
modules/Isolation.cc
r004aa60 rfdfac34 206 206 { 207 207 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 } 212 218 } 213 219 else … … 234 240 235 241 // 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); 238 244 ratioDBeta = sumDBeta/candidateMomentum.Pt(); 239 245 ratioRhoCorr = sumRhoCorr/candidateMomentum.Pt();
Note:
See TracChangeset
for help on using the changeset viewer.