void plotMaker(TH1F *h_axes1, TGraphAsymmErrors *gl, TGraphAsymmErrors *gl0, TGraphAsymmErrors *gl1, TGraphAsymmErrors *gl2, int LOG, double Min, double Max, TString TAG, TString fname, bool ONLYEXPECTED=true ) { TCanvas *c1 = new TCanvas("c","c",800,800); h_axes1->GetYaxis()->SetRangeUser(Min, Max); h_axes1->GetYaxis()->SetTitle("#sigma#times BR(#phi#rightarrow#tau_{lep}#tau_{had}) [pb]"); h_axes1->GetXaxis()->SetTitle("m_{#phi} [GeV]"); h_axes1->Draw("HIST"); gPad->SetLogy(LOG); gl0->SetLineStyle(2); gl1->SetLineStyle(0); gl2->SetLineStyle(0); gl0->SetLineWidth(2); gl1->SetLineWidth(0); gl2->SetLineWidth(0); gl1->SetLineColor(0); gl2->SetLineColor(0); gl1->SetFillColor(kGreen); gl2->SetFillColor(kYellow); gl0->SetLineColor(1); gl0->SetLineWidth(2); gl->SetLineColor(1); gl->SetMarkerColor(1); gl->SetLineWidth(2); gl->SetMarkerStyle(20); gl0->Draw("L,same"); gl2->Draw("E3,same"); gl1->Draw("E3,same"); gl->Draw("PL,same"); gl0->Draw("L,same"); TLegend *leg1 = new TLegend(0.15+0.49, 0.55+0.15, 0.45+0.49, 0.77+0.15); leg1->SetFillColor(kWhite); leg1->SetFillStyle(0); leg1->SetBorderSize(0); if (ONLYEXPECTED) { leg1->AddEntry(gl, "Expected CLs", "LP"); } else { leg1->AddEntry(gl, "Observed CLs", "LP"); leg1->AddEntry(gl0, "Expected CLs", "L"); } leg1->AddEntry(gl1,"1 #sigma", "F"); leg1->AddEntry(gl2,"2 #sigma", "F"); leg1->Draw("same"); TLatex *l = new TLatex; l->SetNDC(); l->SetTextFont(72); l->SetTextColor(1); l->SetTextSize(0.04); //l->DrawLatex(0.15,0.87,"ATLAS Internal, 2012"); TLatex *l1 = new TLatex; l1->SetNDC(); l1->SetTextFont(72); l1->SetTextColor(1); l1->SetTextSize(0.04); //l1->DrawLatex(0.15,0.81,"H#rightarrow #tau_{lep}#tau_{had}, m_{h}^{max}, #mu>0"); TLatex *l2 = new TLatex; l2->SetNDC(); l2->SetTextFont(72); l2->SetTextColor(1); l2->SetTextSize(0.04); //l2->DrawLatex(0.15,0.73,"#int L dt = 20.3 fb^{-1}, #sqrt{s}=8 TeV"); TLatex *l3 = new TLatex; l3->SetNDC(); l3->SetTextFont(72); l3->SetTextColor(1); l3->SetTextSize(0.04); l->DrawLatex(0.19,0.87,"ATLAS Internal, 2012"); // l1->DrawLatex(0.15,0.81,"#phi#rightarrow #tau_{lep}#tau_{had}"+ TAG); l1->DrawLatex(0.19,0.81,TAG); l2->DrawLatex(0.19,0.73,"#int L dt = 20.3 fb^{-1}, #sqrt{s}=8 TeV"); // l3->DrawLatex(0.15,0.65,"low mass selection"); gPad->RedrawAxis(); c1->Print(fname+".png"); //c1->Print(fname+".eps"); } // compare 2 histograms void plotMaker(TH1F *h_axes1, TGraphAsymmErrors *g_comp, TString comp_label, TGraphAsymmErrors *gl, TGraphAsymmErrors *gl0, TGraphAsymmErrors *gl1, TGraphAsymmErrors *gl2, int LOG, double Min, double Max, TString TAG, TString fname, bool ONLYEXPECTED=true ) { TCanvas *c1 = new TCanvas("c","c",800,800); h_axes1->GetYaxis()->SetRangeUser(Min, Max); h_axes1->Draw("HIST"); gPad->SetLogy(LOG); gl0->SetLineStyle(2); gl1->SetLineStyle(0); gl2->SetLineStyle(0); gl0->SetLineWidth(2); gl1->SetLineWidth(0); gl2->SetLineWidth(0); gl1->SetLineColor(0); gl2->SetLineColor(0); gl1->SetFillColor(kGreen); gl2->SetFillColor(kYellow); gl0->SetLineColor(1); gl0->SetLineWidth(2); gl->SetLineColor(1); gl->SetLineStyle(1); gl->SetLineWidth(2); gl->SetMarkerColor(1); gl->SetMarkerStyle(20); gl0->Draw("L,same"); gl2->Draw("E3,same"); gl1->Draw("E3,same"); gl->Draw("PL,same"); gl0->Draw("L,same"); g_comp->SetLineColor(9); g_comp->SetLineStyle(9); g_comp->SetLineWidth(3); g_comp->Draw("L,same"); TLegend *leg1 = new TLegend(0.15+0.49, 0.55+0.15, 0.45+0.49, 0.77+0.15); leg1->SetFillColor(kWhite); leg1->SetFillStyle(0); leg1->SetBorderSize(0); if (ONLYEXPECTED) { leg1->AddEntry(gl, "Expected CLs", "LP"); } else { leg1->AddEntry(gl, "Observed CLs", "LP"); leg1->AddEntry(gl0, "Expected CLs", "L"); } leg1->AddEntry(gl1,"1 #sigma", "F"); leg1->AddEntry(gl2,"2 #sigma", "F"); leg1->AddEntry(g_comp, "Exp "+comp_label, "L"); leg1->Draw("same"); TLatex *l = new TLatex; l->SetNDC(); l->SetTextFont(72); l->SetTextColor(1); l->SetTextSize(0.04); //l->DrawLatex(0.15,0.87,"ATLAS Internal, 2012"); TLatex *l1 = new TLatex; l1->SetNDC(); l1->SetTextFont(72); l1->SetTextColor(1); l1->SetTextSize(0.04); //l1->DrawLatex(0.15,0.81,"H#rightarrow #tau_{lep}#tau_{had}, m_{h}^{max}, #mu>0"); TLatex *l2 = new TLatex; l2->SetNDC(); l2->SetTextFont(72); l2->SetTextColor(1); l2->SetTextSize(0.04); //l2->DrawLatex(0.15,0.73,"#int L dt = 20.3 fb^{-1}, #sqrt{s}=8 TeV"); TLatex *l3 = new TLatex; l3->SetNDC(); l3->SetTextFont(72); l3->SetTextColor(1); l3->SetTextSize(0.04); l->DrawLatex(0.19,0.87,"ATLAS Internal, 2012"); // l1->DrawLatex(0.15,0.81,"#phi#rightarrow #tau_{lep}#tau_{had}"+ TAG); l1->DrawLatex(0.19,0.81,TAG); l2->DrawLatex(0.19,0.73,"#int L dt = 20.3 fb^{-1}, #sqrt{s}=8 TeV"); // l3->DrawLatex(0.15,0.65,"low mass selection"); gPad->RedrawAxis(); c1->Print(fname+".png"); // c1->Print(fname+".eps"); } TGraphAsymmErrors* combineGraphs(TGraphAsymmErrors *low, TGraphAsymmErrors *high, int Nlow, int Nhigh) { TGraphAsymmErrors *h_new = new TGraphAsymmErrors(); for (int i=0; iGetPoint(i, x,y); ey_h = low->GetErrorYhigh(i); ey_l = low->GetErrorYlow(i); ex_h = low->GetErrorXhigh(i); ex_l = low->GetErrorXlow(i); h_new->SetPoint(i, x,y); h_new->SetPointError(i, ex_l,ex_h, ey_l, ey_h); } for (int i=Nhigh; i< high->GetN(); ++i) { double x, y, ex_l, ey_l, ex_h, ey_h; high->GetPoint(i, x,y); ey_h = high->GetErrorYhigh(i); ey_l = high->GetErrorYlow(i); ex_h = high->GetErrorXhigh(i); ex_l = high->GetErrorXlow(i); h_new->SetPoint(i-Nhigh+Nlow, x,y); h_new->SetPointError(i-Nhigh+Nlow, ex_l,ex_h, ey_l, ey_h); } return h_new; } void test() { gROOT->Reset(); gROOT->ProcessLine(".L ~/atlasStyle.C"); gROOT->ProcessLine("setAtlasStyle()"); // this is the conf note result TFile *f = new TFile("./LH_tests.root","read"); // high mass TGraphAsymmErrors *ghc = (TGraphAsymmErrors*) f->Get("COMB_high_pre1/Expected"); TGraphAsymmErrors *ghc0 = (TGraphAsymmErrors*) f->Get("COMB_high_pre1/Expected"); TGraphAsymmErrors *ghc1 = (TGraphAsymmErrors*) f->Get("COMB_high_pre1/Expected1sigma"); TGraphAsymmErrors *ghc2 = (TGraphAsymmErrors*) f->Get("COMB_high_pre1/Expected2sigma"); TGraphAsymmErrors *gh = (TGraphAsymmErrors*) f->Get("LH_high/Expected"); TGraphAsymmErrors *gh0 = (TGraphAsymmErrors*) f->Get("LH_high/Expected"); TGraphAsymmErrors *gh1 = (TGraphAsymmErrors*) f->Get("LH_high/Expected1sigma"); TGraphAsymmErrors *gh2 = (TGraphAsymmErrors*) f->Get("LH_high/Expected2sigma"); TGraphAsymmErrors *ghef = (TGraphAsymmErrors*) f->Get("LH_high_efake/Expected"); TGraphAsymmErrors *ghef0 = (TGraphAsymmErrors*) f->Get("LH_high_efake/Expected"); TGraphAsymmErrors *ghef1 = (TGraphAsymmErrors*) f->Get("LH_high_efake/Expected1sigma"); TGraphAsymmErrors *ghef2 = (TGraphAsymmErrors*) f->Get("LH_high_efake/Expected2sigma"); TGraphAsymmErrors *ghe = (TGraphAsymmErrors*) f->Get("LH_high_el/Expected"); TGraphAsymmErrors *ghe0 = (TGraphAsymmErrors*) f->Get("LH_high_el/Expected"); TGraphAsymmErrors *ghe1 = (TGraphAsymmErrors*) f->Get("LH_high_el/Expected1sigma"); TGraphAsymmErrors *ghe2 = (TGraphAsymmErrors*) f->Get("LH_high_el/Expected2sigma"); TGraphAsymmErrors *ghm = (TGraphAsymmErrors*) f->Get("LH_high_mu/Expected"); TGraphAsymmErrors *ghm0 = (TGraphAsymmErrors*) f->Get("LH_high_mu/Expected"); TGraphAsymmErrors *ghm1 = (TGraphAsymmErrors*) f->Get("LH_high_mu/Expected1sigma"); TGraphAsymmErrors *ghm2 = (TGraphAsymmErrors*) f->Get("LH_high_mu/Expected2sigma"); // TGraphAsymmErrors *ghc = (TGraphAsymmErrors*) f->Get("LH_high_comb/Expected"); // TGraphAsymmErrors *ghc0 = (TGraphAsymmErrors*) f->Get("LH_high_comb/Expected"); // TGraphAsymmErrors *ghc1 = (TGraphAsymmErrors*) f->Get("LH_high_comb/Expected1sigma"); // TGraphAsymmErrors *ghc2 = (TGraphAsymmErrors*) f->Get("LH_high_comb/Expected2sigma"); TGraphAsymmErrors *gv = (TGraphAsymmErrors*) f->Get("LH_veto/Expected"); TGraphAsymmErrors *gv0 = (TGraphAsymmErrors*) f->Get("LH_veto/Expected"); TGraphAsymmErrors *gv1 = (TGraphAsymmErrors*) f->Get("LH_veto/Expected1sigma"); TGraphAsymmErrors *gv2 = (TGraphAsymmErrors*) f->Get("LH_veto/Expected2sigma"); TGraphAsymmErrors *gt = (TGraphAsymmErrors*) f->Get("LH_tag/Expected"); TGraphAsymmErrors *gt0 = (TGraphAsymmErrors*) f->Get("LH_tag/Expected"); TGraphAsymmErrors *gt1 = (TGraphAsymmErrors*) f->Get("LH_tag/Expected1sigma"); TGraphAsymmErrors *gt2 = (TGraphAsymmErrors*) f->Get("LH_tag/Expected2sigma"); TGraphAsymmErrors *gl = (TGraphAsymmErrors*) f->Get("LH_low/Expected"); TGraphAsymmErrors *gl0 = (TGraphAsymmErrors*) f->Get("LH_low/Expected"); TGraphAsymmErrors *gl1 = (TGraphAsymmErrors*) f->Get("LH_low/Expected1sigma"); TGraphAsymmErrors *gl2 = (TGraphAsymmErrors*) f->Get("LH_low/Expected2sigma"); TGraphAsymmErrors *gvef = (TGraphAsymmErrors*) f->Get("LH_veto_efake/Expected"); TGraphAsymmErrors *gvef0 = (TGraphAsymmErrors*) f->Get("LH_veto_efake/Expected"); TGraphAsymmErrors *gvef1 = (TGraphAsymmErrors*) f->Get("LH_veto_efake/Expected1sigma"); TGraphAsymmErrors *gvef2 = (TGraphAsymmErrors*) f->Get("LH_veto_efake/Expected2sigma"); TGraphAsymmErrors *gtef = (TGraphAsymmErrors*) f->Get("LH_tag_efake/Expected"); TGraphAsymmErrors *gtef0 = (TGraphAsymmErrors*) f->Get("LH_tag_efake/Expected"); TGraphAsymmErrors *gtef1 = (TGraphAsymmErrors*) f->Get("LH_tag_efake/Expected1sigma"); TGraphAsymmErrors *gtef2 = (TGraphAsymmErrors*) f->Get("LH_tag_efake/Expected2sigma"); TGraphAsymmErrors *glef = (TGraphAsymmErrors*) f->Get("LH_low_efake/Expected"); TGraphAsymmErrors *glef0 = (TGraphAsymmErrors*) f->Get("LH_low_efake/Expected"); TGraphAsymmErrors *glef1 = (TGraphAsymmErrors*) f->Get("LH_low_efake/Expected1sigma"); TGraphAsymmErrors *glef2 = (TGraphAsymmErrors*) f->Get("LH_low_efake/Expected2sigma"); // June18v1 limits - update of MC samples TGraphAsymmErrors *gh_j18v1 = (TGraphAsymmErrors*) f->Get("LH_high_june18v1/Expected"); TGraphAsymmErrors *gh0_j18v1 = (TGraphAsymmErrors*) f->Get("LH_high_june18v1/Expected"); TGraphAsymmErrors *gh1_j18v1 = (TGraphAsymmErrors*) f->Get("LH_high_june18v1/Expected1sigma"); TGraphAsymmErrors *gh2_j18v1 = (TGraphAsymmErrors*) f->Get("LH_high_june18v1/Expected2sigma"); TGraphAsymmErrors *ghe_j18v1 = (TGraphAsymmErrors*) f->Get("LH_high_el_june18v1/Expected"); TGraphAsymmErrors *ghe0_j18v1 = (TGraphAsymmErrors*) f->Get("LH_high_el_june18v1/Expected"); TGraphAsymmErrors *ghe1_j18v1 = (TGraphAsymmErrors*) f->Get("LH_high_el_june18v1/Expected1sigma"); TGraphAsymmErrors *ghe2_j18v1 = (TGraphAsymmErrors*) f->Get("LH_high_el_june18v1/Expected2sigma"); TGraphAsymmErrors *ghm_j18v1 = (TGraphAsymmErrors*) f->Get("LH_high_mu_june18v1/Expected"); TGraphAsymmErrors *ghm0_j18v1 = (TGraphAsymmErrors*) f->Get("LH_high_mu_june18v1/Expected"); TGraphAsymmErrors *ghm1_j18v1 = (TGraphAsymmErrors*) f->Get("LH_high_mu_june18v1/Expected1sigma"); TGraphAsymmErrors *ghm2_j18v1 = (TGraphAsymmErrors*) f->Get("LH_high_mu_june18v1/Expected2sigma"); TGraphAsymmErrors *gh_j19v1 = (TGraphAsymmErrors*) f->Get("LH_high_June19v1/Expected"); TGraphAsymmErrors *gh0_j19v1 = (TGraphAsymmErrors*) f->Get("LH_high_June19v1/Expected"); TGraphAsymmErrors *gh1_j19v1 = (TGraphAsymmErrors*) f->Get("LH_high_June19v1/Expected1sigma"); TGraphAsymmErrors *gh2_j19v1 = (TGraphAsymmErrors*) f->Get("LH_high_June19v1/Expected2sigma"); // the ntuple has a naming bug, otherwise this is what it is TGraphAsymmErrors *gh_j19v1_bc = (TGraphAsymmErrors*) f->Get("LH_high_june18v1_bc/Expected"); TGraphAsymmErrors *gh0_j19v1_bc = (TGraphAsymmErrors*) f->Get("LH_high_june18v1_bc/Expected"); TGraphAsymmErrors *gh1_j19v1_bc = (TGraphAsymmErrors*) f->Get("LH_high_june18v1_bc/Expected1sigma"); TGraphAsymmErrors *gh2_j19v1_bc = (TGraphAsymmErrors*) f->Get("LH_high_june18v1_bc/Expected2sigma"); // comparison with April 19 limits TFile *f1 = new TFile("../mAtb_Apr19/LH_tests.root","read"); TGraphAsymmErrors *g1h = (TGraphAsymmErrors*) f1->Get("LH_high/Expected"); TGraphAsymmErrors *g1h0 = (TGraphAsymmErrors*) f1->Get("LH_high/Expected"); TGraphAsymmErrors *g1h1 = (TGraphAsymmErrors*) f1->Get("LH_high/Expected1sigma"); TGraphAsymmErrors *g1h2 = (TGraphAsymmErrors*) f1->Get("LH_high/Expected2sigma"); TGraphAsymmErrors *g1v = (TGraphAsymmErrors*) f1->Get("LH_veto/Expected"); TGraphAsymmErrors *g1v0 = (TGraphAsymmErrors*) f1->Get("LH_veto/Expected"); TGraphAsymmErrors *g1v1 = (TGraphAsymmErrors*) f1->Get("LH_veto/Expected1sigma"); TGraphAsymmErrors *g1v2 = (TGraphAsymmErrors*) f1->Get("LH_veto/Expected2sigma"); TGraphAsymmErrors *g1t = (TGraphAsymmErrors*) f1->Get("LH_tag/Expected"); TGraphAsymmErrors *g1t0 = (TGraphAsymmErrors*) f1->Get("LH_tag/Expected"); TGraphAsymmErrors *g1t1 = (TGraphAsymmErrors*) f1->Get("LH_tag/Expected1sigma"); TGraphAsymmErrors *g1t2 = (TGraphAsymmErrors*) f1->Get("LH_tag/Expected2sigma"); TGraphAsymmErrors *g1l = (TGraphAsymmErrors*) f1->Get("LH_low/Expected"); TGraphAsymmErrors *g1l0 = (TGraphAsymmErrors*) f1->Get("LH_low/Expected"); TGraphAsymmErrors *g1l1 = (TGraphAsymmErrors*) f1->Get("LH_low/Expected1sigma"); TGraphAsymmErrors *g1l2 = (TGraphAsymmErrors*) f1->Get("LH_low/Expected2sigma"); // cross section limits TGraphAsymmErrors *xbh = (TGraphAsymmErrors*) f->Get("LH_high_bbA/Expected"); TGraphAsymmErrors *xbh0 = (TGraphAsymmErrors*) f->Get("LH_high_bbA/Expected"); TGraphAsymmErrors *xbh1 = (TGraphAsymmErrors*) f->Get("LH_high_bbA/Expected1sigma"); TGraphAsymmErrors *xbh2 = (TGraphAsymmErrors*) f->Get("LH_high_bbA/Expected2sigma"); TGraphAsymmErrors *xgh = (TGraphAsymmErrors*) f->Get("LH_high_ggA/Expected"); TGraphAsymmErrors *xgh0 = (TGraphAsymmErrors*) f->Get("LH_high_ggA/Expected"); TGraphAsymmErrors *xgh1 = (TGraphAsymmErrors*) f->Get("LH_high_ggA/Expected1sigma"); TGraphAsymmErrors *xgh2 = (TGraphAsymmErrors*) f->Get("LH_high_ggA/Expected2sigma"); TGraphAsymmErrors *xbh_j18v1 = (TGraphAsymmErrors*) f->Get("LH_high_bbA_June18v1/Expected"); TGraphAsymmErrors *xbh0_j18v1 = (TGraphAsymmErrors*) f->Get("LH_high_bbA_June18v1/Expected"); TGraphAsymmErrors *xbh1_j18v1 = (TGraphAsymmErrors*) f->Get("LH_high_bbA_June18v1/Expected1sigma"); TGraphAsymmErrors *xbh2_j18v1 = (TGraphAsymmErrors*) f->Get("LH_high_bbA_June18v1/Expected2sigma"); TGraphAsymmErrors *xgh_j18v1 = (TGraphAsymmErrors*) f->Get("LH_high_ggA_June18v1/Expected"); TGraphAsymmErrors *xgh0_j18v1 = (TGraphAsymmErrors*) f->Get("LH_high_ggA_June18v1/Expected"); TGraphAsymmErrors *xgh1_j18v1 = (TGraphAsymmErrors*) f->Get("LH_high_ggA_June18v1/Expected1sigma"); TGraphAsymmErrors *xgh2_j18v1 = (TGraphAsymmErrors*) f->Get("LH_high_ggA_June18v1/Expected2sigma"); TGraphAsymmErrors *xgv = (TGraphAsymmErrors*) f->Get("LH_veto_ggA/Expected"); TGraphAsymmErrors *xgv0 = (TGraphAsymmErrors*) f->Get("LH_veto_ggA/Expected"); TGraphAsymmErrors *xgv1 = (TGraphAsymmErrors*) f->Get("LH_veto_ggA/Expected1sigma"); TGraphAsymmErrors *xgv2 = (TGraphAsymmErrors*) f->Get("LH_veto_ggA/Expected2sigma"); TGraphAsymmErrors *xbv = (TGraphAsymmErrors*) f->Get("LH_veto_bbA/Expected"); TGraphAsymmErrors *xbv0 = (TGraphAsymmErrors*) f->Get("LH_veto_bbA/Expected"); TGraphAsymmErrors *xbv1 = (TGraphAsymmErrors*) f->Get("LH_veto_bbA/Expected1sigma"); TGraphAsymmErrors *xbv2 = (TGraphAsymmErrors*) f->Get("LH_veto_bbA/Expected2sigma"); TGraphAsymmErrors *xbt = (TGraphAsymmErrors*) f->Get("LH_tag_bbA/Expected"); TGraphAsymmErrors *xbt0 = (TGraphAsymmErrors*) f->Get("LH_tag_bbA/Expected"); TGraphAsymmErrors *xbt1 = (TGraphAsymmErrors*) f->Get("LH_tag_bbA/Expected1sigma"); TGraphAsymmErrors *xbt2 = (TGraphAsymmErrors*) f->Get("LH_tag_bbA/Expected2sigma"); TGraphAsymmErrors *xgt = (TGraphAsymmErrors*) f->Get("LH_tag_ggA/Expected"); TGraphAsymmErrors *xgt0 = (TGraphAsymmErrors*) f->Get("LH_tag_ggA/Expected"); TGraphAsymmErrors *xgt1 = (TGraphAsymmErrors*) f->Get("LH_tag_ggA/Expected1sigma"); TGraphAsymmErrors *xgt2 = (TGraphAsymmErrors*) f->Get("LH_tag_ggA/Expected2sigma"); TGraphAsymmErrors *xbl = (TGraphAsymmErrors*) f->Get("LH_low_bbA/Expected"); TGraphAsymmErrors *xbl0 = (TGraphAsymmErrors*) f->Get("LH_low_bbA/Expected"); TGraphAsymmErrors *xbl1 = (TGraphAsymmErrors*) f->Get("LH_low_bbA/Expected1sigma"); TGraphAsymmErrors *xbl2 = (TGraphAsymmErrors*) f->Get("LH_low_bbA/Expected2sigma"); TGraphAsymmErrors *xgl = (TGraphAsymmErrors*) f->Get("LH_low_ggA/Expected"); TGraphAsymmErrors *xgl0 = (TGraphAsymmErrors*) f->Get("LH_low_ggA/Expected"); TGraphAsymmErrors *xgl1 = (TGraphAsymmErrors*) f->Get("LH_low_ggA/Expected1sigma"); TGraphAsymmErrors *xgl2 = (TGraphAsymmErrors*) f->Get("LH_low_ggA/Expected2sigma"); // ................................................................................... TGraphAsymmErrors *x1bh = (TGraphAsymmErrors*) f1->Get("LH_high_bbA/Expected"); TGraphAsymmErrors *x1bh0 = (TGraphAsymmErrors*) f1->Get("LH_high_bbA/Expected"); TGraphAsymmErrors *x1bh1 = (TGraphAsymmErrors*) f1->Get("LH_high_bbA/Expected1sigma"); TGraphAsymmErrors *x1bh2 = (TGraphAsymmErrors*) f1->Get("LH_high_bbA/Expected2sigma"); TGraphAsymmErrors *x1gh = (TGraphAsymmErrors*) f1->Get("LH_high_ggA/Expected"); TGraphAsymmErrors *x1gh0 = (TGraphAsymmErrors*) f1->Get("LH_high_ggA/Expected"); TGraphAsymmErrors *x1gh1 = (TGraphAsymmErrors*) f1->Get("LH_high_ggA/Expected1sigma"); TGraphAsymmErrors *x1gh2 = (TGraphAsymmErrors*) f1->Get("LH_high_ggA/Expected2sigma"); TGraphAsymmErrors *x1gv = (TGraphAsymmErrors*) f1->Get("LH_veto_ggA/Expected"); TGraphAsymmErrors *x1gv0 = (TGraphAsymmErrors*) f1->Get("LH_veto_ggA/Expected"); TGraphAsymmErrors *x1gv1 = (TGraphAsymmErrors*) f1->Get("LH_veto_ggA/Expected1sigma"); TGraphAsymmErrors *x1gv2 = (TGraphAsymmErrors*) f1->Get("LH_veto_ggA/Expected2sigma"); TGraphAsymmErrors *x1bv = (TGraphAsymmErrors*) f1->Get("LH_veto_bbA/Expected"); TGraphAsymmErrors *x1bv0 = (TGraphAsymmErrors*) f1->Get("LH_veto_bbA/Expected"); TGraphAsymmErrors *x1bv1 = (TGraphAsymmErrors*) f1->Get("LH_veto_bbA/Expected1sigma"); TGraphAsymmErrors *x1bv2 = (TGraphAsymmErrors*) f1->Get("LH_veto_bbA/Expected2sigma"); TGraphAsymmErrors *x1bt = (TGraphAsymmErrors*) f1->Get("LH_tag_bbA/Expected"); TGraphAsymmErrors *x1bt0 = (TGraphAsymmErrors*) f1->Get("LH_tag_bbA/Expected"); TGraphAsymmErrors *x1bt1 = (TGraphAsymmErrors*) f1->Get("LH_tag_bbA/Expected1sigma"); TGraphAsymmErrors *x1bt2 = (TGraphAsymmErrors*) f1->Get("LH_tag_bbA/Expected2sigma"); TGraphAsymmErrors *x1gt = (TGraphAsymmErrors*) f1->Get("LH_tag_ggA/Expected"); TGraphAsymmErrors *x1gt0 = (TGraphAsymmErrors*) f1->Get("LH_tag_ggA/Expected"); TGraphAsymmErrors *x1gt1 = (TGraphAsymmErrors*) f1->Get("LH_tag_ggA/Expected1sigma"); TGraphAsymmErrors *x1gt2 = (TGraphAsymmErrors*) f1->Get("LH_tag_ggA/Expected2sigma"); TGraphAsymmErrors *x1bl = (TGraphAsymmErrors*) f1->Get("LH_low_bbA/Expected"); TGraphAsymmErrors *x1bl0 = (TGraphAsymmErrors*) f1->Get("LH_low_bbA/Expected"); TGraphAsymmErrors *x1bl1 = (TGraphAsymmErrors*) f1->Get("LH_low_bbA/Expected1sigma"); TGraphAsymmErrors *x1bl2 = (TGraphAsymmErrors*) f1->Get("LH_low_bbA/Expected2sigma"); TGraphAsymmErrors *x1gl = (TGraphAsymmErrors*) f1->Get("LH_low_ggA/Expected"); TGraphAsymmErrors *x1gl0 = (TGraphAsymmErrors*) f1->Get("LH_low_ggA/Expected"); TGraphAsymmErrors *x1gl1 = (TGraphAsymmErrors*) f1->Get("LH_low_ggA/Expected1sigma"); TGraphAsymmErrors *x1gl2 = (TGraphAsymmErrors*) f1->Get("LH_low_ggA/Expected2sigma"); // 2011 data limit combined TFile *f1 = new TFile("../limitgraphs_2011paper.root","read"); TGraphAsymmErrors *g0_all2011 = (TGraphAsymmErrors*) f1->Get("All_combined/Expected"); TGraphAsymmErrors *g0_lh2011 = (TGraphAsymmErrors*) f1->Get("LH_combined/Expected"); // CMS Limit 20 fb-1 TFile *f_cms = new TFile("../CMS_HCP_combined.root","read"); TGraphAsymmErrors *g0_cms = (TGraphAsymmErrors*) f_cms->Get("CMS_HCP_combined/Expected"); // LEP limit - from Julian update TCutG* cutg_lep = new TCutG("cutg_lep", 38+2); cutg_lep->SetVarX(""); cutg_lep->SetVarY(""); cutg_lep->SetTitle("H_LEP"); cutg_lep->SetLineColor(1); cutg_lep->SetLineWidth(2); cutg_lep->SetMarkerStyle(20); cutg_lep->SetMarkerSize(1.2); // limit 95% from Julian cutg_lep->SetPoint(0,93.4,60.); cutg_lep->SetPoint(0,93.4,60.); cutg_lep->SetPoint(1,93.4,34.2712); cutg_lep->SetPoint(2,93.4,13.3771); cutg_lep->SetPoint(3,94.6728,9.32892); cutg_lep->SetPoint(4,96.4992,7.37376); cutg_lep->SetPoint(5,98.6301,6.17662); cutg_lep->SetPoint(6,102.588,5.95432); cutg_lep->SetPoint(7,104.718,5.91806); cutg_lep->SetPoint(8,104.11,6.21447); cutg_lep->SetPoint(9,105.327,6.93679); cutg_lep->SetPoint(10,108.98,7.67245); cutg_lep->SetPoint(11,112.329,8.33201); cutg_lep->SetPoint(12,115.068,8.43444); cutg_lep->SetPoint(13,116.895,7.86224); cutg_lep->SetPoint(14,124.505,6.06445); cutg_lep->SetPoint(15,135.464,4.76427); cutg_lep->SetPoint(16,148.25,3.91831); cutg_lep->SetPoint(17,166.819,3.37364); cutg_lep->SetPoint(18,188.432,3.04087); cutg_lep->SetPoint(19,236.834,2.60222); cutg_lep->SetPoint(20,286.454,2.41091); cutg_lep->SetPoint(21,373.516,2.26113); cutg_lep->SetPoint(22,452.664,2.19311); cutg_lep->SetPoint(23,536.986,2.19311); cutg_lep->SetPoint(24,589.954,2.14671); // - just extrapolated - to be checked cutg_lep->SetPoint(25,900,2.14); cutg_lep->SetPoint(26,900,0.68); // lower hole cutg_lep->SetPoint(25+2,585.997,0.680772); cutg_lep->SetPoint(26+2,545.205,0.676626); cutg_lep->SetPoint(27+2,499.543,0.680772); cutg_lep->SetPoint(28+2,394.825,0.654269); cutg_lep->SetPoint(29+2,339.422,0.642388); cutg_lep->SetPoint(30+2,308.98,0.609883); cutg_lep->SetPoint(31+2,289.193,0.584353); cutg_lep->SetPoint(32+2,273.668,0.54471); cutg_lep->SetPoint(33+2,262.709,0.513998); cutg_lep->SetPoint(34+2,254.49,0.471866); cutg_lep->SetPoint(35+2,253.881,0.420156); cutg_lep->SetPoint(36+2,253.881,0.392851); // closing the area cutg_lep->SetPoint(37+2,253.881,-1); cutg_lep->SetPoint(38+2,0,-1); TH1F *h_axes = new TH1F("","", 100, 90, 800); h_axes->Fill(-10, 1000); h_axes->GetYaxis()->SetTitle("tan#beta"); h_axes->GetXaxis()->SetTitle("m_{A} [GeV/c^{2}]"); h_axes->GetYaxis()->SetRangeUser(0,60); h_axes->SetLineColor(0); TH1F *h_axes1 = new TH1F("","", 100, 90, 350); h_axes1->Fill(-10, 1000); h_axes1->GetYaxis()->SetTitle("tan#beta"); h_axes1->GetXaxis()->SetTitle("m_{A} [GeV/c^{2}]"); h_axes1->GetYaxis()->SetRangeUser(0,60); h_axes1->SetLineColor(0); TLatex *l = new TLatex; l->SetNDC(); l->SetTextFont(72); l->SetTextColor(1); l->SetTextSize(0.04); //l->DrawLatex(0.15,0.87,"ATLAS Internal, 2012"); TLatex *l1 = new TLatex; l1->SetNDC(); l1->SetTextFont(72); l1->SetTextColor(1); l1->SetTextSize(0.04); //l1->DrawLatex(0.15,0.81,"H#rightarrow #tau_{lep}#tau_{had}, m_{h}^{max}, #mu>0"); TLatex *l2 = new TLatex; l2->SetNDC(); l2->SetTextFont(72); l2->SetTextColor(1); l2->SetTextSize(0.04); //l2->DrawLatex(0.15,0.73,"#int L dt = 20.3 fb^{-1}, #sqrt{s}=8 TeV"); TLatex *l3 = new TLatex; l3->SetNDC(); l3->SetTextFont(72); l3->SetTextColor(1); l3->SetTextSize(0.04); cutg_lep->SetFillColor(1); cutg_lep->SetFillStyle(3004); // Plot #1: high mass optimization // ============================== TCanvas *c1 = new TCanvas("c","c",800,800); h_axes->Draw("HIST"); gh0->SetLineStyle(2); gh1->SetLineStyle(0); gh2->SetLineStyle(0); gh0->SetLineWidth(2); gh1->SetLineWidth(0); gh2->SetLineWidth(0); gh1->SetLineColor(0); gh2->SetLineColor(0); gh1->SetFillColor(kGreen); gh2->SetFillColor(kYellow); gh0->SetLineColor(1); gh0->SetLineWidth(2); gh->SetLineColor(1); gh->SetMarkerColor(1); gh->SetMarkerStyle(20); gh0->Draw("L,same"); gh2->Draw("E3,same"); gh1->Draw("E3,same"); gh->Draw("PL,same"); gh0->Draw("L,same"); g1h0->SetLineColor(2); g1h0->SetLineStyle(9); g1h0->SetLineWidth(2); g1h0->Draw("L,same"); TLegend *leg1 = new TLegend(0.15+0.47, 0.55-0.3, 0.45+0.47, 0.77-0.3); leg1->SetFillColor(kWhite); leg1->SetFillStyle(0); leg1->SetBorderSize(0); // leg1->AddEntry(gl, "Observed CLs", "LP"); // leg1->AddEntry(gl0, "Expected CLs", "L"); leg1->AddEntry(gh, "Expected CLs", "LP"); leg1->AddEntry(gh1,"1 #sigma", "F"); leg1->AddEntry(gh2,"2 #sigma", "F"); leg1->AddEntry(g1h0, "Apr19", "L"); leg1->AddEntry(cutg_lep,"LEP exclusion", "F"); leg1->Draw("same"); cutg_lep->Draw("fsame"); cutg_lep->Draw("lsame"); l->DrawLatex(0.15,0.87,"ATLAS Internal, 2012"); l1->DrawLatex(0.15,0.81,"H#rightarrow #tau_{lep}#tau_{had}, m_{h}^{max}, #mu>0"); l2->DrawLatex(0.15,0.73,"#int L dt = 20.3 fb^{-1}, #sqrt{s}=8 TeV"); l3->DrawLatex(0.15,0.65,"high mass selection, Aug 3"); gPad->RedrawAxis(); c1->Print("limit_lh_2d_highmass.png"); //c1->Print("limit_lh_2d_hihgmass.eps"); TCanvas *c1 = new TCanvas("c","c",800,800); h_axes->Draw("HIST"); ghc0->SetLineStyle(2); ghc1->SetLineStyle(0); ghc2->SetLineStyle(0); ghc0->SetLineWidth(2); ghc1->SetLineWidth(0); ghc2->SetLineWidth(0); ghc1->SetLineColor(0); ghc2->SetLineColor(0); ghc1->SetFillColor(kGreen); ghc2->SetFillColor(kYellow); ghc0->SetLineColor(1); ghc0->SetLineWidth(2); ghc->SetLineColor(1); ghc->SetMarkerColor(1); ghc->SetMarkerStyle(20); ghc0->Draw("L,same"); ghc2->Draw("E3,same"); ghc1->Draw("E3,same"); ghc->Draw("PL,same"); ghc0->Draw("L,same"); gh0->SetLineColor(2); gh0->SetLineStyle(9); gh0->SetLineWidth(2); gh0->Draw("L,same"); gl0->SetLineColor(4); gl0->SetLineStyle(9); gl0->SetLineWidth(2); gl0->Draw("L,same"); g0_cms->SetLineColor(1); g0_cms->SetLineWidth(3); g0_cms->Draw("L,same"); g0_lh2011->SetLineColor(1); g0_lh2011->SetLineWidth(3); g0_lh2011->SetLineStyle(3); g0_lh2011->Draw("L,same"); g0_all2011->SetLineColor(4); g0_all2011->SetLineWidth(3); g0_all2011->SetLineStyle(4); g0_all2011->Draw("L,same"); TLegend *leg1 = new TLegend(0.15+0.52, 0.55-0.34, 0.45+0.52, 0.77-0.34); leg1->SetFillColor(kWhite); leg1->SetFillStyle(0); leg1->SetBorderSize(0); // leg1->AddEntry(gl, "Observed CLs", "LP"); // leg1->AddEntry(gl0, "Expected CLs", "L"); leg1->AddEntry(ghc, "Expected CLs", "LP"); leg1->AddEntry(ghc1,"1 #sigma", "F"); leg1->AddEntry(ghc2,"2 #sigma", "F"); leg1->AddEntry(gh0, "lephad, high", "L"); leg1->AddEntry(gl0, "lephad, low", "L"); leg1->AddEntry(g0_lh2011, "ATLAS lephad, 2011", "L"); leg1->AddEntry(g0_all2011, "ATLAS comb, 2011", "L"); leg1->AddEntry(g0_cms, "CMS HCP", "L"); leg1->AddEntry(cutg_lep,"LEP exclusion", "F"); leg1->Draw("same"); cutg_lep->Draw("fsame"); cutg_lep->Draw("lsame"); l->DrawLatex(0.15,0.87,"ATLAS Internal, 2012"); l1->DrawLatex(0.15,0.81,"H#rightarrow #tau_{lep}#tau_{had}, m_{h}^{max}, #mu>0"); l2->DrawLatex(0.15,0.73,"#int L dt = 19.5 - 20.3 fb^{-1}, #sqrt{s}=8 TeV"); l3->DrawLatex(0.15,0.65,"high mass combination, Aug 3"); gPad->RedrawAxis(); c1->Print("limit_comb_2d_highmass.png"); // high mass combination and comparison with all the rest TCanvas *c1 = new TCanvas("c","c",800,800); h_axes->Draw("HIST"); gh0->SetLineStyle(2); gh1->SetLineStyle(0); gh2->SetLineStyle(0); gh0->SetLineWidth(2); gh1->SetLineWidth(0); gh2->SetLineWidth(0); gh1->SetLineColor(0); gh2->SetLineColor(0); gh1->SetFillColor(kGreen); gh2->SetFillColor(kYellow); gh0->SetLineColor(1); gh0->SetLineWidth(2); gh->SetLineColor(1); gh->SetMarkerColor(1); gh->SetMarkerStyle(20); gh0->Draw("L,same"); gh2->Draw("E3,same"); gh1->Draw("E3,same"); gh->Draw("PL,same"); gh0->Draw("L,same"); ghef0->SetLineColor(2); ghef0->SetLineStyle(9); ghef0->SetLineWidth(2); ghef0->Draw("L,same"); TLegend *leg1 = new TLegend(0.15+0.47, 0.55-0.3, 0.45+0.47, 0.77-0.3); leg1->SetFillColor(kWhite); leg1->SetFillStyle(0); leg1->SetBorderSize(0); // leg1->AddEntry(gl, "Observed CLs", "LP"); // leg1->AddEntry(gl0, "Expected CLs", "L"); leg1->AddEntry(gh, "Expected CLs", "LP"); leg1->AddEntry(gh1,"1 #sigma", "F"); leg1->AddEntry(gh2,"2 #sigma", "F"); leg1->AddEntry(g1h0, "with EFAKE", "L"); leg1->AddEntry(cutg_lep,"LEP exclusion", "F"); leg1->Draw("same"); cutg_lep->Draw("fsame"); cutg_lep->Draw("lsame"); l->DrawLatex(0.15,0.87,"ATLAS Internal, 2012"); l1->DrawLatex(0.15,0.81,"H#rightarrow #tau_{lep}#tau_{had}, m_{h}^{max}, #mu>0"); l2->DrawLatex(0.15,0.73,"#int L dt = 20.3 fb^{-1}, #sqrt{s}=8 TeV"); l3->DrawLatex(0.15,0.65,"high mass selection, Aug 3"); gPad->RedrawAxis(); c1->Print("limit_lh_2d_highmass_2.png"); // low, tag and veto comparison // veto TCanvas *c2 = new TCanvas("c2","c2",800,800); h_axes1->Draw("HIST"); gv0->SetLineStyle(2); gv1->SetLineStyle(0); gv2->SetLineStyle(0); gv0->SetLineWidth(2); gv1->SetLineWidth(0); gv2->SetLineWidth(0); gv1->SetLineColor(0); gv2->SetLineColor(0); gv1->SetFillColor(kGreen); gv2->SetFillColor(kYellow); gv0->SetLineColor(1); gv0->SetLineWidth(2); gv->SetLineColor(1); gv->SetMarkerColor(1); gv->SetMarkerStyle(20); gv0->Draw("L,same"); gv2->Draw("E3,same"); gv1->Draw("E3,same"); gv->Draw("PL,same"); gv0->Draw("L,same"); g1v0->SetLineColor(2); g1v0->SetLineStyle(9); g1v0->SetLineWidth(2); g1v0->Draw("L,same"); TLegend *leg1 = new TLegend(0.15+0.47, 0.55+0.15, 0.45+0.47, 0.77+0.15); leg1->SetFillColor(kWhite); leg1->SetFillStyle(0); leg1->SetBorderSize(0); // leg1->AddEntry(gl, "Observed CLs", "LP"); // leg1->AddEntry(gl0, "Expected CLs", "L"); leg1->AddEntry(gv, "Expected CLs", "LP"); leg1->AddEntry(gv1,"1 #sigma", "F"); leg1->AddEntry(gv2,"2 #sigma", "F"); leg1->AddEntry(g1v0, "April 19", "L"); leg1->AddEntry(cutg_lep,"LEP exclusion", "F"); leg1->Draw("same"); cutg_lep->Draw("fsame"); cutg_lep->Draw("lsame"); l->DrawLatex(0.15,0.87,"ATLAS Internal, 2012"); l1->DrawLatex(0.15,0.81,"H#rightarrow #tau_{lep}#tau_{had}, m_{h}^{max}, #mu>0"); l2->DrawLatex(0.15,0.73,"#int L dt = 20.3 fb^{-1}, #sqrt{s}=8 TeV"); l3->DrawLatex(0.15,0.65,"veto category, Aug 3"); gPad->RedrawAxis(); c2->Print("limit_lh_2d_veto.png"); // tag // veto TCanvas *c2 = new TCanvas("c2","c2",800,800); h_axes1->Draw("HIST"); gt0->SetLineStyle(2); gt1->SetLineStyle(0); gt2->SetLineStyle(0); gt0->SetLineWidth(2); gt1->SetLineWidth(0); gt2->SetLineWidth(0); gt1->SetLineColor(0); gt2->SetLineColor(0); gt1->SetFillColor(kGreen); gt2->SetFillColor(kYellow); gt0->SetLineColor(1); gt0->SetLineWidth(2); gt->SetLineColor(1); gt->SetMarkerColor(1); gt->SetMarkerStyle(20); gt0->Draw("L,same"); gt2->Draw("E3,same"); gt1->Draw("E3,same"); gt->Draw("PL,same"); gt0->Draw("L,same"); g1t0->SetLineColor(2); g1t0->SetLineStyle(9); g1t0->SetLineWidth(2); g1t0->Draw("L,same"); TLegend *leg1 = new TLegend(0.15+0.47, 0.55+0.15, 0.45+0.47, 0.77+0.15); leg1->SetFillColor(kWhite); leg1->SetFillStyle(0); leg1->SetBorderSize(0); // leg1->AddEntry(gl, "Observed CLs", "LP"); // leg1->AddEntry(gl0, "Expected CLs", "L"); leg1->AddEntry(gt, "Expected CLs", "LP"); leg1->AddEntry(gt1,"1 #sigma", "F"); leg1->AddEntry(gt2,"2 #sigma", "F"); leg1->AddEntry(g1t0, "April 19", "L"); leg1->AddEntry(cutg_lep,"LEP exclusion", "F"); leg1->Draw("same"); cutg_lep->Draw("fsame"); cutg_lep->Draw("lsame"); l->DrawLatex(0.15,0.87,"ATLAS Internal, 2012"); l1->DrawLatex(0.15,0.81,"H#rightarrow #tau_{lep}#tau_{had}, m_{h}^{max}, #mu>0"); l2->DrawLatex(0.15,0.73,"#int L dt = 20.3 fb^{-1}, #sqrt{s}=8 TeV"); l3->DrawLatex(0.15,0.65,"tag category, Aug 3"); gPad->RedrawAxis(); c2->Print("limit_lh_2d_tag.png"); // low category TCanvas *c2 = new TCanvas("c2","c2",800,800); h_axes1->Draw("HIST"); gl0->SetLineStyle(2); gl1->SetLineStyle(0); gl2->SetLineStyle(0); gl0->SetLineWidth(2); gl1->SetLineWidth(0); gl2->SetLineWidth(0); gl1->SetLineColor(0); gl2->SetLineColor(0); gl1->SetFillColor(kGreen); gl2->SetFillColor(kYellow); gl0->SetLineColor(1); gl0->SetLineWidth(2); gl->SetLineColor(1); gl->SetMarkerColor(1); gl->SetMarkerStyle(20); gl0->Draw("L,same"); gl2->Draw("E3,same"); gl1->Draw("E3,same"); gl->Draw("PL,same"); gl0->Draw("L,same"); g1l0->SetLineColor(2); g1l0->SetLineStyle(9); g1l0->SetLineWidth(2); g1l0->Draw("L,same"); TLegend *leg1 = new TLegend(0.15+0.47, 0.55+0.15, 0.45+0.47, 0.77+0.15); leg1->SetFillColor(kWhite); leg1->SetFillStyle(0); leg1->SetBorderSize(0); // leg1->AddEntry(gl, "Observed CLs", "LP"); // leg1->AddEntry(gl0, "Expected CLs", "L"); leg1->AddEntry(gl, "Expected CLs", "LP"); leg1->AddEntry(gl1,"1 #sigma", "F"); leg1->AddEntry(gl2,"2 #sigma", "F"); leg1->AddEntry(g1l0, "April 19", "L"); leg1->AddEntry(cutg_lep,"LEP exclusion", "F"); leg1->Draw("same"); cutg_lep->Draw("fsame"); cutg_lep->Draw("lsame"); l->DrawLatex(0.15,0.87,"ATLAS Internal, 2012"); l1->DrawLatex(0.15,0.81,"H#rightarrow #tau_{lep}#tau_{had}, m_{h}^{max}, #mu>0"); l2->DrawLatex(0.15,0.73,"#int L dt = 20.3 fb^{-1}, #sqrt{s}=8 TeV"); l3->DrawLatex(0.15,0.65,"low, Aug 3"); gPad->RedrawAxis(); c2->Print("limit_lh_2d_low.png"); // low tag and veto, but with efake added on the plot to show the effect TCanvas *c2 = new TCanvas("c2","c2",800,800); h_axes1->Draw("HIST"); gv0->SetLineStyle(2); gv1->SetLineStyle(0); gv2->SetLineStyle(0); gv0->SetLineWidth(2); gv1->SetLineWidth(0); gv2->SetLineWidth(0); gv1->SetLineColor(0); gv2->SetLineColor(0); gv1->SetFillColor(kGreen); gv2->SetFillColor(kYellow); gv0->SetLineColor(1); gv0->SetLineWidth(2); gv->SetLineColor(1); gv->SetMarkerColor(1); gv->SetMarkerStyle(20); gv0->Draw("L,same"); gv2->Draw("E3,same"); gv1->Draw("E3,same"); gv->Draw("PL,same"); gv0->Draw("L,same"); g1v0->SetLineColor(2); g1v0->SetLineStyle(9); g1v0->SetLineWidth(2); g1v0->Draw("L,same"); gvef0->SetLineColor(4); gvef0->SetLineStyle(9); gvef0->SetLineWidth(3); gvef0->Draw("L,same"); TLegend *leg1 = new TLegend(0.15+0.47, 0.55+0.15, 0.45+0.47, 0.77+0.15); leg1->SetFillColor(kWhite); leg1->SetFillStyle(0); leg1->SetBorderSize(0); // leg1->AddEntry(gl, "Observed CLs", "LP"); // leg1->AddEntry(gl0, "Expected CLs", "L"); leg1->AddEntry(gv, "Expected CLs", "LP"); leg1->AddEntry(gv1,"1 #sigma", "F"); leg1->AddEntry(gv2,"2 #sigma", "F"); leg1->AddEntry(g1v0, "April 19", "L"); leg1->AddEntry(gvef0, "efake", "L"); leg1->AddEntry(cutg_lep,"LEP exclusion", "F"); leg1->Draw("same"); cutg_lep->Draw("fsame"); cutg_lep->Draw("lsame"); l->DrawLatex(0.15,0.87,"ATLAS Internal, 2012"); l1->DrawLatex(0.15,0.81,"H#rightarrow #tau_{lep}#tau_{had}, m_{h}^{max}, #mu>0"); l2->DrawLatex(0.15,0.73,"#int L dt = 20.3 fb^{-1}, #sqrt{s}=8 TeV"); l3->DrawLatex(0.15,0.65,"veto category, Aug 3"); gPad->RedrawAxis(); c2->Print("limit_lh_2d_veto_2.png"); TCanvas *c2 = new TCanvas("c2","c2",800,800); h_axes1->Draw("HIST"); gt0->SetLineStyle(2); gt1->SetLineStyle(0); gt2->SetLineStyle(0); gt0->SetLineWidth(2); gt1->SetLineWidth(0); gt2->SetLineWidth(0); gt1->SetLineColor(0); gt2->SetLineColor(0); gt1->SetFillColor(kGreen); gt2->SetFillColor(kYellow); gt0->SetLineColor(1); gt0->SetLineWidth(2); gt->SetLineColor(1); gt->SetMarkerColor(1); gt->SetMarkerStyle(20); gt0->Draw("L,same"); gt2->Draw("E3,same"); gt1->Draw("E3,same"); gt->Draw("PL,same"); gt0->Draw("L,same"); g1t0->SetLineColor(2); g1t0->SetLineStyle(9); g1t0->SetLineWidth(2); g1t0->Draw("L,same"); gtef0->SetLineColor(4); gtef0->SetLineStyle(9); gtef0->SetLineWidth(3); gtef0->Draw("L,same"); TLegend *leg1 = new TLegend(0.15+0.47, 0.55+0.15, 0.45+0.47, 0.77+0.15); leg1->SetFillColor(kWhite); leg1->SetFillStyle(0); leg1->SetBorderSize(0); // leg1->AddEntry(gl, "Observed CLs", "LP"); // leg1->AddEntry(gl0, "Expected CLs", "L"); leg1->AddEntry(gt, "Expected CLs", "LP"); leg1->AddEntry(gt1,"1 #sigma", "F"); leg1->AddEntry(gt2,"2 #sigma", "F"); leg1->AddEntry(g1t0, "April 19", "L"); leg1->AddEntry(gtef0, "efake", "L"); leg1->AddEntry(cutg_lep,"LEP exclusion", "F"); leg1->Draw("same"); cutg_lep->Draw("fsame"); cutg_lep->Draw("lsame"); l->DrawLatex(0.15,0.87,"ATLAS Internal, 2012"); l1->DrawLatex(0.15,0.81,"H#rightarrow #tau_{lep}#tau_{had}, m_{h}^{max}, #mu>0"); l2->DrawLatex(0.15,0.73,"#int L dt = 20.3 fb^{-1}, #sqrt{s}=8 TeV"); l3->DrawLatex(0.15,0.65,"tag category, Aug 3"); gPad->RedrawAxis(); c2->Print("limit_lh_2d_tag_2.png"); TCanvas *c2 = new TCanvas("c2","c2",800,800); h_axes1->Draw("HIST"); gl0->SetLineStyle(2); gl1->SetLineStyle(0); gl2->SetLineStyle(0); gl0->SetLineWidth(2); gl1->SetLineWidth(0); gl2->SetLineWidth(0); gl1->SetLineColor(0); gl2->SetLineColor(0); gl1->SetFillColor(kGreen); gl2->SetFillColor(kYellow); gl0->SetLineColor(1); gl0->SetLineWidth(2); gl->SetLineColor(1); gl->SetMarkerColor(1); gl->SetMarkerStyle(20); gl0->Draw("L,same"); gl2->Draw("E3,same"); gl1->Draw("E3,same"); gl->Draw("PL,same"); gl0->Draw("L,same"); g1l0->SetLineColor(2); g1l0->SetLineStyle(9); g1l0->SetLineWidth(2); g1l0->Draw("L,same"); glef0->SetLineColor(4); glef0->SetLineStyle(9); glef0->SetLineWidth(3); glef0->Draw("L,same"); TLegend *leg1 = new TLegend(0.15+0.47, 0.55+0.15, 0.45+0.47, 0.77+0.15); leg1->SetFillColor(kWhite); leg1->SetFillStyle(0); leg1->SetBorderSize(0); // leg1->AddEntry(gl, "Observed CLs", "LP"); // leg1->AddEntry(gl0, "Expected CLs", "L"); leg1->AddEntry(gl, "Expected CLs", "LP"); leg1->AddEntry(gl1,"1 #sigma", "F"); leg1->AddEntry(gl2,"2 #sigma", "F"); leg1->AddEntry(g1l0, "April 19", "L"); leg1->AddEntry(glef0, "efake", "L"); leg1->AddEntry(cutg_lep,"LEP exclusion", "F"); leg1->Draw("same"); cutg_lep->Draw("fsame"); cutg_lep->Draw("lsame"); l->DrawLatex(0.15,0.87,"ATLAS Internal, 2012"); l1->DrawLatex(0.15,0.81,"H#rightarrow #tau_{lep}#tau_{had}, m_{h}^{max}, #mu>0"); l2->DrawLatex(0.15,0.73,"#int L dt = 20.3 fb^{-1}, #sqrt{s}=8 TeV"); l3->DrawLatex(0.15,0.65,"low, Aug 3"); gPad->RedrawAxis(); c2->Print("limit_lh_2d_low_2.png"); plotMaker(h_axes1, xbt, xbt0, xbt1, xbt2, 1, 0.1, 100., "bb#phi#rightarrow #tau_{lep}#tau_{had} tag", "limit_lh_xs_tag_bbA", true); plotMaker(h_axes1, xbv, xbv0, xbv1, xbv2, 1, 0.1, 100., "bb#phi#rightarrow #tau_{lep}#tau_{had} veto", "limit_lh_xs_veto_bbA", true); plotMaker(h_axes1, xgt, xgt0, xgt1, xgt2, 1, 0.1, 100., "gg#phi#rightarrow #tau_{lep}#tau_{had} tag", "limit_lh_xs_tag_ggA", true); plotMaker(h_axes1, xgv, xgv0, xgv1, xgv2, 1, 0.1, 100., "gg#phi#rightarrow #tau_{lep}#tau_{had} veto", "limit_lh_xs_veto_ggA", true); return; plotMaker(h_axes, xbh, xbh0, xbh1, xbh2, 1, 0.01, 100., "bb#phi#rightarrow #tau_{lep}#tau_{had} high mass", "limit_lh_xs_highmass_bbA", true); plotMaker(h_axes, xgh, xgh0, xgh1, xgh2, 1, 0.01, 100., "gg#phi#rightarrow #tau_{lep}#tau_{had} high mass", "limit_lh_xs_highmass_ggA", true); return; // high mass optimization vs June18v1 TCanvas *c11 = new TCanvas("c11","c11",800,800); h_axes->Draw("HIST"); gh0_j18v1->SetLineStyle(2); gh1_j18v1->SetLineStyle(0); gh2_j18v1->SetLineStyle(0); gh0_j18v1->SetLineWidth(2); gh1_j18v1->SetLineWidth(0); gh2_j18v1->SetLineWidth(0); gh1_j18v1->SetLineColor(0); gh2_j18v1->SetLineColor(0); gh1_j18v1->SetFillColor(kGreen); gh2_j18v1->SetFillColor(kYellow); gh0_j18v1->SetLineColor(1); gh0_j18v1->SetLineWidth(2); gh_j18v1->SetLineColor(1); gh_j18v1->SetMarkerColor(1); gh_j18v1->SetMarkerStyle(20); gh0_j18v1->Draw("L,same"); gh2_j18v1->Draw("E3,same"); gh1_j18v1->Draw("E3,same"); gh_j18v1->Draw("PL,same"); gh0_j18v1->Draw("L,same"); g1h0->SetLineColor(2); g1h0->SetLineStyle(9); g1h0->SetLineWidth(2); g1h0->Draw("L,same"); gh0->SetLineColor(4); gh0->SetLineStyle(2); gh0->SetLineWidth(3); gh0->Draw("L,same"); TLegend *leg1 = new TLegend(0.15+0.47, 0.55-0.3, 0.45+0.47, 0.77-0.3); leg1->SetFillColor(kWhite); leg1->SetFillStyle(0); leg1->SetBorderSize(0); // leg1->AddEntry(gl, "Observed CLs", "LP"); // leg1->AddEntry(gl0, "Expected CLs", "L"); leg1->AddEntry(gh_j18v1, "Expected CLs", "LP"); leg1->AddEntry(gh1_j18v1,"1 #sigma", "F"); leg1->AddEntry(gh2_j18v1,"2 #sigma", "F"); leg1->AddEntry(g1h0, "Apr19", "L"); leg1->AddEntry(gh0, "June17v1", "L"); leg1->AddEntry(cutg_lep,"LEP exclusion", "F"); leg1->Draw("same"); cutg_lep->Draw("fsame"); cutg_lep->Draw("lsame"); l->DrawLatex(0.15,0.87,"ATLAS Internal, 2012"); l1->DrawLatex(0.15,0.81,"H#rightarrow #tau_{lep}#tau_{had}, m_{h}^{max}, #mu>0"); l2->DrawLatex(0.15,0.73,"#int L dt = 20.3 fb^{-1}, #sqrt{s}=8 TeV"); l3->DrawLatex(0.15,0.65,"high mass selection, June 18v1"); gPad->RedrawAxis(); c11->Print("limit_lh_2d_highmass_June18v1.png"); TCanvas *c111 = new TCanvas("c111","c111",800,800); h_axes->Draw("HIST"); gh0_j19v1->SetLineStyle(2); gh1_j19v1->SetLineStyle(0); gh2_j19v1->SetLineStyle(0); gh0_j19v1->SetLineWidth(2); gh1_j19v1->SetLineWidth(0); gh2_j19v1->SetLineWidth(0); gh1_j19v1->SetLineColor(0); gh2_j19v1->SetLineColor(0); gh1_j19v1->SetFillColor(kGreen); gh2_j19v1->SetFillColor(kYellow); gh0_j19v1->SetLineColor(1); gh0_j19v1->SetLineWidth(2); gh_j19v1->SetLineColor(1); gh_j19v1->SetMarkerColor(1); gh_j19v1->SetMarkerStyle(20); gh0_j19v1->Draw("L,same"); gh2_j19v1->Draw("E3,same"); gh1_j19v1->Draw("E3,same"); gh_j19v1->Draw("PL,same"); gh0_j19v1->Draw("L,same"); g1h0->SetLineColor(2); g1h0->SetLineStyle(9); g1h0->SetLineWidth(2); g1h0->Draw("L,same"); gh0->SetLineColor(4); gh0->SetLineStyle(2); gh0->SetLineWidth(3); gh0->Draw("L,same"); gh0_j18v1->SetLineColor(4); gh0_j18v1->SetLineStyle(9); gh0_j18v1->SetLineWidth(3); gh0_j18v1->Draw("L,same"); TLegend *leg1 = new TLegend(0.15+0.47, 0.55-0.3, 0.45+0.47, 0.77-0.3); leg1->SetFillColor(kWhite); leg1->SetFillStyle(0); leg1->SetBorderSize(0); // leg1->AddEntry(gl, "Observed CLs", "LP"); // leg1->AddEntry(gl0, "Expected CLs", "L"); leg1->AddEntry(gh_j18v1, "Expected CLs", "LP"); leg1->AddEntry(gh1_j18v1,"1 #sigma", "F"); leg1->AddEntry(gh2_j18v1,"2 #sigma", "F"); leg1->AddEntry(g1h0, "Apr19", "L"); leg1->AddEntry(gh0, "June17v1", "L"); leg1->AddEntry(gh0_j18v1, "June18v1", "L"); leg1->AddEntry(cutg_lep,"LEP exclusion", "F"); leg1->Draw("same"); cutg_lep->Draw("fsame"); cutg_lep->Draw("lsame"); l->DrawLatex(0.15,0.87,"ATLAS Internal, 2012"); l1->DrawLatex(0.15,0.81,"H#rightarrow #tau_{lep}#tau_{had}, m_{h}^{max}, #mu>0"); l2->DrawLatex(0.15,0.73,"#int L dt = 20.3 fb^{-1}, #sqrt{s}=8 TeV"); l3->DrawLatex(0.15,0.65,"high mass selection, June 19v1"); gPad->RedrawAxis(); c111->Print("limit_lh_2d_highmass_June19v1.png"); TCanvas *c112 = new TCanvas("c112","c112",800,800); h_axes->Draw("HIST"); gh0_j19v1_bc->SetLineStyle(2); gh1_j19v1_bc->SetLineStyle(0); gh2_j19v1_bc->SetLineStyle(0); gh0_j19v1_bc->SetLineWidth(2); gh1_j19v1_bc->SetLineWidth(0); gh2_j19v1_bc->SetLineWidth(0); gh1_j19v1_bc->SetLineColor(0); gh2_j19v1_bc->SetLineColor(0); gh1_j19v1_bc->SetFillColor(kGreen); gh2_j19v1_bc->SetFillColor(kYellow); gh0_j19v1_bc->SetLineColor(1); gh0_j19v1_bc->SetLineWidth(2); gh_j19v1_bc->SetLineColor(1); gh_j19v1_bc->SetMarkerColor(1); gh_j19v1_bc->SetMarkerStyle(20); gh0_j19v1_bc->Draw("L,same"); gh2_j19v1_bc->Draw("E3,same"); gh1_j19v1_bc->Draw("E3,same"); gh_j19v1_bc->Draw("PL,same"); gh0_j19v1_bc->Draw("L,same"); g1h0->SetLineColor(2); g1h0->SetLineStyle(9); g1h0->SetLineWidth(2); g1h0->Draw("L,same"); // gh0->SetLineColor(4); // gh0->SetLineStyle(2); // gh0->SetLineWidth(3); // gh0->Draw("L,same"); gh0_j19v1->SetLineColor(4); gh0_j19v1->SetLineStyle(9); gh0_j19v1->SetLineWidth(3); gh0_j19v1->Draw("L,same"); TLegend *leg1 = new TLegend(0.15+0.47, 0.55-0.3, 0.45+0.47, 0.77-0.3); leg1->SetFillColor(kWhite); leg1->SetFillStyle(0); leg1->SetBorderSize(0); // leg1->AddEntry(gl, "Observed CLs", "LP"); // leg1->AddEntry(gl0, "Expected CLs", "L"); leg1->AddEntry(gh_j18v1, "Expected CLs", "LP"); leg1->AddEntry(gh1_j18v1,"1 #sigma", "F"); leg1->AddEntry(gh2_j18v1,"2 #sigma", "F"); leg1->AddEntry(g1h0, "Apr19", "L"); leg1->AddEntry(gh0_j19v1, "June18v1", "L"); leg1->AddEntry(cutg_lep,"LEP exclusion", "F"); leg1->Draw("same"); cutg_lep->Draw("fsame"); cutg_lep->Draw("lsame"); l->DrawLatex(0.15,0.87,"ATLAS Internal, 2012"); l1->DrawLatex(0.15,0.81,"H#rightarrow #tau_{lep}#tau_{had}, m_{h}^{max}, #mu>0"); l2->DrawLatex(0.15,0.73,"#int L dt = 20.3 fb^{-1}, #sqrt{s}=8 TeV"); l3->DrawLatex(0.15,0.65,"high mass selection, June 19v1 bug corr"); gPad->RedrawAxis(); c112->Print("limit_lh_2d_highmass_June19v1_bc.png"); TCanvas *c12 = new TCanvas("c12","c12",800,800); h_axes->Draw("HIST"); gh0_j18v1->SetLineStyle(2); gh1_j18v1->SetLineStyle(0); gh2_j18v1->SetLineStyle(0); gh0_j18v1->SetLineWidth(2); gh1_j18v1->SetLineWidth(0); gh2_j18v1->SetLineWidth(0); gh1_j18v1->SetLineColor(0); gh2_j18v1->SetLineColor(0); gh1_j18v1->SetFillColor(kGreen); gh2_j18v1->SetFillColor(kYellow); gh0_j18v1->SetLineColor(1); gh0_j18v1->SetLineWidth(2); gh_j18v1->SetLineColor(1); gh_j18v1->SetMarkerColor(1); gh_j18v1->SetMarkerStyle(20); gh0_j18v1->Draw("L,same"); gh2_j18v1->Draw("E3,same"); gh1_j18v1->Draw("E3,same"); gh_j18v1->Draw("PL,same"); gh0_j18v1->Draw("L,same"); ghe0_j18v1->SetLineColor(2); ghe0_j18v1->SetLineStyle(9); ghe0_j18v1->SetLineWidth(2); ghe0_j18v1->Draw("L,same"); ghm0_j18v1->SetLineColor(4); ghm0_j18v1->SetLineStyle(2); ghm0_j18v1->SetLineWidth(3); ghm0_j18v1->Draw("L,same"); TLegend *leg1 = new TLegend(0.15+0.47, 0.55-0.3, 0.45+0.47, 0.77-0.3); leg1->SetFillColor(kWhite); leg1->SetFillStyle(0); leg1->SetBorderSize(0); // leg1->AddEntry(gl, "Observed CLs", "LP"); // leg1->AddEntry(gl0, "Expected CLs", "L"); leg1->AddEntry(gh_j18v1, "Expected CLs", "LP"); leg1->AddEntry(gh1_j18v1,"1 #sigma", "F"); leg1->AddEntry(gh2_j18v1,"2 #sigma", "F"); leg1->AddEntry(ghe0_j18v1, "elec", "L"); leg1->AddEntry(ghm0_j18v1, "muon", "L"); leg1->AddEntry(cutg_lep,"LEP exclusion", "F"); leg1->Draw("same"); cutg_lep->Draw("fsame"); cutg_lep->Draw("lsame"); l->DrawLatex(0.15,0.87,"ATLAS Internal, 2012"); l1->DrawLatex(0.15,0.81,"H#rightarrow #tau_{lep}#tau_{had}, m_{h}^{max}, #mu>0"); l2->DrawLatex(0.15,0.73,"#int L dt = 20.3 fb^{-1}, #sqrt{s}=8 TeV"); l3->DrawLatex(0.15,0.65,"high mass selection, June 18v1"); gPad->RedrawAxis(); c12->Print("limit_lh_2d_highmass_June18v1_emu.png"); // veto TCanvas *c2 = new TCanvas("c2","c2",800,800); h_axes1->Draw("HIST"); gv0->SetLineStyle(2); gv1->SetLineStyle(0); gv2->SetLineStyle(0); gv0->SetLineWidth(2); gv1->SetLineWidth(0); gv2->SetLineWidth(0); gv1->SetLineColor(0); gv2->SetLineColor(0); gv1->SetFillColor(kGreen); gv2->SetFillColor(kYellow); gv0->SetLineColor(1); gv0->SetLineWidth(2); gv->SetLineColor(1); gv->SetMarkerColor(1); gv->SetMarkerStyle(20); gv0->Draw("L,same"); gv2->Draw("E3,same"); gv1->Draw("E3,same"); gv->Draw("PL,same"); gv0->Draw("L,same"); g1v0->SetLineColor(2); g1v0->SetLineStyle(9); g1v0->SetLineWidth(2); g1v0->Draw("L,same"); TLegend *leg1 = new TLegend(0.15+0.47, 0.55+0.15, 0.45+0.47, 0.77+0.15); leg1->SetFillColor(kWhite); leg1->SetFillStyle(0); leg1->SetBorderSize(0); // leg1->AddEntry(gl, "Observed CLs", "LP"); // leg1->AddEntry(gl0, "Expected CLs", "L"); leg1->AddEntry(gv, "Expected CLs", "LP"); leg1->AddEntry(gv1,"1 #sigma", "F"); leg1->AddEntry(gv2,"2 #sigma", "F"); leg1->AddEntry(g1v0, "April 19", "L"); leg1->AddEntry(cutg_lep,"LEP exclusion", "F"); leg1->Draw("same"); cutg_lep->Draw("fsame"); cutg_lep->Draw("lsame"); l->DrawLatex(0.15,0.87,"ATLAS Internal, 2012"); l1->DrawLatex(0.15,0.81,"H#rightarrow #tau_{lep}#tau_{had}, m_{h}^{max}, #mu>0"); l2->DrawLatex(0.15,0.73,"#int L dt = 20.3 fb^{-1}, #sqrt{s}=8 TeV"); l3->DrawLatex(0.15,0.65,"veto category, June 17v1"); gPad->RedrawAxis(); c2->Print("limit_lh_2d_veto.png"); //c1->Print("limit_lh_2d_hihgmass.eps"); //return; // tag TCanvas *c3 = new TCanvas("c3","c3",800,800); h_axes1->Draw("HIST"); gt0->SetLineStyle(2); gt1->SetLineStyle(0); gt2->SetLineStyle(0); gt0->SetLineWidth(2); gt1->SetLineWidth(0); gt2->SetLineWidth(0); gt1->SetLineColor(0); gt2->SetLineColor(0); gt1->SetFillColor(kGreen); gt2->SetFillColor(kYellow); gt0->SetLineColor(1); gt0->SetLineWidth(2); gt->SetLineColor(1); gt->SetMarkerColor(1); gt->SetMarkerStyle(20); gt0->Draw("L,same"); gt2->Draw("E3,same"); gt1->Draw("E3,same"); gt->Draw("PL,same"); gt0->Draw("L,same"); g1t0->SetLineColor(2); g1t0->SetLineStyle(9); g1t0->SetLineWidth(2); g1t0->Draw("L,same"); TLegend *leg1 = new TLegend(0.15+0.47, 0.55+0.15, 0.45+0.47, 0.77+0.15); leg1->SetFillColor(kWhite); leg1->SetFillStyle(0); leg1->SetBorderSize(0); leg1->AddEntry(gt, "Expected CLs", "LP"); leg1->AddEntry(gt1,"1 #sigma", "F"); leg1->AddEntry(gt2,"2 #sigma", "F"); leg1->AddEntry(g1t0, "April 19", "L"); leg1->AddEntry(cutg_lep,"LEP exclusion", "F"); leg1->Draw("same"); cutg_lep->Draw("fsame"); cutg_lep->Draw("lsame"); l->DrawLatex(0.15,0.87,"ATLAS Internal, 2012"); l1->DrawLatex(0.15,0.81,"H#rightarrow #tau_{lep}#tau_{had}, m_{h}^{max}, #mu>0"); l2->DrawLatex(0.15,0.73,"#int L dt = 20.3 fb^{-1}, #sqrt{s}=8 TeV"); l3->DrawLatex(0.15,0.65,"tag category, June 17v1"); gPad->RedrawAxis(); c3->Print("limit_lh_2d_tag.png"); TCanvas *c4 = new TCanvas("c4","c4",800,800); h_axes1->Draw("HIST"); gl0->SetLineStyle(2); gl1->SetLineStyle(0); gl2->SetLineStyle(0); gl0->SetLineWidth(2); gl1->SetLineWidth(0); gl2->SetLineWidth(0); gl1->SetLineColor(0); gl2->SetLineColor(0); gl1->SetFillColor(kGreen); gl2->SetFillColor(kYellow); gl0->SetLineColor(1); gl0->SetLineWidth(2); gl->SetLineColor(1); gl->SetMarkerColor(1); gl->SetMarkerStyle(20); gl0->Draw("L,same"); gl2->Draw("E3,same"); gl1->Draw("E3,same"); gl->Draw("PL,same"); gl0->Draw("L,same"); g1l0->SetLineColor(2); g1l0->SetLineStyle(9); g1l0->SetLineWidth(2); g1l0->Draw("L,same"); TLegend *leg1 = new TLegend(0.15+0.47, 0.55+0.15, 0.45+0.47, 0.77+0.15); leg1->SetFillColor(kWhite); leg1->SetFillStyle(0); leg1->SetBorderSize(0); leg1->AddEntry(gl, "Expected CLs", "LP"); leg1->AddEntry(gl1,"1 #sigma", "F"); leg1->AddEntry(gl2,"2 #sigma", "F"); leg1->AddEntry(g1l0, "April 19", "L"); leg1->AddEntry(cutg_lep,"LEP exclusion", "F"); leg1->Draw("same"); cutg_lep->Draw("fsame"); cutg_lep->Draw("lsame"); l->DrawLatex(0.15,0.87,"ATLAS Internal, 2012"); l1->DrawLatex(0.15,0.81,"H#rightarrow #tau_{lep}#tau_{had}, m_{h}^{max}, #mu>0"); l2->DrawLatex(0.15,0.73,"#int L dt = 20.3 fb^{-1}, #sqrt{s}=8 TeV"); l3->DrawLatex(0.15,0.65,"low mass selection, June 17v1"); gPad->RedrawAxis(); c4->Print("limit_lh_2d_low.png"); // cross section limits TH1F *h_axes = new TH1F("","", 100, 90, 800); h_axes->Fill(-10, 1000); h_axes->GetYaxis()->SetTitle("#sigma BR(#phi) [pb]"); h_axes->GetXaxis()->SetTitle("m_{#phi} [GeV/c^{2}]"); h_axes->GetYaxis()->SetRangeUser(0,40); h_axes->SetLineColor(0); TH1F *h_axes1 = new TH1F("","", 100, 90, 350); h_axes1->Fill(-10, 1000); h_axes1->GetYaxis()->SetTitle("#sigma BR(#phi) [pb]"); h_axes1->GetXaxis()->SetTitle("m_{#phi} [GeV/c^{2}]"); h_axes1->GetYaxis()->SetRangeUser(0,40); h_axes1->SetLineColor(0); // x1bh plotMaker(h_axes, x1bh, "April 19", xbh, xbh0, xbh1, xbh2, 1, 0.01, 100., "bb#phi#rightarrow #tau_{lep}#tau_{had} high June 17v1", "limit_lh_high_bbA", true); plotMaker(h_axes, x1gh, "April 19", xgh, xgh0, xgh1, xgh2, 1, 0.01, 100., "gg#phi#rightarrow #tau_{lep}#tau_{had} high June 17v1", "limit_lh_high_ggA", true); plotMaker(h_axes1, x1bl, "April 19", xbl, xbl0, xbl1, xbl2, 1, 0.1, 1000., "bb#phi#rightarrow #tau_{lep}#tau_{had} low June 17v1", "limit_lh_low_bbA", true); plotMaker(h_axes1, x1gl, "April 19", xgl, xgl0, xgl1, xgl2, 1, 0.1, 1000., "gg#phi#rightarrow #tau_{lep}#tau_{had} low June 17v1", "limit_lh_low_ggA", true); plotMaker(h_axes1, x1bv, "April 19", xbv, xbv0, xbv1, xbv2, 1, 0.1, 1000., "bb#phi#rightarrow #tau_{lep}#tau_{had} veto June 17v1", "limit_lh_veto_bbA", true); plotMaker(h_axes1, x1gv, "April 19", xgv, xgv0, xgv1, xgv2, 1, 0.1, 1000., "gg#phi#rightarrow #tau_{lep}#tau_{had} veto June 17v1", "limit_lh_veto_ggA", true); plotMaker(h_axes1, x1bt, "April 19", xbt, xbt0, xbt1, xbt2, 1, 0.1, 1000., "bb#phi#rightarrow #tau_{lep}#tau_{had} tag June 17v1", "limit_lh_tag_bbA", true); plotMaker(h_axes1, x1gt, "April 19", xgt, xgt0, xgt1, xgt2, 1, 0.1, 1000., "gg#phi#rightarrow #tau_{lep}#tau_{had} tag June 17v1", "limit_lh_tag_ggA", true); plotMaker(h_axes, xbh, "June 17v1", xbh_j18v1, xbh0_j18v1, xbh1_j18v1, xbh2_j18v1, 1, 0.01, 100., "bb#phi#rightarrow #tau_{lep}#tau_{had} high June 18v1", "limit_lh_high_bbA_j18v1", true); plotMaker(h_axes, xgh, "June 17v1", xgh_j18v1, xgh0_j18v1, xgh1_j18v1, xgh2_j18v1, 1, 0.01, 100., "gg#phi#rightarrow #tau_{lep}#tau_{had} high June 18v1", "limit_lh_high_ggA_j18v1", true); return; }