AliPhysics  9b6b435 (9b6b435)

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