]>
Commit | Line | Data |
---|---|---|
c9cbd2f2 | 1 | #ifndef ALITPCSPACECHARGE_H |
2 | #define ALITPCSPACECHARGE_H | |
3 | ||
4 | /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * | |
5 | * See cxx source for full Copyright notice */ | |
6 | ||
7d855b04 | 7 | /// \class AliTPCSpaceCharge |
8 | /// \brief \author Jim Thomas, Stefan Rossegger | |
c9cbd2f2 | 9 | |
10 | #include "AliTPCCorrection.h" | |
11 | ||
12 | ||
7d855b04 | 13 | class TH2; |
276c23d3 | 14 | |
c9cbd2f2 | 15 | class AliTPCSpaceCharge : public AliTPCCorrection { |
16 | public: | |
17 | AliTPCSpaceCharge(); | |
18 | virtual ~AliTPCSpaceCharge(); | |
19 | ||
20 | // initialization and update functions | |
21 | virtual void Init(); | |
22 | virtual void Update(const TTimeStamp &timeStamp); | |
276c23d3 | 23 | void SetInputSpaceCharge(TH2 * hisSpaceCharge){fSpaceChargeHistogram = hisSpaceCharge;} // MI add |
7d855b04 | 24 | const TH2 * GetInputSpaceCharge(){return fSpaceChargeHistogram;} // MI add |
c9cbd2f2 | 25 | |
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;} | |
36 | ||
37 | // setters and getters for conical | |
38 | void SetCorrectionFactor(Float_t correctionFactor) {fCorrectionFactor=correctionFactor;} | |
39 | Float_t GetCorrectionFactor() const {return fCorrectionFactor;} | |
40 | ||
41 | void InitSpaceChargeDistortion(); | |
42 | ||
43 | virtual void Print(const Option_t* option="") const; | |
44 | ||
45 | protected: | |
46 | virtual void GetCorrection(const Float_t x[],const Short_t roc,Float_t dx[]); | |
47 | ||
48 | private: | |
7d855b04 | 49 | Float_t fC0; ///< coefficient C0 (compare Jim Thomas's notes for definitions) |
50 | Float_t fC1; ///< coefficient C1 (compare Jim Thomas's notes for definitions) | |
51 | Float_t fCorrectionFactor; ///< Space Charge Correction factor in comparison to initialized | |
c9cbd2f2 | 52 | // look up table which was created for M_mb = 900 and IR = 3000 |
53 | // compare Internal Note Nr: ??? | |
54 | ||
7d855b04 | 55 | TH2 * fSpaceChargeHistogram; ///< Histogram with the input space charge histogram - used as an optional input |
56 | Bool_t fInitLookUp; ///< flag to check it the Look Up table was created | |
c9cbd2f2 | 57 | |
7d855b04 | 58 | Double_t fLookUpErOverEz[kNZ][kNR]; ///< Array to store electric field integral (int Er/Ez) |
59 | Double_t fLookUpDeltaEz[kNZ][kNR]; ///< Array to store electric field integral (int Delta Ez) | |
c9cbd2f2 | 60 | // basic numbers for the poisson relaxation //can be set individually in each class |
61 | enum {kRows =257}; // grid size in r direction used in the poisson relaxation // ( 2**n + 1 ) eg. 65, 129, 257 etc. | |
62 | enum {kColumns=129}; // grid size in z direction used in the poisson relaxation // ( 2**m + 1 ) eg. 65, 129, 257 etc. | |
7d855b04 | 63 | enum {kIterations=100}; // Number of iterations within the poisson relaxation |
c9cbd2f2 | 64 | |
7d855b04 | 65 | /// \cond CLASSIMP |
66 | ClassDef(AliTPCSpaceCharge,0); | |
67 | /// \endcond | |
c9cbd2f2 | 68 | }; |
69 | ||
70 | #endif |