void Integrated_Eff(){ int oniaMode = 1; int ispPb=1; TFile *f = new TFile("/home/christos/public_html/RpA/eff_pPb/Eff_pPb_1S_eta29_5.root"); TH1D* reco = (TH1D*) f->Get("RecoEventsSumET_HF"); TH1D* gen = (TH1D*) f->Get("GenEventsSumET_HF"); TH1D* eff = (TH1D*) f->Get("EffSumET_HF"); TH1D* mc = (TH1D*) f->Get("Ntracks_MC"); TH1D* data = (TH1D*) f->Get("Ntracks_Data"); //TH1D* intReco = new TH1D("intReco","Integrated Reconstructed Event",1,reco->GetXaxis()->GetXmin(),reco->GetXaxis()->GetXmax()); //TH1D* intGen = new TH1D("intGen","Integrated Gen Event",1,gen->GetXaxis()->GetXmin(),gen->GetXaxis()->GetXmax()); gStyle->SetOptStat(0); TH1F *newreco = (TH1F*)reco->Clone("newreco"); TH1F *newgen = (TH1F*)gen->Clone("newgen"); newgen->SetTitle("#Sigma E_{T}^{HF} Gen & Reco"); TCanvas *c1 = new TCanvas(); newreco->SetLineColor(kGreen+3); newreco->SetLineWidth(2); newgen->SetLineColor(kBlue+2); newgen->SetLineWidth(2); newgen->Draw(); newreco->Draw("SAME"); TLegend *legend1 = new TLegend(0.1,0.7,0.48,0.9); legend1->SetHeader("Gen and Reco #Sigma E_{T}^{HF}"); // option "C" allows to center the header legend1->AddEntry(newreco,"Reco #Sigma E_{T}^{HF}","f"); legend1->AddEntry(newgen,"Gen #Sigma E_{T}^{HF}","f"); legend1->Draw("SAME"); c1->SaveAs(Form("extra_plots/GenAndReco_%dS_%s_HFeta4_5.png",oniaMode,ispPb ? "pPb" : "PP")); //TCanvas *c3 = new TCanvas(); reco->Rebin(reco->GetNbinsX()); gen->Rebin(gen->GetNbinsX()); TGraphAsymmErrors *intEff = new TGraphAsymmErrors(1); intEff->BayesDivide(reco, gen); //TH1D* intEff = new TH1D("intEff","Integrated Efficiency",1,reco->GetXaxis()->GetXmin(),reco->GetXaxis()->GetXmax()); //intEff->Divide(reco,gen,1,1,"B"); TCanvas *c2 = new TCanvas("c2","c2",800,600); eff->SetName("Efficiency"); eff->SetTitle("Efficiency"); //eff->GetXaxis()->SetTitle("#Sigma E_{T}^{HF} (GeV)"); eff->Draw("AP"); intEff->Draw("SAME"); intEff->SetLineStyle(2); intEff->SetLineColor(kRed+3); intEff->SetLineWidth(2); c2->SaveAs(Form("extra_plots/EffWithInt_%dS_%s_HFeta4_5.png",oniaMode,ispPb ? "pPb" : "PP")); TCanvas *c3 = new TCanvas(); data->SetTitle("N_{tracks}"); data->GetXaxis()->SetTitle("N_{tracks}"); data->SetLineColor(kRed+3); data->SetLineWidth(2); mc->SetLineColor(kBlue-2); mc->SetLineWidth(2); //data->Draw(); mc->Draw(); TLegend *legend2 = new TLegend(0.1,0.7,0.48,0.9); legend2->SetHeader("MC and Data N_{tracks}"); // option "C" allows to center the header legend2->AddEntry(mc,"N_{tracks} (MC)","f"); legend2->AddEntry(data,"N_{tracks} (Data)","f"); legend2->Draw("SAME"); c3->SaveAs(Form("extra_plots/MCData_%dS_%s_HFeta4_5.png",oniaMode,ispPb ? "pPb" : "PP")); }