PWGCF/Correlations/macros/underlyingevent/GetTrackingEfficiency.C

   1 void GetTrackingEfficiency(const char* fileName = "PWG4_JetTasksOutput.root")
   2 {
   3   gSystem->Load("libANALYSIS");
   4   gSystem->Load("libANALYSISalice");
   5   gSystem->Load("libCORRFW");
   6   gSystem->Load("libJETAN");
   7   gSystem->Load("libPWG4JetTasks");
   8
   9   TFile::Open(fileName);
  10   list = (TList*) gFile->Get("PWG4_LeadingTrackUE/histosLeadingTrackUE");
  11
  12   AliUEHistograms* corr = (AliUEHistograms*) list->FindObject("AliUEHistograms");
  13   corr->SetEtaRange(-0.79, 0.79);
  14
  15   obj = (TH1*) corr->GetNumberDensitypT()->GetTrackingEfficiency(1)->Clone("trackingefficiency");
  16
  17   obj->Draw();
  18   obj->Fit("pol2", "", "", 0.5, 1.8);
  19   obj->Fit("pol0", "+", "SAME", 5, 15);
  20
  21   Printf("pol2:");
  22   for (Int_t i=0; i<3; i++)
  23     Printf("par %d: %f", i, obj->GetFunction("pol2")->GetParameter(i));
  24
  25   Printf("pol0:");
  26   for (Int_t i=0; i<2; i++)
  27     Printf("par %d: %f", i, obj->GetFunction("pol0")->GetParameter(i));
  28
  29   // extend up to pT 100
  30   for (Int_t bin=obj->GetXaxis()->FindBin(10); bin <= obj->GetNbinsX(); bin++)
  31     obj->SetBinContent(bin, obj->GetFunction("pol0")->Eval(obj->GetXaxis()->GetBinCenter(bin)));
  32
  33   file = TFile::Open("trackingefficiency.root", "RECREATE");
  34   obj->Write();
  35   file->Close();
  36 }