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 » REC » AliTPCClusterParam

class AliTPCClusterParam: public TObject


  TPC cluster error, shape and charge parameterization as function
of drift length, and inclination angle

  Following notation is used in following
  Int_t dim 0 - y direction
            1 - z direction

  Int_t type 0 - short pads
             1 - medium pads
             2 - long pads
  Float_t z    - drift length

  Float_t angle - tangent of inclination angle at given dimension

  Implemented parameterization


  1. Resolution as function of drift length and inclination angle
     1.a) GetError0(Int_t dim, Int_t type, Float_t z, Float_t angle)
          Simple error parameterization as derived from analytical formula
          only linear term in drift length and angle^2
          The formula is valid only with precission +-5%
          Separate parameterization for differnt pad geometry
     1.b) GetError0Par
          Parabolic term correction - better precision

     1.c) GetError1 - JUST FOR Study
          Similar to GetError1
          The angular and diffusion effect is scaling with pad length
          common parameterization for different pad length

  2. Error parameterization using charge
     2.a) GetErrorQ
          GetError0+
          adding 1/Q component to diffusion and angluar part
     2.b) GetErrorQPar
          GetError0Par+
          adding 1/Q component to diffusion and angluar part
     2.c) GetErrorQParScaled - Just for study
          One parameterization for all pad shapes
          Smaller precission as previous one


  Example how to retrieve the paramterization:

      AliCDBManager::Instance()->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
      AliCDBManager::Instance()->SetRun(0)
      AliTPCClusterParam * param = AliTPCcalibDB::Instance()->GetClusterParam();

      //
      //
      AliTPCClusterParam::SetInstance(param);
      TF1 f1("f1","AliTPCClusterParam::SGetError0Par(1,0,x,0)",0,250);

Function Members (Methods)

public:

	AliTPCClusterParam()
	AliTPCClusterParam(const AliTPCClusterParam& param)
