AliPhysics  9c66e61 (9c66e61)

Get fraction of clusters with photon shape and large number of cells. More...

#include <TFile.h>
#include <TDirectoryFile.h>
#include <TList.h>
#include <TString.h>
#include <TH1D.h>
#include <TH2F.h>
#include <TH3F.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 <TGraphErrors.h>
#include <TROOT.h>
#include <TStyle.h>
#include "PlotUtils.C"
Include dependency graph for CalculateShowerShapeFractionPerNCellCut.C:

Go to the source code of this file.

Functions

void GetSMNumber (TString titleName, Int_t &totalSM, Int_t &firstSM, Int_t &ncol, Int_t &nrow)
 
void CalculateAndPlot (TString titleName, TString filePath, Int_t iana=0, TString tm="_TMDep", Bool_t plotRat=kFALSE, Bool_t debug=kFALSE)
 
void CompareProductions (Int_t nProd, TString *fileName, TString *legName, Int_t iana=0, Bool_t debug=kFALSE)
 
void CalculateShowerShapeFractionPerNCellCut (Bool_t doCalc=kTRUE, Bool_t doComp=kTRUE, TString tm="_TMDep", Bool_t plotRat=kFALSE, Bool_t debug=kFALSE)
 

Variables

TString anaName [] = {"NoIso","Iso","All"}
 
TString anaTitle [] = {"NOT Isolated","Isolated","Inclusive"}
 
const Int_t nHisto = 3
 
TString histoLeg [] = {"#it{n}^{w}_{cell}>1", "#it{n}^{w}_{cell} > 4","#it{n}^{w}_{cell} #leq 4"}
 
const Int_t nEBins = 18
 
Double_t binE [] = { 2, 3, 4, 5, 6, 7, 8,10,12,14,16,18,20,22, 25, 30, 35, 40,50}
 
Double_t binEErr [] = {0.5,0.5,0.5,0.5,0.5,0.5,0.5, 1, 1, 1, 1, 1, 1, 1,1.5, 2.5, 2.5, 2.5, 5}
 
const Int_t nShShMaxCut = 3
 
Float_t maxShShCut [] = {0.27,0.3,0.4}
 
Int_t shshRef = 1
 
Int_t color []
 
Int_t lineStyle []
 
Int_t marker [] = {20,20,20,21,24,24,24,24,24,24,24}
 
TString fileFormat = ".eps"
 

Detailed Description

Get fraction of clusters with photon shape and large number of cells.

Plot the shower shape distribution for isolated, non isolated and inclusive clusters depending if the cluster has or not more than 4 cells with weight. Calculate the fraction of such clusters over the total when they have a shower shape like the one of photons, for different selection cuts 0.1<M02<0.27,0.3,0.4 Put the fractions vs E in output files to be processed later for comparison in method CompareProductions()

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 CalculateShowerShapeFractionPerNCellCut.C.

Function Documentation

void CalculateAndPlot ( TString  titleName,
TString  filePath,
Int_t  iana = 0,
TString  tm = "_TMDep",
Bool_t  plotRat = kFALSE,
Bool_t  debug = kFALSE 
)

Open the input histogram file and recover the histogram corresponding to the cluster type selected do the plots and calculate the fraction of clusters and put them in output file.

Parameters
titleName: String with simplified production name, it must contain data period and trigger for data and calorimeter
filePath: String input file path
iana: Selected cluster type: 0-not isolated, 1-isolated, 2-inclusive (AliAnaClusterShapeStudies)
tm: String with track matching option used
plotRat: Bool to activate ratio plots.
debug: Bool to activate printfs.

Definition at line 122 of file CalculateShowerShapeFractionPerNCellCut.C.

Referenced by CalculateShowerShapeFractionPerNCellCut(), and GetSMNumber().

void CalculateShowerShapeFractionPerNCellCut ( Bool_t  doCalc = kTRUE,
Bool_t  doComp = kTRUE,
TString  tm = "_TMDep",
Bool_t  plotRat = kFALSE,
Bool_t  debug = kFALSE 
)

