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
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)
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();
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);
28 mPM->AddHisto(histname);
29 mPM->SetInputFile(finname);
30 mPM->SetOutputFile(foutname);
31 mPM->SetMultLowBound(mult);
32 mPM->SetPercentileCrossSection(percentXsec);
33 mPM->MakePercentiles();