CbmRichRingFitterEllipseTau.h

Go to the documentation of this file.

 
#ifndef CBM_RICH_RING_FITTER_ELLIPSE_TAU
#define CBM_RICH_RING_FITTER_ELLIPSE_TAU
 
#include "CbmRichRingFitterEllipseBase.h"
#include "TMath.h"
#include "TMatrixD.h"
#include "TMatrixDEigen.h"
#include "TVectorD.h"
#include <vector>
 
#include <iostream>
 
using std::vector;
 
 
class CbmRichRingFitterEllipseTau : public CbmRichRingFitterEllipseBase {
public:
  CbmRichRingFitterEllipseTau();
 
  virtual ~CbmRichRingFitterEllipseTau();
 
  virtual void DoFit(CbmRichRingLight* ring);
 
private:
  double fM[36];
  double fP[25];
  double fQ[25];
  double fZ[MAX_NOF_HITS_IN_RING * 6];
  double fZT[MAX_NOF_HITS_IN_RING * 6];
  double fAlgPar[6];
 
  void TransformEllipse(CbmRichRingLight* ring);
 
  void InitMatrices(CbmRichRingLight* ring);
 
  void Taubin();
 
  void Inv5x5();
 
  void AMultB(const double* const ap,
              int na,
              int ncolsa,
              const double* const bp,
              int nb,
              int ncolsb,
              double* cp);
 
  void Jacobi(double a[5][5], double d[5], double v[5][5]);
 
  void Eigsrt(double d[5], double v[5][5]);
};
 
#endif