[u/mrichter/AliRoot.git] / PWGPP / analysisQA / README

--------------------------------------------
Analysis QA 
Authur: Satyajit Jena <sjena@cern.ch>
Date: 23 Nov 2013
Last Update: Mon Mar 31 15:47:10 CEST 2014
Last Update: 27/06/2014 - LF
--------------------------------------------


README - This readme file

%+++++++++++++++++++++++++++++++++++++++++++
process.sh:

Script to run each post-processing macro on 
respective outputs 

Simple script to execute macros and produce output

  Requirement: the macro paths need to be defined
  ex:  CODE=$ALICE_ROOT/PWGPP/analysisQA

  Arguments for each macro should be change by
  collecting proper information from respective
  wagon owner.

  sh process.sh <output.root> <aod-number>
  ex: sh process.sh output.root 145 eps


%+++++++++++++++++++++++++++++++++++++++++++
processCFv2vsPt.C
 
Processing of Flow QA task in Analysis QA train


%+++++++++++++++++++++++++++++++++++++++++++
processDrawAnaCaloTrackQA.C

%+++++++++++++++++++++++++++++++++++++++++++
processDrawUDQA.C

%+++++++++++++++++++++++++++++++++++++++++++
processHFEQAtask.C

%+++++++++++++++++++++++++++++++++++++++++++
processJETrackQA.C

Processes the results of the JE_PWG4HighPtTrackQA wagon. Should be run with the following arguments:
TString strFileIn   = "AnalysisResults.root" : the path to the input file.
TString suffix      = "eps"                  : the suffix that determines the output format. For example: eps, pdf, png, etc.
Int_t cent          = 10                     : the centrality class of the results. Options: 0 (0-10%), 1 (10-30%), 2 (30-50%), 3 (50-80%) and 10 (0-80%). Usually they are published with class 10.
Int_t trig          = 1                      : linked to the trigger. Default depends on data period: lhc11h: 1, containing lhc12 or lhc13: 5, or else: 6. On doubt one can always referr to the table at the bottom of this paragraph in after checking the AnalysisResults.root file with a browser.
Bool_t bESD         = kFALSE                 : an obsolete variable. Might be removed in the future.
Int_t run           = 0                      : enables to add a run number to the plots. Not neccessary for general QA, and can be suppresed by giving the number 0.
const char *outfile ="JETrackQA_output.root" : the path to the place to store the plots as root file. This is in addition to the image files made with in the 'suffix' format.
List of triggers:
  if(trig==1) strTrigger = "kCentral";
  if(trig==2) strTrigger = "kSemiCentral";
  if(trig==3) strTrigger = "kMBkCentralkSemiCentral";
  if(trig==4) strTrigger = "kEMCEJE";
  if(trig==5) strTrigger = "kINT7";
  if(trig==6) strTrigger = "kMB";

%+++++++++++++++++++++++++++++++++++++++++++
processJETriggerQA.C

This is an old version of the macro. Use V2 if possible. Only use this for old results.

%+++++++++++++++++++++++++++++++++++++++++++
processJETriggerQA_V2.C

Processes the results of the JE_TriggerQAFullR020 wagon, and if present, can also be used for the results of the JE_TriggerQAFullR020_EMCEGA and JE_TriggerQAFullR020_EMCEJE wagon. Should be run with the following arguments:
TString strFileIn    = "AnalysisResults.root"      : the path to the input file.
TString suftype      ="eps"                        : the suffix that deternmines the output format. For example: eps, pdf, png, etc.
Float_t jetR         = 0.2                         : the jetR. Default train value is 0.2, but can otherwise be found by opening the AnalysisResults.root file with the browser.
Float_t minTrkPT     = 0.15                        : the lower Track momentum cut-off. Default train value is 0.15, but can otherwise be found by opening the AnalysisResults.root file with the browser.
Float_t minClusterET = 0.3                         : the lower Cluster energy cut-off. Default train value is 0.3, but can otherwise be found by opening the AnalysisResults.root file with the browser.
nt_t run            = 0                           : enables to add a run number to the plots. Not neccessary for general QA, and can be suppresed by giving the number 0.
TString trigsuffix   = ""                          : leave empty to open the results of the JE_TriggerQAFullR020 wagon, but or use "EJE" or "EGA" to use the triggered results otherwise.
const char* outfile  = "JETriggerQA_outfile.root"  : the path to the place to store the plots as root file. This is in addition to the image files made with in the 'suffix' format.