virtual	~AliTPCClusterParam()
void	TObject::AbstractMethod(const char* method) const
virtual void	TObject::AppendPad(Option_t* option = "")
virtual void	TObject::Browse(TBrowser* b)
static TClass*	Class()
virtual const char*	TObject::ClassName() const
virtual void	TObject::Clear(Option_t* = "")
virtual TObject*	TObject::Clone(const char* newname = "") const
virtual Int_t	TObject::Compare(const TObject* obj) const
virtual void	TObject::Copy(TObject& object) const
virtual void	TObject::Delete(Option_t* option = "")MENU
virtual Int_t	TObject::DistancetoPrimitive(Int_t px, Int_t py)
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)
virtual void	TObject::Fatal(const char* method, const char* msgfmt) const
virtual TObject*	TObject::FindObject(const char* name) const
virtual TObject*	TObject::FindObject(const TObject* obj) const
void	FitData(TTree* tree)
void	FitResol(TTree* tree)
void	FitResol0(TTree* tree, Int_t dim, Int_t type, Float_t* param0, Float_t* error)
void	FitResol0Par(TTree* tree, Int_t dim, Int_t type, Float_t* param0, Float_t* error)
void	FitResol1(TTree* tree, Int_t dim, Float_t* param0, Float_t* error)
void	FitResolQ(TTree* tree, Int_t dim, Int_t type, Float_t* param0, Float_t* error)
void	FitResolQPar(TTree* tree, Int_t dim, Int_t type, Float_t* param0, Float_t* error)
void	FitRMS(TTree* tree)
void	FitRMS0(TTree* tree, Int_t dim, Int_t type, Float_t* param0, Float_t* error)
void	FitRMS1(TTree* tree, Int_t dim, Float_t* param0, Float_t* error)
void	FitRMSQ(TTree* tree, Int_t dim, Int_t type, Float_t* param0, Float_t* error)
void	FitRMSSigma(TTree* tree, Int_t dim, Int_t type, Float_t* param0, Float_t* error)
static Double_t	GaussConvolution(Double_t x0, Double_t x1, Double_t k0, Double_t k1, Double_t s0, Double_t s1)
static Double_t	GaussConvolutionGamma4(Double_t x0, Double_t x1, Double_t k0, Double_t k1, Double_t s0, Double_t s1, Double_t tau)
static Double_t	GaussConvolutionTail(Double_t x0, Double_t x1, Double_t k0, Double_t k1, Double_t s0, Double_t s1, Double_t tau)
virtual Option_t*	TObject::GetDrawOption() const
static Long_t	TObject::GetDtorOnly()
Float_t	GetError0(Int_t dim, Int_t type, Float_t z, Float_t angle) const
Float_t	GetError0Par(Int_t dim, Int_t type, Float_t z, Float_t angle) const
Float_t	GetError1(Int_t dim, Int_t type, Float_t z, Float_t angle) const
Float_t	GetErrorQ(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean) const
Float_t	GetErrorQPar(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean) const
Float_t	GetErrorQParScaled(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean) const
virtual const char*	TObject::GetIconName() const
virtual const char*	TObject::GetName() const
virtual char*	TObject::GetObjectInfo(Int_t px, Int_t py) const
static Bool_t	TObject::GetObjectStat()
virtual Option_t*	TObject::GetOption() const
Double_t	GetQnormCorr(Int_t ipad, Int_t itype, Int_t corrType) const
TMatrixD*	GetQnormCorrMatrix()
const THnBase*	GetResolutionYMap() const
Float_t	GetRMS0(Int_t dim, Int_t type, Float_t z, Float_t angle) const
Float_t	GetRMS1(Int_t dim, Int_t type, Float_t z, Float_t angle) const
Float_t	GetRMSQ(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean) const
Float_t	GetRMSSigma(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean) const
Float_t	GetShapeFactor(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean, Float_t rmsL, Float_t rmsM) const
virtual const char*	TObject::GetTitle() const
virtual UInt_t	TObject::GetUniqueID() const
Float_t	GetWaveCorrection(Int_t Type, Float_t Z, Int_t QMax, Float_t Pad, Float_t angleY) const
const THnBase*	GetWaveCorrectionMap() const
virtual Bool_t	TObject::HandleTimer(TTimer* timer)
virtual ULong_t	TObject::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	TObject::Inspect() constMENU
static AliTPCClusterParam*	Instance()
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	TObject::IsSortable() const
Bool_t	TObject::IsZombie() const
virtual void	TObject::ls(Option_t* option = "") const
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)
AliTPCClusterParam&	operator=(const AliTPCClusterParam& param)
virtual void	TObject::Paint(Option_t* option = "")
Float_t	ParamS0Par(Int_t i, Int_t j, Int_t k) const
virtual void	TObject::Pop()
Float_t	PosCorrection(Int_t type, Int_t ipad, Float_t pad, Float_t time, Float_t z, Float_t sy2, Float_t sz2, Float_t qm)
TVectorD*&	PosYcor(Int_t ind)
TVectorD*&	PosZcor(Int_t ind)
virtual void	Print(Option_t* option = "") const
static Double_t	QmaxCorrection(Int_t sector, Int_t row, Float_t cpad, Float_t ctime, Float_t ky, Float_t kz, Float_t rmsy0, Float_t rmsz0, Float_t effLength = 0, Float_t effDiff = 1)
Float_t	Qnorm(Int_t ipad, Int_t itype, Float_t dr, Float_t ty, Float_t tz)
Float_t	QnormHis(Int_t ipad, Int_t itype, Float_t dr, Float_t ty, Float_t tz)
Float_t	QnormPos(Int_t ipad, Bool_t isMax, Float_t pad, Float_t time, Float_t z, Float_t sy2, Float_t sz2, Float_t qm, Float_t qt)
TVectorD*	QpadMnorm() const
TVectorD*	QpadTnorm() const
static Double_t	QtotCorrection(Int_t sector, Int_t row, Float_t cpad, Float_t ctime, Float_t ky, Float_t kz, Float_t rmsy0, Float_t rmsz0, Float_t qtot, Float_t thr, Float_t effLength = 0, Float_t effDiff = 1)
virtual Int_t	TObject::Read(const char* name)
virtual void	TObject::RecursiveRemove(TObject* obj)
void	TObject::ResetBit(UInt_t f)
void	ResetQnormCorr()
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)
virtual void	TObject::SetDrawOption(Option_t* option = "")MENU
static void	TObject::SetDtorOnly(void* obj)
void	SetInstance(AliTPCClusterParam *const param)
static void	TObject::SetObjectStat(Bool_t stat)
void	SetQnorm(Int_t ipad, Int_t itype, const TVectorD *const norm)
void	SetQnormCorr(Int_t ipad, Int_t itype, Int_t corrType, Float_t val, Int_t mode = 1)
void	SetResolutionYMap(THnBase* ResolutionYMap)
virtual void	TObject::SetUniqueID(UInt_t uid)
void	SetWaveCorrectionMap(THnBase* WaveCorrectionMap)
static Float_t	SGetError0(Int_t dim, Int_t type, Float_t z, Float_t angle)
static Float_t	SGetError0Par(Int_t dim, Int_t type, Float_t z, Float_t angle)
static Float_t	SGetError1(Int_t dim, Int_t type, Float_t z, Float_t angle)
static Float_t	SGetErrorQ(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean)
static Float_t	SGetErrorQPar(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean)
static Float_t	SGetErrorQParScaled(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean)
static Float_t	SGetRMS0(Int_t dim, Int_t type, Float_t z, Float_t angle)
static Float_t	SGetRMS1(Int_t dim, Int_t type, Float_t z, Float_t angle)
static Float_t	SGetRMSQ(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean)
static Float_t	SGetRMSSigma(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean)
static Float_t	SGetShapeFactor(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean, Float_t rmsL, Float_t rmsM)
static Float_t	SGetWaveCorrection(Int_t Type, Float_t Z, Int_t QMax, Float_t Pad, Float_t angleY)
virtual void	ShowMembers(TMemberInspector&)
static Float_t	SPosCorrection(Int_t type, Int_t ipad, Float_t pad, Float_t time, Float_t z, Float_t sy2, Float_t sz2, Float_t qm)
static Float_t	SQnorm(Int_t ipad, Int_t itype, Float_t dr, Float_t ty, Float_t tz)
static Float_t	SQnormHis(Int_t ipad, Int_t itype, Float_t dr, Float_t ty, Float_t tz)
static Float_t	SQnormPos(Int_t ipad, Bool_t isMax, Float_t pad, Float_t time, Float_t z, Float_t sy2, Float_t sz2, Float_t qm, Float_t qt)
virtual void	Streamer(TBuffer&)
void	StreamerNVirtual(TBuffer& ClassDef_StreamerNVirtual_b)
virtual void	TObject::SysError(const char* method, const char* msgfmt) const
void	Test(TTree* tree, const char* output = "TestClusterParam.root")
Bool_t	TObject::TestBit(UInt_t f) const
Int_t	TObject::TestBits(UInt_t f) const
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:

