PWGPP/EVCHAR/makeCentralityBy1D.C

   1 void makeCentralityBy1D
   2 //(int run =167987, const char *system ="V0M", float percentXsec=90.06, float mult=77.4) //V0M
   3 //(int run =170162, const char *system ="CL1", float percentXsec=90.04, float mult=20) // CL1
   4 //(int run =167693, const char *system ="TRK", float percentXsec=90.20, float mult=11) // TRK
   5
   6 // for pA pilot run
   7 // for MC use percentXsec=100 and mult=0
   8 // for MC available also NPA
   9 // for pA use V0M, V0A, V0C, CL1, TRK, CND, (tbd FMD: not available in 146079)
  10 (int run =195344, const char *system ="FMD", float percentXsec=100.0, float mult=0.0)
  11 {
  12  //load libraries
  13   // gSystem->SetBuildDir("/tmp/");
  14   // gSystem->Load("libCore.so");
  15   // gSystem->Load("libTree.so");
  16   // gSystem->Load("libGeom.so");
  17   // gSystem->Load("libVMC.so");
  18   // gSystem->Load("libPhysics.so");
  19   // gSystem->Load("libSTEERBase.so");
  20   // gROOT->ProcessLine(".include $ALICE_ROOT/include");
  21   gROOT->LoadMacro("AliCentralityBy1D.cxx+");
  22   AliCentralityBy1D *mPM = new AliCentralityBy1D();
  23
  24   TString finname = Form("/home/atoia/analysis/data2013/multRef/EventStat_temp_%i.root",run);
  25   TString foutname = Form("/home/atoia/analysis/data2013/by1D/AliCentralityBy1D_%i_%s.root",run,system);
  26   const char *histname=Form("fHOutMult%s",system);
  27
  28   mPM->AddHisto(histname);
  29   mPM->SetInputFile(finname);
  30   mPM->SetOutputFile(foutname);
  31   mPM->SetMultLowBound(mult);
  32   mPM->SetPercentileCrossSection(percentXsec);
  33   mPM->MakePercentiles();
  34
  35 }
  36