beamtime/cern-oct12/go4/RICH/CbmRichRingFitterEllipseTau.h (r4864/r3568)

Go to the documentation of this file.

#ifndef CBM_RICH_RING_FITTER_ELLIPSE_TAU
#define CBM_RICH_RING_FITTER_ELLIPSE_TAU

#include <vector>
#include "CbmRichRingFitterEllipseBase.h"
#include "TMath.h"
#include "TVectorD.h"
#include "TMatrixD.h"
#include "TMatrixDEigen.h"

#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