Changeset 9e3f2fb in git for modules

Timestamp:

Jan 26, 2020, 4:52:03 PM (5 years ago)

Author:

michele <michele.selvaggi@…>

Branches:

Timing

Children:

f59b36de

Parents:

c7e90d8

Message:

adapted TimeSmearing to neutral particles

Location:

modules

Files:

• TimeSmearing.cc (modified) (5 diffs)
• TimeSmearing.h (modified) (2 diffs)

modules/TimeSmearing.cc

@@ -53 +53 @@
 
 TimeSmearing::TimeSmearing() :
-  fItInputArray(0), fFormula(0)
+  fItInputArray(0), fResolutionFormula(0)
 {
-        fFormula = new DelphesFormula;
+        fResolutionFormula = new DelphesFormula;
 }
 
@@ -61 +61 @@
 
 TimeSmearing::~TimeSmearing()
-{  
-        if(fFormula) delete fFormula;
+{
+        if(fResolutionFormula) delete fResolutionFormula;
 }
 
@@ -69 +69 @@
 void TimeSmearing::Init()
 {
-  // read resolution formula    
-  fFormula->Compile(GetString("TimeResolution", "0.001"));
+  // read time resolution formula in seconds
+  fResolutionFormula->Compile(GetString("TimeResolution", "30e-12"));
 
   // import input array
-  fEtaMax = GetDouble("EtaMax", 6.);
-  fInputArray = ImportArray(GetString("InputArray", "MuonMomentumSmearing/muons"));
+  fInputArray = ImportArray(GetString("InputArray", "TrackMerger/tracks"));
   fItInputArray = fInputArray->MakeIterator();
+
   // create output array
-
-  fOutputArray = ExportArray(GetString("OutputArray", "muons"));
+  fOutputArray = ExportArray(GetString("OutputArray", "tracks"));
 }
 
@@ -93 +92 @@
 {
   Candidate *candidate, *mother;
-  Double_t ti, tf_smeared, tf;
-  Double_t pt, eta, phi, e, p, l;
-  Double_t beta_particle;
+  Double_t ti,ti_smeared, tf, tf_smeared, dt;
+  Double_t eta, energy;
   Double_t timeResolution;
 
@@ -104 +102 @@
   while((candidate = static_cast<Candidate*>(fItInputArray->Next())))
   {
+    const TLorentzVector &candidatePosition = candidate->Position;
+    const TLorentzVector &candidateMomentum = candidate->Momentum;
+
+    // convert mm in seconds
     ti = candidate->InitialPosition.T()*1.0E-3/c_light;
     tf = candidate->Position.T()*1.0E-3/c_light;
-    
-    // dummy, only need to properly call TFormula
-    const TLorentzVector &candidatePosition = candidate->Position;
-    const TLorentzVector &candidateMomentum = candidate->Momentum;
+
     eta = candidatePosition.Eta();
-    phi = candidatePosition.Phi();
-    pt = candidateMomentum.Pt();
-    e = candidateMomentum.E();
-    p = candidateMomentum.P();
-    beta_particle = p/e;
-    l = candidate->L;
-    timeResolution = fFormula->Eval(e);
-    
+    energy = candidateMomentum.E();
+    timeResolution = fResolutionFormula->Eval(0.0, eta, 0.0, energy);
+
+    dt = timeResolution*gRandom->Gaus(0, 1);
+    tf_smeared = tf + dt;
+    ti_smeared = ti + dt;
+    mother = candidate;
+    candidate = static_cast<Candidate*>(candidate->Clone());
+    candidate->AddCandidate(mother);
+
+    candidate->Position.SetT(tf_smeared*1.0E3*c_light);
+    candidate->ErrorT = timeResolution*1.0E3*c_light;
+
+    // treating charged and neutral differently:
+    // for charged we smear the time after propagation, and put a dummy value for time at Vertex
+    // since the correct value will be computed after Vertexing4D
     if(candidate->Charge != 0)
-    { 
-      tf_smeared = tf + timeResolution*gRandom->Gaus(0, 1);
-      mother = candidate;
-      candidate = static_cast<Candidate*>(candidate->Clone());  // I am not sure that we need these lines !!!
-      candidate->AddCandidate(mother);
-      candidate->InitialPosition.SetT((100+ti)*1.0E3*c_light);  // Dummy Value, correct value will be computed by VertexFinderDA4D
-      candidate->Position.SetT(tf_smeared*1.0E3*c_light);
-      candidate->ErrorT = timeResolution*1.0E3*c_light;
-      fOutputArray->Add(candidate);
-    }
+      // Dummy Value, correct value will be computed by VertexFinderDA4D
+      candidate->InitialPosition.SetT((100+ti)*1.0E3*c_light);
     else
-    {
-      ti = timeResolution - l*1.0E3/(c_light*beta_particle);   
-      candidate->InitialPosition.SetT(ti);
-      candidate->ErrorT = timeResolution*1.0E3*c_light;         // Do we need to sum with 100 like in upside ?
-      fOutputArray->Add(candidate);
-    }
-   
+      candidate->InitialPosition.SetT(ti_smeared*1.0E3*c_light);
+
+    fOutputArray->Add(candidate);
+
   }
 }

modules/TimeSmearing.h

@@ -22 +22 @@
 /** \class TimeSmearing
  *
- *  Performs transverse time smearing.
+ *  Performs time of flight smearing.
  *
  *  \author Michele Selvaggi - UCL, Louvain-la-Neuve
@@ -45 +45 @@
 
 private:
-  DelphesFormula *fFormula;
-  Double_t fEtaMax;
 
+  DelphesFormula *fResolutionFormula;
+  
   TIterator *fItInputArray; //!