]>
Commit | Line | Data |
---|---|---|
61827aa9 | 1 | void makeCentralityFit() |
2 | { | |
3 | //load libraries | |
4 | gSystem->SetBuildDir("/tmp/"); | |
5 | gSystem->Load("libCore.so"); | |
6 | gSystem->Load("libTree.so"); | |
7 | gSystem->Load("libGeom.so"); | |
8 | gSystem->Load("libVMC.so"); | |
9 | gSystem->Load("libPhysics.so"); | |
10 | gSystem->Load("libSTEERBase.so"); | |
11 | gROOT->ProcessLine(".include $ALICE_ROOT/include"); | |
12 | gROOT->LoadMacro("AliCentralityGlauberFit.cxx+"); | |
13 | ||
14 | const char *finnameGlau ="/home/atoia/GlauberNtuple/GlauberMC_PbPb_ntuple_sigma64_mind4_r662_a546.root"; | |
15 | ||
16 | const char * run="167693"; | |
17 | const char * system="TRK"; | |
18 | TString finname = Form("/home/atoia/analysis/data2011/multRef/AnalysisResults_%s.root",run); | |
19 | TString foutname = Form("/home/atoia/analysis/data2011/fit/%s_fit_%s.root",system,run); | |
20 | //TString foutnameGlau = Form("/home/atoia/analysis/data2011/fit/%s_ntuple_%s.root",system,run); | |
21 | TString foutnameGlau = Form("test_%s.root",system,run); | |
22 | //const char *histname="fHOutMultV0M"; | |
23 | //const char *histname="fHOutMultCL1"; | |
24 | const char *histname="fHOutMultTRK"; | |
25 | ||
26 | ||
27 | AliCentralityGlauberFit *mPM = new AliCentralityGlauberFit(finnameGlau); | |
28 | mPM->SetInputFile(finname); | |
29 | mPM->SetInputNtuple(finnameGlau); | |
30 | mPM->SetOutputFile(foutname); | |
31 | mPM->SetOutputNtuple(foutnameGlau); | |
32 | mPM->AddHisto(histname); | |
33 | ||
34 | mPM->SetRebin(1); | |
35 | mPM->SetAncestorMode(2); // If 1 use Npart**alpha, if 2 use alpha*Npart + (1-alpha)*Ncoll | |
36 | mPM->SetFastFitMode(0); // If 1 or 2 use cruder approximation to compute curve faster 1:NBD, 2:Gauss | |
37 | mPM->UseChi2(kTRUE); // If TRUE minimize Chi2 | |
38 | mPM->UseAverage(kFALSE); // If TRUE use Average | |
39 | mPM->SetNtrials(1); | |
40 | mPM->SetNevents(1e5); | |
41 | ||
42 | // ----------range to fit-------------- | |
43 | //mPM->SetRangeToFit(100., 20000.); // V0M | |
44 | //mPM->SetRangeToFit(40., 5200.); // CL1 | |
45 | mPM->SetRangeToFit(10., 2400.); // TRK | |
46 | // ----------range to scale-------------- | |
47 | //mPM->SetRangeToScale(100.); // V0M | |
48 | //mPM->SetRangeToScale(40.); // CL1 | |
49 | mPM->SetRangeToScale(10.); //TRK | |
50 | // ----------initial parameters-------------- | |
51 | //mPM->SetGlauberParam(10,29,31, 10,0.7,1.5, 10,0.85,0.87); | |
52 | //mPM->MakeFitsMinuitNBD(0.801,30.,1.2); // initial parameters | |
53 | //mPM->SetGlauberParam(1,28,28.5, 1,1.291,1.9, 1,0.801,0.82); // V0M | |
54 | //mPM->SetGlauberParam(1,7.25,7.5, 1,1.291,1.9, 1,0.801,0.82); // Cl1 | |
55 | mPM->SetGlauberParam(1,3.5,3.7, 1,1.291,1.5, 1,0.801,0.81); // TRK | |
56 | // ----------done-------------- | |
57 | ||
58 | mPM->MakeFits(); | |
59 | ||
60 | TFile * f = new TFile (foutname); | |
61 | TH1 * hd = (TH1*) gDirectory->Get("fHOutMultTRK"); | |
62 | TH1 * hg = (TH1*) gDirectory->Get("fHOutMultTRK_GLAU"); | |
63 | hg->SetLineColor(kRed); | |
64 | hd->Draw("e"); | |
65 | hg->Draw("same"); | |
66 | ||
67 | ||
68 | ||
69 | } | |
70 |