TPC/AliTPCSpaceCharge.h

   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
   7 ////////////////////////////////////////////////////////////////////////////
   8 //                                                                        //
   9 // AliTPCSpaceCharge class                                                //
  10 // date: 19/06/2010                                                       //
  11 // Authors: Jim Thomas, Stefan Rossegger                                  //
  12 ////////////////////////////////////////////////////////////////////////////
  13
  14 #include "AliTPCCorrection.h"
  15
  16
  17 class AliTPCSpaceCharge : public AliTPCCorrection {
  18 public:
  19   AliTPCSpaceCharge();
  20   virtual ~AliTPCSpaceCharge();
  21
  22   // initialization and update functions
  23   virtual void Init();
  24   virtual void Update(const TTimeStamp &timeStamp);
  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:
  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: ???
  54
  55   Bool_t fInitLookUp;                  // flag to check it the Look Up table was created
  56
  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)
  59
  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
  64
  65   ClassDef(AliTPCSpaceCharge,0);
  66 };
  67
  68 #endif