%+++++++++++++++++++++++++++++++++++++++++++
processJpsi2eeQAplots.C

%+++++++++++++++++++++++++++++++++++++++++++
processMakeQA2pc.C

%+++++++++++++++++++++++++++++++++++++++++++
processMultistrangeQA.C

Processes the results of the wagons based on AddTaskQAMultistrange.C.
A copy is present in the folder: /PWGLF/QATasks/post/PostProcessQAMultistrange.C .
Creates a "summary pdf" file with several pages. 
Arguments:
Int_t   icasType        = 0,                             // 0) Xi-   1) Xi+  2) Omega-  3) Omega+
Int_t   collidingsystem = 0,                             // 0) PbPb  1) pp   2) pPb
Bool_t  isMC            = kFALSE,                        // kTRUE-->MC and kFALSE-->Exp.
Char_t *fileDir         = ".",                           // Input file directory
Char_t *filein          = "AnalysisResults.root"         // Input file name
Further details in the macro.

%+++++++++++++++++++++++++++++++++++++++++++
processResonance.C

Processes the results of the wagons based on AddTaskQAPhi.C.
A copy is present in the folder: /PWGLF/QATasks/post/PostProcessQAPhi.C .
Creates a "summary pdf" file with several pages.
Arguments:
char* system          = "pp276",                                 // string giving the collision system ("pp276","pp7", or "PbPb276")
char* name_fin        = "AnalysisResults.root",                  // name of input file
char* name_fout       = "QAphi",                                 // base name of output files (without suffix)
char* name_list       = "RsnHistMini_Phi_PhiNsigma_KTPCnsig30",  // name of the list in fin that contains histograms (may be different)
char* name_hist_base  = "RsnMini_phi.PhiNsigma_KTPCnsig30"       // base name of the THnSparse histograms (may be different)
Further details in the macro.

%+++++++++++++++++++++++++++++++++++++++++++
processV0.C

Processes the results of the wagons based on AddTaskQAV0.C or AddTaskQAV0AOD.C.
A copy is present in the folder: /PWGLF/QATasks/post/PostProcessQAV0.C .
Creates a "summary pdf" file with several pages. 
Arguments:
Bool_t lAttemptInvMassFit = kTRUE           // if kTRUE attempt rough signal extraction 
The Macro run on the file: ./AnalysisResults.root .

%+++++++++++++++++++++++++++++++++++++++++++
processHighPtDeDx.C

Processes the results of the wagons based on AddTaskQAHighPtDeDx.C.
A copy is present in the folder: /PWGLF/QATasks/post/PostProcessQAHighPtDeDx.C .
How to run (use AliRoot because of AliXRDPROOFtoolkit):

gSystem->AddIncludePath("-I$ALICE_ROOT/TPC/Base")
gSystem->AddIncludePath("-I$ALICE_ROOT/PWGLF/SPECTRA/IdentifiedHighPt/macros")
gSystem->AddIncludePath("-I$ALICE_ROOT/PWGLF/SPECTRA/IdentifiedHighPt/grid")
gSystem->AddIncludePath("-I$ALICE_ROOT/PWGLF/SPECTRA/IdentifiedHighPt/lib")
gROOT->SetMacroPath(".:$ALICE_ROOT/PWGLF/SPECTRA/IdentifiedHighPt/macros:$ALICE_ROOT/PWGLF/SPECTRA/IdentifiedHighPt/grid:$ALICE_ROOT/PWGLF/SPECTRA/IdentifiedHighPt/lib/")
.L my_functions.C+
.L my_tools.C+
.L PostProcessQAHighPtDeDx.C+