virtual void	TObject::DoError(int level, const char* location, const char* fmt, va_list va) const
void	TObject::MakeZombie()

Data Members

public:

enum TObject::EStatusBits {	kCanDelete
	kMustCleanup
	kObjInCanvas
	kIsReferenced
	kHasUUID
	kCannotPick
	kNoContextMenu
	kInvalidObject
};
enum TObject::[unnamed] {	kIsOnHeap
	kNotDeleted
	kZombie
	kBitMask
	kSingleKey
	kOverwrite
	kWriteDelete
};

protected:

Float_t	fErrorRMS0[2][3][4]	shape parameterization coeficients
Float_t	fErrorRMS1[2][5]	shape parameterization coeficients
Float_t	fErrorRMSQ[2][3][6]	shape parameterization coeficients
Float_t	fErrorRMSSys[2]	systematic relative error of the parametererization
Float_t	fErrorS0[2][3][4]	error parameterization coeficients
Float_t	fErrorS0Par[2][3][7]	error parameterization coeficients
Float_t	fErrorS1[2][4]	error parameterization coeficients
Float_t	fErrorSQ[2][3][6]	error parameterization coeficients
Float_t	fErrorSQPar[2][3][9]	error parameterization coeficients
Float_t	fParamRMS0[2][3][4]	shape parameterization coeficients
Float_t	fParamRMS1[2][5]	shape parameterization coeficients
Float_t	fParamRMSQ[2][3][6]	shape parameterization coeficients
Float_t	fParamS0[2][3][4]	error parameterization coeficients
Float_t	fParamS0Par[2][3][7]	error parameterization coeficients
Float_t	fParamS1[2][4]	error parameterization coeficients
Float_t	fParamSQ[2][3][6]	error parameterization coeficients
Float_t	fParamSQPar[2][3][9]	error parameterization coeficients
TVectorD*	fPosQMnorm[3]	q position normalization
TVectorD*	fPosQTnorm[3]	q position normalization
TVectorD*	fPosYcor[3]	position correction parameterization
TVectorD*	fPosZcor[3]	position correction parameterization
TObjArray*	fQNorm	q norm paramters
TMatrixD*	fQNormCorr	q norm correction for analytica correction
TObjArray*	fQNormHis	q norm correction for analytical correction
TVectorD*	fQpadMnorm	q pad normalization - Max charge
TVectorD*	fQpadTnorm	q pad normalization - Total charge
Float_t	fRMSSigmaFit[2][3][2]	mean value of the varation of RMS to RMS
Float_t	fRMSSigmaRatio[2][2]	mean value of the varation of RMS to RMS
Float_t	fRatio	ratio of values constibution to error
THnBase*	fResolutionYMap	Map of resolution in Y
THnBase*	fWaveCorrectionMap	dY with respect to the distance to the center of the pad
Bool_t	fWaveCorrectionMirroredAngle	flag is the Angle axis mirrored at 0
Bool_t	fWaveCorrectionMirroredPad	flag is the cog axis mirrored at 0.5
Bool_t	fWaveCorrectionMirroredZ	flag is the Z axis mirrored at 0
static AliTPCClusterParam*	fgInstance	! Instance of this class (singleton implementation)

