1 #ifndef ALITPCFCVOLTERROR3D_H
2 #define ALITPCFCVOLTERROR3D_H
4 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
5 * See cxx source for full Copyright notice */
7 ////////////////////////////////////////////////////////////////////////////
9 // AliTPCFCVoltError3D class //
10 // date: 01/06/2010 //
11 // Authors: Jim Thomas, Stefan Rossegger //
12 ////////////////////////////////////////////////////////////////////////////
14 #include "AliTPCCorrection.h"
17 class AliTPCFCVoltError3D : public AliTPCCorrection {
19 AliTPCFCVoltError3D();
20 virtual ~AliTPCFCVoltError3D();
22 // initialization and update functions
24 virtual void Update(const TTimeStamp &timeStamp);
26 // common setters and getters for tangled ExB effect
27 virtual void SetOmegaTauT1T2(Float_t omegaTau,Float_t t1,Float_t t2) {
29 const Double_t wt0=t2*omegaTau; fC0=1./(1.+wt0*wt0);
30 const Double_t wt1=t1*omegaTau; fC1=wt1/(1.+wt1*wt1);
32 void SetC0C1(Float_t c0,Float_t c1) {fC0=c0;fC1=c1;} // CAUTION: USE WITH CARE
33 Float_t GetC0() const {return fC0;}
34 Float_t GetC1() const {return fC1;}
36 // setters and getters
38 // Set rod shift in Voltage equivalents (40V ~ 1mm)
39 // rod numbers: 0-17 (IFC), 18-35 (OFC)
40 // note: strips move accordingly
41 void SetRodVoltShiftA(Int_t rod, Float_t voltOffset) {fRodVoltShiftA[rod]=voltOffset; fInitLookUp=kFALSE;}
42 void SetRodVoltShiftC(Int_t rod, Float_t voltOffset) {fRodVoltShiftC[rod]=voltOffset; fInitLookUp=kFALSE;}
43 Float_t GetRodVoltShiftA(Int_t i) const {return fRodVoltShiftA[i];}// 0-17: IFC, 18-35; OFC
44 Float_t GetRodVoltShiftC(Int_t i) const {return fRodVoltShiftC[i];}// 0-17: IFC, 18-35; OFC
46 // Set rotated clip (just at High Voltage RODs) in Voltage equivalents (40V ~ 1mm)
47 // rod number: 0 (IFC), 1 (OFC)
48 void SetRotatedClipVoltA(Int_t rod, Float_t voltOffset) {fRotatedClipVoltA[rod]=voltOffset; fInitLookUp=kFALSE;}
49 void SetRotatedClipVoltC(Int_t rod, Float_t voltOffset) {fRotatedClipVoltC[rod]=voltOffset; fInitLookUp=kFALSE;}
50 Float_t GetRotatedClipVoltA(Int_t i) const {return fRotatedClipVoltA[i];}// (0,1):(IFC,OFC)
51 Float_t GetRotatedClipVoltC(Int_t i) const {return fRotatedClipVoltC[i];}// (0,1):(IFC,OFC)
53 // Set rod shift in Voltage equivalents (40V ~ 1mm)
54 // rod numbers: 0-17 (OFC)
55 // note: strips DO NOT move, only the copper rods do ...
56 void SetCopperRodShiftA(Int_t rod, Float_t voltOffset) {fCopperRodShiftA[rod]=voltOffset; fInitLookUp=kFALSE;}
57 void SetCopperRodShiftC(Int_t rod, Float_t voltOffset) {fCopperRodShiftC[rod]=voltOffset; fInitLookUp=kFALSE;}
58 Float_t GetCopperRodShiftA(Int_t i) const {return fCopperRodShiftA[i];}// 0-17: IFC, 18-35; OFC
59 Float_t GetCopperRodShiftC(Int_t i) const {return fCopperRodShiftC[i];}// 0-17: IFC, 18-35; OFC
62 void InitFCVoltError3D(); // Fill the lookup tables
64 virtual void Print(const Option_t* option="") const;
67 virtual void GetCorrection(const Float_t x[],const Short_t roc,Float_t dx[]);
71 AliTPCFCVoltError3D(const AliTPCFCVoltError3D &); // not implemented
72 AliTPCFCVoltError3D &operator=(const AliTPCFCVoltError3D &); // not implemented
74 Float_t fC0; // coefficient C0 (compare Jim Thomas's notes for definitions)
75 Float_t fC1; // coefficient C1 (compare Jim Thomas's notes for definitions)
76 Float_t fRodVoltShiftA[36]; // Rod (&strips) shift in Volt (40V~1mm)
77 Float_t fRodVoltShiftC[36]; // Rod (&strips) shift in Volt (40V~1mm)
78 Float_t fRotatedClipVoltA[2]; // rotated clips at HV rod
79 Float_t fRotatedClipVoltC[2]; // rotated clips at HV rod
80 Float_t fCopperRodShiftA[36]; // only Rod shift
81 Float_t fCopperRodShiftC[36]; // only Rod shift
83 Bool_t fInitLookUp; // flag to check it the Look Up table was created (SUM)
84 Bool_t fInitLookUpBasic[6]; // ! flag if the basic lookup was created (shifted Rod (IFC,OFC) or rotated clip (IFC,OFC))
87 TMatrixD *fLookUpErOverEz[kNPhi]; // Array to store electric field integral (int Er/Ez)
88 TMatrixD *fLookUpEphiOverEz[kNPhi]; // Array to store electric field integral (int Er/Ez)
89 TMatrixD *fLookUpDeltaEz[kNPhi]; // Array to store electric field integral (int Er/Ez)
91 // basic numbers for the poisson relaxation //can be set individually in each class
92 enum {kRows =257}; // grid size in r direction used in the poisson relaxation // ( 2**n + 1 ) eg. 65, 129, 257 etc.
93 enum {kColumns=129}; // grid size in z direction used in the poisson relaxation // ( 2**m + 1 ) eg. 65, 129, 257 etc.
94 enum {kPhiSlicesPerSector = 10 }; // number of points in phi slices
95 enum {kPhiSlices = 1+kPhiSlicesPerSector*3 }; // number of points in phi for the basic lookup tables
96 enum {kIterations=100}; // Number of iterations within the poisson relaxation
98 // ugly way to store "partial" look up tables
99 // needed for the faster calculation of the final distortion table
101 // for Rod and Strip shift
102 TMatrixD *fLookUpBasic1ErOverEz[kPhiSlices]; // ! Array to store electric field integral (int Er/Ez)
103 TMatrixD *fLookUpBasic1EphiOverEz[kPhiSlices]; // ! Array to store electric field integral (int Ephi/Ez)
104 TMatrixD *fLookUpBasic1DeltaEz[kPhiSlices]; // ! Array to store electric field integral (int Ez)
106 TMatrixD *fLookUpBasic2ErOverEz[kPhiSlices]; // ! Array to store electric field integral
107 TMatrixD *fLookUpBasic2EphiOverEz[kPhiSlices]; // ! Array to store electric field integral
108 TMatrixD *fLookUpBasic2DeltaEz[kPhiSlices]; // ! Array to store electric field integral
111 TMatrixD *fLookUpBasic3ErOverEz[kPhiSlices]; // ! Array to store electric field integral
112 TMatrixD *fLookUpBasic3EphiOverEz[kPhiSlices]; // ! Array to store electric field integral
113 TMatrixD *fLookUpBasic3DeltaEz[kPhiSlices]; // ! Array to store electric field integral
115 TMatrixD *fLookUpBasic4ErOverEz[kPhiSlices]; // ! Array to store electric field integral
116 TMatrixD *fLookUpBasic4EphiOverEz[kPhiSlices]; // ! Array to store electric field integral
117 TMatrixD *fLookUpBasic4DeltaEz[kPhiSlices]; // ! Array to store electric field integral
119 // for (only rod) shift (copper rods)
120 TMatrixD *fLookUpBasic5ErOverEz[kPhiSlices]; // ! Array to store electric field integral
121 TMatrixD *fLookUpBasic5EphiOverEz[kPhiSlices]; // ! Array to store electric field integral
122 TMatrixD *fLookUpBasic5DeltaEz[kPhiSlices]; // ! Array to store electric field integral
124 TMatrixD *fLookUpBasic6ErOverEz[kPhiSlices]; // ! Array to store electric field integral
125 TMatrixD *fLookUpBasic6EphiOverEz[kPhiSlices]; // ! Array to store electric field integral
126 TMatrixD *fLookUpBasic6DeltaEz[kPhiSlices]; // ! Array to store electric field integral
129 ClassDef(AliTPCFCVoltError3D,2); //