Changeset 2a3eb22 in git for modules/TreeWriter.cc

Timestamp:

Jun 24, 2016, 3:33:25 PM (8 years ago)

Author:

Michele Selvaggi <michele.selvaggi@…>

Branches:

ImprovedOutputFile, Timing, dual_readout, llp, master

Children:

734b267, 90ccc5d

Parents:

5c402f7

Message:

fix track momentum in TreeWriter

File:

• modules/TreeWriter.cc (modified) (4 diffs)

modules/TreeWriter.cc

@@ -317 +317 @@
   Candidate *particle = 0;
   Track *entry = 0;
-  Double_t pt, signz, cosTheta, eta, rapidity;
+  Double_t pt, signz, cosTheta, eta, rapidity, p, ctgTheta, phi;
   const Double_t c_light = 2.99792458E8;
 
@@ -354 +354 @@
     entry->DZ            = candidate->DZ;
     entry->ErrorDZ       = candidate->ErrorDZ;
-    entry->P             = candidate->P;
+
     entry->ErrorP        = candidate->ErrorP;
-    entry->PT            = candidate->PT;
     entry->ErrorPT       = candidate->ErrorPT;
-    entry->CtgTheta      = candidate->CtgTheta;
     entry->ErrorCtgTheta = candidate->ErrorCtgTheta;
-    entry->Phi           = candidate->Phi;
     entry->ErrorPhi      = candidate->ErrorPhi;
 
@@ -370 +367 @@
 
     pt = momentum.Pt();
+    p = momentum.P();
+    phi = momentum.Phi();
+    ctgTheta = (TMath::Tan(momentum.Theta()) != 0) ? 1/TMath::Tan(momentum.Theta()) : 1e10;
+
     cosTheta = TMath::Abs(momentum.CosTheta());
     signz = (momentum.Pz() >= 0.0) ? 1.0 : -1.0;
@@ -375 +376 @@
     rapidity = (cosTheta == 1.0 ? signz*999.9 : momentum.Rapidity());
 
+    entry->PT  = pt;
     entry->Eta = eta;
+    entry->Phi = phi;
+    entry->CtgTheta = ctgTheta;
 
     particle = static_cast<Candidate*>(candidate->GetCandidates()->At(0));