Changeset 157fe13 in git for modules/LLPFilter.cc

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

File:

• modules/LLPFilter.cc (modified) (5 diffs)

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);