CbmRichRingFitterEllipseTau.h

Go to the documentation of this file.

/* Copyright (C) 2009-2012 UGiessen/JINR-LIT, Giessen/Dubna
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Alexander Ayriyan, Semen Lebedev [committer] */

 
#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 <iostream>
#include <vector>
 
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