Correction scripts

Collaboration diagram for Correction scripts:

Files
file	CompareCorrs.C
	Utilities for comparing correction objects.
file	DrawAnaELoss.C
	Script to draw the energy loss fits from the output file of AliFMDELossFitter(Task).
file	DrawCorrAcc.C
file	DrawCorrAcc2.C
file	DrawCorrCentralSecMap2.C
file	DrawCorrELoss.C
file	DrawCorrSecMap.C
file	DrawCorrSecMap2.C
file	DrawCorrVtxBias.C

Classes
struct	Canvas

Functions
void	CompareCentralSecMaps (const char *fn1, const char *fn2, const char *n1=0, const char *n2=0, bool load=true)
void	GetObjects (const char *what, const char *fn1, const char *fn2, TObject *&o1, TObject *&o2)
void	GetObjects (UShort_t what, const char *fn1, const char *fn2, TObject *&o1, TObject *&o2)
void	CompareSecMaps (const char *fn1, const char *fn2, const char *n1=0, const char *n2=0, bool load=true)
void	CompareVtxBias (const char *fn1, const char *fn2, const char *n1=0, const char *n2=0, bool load=true)
void	CompELossFits (const char *fname1, const char *fname2)
TList *	OpenFile (const char *fname)
TList *	CheckFitter (const char *fname="AnalysisResults.root")
TCanvas *	CheckCanvas ()
void	CleanStack (THStack *stack)
THStack *	AddToStack (TList *stacks, TList *list, const char *name)
void	DrawSummary (const char *fname="forward_eloss.root", bool onlySummary=true)
void	DrawRings (const char *fname="AnalysisResults.root")
void	DrawEtaBins (const char *fname="AnalysisResults.root")
void	DrawAnaELoss (const char *fname="forward_eloss.root", bool onlySummary=true)
void	DrawCorrAcc (ULong_t runNo, UShort_t sys, UShort_t sNN, const char *fname=0, Bool_t details=true)
void	ClearCanvas (TCanvas *c)
void	DrawCorrAcc2 (const char *fname, const char *option="colz")
void	DrawCorrCentralSecMap2 (const char *fname, const char *option="colz", bool tracklets=true)
void	DrawCorrELoss (ULong_t runNo, UShort_t sys, UShort_t sNN, Short_t field, Bool_t mc=false, Bool_t sat=false, const char *fname=0, Bool_t details=true)
void	DrawCorrSecMap (ULong_t runNo, UShort_t sys, UShort_t sNN, Short_t field, Bool_t sat, const char *fname=0, Bool_t details=true)
void	DrawCorrSecMap2 (const char *fname, const char *option="colz")
void	DrawCorrVtxBias (const char *fname, const char *option="colz")
TH2D *	MakeOneRing (UShort_t d, Char_t r, Double_t vz, Int_t &nDead, std::ostream *deadScript)
void	ExtractAcceptance (Int_t runNo=121526, Int_t nVtxBins=10, Float_t vtxLow=-10, Float_t vtxHigh=10)
void	ExtractELoss (const char *fname="forward_eloss.root", Bool_t mc=false, ULong_t runNo=0, UShort_t sys=0, UShort_t sNN=0, Short_t fld=999, Bool_t sat=false)
void	ExtractMCCorr (const char *fname, ULong_t runNo=0, UShort_t sys=0, UShort_t sNN=0, Short_t fld=999, Bool_t sat=false)
TH2D *	MakeOneRing (UShort_t d, Char_t r, Double_t vz, Int_t &nDead)
void	ExtractAcceptance (Int_t runNo=121526, Int_t system=1, Float_t energy=900, Float_t field=5, Int_t nVtxBins=10, Float_t vtxLow=-10, Float_t vtxHigh=10)
void	MakeCorrRepository ()
Bool_t	MakeSecMap (TList *list, Double_t low, Double_t high, AliFMDCorrSecondaryMap *m)
void	MakeCorrSecMap (const char *filename, const char *sys="pp", UShort_t cms=900, Short_t field=+5)

Detailed Description

Scripts and functions to handle corrections

Function Documentation

THStack* AddToStack	(	TList *	stacks,
		TList *	list,
		const char *	name
	)

Add a stack to list

Parameters

stacks	List to add to
list	Input list
name	Name of stack to add

