Changes in modules/TrackCovariance.cc [3051ea17:c18dca6] in git

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, dC;
+  Double_t dd0, ddz, dphi, dct, dp, dpt;
   
 
@@ -115 +115 @@
     candidate->InitialPosition.SetXYZT(track.GetObsX().X(),track.GetObsX().Y(),track.GetObsX().Z(),candidatePosition.T());
     
-    // save full covariance 5x5 matrix internally (D0, phi, Curvature, dz, ctg(theta))
-    candidate->TrackCovariance = track.GetCov();
+    // save full ACTS friendly covariance matrix internally (6x6: D, z0, phi0, theta, q/p, time)
+    candidate->TrackCovarianceACTS = track.GetCovACTS();
 
     pt = candidate->Momentum.Pt();
@@ -128 +128 @@
     
     candidate->D0 = track.GetObsPar()[0];
+    candidate->DZ = track.GetObsPar()[3];
+    candidate->P  = track.GetObsP().Mag();
+    candidate->CtgTheta = track.GetObsPar()[4];
     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->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;