TPC/AliTPCTransformation.h

   1 #ifndef ALITPCTRANSFORMATION_H
   2 #define ALITPCTRANSFORMATION_H
   3
   4 #include "TNamed.h"
   5 #include "TMatrixD.h"
   6 #include "TVectorD.h"
   7 #include "TObjArray.h"
   8 class TTreeSRedirector;
   9 class AliTrackPointArray;
  10 class AliTrackPoint;
  11 class TFormula;
  12 class TBits;
  13 class TString;
  14
  15 class AliTPCTransformation: public TNamed{
  16 public:
  17   typedef Double_t (*GenFuncG)(const Double_t*,const Double_t*);
  18   AliTPCTransformation();
  19   AliTPCTransformation(const char *name,TBits *mask, const char *fx, const char *fy, const char  *fz, Int_t coord);
  20   AliTPCTransformation(const AliTPCTransformation&trafo);
  21   ~AliTPCTransformation();
  22   //
  23   virtual Double_t GetDeltaXYZ(Int_t coord, Int_t volID, Double_t param, Double_t x, Double_t y, Double_t z);
  24   void SetParams(Double_t param, Double_t sigma, Double_t sigma2Time, TVectorD* fixedParams);
  25   Bool_t Init();
  26   void   SetActive(Bool_t flag){ fIsActive = flag;}
  27   Bool_t IsActive(){return fIsActive;}
  28   //
  29   //
  30   static TBits * BitsSide(Bool_t aside);
  31   static TBits * BitsAll();
  32   static void RegisterFormula(const char * name, GenFuncG formula);
  33   static AliTPCTransformation::GenFuncG FindFormula(const char * name);
  34   static Double_t Eval(const char * name, const Double_t*x,const Double_t*par);
  35 public:
  36   //
  37   TString  * fNameX;         // x formula
  38   TString  * fNameY;         // y formula
  39   TString  * fNameZ;         // z formula
  40   //
  41   TBits    * fBitMask;       // bitmaps - transformation only for specified volID
  42   Int_t      fCoordSystem;   // coord system of  output deltas
  43   Double_t   fParam;         // free parameter of transformation
  44   Double_t   fSigma;         // error of the parameter
  45   Double_t   fSigma2Time;    // change of the error in time (per hour) - (For kalman filter)
  46   TVectorD  *fFixedParam;    // fixed parameters of tranformation
  47   Bool_t     fIsActive;      // switch - is transformation active
  48   //
  49   // predefined formulas
  50   //
  51   static  Int_t          BuildBasicFormulas(); //build list of basic formulas
  52   static  Double_t       TPCscalingRPol(Double_t *xyz, Double_t * param);
  53   static  Double_t       TPCscalingZDr(Double_t *xyz, Double_t * param);
  54   static  Double_t       TPCscalingPhiLocal(Double_t *xyz, Double_t * param);
  55   //
  56   // TPC Field cage + ROC misalingment induced distortion
  57   //
  58   static  Double_t       TPCscalingRIFC(Double_t *xyz, Double_t * param); // inner field cage r distorion
  59   static  Double_t       TPCscalingROFC(Double_t *xyz, Double_t * param); // outer field cage r distorion
  60   //
  61   // TPC local misalignment
  62   //
  63   static  Double_t       TPClocaldLxdGX(Double_t *xyz, Double_t * param);
  64   static  Double_t       TPClocaldLxdGY(Double_t *xyz, Double_t * param);
  65   static  Double_t       TPClocaldLydGX(Double_t *xyz, Double_t * param);
  66   static  Double_t       TPClocaldLydGY(Double_t *xyz, Double_t * param);
  67   //
  68   // TPC  quadrant misalignment
  69   //
  70   //  static  Double_t       TPCQuadrantDr(Double_t *xyz, Double_t * param){return 0;}
  71   //static  Double_t       TPCQuadrantDrphi(Double_t *xyz, Double_t * param){return 0;}
  72   //
  73   // Z shift -
  74   //
  75   static  Double_t       TPCDeltaZ(Double_t *xyz, Double_t * param);
  76   static  Double_t       TPCDeltaZMediumLong(Double_t *xyz, Double_t * param);
  77   static  Double_t       TPCTiltingZ(Double_t *xyz, Double_t * param);
  78   //
  79   Bool_t    fInit;          // initialization flag
  80   GenFuncG  fFormulaX;      //! x formula
  81   GenFuncG  fFormulaY;      //! y formula
  82   GenFuncG  fFormulaZ;      //! z formula
  83   static  GenFuncG    fgFormulas[10000];   //! array of pointers to formula
  84   static  TObjArray*  fgFormulasName;      //! array of formalas name
  85 private:
  86   AliTPCTransformation &operator=(const AliTPCTransformation&);    // not implemented
  87
  88   ClassDef(AliTPCTransformation,1);
  89 };
  90
  91 #endif
  92