Changeset 9e08f27 in git

Timestamp:

Apr 17, 2021, 10:25:04 AM (4 years ago)

Author:

Michele Selvaggi <michele.selvaggi@…>

Branches:

master

Children:

3b3071a

Parents:

08d82fe

Message:

added e/mu mistag rate for taus (need to update partonArray in other than FWLite readers)

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;
-    else if(pdgCode == 24)
+    //if(pdgCode == 11 || pdgCode == 13 || pdgCode == 15)
+    //  return -1;
+    if(pdgCode == 24)
     {
       if(daughter1->D1 < 0) return -1;
@@ -96 +96 @@
     }
   }
-
   return 0;
 }
@@ -190 +189 @@
 void TauTagging::Process()
 {
-  Candidate *jet, *tau, *daughter;
+  Candidate *jet, *tau, *daughter, *part;
   TLorentzVector tauMomentum;
   Double_t pt, eta, phi, e, eff;
@@ -204 +203 @@
   // loop over all input jets
   fItJetInputArray->Reset();
+
   while((jet = static_cast<Candidate *>(fItJetInputArray->Next())))
   {
+
     const TLorentzVector &jetMomentum = jet->Momentum;
     pdgCode = 0;
@@ -243 +244 @@
       }
     }
+    
+    // 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);
-    }
-    else if(pdgCode <= 5 || pdgCode == 21 || pdgCode == 15)
+      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)
     {
       partonOutputArray->Add(candidate);