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 AliTPCCorrection *	GetVisualCorrection (Int_t position)

static AliTPCCorrection *	GetVisualCorrection (const char *corName)

static TObjArray *	GetVisualCorrections ()

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 TObjArray *	fgVisualCorrection =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

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:

TPC/TPCbase/AliTPCSpaceCharge.h
TPC/TPCbase/AliTPCSpaceCharge.cxx

Public Member Functions

Protected Member Functions

Private Types

Private Attributes

Additional Inherited Members

Detailed Description

Member Enumeration Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation