Quick Links: ROOT Homepage Class Index Class Hierarchy

Source: header file source file inheritance tree (.pdf)

Sections: class description function members data members class charts

AliRoot » TPC » BASE » AliTPCComposedCorrection

class AliTPCComposedCorrection: public AliTPCCorrection


AliTPCComposedCorrection class

This class is creating a correction that is composed out of smaller
corrections.
There are two ways the sub-corrections can be combined into this one:
1. kParallel: All corrections are applied at the given position x and
the dx terms are summed up (this commutes).
2. kQueue: The corrections are called in order. The first one at the
given position x resulting in dx1, the second one is called at
the corrected position (x+dx1) resulting in dx2, the third one
is then called at position (x+dx1+dx2) and so forth. dx=dx1+dx2+...
is returned.
3. kQueueResidual: like kQueue with the exception that in case of
a positive weight the 'Distortion' is called and in case of a negative
weight the 'Correction' is called, where the absolute of the weight
    will be applied to the correction
For the inverse of the correction this is taken into account by reversing
the order the corrections are applied in the kQueue case (no issue for
kParallel).

date: 27/04/2010
Authors: Magnus Mager, Stefan Rossegger, Jim Thomas

Example usage:

AliMagF mag("mag","mag");
AliTPCExBBShape exb;             // B field shape distortions
exb.SetBField(&mag);

AliTPCExBTwist twist;            // ExB Twist distortions
twist.SetXTwist(0.001);

TObjArray cs;  cs.Add(&exb); cs.Add(&twist);

AliTPCComposedCorrection cc;
cc.SetCorrections(&cs);
cc.SetOmegaTauT1T2(wt,T1,T2);
cc.Print("DA");
cc.CreateHistoDRPhiinZR(0,100,100)->Draw("surf2");

Function Members (Methods)

public:

	AliTPCComposedCorrection()
	AliTPCComposedCorrection(TCollection* corrections, AliTPCComposedCorrection::CompositionType mode)
