#include <AliFMDDensityCalculator.h>

Inheritance diagram for AliFMDDensityCalculator:

Classes
struct	RingHistos

Public Types
enum	{ kPhiNoCorrect, kPhiCorrectNch, kPhiCorrectELoss }

Public Member Functions
	AliFMDDensityCalculator ()

	AliFMDDensityCalculator (const char *name)

	AliFMDDensityCalculator (const AliFMDDensityCalculator &o)

virtual	~AliFMDDensityCalculator ()

AliFMDDensityCalculator &	operator= (const AliFMDDensityCalculator &o)

virtual void	SetupForData (const TAxis &etaAxis)

virtual Bool_t	Calculate (const AliESDFMD &fmd, AliForwardUtil::Histos &hists, Bool_t lowFlux, Double_t cent=-1, const TVector3 &ip=TVector3(1024, 1024, 0))

virtual void	Terminate (const TList dir, TList output, Int_t nEvents)

virtual void	CreateOutputObjects (TList *dir)

void	SetDebug (Int_t dbg=1)

void	SetDoTiming (Bool_t enable=true)

void	SetMaxParticles (UShort_t m)

void	SetUsePoisson (Bool_t u)

void	SetRecalculatePhi (Bool_t use)

void	SetUsePhiAcceptance (UShort_t u=kPhiCorrectNch)

void	SetLumping (Int_t eta, Int_t phi)

void	SetMinQuality (UShort_t cut=10)

void	SetMaxOutliers (Double_t ratio=0.10)

void	SetOutlierCut (Double_t cut=0.50)

Double_t	GetMultCut (UShort_t d, Char_t r, Double_t eta, Bool_t errors=true) const

Double_t	GetMultCut (UShort_t d, Char_t r, Int_t ieta, Bool_t errors=true) const

UShort_t	GetMinQuality () const

void	Print (Option_t *option="") const

AliFMDMultCuts &	GetCuts ()

void	SetCuts (const AliFMDMultCuts &c)

Static Public Attributes
static const char *	fgkFolderName

Protected Member Functions
Int_t	FindMaxWeight (const AliFMDCorrELossFit *cor, UShort_t d, Char_t r, Int_t iEta) const

Int_t	FindMaxWeight (const AliFMDCorrELossFit *cor, UShort_t d, Char_t r, Double_t iEta) const

void	CacheMaxWeights (const TAxis &axis)

Int_t	GetMaxWeight (UShort_t d, Char_t r, Int_t iEta) const

Int_t	GetMaxWeight (UShort_t d, Char_t r, Float_t eta) const

virtual Float_t	NParticles (Float_t mult, UShort_t d, Char_t r, Float_t eta, Bool_t lowFlux) const

virtual Float_t	Correction (UShort_t d, Char_t r, UShort_t t, Float_t eta, Bool_t lowFlux) const

virtual Float_t	AcceptanceCorrection (Char_t r, UShort_t t) const

virtual TH1D *	GenerateAcceptanceCorrection (Char_t r) const

virtual Bool_t	CheckOutlier (Double_t eloss, Double_t poisson, Double_t cut=0.5) const

RingHistos *	GetRingHistos (UShort_t d, Char_t r) const

	ClassDef (AliFMDDensityCalculator, 14)

Protected Attributes
TList	fRingHistos

TH1D *	fSumOfWeights

TH1D *	fWeightedSum

TH1D *	fCorrections

UShort_t	fMaxParticles

Bool_t	fUsePoisson

UShort_t	fUsePhiAcceptance

TH1D *	fAccI

TH1D *	fAccO

TArrayI	fFMD1iMax

TArrayI	fFMD2iMax

TArrayI	fFMD2oMax

TArrayI	fFMD3iMax

TArrayI	fFMD3oMax

TH2D *	fMaxWeights

TH2D *	fLowCuts

Int_t	fEtaLumping

Int_t	fPhiLumping

Int_t	fDebug

AliFMDMultCuts	fCuts

Bool_t	fRecalculatePhi

UShort_t	fMinQuality

AliForwardUtil::Histos	fCache

Bool_t	fDoTiming

TProfile *	fHTiming

Double_t	fMaxOutliers

Double_t	fOutlierCut

Detailed Description

This class calculates the inclusive charged particle density in each for the 5 FMD rings.

Input:

AliESDFMD object possibly corrected for sharing

Output:

5 RingHistos objects - each with a number of vertex dependent 2D histograms of the inclusive charge particle density

Corrections used:

AliFMDAnaCalibEnergyDistribution
AliFMDDoubleHitCorrection
AliFMDDeadCorrection

Definition at line 47 of file AliFMDDensityCalculator.h.

Member Enumeration Documentation

anonymous enum

