Changeset 0203656 in git for modules/StatusPidFilter.cc

Timestamp:

Mar 12, 2018, 5:13:00 PM (7 years ago)

Author:

GitHub <noreply@…>

Branches:

ImprovedOutputFile, Timing, dual_readout, llp, master

Children:

070a19b, 7429c6a

Parents:

e4cc3f7 (diff), 010117a (diff)
Note: this is a merge changeset, the changes displayed below correspond to the merge itself.
Use the (diff) links above to see all the changes relative to each parent.

git-author:

Michele Selvaggi <michele.selvaggi@…> (03/12/18 17:13:00)

git-committer:

GitHub <noreply@…> (03/12/18 17:13:00)

Message:

Merge pull request #50 from mflechl/master

Saving all gen-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);