Changeset 3051ea17 in git for modules/TrackCovariance.cc

Timestamp:

Jul 13, 2020, 5:19:21 PM (4 years ago)

Author:

Michele Selvaggi <michele.selvaggi@…>

Branches:

master

Children:

a0f5d71

Parents:

c18dca6

Message:

added track curvature to candidate class and switched back to regular covariance matrix

File:

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

modules/TrackCovariance.cc

@@ -96 +96 @@
   Candidate *candidate, *mother;
   Double_t mass, p, pt, q, ct;
-  Double_t dd0, ddz, dphi, dct, dp, dpt;
+  Double_t dd0, ddz, dphi, dct, dp, dpt, dC;
   
 
@@ -115 +115 @@
     candidate->InitialPosition.SetXYZT(track.GetObsX().X(),track.GetObsX().Y(),track.GetObsX().Z(),candidatePosition.T());
     
-    // save full ACTS friendly covariance matrix internally (6x6: D, z0, phi0, theta, q/p, time)
-    candidate->TrackCovarianceACTS = track.GetCovACTS();
+    // save full covariance 5x5 matrix internally (D0, phi, Curvature, dz, ctg(theta))
+    candidate->TrackCovariance = track.GetCov();
 
     pt = candidate->Momentum.Pt();
@@ -128 +128 @@
     
     candidate->D0 = track.GetObsPar()[0];
+    candidate->Phi = track.GetObsPar()[1];
+    candidate->C = track.GetObsPar()[2];
     candidate->DZ = track.GetObsPar()[3];
+    candidate->CtgTheta = track.GetObsPar()[4];
     candidate->P  = track.GetObsP().Mag();
-    candidate->CtgTheta = track.GetObsPar()[4];
-    candidate->Phi = track.GetObsPar()[1];
-
     candidate->PT = pt;
     candidate->Charge = q;
@@ -142 +142 @@
     dpt       = 2 * TMath::Sqrt( track.GetCov()(2, 2))*pt*pt / (0.2998*fBz);
     dp        = TMath::Sqrt((1.+ct*ct)*dpt*dpt + 4*pt*pt*ct*ct*dct*dct/(1.+ct*ct)/(1.+ct*ct));
+    dC        = TMath::Sqrt(track.GetCov()(2, 2));
 
     candidate->ErrorD0 = dd0;
     candidate->ErrorDZ = ddz;
     candidate->ErrorP = dp;
+    candidate->ErrorC = dC;
     candidate->ErrorCtgTheta = dct;
     candidate->ErrorPhi = dphi;