Changeset dd514ae in git for modules/TauTagging.cc
- Timestamp:
- Jun 10, 2015, 4:04:47 PM (9 years ago)
- Branches:
- ImprovedOutputFile, Timing, dual_readout, llp, master
- Children:
- 8707eeb
- Parents:
- 6ec8d18 (diff), 9c491e1 (diff)
Note: this is a merge changeset, the changes displayed below correspond to the merge itself.
Use the(diff)
links above to see all the changes relative to each parent. - File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
modules/TauTagging.cc
r6ec8d18 rdd514ae 73 73 74 74 if(tau->D1 < 0) return -1; 75 76 if(tau->D2 < tau->D1) return -1; 75 77 76 78 if(tau->D1 >= fParticleInputArray->GetEntriesFast() || … … 201 203 tauArray = fFilter->GetSubArray(fClassifier, 0); 202 204 203 if(tauArray == 0) return;204 205 TIter itTauArray(tauArray);206 207 205 // loop over all input jets 208 206 fItJetInputArray->Reset(); … … 217 215 218 216 // loop over all input taus 219 itTauArray.Reset(); 220 while((tau = static_cast<Candidate *>(itTauArray.Next()))) 221 { 222 if(tau->D1 < 0) continue; 223 224 if(tau->D1 >= fParticleInputArray->GetEntriesFast() || 225 tau->D2 >= fParticleInputArray->GetEntriesFast()) 217 if(tauArray){ 218 TIter itTauArray(tauArray); 219 while((tau = static_cast<Candidate *>(itTauArray.Next()))) 226 220 { 227 throw runtime_error("tau's daughter index is greater than the ParticleInputArray size"); 228 } 229 230 tauMomentum.SetPxPyPzE(0.0, 0.0, 0.0, 0.0); 231 232 for(i = tau->D1; i <= tau->D2; ++i) 233 { 234 daughter = static_cast<Candidate *>(fParticleInputArray->At(i)); 235 if(TMath::Abs(daughter->PID) == 16) continue; 236 tauMomentum += daughter->Momentum; 237 } 238 239 if(jetMomentum.DeltaR(tauMomentum) <= fDeltaR) 240 { 241 pdgCode = 15; 242 charge = tau->Charge; 221 if(tau->D1 < 0) continue; 222 223 if(tau->D1 >= fParticleInputArray->GetEntriesFast() || 224 tau->D2 >= fParticleInputArray->GetEntriesFast()) 225 { 226 throw runtime_error("tau's daughter index is greater than the ParticleInputArray size"); 227 } 228 229 tauMomentum.SetPxPyPzE(0.0, 0.0, 0.0, 0.0); 230 231 for(i = tau->D1; i <= tau->D2; ++i) 232 { 233 daughter = static_cast<Candidate *>(fParticleInputArray->At(i)); 234 if(TMath::Abs(daughter->PID) == 16) continue; 235 tauMomentum += daughter->Momentum; 236 } 237 238 if(jetMomentum.DeltaR(tauMomentum) <= fDeltaR) 239 { 240 pdgCode = 15; 241 charge = tau->Charge; 242 } 243 243 } 244 244 }
Note:
See TracChangeset
for help on using the changeset viewer.