PWG2/EVCHAR/makeCentralityBy1D.C

   1 void makeCentralityBy1D
   2 //(int run =167693, const char *histname ="fHOutMultV0M", float percentXsec=90.05, float mult=77.4) //V0M
   3 //(int run =167693, const char *histname ="fHOutMultCL1", float percentXsec=90.14, float mult=20) // CL1
   4 (int run =167693, const char *histname ="fHOutMultTRK", float percentXsec=90.09, float mult=10) // TRK
   5 {
   6  //load libraries
   7   // gSystem->SetBuildDir("/tmp/");
   8   // gSystem->Load("libCore.so");
   9   // gSystem->Load("libTree.so");
  10   // gSystem->Load("libGeom.so");
  11   // gSystem->Load("libVMC.so");
  12   // gSystem->Load("libPhysics.so");
  13   // gSystem->Load("libSTEERBase.so");
  14   // gROOT->ProcessLine(".include $ALICE_ROOT/include");
  15   gROOT->LoadMacro("AliCentralityBy1D.cxx+");
  16   AliCentralityBy1D *mPM = new AliCentralityBy1D();
  17
  18   TString finname = Form("/home/atoia/analysis/data2011/multRef/AnalysisResults_%i.root",run);
  19   TString foutname = Form("/home/atoia/analysis/data2011/by1D/AliCentralityBy1D_%i_TRK.root",run);
  20
  21   mPM->AddHisto(histname);
  22   mPM->SetInputFile(finname);
  23   mPM->SetOutputFile(foutname);
  24   mPM->SetMultLowBound(mult);
  25   mPM->SetPercentileCrossSection(percentXsec);
  26   mPM->MakePercentiles();
  27
  28 }
  29