AUI_ERRCODE aui_Region::Draw( aui_Surface *surface, sint32 x, sint32 y )
{
AUI_ERRCODE errcode = AUI_ERRCODE_OK;
if ( !IsHidden() )
{
errcode = DrawThis( surface, x, y );
if ( errcode == AUI_ERRCODE_OK )
DrawChildren( surface, x, y );
}
m_draw = 0;
return AUI_ERRCODE_OK;
}
void Plotter_sbatb2(){
Int_t mA=248;
TString strmA = "248";
Int_t mH=500;
TString strmH = "500";
gStyle->SetOptStat(0);
gStyle->SetOptTitle(0);
gStyle->SetPadLeftMargin(0.1191275);
gStyle->SetPadRightMargin(0.05944056);
gStyle->SetPadTopMargin(0.05944056);
gStyle->SetPadBottomMargin(0.1206294);
gStyle->SetTitleSize(0.04, "XYZ");
gStyle->SetTitleXOffset(1.1);
gStyle->SetTitleYOffset(1.45);
gStyle->SetPalette(1);
gStyle->SetNdivisions(505);
gStyle->SetLabelSize(0.04,"XYZ");
// Getting the limits from combine root files
TString file_path_sig1 = "/home/fynu/amertens/storage/CMSSW/CMSSW_6_1_1/src/HiggsAnalysis/CombinedLimit/python/v1_MH_"+strmH+"_Asymptotic_Unblind_Combined_/higgsCombineTest.Asymptotic.mA"+strmA+".root";
Double_t limit;
TGraph2D* myTGraph_ratio = new TGraph2D(10);
myTGraph_ratio->SetName("Ratio CMSSW/Delphes");
myTGraph_ratio->SetTitle("Ratio CMSSW/Delphes");
myTGraph_ratio->SetPoint(0,35,142,0.411);
myTGraph_ratio->SetPoint(1,10,200,0.838);
myTGraph_ratio->SetPoint(2,110,200,0.796);
myTGraph_ratio->SetPoint(3,30,329,0.826);
myTGraph_ratio->SetPoint(4,70,329,0.897);
myTGraph_ratio->SetPoint(5,70,575,0.955);
myTGraph_ratio->SetPoint(6,70,875,0.907);
myTGraph_ratio->SetPoint(7,142,329,0.925);
myTGraph_ratio->SetPoint(8,142,875,0.927);
myTGraph_ratio->SetPoint(9,378,575,0.891);
myTGraph_ratio->SetPoint(10,378,875,0.911);
myTGraph_ratio->SetPoint(11,575,875,0.890);
myTGraph_ratio->SetPoint(12,761,875,0.849);
myTGraph_ratio->SetPoint(13,50,1200,0.9);
myTGraph_ratio->SetPoint(14,1200,1200,0.85);
Double_t ratio = myTGraph_ratio->Interpolate(mA,mH);
pathDelphesEff = "/home/fynu/amertens/storage/THDM/eff_v1/"+strmH+"_"+strmA+".root";
DelphesFile = TFile(pathDelphesEff);
TGraph2D* myG2D = DelphesFile.Get("Graph2D");
Double_t eff = myG2D->Interpolate(int(mA),int(mH));
Double_t lumi = 19.7;
cout << "efficiency : " << eff << " ratio : " << ratio << endl;
Double_t LumiEff = eff*ratio*lumi;
TFile* f_Sig1 = new TFile(file_path_sig1);
TTree* tree_Sig1 = (TTree*) f_Sig1->Get("limit");
TBranch *branch = tree_Sig1->GetBranch("limit");
branch->SetAddress(&limit);
branch->GetEntry(5);
Double_t limit_obs = limit /(0.067*LumiEff);
branch->GetEntry(0);
Double_t limit_m2 = limit /(0.067*LumiEff);
branch->GetEntry(1);
Double_t limit_m1 = limit /(0.067*LumiEff);
branch->GetEntry(2);
Double_t limit_exp = limit /(0.067*LumiEff);
branch->GetEntry(3);
Double_t limit_p1 = limit /(0.067*LumiEff);
branch->GetEntry(4);
Double_t limit_p2 = limit /(0.067*LumiEff);
// Limit vs TanBeta cosba
TString Xsec_path = "/home/fynu/amertens/Delphes/delphes/condor/xsec_H_ZA_ll_bb_tautau_tanb_cosba_"+strmH+"_"+strmA+".root";
TFile *f4 = new TFile(Xsec_path);
TCanvas* C_tbcba = new TCanvas("C_tbcba","C_tbcba",600,600);
gPad->SetGrid();
C_tbcba->SetLogy(true);
//C_tbcba->SetRightMargin(0.17);
TH2D* h_tbcba = f4->Get("h_HZbbXsec");
TGraph* g_obs = getContourFilledX(h_tbcba, C_tbcba, 3, 1,3645, limit_obs);
TGraph* g_exp = getContourFilledX(h_tbcba, C_tbcba, 3, 7,0, limit_exp);
Double_t contours[6];
contours[0] = limit_m2;
contours[1] = limit_m1;
contours[2] = limit_exp;
contours[3] = limit_p1;
contours[4] = limit_p2;
h_tbcba->SetContour(5, contours);
TH2D* finalPlot_tbcba = h_tbcba;
Int_t colors[] = {kYellow, kGreen, kGreen, kYellow, kWhite}; // #colors >= #levels - 1
gStyle->SetPalette((sizeof(colors)/sizeof(Int_t)), colors);
DrawThis();
finalPlot_tbcba->Draw("cont list same");
cout << "limits : " << limit_m2 << endl;
cout << "limits : " << limit_m1 << endl;
cout << "limits : " << limit_exp << endl;
cout << "limits : " << limit_p1 << endl;
cout << "limits : " << limit_p2 << endl;
cout << "limits obs : " << limit_obs << endl;
// TGraph* g_obs = getContourFilledX(finalPlot_tbcba, C_tbcba, 3, 1,3004, limit_obs);
// TGraph* g_exp = getContourFilledX(finalPlot_tbcba, C_tbcba, 3, 7,4000, limit_exp);
g_obs->Draw("CL F same");
g_exp->Draw("CL F same");
DrawMasses("M_{H} = 330 GeV", "M_{A} = 100 GeV");
TGraph* g_yellow = new TGraph();
g_yellow->SetFillColor(kYellow);
TGraph* g_green = new TGraph();
g_green->SetFillColor(kGreen);
TLegend* leg = new TLegend(0.7,0.5,0.9,0.7);
//leg->SetHeader("#splitline{THDM type II}{#splitline{M_{H} = 378}{M_{A} = 216}}");
leg->SetLineColor(0);
leg->SetFillColor(0);
leg->AddEntry(g_obs,"Obs. Excl.","F");
leg->AddEntry(g_exp,"Exp. Excl.","l");
leg->AddEntry(g_green,"1-sigma","F");
leg->AddEntry(g_yellow,"2-sigma","F");
leg->Draw();
gPad->RedrawAxis("g"); // for "grid lines"
}