virtual	~AliTPCComposedCorrection()
void	TObject::AbstractMethod(const char* method) const
virtual Bool_t	AddCorrectionCompact(AliTPCCorrection* corr, Double_t weight)
static void	AliTPCCorrection::AddVisualCorrection(AliTPCCorrection* corr, Int_t position)
virtual void	TObject::AppendPad(Option_t* option = "")
virtual void	TObject::Browse(TBrowser* b)
static TClass*	Class()
virtual const char*	TObject::ClassName() const
virtual void	TNamed::Clear(Option_t* option = "")
virtual TObject*	TNamed::Clone(const char* newname = "") const
virtual Int_t	TNamed::Compare(const TObject* obj) const
virtual void	TNamed::Copy(TObject& named) const
void	AliTPCCorrection::CorrectPoint(Float_t* x, Short_t roc)
void	AliTPCCorrection::CorrectPoint(const Float_t* x, Short_t roc, Float_t* xp)
void	AliTPCCorrection::CorrectPointLocal(Float_t* x, Short_t roc)
TTree*	AliTPCCorrection::CreateDistortionTree(Double_t step = 5)
TH2F*	AliTPCCorrection::CreateHistoDRinXY(Float_t z = 10., Int_t nx = 100, Int_t ny = 100)
TH2F*	AliTPCCorrection::CreateHistoDRinZR(Float_t phi = 0., Int_t nZ = 100, Int_t nR = 100)
TH2F*	AliTPCCorrection::CreateHistoDRPhiinXY(Float_t z = 10., Int_t nx = 100, Int_t nphi = 100)
TH2F*	AliTPCCorrection::CreateHistoDRPhiinZR(Float_t phi = 0., Int_t nZ = 100, Int_t nR = 100)
TH2F*	AliTPCCorrection::CreateHistoDZinXY(Float_t z = 10., Int_t nx = 100, Int_t ny = 100)
TH2F*	AliTPCCorrection::CreateHistoDZinZR(Float_t phi = 0., Int_t nZ = 100, Int_t nR = 100)
virtual void	TObject::Delete(Option_t* option = "")MENU
virtual Int_t	TObject::DistancetoPrimitive(Int_t px, Int_t py)
void	AliTPCCorrection::DistortPoint(Float_t* x, Short_t roc)
void	AliTPCCorrection::DistortPoint(const Float_t* x, Short_t roc, Float_t* xp)
void	AliTPCCorrection::DistortPointLocal(Float_t* x, Short_t roc)
virtual void	TObject::Draw(Option_t* option = "")
virtual void	TObject::DrawClass() constMENU
virtual TObject*	TObject::DrawClone(Option_t* option = "") constMENU
virtual void	TObject::Dump() constMENU
virtual void	TObject::Error(const char* method, const char* msgfmt) const
virtual void	TObject::Execute(const char* method, const char* params, Int_t* error = 0)
virtual void	TObject::Execute(TMethod* method, TObjArray* params, Int_t* error = 0)
virtual void	TObject::ExecuteEvent(Int_t event, Int_t px, Int_t py)
void	AliTPCCorrection::FastSimDistortedVertex(Double_t* orgVertex, Int_t nTracks, AliESDVertex& aV, AliESDVertex& avOrg, AliESDVertex& cV, AliESDVertex& cvOrg, TTreeSRedirector *const pcstream, Double_t etaCuts)
virtual void	TObject::Fatal(const char* method, const char* msgfmt) const
virtual void	TNamed::FillBuffer(char*& buffer)
virtual TObject*	TObject::FindObject(const char* name) const
virtual TObject*	TObject::FindObject(const TObject* obj) const
AliExternalTrackParam*	AliTPCCorrection::FitDistortedTrack(AliExternalTrackParam& trackIn, Double_t refX, Int_t dir, TTreeSRedirector* pcstream)
virtual void	GetCorrection(const Float_t* x, const Short_t roc, Float_t* dx)
virtual void	AliTPCCorrection::GetCorrectionDz(const Float_t* x, Short_t roc, Float_t* dx, Float_t delta)
virtual void	AliTPCCorrection::GetCorrectionIntegralDz(const Float_t* x, Short_t roc, Float_t* dx, Float_t delta)
TCollection*	GetCorrections() const
virtual Float_t	AliTPCCorrection::GetCorrScaleFactor() const
static Double_t	AliTPCCorrection::GetCorrSector(Double_t sector, Double_t r, Double_t kZ, Int_t axisType, Int_t corrType = 0)
static Double_t	AliTPCCorrection::GetCorrXYZ(Double_t gx, Double_t gy, Double_t gz, Int_t axisType, Int_t corrType = 0)
static Double_t	AliTPCCorrection::GetCorrXYZDz(Double_t gx, Double_t gy, Double_t gz, Int_t axisType, Int_t corrType = 0, Double_t delta = 5)
static Double_t	AliTPCCorrection::GetCorrXYZIntegrateZ(Double_t gx, Double_t gy, Double_t gz, Int_t axisType, Int_t corrType = 0, Double_t delta = 5)
virtual void	GetDistortion(const Float_t* x, const Short_t roc, Float_t* dx)
virtual void	AliTPCCorrection::GetDistortionDz(const Float_t* x, Short_t roc, Float_t* dx, Float_t delta)
virtual void	AliTPCCorrection::GetDistortionIntegralDz(const Float_t* x, Short_t roc, Float_t* dx, Float_t delta)
static Double_t	AliTPCCorrection::GetDistXYZ(Double_t gx, Double_t gy, Double_t gz, Int_t axisType, Int_t corrType = 0)
static Double_t	AliTPCCorrection::GetDistXYZDz(Double_t gx, Double_t gy, Double_t gz, Int_t axisType, Int_t corrType = 0, Double_t delta = 5)
static Double_t	AliTPCCorrection::GetDistXYZIntegrateZ(Double_t gx, Double_t gy, Double_t gz, Int_t axisType, Int_t corrType = 0, Double_t delta = 5)
virtual Option_t*	TObject::GetDrawOption() const
static Long_t	TObject::GetDtorOnly()
virtual const char*	TObject::GetIconName() const
AliTPCComposedCorrection::CompositionType	GetMode() const
virtual const char*	TNamed::GetName() const
virtual char*	TObject::GetObjectInfo(Int_t px, Int_t py) const
static Bool_t	TObject::GetObjectStat()
virtual Option_t*	TObject::GetOption() const
virtual AliTPCCorrection*	GetSubCorrection(Int_t ipos)
virtual AliTPCCorrection*	GetSubCorrection(const char* cname)
virtual const char*	TNamed::GetTitle() const
virtual UInt_t	TObject::GetUniqueID() const
static AliTPCCorrection*	AliTPCCorrection::GetVisualCorrection(Int_t position)
const TVectorD*	GetWeights() const
virtual Bool_t	TObject::HandleTimer(TTimer* timer)
virtual ULong_t	TNamed::Hash() const
virtual void	TObject::Info(const char* method, const char* msgfmt) const
virtual Bool_t	TObject::InheritsFrom(const char* classname) const
virtual Bool_t	TObject::InheritsFrom(const TClass* cl) const
virtual void	Init()
virtual void	TObject::Inspect() constMENU
void	TObject::InvertBit(UInt_t f)
virtual TClass*	IsA() const
virtual Bool_t	TObject::IsEqual(const TObject* obj) const
virtual Bool_t	TObject::IsFolder() const
Bool_t	TObject::IsOnHeap() const
virtual Bool_t	TNamed::IsSortable() const
Bool_t	TObject::IsZombie() const
virtual void	TNamed::ls(Option_t* option = "") const
static void	AliTPCCorrection::MakeDistortionMap(THnSparse* his0, TTreeSRedirector* pcstream, const char* hname, Int_t run, Float_t refX, Int_t type, Int_t integ = 1)
static void	AliTPCCorrection::MakeDistortionMapCosmic(THnSparse* his0, TTreeSRedirector* pcstream, const char* hname, Int_t run, Float_t refX, Int_t type)
static void	AliTPCCorrection::MakeDistortionMapSector(THnSparse* his0, TTreeSRedirector* pcstream, const char* hname, Int_t run, Int_t type)
static void	AliTPCCorrection::MakeLaserDistortionTree(TTree* tree, TObjArray* corrArray, Int_t itype)
static void	AliTPCCorrection::MakeLaserDistortionTreeOld(TTree* tree, TObjArray* corrArray, Int_t itype)
static void	AliTPCCorrection::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	AliTPCCorrection::MakeTrackDistortionTree(TTree* tinput, Int_t dtype, Int_t ptype, const TObjArray* corrArray, Int_t step = 1, Int_t offset = 0, Bool_t debug = 0)
void	TObject::MayNotUse(const char* method) const
virtual Bool_t	TObject::Notify()
void	TObject::Obsolete(const char* method, const char* asOfVers, const char* removedFromVers) const
static void	TObject::operator delete(void* ptr)
static void	TObject::operator delete(void* ptr, void* vp)
static void	TObject::operator delete[](void* ptr)
static void	TObject::operator delete[](void* ptr, void* vp)
void*	TObject::operator new(size_t sz)
void*	TObject::operator new(size_t sz, void* vp)
void*	TObject::operator new[](size_t sz)
void*	TObject::operator new[](size_t sz, void* vp)
virtual void	TObject::Paint(Option_t* option = "")
virtual void	TObject::Pop()
virtual void	Print(Option_t* option = "") const
virtual Int_t	TObject::Read(const char* name)
virtual void	TObject::RecursiveRemove(TObject* obj)
void	TObject::ResetBit(UInt_t f)
virtual void	TObject::SaveAs(const char* filename = "", Option_t* option = "") constMENU
virtual void	TObject::SavePrimitive(basic_ostream<char,char_traits<char> >& out, Option_t* option = "")
void	TObject::SetBit(UInt_t f)
void	TObject::SetBit(UInt_t f, Bool_t set)
void	SetCorrections(const TCollection* corrections)
virtual void	AliTPCCorrection::SetCorrScaleFactor(Float_t)
virtual void	TObject::SetDrawOption(Option_t* option = "")MENU
static void	TObject::SetDtorOnly(void* obj)
void	SetMode(AliTPCComposedCorrection::CompositionType mode)
virtual void	TNamed::SetName(const char* name)MENU
virtual void	TNamed::SetNameTitle(const char* name, const char* title)
static void	TObject::SetObjectStat(Bool_t stat)
virtual void	SetOmegaTauT1T2(Float_t omegaTau, Float_t t1, Float_t t2)
virtual void	TNamed::SetTitle(const char* title = "")MENU
virtual void	TObject::SetUniqueID(UInt_t uid)
void	SetWeights(TVectorD* weights)
virtual void	ShowMembers(TMemberInspector&)
virtual Int_t	TNamed::Sizeof() const
void	AliTPCCorrection::StoreInOCDB(Int_t startRun, Int_t endRun, const char* comment = 0)
virtual void	Streamer(TBuffer&)
void	StreamerNVirtual(TBuffer& ClassDef_StreamerNVirtual_b)
virtual void	TObject::SysError(const char* method, const char* msgfmt) const
Bool_t	TObject::TestBit(UInt_t f) const
Int_t	TObject::TestBits(UInt_t f) const
virtual void	Update(const TTimeStamp& timeStamp)
virtual void	TObject::UseCurrentStyle()
virtual void	TObject::Warning(const char* method, const char* msgfmt) const
virtual Int_t	TObject::Write(const char* name = 0, Int_t option = 0, Int_t bufsize = 0)
virtual Int_t	TObject::Write(const char* name = 0, Int_t option = 0, Int_t bufsize = 0) const