How to correct for the missing phi coverage at the corners of the sensors

Enumerator
kPhiNoCorrect	No correction
kPhiCorrectNch	Correct the calculated number charged particles
kPhiCorrectELoss	Correct the energy loss

Definition at line 55 of file AliFMDDensityCalculator.h.

Constructor & Destructor Documentation

AliFMDDensityCalculator::AliFMDDensityCalculator ( )

Constructor

AliFMDDensityCalculator::AliFMDDensityCalculator ( const char * name )

Constructor

Parameters

name	Name of object

Definition at line 71 of file AliFMDDensityCalculator.cxx.

AliFMDDensityCalculator::AliFMDDensityCalculator ( const AliFMDDensityCalculator & o )

Copy constructor

Parameters

o	Object to copy from

Definition at line 138 of file AliFMDDensityCalculator.cxx.

AliFMDDensityCalculator::~AliFMDDensityCalculator ( )

virtual

Destructor

Definition at line 182 of file AliFMDDensityCalculator.cxx.

Member Function Documentation

Float_t AliFMDDensityCalculator::AcceptanceCorrection	(	Char_t	r,
		UShort_t	t
	)		const

protectedvirtual

Get the acceptance correction for strip t in an ring of type r

Parameters

r	Ring type ('I' or 'O')
t	Strip number

Returns: Inverse acceptance correction

Definition at line 1018 of file AliFMDDensityCalculator.cxx.

Referenced by Calculate(), and Correction().

void AliFMDDensityCalculator::CacheMaxWeights ( const TAxis & axis )

protected

Find the max weights and cache them

Parameters

axis	Default $\eta$ axis from parent task

Definition at line 713 of file AliFMDDensityCalculator.cxx.

Referenced by SetupForData().

Bool_t AliFMDDensityCalculator::Calculate	(	const AliESDFMD &	fmd,
		AliForwardUtil::Histos &	hists,
		Bool_t	lowFlux,
		Double_t	cent = `-1`,
		const TVector3 &	ip = `TVector3(1024,1024,0)`
	)

virtual

Do the calculations

Parameters

fmd	AliESDFMD object (possibly) corrected for sharing
hists	Histogram cache
lowFlux	Low flux flag.
cent	Centrality
ip	Coordinates of interaction point

Returns: true on successs

Definition at line 354 of file AliFMDDensityCalculator.cxx.

Referenced by AliForwardMultiplicityTask::Event(), AliForwardMCMultiplicityTask::Event(), and AliForwardQATask::Event().

Bool_t AliFMDDensityCalculator::CheckOutlier	(	Double_t	eloss,
		Double_t	poisson,
		Double_t	cut = `0.5`
	)		const

protectedvirtual

Check if, for a given region, whether this is an outlier

The condition for an outlier event are

\[ |N_{ch}^{Poisson} - N_{ch}^{\Delta}| / N_{ch}^{\Delta} > c \]

Parameters

eloss	$ N_{ch}^{\Delta}$ - number of charged particles
poisson	$ N_{ch}^{Poisson}$ - number of charged particles
cut	$ c$ - the cut

Returns: true if the region reflects an outlier event

Definition at line 650 of file AliFMDDensityCalculator.cxx.

Referenced by Calculate().

AliFMDDensityCalculator::ClassDef	(	AliFMDDensityCalculator	,
		14
	)

protected

Float_t AliFMDDensityCalculator::Correction	(	UShort_t	d,
		Char_t	r,
		UShort_t	t,
		Float_t	eta,
		Bool_t	lowFlux
	)		const

protectedvirtual

Get the inverse correction factor. This consist of

acceptance correction (corners of sensors)
double hit correction (for low-flux events)
dead strip correction

Parameters

d	Detector
r	Ring
t	Strip
eta	Pseudo-rapidity
lowFlux	Low-flux flag

Returns: the correction factor

Definition at line 890 of file AliFMDDensityCalculator.cxx.

Referenced by AliFMDMCDensityCalculator::CalculateMC().

void AliFMDDensityCalculator::CreateOutputObjects ( TList * dir )

virtual

Output diagnostic histograms to directory

Parameters

dir	List to write in

Reimplemented in AliFMDMCDensityCalculator.

Definition at line 1083 of file AliFMDDensityCalculator.cxx.

Referenced by AliForwardQATask::Book(), AliForwardMultiplicityBase::Book(), and AliFMDMCDensityCalculator::CreateOutputObjects().

Int_t AliFMDDensityCalculator::FindMaxWeight	(	const AliFMDCorrELossFit *	cor,
		UShort_t	d,
		Char_t	r,
		Int_t	iEta
	)		const

protected

Find the max weight to use for FMDdr in eta bin iEta

Parameters

