Inheritance diagram for AliTPCclusterFast:

Public Member Functions
	AliTPCclusterFast ()

void	Init ()

virtual	~AliTPCclusterFast ()

void	SetParam (Float_t mnprim, Float_t diff, Float_t diffL, Int_t padrow, Float_t y, Float_t z, Float_t ky, Float_t kz, Float_t yCenter, Float_t zCenter)

void	Digitize ()

Double_t	GetQtot (Float_t gain, Float_t thr, Float_t noise, Bool_t rounding=kTRUE, Bool_t addPedestal=kTRUE)

Double_t	GetQmax (Float_t gain, Float_t thr, Float_t noise, Bool_t rounding=kTRUE, Bool_t addPedestal=kTRUE)

Double_t	GetQmaxCorr (Float_t rmsy0, Float_t rmsz0)

Double_t	GetQtotCorr (Float_t rmsy0, Float_t rmsz0, Float_t gain, Float_t thr)

Double_t	GetClusterProperties (TVectorF &param, Int_t addOverlap, Float_t gain=0.8, Float_t thr=2, Float_t noise=0.7, Bool_t rounding=kTRUE, Bool_t addPedestal=kTRUE, Int_t skipSample=0)

Double_t	GetCOG (Int_t returnType, Int_t addOverlap, Float_t gain=0.8, Float_t thr=2, Float_t noise=0.7, Bool_t rounding=kTRUE, Bool_t addPedestal=kTRUE, Int_t skipSample=0)

Double_t	GetCOGHit (Int_t returnType)

Float_t	UnfoldCluster (Float_t &meani, Float_t &meanj, Float_t &sumu, Float_t &overlap, Int_t addOverlap, Float_t gain=0.8, Float_t thr=2, Float_t noise=0.7, Bool_t rounding=kTRUE, Bool_t addPedestal=kTRUE, Int_t skipSample=0)

Float_t	GetCOGUnfolded (Int_t returnType, Int_t addOverlap=kTRUE, Float_t gain=0.8, Float_t thr=2, Float_t noise=0.7, Bool_t rounding=kTRUE, Bool_t addPedestal=kTRUE, Int_t skipSample=0)

Bool_t	MakeDigitization (Int_t addOverlap, Float_t gain, Float_t thr, Float_t noise, Bool_t rounding, Bool_t addPedestal, Int_t skipSample)

Double_t	GetExpectedRMS (Int_t dim)

Double_t	GetNsec ()

Int_t	GetYMaxBin ()

Int_t	GetZMaxBin ()

TMatrixF *	GetRawDigits ()

Float_t	GetRawDigit (Int_t i, Int_t j)

Float_t	GetDigit (Int_t i, Int_t j)

Float_t	GetDigitsRawMax ()

Float_t	GetDigitsMax ()

Static Public Member Functions
static void	GenerElectrons (AliTPCclusterFast cl0, AliTPCclusterFast clm, AliTPCclusterFast *clp)

static void	InitFormulas ()

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	GaussExpConvolution (Double_t x0, Double_t s0, Double_t t1)

static Double_t	GaussGamma4 (Double_t x, Double_t s0, Double_t p1)

static Double_t	Gamma4 (Double_t x, Double_t p0, Double_t p1)

static Double_t	Gamma4Norm (Double_t x)

static Double_t	GEMPRF (Double_t x, Double_t sigma)

static void	SetMetadata (TTree *tree)

static Bool_t	UnitTest ()

Public Attributes
Float_t	fMNprim
	mean number of primary electrons More...

Int_t	fNprim
	mean number of primary electrons More...

Int_t	fNtot
	total number of electrons More...

Float_t	fQtot
	total charge - Gas gain flucuation taken into account More...

Float_t	fDiff
	diffusion sigma More...

Float_t	fDiffLong
	diffusion sigma longitudinal direction More...

Int_t	fPadRow
	cluster pad row number More...

Float_t	fY
	y ideal position - center bin More...

Float_t	fZ
	z ideal position - center bin More...

Float_t	fYCenterBin
	y center bin More...

