Changeset 3c59f98 in git

Timestamp:

Jan 11, 2022, 9:17:09 PM (3 years ago)

Author:

christinaw97 <christina.wang@…>

Children:

c6bce62

Parents:

c61b5ce

Message:

add example/ExampleCscCluster.py and remove decayPosition from GenParticle class

Files:

• classes/DelphesClasses.h (modified) (2 diffs)
• classes/DelphesHepMC2Reader.cc (modified) (2 diffs)
• examples/ExampleCscCluster.py (added)
• modules/LLPFilter.cc (modified) (2 diffs)
• modules/TreeWriter.cc (modified) (2 diffs)

classes/DelphesClasses.h

@@ -670 +670 @@
   const CompBase *GetCompare() const { return fgCompare; }
 
-  ClassDef(CscCluster, 4)
+  ClassDef(CscCluster, 5)
 };
 
@@ -716 +716 @@
   Float_t DeltaPhi;
 
-  TLorentzVector Momentum, Position, InitialPosition, DecayPosition, PositionError, Area;
+  TLorentzVector Momentum, Position, InitialPosition, PositionError, DecayPosition, Area;
 
   Float_t L; // path length

classes/DelphesHepMC2Reader.cc

@@ -481 +481 @@
         const TLorentzVector &decayPosition = candidate->Position;
         candidate->DecayPosition.SetXYZT(decayPosition.X(), decayPosition.Y(), decayPosition.Z(), decayPosition.T());// decay position
-      }
+     }
       else
       {
@@ -489 +489 @@
         const TLorentzVector &decayPosition = candidateDaughter->Position;
         candidate->DecayPosition.SetXYZT(decayPosition.X(), decayPosition.Y(), decayPosition.Z(), decayPosition.T());// decay position
-
-      }
-    }
-  }
-}
-
-//---------------------------------------------------------------------------
+        
+
+      }
+    }
+  }
+}
+
+//---------------------------------------------------------------------------

modules/LLPFilter.cc

@@ -170 +170 @@
       if (abs(daughterPdg)==12 || abs(daughterPdg)==14 || abs(daughterPdg)==16 || abs(daughterPdg)==13)continue; // ignore neutrinos and muons
       if (abs(daughterPdg) > 1000000) continue;//ignore BSM particles
-      
-      const TLorentzVector &daughterProdPosition = daughter->Position;
+
       const TLorentzVector &daughterMomentum = daughter->Momentum;
-      const TLorentzVector distance = daughterProdPosition - candidateDecayPosition;
 
       // look for mother until find LLP or reach the top of the tree
@@ -187 +185 @@
     }
 
-
-
     // 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

modules/TreeWriter.cc

@@ -198 +198 @@
     const TLorentzVector &momentum = candidate->Momentum;
     const TLorentzVector &position = candidate->Position;
-    const TLorentzVector &DecayPosition = candidate->DecayPosition;
 
     entry = static_cast<GenParticle *>(branch->NewEntry());
@@ -241 +240 @@
     entry->T = position.T() * 1.0E-3 / c_light;
 
-    entry->decayX = DecayPosition.X();
-    entry->decayY = DecayPosition.Y();
-    entry->decayZ = DecayPosition.Z();
-    entry->decayT = DecayPosition.T()* 1.0E-3 / c_light;
   }
 }