Changeset 1ad8eca in git

Timestamp:

Aug 24, 2021, 11:53:15 AM (3 years ago)

Author:

michele <michele.selvaggi@…>

Branches:

master

Children:

61dccd3

Parents:

5c257a5

Message:

fixed momentum cut and generalised to neutrals

Location:

modules

Files:

• TimeOfFlight.cc (modified) (6 diffs)
• TimeOfFlight.h (modified) (1 diff)
• TimeSmearing.cc (modified) (4 diffs)
• TimeSmearing.h (modified) (1 diff)

modules/TimeOfFlight.cc

@@ -52 +52 @@
 
 TimeOfFlight::TimeOfFlight() :
-  fItTrackInputArray(0), fItVertexInputArray(0)
+  fItInputArray(0), fItVertexInputArray(0)
 {
 }
@@ -71 +71 @@
 
   // import track input array
-  fTrackInputArray = ImportArray(GetString("TrackInputArray", "MuonMomentumSmearing/muons"));
-  fItTrackInputArray = fTrackInputArray->MakeIterator();
+  fInputArray = ImportArray(GetString("InputArray", "MuonMomentumSmearing/muons"));
+  fItInputArray = fInputArray->MakeIterator();
 
   // import vertex input array
@@ -86 +86 @@
 void TimeOfFlight::Finish()
 {
-  if(fItTrackInputArray) delete fItTrackInputArray;
+  if(fItInputArray) delete fItInputArray;
   if(fItVertexInputArray) delete fItVertexInputArray;
 }
@@ -103 +103 @@
   ComputeVertexMomenta();
 
-  fItTrackInputArray->Reset();
-  while((candidate = static_cast<Candidate *>(fItTrackInputArray->Next())))
+  fItInputArray->Reset();
+  while((candidate = static_cast<Candidate *>(fItInputArray->Next())))
   {
 
@@ -116 +116 @@
     // time at vertex from MC truth
     t_truth = candidateInitialPosition.T() * 1.0E-3 / c_light;
-
-    if (candidate->Position.Vect().Mag() < 5.) continue;
 
     // various options on how to calculate the vertex time
@@ -210 +208 @@
     while((constituent = static_cast<Candidate *>(itGenParts.Next())))
     {
-      fItTrackInputArray->Reset();
-      while((track = static_cast<Candidate *>(fItTrackInputArray->Next())))
+      fItInputArray->Reset();
+      while((track = static_cast<Candidate *>(fItInputArray->Next())))
       {
         // get gen part that generated track

modules/TimeOfFlight.h

@@ -48 +48 @@
   Int_t fVertexTimeMode;
 
-  TIterator *fItTrackInputArray; //!
+  TIterator *fItInputArray; //!
   TIterator *fItVertexInputArray; //!
 
-  const TObjArray *fTrackInputArray; //!
+  const TObjArray *fInputArray; //!
   const TObjArray *fVertexInputArray; //!
 

modules/TimeSmearing.cc

@@ -52 +52 @@
 
 TimeSmearing::TimeSmearing() :
-  fItTrackInputArray(0), fResolutionFormula(0)
+  fItInputArray(0), fResolutionFormula(0)
 {
         fResolutionFormula = new DelphesFormula;
@@ -74 +74 @@
 
   // import track input array
-  fTrackInputArray = ImportArray(GetString("TrackInputArray", "MuonMomentumSmearing/muons"));
-  fItTrackInputArray = fTrackInputArray->MakeIterator();
+  fInputArray = ImportArray(GetString("InputArray", "MuonMomentumSmearing/muons"));
+  fItInputArray = fInputArray->MakeIterator();
 
 
@@ -86 +86 @@
 void TimeSmearing::Finish()
 {
-  if(fItTrackInputArray) delete fItTrackInputArray;
+  if(fItInputArray) delete fItInputArray;
 }
 
@@ -100 +100 @@
   const Double_t c_light = 2.99792458E8;
 
-  fItTrackInputArray->Reset();
-  while((candidate = static_cast<Candidate *>(fItTrackInputArray->Next())))
+  fItInputArray->Reset();
+  while((candidate = static_cast<Candidate *>(fItInputArray->Next())))
   {
 

modules/TimeSmearing.h

@@ -49 +49 @@
   Int_t fVertexTimeMode;
 
-  TIterator *fItTrackInputArray; //!
+  TIterator *fItInputArray; //!
 
-  const TObjArray *fTrackInputArray; //!
+  const TObjArray *fInputArray; //!
 
   TObjArray *fOutputArray; //!