Changes in / [4564cba:7dbc149] in git

Files:

• cards/delphes_card_CMS.tcl (modified) (1 diff)
• modules/PdgCodeFilter.cc (modified) (3 diffs)
• modules/PdgCodeFilter.h (modified) (1 diff)

cards/delphes_card_CMS.tcl

@@ -585 +585 @@
 module PdgCodeFilter NeutrinoFilter {
 
-  #set InputArray Delphes/stableParticles
-  set InputArray Delphes/allParticles
+  set InputArray Delphes/stableParticles
   set OutputArray filteredParticles
 
   set PTMin 0.0
-  set RequireStatus true
-  set RequireKeepGhostBHadron true
 
   add PdgCode {12}

modules/PdgCodeFilter.cc

@@ -85 +85 @@
   fCharge = GetInt("Charge", 1);
 
-  // keep bhadron
-  fRequireKeepGhostBHadron = GetBool("RequireKeepGhostBHadron", false);
-
   // import input array
   fInputArray = ImportArray(GetString("InputArray", "Delphes/allParticles"));
@@ -127 +124 @@
   {
     pdgCode = candidate->PID;
-
-    if (fRequireKeepGhostBHadron) {
-      if (isBHadron(abs(pdgCode)) ){
-        candidate->PT = candidate->PT * 1e-18;
-        if (candidate->PT ==0) candidate->PT = 1e-18;
-        candidate->Momentum.SetPtEtaPhiM(candidate->PT, candidate->Momentum.Eta(), candidate->Momentum.Phi(), candidate->Momentum.M());
-        fOutputArray->Add(candidate);
-        continue;
-      }
-    } 
-
     const TLorentzVector &candidateMomentum = candidate->Momentum;
     pt = candidateMomentum.Pt();
@@ -153 +139 @@
   }
 }
-
-Bool_t PdgCodeFilter::isBHadron(const unsigned int absPdgId) {
-  if (absPdgId <= 100)
-    return false;  // Fundamental particles and MC internals
-  if (absPdgId >= 1000000000)
-    return false;  // Nuclei, +-10LZZZAAAI
-
-  // General form of PDG ID is 7 digit form
-  // +- n nr nL nq1 nq2 nq3 nJ
-  //const int nJ = absPdgId % 10; // Spin
-  const int nq3 = (absPdgId / 10) % 10;
-  const int nq2 = (absPdgId / 100) % 10;
-  const int nq1 = (absPdgId / 1000) % 10;
-
-  if (nq3 == 0)
-    return false;  // Diquarks
-  if (nq1 == 0 and nq2 == 5)
-    return true;  // B mesons
-  if (nq1 == 5)
-    return true;  // B baryons
-
-  return false;
-}
-

modules/PdgCodeFilter.h

@@ -55 +55 @@
   Bool_t fRequireNotPileup; //!
 
-  Bool_t fRequireKeepGhostBHadron; //!
-
   std::vector<Int_t> fPdgCodes;
-
-  Bool_t isBHadron(const unsigned int absPdgId); //!
 
   TIterator *fItInputArray; //!