Class Charts

Inheritance Chart:

TObject

←

AliTPCClusterParam

Function documentation

AliTPCClusterParam* Instance()

 Singleton implementation
 Returns an instance of this class, it is created if neccessary

AliTPCClusterParam()


 Default constructor

AliTPCClusterParam(const AliTPCClusterParam& param)

 copy constructor

AliTPCClusterParam & operator=(const AliTPCClusterParam& param)

 Assignment operator

~AliTPCClusterParam()

 destructor

void FitResol0(TTree* tree, Int_t dim, Int_t type, Float_t* param0, Float_t* error)

 Fit z - angular dependence of resolution

 Int_t dim=0, type=0;

void FitResol0Par(TTree* tree, Int_t dim, Int_t type, Float_t* param0, Float_t* error)

 Fit z - angular dependence of resolution

 Int_t dim=0, type=0;

void FitResol1(TTree* tree, Int_t dim, Float_t* param0, Float_t* error)

 Fit z - angular dependence of resolution - pad length scaling

 Int_t dim=0, type=0;

void FitResolQ(TTree* tree, Int_t dim, Int_t type, Float_t* param0, Float_t* error)

 Fit z - angular dependence of resolution - Q scaling

 Int_t dim=0, type=0;

void FitResolQPar(TTree* tree, Int_t dim, Int_t type, Float_t* param0, Float_t* error)

 Fit z - angular dependence of resolution - Q scaling  - parabolic correction

 Int_t dim=0, type=0;

void FitRMS0(TTree* tree, Int_t dim, Int_t type, Float_t* param0, Float_t* error)

 Fit z - angular dependence of resolution

 Int_t dim=0, type=0;

void FitRMS1(TTree* tree, Int_t dim, Float_t* param0, Float_t* error)

 Fit z - angular dependence of resolution - pad length scaling

 Int_t dim=0, type=0;

void FitRMSQ(TTree* tree, Int_t dim, Int_t type, Float_t* param0, Float_t* error)

 Fit z - angular dependence of resolution - Q scaling

 Int_t dim=0, type=0;

void FitRMSSigma(TTree* tree, Int_t dim, Int_t type, Float_t* param0, Float_t* error)

 Fit z - angular dependence of resolution - Q scaling

 Int_t dim=0, type=0;

Float_t GetError0(Int_t dim, Int_t type, Float_t z, Float_t angle) const

Float_t GetError0Par(Int_t dim, Int_t type, Float_t z, Float_t angle) const

Float_t GetError1(Int_t dim, Int_t type, Float_t z, Float_t angle) const

Float_t GetErrorQ(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean) const

Float_t GetErrorQPar(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean) const

Float_t GetErrorQParScaled(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean) const

Float_t GetRMS0(Int_t dim, Int_t type, Float_t z, Float_t angle) const

 calculate mean RMS of cluster - z,angle - parameters for each pad and dimension separatelly

Float_t GetRMS1(Int_t dim, Int_t type, Float_t z, Float_t angle) const

 calculate mean RMS of cluster - z,angle - pad length scalling

Float_t GetRMSQ(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean) const

 calculate mean RMS of cluster - z,angle, Q dependence

Float_t GetRMSSigma(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean) const

 calculates RMS of signal shape fluctuation

Float_t GetShapeFactor(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean, Float_t rmsL, Float_t rmsM) const

 calculates vallue - sigma distortion contribution

void FitData(TTree* tree)

 make fits for error param and shape param

void FitResol(TTree* tree)

void FitRMS(TTree* tree)

void Test(TTree* tree, const char* output = "TestClusterParam.root")

 Draw standard quality histograms to output file

void Print(Option_t* option = "") const

 Print param Information

