Changeset ebf40fd in git for external/TrackCovariance/TrkUtil.h

Timestamp:

Nov 29, 2021, 3:18:22 PM (3 years ago)

Author:

Franco BEDESCHI <bed@…>

Branches:

master

Children:

bd376e3

Parents:

9a7ea36

Message:

Major update to handle highly displaced tracks

File:

• external/TrackCovariance/TrkUtil.h (modified) (4 diffs)

external/TrackCovariance/TrkUtil.h

@@ -28 +28 @@
         TVectorD XPtoPar(TVector3 x, TVector3 p, Double_t Q);
         TVector3 ParToP(TVectorD Par);
+        TMatrixDSym RegInv(TMatrixDSym& Min);
+        //
+        // Track trajectory derivatives
+        TMatrixD derXdPar(TVectorD par, Double_t s);    // derivatives of position wrt parameters
+        TVectorD derXds(TVectorD par, Double_t s);              // derivatives of position wrt phase
+        TVectorD dsdPar_R(TVectorD par, Double_t R);    // derivatives of phase at constant R
+        TVectorD dsdPar_z(TVectorD par, Double_t z);    // derivatives of phase at constant z
         //
         // Conversion to ACTS parametrization
@@ -54 +61 @@
                 Double_t c = 2.99792458e8;      // speed of light m/sec
                 //return TMath::C()*1.0e-9;     // Incompatible with root5
-                return c*1.0e-9;                // Reduced speed of light       
+                return c*1.0e-9;                        // Reduced speed of light       
         }
         //
@@ -63 +70 @@
         static TVector3 ParToP(TVectorD Par, Double_t Bz);      // Get Momentum from track parameters
         static Double_t ParToQ(TVectorD Par);                           // Get track charge
+        static void LineDistance(TVector3 x0, TVector3 y0, TVector3 dirx, TVector3 diry, Double_t &sx, Double_t &sy, Double_t &distance);
+        //
+        // Track trajectory
+        //
+        static TVector3 Xtrack(TVectorD par, Double_t s);               // Parametric track trajectory
+        TVectorD derRphi_R(TVectorD par, Double_t R);   // Derivatives of R-phi at constant R
+        TVectorD derZ_R(TVectorD par, Double_t R);              // Derivatives of z at constant R
+        TVectorD derRphi_Z(TVectorD par, Double_t z);   // Derivatives of R-phi at constant z
+        TVectorD derR_Z(TVectorD par, Double_t z);              // Derivatives of R at constant z
+        //
+        // Smear with given covariance matrix
+        //
+        static TVectorD CovSmear(TVectorD x, TMatrixDSym C);
         //
         // Conversion from meters to mm
@@ -68 +88 @@
         static TVectorD ParToMm(TVectorD Par);                  // Parameter conversion
         static TMatrixDSym CovToMm(TMatrixDSym Cov);    // Covariance conversion
+        //
+        // Inside cylindrical volume
+        //
+        static Bool_t IsInside(TVector3 x, Double_t Rout, Double_t Zmin, Double_t Zmax)
+        {
+                Bool_t Is = kFALSE;
+                if (x.Pt() <= Rout && x.z() >= Zmin && x.z() <= Zmax)Is = kTRUE;
+                return Is;
+        }
         //
         // Cluster counting in gas