Public Types
enum	{ kColorMeasured = kOrange-2, kColorTruth = kBlue-3, kColorAccepted = kMagenta-3, kColorTrgVtx = kBlack, kColorUnfolded = kOrange+2, kColorCorrected = kRed+2, kColorError = kBlue-10, kColorALICE = kPink+1, kColorCMS = kGreen+2 }

Public Member Functions
	Unfolder ()

void	Run (const TString &measuredFile, const TString &corrFile, const TString &method="Bayes", Double_t regParam=4)

void	SaveInformation (TDirectory *dir, const TString &method, Int_t mId, Double_t regParam, UShort_t sys, UShort_t sNN, UInt_t trigger, Double_t minIpZ, Double_t maxIpZ, Bool_t self) const

void	SaveSummarize ()

void	ProcessType (TCollection measured, TCollection corrections, UInt_t method, Double_t regParam, TDirectory *out, UShort_t sys, UShort_t sNN)

THStack *	ProcessBin (TCollection measured, TCollection corrections, UInt_t method, Double_t regParam, TDirectory *out)

void	Bin2Stack (const THStack bin, Int_t i, THStack measured, THStack truth, THStack accepted, THStack unfolded, THStack corrected, TH1 *&result)

void	Other2Stack (const TString &name, Int_t i, UShort_t sNN, TMultiGraph allALICE, TMultiGraph allCMS, TGraph &alice, TGraph &cms)

void	Ratio2Stack (Int_t ib, TH1 res, TGraph alice, TGraph cms, THStack all)

TLegend *	StackLegend (THStack *stack)

TGraphAsymmErrors *	GetOther (UShort_t type, Double_t eta, UShort_t sNN, Int_t factor)

Static Public Member Functions
static TCollection *	GetTop (const TString &fileName, Bool_t results=false)

static TObject *	GetObject (TCollection c, const TString &name, TClass cl, Bool_t verbose=true)

static TCollection *	GetCollection (TCollection *c, const TString &name, Bool_t verbose=-true)

static TH1 *	GetH1 (TCollection *c, const TString &name, Bool_t verbose=true)

static TH2 *	GetH2 (TCollection *c, const TString &name, Bool_t verbose=true)

static void	GetParameter (TCollection *c, const TString &name, UShort_t &v)

static void	GetParameter (TCollection *c, const TString &name, ULong_t &v)

static void	GetParameter (TCollection *c, const TString &name, Double_t &v)

static UInt_t	MethodId (TString &method)

static void	AppendAnd (TString &trg, const TString &what)

static void	BinAttributes (Int_t i, Int_t &open, Int_t &closed, Float_t &size, Double_t &factor)

Detailed Description

Class to do unfolding of raw histograms produced by AliForwardMultDists

Definition at line 33 of file UnfoldMultDists.C.

Member Enumeration Documentation

anonymous enum

Colours used

Enumerator
kColorMeasured
kColorTruth
kColorAccepted
kColorTrgVtx
kColorUnfolded
kColorCorrected
kColorError
kColorALICE
kColorCMS

Definition at line 39 of file UnfoldMultDists.C.

Constructor & Destructor Documentation

Unfolder::Unfolder ( )

inline

Constructor

Definition at line 54 of file UnfoldMultDists.C.

Member Function Documentation

static void Unfolder::AppendAnd	(	TString &	trg,
		const TString &	what
	)

inlinestatic

Append an & to a string and the next term.

Parameters

trg	Output string
what	Term

Definition at line 362 of file UnfoldMultDists.C.

Referenced by SaveInformation().

void Unfolder::Bin2Stack	(	const THStack *	bin,
		Int_t	i,
		THStack *	measured,
		THStack *	truth,
		THStack *	accepted,
		THStack *	unfolded,
		THStack *	corrected,
		TH1 *&	result
	)

inline

Add the bin histograms to our summary stacks

Parameters

bin	Bin stack
i	Current off-set in the stacks
measured	All measured \( P(N_{ch})\)
truth	All MC truth \( P(N_{ch})\)
accepted	All MC accepted \( P(N_{ch})\)
unfolded	All unfolded \( P(N_{ch})\)
corrected	All corrected \( P(N_{ch})\)
result	The result in this bin

Definition at line 790 of file UnfoldMultDists.C.

Referenced by ProcessType().

static void Unfolder::BinAttributes	(	Int_t	i,
		Int_t &	open,
		Int_t &	closed,
		Float_t &	size,
		Double_t &	factor
	)

inlinestatic

Definition at line 760 of file UnfoldMultDists.C.

Referenced by Bin2Stack(), and Other2Stack().

static TCollection* Unfolder::GetCollection	(	TCollection *	c,
		const TString &	name,
		Bool_t	verbose = `-true`
	)

