AliRoot Core  v5-06-30 (35d6c57)
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
AliTPCSpaceCharge Class Reference

The class calculates the space point distortions due to a rotational. More...

#include <AliTPCSpaceCharge.h>

Inheritance diagram for AliTPCSpaceCharge:

Public Member Functions

 AliTPCSpaceCharge ()
 
virtual ~AliTPCSpaceCharge ()
 
virtual void Init ()
 
virtual void Update (const TTimeStamp &timeStamp)
 
void SetInputSpaceCharge (TH2 *hisSpaceCharge)
 
const TH2 * GetInputSpaceCharge ()
 
virtual void SetOmegaTauT1T2 (Float_t omegaTau, Float_t t1, Float_t t2)
 
void SetC0C1 (Float_t c0, Float_t c1)
 
Float_t GetC0 () const
 
Float_t GetC1 () const
 
void SetCorrectionFactor (Float_t correctionFactor)
 
Float_t GetCorrectionFactor () const
 
void InitSpaceChargeDistortion ()
 
virtual void Print (const Option_t *option="") const
 
- Public Member Functions inherited from AliTPCCorrection
 AliTPCCorrection ()
 
 AliTPCCorrection (const char *name, const char *title)
 
virtual ~AliTPCCorrection ()
 
virtual Bool_t AddCorrectionCompact (AliTPCCorrection *corr, Double_t weight)
 
void CorrectPoint (Float_t x[], Short_t roc)
 
void CorrectPointLocal (Float_t x[], Short_t roc)
 
void CorrectPoint (const Float_t x[], Short_t roc, Float_t xp[])
 
virtual void GetCorrectionDz (const Float_t x[], Short_t roc, Float_t dx[], Float_t delta)
 
virtual void GetCorrectionIntegralDz (const Float_t x[], Short_t roc, Float_t dx[], Float_t delta)
 
void DistortPoint (Float_t x[], Short_t roc)
 
void DistortPointLocal (Float_t x[], Short_t roc)
 
void DistortPoint (const Float_t x[], Short_t roc, Float_t xp[])
 
virtual void GetDistortion (const Float_t x[], Short_t roc, Float_t dx[])
 
virtual void GetDistortionDz (const Float_t x[], Short_t roc, Float_t dx[], Float_t delta)
 
virtual void GetDistortionIntegralDz (const Float_t x[], Short_t roc, Float_t dx[], Float_t delta)
 
virtual void SetCorrScaleFactor (Float_t)
 
virtual Float_t GetCorrScaleFactor () const
 
virtual void Print (Option_t *option="") const
 
TH2F * CreateHistoDRinXY (Float_t z=10., Int_t nx=100, Int_t ny=100)
 
TH2F * CreateHistoDRPhiinXY (Float_t z=10., Int_t nx=100, Int_t nphi=100)
 
TH2F * CreateHistoDZinXY (Float_t z=10., Int_t nx=100, Int_t ny=100)
 
TH2F * CreateHistoDRinZR (Float_t phi=0., Int_t nZ=100, Int_t nR=100)
 
TH2F * CreateHistoDRPhiinZR (Float_t phi=0., Int_t nZ=100, Int_t nR=100)
 
TH2F * CreateHistoDZinZR (Float_t phi=0., Int_t nZ=100, Int_t nR=100)
 
TTree * CreateDistortionTree (Double_t step=5)
 
AliExternalTrackParam * FitDistortedTrack (AliExternalTrackParam &trackIn, Double_t refX, Int_t dir, TTreeSRedirector *pcstream)
 
void StoreInOCDB (Int_t startRun, Int_t endRun, const char *comment=0)
 
void FastSimDistortedVertex (Double_t orgVertex[3], Int_t nTracks, AliESDVertex &aV, AliESDVertex &avOrg, AliESDVertex &cV, AliESDVertex &cvOrg, TTreeSRedirector *const pcstream, Double_t etaCuts)
 

