AliPhysics  5364b50 (5364b50)

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 ()
 
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...
 
Int_t color [] ={kBlack,kRed,kBlue,kOrange+1,kYellow+1,kGreen+2,kCyan+1,kViolet,kMagenta+2,kGray}
 pre-defined colors list 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;

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 778 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 234 of file DrawProductionComparison.C.

Referenced by ProcessTrigger().

void ClusterCells ( Int_t  icalo)

Plot basic shower shape

Parameters
icalo0 EMCal, 1 DCal

Definition at line 958 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 1630 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 124 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 2339 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 2399 of file DrawProductionComparison.C.

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

void Isol ( Int_t  icalo)

isolation cone

To be updated

Parameters
icalo0 EMCal, 1 DCal

Definition at line 1750 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 1150 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 162 of file DrawProductionComparison.C.

Referenced by DrawProductionComparison().

void ShowerShape ( Int_t  icalo)

Plot basic shower shape

Parameters
icalo0 EMCal, 1 DCal

Definition at line 1303 of file DrawProductionComparison.C.

Referenced by ProcessTrigger().

void Vertex ( )

x y z vertex distribution and ratios to different productions.

Definition at line 2238 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}

pre-defined colors list

Definition at line 113 of file DrawProductionComparison.C.

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

TDirectoryFile* dir[nProd]

total number of productions to compare

Definition at line 103 of file DrawProductionComparison.C.

Referenced by GetFileAndList().

TString errType = "B"

output plots format: eps, pdf, etc.

Definition at line 110 of file DrawProductionComparison.C.

TFile* file[nProd]

TList with histograms for a given trigger.

Definition at line 105 of file DrawProductionComparison.C.

Referenced by GetFileAndList(), and GetHisto().

TString format = "eps"

file names tag, basically the trigger and calorimeter combination

Definition at line 109 of file DrawProductionComparison.C.

TList* histArr[nProd]

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

Definition at line 104 of file DrawProductionComparison.C.

TString histoTag = ""

number of events container

Definition at line 108 of file DrawProductionComparison.C.

Float_t nEvents[nProd]

Input train file.

Definition at line 106 of file DrawProductionComparison.C.

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

TString prod[] = {"LHC17l3b_fast_TM","LHC17l4b_fast_TM"}

productions to be compared, directory name

productions name used in legends

Definition at line 95 of file DrawProductionComparison.C.

Referenced by AliMultiDimVector::GetGlobalAddressFromIndices(), AliMultiDimVector::GetIndicesFromGlobalAddress(), and GetProdErr().

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

Definition at line 96 of file DrawProductionComparison.C.