cor	Correction
d	Detector
r	Ring
iEta	Eta bin

Returns: The maximum weight

Definition at line 660 of file AliFMDDensityCalculator.cxx.

Referenced by CacheMaxWeights().

Int_t AliFMDDensityCalculator::FindMaxWeight	(	const AliFMDCorrELossFit *	cor,
		UShort_t	d,
		Char_t	r,
		Double_t	iEta
	)		const

protected

Find the max weight to use for FMDdr in eta eta

Parameters

cor	Correction
d	Detector
r	Ring
iEta	iEta

Returns: The maximum weight

Definition at line 686 of file AliFMDDensityCalculator.cxx.

TH1D * AliFMDDensityCalculator::GenerateAcceptanceCorrection ( Char_t r ) const

protectedvirtual

Generate the acceptance corrections

Parameters

r	Ring to generate for

Returns: Newly allocated histogram of acceptance corrections

Definition at line 939 of file AliFMDDensityCalculator.cxx.

Referenced by AliFMDDensityCalculator().

AliFMDMultCuts& AliFMDDensityCalculator::GetCuts ( )

inline

Get the cuts used

Returns: Reference to cuts object

Definition at line 254 of file AliFMDDensityCalculator.h.

Referenced by AliForwardQATask::PreData().

Int_t AliFMDDensityCalculator::GetMaxWeight	(	UShort_t	d,
		Char_t	r,
		Int_t	iEta
	)		const

protected

Find the (cached) maximum weight for FMDdr in $\eta$ bin iEta

Parameters

d	Detector
r	Ring
iEta	Eta bin

Returns: max weight or <= 0 in case of problems

Definition at line 773 of file AliFMDDensityCalculator.cxx.

Referenced by GetMaxWeight(), and NParticles().

Int_t AliFMDDensityCalculator::GetMaxWeight	(	UShort_t	d,
		Char_t	r,
		Float_t	eta
	)		const

protected

Find the (cached) maximum weight for FMDdr iat $\eta$

Parameters

d	Detector
r	Ring
eta	Eta bin

Returns: max weight or <= 0 in case of problems

Definition at line 813 of file AliFMDDensityCalculator.cxx.

UShort_t AliFMDDensityCalculator::GetMinQuality ( ) const

inline

Set the minimum quality of the energy loss fits

Returns: Cut value

Definition at line 241 of file AliFMDDensityCalculator.h.

Double_t AliFMDDensityCalculator::GetMultCut	(	UShort_t	d,
		Char_t	r,
		Double_t	eta,
		Bool_t	errors = `true`
	)		const

Get the multiplicity cut. If the user has set fMultCut (via SetMultCut) then that value is used. If not, then the lower value of the fit range for the enery loss fits is returned.

Parameters

d	Detector
r	Ring
eta	Psuedo-rapidity
errors	Factor in errors

Returns: Lower cut on multiplicity

Definition at line 326 of file AliFMDDensityCalculator.cxx.

Referenced by Calculate().

Double_t AliFMDDensityCalculator::GetMultCut	(	UShort_t	d,
		Char_t	r,
		Int_t	ieta,
		Bool_t	errors = `true`
	)		const

Get the multiplicity cut. If the user has set fMultCut (via SetMultCut) then that value is used. If not, then the lower value of the fit range for the enery loss fits is returned.

Parameters

d	Detector
r	Ring
ieta	Psuedo-rapidity bin
errors	Factor in errors

Returns: Lower cut on multiplicity

Definition at line 309 of file AliFMDDensityCalculator.cxx.

AliFMDDensityCalculator::RingHistos * AliFMDDensityCalculator::GetRingHistos	(	UShort_t	d,
		Char_t	r
	)		const

protected

Get the ring histogram container

Parameters

d	Detector
r	Ring

Returns: Ring histogram container

Definition at line 264 of file AliFMDDensityCalculator.cxx.

Referenced by Calculate().

Float_t AliFMDDensityCalculator::NParticles	(	Float_t	mult,
		UShort_t	d,
		Char_t	r,
		Float_t	eta,
		Bool_t	lowFlux
	)		const

protectedvirtual

Get the number of particles corresponding to the signal mult

Parameters

mult	Signal
d	Detector
r	Ring
eta	Pseudo-rapidity
lowFlux	Low-flux flag

Returns: The number of particles

Definition at line 835 of file AliFMDDensityCalculator.cxx.

Referenced by Calculate().

AliFMDDensityCalculator & AliFMDDensityCalculator::operator= ( const AliFMDDensityCalculator & o )

Assignement operator

Parameters

o	Object to assign from

Returns: Reference to this object

Definition at line 193 of file AliFMDDensityCalculator.cxx.

