Plot N productions comparison analysis of QA histograms from EMCal PWG-GA wagon. More...

#include <TFile.h>
#include <TDirectoryFile.h>
#include <TList.h>
#include <TString.h>
#include <TROOT.h>
#include <TStyle.h>
#include <TH1F.h>
#include <TH2F.h>
#include <TCanvas.h>
#include <TPad.h>
#include <TLegend.h>
#include <TObject.h>
#include <TAxis.h>
#include <TGaxis.h>
#include <TLine.h>
#include <TF1.h>
#include <TMath.h>

Include dependency graph for DrawProductionComparison.C:

Go to the source code of this file.

Functions
void	ProcessTrigger (TString trigName="default", TString fileName="AnalysisResults.root", TString listName="Pi0IM_GammaTrackCorr_EMCAL")

void	Cluster (Int_t icalo)

void	ClusterCells (Int_t icalo)

void	Cell ()

void	Correl (Int_t icalo)

void	Isol (Int_t icalo)

void	ShowerShape (Int_t icalo)

void	InvMass (Int_t icalo, TString particle, TString fileName)

void	MCParticle (Int_t icalo, TString particle)

void	Track ()

void	Vertex ()

void	Centrality ()

void	PtHard ()
	Pt Hard distributions. More...

TObject *	GetHisto (TString histoName, Int_t iprod)

Bool_t	GetFileAndList (TString fileName, TString listName, TString trigName)

void	DrawProductionComparison (TString listName="Pi0IM_GammaTrackCorr_EMCAL", TString fileName="AnalysisResults.root", TString fileFormat="pdf")

Variables
TString	prod [] = {"LHC17l3b_fast_TM","LHC17l4b_fast_TM"}
	productions to be compared, directory name More...

TString	prodLeg [] = {"Geant3, b FAST","Geant4, b FAST"}

const Int_t	nProd = 2

TDirectoryFile *	dir [nProd]
	total number of productions to compare More...

TList *	histArr [nProd]
	TDirectory file where lists per trigger are stored in train ouput. More...

TFile *	file [nProd]
	TList with histograms for a given trigger. More...

Float_t	nEvents [nProd]
	Input train file. More...

TString	histoTag = ""
	number of events container More...

TString	format = "eps"
	file names tag, basically the trigger and calorimeter combination More...

TString	errType = "B"
	output plots format: eps, pdf, etc. More...

Bool_t	scaled = kFALSE

Int_t	color [] ={kBlack,kRed,kBlue,kOrange+1,kYellow+1,kGreen+2,kCyan+1,kViolet,kMagenta+2,kGray}
	In case of pT hard binned productions where the scaling and cross section is already done, do not normalize by N events or get the Sumw2 of histogram. More...

Detailed Description

Plot N productions comparison analysis of QA histograms from EMCal PWG-GA wagon.

Macro to plot comparison of different distributions (spectra, correlations) produced in QA trains but different data or cuts on the same data Based on the plots provided by DrawAnaCaloTrackQA.C.

To execute: root -q -b -l DrawProductionComparison.C'("Pi0IM_GammaTrackCorr_EMCAL","AnalysisResults.root")'

The trigger name might change depending on the wagon / data type In simulations only the "default" case is available On data, there can be different triggers, depending on the period

default : min bias like triggers, kINT7, kINT1, kMB
EMCAL_L0: kEMC7 L0 EMCal
DCAL_L0 : kEMC7 L0 DCal
EMCAL_L1: kEMCEGA L1 EG1 EMCal
DCAL_L1 : kEMCEGA L1 EG1 DCal
EMCAL_L2: kEMCEGA L1 EG2 EMCal
DCAL_L2 : kEMCEGA L1 EG2 DCal A plot will be produced for each of the triggers, if they existed in the data.

The input files must be placed in different directories, each one defined in the string array "prod"m for ex.: TString prod[] = {"AOD142","AOD115","ESD"}; and provide a name to be used in the legeds TString prodLeg[] = {"AOD 142","AOD 115","ESD"}; Both have to be modified inside the macro.

The number of productions has to be specified const Int_t nProd = 3;