Float_t Qnorm(Int_t ipad, Int_t itype, Float_t dr, Float_t ty, Float_t tz)

 get Q normalization
 type - 0 Qtot 1 Qmax
 ipad - 0 (0.75 cm) ,1 (1 cm), 2 (1.5 cm)

expession formula - TString *strq0 = toolkit.FitPlane(chain,"dedxQ.fElements[2]","dr++ty++tz++dr*ty++dr*tz++++dr*dr++ty*tz++ty^2++tz^2","IPad==0",chi2,npoints,param,covar,0,100000);

Float_t QnormHis(Int_t ipad, Int_t itype, Float_t dr, Float_t ty, Float_t tz)

 get Q normalization
 type - 0 Qtot 1 Qmax
 ipad - 0 (0.75 cm) ,1 (1 cm), 2 (1.5 cm)

void SetQnorm(Int_t ipad, Int_t itype, const TVectorD *const norm)

 set normalization

 type - 0 Qtot 1 Qmax
 ipad - 0 (0.75 cm) ,1 (1 cm), 2 (1.5 cm)

void ResetQnormCorr()

void SetQnormCorr(Int_t ipad, Int_t itype, Int_t corrType, Float_t val, Int_t mode = 1)

 ipad        - pad type
 itype       - 0- qtot 1-qmax
 corrType    - 0 - s0y corr     - eff. PRF corr
             - 1 - s0z corr     - eff. TRF corr
             - 2 - d0y          - eff. diffusion correction y
             - 3 - d0z          - eff. diffusion correction
             - 4 - eff length   - eff.length - wire pitch + x diffsion
             - 5 - pad type normalization

Double_t GetQnormCorr(Int_t ipad, Int_t itype, Int_t corrType) const

 see AliTPCClusterParam::SetQnormCorr

Float_t QnormPos(Int_t ipad, Bool_t isMax, Float_t pad, Float_t time, Float_t z, Float_t sy2, Float_t sz2, Float_t qm, Float_t qt)

 Make Q normalization as function of following parameters
 Fit parameters to be used in corresponding correction function extracted in the AliTPCclaibTracksGain - Taylor expansion
 1 - dp   - relative pad position
 2 - dt   - relative time position
 3 - di   - drift length (norm to 1);
 4 - dq0  - Tot/Max charge
 5 - dq1  - Max/Tot charge
 6 - sy   - sigma y - shape
 7 - sz   - sigma z - shape

Float_t PosCorrection(Int_t type, Int_t ipad, Float_t pad, Float_t time, Float_t z, Float_t sy2, Float_t sz2, Float_t qm)

Double_t GaussConvolution(Double_t x0, Double_t x1, Double_t k0, Double_t k1, Double_t s0, Double_t s1)

 2 D gaus convoluted with angular effect
 See in mathematica:
Simplify[Integrate[Exp[-(x0-k0*xd)*(x0-k0*xd)/(2*s0*s0)-(x1-k1*xd)*(x1-k1*xd)/(2*s1*s1)]/(s0*s1),{xd,-1/2,1/2}]]

TF1 f1("f1","AliTPCClusterParam::GaussConvolution(x,0,1,0,0.1,0.1)",-2,2)
TF2 f2("f2","AliTPCClusterParam::GaussConvolution(x,y,1,1,0.1,0.1)",-2,2,-2,2)

Double_t GaussConvolutionTail(Double_t x0, Double_t x1, Double_t k0, Double_t k1, Double_t s0, Double_t s1, Double_t tau)

 2 D gaus convoluted with angular effect and exponential tail in z-direction
 tail integrated numerically
 Integral normalized to one
 Mean at 0

 TF1 f1t("f1t","AliTPCClusterParam::GaussConvolutionTail(0,x,0,0,0.5,0.5,0.9)",-5,5)

Double_t GaussConvolutionGamma4(Double_t x0, Double_t x1, Double_t k0, Double_t k1, Double_t s0, Double_t s1, Double_t tau)

 2 D gaus convoluted with angular effect and exponential tail in z-direction
 tail integrated numerically
 Integral normalized to one
 Mean at 0

 TF1 f1g4("f1g4","AliTPCClusterParam::GaussConvolutionGamma4(0,x,0,0,0.5,0.2,1.6)",-5,5)
 TF2 f2g4("f2g4","AliTPCClusterParam::GaussConvolutionGamma4(y,x,0,0,0.5,0.2,1.6)",-5,5,-5,5)

