#include <AliTPCFCVoltError3D.h>

Inheritance diagram for AliTPCFCVoltError3D:

Public Member Functions
	AliTPCFCVoltError3D ()

virtual	~AliTPCFCVoltError3D ()

virtual Bool_t	AddCorrectionCompact (AliTPCCorrection *corr, Double_t weight)

virtual void	Init ()

virtual void	Update (const TTimeStamp &timeStamp)

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	SetRodVoltShiftA (Int_t rod, Float_t voltOffset, Bool_t doInit=kTRUE)

void	SetRodVoltShiftC (Int_t rod, Float_t voltOffset, Bool_t doInit=kTRUE)

Float_t	GetRodVoltShiftA (Int_t i) const

Float_t	GetRodVoltShiftC (Int_t i) const

void	SetRotatedClipVoltA (Int_t rod, Float_t voltOffset, Bool_t doInit=kTRUE)

void	SetRotatedClipVoltC (Int_t rod, Float_t voltOffset, Bool_t doInit=kTRUE)

Float_t	GetRotatedClipVoltA (Int_t i) const

Float_t	GetRotatedClipVoltC (Int_t i) const

void	SetCopperRodShiftA (Int_t rod, Float_t voltOffset, Bool_t doInit=kTRUE)

void	SetCopperRodShiftC (Int_t rod, Float_t voltOffset, Bool_t doInit=kTRUE)

Float_t	GetCopperRodShiftA (Int_t i) const

Float_t	GetCopperRodShiftC (Int_t i) const

void	InitFCVoltError3D ()

void	ForceInitFCVoltError3D ()

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

Public Member Functions inherited from AliTPCCorrection
	AliTPCCorrection ()

	AliTPCCorrection (const char name, const char title)

virtual	~AliTPCCorrection ()

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, Int_t type=0)

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, IntegrationType integrationType=kIntegral)

void	SetIsLocal (Bool_t isLocal)

Bool_t	IsLocal () const

Private Types
enum	{ kRows =257 }

enum	{ kColumns =129 }

enum	{ kPhiSlicesPerSector = 10 }

enum	{ kPhiSlices = 1+kPhiSlicesPerSector*3 }

enum	{ kIterations =100 }

Private Member Functions
	AliTPCFCVoltError3D (const AliTPCFCVoltError3D &)

AliTPCFCVoltError3D &	operator= (const AliTPCFCVoltError3D &)

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	fRodVoltShiftA [36]
	Rod (&strips) shift in Volt (40V~1mm) More...

Float_t	fRodVoltShiftC [36]
	Rod (&strips) shift in Volt (40V~1mm) More...

Float_t	fRotatedClipVoltA [2]
	rotated clips at HV rod More...

Float_t	fRotatedClipVoltC [2]
	rotated clips at HV rod More...

Float_t	fCopperRodShiftA [36]
	only Rod shift More...

Float_t	fCopperRodShiftC [36]
	only Rod shift More...

Bool_t	fInitLookUp
	flag to check if the Look Up table was created (SUM) More...

Bool_t	fInitLookUpBasic [6]
	! flag if the basic lookup was created (shifted Rod (IFC,OFC) or rotated clip (IFC,OFC)) More...

TMatrixF *	fLookUpErOverEz [kNPhi]
	Array to store electric field integral (int Er/Ez) More...

TMatrixF *	fLookUpEphiOverEz [kNPhi]
	Array to store electric field integral (int Er/Ez) More...

TMatrixF *	fLookUpDeltaEz [kNPhi]
	Array to store electric field integral (int Er/Ez) More...

TMatrixD *	fLookUpBasic1ErOverEz [kPhiSlices]
	! Array to store electric field integral (int Er/Ez) More...

TMatrixD *	fLookUpBasic1EphiOverEz [kPhiSlices]
	! Array to store electric field integral (int Ephi/Ez) More...

TMatrixD *	fLookUpBasic1DeltaEz [kPhiSlices]
	! Array to store electric field integral (int Ez) More...

TMatrixD *	fLookUpBasic2ErOverEz [kPhiSlices]
	! Array to store electric field integral More...

TMatrixD *	fLookUpBasic2EphiOverEz [kPhiSlices]
	! Array to store electric field integral More...

