Changes in / [3b3071a:64294db] in git

Files:

• modules/TauTagging.cc (modified) (5 diffs)
• readers/DelphesCMSFWLite.cpp (modified) (1 diff)

modules/TauTagging.cc

@@ -83 +83 @@
     daughter1 = static_cast<Candidate *>(fParticleInputArray->At(i));
     pdgCode = TMath::Abs(daughter1->PID);
-    //if(pdgCode == 11 || pdgCode == 13 || pdgCode == 15)
-    //  return -1;
-    if(pdgCode == 24)
+    if(pdgCode == 11 || pdgCode == 13 || pdgCode == 15)
+      return -1;
+    else if(pdgCode == 24)
     {
       if(daughter1->D1 < 0) return -1;
@@ -96 +96 @@
     }
   }
+
   return 0;
 }
@@ -189 +190 @@
 void TauTagging::Process()
 {
-  Candidate *jet, *tau, *daughter, *part;
+  Candidate *jet, *tau, *daughter;
   TLorentzVector tauMomentum;
   Double_t pt, eta, phi, e, eff;
@@ -203 +204 @@
   // loop over all input jets
   fItJetInputArray->Reset();
-
   while((jet = static_cast<Candidate *>(fItJetInputArray->Next())))
   {
-
     const TLorentzVector &jetMomentum = jet->Momentum;
     pdgCode = 0;
@@ -244 +243 @@
       }
     }
-    
-    // fake electrons and muons
-    
-    if (pdgCode == 0)
-    {
-     
-      Double_t drMin = fDeltaR;   
-      fItPartonInputArray->Reset();
-      while((part = static_cast<Candidate *>(fItPartonInputArray->Next())))
-      {
-        if(TMath::Abs(part->PID) == 11 || TMath::Abs(part->PID) == 13) 
-        {
-            tauMomentum = part->Momentum;
-            if (tauMomentum.Pt() < fClassifier->fPTMin) continue;
-            if (TMath::Abs(tauMomentum.Eta()) > fClassifier->fEtaMax) continue;
-
-            Double_t dr = jetMomentum.DeltaR(tauMomentum);
-            if( dr < drMin)
-            {
-               drMin = dr;
-               pdgCode = TMath::Abs(part->PID);
-               charge = part->Charge;
-            }  
-        }
-      }
-    }
-
     // find an efficency formula
     itEfficiencyMap = fEfficiencyMap.find(pdgCode);

readers/DelphesCMSFWLite.cpp

@@ -217 +217 @@
       // Prevent duplicated particle.
       if(!isMiniAOD) stableParticleOutputArray->Add(candidate);
-      if (pdgCode == 11 || pdgCode == 13) partonOutputArray->Add(candidate);
-    }
-    //else if(pdgCode <= 5 || pdgCode == 21 || pdgCode == 15)
-    else if(pdgCode <= 5 || pdgCode == 21 || pdgCode == 11 || pdgCode == 13 || pdgCode == 15)
+    }
+    else if(pdgCode <= 5 || pdgCode == 21 || pdgCode == 15)
     {
       partonOutputArray->Add(candidate);