Double_t QmaxCorrection(Int_t sector, Int_t row, Float_t cpad, Float_t ctime, Float_t ky, Float_t kz, Float_t rmsy0, Float_t rmsz0, Float_t effLength = 0, Float_t effDiff = 1)


 cpad      - pad (y) coordinate
 ctime     - time(z) coordinate
 ky        - dy/dx
 kz        - dz/dx
 rmsy0     - RF width in pad units
 rmsz0     - RF width in time bin  units
 effLength - contibution of PRF and diffusion
 effDiff   - overwrite diffusion

Double_t QtotCorrection(Int_t sector, Int_t row, Float_t cpad, Float_t ctime, Float_t ky, Float_t kz, Float_t rmsy0, Float_t rmsz0, Float_t qtot, Float_t thr, Float_t effLength = 0, Float_t effDiff = 1)


 cpad      - pad (y) coordinate
 ctime     - time(z) coordinate
 ky        - dy/dx
 kz        - dz/dx
 rmsy0     - RF width in pad units
 rmsz0     - RF width in time bin  units
 qtot      - the sum of signal in cluster - without thr correction
 thr       - threshold
 effLength - contibution of PRF and diffusion
 effDiff   - overwrite diffusion

void SetWaveCorrectionMap(THnBase* WaveCorrectionMap)

 Set Correction Map for Y

void SetResolutionYMap(THnBase* ResolutionYMap)

 Set Resolution Map for Y

Float_t GetWaveCorrection(Int_t Type, Float_t Z, Int_t QMax, Float_t Pad, Float_t angleY) const

 Correct Y cluster coordinate using a map

void SetInstance(AliTPCClusterParam *const param)

 Seting functions

{fgInstance = param;}

TMatrixD * GetQnormCorrMatrix()

{return fQNormCorr;}

Float_t SQnormPos(Int_t ipad, Bool_t isMax, Float_t pad, Float_t time, Float_t z, Float_t sy2, Float_t sz2, Float_t qm, Float_t qt)

{ return fgInstance->QnormPos(ipad,isMax,pad,time,z,sy2,sz2,qm,qt);;}

Float_t SPosCorrection(Int_t type, Int_t ipad, Float_t pad, Float_t time, Float_t z, Float_t sy2, Float_t sz2, Float_t qm)

 Error parameterization

{ return fgInstance->PosCorrection(type,ipad,pad,time,z,sy2,sz2,qm);}

const THnBase * GetWaveCorrectionMap() const

 Correction and resolution maps

{ return fWaveCorrectionMap; }

const THnBase * GetResolutionYMap() const

{ return fResolutionYMap; }

Float_t SGetWaveCorrection(Int_t Type, Float_t Z, Int_t QMax, Float_t Pad, Float_t angleY)

{return (fgInstance) ? fgInstance->GetWaveCorrection(Type,Z,QMax,Pad,angleY):0;}

Float_t SGetError0(Int_t dim, Int_t type, Float_t z, Float_t angle)

 static methods equivalents  - use instance of param object - useful for tree draw and TF2 visualization

Float_t SGetError0Par(Int_t dim, Int_t type, Float_t z, Float_t angle)

Float_t SGetError1(Int_t dim, Int_t type, Float_t z, Float_t angle)

Float_t SGetErrorQ(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean)

Float_t SGetErrorQPar(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean)

Float_t SGetErrorQParScaled(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean)

Float_t SGetRMS0(Int_t dim, Int_t type, Float_t z, Float_t angle)

Float_t SGetRMS1(Int_t dim, Int_t type, Float_t z, Float_t angle)

Float_t SGetRMSQ(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean)

Float_t SGetRMSSigma(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean)

Float_t SGetShapeFactor(Int_t dim, Int_t type, Float_t z, Float_t angle, Float_t Qmean, Float_t rmsL, Float_t rmsM)

Float_t SQnorm(Int_t ipad, Int_t itype, Float_t dr, Float_t ty, Float_t tz)

{return fgInstance->Qnorm(ipad, itype, dr,ty,tz);}

Float_t SQnormHis(Int_t ipad, Int_t itype, Float_t dr, Float_t ty, Float_t tz)

{return fgInstance->QnormHis(ipad, itype, dr,ty,tz);}

TVectorD*& PosYcor(Int_t ind)

{return fPosYcor[ind];}

TVectorD*& PosZcor(Int_t ind)

{return fPosZcor[ind];}

Float_t ParamS0Par(Int_t i, Int_t j, Int_t k) const

{return fParamS0Par[i][j][k];}

TVectorD* QpadTnorm() const

{return fQpadTnorm;}

TVectorD* QpadMnorm() const

{return fQpadMnorm;}