#include #include #include "Rtypes.h" #include "AtlasUtils.h" // #ifndef __CINT__ #include "AtlasStyle.C" // #endif #include "TCanvas.h" #include "TFile.h" #include "TROOT.h" #include "TH1F.h" #include "TRandom.h" #include "TGraphErrors.h" using namespace std; void prepare() { TFile *F1= new TFile("../test2010_attilio.root","read"); TFile *outfile= new TFile("test_divide.root", "RECREATE"); //histo useful for efficiency //numerators TH1D *PtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4=(TH1D*)F1->Get("PtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4") ; TH1D *PtMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4=(TH1D*)F1->Get("PtMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4") ; TH1D *EtaMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4=(TH1D*)F1->Get("EtaMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4") ; TH1D *EtaMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4=(TH1D*)F1->Get("EtaMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4") ; TH1D *DeltaR_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4=(TH1D*)F1->Get("DeltaR_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4") ; //Barrel //numerators TH1D *WPtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4b=(TH1D*)F1->Get("WPtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4b") ; TH1D *WPtMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4b=(TH1D*)F1->Get("WPtMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4b") ; //denominators TH1D *PtMajor_Jpsi_EFMU4b=(TH1D*)F1->Get("PtMajor_Jpsi_EFMU4b") ; TH1D *PtMinor_Jpsi_EFMU4b=(TH1D*)F1->Get("PtMinor_Jpsi_EFMU4b") ; //EndCap //numerators TH1D *WPtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4e=(TH1D*)F1->Get("WPtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4e") ; TH1D *WPtMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4e=(TH1D*)F1->Get("WPtMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4e") ; //denominators TH1D *PtMajor_Jpsi_EFMU4e=(TH1D*)F1->Get("PtMajor_Jpsi_EFMU4e") ; TH1D *PtMinor_Jpsi_EFMU4e=(TH1D*)F1->Get("PtMinor_Jpsi_EFMU4e") ; //numerators TH1D *WPtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4=(TH1D*)F1->Get("WPtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4") ; TH1D *WPtMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4=(TH1D*)F1->Get("WPtMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4") ; TH1D *WEtaMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4=(TH1D*)F1->Get("WEtaMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4") ; TH1D *WEtaMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4=(TH1D*)F1->Get("WEtaMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4") ; TH1D *WDeltaR_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4=(TH1D*)F1->Get("WDeltaR_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4") ; TH1D *WPtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4_mis=(TH1D*)F1->Get("WPtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4_mis") ; TH1D *WPtMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4_mis=(TH1D*)F1->Get("WPtMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4_mis") ; //denominators TH1D *PtMajor_Jpsi_EFMU4=(TH1D*)F1->Get("PtMajor_Jpsi_EFMU4") ; TH1D *PtMinor_Jpsi_EFMU4=(TH1D*)F1->Get("PtMinor_Jpsi_EFMU4") ; TH1D *EtaMajor_Jpsi_EFMU4=(TH1D*)F1->Get("EtaMajor_Jpsi_EFMU4") ; TH1D *EtaMinor_Jpsi_EFMU4=(TH1D*)F1->Get("EtaMinor_Jpsi_EFMU4") ; TH1D *DeltaR_Jpsi_EFMU4=(TH1D*)F1->Get("DeltaR_Jpsi_EFMU4") ; TH1D *PtMajor_Jpsi_EFMU4_mis=(TH1D*)F1->Get("PtMajor_Jpsi_EFMU4_mis") ; TH1D *PtMinor_Jpsi_EFMU4_mis=(TH1D*)F1->Get("PtMinor_Jpsi_EFMU4_mis") ; TGraphAsymmErrors *Eff_ptMajor_All_noW=new TGraphAsymmErrors(); Eff_ptMajor_All_noW->BayesDivide(PtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4,PtMajor_Jpsi_EFMU4); TGraphAsymmErrors *Eff_ptMajor_All=new TGraphAsymmErrors(); Eff_ptMajor_All->BayesDivide(WPtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4,PtMajor_Jpsi_EFMU4); TGraphAsymmErrors *Eff_ptMajor_Barrel=new TGraphAsymmErrors(); Eff_ptMajor_Barrel->BayesDivide(WPtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4b,PtMajor_Jpsi_EFMU4b); TGraphAsymmErrors *Eff_ptMajor_Endcap=new TGraphAsymmErrors(); Eff_ptMajor_Endcap->BayesDivide(WPtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4e,PtMajor_Jpsi_EFMU4e); //TGraphAsymmErrors *Eff_ptMajor_mis=new TGraphAsymmErrors(); //Eff_ptMajor_mis->BayesDivide(WPtMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4_mis,PtMajor_Jpsi_EFMU4_mis); TGraphAsymmErrors *Eff_ptMinor_All=new TGraphAsymmErrors(); Eff_ptMinor_All->BayesDivide(WPtMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4,PtMinor_Jpsi_EFMU4); TGraphAsymmErrors *Eff_ptMinor_Barrel=new TGraphAsymmErrors(); Eff_ptMinor_Barrel->BayesDivide(WPtMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4b,PtMinor_Jpsi_EFMU4b); TGraphAsymmErrors *Eff_ptMinor_Endcap=new TGraphAsymmErrors(); Eff_ptMinor_Endcap->BayesDivide(WPtMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4e,PtMinor_Jpsi_EFMU4e); //TGraphAsymmErrors *Eff_ptMinor_mis=new TGraphAsymmErrors(); //Eff_ptMinor_mis->BayesDivide(WPtMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4_mis,PtMinor_Jpsi_EFMU4_mis); TGraphAsymmErrors *Eff_etaMajor_All=new TGraphAsymmErrors(); Eff_etaMajor_All->BayesDivide(WEtaMajor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4,EtaMajor_Jpsi_EFMU4); TGraphAsymmErrors *Eff_etaMinor_All=new TGraphAsymmErrors(); Eff_etaMinor_All->BayesDivide(WEtaMinor_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4,EtaMinor_Jpsi_EFMU4); TGraphAsymmErrors *Eff_DeltaR_All=new TGraphAsymmErrors(); Eff_DeltaR_All->BayesDivide(WDeltaR_Jpsi_EFMU4_Jpsimumupass_vs_EFMU4,DeltaR_Jpsi_EFMU4); Eff_ptMajor_All->Write("Eff_ptMajor_All"); //Eff_ptMajor_mis->Write("Eff_ptMajor_mis"); Eff_ptMajor_Barrel->Write("Eff_ptMajor_Barrel"); Eff_ptMajor_Endcap->Write("Eff_ptMajor_Endcap"); Eff_ptMajor_All_noW->Write("Eff_ptMajor_All_noW"); Eff_ptMinor_All->Write("Eff_ptMinor_All"); //Eff_ptMinor_mis->Write("Eff_ptMinor_mis"); Eff_ptMinor_Barrel->Write("Eff_ptMinor_Barrel"); Eff_ptMinor_Endcap->Write("Eff_ptMinor_Endcap"); Eff_etaMajor_All->Write("Eff_etaMajor_All"); Eff_etaMinor_All->Write("Eff_etaMinor_All"); Eff_DeltaR_All->Write("Eff_DeltaR_All"); outfile->Close(); } void plot () { #ifdef __CINT__ gROOT->LoadMacro("AtlasUtils.C"); #endif SetAtlasStyle(); TFile * file0MC_NW = TFile::Open("test_divide.root"); TFile * fileData = TFile::Open("test_divide.root"); TCanvas *c1 = new TCanvas("c1","",10,10,800,600); TH1F *hr = c1->DrawFrame(0.,0.,50.,1.); hr->GetXaxis()->SetTitle("#font[42]{p_{T} [GeV]}"); hr->GetYaxis()->SetTitle("EF_2mu10 Efficiency Major"); fileData->cd(); //Eff_ptMajor_All->SetLineColor(2); Eff_ptMajor_All->SetLineColor(2); Eff_ptMajor_All->SetMarkerColor(2); Eff_ptMajor_All->SetMarkerStyle(24); Eff_ptMajor_All->SetFillStyle(3001); Eff_ptMajor_All->SetFillColor(2); Eff_ptMajor_All->Draw("P"); file0MC_NW->cd(); Eff_ptMajor_All->SetLineColor(4); Eff_ptMajor_All->SetMarkerColor(4); Eff_ptMajor_All->SetMarkerStyle(22); Eff_ptMajor_All->SetFillStyle(3001); Eff_ptMajor_All->SetFillColor(4); Eff_ptMajor_All->Draw("P"); ATLAS_LABEL(0.55, 0.88); //myText( 0.33,0.88, 1,"Preliminary"); //myText(0.56, 0.82,1,"2010 Data #int L dt = 38 pb^{-1}"); myMarkerText(0.7, 0.40,4,22, "MC Zmumu",1.0); myMarkerText(0.7, 0.35,2,24, "MC H150 ",1.0); //myMarkerText(0.7, 0.34,2,24,"Data",1.0); myText(0.7, 0.25,1,"|#eta|<2.4"); //c1->Print("TurnOn_EF_2mu4_Jpsimumu_PtMajor_D_E_data_all.eps"); ////////// TCanvas *c5 = new TCanvas("c5","",10,10,800,600); TH1F *hr = c5->DrawFrame(0.,0.,50.,1.); hr->GetXaxis()->SetTitle("#font[42]{p_{T} [GeV]}"); hr->GetYaxis()->SetTitle("EF_2mu10 Efficiency Minor"); fileData->cd(); Eff_ptMinor_All->SetLineColor(2); Eff_ptMinor_All->SetMarkerColor(2); Eff_ptMinor_All->SetMarkerStyle(24); Eff_ptMinor_All->SetFillStyle(3001); Eff_ptMinor_All->SetFillColor(2); Eff_ptMinor_All->Draw("P"); file0MC_NW->cd(); Eff_ptMinor_All->SetLineColor(4); Eff_ptMinor_All->SetMarkerColor(4); Eff_ptMinor_All->SetMarkerStyle(22); Eff_ptMinor_All->SetFillStyle(3001); Eff_ptMinor_All->SetFillColor(4); Eff_ptMinor_All->Draw("P"); ATLAS_LABEL(0.5, 0.88); myMarkerText(0.7, 0.40,4,22, "MC Zmumu",1.0); myMarkerText(0.7, 0.35,2,24, "MC H150 ",1.0); //myText( 0.33,0.88, 1,"Preliminary"); //myText(0.51, 0.82,1,"2010 Data #int L dt = 38 pb^{-1}"); //myMarkerText(0.7, 0.49,4,22, "MC",1.0); //myMarkerText(0.6, 0.55,3,23,"MC emul : #Delta R ",1.0); //myMarkerText(0.7, 0.44,2,24,"Data",1.0); myText(0.7, 0.25,1,"|#eta_{mu1-mu2}|<2.4"); //c5->Print("TurnOn_EF_2mu4_Jpsimumu_PtMinor_D_E_data_all.eps"); TCanvas *c8 = new TCanvas("c8","",10,10,800,600); TH1F *hr = c8->DrawFrame(-3.,0.,3.,1.); hr->GetXaxis()->SetTitle("#font[42]{#eta}"); hr->GetYaxis()->SetTitle("EF_2mu10 Eta major Efficiency"); fileData->cd(); Eff_etaMajor_All->SetLineColor(2); Eff_etaMajor_All->SetMarkerColor(2); Eff_etaMajor_All->SetMarkerStyle(24); Eff_etaMajor_All->SetFillStyle(3001); Eff_etaMajor_All->SetFillColor(2); Eff_etaMajor_All->Draw("P"); file0MC_NW->cd(); Eff_etaMajor_All->SetLineColor(4); Eff_etaMajor_All->SetMarkerColor(4); Eff_etaMajor_All->SetMarkerStyle(22); Eff_etaMajor_All->SetFillStyle(3001); Eff_etaMajor_All->SetFillColor(4); Eff_etaMajor_All->Draw("P"); ATLAS_LABEL(0.2, 0.88); myMarkerText(0.7, 0.40,4,22, "MC Zmumu",1.0); myMarkerText(0.7, 0.35,2,24, "MC H150 ",1.0); //myText( 0.33,0.88, 1,"Preliminary"); //myText(0.21, 0.82,1,"2010 Data #int L dt = 38 pb^{-1}"); //myMarkerText(0.7, 0.89,4,22, "MC",1.0); //myMarkerText(0.6, 0.55,3,23,"MC emul : #Delta R ",1.0); //myMarkerText(0.7, 0.84,2,24,"Data",1.0); //myText(0.7, 0.75,1,"|#eta_{mu1-mu2}|<2.4"); //c8->Print("TurnOn_EF_2mu4_Jpsimumu_EtaMajor_D_E_data.eps"); TCanvas *c9 = new TCanvas("c9","",10,10,800,600); TH1F *hr = c9->DrawFrame(-3.,0.,3.,1.); hr->GetXaxis()->SetTitle("#font[42]{#eta}"); hr->GetYaxis()->SetTitle("EF_2mu10 eta minor Efficiency"); fileData->cd(); Eff_etaMinor_All->SetLineColor(2); Eff_etaMinor_All->SetMarkerColor(2); Eff_etaMinor_All->SetMarkerStyle(24); Eff_etaMinor_All->SetFillStyle(3001); Eff_etaMinor_All->SetFillColor(2); Eff_etaMinor_All->Draw("P"); file0MC_NW->cd(); Eff_etaMinor_All->SetLineColor(4); Eff_etaMinor_All->SetMarkerColor(4); Eff_etaMinor_All->SetMarkerStyle(22); Eff_etaMinor_All->SetFillStyle(3001); Eff_etaMinor_All->SetFillColor(4); Eff_etaMinor_All->Draw("P"); ATLAS_LABEL(0.2, 0.88); myMarkerText(0.7, 0.40,4,22, "MC Zmumu",1.0); myMarkerText(0.7, 0.35,2,24, "MC H150 ",1.0); //myText( 0.33,0.88, 1,"Preliminary"); //myText(0.21, 0.82,1,"2010 Data #int L dt = 38 pb^{-1}"); //myMarkerText(0.7, 0.89,4,22, "MC",1.0); //myMarkerText(0.6, 0.55,3,23,"MC emul : #Delta R ",1.0); //myMarkerText(0.7, 0.84,2,24,"Data",1.0); //myText(0.7, 0.75,1,"|#eta_{mu1-mu2}|<2.4"); //c9->Print("TurnOn_EF_2mu4_Jpsimumu_EtaMinor_D_E_data.eps"); ////////// TCanvas *c10 = new TCanvas("c10","",10,10,800,600); TH1F *hr = c10->DrawFrame(0.,0.,2.0,1.); hr->GetXaxis()->SetTitle("#font[42]{#Delta R}"); hr->GetYaxis()->SetTitle("EF_2mu10 DR Efficiency"); fileData->cd(); Eff_DeltaR_All->SetLineColor(2); Eff_DeltaR_All->SetMarkerColor(2); Eff_DeltaR_All->SetMarkerStyle(24); Eff_DeltaR_All->SetFillStyle(3001); Eff_DeltaR_All->SetFillColor(2); Eff_DeltaR_All->Draw("P"); file0MC_NW->cd(); Eff_DeltaR_All->SetLineColor(4); Eff_DeltaR_All->SetMarkerColor(4); Eff_DeltaR_All->SetMarkerStyle(22); Eff_DeltaR_All->SetFillStyle(3001); Eff_DeltaR_All->SetFillColor(4); Eff_DeltaR_All->Draw("P"); ATLAS_LABEL(0.2, 0.88); myMarkerText(0.7, 0.40,4,22, "MC Zmumu",1.0); myMarkerText(0.7, 0.35,2,24, "MC H150 ",1.0); //myText( 0.33,0.88, 1,"Preliminary"); //myText(0.21, 0.82,1,"2010 Data #int L dt = 38 pb^{-1}"); //myMarkerText(0.7, 0.89,4,22, "MC",1.0); //myMarkerText(0.6, 0.55,3,23,"MC emul : #Delta R ",1.0); //myMarkerText(0.7, 0.84,2,24,"Data",1.0); //myText(0.7, 0.75,1,"|#eta_{mu1-mu2}|<2.4"); //c10->Print("TurnOn_EF_2mu4_Jpsimumu_DR_D_E_data.eps"); }