- Timestamp:
- Aug 19, 2013, 3:40:54 PM (11 years ago)
- Branches:
- ImprovedOutputFile, Timing, dual_readout, llp, master
- Children:
- 76d3973
- Parents:
- 82e07cf
- Location:
- modules
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
modules/Isolation.cc
r82e07cf rb286067 90 90 fPTRatioMax = GetDouble("PTRatioMax", 0.1); 91 91 92 fPTSumMax = GetDouble("PTSumMax", 5.0); 93 94 fUsePTSum = GetBool("UsePTSum", false); 95 92 96 fClassifier->fPTMin = GetDouble("PTMin", 0.5); 93 97 … … 132 136 Candidate *candidate, *isolation; 133 137 TObjArray *isolationArray; 134 Double_t sum PT, ratio;138 Double_t sum, ratio; 135 139 Int_t counter; 136 140 Double_t rho = 0.0; … … 157 161 158 162 // loop over all input tracks 159 sum PT= 0.0;163 sum = 0.0; 160 164 counter = 0; 161 165 itIsolationArray.Reset(); … … 167 171 !candidate->Overlaps(isolation)) 168 172 { 169 sum PT+= isolationMomentum.Pt();173 sum += isolationMomentum.Pt(); 170 174 ++counter; 171 175 } 172 176 } 173 177 174 // correct sum PTfor pile-up contamination175 sum PT = sumPT- rho*fDeltaRMax*fDeltaRMax*TMath::Pi();178 // correct sum for pile-up contamination 179 sum = sum - rho*fDeltaRMax*fDeltaRMax*TMath::Pi(); 176 180 177 ratio = sum PT/candidateMomentum.Pt();178 if( ratio > fPTRatioMax) continue;181 ratio = sum/candidateMomentum.Pt(); 182 if((fUsePTSum && sum > fPTSumMax) || ratio > fPTRatioMax) continue; 179 183 180 184 fOutputArray->Add(candidate); -
modules/Isolation.h
r82e07cf rb286067 41 41 Double_t fPTRatioMax; 42 42 43 Double_t fPTSumMax; 44 45 Bool_t fUsePTSum; 46 43 47 IsolationClassifier *fClassifier; //! 44 48
Note:
See TracChangeset
for help on using the changeset viewer.