PWG2/SPECTRA/macros/configProtonAnalysis.C

   1 //__________________________________________________//
   2 AliProtonAnalysis *GetProtonAnalysisObject(const char* analysisLevel = "ESD",
   3                                            Bool_t kAnalyzeMC = kTRUE,
   4                                            const char* esdAnalysisType = "Hybrid",
   5                                            const char* pidMode = "Bayesian") {
   6   //                                       UInt_t runNumberForOfflineTrigger = -1) {
   7   gROOT->LoadMacro("$ALICE_ROOT/PWG2/SPECTRA/macros/configProtonAnalysisBaseObject.C");
   8   //Function to setup the AliProtonAnalysis object and return it
   9   AliProtonAnalysisBase *baseAnalysis = GetProtonAnalysisBaseObject(analysisLevel,kAnalyzeMC,esdAnalysisType,pidMode);
  10
  11   AliProtonAnalysis *analysis = new AliProtonAnalysis();
  12   analysis->SetBaseAnalysis(baseAnalysis);
  13   if(baseAnalysis->GetAnalysisMode() == AliProtonAnalysisBase::kGlobal) {
  14     Double_t gY[17] = {-0.9,-0.75,-0.6,-0.5,-0.4,-0.3,-0.2,-0.1,0.0,0.1,0.2,0.3,0.4,0.5,0.6,0.75,0.9};
  15     Double_t gPt[18] = {0.4,0.5,0.6,0.7,0.8,0.9,1.0,1.1,1.2,1.3,1.4,1.5,1.7,1.9,2.1,2.4,2.7,3.1};
  16     analysis->InitAnalysisHistograms(16,gY,17,gPt);
  17     //Double_t gPt[27] = {0.4,0.5,0.6,0.7,0.8,0.9,1.0,1.1,1.2,1.3,1.4,1.5,1.7,1.9,2.1,2.4,2.7,3.0,3.4,3.8,4.2,4.7,5.1,6.0,7.0,9.0,12.0};
  18     //analysis->InitAnalysisHistograms(16,gY,26,gPt);
  19   }
  20   else
  21     analysis->InitAnalysisHistograms(baseAnalysis->GetNBinsX(),
  22                                      baseAnalysis->GetMinX(),
  23                                      baseAnalysis->GetMaxX(),
  24                                      baseAnalysis->GetNBinsY(),
  25                                      baseAnalysis->GetMinY(),
  26                                      baseAnalysis->GetMaxY());
  27
  28   return analysis;
  29 }
  30
  31 //__________________________________________________//
  32 AliProtonQAAnalysis *GetProtonQAAnalysisObject(const char* analysisLevel = "ESD",
  33                                                const char* esdAnalysisType = "Hybrid",
  34                                                const char* pidMode = "Bayesian") {
  35   gROOT->LoadMacro("$ALICE_ROOT/PWG2/SPECTRA/macros/configProtonAnalysisBaseObject.C");
  36   //Function to setup the AliProtonQAAnalysis object and return it
  37   AliProtonAnalysisBase *baseAnalysis = GetProtonAnalysisBaseObject(analysisLevel,kTRUE,esdAnalysisType,pidMode);
  38
  39   AliProtonQAAnalysis *analysis = new AliProtonQAAnalysis();
  40   analysis->SetBaseAnalysis(baseAnalysis);
  41   analysis->SetRunMCAnalysis();
  42   //analysis->SetMCProcessId(4);//4: weak decay - 13: hadronic interaction
  43   //analysis->SetMotherParticlePDGCode(3122);//3122: Lambda
  44   analysis->SetRunEfficiencyAnalysis(kFALSE);//use cuts in the eff. analysis
  45   analysis->SetQAYPtBins(baseAnalysis->GetNBinsX(),
  46                          baseAnalysis->GetMinX(),
  47                          baseAnalysis->GetMaxX(),
  48                          baseAnalysis->GetNBinsY(),
  49                          baseAnalysis->GetMinY(),
  50                          baseAnalysis->GetMaxY());
  51
  52   return analysis;
  53 }
  54
  55 //__________________________________________________//
  56