Changes in external/TrackCovariance/VertexFit.h [82db145:127644a] in git
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
external/TrackCovariance/VertexFit.h
r82db145 r127644a 7 7 #include <TMatrixDSym.h> 8 8 #include "ObsTrk.h" 9 #include <vector>10 9 #include <iostream> 11 10 // … … 22 21 // Inputs 23 22 Int_t fNtr; // Number of tracks 24 std::vector<TVectorD*> fPar;// Input parameter array25 std::vector<TMatrixDSym*> fCov;// Input parameter covariances23 TVectorD** fPar; // Input parameter array 24 TMatrixDSym** fCov; // Input parameter covariances 26 25 // Constraints 27 26 Bool_t fVtxCst; // Vertex constraint flag 28 27 TVectorD fxCst; // Constraint value 29 TMatrixDSym fCovCst; 28 TMatrixDSym fCovCst; // Constraint covariance 30 29 // 31 30 // Results 32 Bool_t fVtxDone; // Flag vertex fit completed 33 Double_t fRold; // Current value of vertex radius 31 Bool_t fVtxDone; // Flag vertex fit completed 34 32 TVectorD fXv; // Found vertex 35 33 TMatrixDSym fcovXv; // Vertex covariance … … 37 35 TVectorD fChi2List; // List of Chi2 contributions 38 36 // 39 // Workarrays40 std::vector<Double_t>ffi; // Fit phases41 std::vector<TVectorD*> fx0i;// Track expansion points42 std::vector<TVectorD*>fai; // dx/dphi43 std::vector<Double_t>fa2i; // a'Wa44 std::vector<TMatrixDSym*>fDi; // W-WBW45 std::vector<TMatrixDSym*>fWi; // (ACA')^-146 std::vector<TMatrixDSym*>fWinvi; // ACA'37 // Transient arrays 38 Double_t* ffi; // Fit phases 39 TVectorD** fx0i; // Track expansion points 40 TVectorD** fai; // dx/dphi 41 Double_t* fa2i; // a'Wa 42 TMatrixDSym** fDi; // W-WBW 43 TMatrixDSym** fWi; // (ACA')^-1 44 TMatrixDSym** fWinvi; // ACA' 47 45 // 48 46 // Service routines 49 //void InitWrkArrays(); // Initializations 50 void ResetWrkArrays(); // Clear work arrays 51 Double_t StartRadius(); // Starting vertex radius determination 47 Double_t FastRv1(TVectorD p1, TVectorD p2); // Fast vertex radius determination 52 48 Double_t FastRv(TVectorD p1, TVectorD p2); // Fast vertex radius determination 53 TMatrixDSym RegInv (TMatrixDSym& Smat0);// Regularized 3D matrix inversion54 TMatrixD Fill_A(TVectorD par, Double_t phi); 55 TVectorD Fill_a(TVectorD par, Double_t phi); 49 TMatrixDSym RegInv3(TMatrixDSym& Smat0); // Regularized 3D matrix inversion 50 TMatrixD Fill_A(TVectorD par, Double_t phi); // Derivative of track position wrt track parameters 51 TVectorD Fill_a(TVectorD par, Double_t phi); // Derivative of track position wrt track phase 56 52 TVectorD Fill_x0(TVectorD par); // Track position at dma to z-axis 57 TVectorD Fill_x(TVectorD par, Double_t phi); // Track position at given phase 58 void UpdateTrkArrays(Int_t i); // Fill track realted arrays 59 void VertexFitter(); // Vertex finder routine 53 TVectorD Fill_x(TVectorD par, Double_t phi); // Track position at given phase 54 void VertexFinder(); // Vertex finder routine 60 55 public: 61 56 // 62 57 // Constructors 63 58 VertexFit(); // Initialize waiting for tracks 64 VertexFit(Int_t Ntr, ObsTrk** tracks); // Initialize with ObsTrk tracks59 VertexFit(Int_t Ntr, ObsTrk** tracks); // Initialize with ObsTrk tracks 65 60 VertexFit(Int_t Ntr, TVectorD** trkPar, TMatrixDSym** trkCov); // Initialize with parameters and covariances 66 61 // Destructor … … 76 71 // Handle tracks/constraints 77 72 void AddVtxConstraint(TVectorD xv, TMatrixDSym cov); // Add gaussian vertex constraint 78 void AddTrk(TVectorD *par, TMatrixDSym *Cov);// Add track to input list79 void RemoveTrk(Int_t iTrk); 73 void AddTrk(TVectorD par, TMatrixDSym Cov); // Add track to input list 74 void RemoveTrk(Int_t iTrk); // Remove iTrk track 80 75 // 81 76 };
Note:
See TracChangeset
for help on using the changeset viewer.