Analysis/VertexFit/include/VertexFit/FastVertexFit.h

Go to the documentation of this file.

#ifndef VertexFit_FastVertexFit_H
#define VertexFit_FastVertexFit_H
 
#include "CLHEP/Matrix/Matrix.h"
#include "CLHEP/Matrix/SymMatrix.h"
#include "CLHEP/Vector/ThreeVector.h"
using CLHEP::Hep3Vector;
using CLHEP::HepMatrix;
using CLHEP::HepSymMatrix;
using CLHEP::HepVector;
#include "CLHEP/Geometry/Point3D.h"
#ifndef ENABLE_BACKWARDS_COMPATIBILITY
typedef HepGeom::Point3D<double> HepPoint3D;
#endif
#include "VertexFit/VertexParameter.h"
#include <vector>
 
class FastVertexFit {
 
public:
  // constructor & deconstructor
  static FastVertexFit* instance();
  ~FastVertexFit();
 
  // initialization, must be called before fit
  void init();
 
  // void initVertex(); //FIXME need?
 
  // add track information
  void addTrack( const int number, const HepVector helix, const HepSymMatrix err );
 
  // fit function
  bool Fit();
 
  // Output
  double          chisq() const { return m_chisq; }
  VertexParameter vtx() const { return m_vtx_infit; }
  HepVector       Vx() const { return m_Vx; }
  HepSymMatrix    Evx() const { return m_Evx; }
  HepVector       Pull() const;
 
private:
  void updateMatrices( const HepVector helix, const HepSymMatrix err );
  // std::vector<HepVector> m_helix;
  // std::vector<HepSymMatrix> m_err;
  std::vector<HepMatrix>    m_D;
  std::vector<HepSymMatrix> m_W;
  std::vector<HepSymMatrix> m_DTWD;
  std::vector<HepVector>    m_xp;
 
private:
  double          m_chisq;     // chisq in fitting
  VertexParameter m_vtx_infit; // updated vertex parameters
  HepVector       m_Vx;        // updated vertex position
  HepSymMatrix    m_Evx;       // coviance matrix of updated vertex
 
private:
  FastVertexFit();                 // constructor
  static FastVertexFit* m_pointer; // pointer
};
 
#endif