Returns

The found stack

Definition at line 181 of file DrawAnaELoss.C.

Referenced by DrawSummary().

TCanvas* CheckCanvas

(

)

Make canvas if not done already

Returns

Canvas

Definition at line 112 of file DrawAnaELoss.C.

Referenced by DrawAnaELoss(), DrawEtaBins(), DrawRings(), and DrawSummary().

TList* CheckFitter

(

const char *

fname = "AnalysisResults.root"

)

Open file if not done already

Parameters

fname

File to open

Returns

List of fits

Definition at line 97 of file DrawAnaELoss.C.

Referenced by DrawEtaBins(), DrawRings(), and DrawSummary().

void CleanStack

(

THStack *

stack

)

Clean stack for full ring histograms

Parameters

stack

Input stack

Definition at line 157 of file DrawAnaELoss.C.

Referenced by AddToStack().

void ClearCanvas

(

TCanvas *

c

)

Clear canvas

Parameters

c	Canvas to clear

Definition at line 16 of file DrawCorrAcc2.C.

Referenced by DrawCorrAcc2(), DrawCorrCentralSecMap2(), DrawCorrSecMap2(), DrawCorrVtxBias(), AliFMDCorrAcceptance::ls(), DetailsComparer::MakeTitlePage(), and DetailsComparer::ProcessOne().

void CompareCentralSecMaps	(	const char *	fn1,
		const char *	fn2,
		const char *	n1 = 0,
		const char *	n2 = 0,
		bool	load = true
	)

Compare secondary maps

Parameters

fn1	File 1
fn2	File 2
n1	Name 1
n2	Name 2
load

Definition at line 22 of file CompareCentralSecMaps.C.

void CompareSecMaps	(	const char *	fn1,
		const char *	fn2,
		const char *	n1 = 0,
		const char *	n2 = 0,
		bool	load = true
	)

Parameters

fn1
fn2
n1
n2
load

Definition at line 14 of file CompareSecMaps.C.

void CompareVtxBias	(	const char *	fn1,
		const char *	fn2,
		const char *	n1 = 0,
		const char *	n2 = 0,
		bool	load = true
	)

Parameters

fn1
fn2
n1
n2
load

Definition at line 14 of file CompareVtxBias.C.

void CompELossFits	(	const char *	fname1,
		const char *	fname2
	)

Parameters

fname1
fname2

Definition at line 10 of file CompELossFits.C.

void DrawAnaELoss	(	const char *	fname = "forward_eloss.root",
		bool	onlySummary = true
	)

Draw energy loss fits to a multi-page PDF

The input file is the result of running AliFMDELossFitter - either directly via AliFMDELossFitterTask or as part of a larger train (AliForwardMultiplicityTask or AliForwardMCMultiplicityTask).

*  file
*   +- ForwardResults 
*       +- fmdEnergyFitter 
*           +- chi2   (THStack)
*           +- c      (THStack)
*           +- delta  (THStack)
*           +- xi     (THStack)
*           +- sigma  (THStack)
*           +- sigman (THStack)
*           +- n      (THStack)
*           +- a2     (THStack)
*           +- ...    (THStack)
*           +- an     (THStack)
*           +- FMD1I (TList)
*           |   +- FMD1I_edist (TH1)
*           |   +- EDists      (TList)
*           ...
* 

Parameters

fname	File name to read
onlySummary	Only plot summary pae

Definition at line 489 of file DrawAnaELoss.C.

Referenced by DrawQA().

void DrawCorrAcc	(	ULong_t	runNo,
		UShort_t	sys,
		UShort_t	sNN,
		const char *	fname = 0,
		Bool_t	details = true
	)

Draw energy loss fits to a multi-page PDF.

Input:

The input file is expected to contain a AliFMDCorrELossFit object named elossfits in the top level directory.

Output:

A multi-page PDF. Note, that the PDF generated by ROOT in this way is broken (cannot be read by Acrobat Reader on Windows and MacOSX) and one should pass it through a filter to correct these problems.

Parameters

runNo	Run number
sys	Collision system
sNN	Collision energy in GeV
fname	File name
details	Drawing options

Definition at line 29 of file DrawCorrAcc.C.

void DrawCorrAcc2	(	const char *	fname,
		const char *	option = "colz"
	)

