Changeset 010117a in git for modules/StatusPidFilter.cc

Timestamp:

Mar 9, 2018, 9:07:26 AM (7 years ago)

Author:

Martin <martin.flechl@…>

Branches:

ImprovedOutputFile, Timing, dual_readout, llp, master

Children:

0203656

Parents:

e4cc3f7

Message:

Saving all generator-level tau decay products to allow visible momentum determination.

File:

• modules/StatusPidFilter.cc (modified) (2 diffs)

modules/StatusPidFilter.cc

@@ -104 +104 @@
       return false;
   }
+
+  bool isTauDaughter(int pdgCode, int M1, const TObjArray *fInputArray){
+    //not needed, just to speed up the code - can be further refined but gives only negligible improvement:
+    if ( pdgCode==15 || pdgCode<11 || (pdgCode > 22 && pdgCode < 100) || pdgCode>1000 )
+      return false;
+
+    if ( M1 < 0 ) 
+      return false;
+
+    Candidate *mother;
+    mother = static_cast<Candidate*>(fInputArray->At( M1 ));
+    if ( TMath::Abs(mother->PID) == 15 )
+      return true;
+
+    return false;
+  }
+
 }
 
@@ -177 +194 @@
     bool is_b_quark  = (pdgCode == 5);
 
+    bool is_tau_daughter = isTauDaughter(pdgCode, candidate->M1, fInputArray);
+
     if (is_b_hadron)
       pass = kTRUE;
 
+    if (is_tau_daughter)
+      pass = kTRUE;
+
     // fPTMin not applied to b_hadrons / b_quarks to allow for b-enriched sample stitching
-    if(!pass || (candidate->Momentum.Pt() < fPTMin && !(is_b_hadron || is_b_quark)) ) continue;
+    // fPTMin not applied to tau decay products to allow visible-tau four momentum determination
+    if(!pass || (candidate->Momentum.Pt() < fPTMin && !(is_b_hadron || is_b_quark || is_tau_daughter)) ) continue;
 
     fOutputArray->Add(candidate);