#include <AliFMDHistCollector.h>

Inheritance diagram for AliFMDHistCollector:

Classes
struct	VtxBin

Public Types
enum	MergeMethod { kStraightMean, kStraightMeanNoZero, kWeightedMean, kLeastError, kSum, kPreferInner, kPreferOuter }

enum	FiducialMethod { kByCut, kDistance }

enum	FMDRingBits { kFMD1I =0x01, kFMD1 =kFMD1I, kFMD2I =0x02, kFMD2O =0x04, kFMD2 =kFMD2I\|kFMD2O, kFMD3I =0x08, kFMD3O =0x10, kFMD3 =kFMD3I\|kFMD3O }

Public Member Functions
	AliFMDHistCollector ()

	AliFMDHistCollector (const char *title)

	AliFMDHistCollector (const AliFMDHistCollector &o)

virtual	~AliFMDHistCollector ()

AliFMDHistCollector &	operator= (const AliFMDHistCollector &)

virtual void	SetupForData (const TAxis &vtxAxis, const TAxis &etaAxis)

virtual Bool_t	Collect (const AliForwardUtil::Histos &hists, AliForwardUtil::Histos &sums, UShort_t vtxBin, TH2D &out, Double_t cent=-1.0, Bool_t eta2phi=false, Bool_t add=true)

virtual void	CreateOutputObjects (TList *dir)

void	SetMergeMethod (MergeMethod m)

MergeMethod	GetMergeMethod () const

void	SetFiducialMethod (FiducialMethod m)

void	SetNCutBins (UInt_t n=2)

void	SetCorrectionCut (Float_t cut=0.5)

void	SetFMDRingsToSkip (UShort_t mask)

void	SetMakeBGHitMaps (Bool_t use)

void	SetMakeCentralitySums (Bool_t use)

void	SetDebug (Int_t dbg=1)

void	Print (Option_t *option="") const

Protected Member Functions
VtxBin *	GetVtxBin (Int_t ivtx)

const VtxBin *	GetVtxBin (Int_t ivtx) const

Static Protected Member Functions
static void	GetDetRing (Int_t idx, UShort_t &d, Char_t &r)

static Int_t	GetIdx (UShort_t d, Char_t r)

static Bool_t	CheckSkip (UShort_t d, Char_t r, UShort_t skips)

static Bool_t	CheckCorrection (FiducialMethod m, Double_t cut, const TH2D *bg, Int_t ie, Int_t ip)

static void	MergeBins (MergeMethod m, Double_t c, Double_t e, Double_t oc, Double_t oe, Double_t &rc, Double_t &re)

Protected Attributes
Int_t	fNCutBins

Float_t	fCorrectionCut

Int_t	fDebug

TList *	fList

TH2D *	fSumRings

TH2D *	fCoverage

TH1D *	fSkipped

MergeMethod	fMergeMethod

FiducialMethod	fFiducialMethod

UShort_t	fSkipFMDRings

Bool_t	fBgAndHitMaps

TObjArray *	fVtxList

TList *	fByCent
	Per-vertex list. More...

Bool_t	fDoByCent

Detailed Description

This class collects the event histograms into single histograms, one for each ring in each vertex bin.

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

HistCollector used:

AliFMDCorrSecondaryMap

Definition at line 43 of file AliFMDHistCollector.h.

Member Enumeration Documentation

enum AliFMDHistCollector::FiducialMethod

How to obtain the fiducial cuts

Enumerator
kByCut	Select bins by fixed cut. Bins with a secondary correction less than the cut is considered as non-valid
kDistance	A bin is considered non-valid, if it is less then twice as large as it's neighbors (in eta)

Definition at line 102 of file AliFMDHistCollector.h.

enum AliFMDHistCollector::FMDRingBits

FMD ring bits for skipping

Enumerator
kFMD1I
kFMD1
kFMD2I
kFMD2O
kFMD2
kFMD3I
kFMD3O
kFMD3