protected:

TH2F*	AliTPCCorrection::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)
virtual void	TObject::DoError(int level, const char* location, const char* fmt, va_list va) const
Double_t	AliTPCCorrection::Interpolate(const Double_t* xArray, const Double_t* yArray, Int_t order, Double_t x)
Float_t	AliTPCCorrection::Interpolate(const Double_t* xArray, const Float_t* yArray, Int_t order, Double_t x)
void	AliTPCCorrection::Interpolate2DEdistortion(Int_t order, Double_t r, Double_t z, const Double_t er[][kNR], Double_t& erValue)
Double_t	AliTPCCorrection::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)
Float_t	AliTPCCorrection::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)
void	AliTPCCorrection::Interpolate3DEdistortion(Int_t order, Double_t r, Float_t phi, Double_t z, const Double_t er[][kNPhi][kNR], const Double_t ephi[][kNPhi][kNR], const Double_t ez[][kNPhi][kNR], Double_t& erValue, Double_t& ephiValue, Double_t& ezValue)
Double_t	AliTPCCorrection::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)
Float_t	AliTPCCorrection::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)
Bool_t	AliTPCCorrection::IsLocal() const
virtual Int_t	AliTPCCorrection::IsPowerOfTwo(Int_t i) const
void	TObject::MakeZombie()
void	AliTPCCorrection::PoissonRelaxation2D(TMatrixD& arrayV, TMatrixD& chargeDensity, TMatrixD& arrayErOverEz, TMatrixD& arrayDeltaEz, Int_t rows, Int_t columns, Int_t iterations, Bool_t rocDisplacement = kTRUE)
void	AliTPCCorrection::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	AliTPCCorrection::Search(Int_t n, const Double_t* xArray, Double_t x, Int_t& low)
void	AliTPCCorrection::SetIsLocal(Bool_t isLocal)

