PWGPP/analysisQA/README

   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
   6 Last Update: 27/06/2014 - LF
   7 --------------------------------------------
   8
   9
  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
  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 strFileIn2  = ""                     : the path to a second input file (the train output file), that will be used to compare the graphs with. In the plots this data will be labelled 'previous production. An empty string will give the bare plots without comparisons.
  54 TString suffix      = "eps"                  : the suffix that determines the output format. For example: eps, pdf, png, etc.
  55 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.
  56 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.
  57 Bool_t bESD         = kFALSE                 : an obsolete variable. Might be removed in the future.
  58 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.
  59 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.
  60 List of triggers:
  61   if(trig==1) strTrigger = "kCentral";
  62   if(trig==2) strTrigger = "kSemiCentral";
  63   if(trig==3) strTrigger = "kMBkCentralkSemiCentral";
  64   if(trig==4) strTrigger = "kEMCEJE";
  65   if(trig==5) strTrigger = "kINT7";
  66   if(trig==6) strTrigger = "kMB";
  67
  68 %+++++++++++++++++++++++++++++++++++++++++++
  69 processJETriggerQA.C
  70
  71 This is an old version of the macro. Use V2 if possible. Only use this for old results.
  72
  73 %+++++++++++++++++++++++++++++++++++++++++++
  74 processJETriggerQA_V2.C
  75
  76 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:
  77 TString strFileIn    = "AnalysisResults.root"      : the path to the input file.
  78 TString strFileIn2   = ""                          : the path to a second input file (the train output file), that will be used to compare the graphs with. In the plots this data will be labelled 'previous production. An empty string will give the bare plots without comparisons.
  79 TString suftype      ="eps"                        : the suffix that deternmines the output format. For example: eps, pdf, png, etc.
  80 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.
  81 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.
  82 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.
  83 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.
  84 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.
  85 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.
  86
  87 %+++++++++++++++++++++++++++++++++++++++++++
  88 processJpsi2eeQAplots.C
  89
  90 %+++++++++++++++++++++++++++++++++++++++++++
  91 processMakeQA2pc.C
  92
  93 %+++++++++++++++++++++++++++++++++++++++++++
  94 processMultistrangeQA.C
  95
  96 Processes the results of the wagons based on AddTaskQAMultistrange.C.
  97 A copy is present in the folder: /PWGLF/QATasks/post/PostProcessQAMultistrange.C .
  98 Creates a "summary pdf" file with several pages.
  99 Arguments:
 100 Int_t   icasType        = 0,                             // 0) Xi-   1) Xi+  2) Omega-  3) Omega+
 101 Int_t   collidingsystem = 0,                             // 0) PbPb  1) pp   2) pPb
 102 Bool_t  isMC            = kFALSE,                        // kTRUE-->MC and kFALSE-->Exp.
 103 Char_t *fileDir         = ".",                           // Input file directory
 104 Char_t *filein          = "AnalysisResults.root"         // Input file name
 105 Further details in the macro.
 106
 107 %+++++++++++++++++++++++++++++++++++++++++++
 108 processResonance.C
 109
 110 Processes the results of the wagons based on AddTaskQAPhi.C.
 111 A copy is present in the folder: /PWGLF/QATasks/post/PostProcessQAPhi.C .
 112 Creates a "summary pdf" file with several pages.
 113 Arguments:
 114 char* system          = "pp276",                                 // string giving the collision system ("pp276","pp7", or "PbPb276")
 115 char* name_fin        = "AnalysisResults.root",                  // name of input file
 116 char* name_fout       = "QAphi",                                 // base name of output files (without suffix)
 117 char* name_list       = "RsnHistMini_Phi_PhiNsigma_KTPCnsig30",  // name of the list in fin that contains histograms (may be different)
 118 char* name_hist_base  = "RsnMini_phi.PhiNsigma_KTPCnsig30"       // base name of the THnSparse histograms (may be different)
 119 Further details in the macro.
 120
 121 %+++++++++++++++++++++++++++++++++++++++++++
 122 processV0.C
 123
 124 Processes the results of the wagons based on AddTaskQAV0.C or AddTaskQAV0AOD.C.
 125 A copy is present in the folder: /PWGLF/QATasks/post/PostProcessQAV0.C .
 126 Creates a "summary pdf" file with several pages.
 127 Arguments:
 128 Bool_t lAttemptInvMassFit = kTRUE           // if kTRUE attempt rough signal extraction
 129 The Macro run on the file: ./AnalysisResults.root .
 130
 131 %+++++++++++++++++++++++++++++++++++++++++++
 132 processHighPtDeDx.C
 133
 134 Processes the results of the wagons based on AddTaskQAHighPtDeDx.C.
 135 A copy is present in the folder: /PWGLF/QATasks/post/PostProcessQAHighPtDeDx.C .
 136 How to run (use AliRoot because of AliXRDPROOFtoolkit):
 137
 138 gSystem->AddIncludePath("-I$ALICE_ROOT/TPC/Base")
 139 gSystem->AddIncludePath("-I$ALICE_ROOT/PWGLF/SPECTRA/IdentifiedHighPt/macros")
 140 gSystem->AddIncludePath("-I$ALICE_ROOT/PWGLF/SPECTRA/IdentifiedHighPt/grid")
 141 gSystem->AddIncludePath("-I$ALICE_ROOT/PWGLF/SPECTRA/IdentifiedHighPt/lib")
 142 gROOT->SetMacroPath(".:$ALICE_ROOT/PWGLF/SPECTRA/IdentifiedHighPt/macros:$ALICE_ROOT/PWGLF/SPECTRA/IdentifiedHighPt/grid:$ALICE_ROOT/PWGLF/SPECTRA/IdentifiedHighPt/lib/")
 143 .L my_functions.C+
 144 .L my_tools.C+
 145 .L PostProcessQAHighPtDeDx.C+
 146
 147 PlotQA("FileRoot/AnalysisResults.root")
 148
 149 MakeFitsExternalData("FileRoot/AnalysisResults.root", "HistosForBB")
 150 MakeFitsV0s("FileRoot/AnalysisResults.root", "HistosForBB/PrimaryElectrons.root", "HistosForBB",0)
 151 MakeFitsV0s("FileRoot/AnalysisResults.root", "HistosForBB/PrimaryElectrons.root", "HistosForBB",1)
 152 MakeFitsV0s("FileRoot/AnalysisResults.root", "HistosForBB/PrimaryElectrons.root", "HistosForBB",2)
 153 MakeFitsV0s("FileRoot/AnalysisResults.root", "HistosForBB/PrimaryElectrons.root", "HistosForBB",3)
 154 PlotParametrizations("HistosForBB")
 155
 156 FitDeDxVsP("FileRoot/AnalysisResults.root", 3.0, 10.0, 0, 6, 13, kTRUE,  0,  2, 0,1, 0, "HistosForBB/hres_0_5_02.root",27);
 157 FitDeDxVsP("FileRoot/AnalysisResults.root", 3.0, 10.0, 0, 6, 13, kTRUE,  2,  4, 0,1, 0, "HistosForBB/hres_0_5_24.root",27);
 158 FitDeDxVsP("FileRoot/AnalysisResults.root", 3.0, 10.0, 0, 6, 13, kTRUE,  4,  6, 0,1, 0, "HistosForBB/hres_0_5_46.root",27);
 159 FitDeDxVsP("FileRoot/AnalysisResults.root", 3.0, 10.0, 0, 6, 13, kTRUE,  6,  8, 0,1, 0, "HistosForBB/hres_0_5_68.root",27);
 160 MakeNSigmaPlot("FileRoot/AnalysisResults.root","fitparameters/MB/02_dataPbPb.root",2,50, 0,  "02_dataPbPb.root");
 161 MakeNSigmaPlot("FileRoot/AnalysisResults.root","fitparameters/MB/24_dataPbPb.root",2,50, 1,  "24_dataPbPb.root");
 162 MakeNSigmaPlot("FileRoot/AnalysisResults.root","fitparameters/MB/46_dataPbPb.root",2,50, 2,  "46_dataPbPb.root");
 163 MakeNSigmaPlot("FileRoot/AnalysisResults.root","fitparameters/MB/68_dataPbPb.root",2,50, 3,  "68_dataPbPb.root");
 164 PlotNSigma("nsigma_results")
 165
 166 ExtractUncPartFractvsP("FileRoot/AnalysisResults.root", 2, 50,0,0, "results/eta02","fractions");
 167 ExtractUncPartFractvsP("FileRoot/AnalysisResults.root", 2, 50,0,1, "results/eta24","fractions");
 168 ExtractUncPartFractvsP("FileRoot/AnalysisResults.root", 2, 50,0,2, "results/eta46","fractions");
 169 ExtractUncPartFractvsP("FileRoot/AnalysisResults.root", 2, 50,0,3, "results/eta68","fractions");
 170
 171 Further details in the macro.
 172
 173
 174 %+++++++++++++++++++++++++++++++++++++++++++
 175 processProduceFastQA.C
 176
 177 %+++++++++++++++++++++++++++++++++++++++++++