Draw energy loss fits to a multi-page PDF.

Input:

The input file is expected to contain a AliFMDCorrELossFit object named elossfits in the top level directory.

Output:

A multi-page PDF. Note, that the PDF generated by ROOT in this way is broken (cannot be read by Acrobat Reader on Windows and MacOSX) and one should pass it through a filter to correct these problems.

Parameters

fname	File name
option	Drawing options

Definition at line 43 of file DrawCorrAcc2.C.

void DrawCorrCentralSecMap2	(	const char *	fname,
		const char *	option = "colz",
		bool	tracklets = true
	)

Draw energy loss fits to a multi-page PDF.

Input:

The input file is expected to contain a AliFMDCorrELossFit object named elossfits in the top level directory.

Output:

A multi-page PDF. Note, that the PDF generated by ROOT in this way is broken (cannot be read by Acrobat Reader on Windows and MacOSX) and one should pass it through a filter to correct these problems.

Parameters

fname	File name
option	Drawing options
tracklets	Tracklets flag

Definition at line 56 of file DrawCorrCentralSecMap2.C.

void DrawCorrELoss	(	ULong_t	runNo,
		UShort_t	sys,
		UShort_t	sNN,
		Short_t	field,
		Bool_t	mc = false,
		Bool_t	sat = false,
		const char *	fname = 0,
		Bool_t	details = true
	)

Draw energy loss fits to a multi-page PDF.

Input:

The input file is expected to contain a AliFMDCorrELossFit object named elossfits in the top level directory.

Output:

A multi-page PDF. Note, that the PDF generated by ROOT in this way is broken (cannot be read by Acrobat Reader on Windows and MacOSX) and one should pass it through a filter to correct these problems.

Parameters

runNo	Run number
sys	Collision system
sNN	Collision energy in GeV
field	L3 field strength
mc	For simulation input
sat	For satellite interactions
fname	File name
details	Drawing options

Definition at line 46 of file DrawCorrELoss.C.

void DrawCorrSecMap	(	ULong_t	runNo,
		UShort_t	sys,
		UShort_t	sNN,
		Short_t	field,
		Bool_t	sat,
		const char *	fname = 0,
		Bool_t	details = true
	)

Draw secondary maps fits to a multi-page PDF.

Input:

The input is expected to be parameters to obtain a AliFMDCorrSecondaryMap object from the correction manager

Output:

A multi-page PDF. Note, that the PDF generated by ROOT in this way is broken (cannot be read by Acrobat Reader on Windows and MacOSX) and one should pass it through a filter to correct these problems.

Parameters

runNo	Run number
sys	Collision system
sNN	Collision energy in GeV
field	L3 field strength
sat	If true, assume satellite-main collisions
fname	File name
details	Drawing options

Definition at line 31 of file DrawCorrSecMap.C.

void DrawCorrSecMap2	(	const char *	fname,
		const char *	option = "colz"
	)

Draw energy loss fits to a multi-page PDF.

Input:

The input file is expected to contain a AliFMDCorrELossFit object named elossfits in the top level directory.

Output:

A multi-page PDF. Note, that the PDF generated by ROOT in this way is broken (cannot be read by Acrobat Reader on Windows and MacOSX) and one should pass it through a filter to correct these problems.

Parameters

fname	File name
option	Drawing options

Definition at line 43 of file DrawCorrSecMap2.C.

void DrawCorrVtxBias	(	const char *	fname,
		const char *	option = "colz"
	)

Draw energy loss fits to a multi-page PDF.

Input:

The input file is expected to contain a AliFMDCorrELossFit object named elossfits in the top level directory.

Output:

A multi-page PDF. Note, that the PDF generated by ROOT in this way is broken (cannot be read by Acrobat Reader on Windows and MacOSX) and one should pass it through a filter to correct these problems.

Parameters

fname	File name
option	Drawing options

Definition at line 43 of file DrawCorrVtxBias.C.

void DrawEtaBins

(

const char *

fname = "AnalysisResults.root"

)

Draw fits in eta bins

Parameters

fname

Definition at line 389 of file DrawAnaELoss.C.

Referenced by DrawAnaELoss().

void DrawRings

(

const char *

fname = "AnalysisResults.root"

)

Draw ring fits

Parameters

fname

