Changeset 41da326 in git for modules/StatusPidFilter.cc

Timestamp:

May 20, 2013, 10:10:25 PM (11 years ago)

Author:

pavel <pavel@…>

Branches:

ImprovedOutputFile, Timing, dual_readout, llp, master

Children:

d5104a4

Parents:

1ba5c52

Message:

add more particles to StatusPidFilter

File:

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

modules/StatusPidFilter.cc

@@ -80 +80 @@
   Candidate *candidate;
   Int_t status, pdgCode;
+  Bool_t pass;
 
   fItInputArray->Reset();
@@ -87 +88 @@
     pdgCode = TMath::Abs(candidate->PID);
 
-    // Write all electrons, muons, taus and status == 3;
-    if(pdgCode != 11 && pdgCode != 13 && pdgCode != 15 && status != 3) continue;
+    pass = kFALSE;
 
-    if(candidate->Momentum.Pt() <= fPTMin) continue;
+    // status == 3
+    if(status == 3) pass = kTRUE;
+
+    // electrons, muons, taus
+    if(pdgCode == 11 || pdgCode == 13 || pdgCode == 15) pass = kTRUE;
+
+    // neutrinos
+    if(pdgCode == 12 || pdgCode == 14 || pdgCode == 16) pass = kTRUE;
+
+    // heavy quarks
+    if(pdgCode == 5 || pdgCode == 6) pass = kTRUE;
+
+    // Gauge bosons and other fundamental bosons
+    if(pdgCode > 22 &&  pdgCode < 43) pass = kTRUE;
+
+    if(!pass || candidate->Momentum.Pt() <= fPTMin) continue;
 
     fOutputArray->Add(candidate);