PlotQA("FileRoot/AnalysisResults.root")

MakeFitsExternalData("FileRoot/AnalysisResults.root", "HistosForBB")
MakeFitsV0s("FileRoot/AnalysisResults.root", "HistosForBB/PrimaryElectrons.root", "HistosForBB",0)
MakeFitsV0s("FileRoot/AnalysisResults.root", "HistosForBB/PrimaryElectrons.root", "HistosForBB",1)
MakeFitsV0s("FileRoot/AnalysisResults.root", "HistosForBB/PrimaryElectrons.root", "HistosForBB",2)
MakeFitsV0s("FileRoot/AnalysisResults.root", "HistosForBB/PrimaryElectrons.root", "HistosForBB",3)
PlotParametrizations("HistosForBB")

FitDeDxVsP("FileRoot/AnalysisResults.root", 3.0, 10.0, 0, 6, 13, kTRUE,  0,  2, 0,1, 0, "HistosForBB/hres_0_5_02.root",27);
FitDeDxVsP("FileRoot/AnalysisResults.root", 3.0, 10.0, 0, 6, 13, kTRUE,  2,  4, 0,1, 0, "HistosForBB/hres_0_5_24.root",27);
FitDeDxVsP("FileRoot/AnalysisResults.root", 3.0, 10.0, 0, 6, 13, kTRUE,  4,  6, 0,1, 0, "HistosForBB/hres_0_5_46.root",27);
FitDeDxVsP("FileRoot/AnalysisResults.root", 3.0, 10.0, 0, 6, 13, kTRUE,  6,  8, 0,1, 0, "HistosForBB/hres_0_5_68.root",27);
MakeNSigmaPlot("FileRoot/AnalysisResults.root","fitparameters/MB/02_dataPbPb.root",2,50, 0,  "02_dataPbPb.root");
MakeNSigmaPlot("FileRoot/AnalysisResults.root","fitparameters/MB/24_dataPbPb.root",2,50, 1,  "24_dataPbPb.root");
MakeNSigmaPlot("FileRoot/AnalysisResults.root","fitparameters/MB/46_dataPbPb.root",2,50, 2,  "46_dataPbPb.root");
MakeNSigmaPlot("FileRoot/AnalysisResults.root","fitparameters/MB/68_dataPbPb.root",2,50, 3,  "68_dataPbPb.root");
PlotNSigma("nsigma_results")

ExtractUncPartFractvsP("FileRoot/AnalysisResults.root", 2, 50,0,0, "results/eta02","fractions");
ExtractUncPartFractvsP("FileRoot/AnalysisResults.root", 2, 50,0,1, "results/eta24","fractions");
ExtractUncPartFractvsP("FileRoot/AnalysisResults.root", 2, 50,0,2, "results/eta46","fractions");
ExtractUncPartFractvsP("FileRoot/AnalysisResults.root", 2, 50,0,3, "results/eta68","fractions");

Further details in the macro.


%+++++++++++++++++++++++++++++++++++++++++++
processProduceFastQA.C

%+++++++++++++++++++++++++++++++++++++++++++
Commit	Line	Data
c722c8a3	1	--------------------------------------------
	2	Analysis QA
	3	Authur: Satyajit Jena <sjena@cern.ch>
	4	Date: 23 Nov 2013
	5	Last Update: Mon Mar 31 15:47:10 CEST 2014