Definition at line 325 of file DrawAnaELoss.C.

Referenced by DrawAnaELoss().

void DrawSummary	(	const char *	fname = "forward_eloss.root",
		bool	onlySummary = true
	)

Draw summary

Parameters

fname	File name to read
onlySummary	Only plot summary pae

Definition at line 207 of file DrawAnaELoss.C.

Referenced by DrawAnaELoss().

void ExtractAcceptance	(	Int_t	runNo = 121526,
		Int_t	system = 1,
		Float_t	energy = 900,
		Float_t	field = 5,
		Int_t	nVtxBins = 10,
		Float_t	vtxLow = -10,
		Float_t	vtxHigh = 10
	)

Parameters

runNo
system
energy
field
nVtxBins
vtxLow
vtxHigh

Definition at line 101 of file MakeAcceptanceCorrection.C.

void ExtractAcceptance	(	Int_t	runNo = 121526,
		Int_t	nVtxBins = 10,
		Float_t	vtxLow = -10,
		Float_t	vtxHigh = 10
	)

Parameters

runNo	Run number
nVtxBins	Number of \(IP_{z}\) bins
vtxLow	Least \(IP_{z}\)
vtxHigh	Largest \(IP_{z}\)

Definition at line 166 of file ExtractAcceptance.C.

void ExtractELoss	(	const char *	fname = "forward_eloss.root",
		Bool_t	mc = false,
		ULong_t	runNo = 0,
		UShort_t	sys = 0,
		UShort_t	sNN = 0,
		Short_t	fld = 999,
		Bool_t	sat = false
	)

Script to draw the energy loss fits Extract the energy loss correction object from file and rename it according to the settings

Parameters

fname	File to extract from
mc	Whether this is from MC data or not
runNo	Run number
sys	Collision system (pp, PbPb)
sNN	Center of mass energy (in GeV) per nucleon
fld	L3 magnetic field (-5,0,5) in kGaus
sat	For satellite interactions

Definition at line 23 of file ExtractELoss.C.

void ExtractMCCorr	(	const char *	fname,
		ULong_t	runNo = 0,
		UShort_t	sys = 0,
		UShort_t	sNN = 0,
		Short_t	fld = 999,
		Bool_t	sat = false
	)

Script to draw the energy loss fits Extract the energy loss correction object from file and rename it according to the settings

Parameters

fname	File to extract from
runNo	Run number
sys	Collision system (pp, PbPb)
sNN	Center of mass energy (in GeV) per nucleon
fld	L3 magnetic field (-5,0,5) in kGaus
sat	For satellite interactions

Definition at line 22 of file ExtractMCCorr.C.

void GetObjects	(	const char *	what,
		const char *	fn1,
		const char *	fn2,
		TObject *&	o1,
		TObject *&	o2
	)

Parameters

what
fn1
fn2
o1
o2

Definition at line 174 of file CompareCorrs.C.

Referenced by CompareSecMaps(), CompareVtxBias(), and GetObjects().

void GetObjects	(	UShort_t	what,
		const char *	fn1,
		const char *	fn2,
		TObject *&	o1,
		TObject *&	o2
	)

Parameters

what
fn1
fn2
o1
o2

Definition at line 229 of file CompareCorrs.C.

void MakeCorrRepository

(

)

Make the repository for corrections

Definition at line 8 of file MakeCorrRepository.C.

void MakeCorrSecMap	(	const char *	filename,
		const char *	sys = "pp",
		UShort_t	cms = 900,
		Short_t	field = +5
	)

Parameters

filename
sys
cms
field

Definition at line 86 of file MakeCorrSecMap.C.

TH2D* MakeOneRing	(	UShort_t	d,
		Char_t	r,
		Double_t	vz,
		Int_t &	nDead
	)

Parameters

d
r
vz
nDead

Returns

Definition at line 14 of file MakeAcceptanceCorrection.C.

Referenced by ExtractAcceptance().

TH2D* MakeOneRing	(	UShort_t	d,
		Char_t	r,
		Double_t	vz,
		Int_t &	nDead,
		std::ostream *	deadScript
	)

Parameters

d	Detector
r	Ring
vz	Z–coordinate of interaction point
nDead	On returm the number of dead strips
deadScript	Output stream for dead strips

Returns