private:

	AliTPCComposedCorrection(const AliTPCComposedCorrection&)
AliTPCComposedCorrection&	operator=(const AliTPCComposedCorrection&)

Data Members

public:

enum CompositionType {	kParallel
	kQueue
	kQueueResidual
};
enum AliTPCCorrection::CompositionType {	kParallel
	kQueue
};
enum AliTPCCorrection::[unnamed] {	kNR
	kNPhi
	kNZ
};
enum TObject::EStatusBits {	kCanDelete
	kMustCleanup
	kObjInCanvas
	kIsReferenced
	kHasUUID
	kCannotPick
	kNoContextMenu
	kInvalidObject
};
enum TObject::[unnamed] {	kIsOnHeap
	kNotDeleted
	kZombie
	kBitMask
	kSingleKey
	kOverwrite
	kWriteDelete
};

protected:

Int_t	AliTPCCorrection::fILow
Bool_t	AliTPCCorrection::fIsLocal	switch to indicate that the distortion is a local vector drphi/dz, dr/dz
Int_t	AliTPCCorrection::fJLow
Int_t	AliTPCCorrection::fKLow	variable to help in the interpolation
TString	TNamed::fName	object identifier
Double_t	AliTPCCorrection::fT1	tensor term of wt - T1
Double_t	AliTPCCorrection::fT2	tensor term of wt - T2
TString	TNamed::fTitle	object title
static TObjArray*	AliTPCCorrection::fgVisualCorrection	array of orrection for visualization
static const Double_t	AliTPCCorrection::fgkCathodeV	Cathode Voltage (volts)
static const Double_t	AliTPCCorrection::fgkEM	charge/mass in [C/kg]
static const Double_t	AliTPCCorrection::fgkGG	Gating Grid voltage (volts)
static const Double_t	AliTPCCorrection::fgkIFCRadius	Mean Radius of the Inner Field Cage ( 82.43 min, 83.70 max) (cm)
static const Double_t	AliTPCCorrection::fgkOFCRadius	Mean Radius of the Outer Field Cage (252.55 min, 256.45 max) (cm)
Double_t	AliTPCCorrection::fgkPhiList[181]	points in the phi direction (for the lookup table)
Double_t	AliTPCCorrection::fgkRList[72]	points in the radial direction (for the lookup table)
static const Double_t	AliTPCCorrection::fgkTPCZ0	nominal gating grid position
Double_t	AliTPCCorrection::fgkZList[166]	points in the z direction (for the lookup table)
static const Double_t	AliTPCCorrection::fgkZOffSet	Offset from CE: calculate all distortions closer to CE as if at this point
static const Double_t	AliTPCCorrection::fgkdvdE	[cm/V] drift velocity dependency on the E field (from Magboltz for NeCO2N2 at standard environment)
static const Double_t	AliTPCCorrection::fgke0	vacuum permittivity [A·s/(V·m)]