0b6cb8ab	6	Last Update: 27/06/2014 - LF
c722c8a3	7	--------------------------------------------
c722c8a3	8
ca14f2e9	9
c722c8a3	10	README - This readme file
	11
	12	%+++++++++++++++++++++++++++++++++++++++++++
	13	process.sh:
	14
	15	Script to run each post-processing macro on
	16	respective outputs
	17
	18	Simple script to execute macros and produce output
	19
	20	Requirement: the macro paths need to be defined
	21	ex: CODE=$ALICE_ROOT/PWGPP/analysisQA
	22
	23	Arguments for each macro should be change by
	24	collecting proper information from respective
	25	wagon owner.
	26
	27	sh process.sh <output.root> <aod-number>
	28	ex: sh process.sh output.root 145 eps
	29
	30
	31
	32	%+++++++++++++++++++++++++++++++++++++++++++
	33	processCFv2vsPt.C
	34
	35	Processing of Flow QA task in Analysis QA train
	36
	37
	38
	39	%+++++++++++++++++++++++++++++++++++++++++++
	40	processDrawAnaCaloTrackQA.C
	41
	42	%+++++++++++++++++++++++++++++++++++++++++++
	43	processDrawUDQA.C
	44
	45	%+++++++++++++++++++++++++++++++++++++++++++
	46	processHFEQAtask.C
	47
	48	%+++++++++++++++++++++++++++++++++++++++++++
	49	processJETrackQA.C
	50
78c54fd7	51	Processes the results of the JE_PWG4HighPtTrackQA wagon. Should be run with the following arguments:
	52	TString strFileIn = "AnalysisResults.root" : the path to the input file.
	53	TString suffix = "eps" : the suffix that determines the output format. For example: eps, pdf, png, etc.
	54	Int_t cent = 10 : the centrality class of the results. Options: 0 (0-10%), 1 (10-30%), 2 (30-50%), 3 (50-80%) and 10 (0-80%). Usually they are published with class 10.
	55	Int_t trig = 1 : linked to the trigger. Default depends on data period: lhc11h: 1, containing lhc12 or lhc13: 5, or else: 6. On doubt one can always referr to the table at the bottom of this paragraph in after checking the AnalysisResults.root file with a browser.
	56	Bool_t bESD = kFALSE : an obsolete variable. Might be removed in the future.
	57	Int_t run = 0 : enables to add a run number to the plots. Not neccessary for general QA, and can be suppresed by giving the number 0.
	58	const char *outfile ="JETrackQA_output.root" : the path to the place to store the plots as root file. This is in addition to the image files made with in the 'suffix' format.
	59	List of triggers:
	60	if(trig==1) strTrigger = "kCentral";
	61	if(trig==2) strTrigger = "kSemiCentral";
	62	if(trig==3) strTrigger = "kMBkCentralkSemiCentral";
	63	if(trig==4) strTrigger = "kEMCEJE";
	64	if(trig==5) strTrigger = "kINT7";
	65	if(trig==6) strTrigger = "kMB";
	66
c722c8a3	67	%+++++++++++++++++++++++++++++++++++++++++++
	68	processJETriggerQA.C
	69
78c54fd7	70	This is an old version of the macro. Use V2 if possible. Only use this for old results.
78c54fd7	71
c722c8a3	72	%+++++++++++++++++++++++++++++++++++++++++++
	73	processJETriggerQA_V2.C
	74
78c54fd7	75	Processes the results of the JE_TriggerQAFullR020 wagon, and if present, can also be used for the results of the JE_TriggerQAFullR020_EMCEGA and JE_TriggerQAFullR020_EMCEJE wagon. Should be run with the following arguments:
	76	TString strFileIn = "AnalysisResults.root" : the path to the input file.
	77	TString suftype ="eps" : the suffix that deternmines the output format. For example: eps, pdf, png, etc.
	78	Float_t jetR = 0.2 : the jetR. Default train value is 0.2, but can otherwise be found by opening the AnalysisResults.root file with the browser.
	79	Float_t minTrkPT = 0.15 : the lower Track momentum cut-off. Default train value is 0.15, but can otherwise be found by opening the AnalysisResults.root file with the browser.
	80	Float_t minClusterET = 0.3 : the lower Cluster energy cut-off. Default train value is 0.3, but can otherwise be found by opening the AnalysisResults.root file with the browser.
