Changeset 1216 in svn for trunk

Timestamp:

Jul 12, 2013, 12:12:45 AM (11 years ago)

Author:

Pavel Demin

Message:

calculate tau momentum without neutrinos

File:

• trunk/modules/TauTagging.cc (modified) (2 diffs)

trunk/modules/TauTagging.cc

@@ -180 +180 @@
 void TauTagging::Process()
 {
-  Candidate *jet, *tau;
+  Candidate *jet, *tau, *daughter;
+  TLorentzVector tauMomentum;
   Double_t pt, eta, phi;
   TObjArray *tauArray;
   map< Int_t, DelphesFormula * >::iterator itEfficiencyMap;
   DelphesFormula *formula;
-  Int_t pdgCode, charge;
+  Int_t pdgCode, charge, i;
 
   // select taus
@@ -210 +211 @@
     while((tau = static_cast<Candidate *>(itTauArray.Next())))
     {
-      if(jetMomentum.DeltaR(tau->Momentum) <= fDeltaR)
+      tauMomentum.SetPxPyPzE(0.0, 0.0, 0.0, 0.0);
+      for(i = tau->D1; i <= tau->D2; ++i)
+      {
+        daughter = static_cast<Candidate *>(fParticleInputArray->At(i));
+        if(TMath::Abs(daughter->PID) == 16) continue;
+        tauMomentum += daughter->Momentum;
+      }
+      if(jetMomentum.DeltaR(tauMomentum) <= fDeltaR)
       {
         pdgCode = 15;