inlinestatic

Get a collection

Parameters

c	Parent collection
name	Name of object to findf
verbose	Be verbose

Returns

Definition at line 115 of file UnfoldMultDists.C.

Referenced by ProcessType(), and Run().

static TH1* Unfolder::GetH1	(	TCollection *	c,
		const TString &	name,
		Bool_t	verbose = `true`
	)

inlinestatic

Get a 1D histogram from a collection

Parameters

c	Collection
name	Nanme of histogram
verbose	Be verbose

Returns: Pointer to object or null

Definition at line 132 of file UnfoldMultDists.C.

Referenced by ProcessBin().

static TH2* Unfolder::GetH2	(	TCollection *	c,
		const TString &	name,
		Bool_t	verbose = `true`
	)

inlinestatic

Get a 2D histogram from a collection

Parameters

c	Collection
name	Nanme of histogram
verbose	Be verbose

Returns: Pointer to object or null

Definition at line 145 of file UnfoldMultDists.C.

Referenced by ProcessBin().

static TObject* Unfolder::GetObject	(	TCollection *	c,
		const TString &	name,
		TClass *	cl,
		Bool_t	verbose = `true`
	)

inlinestatic

Get an object from a collection

Parameters

c	Collection
name	Name of object
cl	Possible class to check against
verbose	Be verbose

Returns: Pointer to object or null

Definition at line 89 of file UnfoldMultDists.C.

Referenced by GetCollection(), GetH1(), GetH2(), and GetParameter().

TGraphAsymmErrors* Unfolder::GetOther	(	UShort_t	type,
		Double_t	eta,
		UShort_t	sNN,
		Int_t	factor
	)

inline

Definition at line 1021 of file UnfoldMultDists.C.

Referenced by Other2Stack().

static void Unfolder::GetParameter	(	TCollection *	c,
		const TString &	name,
		UShort_t &	v
	)

inlinestatic

Get an unsigned short parameter from the collection

Parameters

c	Collection
name	Parameter name
v	Value

Returns: Value

Definition at line 158 of file UnfoldMultDists.C.

Referenced by Run().

static void Unfolder::GetParameter	(	TCollection *	c,
		const TString &	name,
		ULong_t &	v
	)

inlinestatic

Get an unsigned short parameter from the collection

Parameters

c	Collection
name	Parameter name
v	Value

Returns: Value

Definition at line 172 of file UnfoldMultDists.C.

static void Unfolder::GetParameter	(	TCollection *	c,
		const TString &	name,
		Double_t &	v
	)

inlinestatic

Get an unsigned short parameter from the collection

Parameters

c	Collection
name	Parameter name
v	Value

Returns: Value

Definition at line 186 of file UnfoldMultDists.C.

static TCollection* Unfolder::GetTop	(	const TString &	fileName,
		Bool_t	results = `false`
	)

inlinestatic

Get a top collection from a file

Parameters

fileName	Name of file
results	Wheter it's the results collection or not

Returns: Collection or null

Definition at line 63 of file UnfoldMultDists.C.

Referenced by Run().

static UInt_t Unfolder::MethodId ( TString & method )

inlinestatic

Get the method identifier

Parameters

method Method

Returns: Method identifier

Definition at line 198 of file UnfoldMultDists.C.

Referenced by Run().

void Unfolder::Other2Stack	(	const TString &	name,
		Int_t	i,
		UShort_t	sNN,
		TMultiGraph *	allALICE,
		TMultiGraph *	allCMS,
		TGraph *&	alice,
		TGraph *&	cms
	)

inline

Add distributions from other experiments to stacks

Parameters

name	Name of current bin
i	Index of current bin
sNN	Center of mass energy
allALICE	Stack of ALICE data
allCMS	Stack of CMS data
alice	Possible ALICE result on return
cms	Possible CMS result on return

Definition at line 865 of file UnfoldMultDists.C.

Referenced by ProcessType().

THStack* Unfolder::ProcessBin	(	TCollection *	measured,
		TCollection *	corrections,
		UInt_t	method,
		Double_t	regParam,
		TDirectory *	out
	)

inline

Process a single eta bin

Parameters

measured	Input collection of measured data
corrections	Input collection of correction data
method	Unfolding method to use
regParam	Regularisation parameter
out	Output directory.

Returns: Stack of histograms or null

Definition at line 604 of file UnfoldMultDists.C.

Referenced by ProcessType().

void Unfolder::ProcessType	(	TCollection *	measured,
		TCollection *	corrections,
		UInt_t	method,
		Double_t	regParam,
		TDirectory *	out,
		UShort_t	sys,
		UShort_t	sNN
	)

inline

