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 binE[] = { 2, 3, 4, 5, 6, 7, 8,10,12,14,16,18,20,22, 25, 30, 35, 40,50}

Definition at line 53 of file CalculateShowerShapeFractionPerNCellCut.C.

Referenced by CalculateAndPlot(), CalculateParamChi2MCxTalkDataPerSM(), CompareTH3DataAndMCProd(), MakeDataMCComparisonPerSMClusterEbin(), MakeDataMCComparisonPerSMClusterEbinAndChi2(), MakeDoubleRatios(), and MakeNDataNMCComparisonPerSMClusterEbin().

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}

Definition at line 54 of file CalculateShowerShapeFractionPerNCellCut.C.

Referenced by CalculateParamChi2MCxTalkDataPerSM(), MakeDataMCComparisonPerSMClusterEbin(), and MakeDataMCComparisonPerSMClusterEbinAndChi2().

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"

Definition at line 74 of file CalculateShowerShapeFractionPerNCellCut.C.

Referenced by CalculateAndPlot(), and CompareProductions().

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}

Definition at line 58 of file CalculateShowerShapeFractionPerNCellCut.C.

Referenced by CalculateAndPlot(), and CompareProductions().

const Int_t nEBins = 18

Definition at line 52 of file CalculateShowerShapeFractionPerNCellCut.C.

Referenced by AddChi2JJDecLowHighGraphs(), CalculateAndPlot(), CalculateParamChi2MCxTalkDataPerSM(), CompareTH3DataAndMCProd(), MakeDataMCComparisonPerSMClusterEbin(), MakeDataMCComparisonPerSMClusterEbinAndChi2(), and MakeNDataNMCComparisonPerSMClusterEbin().

const Int_t nHisto = 3

Definition at line 48 of file CalculateShowerShapeFractionPerNCellCut.C.

Referenced by CalculateAndPlot(), MakeDataMCComparisonPerSMClusterEbinAndChi2(), MakeNDataNMCComparisonPerSMClusterEbin(), and ScalerTrending().

const Int_t nShShMaxCut = 3

Definition at line 57 of file CalculateShowerShapeFractionPerNCellCut.C.

Referenced by CalculateAndPlot(), and CompareProductions().

Int_t shshRef = 1

Definition at line 63 of file CalculateShowerShapeFractionPerNCellCut.C.

Referenced by CalculateAndPlot().

Functions

Variables

Detailed Description

Function Documentation

Variable Documentation