CompareTH3DataAndMCProd.C File Reference

Compare a cluster parameter (shower shape) between different productions. More...

#include <TFile.h>
#include <TDirectoryFile.h>
#include <TList.h>
#include <TString.h>
#include <TROOT.h>
#include <TStyle.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 <TArrayD.h>
#include <TGraphErrors.h>
#include "PlotUtils.C"

Include dependency graph for CompareTH3DataAndMCProd.C:

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Functions
void	DefineHistogramSettings (TString histoName, Float_t &xmin, Float_t &xmax, Float_t &xNormMin, Float_t &xNormMax, Bool_t &normtomax, Int_t &rebin)
Bool_t	SetHistogramRangesAndTitles (TH1D *histo, TString histoName, Int_t iprod, Int_t rebin, Float_t xmin, Float_t xmax, Bool_t normtomax, Float_t xNormMin, Float_t xNormMax, Float_t lowEBin, Float_t highEBin, Bool_t debug)
void	CompareTH3DataAndMCProd (TString &outputFileName, TString histoName="SMM02NoCut", const Int_t nProd=3, TString *prod=0x0, TString *prodLeg=0x0, TString clusterization="", TString tm="_TMDep", TString pid="_Neutral", TString titleMC="JJDecLow", TString mcLeg="MC: #gammaJ+JJ({p^{EMCal}_{T,#gamma}>3.5 GeV/#it{c})", TString titleData="LHC11cd_EMC7", TString daLeg="pp@7 TeV, LHC11c+d EMC7", Int_t firstP=0, Int_t lastP=9, TArrayD binE=0, Int_t firstMC=1, Bool_t plotRatio=kFALSE, Bool_t saveHisto=kFALSE, TString opt="", Bool_t debug=kFALSE)

Detailed Description

Compare a cluster parameter (shower shape) between different productions.

Compare different shape parameter histograms in data and MC productions depending on the SM per cluster energy bis. Treatment of the output of the class AliAnaClusterShapeCorrelStudies Input are TH3 histograms where x=energy, y= SM number or T-Card index and z is one cluster parameter See the different options in MakeDataMCComparisonPerSMClusterEbin.C Here plots are made or projected distributions are stored on another file for further processing like in CalculateParamChi2MCxTalkDataPerSM.C

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

Function Documentation

void CompareTH3DataAndMCProd	(	TString &	outputFileName,
		TString	histoName = "SMM02NoCut",
		const Int_t	nProd = 3,
		TString *	prod = 0x0,
		TString *	prodLeg = 0x0,
		TString	clusterization = "",
		TString	tm = "_TMDep",
		TString	pid = "_Neutral",
		TString	titleMC = "JJDecLow",
		TString	mcLeg = "MC: #gammaJ+JJ({p^{EMCal}_{T,#gamma}>3.5 GeV/#it{c})",
		TString	titleData = "LHC11cd_EMC7",
		TString	daLeg = "pp@7 TeV, LHC11c+d EMC7",
		Int_t	firstP = 0,
		Int_t	lastP = 9,
		TArrayD	binE = 0,
		Int_t	firstMC = 1,
		Bool_t	plotRatio = kFALSE,
		Bool_t	saveHisto = kFALSE,
		TString	opt = "",
		Bool_t	debug = kFALSE
	)

Main method

Parameters

outputFileName	: String with part of output file name to be recovered and used by another macro, see MakeDataMCComparisonPerSMClusterEbinAndChi2.C
histoName	: Name of histogram with the parameter to study = "SMM02NoCut", "SMM02NoCut","SMM02", "SMM20LowM02NoCut","SMM20LowM02","SMM20HighM02NoCut","SMM20HighM02","SMNCell" and other see MakeDataMCComparisonPerSMClusterEbin.C
nProd	: Number of analyzed productions
prod	: Strings array with input files location
prodLeg	: Strings array with title name of analized productions
clusterization	: String with clusterizer used, in case it is included in the input data file name
tm	: String with track matching option used
pid	: String with PID cut used, see AliAnaClusterShapeAnalysis
titleMC	: Simplified string acronym of input MC
mcLeg	: Full title string of input MC
titleData	: Simplified string acronym of input data
daLeg	: Full title string of input data
firstP	: First SM or T-Card index to be studied
lastP	: Last SM or T-Card index to be studied
binE	: TArrayD with list of energy bin limits
firstMC	: prod should contain first data and then MC arrays, this int indicates which one is the first MC
plotRatio	: Optional plotting of different ratio histograms
saveHisto	: Optional storing of projected histograms into output file, if active, no plots are done
opt	: Additional string to differenciate output histograms for different analysis
debug	: Bool to activate printfs.

Definition at line 81 of file CompareTH3DataAndMCProd.C.

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

void DefineHistogramSettings	(	TString	histoName,
		Float_t &	xmin,
		Float_t &	xmax,
		Float_t &	xNormMin,
		Float_t &	xNormMax,
		Bool_t &	normtomax,
		Int_t &	rebin
	)

Define the histogram x axis range, normalization and rebin, depending the cluster parameter

Definition at line 982 of file CompareTH3DataAndMCProd.C.

Referenced by CompareTH3DataAndMCProd().

Bool_t SetHistogramRangesAndTitles	(	TH1D *	histo,
		TString	histoName,
		Int_t	iprod,
		Int_t	rebin,
		Float_t	xmin,
		Float_t	xmax,
		Bool_t	normtomax,
		Float_t	xNormMin,
		Float_t	xNormMax,
		Float_t	lowEBin,
		Float_t	highEBin,
		Bool_t	debug
	)

Prepare the histogram for plotting, set the axis ranges, titles, line and markers style colors, etc.

Definition at line 1074 of file CompareTH3DataAndMCProd.C.

Referenced by CompareTH3DataAndMCProd(), and DefineHistogramSettings().