Referenced by AliFMDMCDensityCalculator::operator=().

void AliFMDDensityCalculator::Print ( Option_t * option = "" ) const

Print information

Parameters

option

Print options

max Print max weights

Definition at line 1170 of file AliFMDDensityCalculator.cxx.

Referenced by AliForwardMultiplicityBase::Print(), and AliForwardQATask::Print().

void AliFMDDensityCalculator::SetCuts ( const AliFMDMultCuts & c )

inline

Set the cuts to use

Parameters

c	Cuts to use

Definition at line 260 of file AliFMDDensityCalculator.h.

Referenced by AddTaskForwardQA(), ForwardAODConfig(), and TunedConfig().

void AliFMDDensityCalculator::SetDebug ( Int_t dbg = 1 )

inline

Set the debug level. The higher the value the more output

Parameters

dbg	Debug level

Definition at line 136 of file AliFMDDensityCalculator.h.

Referenced by AliForwardQATask::SetDebug(), and AliForwardMultiplicityBase::SetDebug().

void AliFMDDensityCalculator::SetDoTiming ( Bool_t enable = true )

inline

Definition at line 137 of file AliFMDDensityCalculator.h.

void AliFMDDensityCalculator::SetLumping	(	Int_t	eta,
		Int_t	phi
	)

inline

Set the luming factors used in the Poisson method

Parameters

eta	Must be 1 or larger
phi	Must be 1 or larger

Definition at line 184 of file AliFMDDensityCalculator.h.

Referenced by AddTaskForwardQA(), ForwardAODConfig(), and TunedConfig().

void AliFMDDensityCalculator::SetMaxOutliers ( Double_t ratio = 0.10 )

inline

Set the maximum ratio of outlier bins to the total number of bins with data.

Parameters

ratio Maximum ratio (number between 0 and 1)

Definition at line 200 of file AliFMDDensityCalculator.h.

Referenced by AddTaskForwardQA(), ForwardAODConfig(), and TunedConfig().

void AliFMDDensityCalculator::SetMaxParticles ( UShort_t m )

inline

Maximum particle weight to use

Parameters

m

Definition at line 143 of file AliFMDDensityCalculator.h.

Referenced by AddTaskForwardQA(), ForwardAODConfig(), and TunedConfig().

void AliFMDDensityCalculator::SetMinQuality ( UShort_t cut = 10 )

inline

Set the minimum quality of the energy loss fits

Parameters

cut Cut value

Definition at line 193 of file AliFMDDensityCalculator.h.

Referenced by AddTaskForwardQA(), ForwardAODConfig(), and TunedConfig().

void AliFMDDensityCalculator::SetOutlierCut ( Double_t cut = 0.50 )

inline

Set the maximum relative diviation between $N_{ch}^{Poisson}$ and $N_{ch}^{\Delta}$

Parameters

cut	Relative cut (number between 0 and 1)

Definition at line 207 of file AliFMDDensityCalculator.h.

Referenced by ForwardAODConfig(), and TunedConfig().

void AliFMDDensityCalculator::SetRecalculatePhi ( Bool_t use )

inline

In case of a displaced vertices recalculate eta and angle correction

Parameters

use	recalculate or not

Definition at line 165 of file AliFMDDensityCalculator.h.

Referenced by ForwardAODConfig(), and TunedConfig().

void AliFMDDensityCalculator::SetupForData ( const TAxis & etaAxis )

virtual

Initialize this sub-algorithm

Parameters

etaAxis Not used

Reimplemented in AliFMDMCDensityCalculator.

Definition at line 242 of file AliFMDDensityCalculator.cxx.

Referenced by AliForwardQATask::PreData(), AliForwardMultiplicityBase::PreData(), and AliFMDMCDensityCalculator::SetupForData().

void AliFMDDensityCalculator::SetUsePhiAcceptance ( UShort_t u = kPhiCorrectNch )

inline

Set whether to use the phi acceptance correction.

How the phi acceptance is used depends on the value passed.

0: No phi acceptance
1: Phi acceptance correction done to estimate of particles
2: Phi acceptance correction done to energy deposited

Parameters

u	If >0, use the phi acceptance (default is false)

Definition at line 177 of file AliFMDDensityCalculator.h.

Referenced by AddTaskForwardQA(), ForwardAODConfig(), and TunedConfig().

void AliFMDDensityCalculator::SetUsePoisson ( Bool_t u )

inline

Set whether to use poisson statistics to estimate the number of particles that has hit within a region. If this is true, then the average charge particle density is given by

\[ \lambda = -\log\left(\frac{N_e}{N_t}\right) \]

where $N_e$ is the number of strips within the region that has no hits over threshold, and $N_t$ is the total number of strips in the region/

Parameters