Definition at line 15 of file ExtractAcceptance.C.

Referenced by ExtractAcceptance().

Bool_t MakeSecMap	(	TList *	list,
		Double_t	low,
		Double_t	high,
		AliFMDCorrSecondaryMap *	m
	)

Parameters

list
low
high
m

Returns

Definition at line 15 of file MakeCorrSecMap.C.

Referenced by MakeCorrSecMap().

TList* OpenFile

(

const char *

fname

)

Open a file. The file is expected to contain the directory structure

*  file
*   +- ForwardResults
*       +- fmdEnergyFitter 
*           +- chi2   (THStack)
*           +- c      (THStack)
*           +- delta  (THStack)
*           +- xi     (THStack)
*           +- sigma  (THStack)
*           +- sigman (THStack)
*           +- n      (THStack)
*           +- a2     (THStack)
*           +- ...    (THStack)
*           +- an     (THStack)
*           +- FMD1I (TList)
*           |   +- FMD1I_edist (TH1)
*           |   +- EDists      (TList)
*           ...
* 

Parameters

fname

File to open

Returns

Pointer to the list of objects

Definition at line 65 of file DrawAnaELoss.C.

Referenced by CheckFitter(), AliTrackComparisonESD::CreateOutputObjects(), AliEmcalPhysicsSelectionTask::Terminate(), AliEsdSkimTask::UserCreateOutputObjects(), AliAnalysisTaskEmcalTriggerTreeWriter::UserCreateOutputObjects(), AliAnalysisTaskSOH::UserCreateOutputObjects(), EMCalTriggerPtAnalysis::AliAnalysisTaskPtEfficiencyJets::UserCreateOutputObjects(), AliEmcalMCTreeWriter::UserCreateOutputObjects(), AliAnalysisTaskCmeEse::UserCreateOutputObjects(), AliAnalysisTaskPiKpK0Lamba::UserCreateOutputObjects(), EMCalTriggerPtAnalysis::AliAnalysisTaskTracksInJet::UserCreateOutputObjects(), AliAnalysisTaskJetsEvshape::UserCreateOutputObjects(), AliJetEmbeddingTask::UserCreateOutputObjects(), AliAnalysisTaskHJetEmbed::UserCreateOutputObjects(), AliAnalysisTaskSECompareHF::UserCreateOutputObjects(), AliAnalysisTaskPSHFE::UserCreateOutputObjects(), AliAnalysisTaskHJetDphi::UserCreateOutputObjects(), AliAnalysisTaskCLQA::UserCreateOutputObjects(), AliAnalysisTaskIDFragmentationFunction::UserCreateOutputObjects(), AliAnalysisTaskSEDStarJets::UserCreateOutputObjects(), AliAnalysisTaskCaloTrackCorrelation::UserCreateOutputObjects(), AliAnalysisTaskFilteredTree::UserCreateOutputObjects(), AliAnalysisTaskPID::UserCreateOutputObjects(), AliAnalysisTaskCaloTrackCorrelationM::UserCreateOutputObjects(), AliAnalysisTaskFullppJet::UserCreateOutputObjects(), AliAnalysisTaskPi0v2::UserCreateOutputObjects(), AliCFTaskForDStarAnalysis::UserCreateOutputObjects(), HighPtTracks::AliHighPtReconstructionEfficiency::UserCreateOutputObjects(), AliJetModelBaseTask::UserCreateOutputObjects(), AliCFTaskVertexingHFCutVarFDSub::UserCreateOutputObjects(), AliAnalysisTaskEmcalEmbeddingHelper::UserCreateOutputObjects(), AliAnalysisTaskSEDs::UserCreateOutputObjects(), AliAnalysisTaskEmcalJetHadEPpid::UserCreateOutputObjects(), AliCFTaskVertexingHF::UserCreateOutputObjects(), AliAnalysisTaskSEDplus::UserCreateOutputObjects(), AliEmcalCorrectionTask::UserCreateOutputObjects(), AliTrackletTaskMulti::UserCreateOutputObjects(), AliAnalysisTaskEmcalLight::UserCreateOutputObjects(), EmcalTriggerJets::AliAnalysisTaskEmcalJetSubstructureTree::UserCreateOutputObjects(), and AliAnalysisTaskEmcal::UserCreateOutputObjects().