Main steering method. Execute here the different cluster selection options for not isolated, isolated or inclusive. Set also here the productions to be analyzed and call the methods doing the calculation of the fractions and previous plotting and/or the comparison of the fractions of different productions.

Parameters
doCalc: Execute method CalculateAndPlot()
doComp: Execute method CompareProductions()
tm: String with track matching option used
plotRat: Bool to activate ratio plots.
debug: Bool to activate printfs.

Definition at line 1343 of file CalculateShowerShapeFractionPerNCellCut.C.

Referenced by CompareProductions().

void CompareProductions ( Int_t  nProd,
TString fileName,
TString legName,
Int_t  iana = 0,
Bool_t  debug = kFALSE 
)

Open the files with output graphs in different productions from CalculateAndPlot and recover the graph corresponding to the cluster type selected and plot the production comparison

Parameters
nProd: Number of productions
fileName: String input file path, file with graphs
legName: String input production title
iana: Selected cluster type: 0-not isolated, 1-isolated, 2-inclusive (AliAnaClusterShapeStudies)
debug: Bool to activate printfs.

Definition at line 1065 of file CalculateShowerShapeFractionPerNCellCut.C.

Referenced by CalculateAndPlot(), and CalculateShowerShapeFractionPerNCellCut().

void GetSMNumber ( TString  titleName,
Int_t totalSM,
Int_t firstSM,
Int_t ncol,
Int_t nrow 
)

Get SM total number, first SM, and n col/rows depending data string name

Parameters
titleName: String with simplified production name, it must contain data period and trigger for data and calorimeter
totalSM: Number of SM
firstSM: 0 for EMCal, 12 for DCal
ncol: Number of canvas columns
nrow: Number of canvas rows

Definition at line 85 of file CalculateShowerShapeFractionPerNCellCut.C.

Referenced by CalculateAndPlot(), and CompareProductions().

Variable Documentation

TString anaName[] = {"NoIso","Iso","All"}

Definition at line 44 of file CalculateShowerShapeFractionPerNCellCut.C.

TString anaTitle[] = {"NOT Isolated","Isolated","Inclusive"}

Definition at line 45 of file CalculateShowerShapeFractionPerNCellCut.C.

Double_t binEErr[] = {0.5,0.5,0.5,0.5,0.5,0.5,0.5, 1, 1, 1, 1, 1, 1, 1,1.5, 2.5, 2.5, 2.5, 5}
Int_t color[]
Initial value:
= {1,4,2,6,8,9,kYellow-6,kCyan,kOrange+2,kViolet,kOrange-2,
4,2,6,8,9,kYellow-6,kCyan,kOrange+2,kViolet,kOrange-2}

Definition at line 66 of file CalculateShowerShapeFractionPerNCellCut.C.

Referenced by CalculateAndPlot().

TString fileFormat = ".eps"
TString histoLeg[] = {"#it{n}^{w}_{cell}>1", "#it{n}^{w}_{cell} > 4","#it{n}^{w}_{cell} #leq 4"}

Definition at line 49 of file CalculateShowerShapeFractionPerNCellCut.C.

Int_t lineStyle[]
Initial value:
= {1,1,1,1,1,1, 1, 1, 1,
1, 1,2,2,2,2,2, 2, 2, 2, 2,
2,2,2,2,2,2,2,2,2,2,2,2}

Definition at line 68 of file CalculateShowerShapeFractionPerNCellCut.C.

Referenced by CalculateAndPlot().

Int_t marker[] = {20,20,20,21,24,24,24,24,24,24,24}

Definition at line 71 of file CalculateShowerShapeFractionPerNCellCut.C.

Referenced by CalculateAndPlot().

Float_t maxShShCut[] = {0.27,0.3,0.4}
const Int_t nShShMaxCut = 3
Int_t shshRef = 1

Definition at line 63 of file CalculateShowerShapeFractionPerNCellCut.C.

Referenced by CalculateAndPlot().