]>
Commit | Line | Data |
---|---|---|
e6f3f2fe | 1 | { |
2 | //load libraries | |
3 | //gSystem->SetBuildDir("/tmp"); | |
d20baea6 | 4 | gSystem->Load("libCore.so"); |
5 | gSystem->Load("libTree.so"); | |
6 | gSystem->Load("libGeom.so"); | |
7 | gSystem->Load("libVMC.so"); | |
8 | gSystem->Load("libPhysics.so"); | |
9 | gSystem->Load("libSTEERBase"); | |
10 | gROOT->ProcessLine(".include $ALICE_ROOT/include"); | |
e6f3f2fe | 11 | gROOT->LoadMacro("AliCentralityGlauberFit.cxx+"); |
d20baea6 | 12 | |
13 | const char *finname ="../centrality/cvetan/VZero_366.root"; // name input file | |
14 | const char *foutname="test.root"; // name output file | |
e6f3f2fe | 15 | |
d20baea6 | 16 | float percentXsec=100.0; |
17 | ||
e6f3f2fe | 18 | AliCentralityGlauberFit *mPM = new AliCentralityGlauberFit(); |
e6f3f2fe | 19 | mPM.SetOutputFile(foutname); |
20 | ||
d20baea6 | 21 | // latest V0 cvetan corrected |
22 | mPM.AddHisto("fCorrMult1d"); | |
23 | //mPM.SetGlauberParam(1,27.550,28., 10,0.4,0.8, 1,1.116,0.825, 1,0.974,0.98); //minuit Npart**alpha | |
24 | mPM.SetGlauberParam(1,23.5,25., 1,0.3,1.5, 1,1.115,0.83, 1,0.99,0.99); | |
25 | mPM.SetRebin(10); | |
26 | mPM.SetRangeToFit(100., 21000.); | |
27 | mPM.MakeFitsMinuitNBD(finname); | |
28 | //mPM->MakeFits(finname); | |
e6f3f2fe | 29 | |
d20baea6 | 30 | TFile * f = new TFile (foutname); |
31 | gDirectory->Get("fCorrMult1d")->Draw();; | |
32 | TH1 * hg = (TH1*) gDirectory->Get("fCorrMult1d_GLAU"); | |
33 | hg->SetLineColor(kRed); | |
34 | hg->Draw("same"); | |
e6f3f2fe | 35 | } |
d20baea6 | 36 |