1 #ifndef ALITPCSPACECHARGE_H
2 #define ALITPCSPACECHARGE_H
4 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
5 * See cxx source for full Copyright notice */
7 ////////////////////////////////////////////////////////////////////////////
9 // AliTPCSpaceCharge class //
10 // date: 19/06/2010 //
11 // Authors: Jim Thomas, Stefan Rossegger //
12 ////////////////////////////////////////////////////////////////////////////
14 #include "AliTPCCorrection.h"
17 class AliTPCSpaceCharge : public AliTPCCorrection {
20 virtual ~AliTPCSpaceCharge();
22 // initialization and update functions
24 virtual void Update(const TTimeStamp &timeStamp);
27 // common setters and getters for tangled ExB effect
28 virtual void SetOmegaTauT1T2(Float_t omegaTau,Float_t t1,Float_t 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);
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;}
37 // setters and getters for conical
38 void SetCorrectionFactor(Float_t correctionFactor) {fCorrectionFactor=correctionFactor;}
39 Float_t GetCorrectionFactor() const {return fCorrectionFactor;}
41 void InitSpaceChargeDistortion();
43 virtual void Print(const Option_t* option="") const;
46 virtual void GetCorrection(const Float_t x[],const Short_t roc,Float_t dx[]);
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
52 // look up table which was created for M_mb = 900 and IR = 3000
53 // compare Internal Note Nr: ???
55 Bool_t fInitLookUp; // flag to check it the Look Up table was created
57 Double_t fLookUpErOverEz[kNZ][kNR]; // Array to store electric field integral (int Er/Ez)
58 Double_t fLookUpDeltaEz[kNZ][kNR]; // Array to store electric field integral (int Delta Ez)
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.
63 enum {kIterations=100}; // Number of iterations within the poisson relaxation
65 ClassDef(AliTPCSpaceCharge,0);