Changeset 1ad8eca in git for modules/TimeOfFlight.cc

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

File:

• modules/TimeOfFlight.cc (modified) (6 diffs)

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