Float_t	fZCenterBin
	z center bin More...

Int_t	fYMaxBin

Int_t	fZMaxBin
	y maximum position as position as calculated in the MakeDigitization (should be 0 in mean ) More...

Float_t	fAngleY
	z maximum position as position as calculated in the MakeDigitization (should be 0 in mean ) More...

Float_t	fAngleZ
	z angle - tan z More...

AliTPCclusterFast *	fOverlapCluster

TVectorF	fSec
	number of secondary electrons More...

TVectorF	fPosY
	position y for each electron More...

TVectorF	fPosZ
	position z for each electron More...

TVectorF	fGain
	gg for each electron More...

TVectorF	fStatY

TVectorF	fStatZ

TMatrixF	fDigits
	response matrix (ideal signal without noise, trheshold rounding, pile-up) More...

TMatrixF	fRawDigits

Float_t	fPRFRMS
	"on the fly caclualted digits matrix" for current version of setup - gain, noise, thr. More...

Float_t	fTRFRMS
	pad respons function width More...

Static Public Attributes
static TF1 *	fPRF =0
	time rsponsefunction width More...

static TF1 *	fTRF =0
	Time response function. More...

static Float_t	fgZSamplingFactor =2

static Int_t	fgDebugLevel =0

Detailed Description

Code for the fast MC simulation of the TPC response

Microscopic simulation includes all features from the FULL MC simulation:

primary ionization
secondary inization
diffusion
gas gain with fluctuation
pad and time response fucntion

Usage:

fast sensitivity studies - dEdx
development of new algorithm
resolution studies
NEW: double track resolutuion studies - are currently implemented