cd65ede4	81	nt_t run = 0 : enables to add a run number to the plots. Not neccessary for general QA, and can be suppresed by giving the number 0.
78c54fd7	82	TString trigsuffix = "" : leave empty to open the results of the JE_TriggerQAFullR020 wagon, but or use "EJE" or "EGA" to use the triggered results otherwise.
	83	const char* outfile = "JETriggerQA_outfile.root" : the path to the place to store the plots as root file. This is in addition to the image files made with in the 'suffix' format.
	84
c722c8a3	85	%+++++++++++++++++++++++++++++++++++++++++++
	86	processJpsi2eeQAplots.C
	87
	88	%+++++++++++++++++++++++++++++++++++++++++++
	89	processMakeQA2pc.C
	90
	91	%+++++++++++++++++++++++++++++++++++++++++++
	92	processMultistrangeQA.C
	93
0b6cb8ab	94	Processes the results of the wagons based on AddTaskQAMultistrange.C.
	95	A copy is present in the folder: /PWGLF/QATasks/post/PostProcessQAMultistrange.C .
	96	Creates a "summary pdf" file with several pages.
	97	Arguments:
	98	Int_t icasType = 0, // 0) Xi- 1) Xi+ 2) Omega- 3) Omega+
	99	Int_t collidingsystem = 0, // 0) PbPb 1) pp 2) pPb
	100	Bool_t isMC = kFALSE, // kTRUE-->MC and kFALSE-->Exp.
	101	Char_t *fileDir = ".", // Input file directory
	102	Char_t *filein = "AnalysisResults.root" // Input file name
	103	Further details in the macro.
	104
	105	%+++++++++++++++++++++++++++++++++++++++++++
	106	processResonance.C
	107
	108	Processes the results of the wagons based on AddTaskQAPhi.C.
	109	A copy is present in the folder: /PWGLF/QATasks/post/PostProcessQAPhi.C .
	110	Creates a "summary pdf" file with several pages.
	111	Arguments:
	112	char* system = "pp276", // string giving the collision system ("pp276","pp7", or "PbPb276")
	113	char* name_fin = "AnalysisResults.root", // name of input file
	114	char* name_fout = "QAphi", // base name of output files (without suffix)
	115	char* name_list = "RsnHistMini_Phi_PhiNsigma_KTPCnsig30", // name of the list in fin that contains histograms (may be different)
	116	char* name_hist_base = "RsnMini_phi.PhiNsigma_KTPCnsig30" // base name of the THnSparse histograms (may be different)
	117	Further details in the macro.
	118
	119	%+++++++++++++++++++++++++++++++++++++++++++
	120	processV0.C
	121
	122	Processes the results of the wagons based on AddTaskQAV0.C or AddTaskQAV0AOD.C.
	123	A copy is present in the folder: /PWGLF/QATasks/post/PostProcessQAV0.C .
	124	Creates a "summary pdf" file with several pages.
	125	Arguments:
	126	Bool_t lAttemptInvMassFit = kTRUE // if kTRUE attempt rough signal extraction
	127	The Macro run on the file: ./AnalysisResults.root .
	128
	129	%+++++++++++++++++++++++++++++++++++++++++++
	130	processHighPtDeDx.C
	131
	132	Processes the results of the wagons based on AddTaskQAHighPtDeDx.C.
	133	A copy is present in the folder: /PWGLF/QATasks/post/PostProcessQAHighPtDeDx.C .
	134	How to run (use AliRoot because of AliXRDPROOFtoolkit):
	135
	136	gSystem->AddIncludePath("-I$ALICE_ROOT/TPC/Base")
	137	gSystem->AddIncludePath("-I$ALICE_ROOT/PWGLF/SPECTRA/IdentifiedHighPt/macros")
	138	gSystem->AddIncludePath("-I$ALICE_ROOT/PWGLF/SPECTRA/IdentifiedHighPt/grid")
	139	gSystem->AddIncludePath("-I$ALICE_ROOT/PWGLF/SPECTRA/IdentifiedHighPt/lib")
	140	gROOT->SetMacroPath(".:$ALICE_ROOT/PWGLF/SPECTRA/IdentifiedHighPt/macros:$ALICE_ROOT/PWGLF/SPECTRA/IdentifiedHighPt/grid:$ALICE_ROOT/PWGLF/SPECTRA/IdentifiedHighPt/lib/")
	141	.L my_functions.C+
	142	.L my_tools.C+
	143	.L PostProcessQAHighPtDeDx.C+
	144
	145	PlotQA("FileRoot/AnalysisResults.root")
	146
	147	MakeFitsExternalData("FileRoot/AnalysisResults.root", "HistosForBB")
	148	MakeFitsV0s("FileRoot/AnalysisResults.root", "HistosForBB/PrimaryElectrons.root", "HistosForBB",0)
	149	MakeFitsV0s("FileRoot/AnalysisResults.root", "HistosForBB/PrimaryElectrons.root", "HistosForBB",1)
	150	MakeFitsV0s("FileRoot/AnalysisResults.root", "HistosForBB/PrimaryElectrons.root", "HistosForBB",2)
	151	MakeFitsV0s("FileRoot/AnalysisResults.root", "HistosForBB/PrimaryElectrons.root", "HistosForBB",3)
	152	PlotParametrizations("HistosForBB")
	153
	154	FitDeDxVsP("FileRoot/AnalysisResults.root", 3.0, 10.0, 0, 6, 13, kTRUE, 0, 2, 0,1, 0, "HistosForBB/hres_0_5_02.root",27);
	155	FitDeDxVsP("FileRoot/AnalysisResults.root", 3.0, 10.0, 0, 6, 13, kTRUE, 2, 4, 0,1, 0, "HistosForBB/hres_0_5_24.root",27);
	156	FitDeDxVsP("FileRoot/AnalysisResults.root", 3.0, 10.0, 0, 6, 13, kTRUE, 4, 6, 0,1, 0, "HistosForBB/hres_0_5_46.root",27);
	157	FitDeDxVsP("FileRoot/AnalysisResults.root", 3.0, 10.0, 0, 6, 13, kTRUE, 6, 8, 0,1, 0, "HistosForBB/hres_0_5_68.root",27);