Process a single type - i.e., one of symmetric, positive, negative, or other - by looping over all contained objects and call ProcessBin for each found bin.

Parameters

measured	Input collection of measured data
corrections	Input collection of correction data
method	Unfolding method to use
regParam	Regularisation parameter
out	Output directory.
sys	Collision system
sNN	Collision energy

Definition at line 503 of file UnfoldMultDists.C.

Referenced by Run().

void Unfolder::Ratio2Stack	(	Int_t	ib,
		TH1 *	res,
		TGraph *	alice,
		TGraph *	cms,
		THStack *	all
	)

inline

Create ratios to other data

Parameters

ib	Bin number
res	Result
alice	ALICE result if any
cms	CMS result if any
all	Stack to add ratio to

Definition at line 929 of file UnfoldMultDists.C.

Referenced by ProcessType().

void Unfolder::Run	(	const TString &	measuredFile,
		const TString &	corrFile,
		const TString &	method = `"Bayes"`,
		Double_t	regParam = `4`
	)

inline

Run the unfolding and correction task.

The measuredFile is assumed to have the structure

*     /- ForwardMultSums         (TCollection)
*          |- [type]             (TCollection)
*          |    |- [bin]         (TCollection) 
*          |    |    `- rawDist  (TH1)
*          |    |- [bin]
*          |    ...
*          |- [type]
*          ...
*

and corrFile is assumed to have the structure

*     /- ForwardMultResults            (TCollection)
*          |- [type]                   (TCollection)
*          |    |- [bin]               (TCollection) 
*          |    |    |- truth          (TH1)
*          |    |    |- truthAccepted  (TH1)
*          |    |    |- triggerVertex  (TH1)
*          |    |    `- response       (TH2)
*          |    |- [bin]
*          |    ...
*          |- [type]
*          ...
*

where [type] is one of symmetric, positive, negative, or other, and [bin] is the \( \eta\) bin named like

*   [bin]          := [eta_spec] _ [eta_spec]
*   [eta_spec]     := [sign_char] [integer_part] d [decimal_part]
*   [sign_part]    := p          positive eta 
*                  |  m          negative eta 
*   [integer_part] := [number]
*   [decimal_part] := [number] [number]
*   [number]       := 0 | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 
*

That is, the bin \( -3\le\eta\ge3\) is labeled m3d00_p3d00, \( 0\le\eta\ge2.5\) is p0d00_p2d50

measuredFile and corrFile can point to the same file. If corrFile is not specified, it is assumed that measuredFile has the expected corrFile in addition to the expected content of that file.

Parameters

measuredFile	Name of file containing measured data
corrFile	Name of file containing correction data
method	Unfolding method to use
regParam	Regularization parameter

Definition at line 285 of file UnfoldMultDists.C.

Referenced by UnfoldMultDists().

void Unfolder::SaveInformation	(	TDirectory *	dir,
		const TString &	method,
		Int_t	mId,
		Double_t	regParam,
		UShort_t	sys,
		UShort_t	sNN,
		UInt_t	trigger,
		Double_t	minIpZ,
		Double_t	maxIpZ,
		Bool_t	self
	)		const

inline

Store some information on the output

Parameters

dir	Where to store
method	Method used
mId	Method identifier
regParam	Regularization parameter
sys	Collision system
sNN	Center of mass energy
trigger	Trigger mask
minIpZ	Least z coordinate of interaction point
maxIpZ	Largest z coordinate of interaction point
self	Self-consistency check

Bits of the trigger pattern

In-elastic collision

In-elastic collision with at least one SPD tracklet

Non-single diffractive collision

Empty bunch crossing

A-side trigger

B(arrel) trigger

C-side trigger

Empty trigger

pileup from SPD

true NSD from MC

Offline MB triggered

At least one SPD cluster

V0-AND trigger

Satellite event

Definition at line 381 of file UnfoldMultDists.C.

Referenced by Run().

void Unfolder::SaveSummarize ( )

inline

Save a script to do a summary of this step

Definition at line 477 of file UnfoldMultDists.C.

Referenced by Run().

TLegend* Unfolder::StackLegend ( THStack * stack )

inline

Get or create a stack legend. This is done by adding a TLegend object to the list of functions for the first histogram in the stack.

Parameters

stack Stack to get the legend from/modify

Returns: The legend object or null

Definition at line 993 of file UnfoldMultDists.C.

Referenced by Bin2Stack(), ProcessBin(), and Ratio2Stack().

The documentation for this struct was generated from the following file:

PWGLF/FORWARD/analysis2/scripts/UnfoldMultDists.C

Public Types

Public Member Functions

Static Public Member Functions

Detailed Description

Member Enumeration Documentation

Constructor & Destructor Documentation

Member Function Documentation