Changeset 157fe13 in git

Timestamp:

Jan 12, 2022, 8:20:47 PM (3 years ago)

Author:

christinaw97 <christina.wang@…>

Children:

4aec383

Parents:

d1ab205

Message:

LLPFilter hard-coded decay region requirements are moved to configuration cards

Files:

• cards/delphes_card_CMS_CSCCluster.tcl (modified) (2 diffs)
• modules/LLPFilter.cc (modified) (5 diffs)
• modules/LLPFilter.h (modified) (1 diff)

cards/delphes_card_CMS_CSCCluster.tcl

@@ -411 +411 @@
   set InputArray Delphes/allParticles
   set OutputArray LLP
-  # DecayRegion = 0: no cuts on decay region
-  # DecayRegion = 1: select LLP that decays in CSC volume
-  # DecayRegion = 2: select LLP that decays outside of calorimeters, for genMET calculation
-  set DecayRegion 1
+
+  # CMS CSC region
+  # used detector geometry in Figure 4.1.1, page141 from CERN-LHCC-97-032: https://cds.cern.ch/record/343814?ln=en
+
+  set RequireDecayRegion true
+  set DecayRegionRMax 6955
+  set DecayRegionRMin 0
+  set DecayRegionZMax 11000
+  set DecayRegionZMin 4000
+  set DecayRegionEtaMax 2
+  set DecayRegionEtaMin 0
+
   set RequireStatus false
   add PdgCode {1500001}
@@ -424 +432 @@
   set InputArray Delphes/allParticles
   set OutputArray LLP
-  # DecayRegion = 0: no cuts on decay region
-  # DecayRegion = 1: select LLP that decays in CSC volume
-  # DecayRegion = 2: select LLP that decays outside of calorimeters, for genMET calculation
-  set DecayRegion 0
+  set RequireDecayRegion false
+
+
   set RequireStatus false
   add PdgCode {1500001}

modules/LLPFilter.cc

@@ -22 +22 @@
  *  The classification of EM and hadronic energy of LLP is based on instructions from the HEPData entry for the CMS paper searching
  *  for neutral LLPs in the CMS endcap muon detectors: https://www.hepdata.net/record/104408
+ *  Muons and neutrinos are ignored. Photons, electrons, and pi0 are EM energy and everything else is hadronic energy.
  *
  *  \author Christina Wang
@@ -77 +78 @@
 
   fInvert = GetBool("Invert", false);
-  fDecayRegion = GetInt("DecayRegion", 0);
   fDaughterNumber = GetInt("DaughterNumber", 0);
-
+  fRequireDecayRegion = GetBool("RequireDecayRegion", 0);
+
+  fDecayRegionRMax = GetInt("DecayRegionRMax", 0); //mm
+  fDecayRegionRMin = GetInt("DecayRegionRMin", 0); //mm
+  fDecayRegionZMax = GetInt("DecayRegionZMax", 0); //mm
+  fDecayRegionZMin = GetInt("DecayRegionZMin", 0); //mm
+  fDecayRegionEtaMax = GetInt("DecayRegionEtaMax", 0); // requirement on abs(eta)
+  fDecayRegionEtaMin = GetInt("DecayRegionEtaMin", 0); //requirement on abs(eta)
 
 
@@ -148 +155 @@
     const TLorentzVector &candidateProdPosition = candidate->Position;
     const TLorentzVector &candidateDecayPosition = candidate->DecayPosition;
-    // TLorentzVector candidateDecayPosition;
     pt = candidateMomentum.Pt();
     eta = candidateMomentum.Eta();
@@ -184 +190 @@
       if (tempCandidate->M1 == -1) continue;
 
-      // candidateDecayPosition = daughter->Position;
-
+      // assign LLP EM or hadronic energy, depending on the daughter ID
       if (abs(daughterPdg)==11 || abs(daughterPdg)==22 || abs(daughterPdg)==111)candidate->Eem += daughterMomentum.E();
       else candidate->Ehad += daughterMomentum.E();
     }
 
-    // used detector geometry in Figure 4.1.1, page141 from CERN-LHCC-97-032: https://cds.cern.ch/record/343814?ln=en
-    // decayRegion = 0: no cuts on decay region
-    // decayRegion = 1: select LLP that decays in CSC volume
-    // decayRegion = 2: select LLP that decays outside of calorimeters
-    if (fDecayRegion == 1)
+    if (fRequireDecayRegion)
     {
-      if (abs(eta) < 2
-         && abs(candidateDecayPosition.Z())<11000 && abs(candidateDecayPosition.Z())>4000
-         && sqrt(pow(candidateDecayPosition.X(),2)+pow(candidateDecayPosition.Y(),2)) < 6955)
+      if (abs(eta) < fDecayRegionEtaMax && abs(eta) > fDecayRegionEtaMin
+         && abs(candidateDecayPosition.Z()) < fDecayRegionZMax && abs(candidateDecayPosition.Z()) > fDecayRegionZMin
+         && sqrt(pow(candidateDecayPosition.X(),2)+pow(candidateDecayPosition.Y(),2)) < fDecayRegionRMax
+         && sqrt(pow(candidateDecayPosition.X(),2)+pow(candidateDecayPosition.Y(),2)) > fDecayRegionRMin)
       {
         fOutputArray->Add(candidate);
@@ -204 +206 @@
 
     }
-    else if(fDecayRegion == 2)
-    {
-      if (abs(candidateDecayPosition.Z()) > 5680 && sqrt(pow(candidateDecayPosition.X(),2)+pow(candidateDecayPosition.Y(),2)) > 3000)
-      {
-        fOutputArray->Add(candidate);
-      }
-    }
     else{
       fOutputArray->Add(candidate);

modules/LLPFilter.h

@@ -51 +51 @@
 private:
   Double_t fPTMin; //!
-  Int_t fDecayRegion;
+  Bool_t fRequireDecayRegion;
+  Int_t fDecayRegionRMax;
+  Int_t fDecayRegionRMin;
+  Int_t fDecayRegionZMax;
+  Int_t fDecayRegionZMin;
+  Int_t fDecayRegionEtaMax;
+  Int_t fDecayRegionEtaMin;
   Int_t fDaughterNumber;
   Bool_t fInvert; //!