158	MakeNSigmaPlot("FileRoot/AnalysisResults.root","fitparameters/MB/02_dataPbPb.root",2,50, 0, "02_dataPbPb.root");
159	MakeNSigmaPlot("FileRoot/AnalysisResults.root","fitparameters/MB/24_dataPbPb.root",2,50, 1, "24_dataPbPb.root");
160	MakeNSigmaPlot("FileRoot/AnalysisResults.root","fitparameters/MB/46_dataPbPb.root",2,50, 2, "46_dataPbPb.root");
161	MakeNSigmaPlot("FileRoot/AnalysisResults.root","fitparameters/MB/68_dataPbPb.root",2,50, 3, "68_dataPbPb.root");
162	PlotNSigma("nsigma_results")
163
164	ExtractUncPartFractvsP("FileRoot/AnalysisResults.root", 2, 50,0,0, "results/eta02","fractions");
165	ExtractUncPartFractvsP("FileRoot/AnalysisResults.root", 2, 50,0,1, "results/eta24","fractions");
166	ExtractUncPartFractvsP("FileRoot/AnalysisResults.root", 2, 50,0,2, "results/eta46","fractions");
167	ExtractUncPartFractvsP("FileRoot/AnalysisResults.root", 2, 50,0,3, "results/eta68","fractions");
168
169	Further details in the macro.
170
171
c722c8a3	172	%+++++++++++++++++++++++++++++++++++++++++++
	173	processProduceFastQA.C
	174
	175	%+++++++++++++++++++++++++++++++++++++++++++