Changeset 7dac4ea in git for external/TrackCovariance/VertexFit.h

Timestamp:

May 17, 2021, 6:05:45 PM (3 years ago)

Author:

GitHub <noreply@…>

Branches:

master

Children:

1c1c9c2, 33a6b3a, f8e61b2

Parents:

4afb18d (diff), 4acf2fd (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@…> (05/17/21 18:05:45)

git-committer:

GitHub <noreply@…> (05/17/21 18:05:45)

Message:

Merge pull request #95 from fbedesch/master

Fix for backward tracks and vertex fitting improvements

File:

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

external/TrackCovariance/VertexFit.h

@@ -23 +23 @@
         Int_t fNtr;                                     // Number of tracks
         std::vector<TVectorD*> fPar;            // Input parameter array
-        std::vector<TMatrixDSym*> fCov;// Input parameter covariances
+        std::vector<TVectorD*> fParNew;         // Updated parameter array
+        std::vector<TMatrixDSym*> fCov;         // Input parameter covariances
+        std::vector<TMatrixDSym*> fCovNew;      // Updated parameter covariances
         // Constraints
         Bool_t fVtxCst;                         // Vertex constraint flag
         TVectorD fxCst;                         // Constraint value
-        TMatrixDSym fCovCst;                    // Constraint covariance
+        TMatrixDSym fCovCst;                    // Constraint 
+        TMatrixDSym fCovCstInv;                 // Inverse of constraint covariance
         //
         // Results
@@ -39 +42 @@
         // Work arrays
         std::vector<Double_t> ffi;                              // Fit phases
-        std::vector<TVectorD*> fx0i;                            // Track expansion points
+        std::vector<TVectorD*> fx0i;                    // Track expansion points
         std::vector<TVectorD*> fai;                             // dx/dphi
+        std::vector<TVectorD*> fdi;                             // x-shift
         std::vector<Double_t> fa2i;                             // a'Wa
+        std::vector<TMatrixD*> fAti;                    // A transposed
         std::vector<TMatrixDSym*> fDi;                  // W-WBW
         std::vector<TMatrixDSym*> fWi;                  // (ACA')^-1