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