Changeset 0b8551f in git for external/TrackCovariance/TrkUtil.h

Timestamp:

Nov 29, 2021, 4:04:38 PM (3 years ago)

Author:

GitHub <noreply@…>

Branches:

master

Children:

0c0c9af

Parents:

9a7ea36 (diff), bd376e3 (diff)
Note: this is a merge changeset, the changes displayed below correspond to the merge itself.
Use the (diff) links above to see all the changes relative to each parent.

git-author:

Michele Selvaggi <michele.selvaggi@…> (11/29/21 16:04:38)

git-committer:

GitHub <noreply@…> (11/29/21 16:04:38)

Message:

Merge pull request #102 from fbedesch/master

Major update to handle highly displaced tracks

File:

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

external/TrackCovariance/TrkUtil.h

@@ -7 +7 @@
 #include <TMatrixDSym.h>
 #include <TRandom.h>
+#include <TMath.h>
 //
 //
@@ -28 +29 @@
         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 +62 @@
                 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 +71 @@
         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 +89 @@
         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