]>
Commit | Line | Data |
---|---|---|
c9cbd2f2 | 1 | #ifndef ALITPCROCVOLTERROR3D_H |
2 | #define ALITPCROCVOLTERROR3D_H | |
3 | ||
4 | /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * | |
5 | * See cxx source for full Copyright notice */ | |
6 | ||
7 | //////////////////////////////////////////////////////////////////////////// | |
8 | // // | |
9 | // AliTPCROCVoltError3D class // | |
10 | // date: 01/06/2010 // | |
11 | // Authors: Jim Thomas, Stefan Rossegger // | |
12 | //////////////////////////////////////////////////////////////////////////// | |
13 | ||
14 | #include "AliTPCCorrection.h" | |
15 | #include "TH2F.h" | |
16 | ||
17 | ||
18 | class AliTPCROCVoltError3D : public AliTPCCorrection { | |
19 | public: | |
20 | AliTPCROCVoltError3D(); | |
21 | virtual ~AliTPCROCVoltError3D(); | |
22 | ||
23 | // initialization and update functions | |
24 | virtual void Init(); | |
25 | virtual void Update(const TTimeStamp &timeStamp); | |
26 | ||
27 | // common setters and getters for tangled ExB effect | |
28 | virtual void SetOmegaTauT1T2(Float_t omegaTau,Float_t t1,Float_t t2) { | |
29 | fT1=t1; fT2=t2; | |
30 | const Double_t wt0=t2*omegaTau; fC0=1./(1.+wt0*wt0); | |
31 | const Double_t wt1=t1*omegaTau; fC1=wt1/(1.+wt1*wt1); | |
32 | }; | |
33 | void SetC0C1(Float_t c0,Float_t c1) {fC0=c0;fC1=c1;} // CAUTION: USE WITH CARE | |
34 | Float_t GetC0() const {return fC0;} | |
35 | Float_t GetC1() const {return fC1;} | |
acf5907b | 36 | void SetROCData(TMatrixD * matrix); |
c9cbd2f2 | 37 | // setters and getters |
acf5907b | 38 | void SetROCDataFileName(const char * fname); |
39 | const Char_t* GetROCDataFileName() const {return fROCDataFileName.Data();} | |
c9cbd2f2 | 40 | |
41 | // flag to wheter or not include the z aligment in the dz calculation | |
42 | // if FALSE, the dz offset is purely due to the electric field change | |
7a348589 | 43 | void SetROCDisplacement(Bool_t flag) { |
44 | if (flag!=fROCdisplacement) { fROCdisplacement = flag; fInitLookUp=kFALSE; } | |
45 | } | |
c9cbd2f2 | 46 | Bool_t GetROCDisplacement() const { return fROCdisplacement; } |
7a348589 | 47 | |
48 | // flag on wheter to consider the difference in the electron arrival between IROC and OROC | |
49 | // due to the different position of the Anode wires | |
50 | void SetElectronArrivalCorrection(Bool_t flag) { | |
51 | if (flag!=fElectronArrivalCorrection) { fElectronArrivalCorrection = flag; fInitLookUp=kFALSE; } | |
52 | } | |
53 | Bool_t GetElectronArrivalCorrection() const { return fElectronArrivalCorrection; } | |
c9cbd2f2 | 54 | |
55 | ||
56 | void InitROCVoltError3D(); // Fill the lookup tables | |
57 | ||
58 | Float_t GetROCVoltOffset(Int_t side, Float_t r0, Float_t phi0); | |
7a348589 | 59 | TH2F* CreateHistoOfZAlignment(Int_t side, Int_t nx=250, Int_t ny=250); |
c9cbd2f2 | 60 | |
61 | virtual void Print(const Option_t* option="") const; | |
62 | ||
63 | protected: | |
64 | virtual void GetCorrection(const Float_t x[],const Short_t roc,Float_t dx[]); | |
65 | ||
66 | private: | |
67 | ||
68 | AliTPCROCVoltError3D(const AliTPCROCVoltError3D &); // not implemented | |
69 | AliTPCROCVoltError3D &operator=(const AliTPCROCVoltError3D &); // not implemented | |
70 | ||
71 | Float_t fC0; // coefficient C0 (compare Jim Thomas's notes for definitions) | |
72 | Float_t fC1; // coefficient C1 (compare Jim Thomas's notes for definitions) | |
73 | ||
74 | Bool_t fROCdisplacement; // flag on wheter to consider the ROC displacement | |
75 | // when calculating the z distortions | |
7a348589 | 76 | Bool_t fElectronArrivalCorrection; // flag on wheter to consider the difference |
77 | // in the electron arrival between IROC and OROC | |
78 | // due to the different position of the Anode wires | |
79 | ||
80 | ||
c9cbd2f2 | 81 | Bool_t fInitLookUp; // flag to check it the Look Up table was created (SUM) |
82 | ||
83 | TMatrixD *fLookUpErOverEz[kNPhi]; // Array to store electric field integral (int Er/Ez) | |
84 | TMatrixD *fLookUpEphiOverEz[kNPhi]; // Array to store electric field integral (int Er/Ez) | |
85 | TMatrixD *fLookUpDeltaEz[kNPhi]; // Array to store electric field integral (int Er/Ez) | |
86 | ||
acf5907b | 87 | TString fROCDataFileName; // filename of the survey data containing the lin Fit values |
c9cbd2f2 | 88 | TMatrixD *fdzDataLinFit; // Linear fits of dz survey points (each sector=72) (z0,slopeX,slopeY) |
89 | ||
90 | // basic numbers for the poisson relaxation //can be set individually in each class | |
91 | enum {kRows =257}; // grid size in r direction used in the poisson relaxation // ( 2**n + 1 ) eg. 65, 129, 257 etc. | |
92 | enum {kColumns=129}; // grid size in z direction used in the poisson relaxation // ( 2**m + 1 ) eg. 65, 129, 257 etc. | |
8847ede1 | 93 | enum {kPhiSlicesPerSector=10}; // phi slices per sector |
c9cbd2f2 | 94 | enum {kPhiSlices = 18*kPhiSlicesPerSector }; // number of points in phi for the basic lookup tables |
95 | enum {kIterations=100}; // Number of iterations within the poisson relaxation | |
96 | ||
8847ede1 | 97 | ClassDef(AliTPCROCVoltError3D,1); |
c9cbd2f2 | 98 | }; |
99 | ||
100 | #endif |