Changeset 84a097e in git for modules/TrackSmearing.cc

Timestamp:

Dec 2, 2019, 6:55:22 PM (5 years ago)

Author:

Michele Selvaggi <michele.selvaggi@…>

Branches:

Timing

Children:

2ad823e, 6fc566b

Parents:

c614dd7

Message:

first iteration of adding timing in Delphes

File:

• modules/TrackSmearing.cc (modified) (3 diffs)

modules/TrackSmearing.cc

@@ -158 +158 @@
   TLorentzVector beamSpotPosition;
   Candidate *candidate, *mother;
-  Double_t pt, eta, d0, d0Error, trueD0, dz, dzError, trueDZ, p, pError, trueP, ctgTheta, ctgThetaError, trueCtgTheta, phi, phiError, truePhi;
+  Double_t pt, eta, d0, d0Error, trueD0, dz, dzError, trueDZ, p, pError, trueP, ctgTheta, ctgThetaError, trueCtgTheta, phi, phiError, truePhi, m;
   Double_t x, y, z, t, px, py, pz, theta;
   Double_t q, r;
@@ -328 +328 @@
 
     theta = TMath::ACos(ctgTheta / TMath::Sqrt(1.0 + ctgTheta * ctgTheta));
+    m = candidate->Momentum.M();
+    
     candidate->Momentum.SetPx(p * TMath::Cos(phi) * TMath::Sin(theta));
     candidate->Momentum.SetPy(p * TMath::Sin(phi) * TMath::Sin(theta));
     candidate->Momentum.SetPz(p * TMath::Cos(theta));
-    candidate->Momentum.SetE(candidate->Momentum.Pt() * TMath::CosH(eta));
+    //candidate->Momentum.SetE(candidate->Momentum.Pt() * TMath::CosH(eta));
+    candidate->Momentum.SetE( TMath::Sqrt(m*m + p*p) );
     candidate->PT = candidate->Momentum.Pt();
 
@@ -382 +385 @@
     candidate->Yd = yd * 1.0E3;
     candidate->Zd = zd * 1.0E3;
-
+    candidate->Td = -9999*c_light*1E-3;
+    
     if(fApplyToPileUp || !candidate->IsPU)
     {