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