Fork me on GitHub

Changeset 1216 in svn for trunk/modules


Ignore:
Timestamp:
Jul 12, 2013, 12:12:45 AM (11 years ago)
Author:
Pavel Demin
Message:

calculate tau momentum without neutrinos

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/modules/TauTagging.cc

    r1156 r1216  
    180180void TauTagging::Process()
    181181{
    182   Candidate *jet, *tau;
     182  Candidate *jet, *tau, *daughter;
     183  TLorentzVector tauMomentum;
    183184  Double_t pt, eta, phi;
    184185  TObjArray *tauArray;
    185186  map< Int_t, DelphesFormula * >::iterator itEfficiencyMap;
    186187  DelphesFormula *formula;
    187   Int_t pdgCode, charge;
     188  Int_t pdgCode, charge, i;
    188189
    189190  // select taus
     
    210211    while((tau = static_cast<Candidate *>(itTauArray.Next())))
    211212    {
    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)
    213221      {
    214222        pdgCode = 15;
Note: See TracChangeset for help on using the changeset viewer.