Protected Member Functions

virtual void GetCorrection (const Float_t x[], const Short_t roc, Float_t dx[])
 
- Protected Member Functions inherited from AliTPCCorrection
TH2F * CreateTH2F (const char *name, const char *title, const char *xlabel, const char *ylabel, const char *zlabel, Int_t nbinsx, Double_t xlow, Double_t xup, Int_t nbinsy, Double_t ylow, Double_t yup)
 
void Interpolate2DEdistortion (Int_t order, Double_t r, Double_t z, const Double_t er[kNZ][kNR], Double_t &erValue)
 
void Interpolate3DEdistortion (Int_t order, Double_t r, Float_t phi, Double_t z, const Double_t er[kNZ][kNPhi][kNR], const Double_t ephi[kNZ][kNPhi][kNR], const Double_t ez[kNZ][kNPhi][kNR], Double_t &erValue, Double_t &ephiValue, Double_t &ezValue)
 
Double_t Interpolate2DTable (Int_t order, Double_t x, Double_t y, Int_t nx, Int_t ny, const Double_t xv[], const Double_t yv[], const TMatrixD &array)
 
Double_t Interpolate3DTable (Int_t order, Double_t x, Double_t y, Double_t z, Int_t nx, Int_t ny, Int_t nz, const Double_t xv[], const Double_t yv[], const Double_t zv[], TMatrixD **arrayofArrays)
 
Double_t Interpolate (const Double_t xArray[], const Double_t yArray[], Int_t order, Double_t x)
 
void Search (Int_t n, const Double_t xArray[], Double_t x, Int_t &low)
 
Float_t Interpolate2DTable (Int_t order, Double_t x, Double_t y, Int_t nx, Int_t ny, const Double_t xv[], const Double_t yv[], const TMatrixF &array)
 
Float_t Interpolate3DTable (Int_t order, Double_t x, Double_t y, Double_t z, Int_t nx, Int_t ny, Int_t nz, const Double_t xv[], const Double_t yv[], const Double_t zv[], TMatrixF **arrayofArrays)
 
Float_t Interpolate (const Double_t xArray[], const Float_t yArray[], Int_t order, Double_t x)
 
virtual Int_t IsPowerOfTwo (Int_t i) const
 
void PoissonRelaxation2D (TMatrixD &arrayV, TMatrixD &chargeDensity, TMatrixD &arrayErOverEz, TMatrixD &arrayDeltaEz, Int_t rows, Int_t columns, Int_t iterations, Bool_t rocDisplacement=kTRUE)
 
void PoissonRelaxation3D (TMatrixD **arrayofArrayV, TMatrixD **arrayofChargeDensities, TMatrixD **arrayofEroverEz, TMatrixD **arrayofEPhioverEz, TMatrixD **arrayofEz, Int_t rows, Int_t columns, Int_t phislices, Float_t deltaphi, Int_t iterations, Int_t summetry, Bool_t rocDisplacement=kTRUE)
 
void SetIsLocal (Bool_t isLocal)
 
Bool_t IsLocal () const
 

Private Types

enum  { kRows =257 }
 
enum  { kColumns =129 }
 
enum  { kIterations =100 }
 

Private Attributes