Documentation for the reconstruction part - we will refer to the CHEP paper - (http://arxiv.org/pdf/physics/0306108.pdf

Definition at line 73 of file AliTPCclusterFast.cxx.

Constructor & Destructor Documentation

AliTPCclusterFast::AliTPCclusterFast ( )

Definition at line 593 of file AliTPCclusterFast.cxx.

Referenced by AliTPCtrackFast::MakeTrack().

AliTPCclusterFast::~AliTPCclusterFast ( )

virtual

Definition at line 641 of file AliTPCclusterFast.cxx.

Member Function Documentation

void AliTPCclusterFast::Digitize ( )

Clear digits

Definition at line 719 of file AliTPCclusterFast.cxx.

Referenced by AliTPCtrackFast::MakeTrack().

Double_t AliTPCclusterFast::Gamma4	(	Double_t	x,
		Double_t	p0,
		Double_t	p1
	)

static

Gamma 4 Time response function of ALTRO Gamma 4 Time response function of ALTRO

Definition at line 991 of file AliTPCclusterFast.cxx.

Referenced by Gamma4Norm(), GetDigitsMax(), and InitFormulas().

Double_t AliTPCclusterFast::Gamma4Norm ( Double_t x )

static

Gamma 4 Time response function of ALTRO with hardwired normalization parameters X in bin unit=100 ns Maximum of the TRF ==1 - sum of all time bins ~ 1.9

Definition at line 1001 of file AliTPCclusterFast.cxx.

Referenced by GetDigitsMax(), and InitFormulas().

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

static

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","AliTPCclusterFast::GaussConvolution(x,0,1,0,0.1,0.1)",-2,2) TF2 f2("f2","AliTPCclusterFast::GaussConvolution(x,y,1,1,0.1,0.1)",-2,2,-2,2)

Definition at line 936 of file AliTPCclusterFast.cxx.

Referenced by GetDigitsMax(), GetQmaxCorr(), and GetQtotCorr().

Double_t AliTPCclusterFast::GaussExpConvolution	(	Double_t	x0,
		Double_t	s0,
		Double_t	t1
	)

static

2 D gaus convoluted with exponential Integral nomalized to 1 See in mathematica: Simplify[Integrate[Exp[-(x0-x1)*(x0-x1)/(2*s0*s0)]*Exp[-x1*t1],{x1,0,Infinity}]] TF1 fgexp("fgexp","AliTPCclusterFast::GaussExpConvolution(x,0.5,1)",-2,2)

Definition at line 965 of file AliTPCclusterFast.cxx.

Referenced by GetDigitsMax().

Double_t AliTPCclusterFast::GaussGamma4	(	Double_t	x,
		Double_t	s0,
		Double_t	p1
	)

static

Gamma 4 Time response function of ALTRO convoluted with Gauss Simplify[Integrate[Exp[-(x0-x1)*(x0-x1)/(2*s0*s0)]*Exp[-4*x1/p1]*(x/p1)^4/s0,{x1,0,Infinity}]] TF1 fgg4("fgg4","AliTPCclusterFast::GaussGamma4(x,0.5,0.5)",-2,2)

Definition at line 1023 of file AliTPCclusterFast.cxx.

Referenced by GetDigitsMax().

Double_t AliTPCclusterFast::GEMPRF	(	Double_t	x,
		Double_t	sigma
	)

static

Definition at line 1013 of file AliTPCclusterFast.cxx.

Referenced by GetDigitsMax(), and InitFormulas().

void AliTPCclusterFast::GenerElectrons	(	AliTPCclusterFast *	cl0,
		AliTPCclusterFast *	clm,
		AliTPCclusterFast *	clp
	)

static

Definition at line 670 of file AliTPCclusterFast.cxx.

Referenced by AliTPCtrackFast::MakeTrack().

Double_t AliTPCclusterFast::GetClusterProperties	(	TVectorF &	param,
		Int_t	addOverlap,
		Float_t	gain = `0.8`,
		Float_t	thr = `2`,
		Float_t	noise = `0.7`,
		Bool_t	rounding = `kTRUE`,
		Bool_t	addPedestal = `kTRUE`,
		Int_t	skipSample = `0`
	)

Definition at line 828 of file AliTPCclusterFast.cxx.

Referenced by GetCOG().

Double_t AliTPCclusterFast::GetCOG	(	Int_t	returnType,
		Int_t	addOverlap,
		Float_t	gain = `0.8`,
		Float_t	thr = `2`,
		Float_t	noise = `0.7`,
		Bool_t	rounding = `kTRUE`,
		Bool_t	addPedestal = `kTRUE`,
		Int_t	skipSample = `0`
	)

Definition at line 910 of file AliTPCclusterFast.cxx.

Double_t AliTPCclusterFast::GetCOGHit ( Int_t returnType )

Definition at line 919 of file AliTPCclusterFast.cxx.

Float_t AliTPCclusterFast::GetCOGUnfolded	(	Int_t	returnType,
		Int_t	addOverlap = `kTRUE`,
		Float_t	gain = `0.8`,
		Float_t	thr = `2`,
		Float_t	noise = `0.7`,
		Bool_t	rounding = `kTRUE`,
		Bool_t	addPedestal = `kTRUE`,
		Int_t	skipSample = `0`
	)

Definition at line 1036 of file AliTPCclusterFast.cxx.

Float_t AliTPCclusterFast::GetDigit	(	Int_t	i,
		Int_t	j
	)

inline

Definition at line 102 of file AliTPCclusterFast.cxx.

Float_t AliTPCclusterFast::GetDigitsMax ( )

inline

Definition at line 104 of file AliTPCclusterFast.cxx.

Float_t AliTPCclusterFast::GetDigitsRawMax ( )

inline

Definition at line 103 of file AliTPCclusterFast.cxx.

Double_t AliTPCclusterFast::GetExpectedRMS ( Int_t dim )

Definition at line 896 of file AliTPCclusterFast.cxx.

Double_t AliTPCclusterFast::GetNsec ( )

Generate number of secondary electrons copy of procedure implemented in geant

Definition at line 656 of file AliTPCclusterFast.cxx.

Referenced by GenerElectrons().

Double_t AliTPCclusterFast::GetQmax	(	Float_t	gain,
		Float_t	thr,
		Float_t	noise,
		Bool_t	rounding = `kTRUE`,
		Bool_t	addPedestal = `kTRUE`
	)

Definition at line 769 of file AliTPCclusterFast.cxx.

Referenced by AliTPCtrackFast::CookdEdxDmax(), and GetQtotCorr().

Double_t AliTPCclusterFast::GetQmaxCorr	(	Float_t	rmsy0,
		Float_t	rmsz0
	)

Gaus distribution convolueted with rectangular Gaus width sy and sz is determined by RF width and diffusion Integral of Q is equal 1 Q max is calculated at position fY,fX

Definition at line 787 of file AliTPCclusterFast.cxx.

Referenced by AliTPCtrackFast::CookdEdxDmax(), and GetQtotCorr().

Double_t AliTPCclusterFast::GetQtot	(	Float_t	gain,
		Float_t	thr,
		Float_t	noise,
		Bool_t	rounding = `kTRUE`,
		Bool_t	addPedestal = `kTRUE`
	)

Definition at line 753 of file AliTPCclusterFast.cxx.

Referenced by AliTPCtrackFast::CookdEdxDtot(), and GetQtotCorr().

Double_t AliTPCclusterFast::GetQtotCorr	(	Float_t	rmsy0,
		Float_t	rmsz0,
		Float_t	gain,
		Float_t	thr
	)

Calculates the fraction of the charge over threshol to total charge The response function

Definition at line 799 of file AliTPCclusterFast.cxx.

Referenced by AliTPCtrackFast::CookdEdxDtot().

Float_t AliTPCclusterFast::GetRawDigit	(	Int_t	i,
		Int_t	j
	)

inline

Definition at line 101 of file AliTPCclusterFast.cxx.

TMatrixF* AliTPCclusterFast::GetRawDigits ( )

inline

Definition at line 100 of file AliTPCclusterFast.cxx.

Int_t AliTPCclusterFast::GetYMaxBin ( )

inline

Definition at line 98 of file AliTPCclusterFast.cxx.

Int_t AliTPCclusterFast::GetZMaxBin ( )

inline

Definition at line 99 of file AliTPCclusterFast.cxx.

void AliTPCclusterFast::Init ( )

reset all counters

< diffusion sigma

< diffusion sigma longitudinal direction

< y ideal position - center bin

< z ideal position - center bin

< y center bin

< z center bin

< y angle - tan(y)

< z angle - tan z

Definition at line 604 of file AliTPCclusterFast.cxx.

Referenced by AliTPCtrackFast::MakeTrack().

void AliTPCclusterFast::InitFormulas ( )

static

Definition at line 234 of file AliTPCclusterFast.cxx.

Referenced by GetDigitsMax().

Bool_t AliTPCclusterFast::MakeDigitization	(	Int_t	addOverlap,
		Float_t	gain,
		Float_t	thr,
		Float_t	noise,
		Bool_t	rounding,
		Bool_t	addPedestal,
		Int_t	skipSample
	)

Definition at line 1049 of file AliTPCclusterFast.cxx.

Referenced by UnfoldCluster().

void AliTPCclusterFast::SetMetadata ( TTree * tree )

static

Definition at line 1226 of file AliTPCclusterFast.cxx.

Referenced by GetDigitsMax().

void AliTPCclusterFast::SetParam	(	Float_t	mnprim,
		Float_t	diff,
		Float_t	diffL,
		Int_t	padrow,
		Float_t	y,
		Float_t	z,
		Float_t	ky,
		Float_t	kz,
		Float_t	yCenter,
		Float_t	zCenter
	)

Definition at line 645 of file AliTPCclusterFast.cxx.

Referenced by AliTPCtrackFast::MakeTrack().

Float_t AliTPCclusterFast::UnfoldCluster	(	Float_t &	meani,
		Float_t &	meanj,
		Float_t &	sumu,
		Float_t &	overlap,
		Int_t	addOverlap,
		Float_t	gain = `0.8`,
		Float_t	thr = `2`,
		Float_t	noise = `0.7`,
		Bool_t	rounding = `kTRUE`,
		Bool_t	addPedestal = `kTRUE`,
		Int_t	skipSample = `0`
	)