Script supersceeded by AliForwarddNdetaTask.C and
[u/mrichter/AliRoot.git] / 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, const TVectorD *const fixedParams);
25   Bool_t Init();
26   void   SetActive(Bool_t flag){ fIsActive = flag;}
27   Bool_t IsActive() const {return fIsActive;}
28   //
29   Double_t GetParam() const {return fParam;}
30   void SetParam(Double_t param) {fParam=param;}
31   Double_t GetSigma() const {return fSigma;}
32   Double_t GetSigmaMax() const {return fSigmaMax;}
33   Double_t GetSigma2Time() const {return fSigma2Time;}
34   //
35   static TBits * BitsSide(Bool_t aside);
36   static TBits * BitsAll();
37   static void RegisterFormula(const char * name, GenFuncG formula);
38   static AliTPCTransformation::GenFuncG FindFormula(const char * name);
39   static Double_t Eval(const char * name, const Double_t*x,const Double_t*par);
40
41  private:
42   //
43   TString  * fNameX;         // x formula
44   TString  * fNameY;         // y formula
45   TString  * fNameZ;         // z formula  
46   //  
47   TBits    * fBitMask;       // bitmaps - transformation only for specified volID
48   Int_t      fCoordSystem;   // coord system of  output deltas 
49   Double_t   fParam;         // free parameter of transformation
50   Double_t   fSigma;         // error of the parameter
51   Double_t   fSigmaMax;      // maximal sigma (Not allowed to increase in propagate time by bigger factor)
52   Double_t   fSigma2Time;    // change of the error in time (per hour) - (For kalman filter) 
53   TVectorD  *fFixedParam;    // fixed parameters of tranformation
54   Bool_t     fIsActive;      // switch - is transformation active
55   //
56   // predefined formulas
57   //
58   static  Int_t          BuildBasicFormulas(); //build list of basic formulas
59   static  Double_t       TPCscalingRPol(Double_t *xyz, Double_t * param);
60   static  Double_t       TPCscalingZDrift(Double_t *xyz, Double_t * param);
61   static  Double_t       TPCscalingZDriftGy(Double_t *xyz, Double_t * param);
62   static  Double_t       TPCscalingZDriftT0(Double_t *xyz, Double_t * param);
63   static  Double_t       TPCscalingPhiLocal(Double_t *xyz, Double_t * param);
64   static  Double_t       TPClocalRPhiEdge(Double_t *xyz, const Double_t *const param);
65   //
66   // TPC Field cage + ROC misalingment induced distortion
67   //
68   static  Double_t       TPCscalingRIFC(Double_t *xyz, Double_t * param); // inner field cage r distorion
69   static  Double_t       TPCscalingROFC(Double_t *xyz, Double_t * param); // outer field cage r distorion
70   //
71   // TPC field cage + ROC misalignemnt induced distortion
72   //
73   static  Double_t       TPCdeltaFCROC(Double_t *xyz, const Double_t *const param); 
74   static  Double_t       TPCdeltaFCCE(Double_t *xyz, const Double_t *const param); 
75
76   //
77   // TPC local misalignment
78   //
79   static  Double_t       TPClocaldLxdGX(Double_t *xyz, const Double_t *const param);
80   static  Double_t       TPClocaldLxdGY(Double_t *xyz, const Double_t *const param);
81   static  Double_t       TPClocaldLydGX(Double_t *xyz, const Double_t *const param);
82   static  Double_t       TPClocaldLydGY(Double_t *xyz, const Double_t *const param);
83   static  Double_t       TPClocaldRzdGX(Double_t *xyz, const Double_t *const param);
84   static  Double_t       TPClocaldRzdGY(Double_t *xyz, const Double_t *const param);
85
86   //
87   // TPC  quadrant misalignment
88   //
89   //  static  Double_t       TPCQuadrantDr(Double_t *xyz, Double_t * param){return 0;}
90   //static  Double_t       TPCQuadrantDrphi(Double_t *xyz, Double_t * param){return 0;}
91   //
92   // Z shift -
93   //
94   static  Double_t       TPCDeltaZ(Double_t *xyz, const Double_t *const param);
95   static  Double_t       TPCDeltaZMediumLong(Double_t *xyz, Double_t * param);
96   static  Double_t       TPCTiltingZ(Double_t *xyz, const Double_t *const param);
97   //
98   Bool_t    fInit;          // initialization flag
99   GenFuncG  fFormulaX;      //! x formula
100   GenFuncG  fFormulaY;      //! y formula
101   GenFuncG  fFormulaZ;      //! z formula
102   static  GenFuncG    fgFormulas[10000];   //! array of pointers to formula
103   static  TObjArray*  fgFormulasName;      //! array of formalas name
104
105   AliTPCTransformation &operator=(const AliTPCTransformation&);    // not implemented
106
107   ClassDef(AliTPCTransformation,2);
108 };
109
110 #endif
111