Float_t fC0
 coefficient C0 (compare Jim Thomas's notes for definitions) More...
 
Float_t fC1
 coefficient C1 (compare Jim Thomas's notes for definitions) More...
 
Float_t fCorrectionFactor
 Space Charge Correction factor in comparison to initialized. More...
 
TH2 * fSpaceChargeHistogram
 Histogram with the input space charge histogram - used as an optional input. More...
 
Bool_t fInitLookUp
 flag to check it the Look Up table was created More...
 
Double_t fLookUpErOverEz [kNZ][kNR]
 Array to store electric field integral (int Er/Ez) More...
 
Double_t fLookUpDeltaEz [kNZ][kNR]
 Array to store electric field integral (int Delta Ez) More...
 

Additional Inherited Members

- Public Types inherited from AliTPCCorrection
enum  CompositionType { kParallel, kQueue }
 
- Static Public Member Functions inherited from AliTPCCorrection
static void MakeDistortionMap (THnSparse *his0, TTreeSRedirector *pcstream, const char *hname, Int_t run, Float_t refX, Int_t type, Int_t integ=1)
 
static void MakeDistortionMapCosmic (THnSparse *his0, TTreeSRedirector *pcstream, const char *hname, Int_t run, Float_t refX, Int_t type)
 
static void MakeDistortionMapSector (THnSparse *his0, TTreeSRedirector *pcstream, const char *hname, Int_t run, Int_t type)
 
static void MakeTrackDistortionTree (TTree *tinput, Int_t dtype, Int_t ptype, const TObjArray *corrArray, Int_t step=1, Int_t offset=0, Bool_t debug=0)
 
static void MakeSectorDistortionTree (TTree *tinput, Int_t dtype, Int_t ptype, const TObjArray *corrArray, Int_t step=1, Int_t offset=0, Bool_t debug=0)
 
static void MakeLaserDistortionTreeOld (TTree *tree, TObjArray *corrArray, Int_t itype)
 
static void MakeLaserDistortionTree (TTree *tree, TObjArray *corrArray, Int_t itype)
 
static void AddVisualCorrection (AliTPCCorrection *corr, Int_t position)
 
static AliTPCCorrectionGetVisualCorrection (Int_t position)
 
static AliTPCCorrectionGetVisualCorrection (const char *corName)
 
static TObjArrayGetVisualCorrections ()
 
static Double_t GetCorrSector (Double_t sector, Double_t r, Double_t kZ, Int_t axisType, Int_t corrType=0)
 
static Double_t GetCorrXYZ (Double_t gx, Double_t gy, Double_t gz, Int_t axisType, Int_t corrType=0)
 
static Double_t GetCorrXYZDz (Double_t gx, Double_t gy, Double_t gz, Int_t axisType, Int_t corrType=0, Double_t delta=5)
 
static Double_t GetCorrXYZIntegrateZ (Double_t gx, Double_t gy, Double_t gz, Int_t axisType, Int_t corrType=0, Double_t delta=5)
 
static Double_t GetDistXYZ (Double_t gx, Double_t gy, Double_t gz, Int_t axisType, Int_t corrType=0)
 
static Double_t GetDistXYZDz (Double_t gx, Double_t gy, Double_t gz, Int_t axisType, Int_t corrType=0, Double_t delta=5)
 
static Double_t GetDistXYZIntegrateZ (Double_t gx, Double_t gy, Double_t gz, Int_t axisType, Int_t corrType=0, Double_t delta=5)
 
- Protected Types inherited from AliTPCCorrection
enum  { kNR = 72 }
 
enum  { kNPhi = 18*10+1 }
 
enum  { kNZ = 166 }
 
- Protected Attributes inherited from AliTPCCorrection
Double_t fgkRList [kNR]
 points in the radial direction (for the lookup table) More...
 
Double_t fgkPhiList [kNPhi]
 points in the phi direction (for the lookup table) More...
 
Double_t fgkZList [kNZ]
 points in the z direction (for the lookup table) More...
 
Int_t fILow
 
Int_t fJLow
 
Int_t fKLow
 variable to help in the interpolation More...
 
Double_t fT1
 tensor term of wt - T1 More...
 
Double_t fT2
 tensor term of wt - T2 More...
 
Bool_t fIsLocal
 switch to indicate that the distortion is a local vector drphi/dz, dr/dz More...
 
- Static Protected Attributes inherited from AliTPCCorrection
static const Double_t fgkTPCZ0 = 249.7
 nominal gating grid position More...
 
static const Double_t fgkIFCRadius = 83.5
 Mean Radius of the Inner Field Cage ( 82.43 min, 83.70 max) (cm) More...
 
static const Double_t fgkOFCRadius = 254.5
 Mean Radius of the Outer Field Cage (252.55 min, 256.45 max) (cm) More...
 
static const Double_t fgkZOffSet = 0.2
 Offset from CE: calculate all distortions closer to CE as if at this point. More...
 
static const Double_t fgkCathodeV = -100000.0
 Cathode Voltage (volts) More...
 
static const Double_t fgkGG = -70.0
 Gating Grid voltage (volts) More...
 
static const Double_t fgkdvdE = 0.0024
 [cm/V] drift velocity dependency on the E field (from Magboltz for NeCO2N2 at standard environment) More...
 
static const Double_t fgkEM = -1.602176487e-19/9.10938215e-31
 charge/mass in [C/kg] More...
 
static const Double_t fgke0 = 8.854187817e-12
 vacuum permittivity [A·s/(V·m)] More...
 
static TObjArrayfgVisualCorrection =0
 array of orrection for visualization More...
 

Detailed Description

The class calculates the space point distortions due to a rotational.

symmetric space charge distribution with the TPC drift volume.

The class uses the PoissonRelaxation2D to calculate the resulting electrical field inhomogeneities in the (r,z)-plane. Then, the Langevin-integral formalism is used to calculate the space point distortions.

The class assumes, that the distortions scales linearly with the magnitude of the space charge distribution $(r,z)$. The in here assumed distribution is $$(r,z) = {(A-B\,z)}{r^2} $$ wherein the factors A and B scale with the event multiplicity and the interaction rate.

The scaling factor can be set via the function SetCorrectionFactor. An example of the shape of the distortions is given below.

MI modification - 22.05.2013 As an optional input the Space charge histogram RZ is used in case it is provided

  • using the SetInputSpaceCharge function
    AliTPCSpaceCharge_cxx_82e9c78.png
    Picture from ROOT macro
    Author
    Jim Thomas, Stefan Rossegger
    Date
    23/08/2010
    Author
    Jim Thomas, Stefan Rossegger

Definition at line 15 of file AliTPCSpaceCharge.h.

Member Enumeration Documentation

anonymous enum
private
Enumerator
kRows 

Definition at line 61 of file AliTPCSpaceCharge.h.

anonymous enum
private
Enumerator
kColumns 

Definition at line 62 of file AliTPCSpaceCharge.h.

anonymous enum
private
Enumerator
kIterations 

Definition at line 63 of file AliTPCSpaceCharge.h.

Constructor & Destructor Documentation

AliTPCSpaceCharge::AliTPCSpaceCharge ( )

Definition at line 59 of file AliTPCSpaceCharge.cxx.

AliTPCSpaceCharge::~AliTPCSpaceCharge ( )
virtual

default destructor

Definition at line 70 of file AliTPCSpaceCharge.cxx.

Member Function Documentation

Float_t AliTPCSpaceCharge::GetC0 ( ) const
inline

Definition at line 34 of file AliTPCSpaceCharge.h.

Float_t AliTPCSpaceCharge::GetC1 ( ) const
inline

Definition at line 35 of file AliTPCSpaceCharge.h.

void AliTPCSpaceCharge::GetCorrection ( const Float_t  x[],
const Short_t  roc,
Float_t  dx[] 
)
protectedvirtual

Calculates the correction due the Space Charge effect within the TPC drift volume

Reimplemented from AliTPCCorrection.

Definition at line 114 of file AliTPCSpaceCharge.cxx.

Float_t AliTPCSpaceCharge::GetCorrectionFactor ( ) const
inline

Definition at line 39 of file AliTPCSpaceCharge.h.

const TH2* AliTPCSpaceCharge::GetInputSpaceCharge ( )
inline

Definition at line 24 of file AliTPCSpaceCharge.h.

void AliTPCSpaceCharge::Init ( )
virtual

Initialization funtion

Reimplemented from AliTPCCorrection.

Definition at line 77 of file AliTPCSpaceCharge.cxx.

void AliTPCSpaceCharge::InitSpaceChargeDistortion ( )

Initialization of the Lookup table which contains the solutions of the poisson problem

Definition at line 168 of file AliTPCSpaceCharge.cxx.

Referenced by GetCorrection(), and Init().

void AliTPCSpaceCharge::Print ( const Option_t *  option = "") const
virtual

Print function to check the settings of the boundary vectors option=="a" prints the C0 and C1 coefficents for calibration purposes

Definition at line 262 of file AliTPCSpaceCharge.cxx.

void AliTPCSpaceCharge::SetC0C1 ( Float_t  c0,
Float_t  c1 
)
inline

Definition at line 33 of file AliTPCSpaceCharge.h.

void AliTPCSpaceCharge::SetCorrectionFactor ( Float_t  correctionFactor)
inline

Definition at line 38 of file AliTPCSpaceCharge.h.

void AliTPCSpaceCharge::SetInputSpaceCharge ( TH2 *  hisSpaceCharge)
inline

Definition at line 23 of file AliTPCSpaceCharge.h.

virtual void AliTPCSpaceCharge::SetOmegaTauT1T2 ( Float_t  omegaTau,
Float_t  t1,
Float_t  t2 
)
inlinevirtual

Reimplemented from AliTPCCorrection.

Definition at line 28 of file AliTPCSpaceCharge.h.

Referenced by Init(), and Update().

void AliTPCSpaceCharge::Update ( const TTimeStamp &  timeStamp)
virtual

Update function

Reimplemented from AliTPCCorrection.

Definition at line 94 of file AliTPCSpaceCharge.cxx.

Member Data Documentation

Float_t AliTPCSpaceCharge::fC0
private

coefficient C0 (compare Jim Thomas's notes for definitions)

Definition at line 49 of file AliTPCSpaceCharge.h.

Referenced by GetC0(), GetCorrection(), Print(), SetC0C1(), and SetOmegaTauT1T2().

Float_t AliTPCSpaceCharge::fC1
private

coefficient C1 (compare Jim Thomas's notes for definitions)

Definition at line 50 of file AliTPCSpaceCharge.h.

Referenced by GetC1(), GetCorrection(), Print(), SetC0C1(), and SetOmegaTauT1T2().

Float_t AliTPCSpaceCharge::fCorrectionFactor
private

Space Charge Correction factor in comparison to initialized.

Definition at line 51 of file AliTPCSpaceCharge.h.

Referenced by GetCorrection(), GetCorrectionFactor(), Print(), and SetCorrectionFactor().

Bool_t AliTPCSpaceCharge::fInitLookUp
private

flag to check it the Look Up table was created

Definition at line 56 of file AliTPCSpaceCharge.h.

Referenced by GetCorrection(), InitSpaceChargeDistortion(), and Print().

Double_t AliTPCSpaceCharge::fLookUpDeltaEz[kNZ][kNR]
private

Array to store electric field integral (int Delta Ez)

Definition at line 59 of file AliTPCSpaceCharge.h.

Referenced by GetCorrection(), and InitSpaceChargeDistortion().

Double_t AliTPCSpaceCharge::fLookUpErOverEz[kNZ][kNR]
private

Array to store electric field integral (int Er/Ez)

Definition at line 58 of file AliTPCSpaceCharge.h.

Referenced by GetCorrection(), and InitSpaceChargeDistortion().

TH2* AliTPCSpaceCharge::fSpaceChargeHistogram
private

Histogram with the input space charge histogram - used as an optional input.

Definition at line 55 of file AliTPCSpaceCharge.h.

Referenced by GetInputSpaceCharge(), InitSpaceChargeDistortion(), and SetInputSpaceCharge().


The documentation for this class was generated from the following files: