AliPhysics  e6d2b2b (e6d2b2b)

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...
 
TObjectGetHisto (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...
 
TListhistArr [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 ( )
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
icalo0 EMCal, 1 DCal

Definition at line 238 of file DrawProductionComparison.C.

Referenced by ProcessTrigger().

void ClusterCells ( Int_t  icalo)

Plot basic shower shape

Parameters
icalo0 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
icalo0 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
listNameName of list with histograms in file (same for all productions)
fileNameFile name (same for all productions)
fileFormatdefine 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
fileNameFile name (same for all productions)
listNameName of list with histograms in file
trigNametrigger 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
histoNamehistogram name
iprodproduction index

Definition at line 2548 of file DrawProductionComparison.C.

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

void Isol ( Int_t  icalo)

isolation cone

To be updated

Parameters
icalo0 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
icalo0 EMCal, 1 DCal
particleMC 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
trigNametrigger case name
fileNameFile name (same for all productions)
listNameName 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
icalo0 EMCal, 1 DCal

Definition at line 1359 of file DrawProductionComparison.C.

Referenced by ProcessTrigger().

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().

TDirectoryFile* dir[nProd]

total number of productions to compare

Definition at line 104 of file DrawProductionComparison.C.

Referenced by GetFileAndList().

TString errType = "B"

output plots format: eps, pdf, etc.

Definition at line 111 of file DrawProductionComparison.C.

TFile* file[nProd]

TList with histograms for a given trigger.

Definition at line 106 of file DrawProductionComparison.C.

Referenced by GetFileAndList(), and GetHisto().

TString format = "eps"

file names tag, basically the trigger and calorimeter combination

Definition at line 110 of file DrawProductionComparison.C.

TList* histArr[nProd]

TDirectory file where lists per trigger are stored in train ouput.

Definition at line 105 of file DrawProductionComparison.C.

TString histoTag = ""

number of events container

Definition at line 109 of file DrawProductionComparison.C.

Float_t nEvents[nProd]

Input train file.

Definition at line 107 of file DrawProductionComparison.C.

Referenced by Cell(), Centrality(), Cluster(), ClusterCells(), GetFileAndList(), MCParticle(), ShowerShape(), Track(), and Vertex().