Changeset 7a86a5c in git

Timestamp:

Nov 23, 2018, 2:37:05 AM (6 years ago)

Author:

chekanov <chekanov@…>

Branches:

ImprovedOutputFile, Timing, llp, master

Children:

7c7fe5e

Parents:

e5fa629

Message:

Added new option to remove pileup particles from truth record

Location:

modules

Files:

• PdgCodeFilter.cc (modified) (2 diffs)
• PdgCodeFilter.h (modified) (1 diff)
• StatusPidFilter.cc (modified) (2 diffs)
• StatusPidFilter.h (modified) (1 diff)

modules/PdgCodeFilter.cc

@@ -76 +76 @@
   fInvert = GetBool("Invert", false);
 
+  // no pileup 
+  fRequireNotPileup = GetBool("RequireNotPileup", false);
+
   fRequireStatus = GetBool("RequireStatus", false);
   fStatus = GetInt("Status", 1);
@@ -127 +130 @@
     if(fRequireStatus && (candidate->Status != fStatus)) continue;
     if(fRequireCharge && (candidate->Charge != fCharge)) continue;
+    if(fRequireNotPileup && (candidate->IsPU   >0     )) continue;
 
     pass = kTRUE;

modules/PdgCodeFilter.h

@@ -55 +55 @@
   Bool_t fRequireCharge; //!
   Int_t fCharge; //!
-
+  Bool_t fRequireNotPileup; //!
 
   std::vector<Int_t> fPdgCodes;

modules/StatusPidFilter.cc

@@ -153 +153 @@
 {
   // PT threshold
-
   fPTMin = GetDouble("PTMin", 0.5);
+
+  // keep or remove pileup particles
+  fRequireNotPileup = GetBool("RequireNotPileup", false);
 
   // import input array
@@ -227 +229 @@
     if(!pass || (candidate->Momentum.Pt() < fPTMin && !(is_b_hadron || is_b_quark || is_tau_daughter || is_W_daughter)) ) continue;
 
+    // not pileup particles
+    if(fRequireNotPileup && (candidate->IsPU >0)) continue;
+
     fOutputArray->Add(candidate);
   }

modules/StatusPidFilter.h

@@ -51 +51 @@
   Double_t fPTMin; //!
 
+  Bool_t fRequireNotPileup; //!
+
   TIterator *fItInputArray; //!