- Timestamp:
- Jun 6, 2016, 4:35:02 PM (9 years ago)
- Branches:
- ImprovedOutputFile, Timing, dual_readout, llp, master
- Children:
- 740d430
- Parents:
- f8299bc
- Location:
- modules
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
modules/Isolation.cc
rf8299bc re2339af 110 110 111 111 fUseRhoCorrection = GetBool("UseRhoCorrection", true); 112 113 fDeltaRMin = GetDouble("DeltaRMin", 0.01); 114 fUseMiniCone = GetBool("UseMiniCone", false); 112 115 113 116 fClassifier->fPTMin = GetDouble("PTMin", 0.5); … … 157 160 Double_t sumChargedNoPU, sumChargedPU, sumNeutral, sumAllParticles; 158 161 Double_t sumDBeta, ratioDBeta, sumRhoCorr, ratioRhoCorr, sum, ratio; 159 Int_t counter;162 Bool_t pass = kFALSE; 160 163 Double_t eta = 0.0; 161 164 Double_t rho = 0.0; … … 197 200 sumAllParticles = 0.0; 198 201 199 counter = 0;200 202 itIsolationArray.Reset(); 201 202 203 while((isolation = static_cast<Candidate*>(itIsolationArray.Next()))) 203 204 { 204 205 const TLorentzVector &isolationMomentum = isolation->Momentum; 205 206 206 if(candidateMomentum.DeltaR(isolationMomentum) <= fDeltaRMax && 207 candidate->GetUniqueID() != isolation->GetUniqueID()) 208 { 207 if(fUseMiniCone) 208 { 209 pass = candidateMomentum.DeltaR(isolationMomentum) <= fDeltaRMax && 210 candidateMomentum.DeltaR(isolationMomentum) > fDeltaRMin; 211 } 212 else 213 { 214 pass = candidateMomentum.DeltaR(isolationMomentum) <= fDeltaRMax && 215 candidate->GetUniqueID() != isolation->GetUniqueID(); 216 } 217 218 if(pass) 219 { 220 209 221 sumAllParticles += isolationMomentum.Pt(); 210 222 if(isolation->Charge != 0) … … 223 235 sumNeutral += isolationMomentum.Pt(); 224 236 } 225 ++counter;226 } 237 } 238 227 239 } 228 240 229 241 // find rho 230 242 rho = 0.0; 231 243 if(fRhoInputArray) … … 240 252 } 241 253 } 254 255 242 256 243 257 // correct sum for pile-up contamination -
modules/Isolation.h
rf8299bc re2339af 57 57 Double_t fPTSumMax; 58 58 59 Double_t fDeltaRMin; 60 59 61 Bool_t fUsePTSum; 60 62 61 63 Bool_t fUseRhoCorrection; 64 65 Bool_t fUseMiniCone; 62 66 63 67 IsolationClassifier *fClassifier; //!
Note:
See TracChangeset
for help on using the changeset viewer.