TMatrixD *	fLookUpBasic2DeltaEz [kPhiSlices]
	! Array to store electric field integral More...

TMatrixD *	fLookUpBasic3ErOverEz [kPhiSlices]
	! Array to store electric field integral More...

TMatrixD *	fLookUpBasic3EphiOverEz [kPhiSlices]
	! Array to store electric field integral More...

TMatrixD *	fLookUpBasic3DeltaEz [kPhiSlices]
	! Array to store electric field integral More...

TMatrixD *	fLookUpBasic4ErOverEz [kPhiSlices]
	! Array to store electric field integral More...

TMatrixD *	fLookUpBasic4EphiOverEz [kPhiSlices]
	! Array to store electric field integral More...

TMatrixD *	fLookUpBasic4DeltaEz [kPhiSlices]
	! Array to store electric field integral More...

TMatrixD *	fLookUpBasic5ErOverEz [kPhiSlices]
	! Array to store electric field integral More...

TMatrixD *	fLookUpBasic5EphiOverEz [kPhiSlices]
	! Array to store electric field integral More...

TMatrixD *	fLookUpBasic5DeltaEz [kPhiSlices]
	! Array to store electric field integral More...

TMatrixD *	fLookUpBasic6ErOverEz [kPhiSlices]
	! Array to store electric field integral More...

TMatrixD *	fLookUpBasic6EphiOverEz [kPhiSlices]
	! Array to store electric field integral More...

TMatrixD *	fLookUpBasic6DeltaEz [kPhiSlices]
	! Array to store electric field integral More...

Additional Inherited Members
Public Types inherited from AliTPCCorrection
enum	CompositionType { kParallel, kQueue }

enum	IntegrationType { kIntegral, kDifferential }

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)

static Double_t	GetCorrectionSector (Double_t sector, Double_t r, Double_t kZ, Int_t axisType, Int_t corrType=0)

static Double_t	GetDistortionSector (Double_t sector, Double_t r, Double_t kZ, Int_t axisType, Int_t corrType=0)

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...

IntegrationType	fIntegrationType
	Presentation of the underlying corrections, integrated, or differential. 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)/hera/alice/wiechula/calib/guiTrees. 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

AliTPCFCVoltError3D class.

The class calculates the space point distortions due to residual voltage errors on the Field Cage (FC) boundaries (rods and strips) of the TPC in 3D. It uses the Poisson relaxation technique in three dimension as implemented in the parent class.

Although the calculation is performed in 3D, the calculation time was significantly reduced by using certain symmetry conditions within the calculation.

The input parameters can be set via the functions (e.g.) SetRodVoltShift(rod,dV) where rod is the number of the rod on which the voltage offset dV is set. The corresponding shift in z direction would be $dz=dV/40$ with an opposite sign for the C side. The rods are numbered in anti-clock-wise direction starting at $=0$. Rods in the IFC are from 0 to 17, rods on the OFC go from 18 to 35. This convention is following the offline numbering scheme of the ROCs.

Different misalignment scenarios can be modeled:

A rotated clip scenario is only possible at two positions (Rod 11 at IFC, rod 3(+18) at OFC) and can be set via SetRotatedClipVolt. The key feature is that at the mentioned positions, the strip-ends were combined. At these positions, a systematic offset of one strip-end in respect to the other is possible.
A normal rod offset, where the strips follow the movement of the rod, can be set via SetRodVoltShift. It has a anti-mirrored signature in terms of distortions when compared to the rotated clip. This misalignment is possible at each single rod of the FC.
A simple rod offset, where the strips do not follow the shift, results in an even more localized distortion close to the rod. The difference to a rod shift, where the strips follow, is more dominant on the OFC. This effect can be set via the function SetCopperRodShift.

Picture from ROOT macro

Author: Jim Thomas, Stefan Rossegger

Date: 08/08/2010

Author: Jim Thomas, Stefan Rossegger

Definition at line 15 of file AliTPCFCVoltError3D.h.

Member Enumeration Documentation

anonymous enum

private

Enumerator
kRows

Definition at line 93 of file AliTPCFCVoltError3D.h.

anonymous enum

private

Enumerator
kColumns

Definition at line 94 of file AliTPCFCVoltError3D.h.