It is also possible to run comparisons done for pT hard binned productions. Since the scaling is usually already done, change the bool scaled to "true" to avoir overnormalization.

Author: Gustavo Conesa Balbastre Gusta.nosp@m.vo.C.nosp@m.onesa.nosp@m..Bal.nosp@m.bastr.nosp@m.e@ce.nosp@m.rn.ch, (LPSC-CNRS)

Definition in file DrawProductionComparison.C.

Function Documentation

void Cell ( )

Plot basic calorimeter Cell QA histograms. 2 canvases with 2-4 pads

cell spectra
number of cells distribution

Definition at line 793 of file DrawProductionComparison.C.

Referenced by ProcessTrigger().

void Centrality ( )

Centrality To be updated

Definition at line 2299 of file DrawProductionComparison.C.

Referenced by AliAnalysisTaskCRCZDC::GetZDCSpectraCorrList(), ProcessTrigger(), and AliAnalysisTaskZDCEP::UserExec().

void Cluster ( Int_t icalo )

Plot basic calorimeter QA histograms. 2 canvases with 2-4 pads

cluster spectra after matching and PID cuts, and their ratio to different productions
cluster track matching residuals in eta and phi directions, and their ratio to different productions

To be updated

Parameters

icalo 0 EMCal, 1 DCal

Definition at line 238 of file DrawProductionComparison.C.

Referenced by ProcessTrigger().

void ClusterCells ( Int_t icalo )

Plot basic shower shape

Parameters

icalo 0 EMCal, 1 DCal

Definition at line 976 of file DrawProductionComparison.C.

Referenced by ProcessTrigger().

void Correl ( Int_t icalo )

cluster-track correlation

To be updated

Parameters

icalo 0 EMCal, 1 DCal

Definition at line 1773 of file DrawProductionComparison.C.

Referenced by ProcessTrigger().

void DrawProductionComparison	(	TString	listName = `"Pi0IM_GammaTrackCorr_EMCAL"`,
		TString	fileName = `"AnalysisResults.root"`,
		TString	fileFormat = `"pdf"`
	)

Main method, produce the plots with the comparisons

Parameters

listName	Name of list with histograms in file (same for all productions)
fileName	File name (same for all productions)
fileFormat	define the type of figures: eps, pdf, etc.

Definition at line 126 of file DrawProductionComparison.C.

Bool_t GetFileAndList	(	TString	fileName,
		TString	listName,
		TString	trigName
	)

Access the file and list with histograms and number of analyzed events per each production.

Parameters

fileName	File name (same for all productions)
listName	Name of list with histograms in file
trigName	trigger case name

Definition at line 2486 of file DrawProductionComparison.C.

Referenced by ProcessTrigger().

TObject * GetHisto	(	TString	histoName,
		Int_t	iprod
	)

Check if the list is available, if not get the histo directly from file

Returns: the histogram with the provided name

Parameters

histoName	histogram name
iprod	production index

Definition at line 2548 of file DrawProductionComparison.C.

Referenced by Cell(), Centrality(), Cluster(), ClusterCells(), Correl(), Isol(), MCParticle(), PtHard(), ShowerShape(), Track(), and Vertex().

void InvMass	(	Int_t	icalo,
		TString	particle,
		TString	fileName
	)

invariant mass plots. Before InvMassFit must run.

To be updated

Parameters

icalo	0 EMCal, 1 DCal
particle	Pi0 or Eta

Definition at line 2021 of file DrawProductionComparison.C.

