Changeset b750b0a in git for external/TrackCovariance/VertexFit.h

Timestamp:

Jan 19, 2022, 10:11:39 AM (3 years ago)

Author:

Franco BEDESCHI <bed@…>

Branches:

master

Children:

78ce8d1, 7bca620

Parents:

bd376e3

Message:

Fix CovarianceMatrix scaling error. Vertexing improvements.

File:

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

external/TrackCovariance/VertexFit.h

@@ -6 +6 @@
 #include <TVectorD.h>
 #include <TMatrixDSym.h>
+#include "TrkUtil.h"
 #include "ObsTrk.h"
 #include <vector>
@@ -12 +13 @@
 // Class for vertex fitting
 
-class VertexFit {
+class VertexFit: public TrkUtil
+{
         //
         // Vertex fitting with track parameters steering
@@ -21 +23 @@
         //
         // Inputs
-        Int_t fNtr;                                     // Number of tracks
+        Int_t fNtr;                                                     // Number of tracks
         std::vector<TVectorD*> fPar;            // Input parameter array
         std::vector<TVectorD*> fParNew;         // Updated parameter array
@@ -27 +29 @@
         std::vector<TMatrixDSym*> fCovNew;      // Updated parameter covariances
         // Constraints
-        Bool_t fVtxCst;                         // Vertex constraint flag
-        TVectorD fxCst;                         // Constraint value
+        Bool_t fVtxCst;                                 // Vertex constraint flag
+        TVectorD fxCst;                                 // Constraint value
         TMatrixDSym fCovCst;                    // Constraint 
         TMatrixDSym fCovCstInv;                 // Inverse of constraint covariance
         //
         // Results
-        Bool_t fVtxDone;                                // Flag vertex fit completed
+        Bool_t fVtxDone;                        // Flag vertex fit completed
         Double_t fRold;                         // Current value of vertex radius
         TVectorD fXv;                           // Found vertex
@@ -52 +54 @@
         //
         // Service routines
-        //void InitWrkArrays();                                                 // Initializations
         void ResetWrkArrays();                                                  // Clear work arrays
-        Double_t StartRadius();                                                 // Starting vertex radius determination
-        Double_t FastRv(TVectorD p1, TVectorD p2);              // Fast vertex radius determination
-        TMatrixDSym RegInv(TMatrixDSym& Smat0);                 // Regularized 3D matrix inversion 
-        TMatrixD Fill_A(TVectorD par, Double_t phi);            // Derivative of track position wrt track parameters
-        TVectorD Fill_a(TVectorD par, Double_t phi);            // Derivative of track position wrt track phase
+        //Double_t StartRadius();                                                       // Starting vertex radius determination
+        //Double_t FastRv(TVectorD p1, TVectorD p2);            // Fast vertex radius determination
+        //TMatrixDSym RegInv(TMatrixDSym& Smat0);               // Regularized 3D matrix inversion 
+        //TMatrixD Fill_A(TVectorD par, Double_t phi);  // Derivative of track position wrt track parameters
+        //TVectorD Fill_a(TVectorD par, Double_t phi);  // Derivative of track position wrt track phase
         TVectorD Fill_x0(TVectorD par);                                 // Track position at dma to z-axis
-        TVectorD Fill_x(TVectorD par, Double_t phi);            // Track position at given phase
+        TVectorD Fill_x(TVectorD par, Double_t phi);    // Track position at given phase
         void UpdateTrkArrays(Int_t i);                                  // Fill track realted arrays
-        void VertexFitter();                                                            // Vertex finder routine
+        void VtxFitNoSteer();                                                   // Vertex fitter routine w/o parameter steering
+        void VertexFitter();                                                    // Vertex fitter routine w/  parameter steering
 public:
         //