Changeset c61b5ce in git

Timestamp:

Dec 10, 2021, 2:11:26 AM (3 years ago)

Author:

christinaw97 <christina.wang@…>

Children:

3c59f98

Parents:

a09b75f

Message:

cleaning

Files:

• Makefile (modified) (2 diffs)
• cards/delphes_card_CMS_CSCCluster.tcl (modified) (5 diffs)
• classes/DelphesClasses.h (modified) (4 diffs)
• modules/TreeWriter.cc (modified) (3 diffs)

Makefile

@@ -523 +523 @@
         classes/SortableObject.h
 tmp/classes/DelphesCscClusterFormula.$(ObjSuf): \
-                classes/DelphesCscClusterFormula.$(SrcSuf) \
-                classes/DelphesCscClusterFormula.h
+        classes/DelphesCscClusterFormula.$(SrcSuf) \
+        classes/DelphesCscClusterFormula.h
 tmp/classes/DelphesCylindricalFormula.$(ObjSuf): \
         classes/DelphesCylindricalFormula.$(SrcSuf) \
@@ -762 +762 @@
         modules/CscClusterId.h \
         classes/DelphesClasses.h \
-  classes/DelphesFactory.h \
+        classes/DelphesFactory.h \
         classes/DelphesCscClusterFormula.h \
         external/ExRootAnalysis/ExRootClassifier.h \

cards/delphes_card_CMS_CSCCluster.tcl

@@ -5 +5 @@
 set ExecutionPath {
   ParticlePropagator
-
-
 
   ChargedHadronTrackingEfficiency
@@ -61 +59 @@
   CutBasedIDEfficiency
   ClusterEfficiency
-
 
   TreeWriter
@@ -414 +411 @@
   set InputArray Delphes/allParticles
   set OutputArray LLP
-  set DecayRegion 1
   # 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
   set RequireStatus false
   add PdgCode {1500001}
@@ -427 +424 @@
   set InputArray Delphes/allParticles
   set OutputArray LLP
-  set DecayRegion 0
-  set RequireStatus false
-
   # 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 RequireStatus false
   add PdgCode {1500001}
 
@@ -1034 +1030 @@
 module TreeWriter TreeWriter {
 # add Branch InputArray BranchName BranchClass
-#add Branch Delphes/allParticles Particle GenParticle
-  #add Branch Delphes/stableParticles Particle GenParticle
-  #add Branch TrackMerger/tracks Track Track
-  #add Branch Calorimeter/towers Tower Tower
-
-  #add Branch HCal/eflowTracks EFlowTrack Track
-  #add Branch ECal/eflowPhotons EFlowPhoton Tower
-  #add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
-
-  #add Branch GenJetFinder/jets GenJet Jet
-  #add Branch GenMissingET/momentum GenMissingET MissingET
-
-
+# add Branch Delphes/allParticles Particle GenParticle
 
   add Branch UniqueObjectFinder/jets Jet Jet
   add Branch UniqueObjectFinder/electrons Electron Electron
-  #add Branch UniqueObjectFinder/photons Photon Photon
   add Branch UniqueObjectFinder/muons Muon Muon
 
-  #add Branch FatJetFinder/jets FatJet Jet
-
   add Branch MissingET/momentum MissingET MissingET
-  #add Branch ScalarHT/energy ScalarHT ScalarHT
   add Branch llpFilter/LLP llp CscCluster
   add Branch CSCFilter/LLP Cscllp CscCluster
-
   add Branch ClusterEfficiency/cluster CscCluster CscCluster
 }

classes/DelphesClasses.h

@@ -163 +163 @@
   Float_t decayT;
 
-  Float_t ctau;
-
   static CompBase *fgCompare; //!
   const CompBase *GetCompare() const { return fgCompare; }
@@ -650 +648 @@
 {
 public:
-  Int_t NHits; //nCSC hits
   Float_t Eta; // eta of LLP
   Float_t Phi; // phi of LLP
@@ -660 +657 @@
   Float_t Ehad; // had energy of LLP
   Float_t Eem; // em energy of LLP
-
   Float_t pid; // LLP pid
-
   Float_t T; // LLP decay time-photon travel time
   Float_t X; // LLP decay x
@@ -668 +663 @@
   Float_t Z; //  LLP decay z
   Float_t R; //  LLP decay z
-  Float_t beta;
-  Float_t ctau;
+  Float_t beta; // LLP beta
+  Float_t ctau; //LLP ctau
 
 

modules/TreeWriter.cc

@@ -245 +245 @@
     entry->decayZ = DecayPosition.Z();
     entry->decayT = DecayPosition.T()* 1.0E-3 / c_light;
-    float beta = entry->P/momentum.E();
-    float gamma = 1./sqrt(1-beta*beta);
-    entry->ctau = sqrt(pow(entry->decayX-entry->X,2)+pow(entry->decayY-entry->Y,2)+pow(entry->decayZ-entry->Z,2))/(beta*gamma);// in millimeter
   }
 }
@@ -930 +927 @@
     entry->Phi = momentum.Phi();
 
-
-    // entry->Eta = position.Eta();
-    // entry->Phi = position.Phi();
-
     entry->PT = momentum.Pt(); // pt of LLP
     entry->Px = momentum.Px();// px of LLP
@@ -940 +933 @@
     entry->E = momentum.E(); // E of LLP
     entry->pid = candidate->PID; // LLP pid
-    entry->Eem = candidate->Eem; // LLP pid
-    entry->Ehad = candidate->Ehad; // LLP pid
+    entry->Eem = candidate->Eem; // LLP Eem
+    entry->Ehad = candidate->Ehad; // LLP Ehad
     Double_t beta = momentum.P()/momentum.E();
     Double_t gamma = 1.0/sqrt(1-beta*beta);
     Double_t decayDistance = sqrt(pow(position.X(),2)+pow(position.Y(),2)+pow(position.Z(),2)); // mm
     entry->beta = beta; // LLP pid
-    entry->ctau = decayDistance/(beta * gamma);; // LLP pid
-
-    // entry->T = (position.T()-sqrt(pow(position.X(),2)+pow(position.Y(),2)+pow(position.Z(),2)))* 1.0E-3 / c_light; // LLP decay time-photon travel time
-
+    entry->ctau = decayDistance/(beta * gamma); // LLP travel time in its rest frame
     entry->T = decayDistance*(1./beta-1)* 1.0E-3/c_light*1e9; // ns
     entry->X = position.X(); // LLP decay x
     entry->Y = position.Y(); //  LLP decay y
     entry->Z = position.Z(); //  LLP decay z
-    // entry->Size = 100;
   }
 }