Referenced by AliAODRecoDecayHF4Prong::CutRhoMass(), AliRDHFCutsB0toDStarPi::DeltaInvMassB0Kpipipi(), AliAnalysisTaskSEB0toDStarPi::DeltaInvMassB0Kpipipi(), AliRDHFCutsB0toDStarPi::DeltaInvMassDStarKpipi(), AliAnalysisTaskSEB0toDStarPi::DeltaInvMassDStarKpipi(), AliAnalysisTaskPSHFE::FillDPhiHistos(), AliAnalysisTaskPSHFE::FillPhotoElecHistos(), AliAODRecoDecayHF2Prong::InvMassD0(), AliAODRecoDecayHF4Prong::InvMassD0(), AliAODRecoDecayHF2Prong::InvMassD0bar(), AliAODRecoDecayHF4Prong::InvMassD0bar(), AliAODRecoDecayHF3Prong::InvMassDplus(), AliAODRecoCascadeHF::InvMassDplustoK0spi(), AliAODRecoDecayHF3Prong::InvMassDsKKpi(), AliAODRecoDecayHF3Prong::InvMassDspiKK(), AliAODRecoCascadeHF::InvMassDstoK0sK(), AliAODRecoDecayHF2Prong::InvMassJPSIee(), AliAODRecoDecayHF3Prong::InvMassLcpiKp(), AliAODRecoDecayHF3Prong::InvMassLcpKpi(), AliAODRecoCascadeHF::InvMassLctoK0sP(), AliAODRecoCascadeHF::InvMassLctoLambdaPi(), and AliAODRecoCascadeHF3Prong::InvMassPiXiPi().

void Isol ( Int_t icalo )

isolation cone

To be updated

Parameters

icalo 0 EMCal, 1 DCal

Definition at line 1895 of file DrawProductionComparison.C.

Referenced by ProcessTrigger().

void MCParticle	(	Int_t	icalo,
		TString	particle
	)

Plot basic calorimeter QA histograms depending on MC origin of cluster

cluster spectra from particle
reconstructed minus generated energy of cluster particle

Parameters

icalo	0 EMCal, 1 DCal
particle	MC origin of cluster: Photon, PhotonPi0Decay, Pi0 (Merged), Electron, ...

Definition at line 1174 of file DrawProductionComparison.C.

void ProcessTrigger	(	TString	trigName,
		TString	fileName,
		TString	listName
	)

Produce the plots per trigger, options are:

default : min bias like triggers, kINT7, kINT1, kMB
EMCAL_L0: kEMC7 L0 EMCal
DCAL_L0 : kEMC7 L0 DCal
EMCAL_L1: kEMCEGA L1 EG1 EMCal
DCAL_L1 : kEMCEGA L1 EG1 DCal
EMCAL_L2: kEMCEGA L1 EG2 EMCal
DCAL_L2 : kEMCEGA L1 EG2 DCal

Input:

Parameters

trigName	trigger case name
fileName	File name (same for all productions)
listName	Name of list with histograms in file

Definition at line 164 of file DrawProductionComparison.C.

Referenced by DrawProductionComparison().

void PtHard ( )

Pt Hard distributions.

Definition at line 1691 of file DrawProductionComparison.C.

void ShowerShape ( Int_t icalo )

Plot basic shower shape

Parameters

icalo 0 EMCal, 1 DCal

Definition at line 1359 of file DrawProductionComparison.C.

Referenced by ProcessTrigger().

void Track ( )

Hybrid Tracks distributions To be updated

Definition at line 1494 of file DrawProductionComparison.C.

Referenced by AliAnalysisTaskEmcalQGTagging::Coronna(), AliAnalysisTaskEmcalHfeTagging::Coronna(), AliAnalysisTaskEmcalQGTagging::FillHistograms(), ProcessTrigger(), AliAnalysisTaskEmcalQGTagging::SelectTrigger(), AliAnalysisTaskEmcalHfeTagging::SelectTrigger(), AliAnalysisTaskSubJetFraction::SelectTriggerHadron(), AliAnalysisTaskRecoilJetYield::SelectTriggerHadron(), and Cumulants::SetTracks().

void Vertex ( )

x y z vertex distribution and ratios to different productions.

Definition at line 2387 of file DrawProductionComparison.C.

Referenced by ProcessTrigger().

Variable Documentation

Int_t color[] ={kBlack,kRed,kBlue,kOrange+1,kYellow+1,kGreen+2,kCyan+1,kViolet,kMagenta+2,kGray}

In case of pT hard binned productions where the scaling and cross section is already done, do not normalize by N events or get the Sumw2 of histogram.

pre-defined colors list

Definition at line 115 of file DrawProductionComparison.C.

Referenced by Cell(), Centrality(), Cluster(), ClusterCells(), Correl(), InvMass(), Isol(), MCParticle(), PtHard(), ShowerShape(), Track(), and Vertex().