Definition at line 117 of file AliFMDHistCollector.h.

enum AliFMDHistCollector::MergeMethod

Methods to use when merging overlapping bins $b_1$, $b_2$ with content $c_1$, $c_2$, and errors $e_1$, $e_2$ into bin $ b$ with content $c$ and error $e$

Enumerator
kStraightMean	\[ c = \frac{1}{2}(c_1+c_2) \] \[ e = \sqrt{e_1^2+e_2^2} \]
kStraightMeanNoZero	As above, exept zero's are ignored
kWeightedMean	\[ c = \frac{\frac{c_1}{e_1^2}+\frac{c_2}{e_2^2}}{ \frac{1}{e_1^2}+\frac{1}{e_2^2}} \] \[ e = \sqrt{\frac{1}{\frac{1}{e_1^2}+\frac{1}{e_2^2}}} \]
kLeastError	\[ c = \left\{\begin{array}{cl} c_1 & \mbox{if $e_1 < e_2$} \\ c_2 & \mbox{otherwise}\end{array}\right. \]
kSum	Just sum the signals
kPreferInner	In overlaps, prefer inners, or if both are inners, do the straight mean
kPreferOuter	In overlaps, prefer outers, or if both are outers (doesn't happen), do the straight mean

Definition at line 51 of file AliFMDHistCollector.h.

Constructor & Destructor Documentation

AliFMDHistCollector::AliFMDHistCollector ( )

Constructor

Definition at line 40 of file AliFMDHistCollector.cxx.

AliFMDHistCollector::AliFMDHistCollector ( const char * title )

Constructor

Parameters

title Name of object

Definition at line 60 of file AliFMDHistCollector.cxx.

AliFMDHistCollector::AliFMDHistCollector ( const AliFMDHistCollector & o )

Copy constructor

Parameters

o	Object to copy from

Definition at line 80 of file AliFMDHistCollector.cxx.

AliFMDHistCollector::~AliFMDHistCollector ( )

virtual

Destructor

Definition at line 101 of file AliFMDHistCollector.cxx.

Member Function Documentation

Bool_t AliFMDHistCollector::CheckCorrection	(	FiducialMethod	m,
		Double_t	cut,
		const TH2D *	bg,
		Int_t	ie,
		Int_t	ip
	)

staticprotected

Check the correction

Parameters

m	Fiducial method used
cut	Cut value
bg	Secondary map
ie	$\eta$ bin
ip	$\varphi$ bin

Returns: true if OK.

Definition at line 244 of file AliFMDHistCollector.cxx.

Referenced by SetDebug(), and AliFMDHistCollector::VtxBin::SetupForData().

Bool_t AliFMDHistCollector::CheckSkip	(	UShort_t	d,
		Char_t	r,
		UShort_t	skips
	)

staticprotected

Check if the detector d, ring r is listed in the skips bit mask. If the detector/ring is in the mask, return true.

That is, use case is

for (UShort_t d=1. d<=3, d++) {
  UShort_t nr = (d == 1 ? 1 : 2);
  for (UShort_t q = 0; q < nr; q++) { 
    Char_t r = (q == 0 ? 'I' : 'O');
    if (CheckSkips(d, r, skips)) continue; 
    // Process detector/ring 
  }
}

Parameters

d	Detector
r	Ring
skips	Mask of detector/rings to skip

Returns: True if detector d, ring r is in the mask skips

Definition at line 295 of file AliFMDHistCollector.cxx.

Referenced by AliFMDHistCollector::VtxBin::Collect(), SetDebug(), and AliFMDHistCollector::VtxBin::SetupForData().

Bool_t AliFMDHistCollector::Collect	(	const AliForwardUtil::Histos &	hists,
		AliForwardUtil::Histos &	sums,
		UShort_t	vtxBin,
		TH2D &	out,
		Double_t	cent = `-1.0`,
		Bool_t	eta2phi = `false`,
		Bool_t	add = `true`
	)

virtual

Do the calculations

Parameters

hists	Cache of histograms
sums	Cache to sum ring histograms in
vtxBin	Vertex bin (1 based)
out	Output histogram
cent	Centrality
eta2phi	Copy eta coverage to phi acceptance
add	If true, add to internal caches

Returns: true on successs

Definition at line 467 of file AliFMDHistCollector.cxx.

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

void AliFMDHistCollector::CreateOutputObjects ( TList * dir )

virtual

Output diagnostic histograms to directory

Parameters

dir	List to write in

Definition at line 277 of file AliFMDHistCollector.cxx.

Referenced by AliForwardMultiplicityBase::Book().

void AliFMDHistCollector::GetDetRing	(	Int_t	idx,
		UShort_t &	d,
		Char_t &	r
	)

staticprotected

Get the detector and ring from the ring index

Parameters

idx	Ring index
d	On return, the detector or 0 in case of errors
r	On return, the ring id or '0' in case of errors

Definition at line 328 of file AliFMDHistCollector.cxx.

Referenced by Print(), SetDebug(), SetupForData(), and AliFMDHistCollector::VtxBin::SetupForData().

Int_t AliFMDHistCollector::GetIdx	(	UShort_t	d,
		Char_t	r
	)

staticprotected

Get the ring index from detector number and ring identifier

Parameters

d	Detector
r	Ring identifier

Returns: ring index or -1 in case of problems

Definition at line 306 of file AliFMDHistCollector.cxx.

Referenced by AliFMDHistCollector::VtxBin::GetFirst(), AliFMDHistCollector::VtxBin::GetFirstAndLast(), AliFMDHistCollector::VtxBin::GetLast(), AliFMDHistCollector::VtxBin::GetOverlap(), and SetDebug().

MergeMethod AliFMDHistCollector::GetMergeMethod ( ) const

inline

Definition at line 196 of file AliFMDHistCollector.h.

AliFMDHistCollector::VtxBin * AliFMDHistCollector::GetVtxBin ( Int_t ivtx )

protected

Get a vertex bin

Parameters

ivtx	Bin number (1-nVz)

Returns: Bin or null

Definition at line 351 of file AliFMDHistCollector.cxx.

Referenced by Collect(), and Print().

const AliFMDHistCollector::VtxBin * AliFMDHistCollector::GetVtxBin ( Int_t ivtx ) const

protected

Get a vertex bin

Parameters

ivtx	Bin number (1-nVz)

Returns: Bin or null

Definition at line 362 of file AliFMDHistCollector.cxx.

void AliFMDHistCollector::MergeBins	(	MergeMethod	m,
		Double_t	c,
		Double_t	e,
		Double_t	oc,
		Double_t	oe,
		Double_t &	rc,
		Double_t &	re
	)

staticprotected

Merge bins accoring to set method

Parameters

m	Merging method
c	Current content
e	Current error
oc	Old content
oe	Old error
rc	On return, the new content
re	On return, tne new error

Definition at line 374 of file AliFMDHistCollector.cxx.

Referenced by AliFMDHistCollector::VtxBin::Collect(), and SetDebug().

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

Assignement operator

Parameters

o	Object to assign from

Returns: Reference to this object

Definition at line 108 of file AliFMDHistCollector.cxx.

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

Print information

Parameters

option Not used

Definition at line 514 of file AliFMDHistCollector.cxx.

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

void AliFMDHistCollector::SetCorrectionCut ( Float_t cut = 0.5 )

inline

Set the correction cut, that is, when bins in the secondary correction maps have a value lower than this cut, they are considered uncertain and not used

Parameters

cut Cut-off

Definition at line 217 of file AliFMDHistCollector.h.

Referenced by ForwardAODConfig(), and TunedConfig().

void AliFMDHistCollector::SetDebug ( Int_t dbg = 1 )

inline

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

Parameters