PndTrkLegendreFits.h

Go to the documentation of this file.

#ifndef PndTrkLegendreFits_H
#define PndTrkLegendreFits_H 1

// Root includes
#include "TROOT.h"

class PndTrkLegendreFits : public TObject
{

 private:

 int fIcounter;

 int    fNThetaDiv,
        fNRDiv;

 Double_t
        fDeltaT,
        fRMin,
        fThetaMax,
        fThetaMin;


 public:

  PndTrkLegendreFits();


  ~PndTrkLegendreFits(){};


 int FindMaximumInMatrix(
                int nRDiv,              // input, n. division in the R axis;
                int nThetaD,    // input  n. division in the Theta axis;
                UShort_t *Mat,  // input Matrix of dimensions [fNRDiv][fNThetaDiv] of
                                // which the Maximum has to be searched;

                int *iRMax,             // output, bin in R of the location of the Maximum;
                int *iTMax              // output, bin in Theta of the location of the Maximum;
                                                );


  Short_t FitHelixCylinder(
        Short_t nHitsinTrack,
        Double_t *Xconformal,
        Double_t *Yconformal,
        Double_t *DriftRadiusconformal,
        Double_t *ErrorDriftRadiusconformal,
        Double_t rotationangle,
        Double_t trajectory_vertex[2],
        Short_t NMAX,
        Double_t *m,
        Double_t *q,
        Double_t *pAlfa,
        Double_t *pBeta,
        Double_t *pGamma,
        bool *Type,
        int istampa,
        int IVOLTE
        );


  Short_t FitHelixCylinder2(
        Double_t * Cosine,
        Short_t LEGIANDRE_NTHETADIV,    // input, the theta divisions;
        Short_t LEGIANDRE_NRADIUSDIV,   // input, the theta divisions;
        Short_t nHitsinTrack,
        Double_t *Xconformal,
        Double_t *Yconformal,
        Double_t *DriftRadiusconformal,
        Double_t *ErrorDriftRadiusconformal,
        Double_t rotationangle,
        Double_t *Sinus,
        Double_t THETAMAX,
        Double_t THETAMIN,
        Double_t trajectory_vertex[2],
        Short_t NMAX,
        Double_t *m,
        Double_t *q,
        Double_t *pAlfa,
        Double_t *pBeta,
        Double_t *pGamma,
        bool *Type,
        int istampa,
        int IVOLTE
        );

  Short_t FitSZspace(
        Short_t nHitsinTrack,
        Double_t *S,
        Double_t *Z,
        Double_t *DriftRadius,
        Double_t *ErrorDriftRadius,
        Double_t FInot,
        Short_t NMAX,
        Double_t *emme,
        int PlotNumber
        );


  int LoadMatrix_FindMaximum(
        Short_t nHitsinTrack,                   // input
        Double_t *X,                    // X position (in conformal or SZ or whatever);
        Double_t *Y,                    // Y position (in conformal or SZ or whatever);
        Double_t *DriftRadius,          // negative if Mvd hit or similar;
        Double_t *ErrorDriftRadiusconformal,    // for the Mvd this is the Radius of the circumference
                                                // translated with the Conformal transformation,
                                                // which is, in the XY space : a) centered on the Pixel
                                                // or Strip; with radius = 0.01 cm --> therefore encompassing
                                                // completely the Pixel or Stip hit.

        Double_t *Rout, // output parameter of the straight line; Xcos(Theta)+Y*sin(Theta)=R;
        Double_t *Thetaout  // output parameter of the straight line; Xcos(Theta)+Y*sin(Theta)=R;
                                        );



  void LoadMatrix_FindMaximum2(
        Double_t * Cosine,              // input, the precalculated values of cosine;
        Short_t LEGIANDRE_NTHETADIV,    // input, the theta divisions;
        Short_t LEGIANDRE_NRADIUSDIV,   // input, the theta divisions;
        Short_t nHitsinTrack,           // input
        Double_t * Sinus,               // input, the precalculated values of sinus;
        Double_t THETAMAX,              // input, maximum of Theta range (usually 2PI radians);
        Double_t THETAMIN,              // input, minimum of Theta range (usually 0 radians);
        Double_t *X,                    // X position (in conformal or SZ or whatever);
        Double_t *Y,                    // Y position (in conformal or SZ or whatever);
        Double_t *DriftRadius,          // negative if Mvd hit or similar;
        Double_t *ErrorDriftRadiusconformal,    // for the Mvd this is the Radius of the circumference
                                                // translated with the Conformal transformation,
                                                // which is, in the XY space : a) centered on the Pixel
                                                // or Strip; with radius = 0.01 cm --> therefore encompassing
                                                // completely the Pixel or Stip hit.

        Double_t *Rout, // output parameter of the straight line; Xcos(Theta)+Y*sin(Theta)=R;
        Double_t *Thetaout  // output parameter of the straight line; Xcos(Theta)+Y*sin(Theta)=R;
                                        );

  ClassDef(PndTrkLegendreFits,1);

};

#endif