#include <vector>

#include <cstdio>

#include <iostream> // std::cout

#include <string> // std::string, std::to_string

#include <tdrstyle.C>

#include <CMS_lumi.C>

using namespace std;

void shift_bin(TH1* input, TH1* output){

char tempname[200];

char temptitle[200];

output->SetName(tempname);

output->SetTitle(temptitle);

output->SetBins(input->GetNbinsX(),input->GetBinLowEdge(1)-0.5,input->GetBinLowEdge(input->GetNbinsX()+1)-0.5);

//input->Print("all");

//output = new TH1D(tempname,temptitle,input->GetNbinsX(),input->GetBinLowEdge(1)-0.5,input->GetBinLowEdge(input->GetNbinsX()+1)-0.5);

// 0: underflow

// 1: first bin [Use the lowedge of this bin]

// input->GetNbinsX(): highest bin

// input->GetNbinsX()+1: overflow bin [use the lowedge of this bin]

//

for (int ibin=1;ibin<=input->GetNbinsX();ibin++){

output->SetBinContent(ibin,input->GetBinContent(ibin));

output->SetBinError(ibin,input->GetBinError(ibin));

//std::cout << input->GetBinContent(ibin) << std::endl;

}

}

void Plot_searchBin_full_ICHEP_wQCD(string option="", int pull=0){ // string option="QCD"

// Use option="QCD" to produce plots in QCD binning

char tempname[200];

// Open root file

sprintf(tempname,"LLPrediction_QCD.root");

// true: do closure test (MC prediction vs MC truth)

// false: do data driven prediction and compare to MC truth

bool doDataVsMC = false;

// Add systematics in quadrature to stat. uncertainty on prediction

// Non-closure systematic not included yet!

bool showSystematics = false;

bool doClosurewoIsoTrackVeto = false;

///////////////////////////////////////////////////////////////////////////////////////////

////Some cosmetic work for official documents.

//

// Set basic style

//

setTDRStyle();

gStyle->SetPalette(1) ; // for better color output

//

// Canvas size

int W = 1200;

int H = 740;

int H_ref = 740;

int W_ref = 800;

float T = 0.10*H_ref;

float B = 0.06*H_ref;

float L = 0.16*W_ref;

float R = 0.04*W_ref;

//

// Various vertical line coordinates

float ymax_top = 40000.;

float ymin_top = 0.015;

float ymax2_top = 1000.;

float ymax3_top = 200.;

float ymax4_top = 30.;

float ymax5_top = 5.;

float ymax_bottom = 1.99;

float ymin_bottom = 0.01;

float ymax2_bottom = 2.15;

float ymax3_bottom = 2.15;

float ymax4_bottom = 2.15;

//

// Luminosity information for scaling

double lumi = 7.6; // normaliza to this lumi (fb-1)

double lumi_ref = 7.6; // normaliza to 3 (fb-1)

///////////////////////////////////////////////////////////////////////////////////////////

//

// More specific style set, opening input files etc

gStyle->SetOptStat(0); ///to avoid the stat. on the plots

//gStyle->SetErrorX(0);

char xtitlename[200];

char ytitlename[200];

TFile * LLFile = new TFile(tempname,"R");

printf("Opened %s\n",tempname);

//

// Define legend

//

Float_t legendX1 = .655; //.50;

Float_t legendX2 = .955; //.70;

Float_t legendY1 = .53; //.65;

Float_t legendY2 = .78;

TLegend* catLeg1 = new TLegend(legendX1,legendY1,legendX2,legendY2);

//catLeg1->SetTextSize(0.060);

catLeg1->SetTextSize(0.044);

catLeg1->SetTextFont(42);

catLeg1->SetFillColor(0);

catLeg1->SetLineColor(1);

catLeg1->SetBorderSize(1);

//

// Define canvas

//

TCanvas *canvas = new TCanvas("canvas","canvas",10,10,W,H);

canvas->SetFillColor(0);

canvas->SetBorderMode(0);

canvas->SetFrameFillStyle(0);

canvas->SetFrameBorderMode(0);

canvas->SetLeftMargin( L/W );

canvas->SetRightMargin( R/W );

canvas->SetTopMargin( T/H );

canvas->SetBottomMargin( B/H );

canvas->SetTickx(0);

canvas->SetTicky(0);

canvas->Divide(1, 2);

//

// Define pads

//

TPad* canvas_up = (TPad*) canvas->GetListOfPrimitives()->FindObject("canvas_1");

TPad* canvas_dw = (TPad*) canvas->GetListOfPrimitives()->FindObject("canvas_2");

//

// define the size

double up_height = 0.8; // please tune so that the upper figures size will meet your requirement

double dw_correction = 1.30; // please tune so that the smaller canvas size will work in your environment

double font_size_dw = 0.1; // please tune the font size parameter for bottom figure

double dw_height = (1. - up_height) * dw_correction;

double dw_height_offset = 0.04; // KH, added to put the bottom one closer to the top panel

//

// set pad size

canvas_up->SetPad(0., 1 - up_height, 1., 1.00);

canvas_dw->SetPad(0., 0., 1., dw_height+dw_height_offset);

//

canvas_up->SetFrameFillColor(0);

canvas_up->SetFillColor(0);

canvas_up->SetTopMargin(0.12);

canvas_up->SetLeftMargin(0.1);

//

canvas_dw->SetFillColor(0);

canvas_dw->SetFrameFillColor(0);

canvas_dw->SetBottomMargin(0.35);

canvas_dw->SetTopMargin(0);

canvas_dw->SetLeftMargin(0.1);

//

// draw top figure

canvas_up->cd();

TH1D * GenHist, * EstHist,* thist;

TH1D * GenHistTemp, * EstHistTemp;

TH1D * GenHistD, * EstHistD;

TH1D * GenHistDTemp, * EstHistDTemp;

TH1D * EstSystematics;

TH1D * histTemplate;

double HT_x_max=2500.;

double HT_x_min=400.;

double MHT_x_max=1000.;

double NJet_x_max=15.;

double NBtag_x_max=4.;

double search_x_max=208.-0.5;

if(option.find("QCD")!=string::npos)search_x_max=208.-0.5;

double search_x_min=1.-0.5;

TDirectory *dPre = 0;

TDirectory *dExp = (TDirectory*)LLFile->Get("Expectation");

if(doDataVsMC){

dPre = (TDirectory*)LLFile->Get("Prediction_data");

}else{

dPre = (TDirectory*)LLFile->Get("Prediction_MC");

}

if(doDataVsMC){

EstHistTemp=(TH1D*) dPre->Get("totalPred_LL")->Clone();

EstHistDTemp=(TH1D*) dPre->Get("totalPred_LL")->Clone();

}else{

if(doClosurewoIsoTrackVeto){

EstHistTemp=(TH1D*) dPre->Get("totalPred_woIsoTrack_LL_MC")->Clone();

EstHistDTemp=(TH1D*) dPre->Get("totalPred_woIsoTrack_LL_MC")->Clone();

}else{

EstHistTemp=(TH1D*) dPre->Get("totalPred_LL_MC")->Clone();

EstHistDTemp=(TH1D*) dPre->Get("totalPred_LL_MC")->Clone();

}

}

if(doClosurewoIsoTrackVeto){

GenHistTemp=(TH1D*) dExp->Get("totalExp_woIsoTrack_LL")->Clone();

GenHistDTemp=(TH1D*) dExp->Get("totalExp_woIsoTrack_LL")->Clone();;

}else{

GenHistTemp=(TH1D*) dExp->Get("totalExp_LL")->Clone();

GenHistDTemp=(TH1D*) dExp->Get("totalExp_LL")->Clone();

}

if(showSystematics){

TDirectory *dSys = (TDirectory*)LLFile->Get("AdditionalContent");

if(doDataVsMC){

EstSystematics=(TH1D*) dSys->Get("totalPropSysUp_LL")->Clone();

}else{

EstSystematics=(TH1D*) dSys->Get("totalPropSysUp_LL_MC")->Clone();

}

}

if(EstHistTemp->GetNbinsX() != GenHistTemp->GetNbinsX()) std::cout<<"NbinsX of Expectation and Prediction don't agree!"<<std::endl;

EstHist = new TH1D("Exp", "Exp", EstHistTemp->GetNbinsX(), 0.5, EstHistTemp->GetNbinsX()+0.5);

EstHistD = new TH1D("ExpD", "ExpD", EstHistDTemp->GetNbinsX(), 0.5, EstHistDTemp->GetNbinsX()+0.5);

GenHist = new TH1D("Pred", "Pred", GenHistTemp->GetNbinsX(), 0.5, GenHistTemp->GetNbinsX()+0.5);

GenHistD = new TH1D("PredD", "PredD", GenHistDTemp->GetNbinsX(), 0.5, GenHistDTemp->GetNbinsX()+0.5);

for(int i = 0; i <= EstHistTemp->GetNbinsX()+1; i++){

EstHist->SetBinContent(i, EstHistTemp->GetBinContent(i));

EstHistD->SetBinContent(i, EstHistDTemp->GetBinContent(i));

if(showSystematics){

EstHist->SetBinError(i, std::sqrt(EstHistTemp->GetBinError(i)*EstHistTemp->GetBinError(i)+EstSystematics->GetBinContent(i)*EstSystematics->GetBinContent(i)*EstHistTemp->GetBinError(i)*EstHistTemp->GetBinError(i)));

EstHistD->SetBinError(i, std::sqrt(EstHistDTemp->GetBinError(i)*EstHistDTemp->GetBinError(i)+EstSystematics->GetBinContent(i)*EstSystematics->GetBinContent(i)*EstHistDTemp->GetBinError(i)*EstHistDTemp->GetBinError(i)));

}else{

EstHist->SetBinError(i, EstHistTemp->GetBinError(i));

EstHistD->SetBinError(i, EstHistDTemp->GetBinError(i));

}

GenHist->SetBinContent(i, GenHistTemp->GetBinContent(i));

GenHist->SetBinError(i, GenHistTemp->GetBinError(i));

GenHistD->SetBinContent(i, GenHistDTemp->GetBinContent(i));

GenHistD->SetBinError(i, GenHistDTemp->GetBinError(i));

}

GenHist->SetLineColor(4);

EstHist->SetLineColor(4);

//GenHist->GetXaxis()->SetLabelFont(42);

//GenHist->GetXaxis()->SetLabelOffset(0.007);

//GenHist->GetXaxis()->SetLabelSize(0.04);

//GenHist->GetXaxis()->SetTitleSize(0.05);

//GenHist->GetXaxis()->SetTitleOffset(0.9);

//GenHist->GetXaxis()->SetTitleOffset(0.5);

//GenHist->GetXaxis()->SetTitleFont(42);

//GenHist->GetYaxis()->SetLabelFont(42);

//GenHist->GetYaxis()->SetLabelOffset(0.007);

//GenHist->GetYaxis()->SetLabelSize(0.04);

GenHist->GetYaxis()->SetLabelSize(0.045*1.15);

GenHist->GetYaxis()->SetTitleSize(0.06*1.15);

GenHist->GetYaxis()->SetTitleOffset(0.6);

GenHist->GetYaxis()->SetTitleFont(42);

//EstHist->GetXaxis()->SetLabelFont(42);

//EstHist->GetXaxis()->SetLabelOffset(0.007);

//EstHist->GetXaxis()->SetLabelSize(0.04);

//EstHist->GetXaxis()->SetTitleSize(0.05);

//EstHist->GetXaxis()->SetTitleOffset(0.9);

//EstHist->GetXaxis()->SetTitleFont(42);

//EstHist->GetYaxis()->SetLabelFont(42);

//EstHist->GetYaxis()->SetLabelOffset(0.007);

//EstHist->GetYaxis()->SetLabelSize(0.04);

//EstHist->GetYaxis()->SetTitleSize(0.08);

//EstHist->GetYaxis()->SetTitleOffset(2.0);

//EstHist->GetYaxis()->SetTitleFont(42);

sprintf(xtitlename,"Search region bin number");

sprintf(ytitlename,"Events");

gPad->SetLogy();

GenHist->SetMaximum(ymax_top);

GenHist->SetMinimum(ymin_top);

GenHist->GetXaxis()->SetRangeUser(search_x_min,search_x_max);

//GenHist->GetYaxis()->SetTickLength(0.015);

//GenHist->GetXaxis()->SetTickLength(0.02);

//gPad->SetGridx(1);

TExec *ex1 = new TExec("ex1","gStyle->SetErrorX(0);");

TExec *ex2 = new TExec("ex2","gStyle->SetErrorX(0.5);");

GenHist->SetTitle("");

GenHist->SetMarkerStyle(20);

GenHist->SetMarkerSize(1.2);

GenHist->SetLineColor(1);

GenHist->GetXaxis()->SetTitle(xtitlename);

GenHist->GetYaxis()->SetTitle(ytitlename);

GenHist->Scale(lumi/lumi_ref);

EstHist->Scale(lumi/lumi_ref);

TH1D * GenHist_Normalize = static_cast<TH1D*>(GenHist->Clone("GenHist_Normalize"));

GenHist_Normalize->SetMaximum(ymax_top);

GenHist_Normalize->SetMinimum(ymin_top);

ex1->Draw();

//GenHist_Normalize->GetListOfFunctions()->Add(ex1);

GenHist_Normalize->DrawCopy("e");

EstHist->SetFillStyle(3144);

EstHist->SetFillColor(kRed-10);

EstHist->SetMarkerStyle(20);

EstHist->SetMarkerSize(0.0001);

TH1D * EstHist_Normalize = static_cast<TH1D*>(EstHist->Clone("EstHist_Normalize"));

ex2->Draw();

//EstHist_Normalize->GetListOfFunctions()->Add(ex2);

EstHist_Normalize->DrawCopy("e2same");

//EstHist_Normalize->DrawCopy("esame");

TH1D *EstHist_Normalize_Clone = (TH1D*)EstHist_Normalize->Clone();

for(int i=1; i<208; i++) {

EstHist_Normalize_Clone->SetBinError(i,0);

}

EstHist_Normalize_Clone->SetFillColor(kWhite);

EstHist_Normalize_Clone->Draw("esame");

GenHist->Print("all");

EstHist->Print("all");

//

// Re-draw to have "expectation" on top of "prediction"

ex1->Draw();

GenHist_Normalize->DrawCopy("esame");

//

TString line = "";

sprintf(tempname,"%8.1f",lumi);

line+=tempname;

line+=" fb^{-1} (13 TeV)";

int iPeriod = 0; // 1=7TeV, 2=8TeV, 3=7+8TeV, 7=7+8+13TeV, 0=free form (uses lumi_sqrtS)

int iPos=0;

writeExtraText = true;

if(doDataVsMC) extraText = " Preliminary";

else extraText = " Simulation";

//float extraTextFont = 52; // default is helvetica-italics

// text sizes and text offsets with respect to the top frame

// in unit of the top margin size

//lumiTextSize = 0.5;

//float lumiTextOffset = 0.2;

//cmsTextSize = 0.65;

//float cmsTextOffset = 0.1; // only used in outOfFrame version

//relPosX = 0.045;

//relPosY = 0.035;

//relExtraDY = 1.2;

// ratio of "CMS" and extra text size

//float extraOverCmsTextSize = 0.76;

//TString lumi_13TeV = "20.1 fb^{-1}";

//TString lumi_8TeV = "19.7 fb^{-1}";

//TString lumi_7TeV = "5.1 fb^{-1}";

TString lumi_sqrtS = line;

//

if(option.find("QCD")==string::npos ){

//-----------------------------------------------------------

// Putting lines and labels explaining search region definitions

//-----------------------------------------------------------

//TString CMSlabel = "";

//cmsText = "#bf{CMS} #it{Simulation}";

//CMSlabel += "#splitline{#bf{CMS}}{#scale[0.6]{#it{Simulation}}}";

/*

double x0 = gStyle->GetPadLeftMargin();

double x1 = 1.-gStyle->GetPadRightMargin();

double y0 = 1.005-gStyle->GetPadTopMargin();

double y1 = 0.96;

TPaveText *Lumitxt = new TPaveText(x0,y0,x1,y1,"NDC");

Lumitxt->SetBorderSize(0);

Lumitxt->SetFillColor(0);

Lumitxt->SetTextFont(42);

Lumitxt->SetTextAlign(31);

Lumitxt->SetTextSize(1.2*gStyle->GetPadTopMargin());

Lumitxt->SetMargin(0.);

Lumitxt->AddText(line);

//Lumitxt->Draw("same");

x0 = gStyle->GetPadLeftMargin()+0.03;

x1 = gStyle->GetPadLeftMargin()+0.13;

y0 = 0.905-gStyle->GetPadTopMargin();

y1 = 0.88;

TPaveText *CMStxt = new TPaveText(x0,y0,x1,y1,"NDC");

CMStxt->SetBorderSize(0);

CMStxt->SetFillColor(0);

CMStxt->SetTextFont(42);

CMStxt->SetTextAlign(11);

CMStxt->SetTextSize(1.2*gStyle->GetPadTopMargin());

CMStxt->SetMargin(0.);

CMStxt->AddText(CMSlabel);

//CMStxt->Draw("same");

*/

// Njet separation lines

TLine *tl_njet = new TLine();

tl_njet->SetLineStyle(2);

tl_njet->DrawLine(53.-0.5,ymin_top,53.-0.5,ymax_top);

tl_njet->DrawLine(105.-0.5,ymin_top,105.-0.5,ymax_top);

tl_njet->DrawLine(157.-0.5,ymin_top,157.-0.5,ymax_top);

// Njet labels

TLatex * ttext_njet = new TLatex();

ttext_njet->SetTextFont(42);

ttext_njet->SetTextSize(0.04);

ttext_njet->SetTextAlign(22);

ttext_njet->DrawLatex(26.-0.5 , ymax_top/4. , "3 #leq N_{#scale[0.2]{ }jet} #leq 4");

ttext_njet->DrawLatex(78.-0.5 , ymax_top/4. , "5 #leq N_{#scale[0.2]{ }jet} #leq 6");

ttext_njet->DrawLatex(130.-0.5 , ymax_top/4. , "7 #leq N_{#scale[0.2]{ }jet} #leq 8");

ttext_njet->DrawLatex(182.-0.5 , ymax_top/4. , "N_{#scale[0.2]{ }jet} #geq 9");

// Nb separation lines

TLine *tl_nb = new TLine();

tl_nb->SetLineStyle(3);

tl_nb->DrawLine(14.-0.5,ymin_top,14.-0.5,ymax2_top);

tl_nb->DrawLine(27.-0.5,ymin_top,27.-0.5,ymax2_top);

tl_nb->DrawLine(40.-0.5,ymin_top,40.-0.5,ymax2_top);

tl_nb->DrawLine(53.-0.5,ymin_top,53.-0.5,ymax2_top);

tl_nb->DrawLine(66.-0.5,ymin_top,66.-0.5,ymax3_top);

tl_nb->DrawLine(79.-0.5,ymin_top,79.-0.5,ymax3_top);

tl_nb->DrawLine(92.-0.5,ymin_top,92.-0.5,ymax3_top);

tl_nb->DrawLine(105.-0.5,ymin_top,105.-0.5,ymax3_top);

tl_nb->DrawLine(118.-0.5,ymin_top,118.-0.5,ymax4_top);

tl_nb->DrawLine(131.-0.5,ymin_top,131.-0.5,ymax4_top);

tl_nb->DrawLine(144.-0.5,ymin_top,144.-0.5,ymax4_top);

tl_nb->DrawLine(157.-0.5,ymin_top,157.-0.5,ymax4_top);

tl_nb->DrawLine(170.-0.5,ymin_top,170.-0.5,ymax5_top);

tl_nb->DrawLine(183.-0.5,ymin_top,183.-0.5,ymax5_top);

tl_nb->DrawLine(196.-0.5,ymin_top,196.-0.5,ymax5_top);

// Nb labels

TLatex * ttext_nb = new TLatex();

ttext_nb->SetTextFont(42);

ttext_nb->SetTextSize(0.04);

ttext_nb->SetTextAlign(22);

ttext_nb->DrawLatex(9.-0.5 , ymax_top/12. , "N_{#scale[0.2]{ }b-jet}");

ttext_nb->DrawLatex(8.-0.5 , ymax_top/40. , "0");

ttext_nb->DrawLatex(20.-0.5 , ymax_top/40. , "1");

ttext_nb->DrawLatex(33.-0.5 , ymax_top/40. , "2");

ttext_nb->DrawLatex(46.-0.5 , ymax_top/40. , "#geq 3");

//

} else {

//-----------------------------------------------------------

// Putting lines and labels explaining search region definitions

//-----------------------------------------------------------

// Njet separation lines

TLine *tl_njet = new TLine();

tl_njet->SetLineStyle(2);

tl_njet->DrawLine( 45.,ymin_top, 45.,ymax_top);

tl_njet->DrawLine( 89.,ymin_top, 89.,ymax_top);

tl_njet->DrawLine(133.,ymin_top,133.,ymax_top);

tl_njet->DrawLine(177.,ymin_top,177.,ymax_top);

// Njet labels

TLatex * ttext_njet = new TLatex();

ttext_njet->SetTextFont(42);

ttext_njet->SetTextSize(0.04);

ttext_njet->SetTextAlign(22);

ttext_njet->DrawLatex(23. , ymax_top/4. , "N_{jets} = 4");

ttext_njet->DrawLatex(67. , ymax_top/4. , "N_{jets} = 5");

ttext_njet->DrawLatex(111., ymax_top/4. , "N_{jets} = 6");

ttext_njet->DrawLatex(155., ymax_top/4. , "7 #leq N_{jets} #leq 8");

ttext_njet->DrawLatex(199., ymax_top/4. , "N_{jets} #geq 9");

// Nb separation lines

TLine *tl_nb = new TLine();

tl_nb->SetLineStyle(3);

tl_nb->DrawLine(12.,ymin_top,12.,ymax2_top);

tl_nb->DrawLine(23.,ymin_top,23.,ymax2_top);

tl_nb->DrawLine(34.,ymin_top,34.,ymax2_top);

tl_nb->DrawLine(56.,ymin_top,56.,ymax2_top);

tl_nb->DrawLine(67.,ymin_top,67.,ymax2_top);

tl_nb->DrawLine(78.,ymin_top,78.,ymax2_top);

tl_nb->DrawLine(100.,ymin_top,100.,ymax2_top);

tl_nb->DrawLine(111.,ymin_top,111.,ymax2_top);

tl_nb->DrawLine(122.,ymin_top,122.,ymax2_top);

tl_nb->DrawLine(144.,ymin_top,144.,ymax3_top);

tl_nb->DrawLine(155.,ymin_top,155.,ymax3_top);

tl_nb->DrawLine(166.,ymin_top,166.,ymax3_top);

tl_nb->DrawLine(188.,ymin_top,188.,ymax4_top);

tl_nb->DrawLine(199.,ymin_top,199.,ymax4_top);

tl_nb->DrawLine(210.,ymin_top,210.,ymax4_top);

// Nb labels

TLatex * ttext_nb = new TLatex();

ttext_nb->SetTextFont(42);

ttext_nb->SetTextSize(0.04);

ttext_nb->SetTextAlign(22);

ttext_nb->SetTextAngle(90);

ttext_nb->DrawLatex( 6. , ymax_top/50. , "N_{b} = 0");

ttext_nb->DrawLatex(17. , ymax_top/50. , "N_{b} = 1");

ttext_nb->DrawLatex(28. , ymax_top/50. , "N_{b} = 2");

ttext_nb->DrawLatex(39. , ymax_top/50. , "N_{b} #geq 3");

TText * ttext = new TLatex(160. , ymax_top/50. , "Normalized to 10 fb^{-1}");

ttext->SetTextFont(42);

ttext->SetTextSize(0.045);

ttext->SetTextAlign(22);

ttext->Draw();

}

// Legend & texts

sprintf(tempname,"Lost-lepton background");

catLeg1->SetHeader(tempname);

//sprintf(tempname,"#tau_{hadronic} BG expectation (MC truth)");

sprintf(tempname,"Direct from simulation");

catLeg1->AddEntry(GenHist,tempname,"p");

//sprintf(tempname,"Prediction from MC");

if(doDataVsMC) sprintf(tempname,"Prediction from data");

else sprintf(tempname,"Treat simulation like data");

catLeg1->AddEntry(EstHist,tempname);

catLeg1->Draw();

gPad->RedrawAxis();

//

// Bottom ratio plot

//

// ----------

//

// Preparing ratio histograms

TH1D * numerator = static_cast<TH1D*>(GenHist->Clone("numerator"));

TH1D * numerator_fullstaterr = static_cast<TH1D*>(GenHist->Clone("numerator_fullstaterr"));

TH1D * denominator = static_cast<TH1D*>(EstHist->Clone("denominator"));

TH1D * GenHist_Clone = static_cast<TH1D*>(GenHist->Clone("GenHist_Clone"));

TH1D * EstHist_Clone = static_cast<TH1D*>(EstHist->Clone("EstHist_Clone"));

TH1D * EstHist_NoError = static_cast<TH1D*>(EstHist->Clone("EstHist_NoError"));

TH1D * One_NoError = static_cast<TH1D*>(EstHist->Clone("EstHist_NoError"));

for (int ibin=0; ibin<EstHist_NoError->GetNbinsX()+2; ibin++){ // scan including underflow and overflow bins

EstHist_NoError->SetBinError(ibin,0.);

One_NoError->SetBinContent(ibin,1.);

One_NoError->SetBinError(ibin,0.);

}

//EstHistD->Add(GenHistD,-1);

numerator->Divide(GenHist_Clone,EstHist_NoError,1,1,"");

denominator->Divide(EstHist_Clone,EstHist_NoError,1,1,"");

numerator_fullstaterr->Divide(GenHist_Clone,EstHist_Clone,1,1,""); // Expectation/Prediction

numerator_fullstaterr->Add(One_NoError,-1.); // Expectation/Prediction-1

// draw bottom figure

canvas_dw->cd();

// font size

numerator->GetXaxis()->SetLabelSize(font_size_dw);

numerator->GetXaxis()->SetTitleSize(font_size_dw);

numerator->GetYaxis()->SetLabelSize(font_size_dw);

numerator->GetYaxis()->SetTitleSize(font_size_dw);

//

// Horizontal Lines

TLine *tline = new TLine(search_x_min,1.,search_x_max,1.);

TLine *tline0 = new TLine(search_x_min,0.,search_x_max,0.);

//

// Common to all bottom plots

//

//sprintf(ytitlename,"#frac{Estimate - #tau_{had} BG}{#tau_{had} BG} ");

sprintf(ytitlename,"#frac{Direct}{Prediction} ");

numerator->SetMaximum(ymax_bottom);

numerator->SetMinimum(ymin_bottom);

//

// Specific to each bottom plot

//

// Setting style

//numerator->SetMaximum(1.4);

//numerator->GetXaxis()->SetLabelFont(42);

//numerator->GetXaxis()->SetLabelOffset(0.007);

numerator->GetXaxis()->SetLabelSize(0.18*0.045/0.06);

numerator->GetXaxis()->SetTitleSize(0.18);

numerator->GetXaxis()->SetTitleOffset(0.9);

numerator->GetXaxis()->SetTitleFont(42);

//numerator->GetYaxis()->SetLabelFont(42);

//numerator->GetYaxis()->SetLabelOffset(0.007);

numerator->GetYaxis()->SetLabelSize(0.18*0.045/0.06);

numerator->GetYaxis()->SetTitleSize(0.18);

//numerator->GetYaxis()->SetTitleOffset(0.5);

numerator->GetYaxis()->SetTitleOffset(0.25);

numerator->GetYaxis()->SetTitleFont(42);

numerator->GetXaxis()->SetTitle(xtitlename);

numerator->GetYaxis()->SetTitle(ytitlename);

//gPad->SetGridx(1);

if (pull==1){

sprintf(ytitlename,"#frac{Exp - Pre}{Stat Error} ");

numerator->SetMaximum(8.);

numerator->SetMinimum(-8.);

//

// Specific to each bottom plot

//

// Setting style

for (int ibin=0; ibin<numerator_fullstaterr->GetNbinsX()+2; ibin++){ // scan including underflow and overflow bins

numerator_fullstaterr->SetBinContent(ibin,numerator_fullstaterr->GetBinContent(ibin)/numerator_fullstaterr->GetBinError(ibin));

numerator_fullstaterr->SetBinError(ibin,0.);

}

numerator_fullstaterr->Print("all");

numerator_fullstaterr->GetXaxis()->SetLabelSize(font_size_dw);

numerator_fullstaterr->GetXaxis()->SetTitleSize(font_size_dw);

numerator_fullstaterr->GetYaxis()->SetLabelSize(font_size_dw);

numerator_fullstaterr->GetYaxis()->SetTitleSize(font_size_dw);

numerator_fullstaterr->GetXaxis()->SetTitleSize(0.12);

numerator_fullstaterr->GetXaxis()->SetTitleOffset(0.9);

numerator_fullstaterr->GetXaxis()->SetTitleFont(42);

numerator_fullstaterr->GetYaxis()->SetTitleSize(0.13);

numerator_fullstaterr->GetYaxis()->SetTitleOffset(0.5);

numerator_fullstaterr->GetYaxis()->SetTitleFont(42);

numerator_fullstaterr->GetXaxis()->SetTitle(xtitlename);

numerator_fullstaterr->GetYaxis()->SetTitle(ytitlename);

//numerator_fullstaterr->SetFillColor(kGreen-3);

numerator_fullstaterr->SetFillColor(kRed-10);

numerator_fullstaterr->DrawCopy();

//

// Drawing lines

tline0->SetLineStyle(2);

//tline0->Draw();

}

else {

//

// Plotting

numerator->GetYaxis()->SetNdivisions(505);

numerator->GetYaxis()->SetTickLength(0.015);

numerator->GetXaxis()->SetTickLength(0.08);

numerator->SetTitle("");

ex1->Draw();

numerator->DrawCopy();

ex2->Draw();

denominator->DrawCopy("e2same");

//denominator->DrawCopy("same");

TH1D *denominator_Clone = (TH1D*)denominator->Clone();

denominator_Clone->SetFillColor(kWhite);

denominator_Clone->Draw("hist same");

ex1->Draw();

numerator->DrawCopy("same");

numerator->Print("all");

denominator->Print("all");

numerator_fullstaterr->Print("all");

//

// Drawing lines

tline->SetLineStyle(2);

//tline->Draw();

}

//

if(option.find("QCD")==string::npos ){

// Njet separation lines

TLine *tl_njet = new TLine();

tl_njet->SetLineStyle(2);

tl_njet->DrawLine(53.-0.5,ymin_bottom,53.-0.5,ymax_bottom);

tl_njet->DrawLine(105.-0.5,ymin_bottom,105.-0.5,ymax_bottom);

tl_njet->DrawLine(157.-0.5,ymin_bottom,157.-0.5,ymax_bottom);

// Nb separation lines

TLine *tl_nb = new TLine();

tl_nb->SetLineStyle(3);

tl_nb->DrawLine(14.-0.5,ymin_bottom,14.-0.5,ymax2_bottom);

tl_nb->DrawLine(27.-0.5,ymin_bottom,27.-0.5,ymax2_bottom);

tl_nb->DrawLine(40.-0.5,ymin_bottom,40.-0.5,ymax2_bottom);

tl_nb->DrawLine(53.-0.5,ymin_bottom,53.-0.5,ymax2_bottom);

tl_nb->DrawLine(66.-0.5,ymin_bottom,66.-0.5,ymax2_bottom);

tl_nb->DrawLine(79.-0.5,ymin_bottom,79.-0.5,ymax2_bottom);

tl_nb->DrawLine(92.-0.5,ymin_bottom,92.-0.5,ymax2_bottom);

tl_nb->DrawLine(105.-0.5,ymin_bottom,105.-0.5,ymax2_bottom);

tl_nb->DrawLine(118.-0.5,ymin_bottom,118.-0.5,ymax2_bottom);

tl_nb->DrawLine(131.-0.5,ymin_bottom,131.-0.5,ymax2_bottom);

tl_nb->DrawLine(144.-0.5,ymin_bottom,144.-0.5,ymax2_bottom);

tl_nb->DrawLine(157.-0.5,ymin_bottom,157.-0.5,ymax2_bottom);

tl_nb->DrawLine(170.-0.5,ymin_bottom,170.-0.5,ymax2_bottom);

tl_nb->DrawLine(183.-0.5,ymin_bottom,183.-0.5,ymax2_bottom);

tl_nb->DrawLine(196.-0.5,ymin_bottom,196.-0.5,ymax2_bottom);

} else {

// Njet separation lines

TLine *tl_njet = new TLine();

tl_njet->SetLineStyle(2);

tl_njet->DrawLine( 45.,ymin_bottom, 45.,ymax_bottom);

tl_njet->DrawLine( 89.,ymin_bottom, 89.,ymax_bottom);

tl_njet->DrawLine(133.,ymin_bottom,133.,ymax_bottom);

tl_njet->DrawLine(177.,ymin_bottom,177.,ymax_bottom);

// Nb separation lines

TLine *tl_nb = new TLine();

tl_nb->SetLineStyle(3);

tl_nb->DrawLine(12.,ymin_bottom,12.,ymax2_bottom);

tl_nb->DrawLine(23.,ymin_bottom,23.,ymax2_bottom);

tl_nb->DrawLine(34.,ymin_bottom,34.,ymax2_bottom);

tl_nb->DrawLine(56.,ymin_bottom,56.,ymax2_bottom);

tl_nb->DrawLine(67.,ymin_bottom,67.,ymax2_bottom);

tl_nb->DrawLine(78.,ymin_bottom,78.,ymax2_bottom);

tl_nb->DrawLine(100.,ymin_bottom,100.,ymax2_bottom);

tl_nb->DrawLine(111.,ymin_bottom,111.,ymax2_bottom);

tl_nb->DrawLine(122.,ymin_bottom,122.,ymax2_bottom);

tl_nb->DrawLine(144.,ymin_bottom,144.,ymax2_bottom);

tl_nb->DrawLine(155.,ymin_bottom,155.,ymax2_bottom);

tl_nb->DrawLine(166.,ymin_bottom,166.,ymax2_bottom);

tl_nb->DrawLine(188.,ymin_bottom,188.,ymax2_bottom);

tl_nb->DrawLine(199.,ymin_bottom,199.,ymax2_bottom);

tl_nb->DrawLine(210.,ymin_bottom,210.,ymax2_bottom);

}

gPad->RedrawAxis();

//

//

CMS_lumi(canvas, iPeriod, iPos, lumi_sqrtS);

if(doDataVsMC){

sprintf(tempname,"DataMC_Full_Plot.pdf");

if (pull==1) sprintf(tempname,"DataMCPull_Full_Plot.pdf");

}else{

if(doClosurewoIsoTrackVeto){

if(option.find("QCD")!=string::npos) sprintf(tempname,"Closure_QCD_HDP_woIsoTrack_Full_Plot.pdf");

else sprintf(tempname,"Closure_woIsoTrack_Full_Plot.pdf");

if (pull==1) sprintf(tempname,"ClosurePull_woIsoTrack_Full_Plot.pdf");

}else{

if(option.find("QCD")!=string::npos) sprintf(tempname,"Closure_QCD_HDP_Full_Plot.pdf");

else sprintf(tempname,"Closure_Full_Plot.pdf");

if (pull==1) sprintf(tempname,"ClosurePull_Full_Plot.pdf");

}

}

canvas->Print(tempname);

}