anonymous enum

private

Enumerator
kPhiSlicesPerSector

Definition at line 95 of file AliTPCFCVoltError3D.h.

anonymous enum

private

Enumerator
kPhiSlices

Definition at line 96 of file AliTPCFCVoltError3D.h.

anonymous enum

private

Enumerator
kIterations

Definition at line 97 of file AliTPCFCVoltError3D.h.

Constructor & Destructor Documentation

AliTPCFCVoltError3D::AliTPCFCVoltError3D ( )

Definition at line 68 of file AliTPCFCVoltError3D.cxx.

AliTPCFCVoltError3D::~AliTPCFCVoltError3D ( )

virtual

destructor

Definition at line 134 of file AliTPCFCVoltError3D.cxx.

AliTPCFCVoltError3D::AliTPCFCVoltError3D ( const AliTPCFCVoltError3D & )

private

Member Function Documentation

Bool_t AliTPCFCVoltError3D::AddCorrectionCompact	(	AliTPCCorrection *	corr,
		Double_t	weight
	)

virtual

Add correction and make them compact Assumptions:

origin of distortion/correction are additive
only correction ot the same type supported ()

Reimplemented from AliTPCCorrection.

Definition at line 172 of file AliTPCFCVoltError3D.cxx.

void AliTPCFCVoltError3D::ForceInitFCVoltError3D ( )

inline

Definition at line 61 of file AliTPCFCVoltError3D.h.

Referenced by GetCorrection().

Float_t AliTPCFCVoltError3D::GetC0 ( ) const

inline

Definition at line 31 of file AliTPCFCVoltError3D.h.

Float_t AliTPCFCVoltError3D::GetC1 ( ) const

inline

Definition at line 32 of file AliTPCFCVoltError3D.h.

Float_t AliTPCFCVoltError3D::GetCopperRodShiftA ( Int_t i ) const

inline

Definition at line 56 of file AliTPCFCVoltError3D.h.

Referenced by TestCorrection_AliTPCFCVoltError3DAddCorrectionCompact().

Float_t AliTPCFCVoltError3D::GetCopperRodShiftC ( Int_t i ) const

inline

Definition at line 57 of file AliTPCFCVoltError3D.h.

Referenced by TestCorrection_AliTPCFCVoltError3DAddCorrectionCompact().

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

protectedvirtual

Calculates the correction due e.g. residual voltage errors on the TPC boundaries

Reimplemented from AliTPCCorrection.

Definition at line 240 of file AliTPCFCVoltError3D.cxx.

Float_t AliTPCFCVoltError3D::GetRodVoltShiftA ( Int_t i ) const

inline

Definition at line 41 of file AliTPCFCVoltError3D.h.

Referenced by TestCorrection_AliTPCFCVoltError3DAddCorrectionCompact().

Float_t AliTPCFCVoltError3D::GetRodVoltShiftC ( Int_t i ) const

inline

Definition at line 42 of file AliTPCFCVoltError3D.h.

Referenced by TestCorrection_AliTPCFCVoltError3DAddCorrectionCompact().

Float_t AliTPCFCVoltError3D::GetRotatedClipVoltA ( Int_t i ) const

inline

Definition at line 48 of file AliTPCFCVoltError3D.h.

Float_t AliTPCFCVoltError3D::GetRotatedClipVoltC ( Int_t i ) const

inline

Definition at line 49 of file AliTPCFCVoltError3D.h.

void AliTPCFCVoltError3D::Init ( )

virtual

Initialization funtion

Reimplemented from AliTPCCorrection.

Definition at line 204 of file AliTPCFCVoltError3D.cxx.

void AliTPCFCVoltError3D::InitFCVoltError3D ( )

Initialization of the Lookup table which contains the solutions of the Dirichlet boundary problem Calculation of the single 3D-Poisson solver is done just if needed (see basic lookup tables in header file)

Definition at line 307 of file AliTPCFCVoltError3D.cxx.

Referenced by ForceInitFCVoltError3D(), GetCorrection(), Init(), LoadModels(), MakeEfieldCorrection(), RegisterAliTPCFCVoltError3D(), and RegisterAliTPCFCVoltError3DRodFCSideRadiusType().