private:

TCollection*	fCorrections	The corrections this one is composed of.
AliTPCComposedCorrection::CompositionType	fMode	The way to apply the corrections (see general class documentation)
TVectorD*	fWeights	optional vector with weights - used for fit benchmarking

Class Charts

Inheritance Chart:

TNamed

←

AliTPCCorrection

←

AliTPCComposedCorrection

Function documentation

AliTPCComposedCorrection()

 default constructor

AliTPCComposedCorrection(TCollection* corrections, AliTPCComposedCorrection::CompositionType mode)

 Constructor that defines the set of corrections, this one is composed of.

~AliTPCComposedCorrection()

 destructor

Bool_t AddCorrectionCompact(AliTPCCorrection* corr, Double_t weight)

 Add correction  - better name needed (left/right) - for the moment I assumme they commute
 Why not to just use array of corrections - CPU consideration
 Assumptions:
  - origin of distortion/correction are additive
  - corrections/distortion are small   and they commute
  - only correction ot the same type supported

AliTPCCorrection * GetSubCorrection(Int_t ipos)

AliTPCCorrection * GetSubCorrection(const char* cname)

void GetCorrection(const Float_t* x, const Short_t roc, Float_t* dx)

 This applies all correction and the specified manner (see general
 class description for details).

void GetDistortion(const Float_t* x, const Short_t roc, Float_t* dx)

 This applies all distortions and the specified manner (see general
 class descxiption for details).

void Print(Option_t* option = "") const

 Print function to check which correction classes are used
 option=="d" prints details regarding the setted magnitude
 option=="a" prints the C0 and C1 coefficents for calibration purposes

void Init()

 Initialization funtion (not used at the moment)

void Update(const TTimeStamp& timeStamp)

 Update function

void SetOmegaTauT1T2(Float_t omegaTau, Float_t t1, Float_t t2)

 Gives the possibility to set the OmegaTau plus Tensor corrections T1 and T2 (effective omega Tau)
 to each subcorrection (since they might become event specific due to changing drift velocity)

 The omegaTau comes idealy from the Database, since it is a function of drift velocity, B and E field
 e.g. omegaTau  = -10.0 * Bz * vdrift / Ez ; // with Bz in kG and Ez in V/cm
      omegaTau  = -0.325 for Bz=5kG, Ez=400V/cm and vdrift = 2.6cm/muSec
 The T1 and T2 tensors were measured in a dedicated calibration run

 Note: overwrites previously set values!

AliTPCComposedCorrection()

TCollection* GetCorrections() const

{return fCorrections;}

void SetCorrections(const TCollection* corrections)

{fCorrections=(TCollection*)corrections;}

CompositionType GetMode() const

{return fMode;}

void SetMode(AliTPCComposedCorrection::CompositionType mode)

{fMode=mode;}

void SetWeights(TVectorD* weights)

{fWeights= (TVectorD*) weights->Clone();}

const TVectorD * GetWeights() const

{return fWeights;}