1 #ifndef ALITPCTRANSFORMATION_H
2 #define ALITPCTRANSFORMATION_H
4 //-------------------------------------------------------
8 // Origin: marian.ivanov@cern.ch
9 // Code is not used anymore for the TPC corrections
10 // Obsolete - will be removed soon
12 //-------------------------------------------------------
18 #include "TObjArray.h"
19 class TTreeSRedirector;
20 class AliTrackPointArray;
26 class AliTPCTransformation: public TNamed{
28 typedef Double_t (*GenFuncG)(const Double_t*,const Double_t*);
29 AliTPCTransformation();
30 AliTPCTransformation(const char *name,TBits *mask, const char *fx, const char *fy, const char *fz, Int_t coord);
31 AliTPCTransformation(const AliTPCTransformation&trafo);
32 ~AliTPCTransformation();
34 virtual Double_t GetDeltaXYZ(Int_t coord, Int_t volID, Double_t param, Double_t x, Double_t y, Double_t z);
35 void SetParams(Double_t param, Double_t sigma, Double_t sigma2Time, const TVectorD *const fixedParams);
37 void SetActive(Bool_t flag){ fIsActive = flag;}
38 Bool_t IsActive() const {return fIsActive;}
40 Double_t GetParam() const {return fParam;}
41 void SetParam(Double_t param) {fParam=param;}
42 Double_t GetSigma() const {return fSigma;}
43 Double_t GetSigmaMax() const {return fSigmaMax;}
44 Double_t GetSigma2Time() const {return fSigma2Time;}
46 static TBits * BitsSide(Bool_t aside);
47 static TBits * BitsAll();
48 static void RegisterFormula(const char * name, GenFuncG formula);
49 static AliTPCTransformation::GenFuncG FindFormula(const char * name);
50 static Double_t Eval(const char * name, const Double_t*x,const Double_t*par);
54 TString * fNameX; // x formula
55 TString * fNameY; // y formula
56 TString * fNameZ; // z formula
58 TBits * fBitMask; // bitmaps - transformation only for specified volID
59 Int_t fCoordSystem; // coord system of output deltas
60 Double_t fParam; // free parameter of transformation
61 Double_t fSigma; // error of the parameter
62 Double_t fSigmaMax; // maximal sigma (Not allowed to increase in propagate time by bigger factor)
63 Double_t fSigma2Time; // change of the error in time (per hour) - (For kalman filter)
64 TVectorD *fFixedParam; // fixed parameters of tranformation
65 Bool_t fIsActive; // switch - is transformation active
67 // predefined formulas
69 static Int_t BuildBasicFormulas(); //build list of basic formulas
70 static Double_t TPCscalingRPol(Double_t *xyz, const Double_t * const param);
71 static Double_t TPCscalingZDrift(Double_t *xyz, const Double_t * const param);
72 static Double_t TPCscalingZDriftGy(Double_t *xyz, const Double_t * const param);
73 static Double_t TPCscalingZDriftT0(Double_t *xyz, const Double_t * const param);
74 static Double_t TPCscalingPhiLocal(Double_t *xyz, const Double_t * const param);
75 static Double_t TPClocalRPhiEdge(Double_t *xyz, const Double_t *const param);
77 // TPC Field cage + ROC misalingment induced distortion
79 static Double_t TPCscalingRIFC(Double_t *xyz, const Double_t * const param); // inner field cage r distorion
80 static Double_t TPCscalingROFC(Double_t *xyz, const Double_t * const param); // outer field cage r distorion
82 // TPC field cage + ROC misalignemnt induced distortion
84 static Double_t TPCdeltaFCROC(Double_t *xyz, const Double_t *const param);
85 static Double_t TPCdeltaFCCE(Double_t *xyz, const Double_t *const param);
88 // TPC local misalignment
90 static Double_t TPClocaldLxdGX(Double_t *xyz, const Double_t *const param);
91 static Double_t TPClocaldLxdGY(Double_t *xyz, const Double_t *const param);
92 static Double_t TPClocaldLydGX(Double_t *xyz, const Double_t *const param);
93 static Double_t TPClocaldLydGY(Double_t *xyz, const Double_t *const param);
94 static Double_t TPClocaldRzdGX(Double_t *xyz, const Double_t *const param);
95 static Double_t TPClocaldRzdGY(Double_t *xyz, const Double_t *const param);
98 // TPC quadrant misalignment
100 // static Double_t TPCQuadrantDr(Double_t *xyz, Double_t * param){return 0;}
101 //static Double_t TPCQuadrantDrphi(Double_t *xyz, Double_t * param){return 0;}
105 static Double_t TPCDeltaZ(Double_t *xyz, const Double_t *const param);
106 static Double_t TPCDeltaZMediumLong(Double_t *xyz, Double_t * param);
107 static Double_t TPCTiltingZ(Double_t *xyz, const Double_t *const param);
109 Bool_t fInit; // initialization flag
110 GenFuncG fFormulaX; //! x formula
111 GenFuncG fFormulaY; //! y formula
112 GenFuncG fFormulaZ; //! z formula
113 static GenFuncG fgFormulas[10000]; //! array of pointers to formula
114 static TObjArray* fgFormulasName; //! array of formalas name
116 AliTPCTransformation &operator=(const AliTPCTransformation&); // not implemented
118 ClassDef(AliTPCTransformation,2);