Changeset f229d8a in git
- Timestamp:
- Jul 12, 2013, 12:12:45 AM (11 years ago)
- Branches:
- ImprovedOutputFile, Timing, dual_readout, llp, master
- Children:
- 7eccbe7
- Parents:
- d6eccdb
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
modules/TauTagging.cc
rd6eccdb rf229d8a 180 180 void TauTagging::Process() 181 181 { 182 Candidate *jet, *tau; 182 Candidate *jet, *tau, *daughter; 183 TLorentzVector tauMomentum; 183 184 Double_t pt, eta, phi; 184 185 TObjArray *tauArray; 185 186 map< Int_t, DelphesFormula * >::iterator itEfficiencyMap; 186 187 DelphesFormula *formula; 187 Int_t pdgCode, charge ;188 Int_t pdgCode, charge, i; 188 189 189 190 // select taus … … 210 211 while((tau = static_cast<Candidate *>(itTauArray.Next()))) 211 212 { 212 if(jetMomentum.DeltaR(tau->Momentum) <= fDeltaR) 213 tauMomentum.SetPxPyPzE(0.0, 0.0, 0.0, 0.0); 214 for(i = tau->D1; i <= tau->D2; ++i) 215 { 216 daughter = static_cast<Candidate *>(fParticleInputArray->At(i)); 217 if(TMath::Abs(daughter->PID) == 16) continue; 218 tauMomentum += daughter->Momentum; 219 } 220 if(jetMomentum.DeltaR(tauMomentum) <= fDeltaR) 213 221 { 214 222 pdgCode = 15;
Note:
See TracChangeset
for help on using the changeset viewer.