PWG2/EVCHAR/makeCentralityFit.C

   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