void PostFitCombine(TString Plots = "fit_s", TString InpDir = "FitResults_DataCardFixbtagSysVisPhSp_hSF-PreApp-v0_Tree_LepJets_NewJEC-OldKinFit_v8-0-6_Spring16-80X_36814pb-1_2btag", TString FitDir = "OBSERVED"){
TString PlotsFileName;
if(FitDir == "OBSERVED") PlotsFileName = "obs";
if(FitDir == "EXPECTED") PlotsFileName = "exp";
TString inputfile = "CombineResults/" + InpDir + "/" + FitDir + "/fitDiagnostics" + PlotsFileName + "MLF.root";
setTDRStyle();
gROOT->SetStyle("Plain");
gStyle->SetOptFit(1000);
gStyle->SetOptStat("emruo");
gStyle->SetOptStat(kFALSE);
gStyle->SetPadTickY(1);
gStyle->SetPadTickX(1);
int col_ttbb = TColor::GetColor("#660000");
int col_ttb = TColor::GetColor("#ffcc00");
int col_ttcc = TColor::GetColor("#cc6600");
int col_ttc = TColor::GetColor("#cc6600");
int col_ttLF = TColor::GetColor("#ff0000");
int col_tt = TColor::GetColor("#FF7F7F");
int col_ttbarBkg = TColor::GetColor("#ff6565");
int col_SingleTop = TColor::GetColor("#ff00ff");
int col_WJets = TColor::GetColor("#33cc33");
int col_ZJets = TColor::GetColor("#3366ff");
int col_QCD = TColor::GetColor("#ffff00");
int col_ttbarV = TColor::GetColor("#e75c8d");
int col_ttbarH = TColor::GetColor("#e5c4f4");
int col_VV = TColor::GetColor("#ffffff");
TFile *hfile = NULL;
hfile = TFile::Open(inputfile);
cout << "file loaded: " << inputfile << endl;
std::vector<TString> hNamefile;
std::vector<int> hColor;
hNamefile.push_back("ttbar_LepJetsPowhegPythiattbb"); // 0
hColor.push_back(TColor::GetColor("#660000"));
hNamefile.push_back("ttbar_LepJetsPowhegPythiattbj"); // 1
hColor.push_back(TColor::GetColor("#ffcc00"));
hNamefile.push_back("ttbar_LepJetsPowhegPythiattcc"); // 2
hColor.push_back(TColor::GetColor("#cc6600"));
hNamefile.push_back("ttbar_LepJetsPowhegPythiattLF"); // 3
hColor.push_back(TColor::GetColor("#ff0000"));
hNamefile.push_back("ttbar_PowhegPythiaBkgtt");
hColor.push_back(TColor::GetColor("#FF7F7F"));
hNamefile.push_back("ttHbb_PowhegPythia");
hColor.push_back(TColor::GetColor("#e5c4f4"));
hNamefile.push_back("ttV_Madgraph");
hColor.push_back(TColor::GetColor("#e75c8d"));
hNamefile.push_back("WJets_aMCatNLO");
hColor.push_back(TColor::GetColor("#33cc33"));
hNamefile.push_back("QCD");
hColor.push_back(TColor::GetColor("#ffff00"));
hNamefile.push_back("SingleTop");
hColor.push_back(TColor::GetColor("#ff00ff"));
hNamefile.push_back("VV");
hColor.push_back(TColor::GetColor("#ffffff"));
hNamefile.push_back("ZJets_aMCatNLO");
hColor.push_back(TColor::GetColor("#3366ff"));
hNamefile.push_back("total");
hColor.push_back(1);
hNamefile.push_back("total_signal");
hColor.push_back(1);
hNamefile.push_back("total_background");
hColor.push_back(1);
//TString dirname[2] = {"Name1","Name2"}; //{mu,e}
TString dirname[2] = {"ch1","ch2"}; //{mu,e}
TString titlechname[2] = {"#mu+Jets","e+Jets"};
TString chname[2] = {"mujets","ejets"};
std::vector<TH1D*> hInput[2];
TH1D *hData[2], *hData_reg[2][20];
THStack *AllMC[2], *AllMC_reg[2][20];
THStack *AllMC_CSV1[2], *AllMC_CSV2[2];
// THStack Initialization
for(int ich=0;ich<2;ich++){
AllMC[ich] = new THStack("PostFit_"+chname[ich], "CSV Distribution Post-Fit ("+titlechname[ich]+")");
for(int ireg=0;ireg<20;ireg++){
TString RegNum;
RegNum.Form("%i",ireg);
AllMC_reg[ich][ireg] = new THStack("PostFit_"+chname[ich]+RegNum, "CSV Distribution Post-Fit ("+titlechname[ich]+") for "+RegNum);
} // for(ireg)
} // for(ich)
cout << "Loading histograms... " << endl;
// MC Profiles
for(int ich=0;ich<2;ich++){
for(int ih=0;ih<hNamefile.size();ih++){
cout << "shapes_"+Plots+"/"+dirname[ich]+"/" + hNamefile.at(ih) << endl;
TH1D *htemp = (TH1D*) hfile->Get("shapes_"+Plots+"/"+dirname[ich]+"/" + hNamefile.at(ih))->Clone("c" + hNamefile.at(ih));
htemp->SetFillColor(hColor.at(ih));
htemp->SetLineColor(1);
hInput[ich].push_back(htemp);
cout << hNamefile.at(ih) << " = " << htemp->Integral() << endl;
if(ih==0) AllMC[ich]->SetHistogram((TH1D*)htemp->Clone("FirstStack"));
if (!hNamefile.at(ih).Contains("total")) AllMC[ich] -> Add(htemp);
// Data
// Clone Histo Structure only once
if(ih==0){
hData[ich] = (TH1D *) htemp->Clone("data_"+dirname[ich]);
hData[ich]->Reset();
}
} // for(ih)
// Data TGraph
TGraphAsymmErrors *DataFull = (TGraphAsymmErrors *) hfile->Get("shapes_"+Plots+"/"+dirname[ich]+"/data")->Clone("data_"+dirname[ich]);
for(int ibin=1;ibin<=hData[ich]->GetNbinsX();ibin++){
double igb, EvtBinData, EvtErrBinData;
DataFull->GetPoint((ibin-1),igb,EvtBinData);
hData[ich]->SetBinContent(ibin,EvtBinData);
EvtErrBinData = DataFull->GetErrorY((ibin-1));
hData[ich]->SetBinError(ibin,EvtErrBinData);
}
// Histograms for each region
for(int ih=0;ih<hNamefile.size();ih++){
TH1D *htemp = (TH1D*) hfile->Get("shapes_"+Plots+"/"+dirname[ich]+"/" + hNamefile.at(ih))->Clone("c_reg" + hNamefile.at(ih));
for(int ireg=0;ireg<20;ireg++){
TH1D *htempreg = new TH1D ("","",20,0,20);
for(int ibin=1;ibin<=20;ibin++){
//if ( (ibin+20*ireg) == 181) cout << hNamefile.at(ih) << " " << ibin+20*ireg << " : " << htemp->GetBinContent(ibin+20*ireg) << " ; " << htemp->GetBinError(ibin+20*ireg) << endl;
htempreg->SetBinContent(ibin,htemp->GetBinContent(ibin+20*ireg));
htempreg->SetBinError(ibin,htemp->GetBinError(ibin+20*ireg));
htempreg->SetFillColor(hColor.at(ih));
} // for(ibin)
if(ih==0) AllMC_reg[ich][ireg]->SetHistogram((TH1D*)htempreg->Clone("FirstStack"));
if (!hNamefile.at(ih).Contains("total")) AllMC_reg[ich][ireg] -> Add(htempreg);
} // for(ireg)
} // for(ih)
// Data
for(int ireg=0;ireg<20;ireg++){
TString RegNum;
RegNum.Form("%i",ireg);
hData_reg[ich][ireg] = new TH1D ("hData_"+chname[ich]+RegNum,"Data Histogram "+titlechname[ich]+" for "+RegNum,20,0,20);
// Data
for(int ibin=1;ibin<=20;ibin++){
double igb, EvtBinData, EvtErrBinData;
DataFull->GetPoint(((ibin-1)+(20*ireg)),igb,EvtBinData);
hData_reg[ich][ireg]->SetBinContent(ibin,EvtBinData);
EvtErrBinData = DataFull->GetErrorY((ibin-1));
hData_reg[ich][ireg]->SetBinError(ibin,EvtErrBinData);
hData_reg[ich][ireg]->SetMarkerStyle(20);
hData_reg[ich][ireg]->SetMarkerSize(0.5);
} // for(ibin)
} // for(ireg)
// Recover basic CSV plots
// -- MC
TH1D *GlobalCSVJet[12][2]; // 12 components
AllMC_CSV1[ich] = new THStack("CSV1_"+chname[ich], "CSV-AddJet1 Distribution Post-Fit ("+titlechname[ich]+")");
AllMC_CSV2[ich] = new THStack("CSV2_"+chname[ich], "CSV-AddJet2 Distribution Post-Fit ("+titlechname[ich]+")");
for(int ih=0;ih<12;ih++){
GlobalCSVJet[ih][0] = new TH1D("GlobalCSVJet1_" + hNamefile.at(ih), "CSV distribution for AddJet-1", 20, 0.0, 1.0);
GlobalCSVJet[ih][1] = new TH1D("GlobalCSVJet2_" + hNamefile.at(ih), "CSV distribution for AddJet-2", 20, 0.0, 1.0);
RecoverCSVHisto (hInput[ich].at(ih), GlobalCSVJet[ih][0], GlobalCSVJet[ih][1]);
GlobalCSVJet[ih][0]->SetFillColor(hColor.at(ih));
GlobalCSVJet[ih][1]->SetFillColor(hColor.at(ih));
if(ih==0) AllMC_CSV1[ich]->SetHistogram((TH1D*)GlobalCSVJet[ih][0]->Clone("FirstStack"));
if(ih==0) AllMC_CSV2[ich]->SetHistogram((TH1D*)GlobalCSVJet[ih][1]->Clone("FirstStack"));
AllMC_CSV1[ich] -> Add(GlobalCSVJet[ih][0]);
AllMC_CSV2[ich] -> Add(GlobalCSVJet[ih][1]);
}
// -- Data
TH1D *GlobalCSVJet_Data[2]; // 12 components
GlobalCSVJet_Data[0] = new TH1D("GlobalCSVJet1_Data", "CSV distribution for AddJet-1", 20, 0.0, 1.0);
GlobalCSVJet_Data[1] = new TH1D("GlobalCSVJet2_Data", "CSV distribution for AddJet-2", 20, 0.0, 1.0);
RecoverCSVHisto (hData[ich], GlobalCSVJet_Data[0], GlobalCSVJet_Data[1]);
// -----------
// Plotting
// -----------
TH1D *hstyle = new TH1D ("","",
hData[ich]->GetNbinsX(),
hData[ich]->GetBinLowEdge (1),
hData[ich]->GetBinLowEdge (hData[ich]->GetNbinsX()+1));
hstyle -> SetMaximum(1.1*hInput[ich].at(12)->GetMaximum());
hstyle -> GetYaxis()->SetTitleFont(42);
hstyle -> GetYaxis()->SetTitleOffset(0.7);
hstyle -> GetYaxis()->SetTitleSize(0.05);
hstyle -> GetYaxis()->SetLabelFont(42);
hstyle -> GetYaxis()->SetLabelSize(0.045);
hstyle -> GetYaxis()->SetNdivisions(607);
hstyle -> GetYaxis()->SetTitle("Events / unit");
hData[ich] -> SetMarkerStyle(20);
hData[ich] -> SetMarkerSize(0.4);
hData[ich] -> SetLineWidth(1);
hData[ich] -> SetTitle("");
TCanvas *cPlots;//histos
cPlots = new TCanvas("cPlots"+dirname[ich] ,"Plots");
cPlots->Divide(1,2);
TPad *pad[4], *glpad[2];
// Global Pad
glpad[0] = (TPad*)cPlots->GetPad(1);
glpad[0]->Divide(1,2);
glpad[1] = (TPad*)cPlots->GetPad(2);
glpad[1]->Divide(1,2);
//Plot Pad
pad[0] = (TPad*)glpad[0]->GetPad(1);
pad[0]->SetPad(0.01, 0.23, 0.99, 0.99);
pad[0]->SetTopMargin(0.1);
pad[0]->SetRightMargin(0.04);
//Ratio Pad
pad[1] = (TPad*)glpad[0]->GetPad(2);
pad[1]->SetPad(0.01, 0.02, 0.99, 0.3);
gStyle->SetGridWidth(1);
gStyle->SetGridColor(14);
pad[1]->SetGridx();
pad[1]->SetGridy();
pad[1]->SetTopMargin(0.05);
pad[1]->SetBottomMargin(0.4);
pad[1]->SetRightMargin(0.04);
//Plot Pad
pad[2] = (TPad*)glpad[1]->GetPad(1);
pad[2]->SetPad(0.01, 0.23, 0.99, 0.99);
pad[2]->SetTopMargin(0.1);
pad[2]->SetRightMargin(0.04);
//Ratio Pad
pad[3] = (TPad*)glpad[1]->GetPad(2);
pad[3]->SetPad(0.01, 0.02, 0.99, 0.3);
gStyle->SetGridWidth(1);
gStyle->SetGridColor(14);
pad[3]->SetGridx();
pad[3]->SetGridy();
pad[3]->SetTopMargin(0.05);
pad[3]->SetBottomMargin(0.4);
pad[3]->SetRightMargin(0.04);
pad[0]->cd();
hstyle->Draw();
AllMC[ich] -> Draw("HISTSAME");
hData[ich] -> Draw("PSAME");
TH1D *RatioFull = HistoRatio (hData[ich] , (TH1D*) AllMC[ich]->GetStack()->Last());
TGraphErrors *gRatioFull = new TGraphErrors(RatioFull);
gRatioFull->SetFillStyle(1001);
gRatioFull->SetFillColor(chatch);
gRatioFull->SetName("gRatioFull");
TLegend *leg;
float legPos[4] = {0.70, // x_o
0.40, // y_o
0.94, // x_f
0.87}; // y_f
leg = new TLegend(legPos[0],legPos[1],legPos[2],legPos[3]);
leg->SetFillColor(0);
leg->SetLineColor(0);
leg->SetLineWidth(0.0);
leg->SetTextFont(62);
leg->SetTextSize(0.03);
leg->SetNColumns(2);
leg->AddEntry(hData[ich], "Data","PL");
leg->AddEntry(hInput[ich].at(11), "Z+Jets","F");
leg->AddEntry(hInput[ich].at(10), "VV","F");
leg->AddEntry(hInput[ich].at(9), "Single t","F");
leg->AddEntry(hInput[ich].at(8), "QCD","F");
leg->AddEntry(hInput[ich].at(7), "W+Jets","F");
leg->AddEntry(hInput[ich].at(6), "t#bar{t}+V","F");
leg->AddEntry(hInput[ich].at(5), "t#bar{t}+H","F");
leg->AddEntry(hInput[ich].at(4), "t#bar{t}+other","F");
leg->AddEntry(hInput[ich].at(3), "t#bar{t}+LF","F");
leg->AddEntry(hInput[ich].at(2), "t#bar{t}+cc","F");
leg->AddEntry(hInput[ich].at(1), "t#bar{t}+bj","F");
leg->AddEntry(hInput[ich].at(0), "t#bar{t}+bb","F");
leg->AddEntry(gRatioFull, "Stat. Unc.","F");
leg->Draw("SAME");
TLatex *titlePr;
titlePr = new TLatex(-20.,50.,"35.9 fb^{-1} (13TeV)");
titlePr->SetNDC();
titlePr->SetTextAlign(12);
titlePr->SetX(0.78);
titlePr->SetY(0.935);
titlePr->SetTextFont(42);
titlePr->SetTextSize(0.05);
titlePr->SetTextSizePixels(24);
titlePr->Draw("SAME");
TLatex *title;
//title = new TLatex(-20.,50.,"CMS(2016) #sqrt{s} = 13TeV, L = 35.9 fb^{-1}");
title = new TLatex(-20.,50.,"CMS");
title->SetNDC();
title->SetTextAlign(12);
title->SetX(0.13);
title->SetY(0.84);
title->SetTextFont(61);
title->SetTextSize(0.06);
title->SetTextSizePixels(24);
title->Draw("SAME");
TLatex *chtitle;
chtitle = new TLatex(-20.,50.,titlechname[ich]+"");
chtitle->SetNDC();
chtitle->SetTextAlign(12);
chtitle->SetX(0.14);
chtitle->SetY(0.74);
chtitle->SetTextFont(42);
chtitle->SetTextSize(0.05);
chtitle->SetTextSizePixels(24);
chtitle->Draw("SAME");
pad[2]->cd();
pad[2]->cd()->SetLogy();
TH1D *hstyleLog = (TH1D *)hstyle->Clone();
hstyleLog -> SetMaximum(10.0*hInput[ich].at(12)->GetMaximum());
hstyleLog -> SetMinimum(0.8);
hstyleLog -> Draw();
AllMC[ich] -> Draw("HISTSAME");
hData[ich] -> Draw("PSAME");
titlePr->Draw("SAME");
title->Draw("SAME");
chtitle->Draw("SAME");
pad[1]->cd();
RatioFull ->Draw("HIST");
gRatioFull->Draw("e2");
RatioFull ->Draw("HISTSAME");
pad[3]->cd();
RatioFull ->Draw("HIST");
gRatioFull->Draw("e2");
RatioFull ->Draw("HISTSAME");
TString dirfigname_pdf = "CombineResults/Figures_" + InpDir + Plots + FitDir + "/";
// make a dir if it does not exist!!
gSystem->mkdir(dirfigname_pdf, kTRUE);
cPlots->SaveAs(dirfigname_pdf + "FullHisto_" + dirname[ich] + "_NormLog.pdf");
// -----------------------------------------------------------------------------
// Only Log Plots
// -----------------------------------------------------------------------------
TCanvas *cPlotsLog;//histos
cPlotsLog = new TCanvas("cPlotsLog"+dirname[ich] ,"Plots");
cPlotsLog->Divide(1,2);
TPad *padLog[2];
//Plot Pad
padLog[0] = (TPad*)cPlotsLog->GetPad(1);
padLog[0]->SetPad(0.01, 0.23, 0.99, 0.99);
padLog[0]->SetTopMargin(0.1);
padLog[0]->SetRightMargin(0.04);
//Ratio Pad
padLog[1] = (TPad*)cPlotsLog->GetPad(2);;
padLog[1]->SetPad(0.01, 0.02, 0.99, 0.3);
gStyle->SetGridWidth(1);
gStyle->SetGridColor(14);
padLog[1]->SetGridx();
padLog[1]->SetGridy();
padLog[1]->SetTopMargin(0.05);
padLog[1]->SetBottomMargin(0.4);
padLog[1]->SetRightMargin(0.04);
padLog[0]->cd();
padLog[0]->cd()->SetLogy();
hstyleLog -> Draw();
AllMC[ich] -> Draw("HISTSAME");
hData[ich] -> Draw("PSAME");
titlePr->Draw("SAME");
title->Draw("SAME");
chtitle->Draw("SAME");
TLegend *legLog;
legLog = new TLegend(0.70,0.60,0.94,0.87);
legLog->SetFillColor(0);
legLog->SetLineColor(0);
legLog->SetLineWidth(0.0);
legLog->SetTextFont(62);
legLog->SetTextSize(0.03);
legLog->SetNColumns(2);
legLog->AddEntry(hData[ich], "Data","PL");
legLog->AddEntry(hInput[ich].at(11), "Z+Jets","F");
legLog->AddEntry(hInput[ich].at(10), "VV","F");
legLog->AddEntry(hInput[ich].at(9), "Single t","F");
legLog->AddEntry(hInput[ich].at(8), "QCD","F");
legLog->AddEntry(hInput[ich].at(7), "W+Jets","F");
legLog->AddEntry(hInput[ich].at(6), "t#bar{t}+V","F");
legLog->AddEntry(hInput[ich].at(5), "t#bar{t}+H","F");
legLog->AddEntry(hInput[ich].at(4), "t#bar{t}+other","F");
legLog->AddEntry(hInput[ich].at(3), "t#bar{t}+LF","F");
legLog->AddEntry(hInput[ich].at(2), "t#bar{t}+cc","F");
legLog->AddEntry(hInput[ich].at(1), "t#bar{t}+bj","F");
legLog->AddEntry(hInput[ich].at(0), "t#bar{t}+bb","F");
legLog->AddEntry(gRatioFull, "Stat. Unc.","F");
legLog->Draw("SAME");
padLog[1]->cd();
RatioFull ->Draw("HIST");
gRatioFull->Draw("e2");
RatioFull ->Draw("HISTSAME");
cPlotsLog->SaveAs(dirfigname_pdf + "FullHisto_" + dirname[ich] + "_OnlyLog.pdf");
// -----------------------------------------------------------------------------
// -----------------------------------------------------------------------------
// -----------------------------------------------------------------------------
// Plots by Regions
// -----------------------------------------------------------------------------
// -----------------------------------------------------------------------------
// -----------------------------------------------------------------------------
int ireg = 0;
for(int ican=0;ican<5;ican++){
TString CanNum;
CanNum.Form("%i",ican);
TCanvas *cPlots_reg;//histos
cPlots_reg = new TCanvas("cPlots_reg" + CanNum + dirname[ich] ,"Plots By regions");
cPlots_reg->Divide(2,2);
for(int icr=1;icr<=4;icr++){
//cPlots_reg->cd(icr);
TPad *glpad_reg = (TPad*)cPlots_reg->cd(icr);
glpad_reg->Divide(2,1);
TPad *pad_reg[2];
//Plot Pad
pad_reg[0] = (TPad*)glpad_reg->GetPad(1);
pad_reg[0]->SetPad(0.01, 0.23, 0.99, 0.99);
pad_reg[0]->SetTopMargin(0.1);
pad_reg[0]->SetRightMargin(0.04);
//Ratio Pad
pad_reg[1] = (TPad*)glpad_reg->GetPad(2);
pad_reg[1]->SetPad(0.01, 0.02, 0.99, 0.3);
gStyle->SetGridWidth(1);
gStyle->SetGridColor(14);
pad_reg[1]->SetGridx();
pad_reg[1]->SetGridy();
pad_reg[1]->SetTopMargin(0.05);
pad_reg[1]->SetBottomMargin(0.4);
pad_reg[1]->SetRightMargin(0.04);
pad_reg[0]->cd();
pad_reg[0]->SetLogy();
TH1D *hstyle_reg = (TH1D*)AllMC_reg[ich][ireg]->GetHistogram();
hstyle_reg -> Reset();
hstyle_reg -> SetMaximum(10.0*hData_reg[ich][ireg]->GetMaximum());
hstyle_reg -> SetMinimum(0.7);
hstyle_reg -> GetYaxis()->SetTitleOffset(0.9);
hstyle_reg -> GetYaxis()->SetTitleSize(0.05);
hstyle_reg -> GetYaxis()->SetLabelSize(0.05);
hstyle_reg -> GetYaxis()->SetTitle("Events");
pad_reg[0]->cd();
hstyle_reg->Draw();
AllMC_reg[ich][ireg] -> Draw("HISTSAME");
hData_reg[ich][ireg] -> Draw("E1SAME");
titlePr->Draw("SAME");
title->Draw("SAME");
chtitle->Draw("SAME");
TH1D *RatioFull_reg = HistoRatio (hData_reg[ich][ireg] , (TH1D*) AllMC_reg[ich][ireg]->GetStack()->Last());
TGraphErrors *gRatioFull_reg = new TGraphErrors(RatioFull_reg);
gRatioFull_reg->SetFillStyle(1001);
gRatioFull_reg->SetFillColor(chatch);
gRatioFull_reg->SetName("gRatioFull");
pad_reg[1]->cd();
RatioFull_reg->GetYaxis()->SetTitleOffset(0.25);
RatioFull_reg->Draw("HIST");
gRatioFull_reg->Draw("e2");
RatioFull_reg->Draw("HISTSAME");
ireg++;
} // for(ireg)
cPlots_reg->SaveAs(dirfigname_pdf + "RegionHisto_" + CanNum + "_" + dirname[ich] + "_NormLog.pdf");
} // for(ican)
TCanvas *cPlotsI;
cPlotsI = new TCanvas("cPlotsI"+dirname[ich] ,"Plots");
cPlotsI->Divide(1,2);
TPad *padI[4];
//Plot Pad
padI[0] = (TPad*)cPlotsI->GetPad(1);
padI[0]->SetPad(0.01, 0.23, 0.99, 0.99);
padI[0]->SetTopMargin(0.1);
padI[0]->SetRightMargin(0.04);
//Ratio Pad
padI[1] = (TPad*)cPlotsI->GetPad(2);
padI[1]->SetPad(0.01, 0.02, 0.99, 0.3);
gStyle->SetGridWidth(1);
gStyle->SetGridColor(14);
padI[1]->SetGridx();
padI[1]->SetGridy();
padI[1]->SetTopMargin(0.05);
padI[1]->SetBottomMargin(0.4);
padI[1]->SetRightMargin(0.04);
TCanvas *cPlotsII;
cPlotsII = new TCanvas("cPlotsII"+dirname[ich] ,"Plots");
cPlotsII->Divide(1,2);
//Plot Pad
padI[2] = (TPad*)cPlotsII->GetPad(1);
padI[2]->SetPad(0.01, 0.23, 0.99, 0.99);
padI[2]->SetTopMargin(0.1);
padI[2]->SetRightMargin(0.04);
//Ratio Pad
padI[3] = (TPad*)cPlotsII->GetPad(2);
padI[3]->SetPad(0.01, 0.02, 0.99, 0.3);
gStyle->SetGridWidth(1);
gStyle->SetGridColor(14);
padI[3]->SetGridx();
padI[3]->SetGridy();
padI[3]->SetTopMargin(0.05);
padI[3]->SetBottomMargin(0.4);
padI[3]->SetRightMargin(0.04);
TH1D *hstyleI = new TH1D ("hstyleI","",
GlobalCSVJet_Data[1]->GetNbinsX(),
GlobalCSVJet_Data[1]->GetBinLowEdge (1),
GlobalCSVJet_Data[1]->GetBinLowEdge (GlobalCSVJet_Data[1]->GetNbinsX()+1));
hstyleI -> SetMaximum(1.1*GlobalCSVJet_Data[1]->GetMaximum());
hstyleI -> GetYaxis()->SetTitleOffset(0.8);
hstyleI -> GetYaxis()->SetTitleSize(0.05);
hstyleI -> GetYaxis()->SetLabelSize(0.05);
hstyleI -> GetYaxis()->SetTitle("Events");
GlobalCSVJet_Data[0] -> SetMarkerStyle(20);
GlobalCSVJet_Data[0] -> SetMarkerSize(0.4);
GlobalCSVJet_Data[0] -> SetLineWidth(1);
GlobalCSVJet_Data[0] -> SetTitle("");
padI[0]->cd();
padI[0]->cd()->SetLogy();
hstyleI -> SetMaximum(100.0*GlobalCSVJet_Data[0]->GetMaximum());
hstyleI -> SetMinimum(0.7);
hstyleI -> Draw();
AllMC_CSV1[ich] -> Draw("HISTSAME");
GlobalCSVJet_Data[0] -> SetMarkerStyle(20);
GlobalCSVJet_Data[0] -> SetMarkerSize(0.6);
GlobalCSVJet_Data[0] -> SetLineWidth(1);
GlobalCSVJet_Data[0] -> SetTitle("");
GlobalCSVJet_Data[0] -> Draw("SAME");
titlePr->Draw("SAME");
title->Draw("SAME");
chtitle->Draw("SAME");
TLegend *legI;
legI = new TLegend(0.70,0.64,0.93,0.87);
legI->SetFillColor(0);
legI->SetLineColor(0);
legI->SetLineWidth(0.0);
legI->SetTextFont(62);
legI->SetTextSize(0.03);
legI->SetNColumns(2);
legI->AddEntry(hData[ich], "Data","PL");
legI->AddEntry(hInput[ich].at(11), "Z+Jets","F");
legI->AddEntry(hInput[ich].at(10), "VV","F");
legI->AddEntry(hInput[ich].at(9), "Single t","F");
legI->AddEntry(hInput[ich].at(8), "QCD","F");
legI->AddEntry(hInput[ich].at(7), "W+Jets","F");
legI->AddEntry(hInput[ich].at(6), "t#bar{t}+V","F");
legI->AddEntry(hInput[ich].at(5), "t#bar{t}+H","F");
legI->AddEntry(hInput[ich].at(4), "t#bar{t}+other","F");
legI->AddEntry(hInput[ich].at(3), "t#bar{t}+LF","F");
legI->AddEntry(hInput[ich].at(2), "t#bar{t}+cc","F");
legI->AddEntry(hInput[ich].at(1), "t#bar{t}+bj","F");
legI->AddEntry(hInput[ich].at(0), "t#bar{t}+bb","F");
legI->AddEntry(gRatioFull, "Stat. Unc.","F");
legI->Draw("SAME");
TH1D *RatioFullJet1 = HistoRatio (GlobalCSVJet_Data[0] , (TH1D*) AllMC_CSV1[ich]->GetStack()->Last());
TGraphErrors *gRatioFullJet1 = new TGraphErrors(RatioFullJet1);
gRatioFullJet1->SetFillStyle(1001);
gRatioFullJet1->SetFillColor(chatch);
gRatioFullJet1->SetName("gRatioFullJet2");
padI[1]->cd();
RatioFullJet1 ->GetYaxis()->SetTitleOffset(0.25);
RatioFullJet1 ->GetXaxis()->SetTitle("CSVv2");
RatioFullJet1 ->Draw("HIST");
gRatioFullJet1->Draw("e2");
RatioFullJet1 ->Draw("HISTSAME");
padI[2]->cd();
padI[2]->cd()->SetLogy();
// No Log
TH1D *hstyleII = (TH1D *)hstyleI -> Clone("StyleII");
hstyleII -> SetMaximum(100.0*GlobalCSVJet_Data[1]->GetMaximum());
hstyleII -> SetMinimum(0.7);
hstyleII -> Draw();
AllMC_CSV2[ich] -> Draw("HISTSAME");
GlobalCSVJet_Data[1] -> SetMarkerStyle(20);
GlobalCSVJet_Data[1] -> SetMarkerSize(0.5);
GlobalCSVJet_Data[1] -> SetLineWidth(1);
GlobalCSVJet_Data[1] -> SetTitle("");
GlobalCSVJet_Data[1] -> Draw("SAME");
titlePr->Draw("SAME");
title->Draw("SAME");
chtitle->Draw("SAME");
legI->Draw("SAME");
TH1D *RatioFullJet2 = HistoRatio (GlobalCSVJet_Data[1] , (TH1D*) AllMC_CSV2[ich]->GetStack()->Last());
TGraphErrors *gRatioFullJet2 = new TGraphErrors(RatioFullJet2);
gRatioFullJet2->SetFillStyle(1001);
gRatioFullJet2->SetFillColor(chatch);
gRatioFullJet2->SetName("gRatioFullJet2");
padI[3]->cd();
RatioFullJet2 ->GetYaxis()->SetTitleOffset(0.25);
RatioFullJet2 ->GetXaxis()->SetTitle("CSVv2");
RatioFullJet2 ->Draw("HIST");
gRatioFullJet2->Draw("e2");
RatioFullJet2 ->Draw("HISTSAME");
cPlotsI ->SaveAs(dirfigname_pdf + "CSVHistosJet1_" + dirname[ich] + "_Log.pdf");
cPlotsII->SaveAs(dirfigname_pdf + "CSVHistosJet2_" + dirname[ich] + "_Log.pdf");
// ttbb Shape
}// for(ich)
}
void mk_sigaccanplots(string flavor = "112", bool tdrstyle = false)
{
setTDRStyle(tdrstyle);
// gStyle->SetOptFit(1100); // chi2 and prob, not parameters
gStyle->SetOptFit(0); // chi2 and prob, not parameters
// gStyle->SetOptStat("irme"); // integral, RMS, mean, # of entries
gStyle->SetStatFontSize(0.005);
gStyle->SetStatY(0.4);
float unused = 0.9, simtxt_y = 0.97;
if (tdrstyle == false) {
simtxt_y = 0.93;
}
TLatex *tex = new TLatex(unused, simtxt_y, "CMS simulation at #sqrt{s} = 8 TeV");
tex->SetNDC();
tex->SetTextAlign(12); // Left-adjusted
tex->SetTextFont(42);
float textsiz = 0.06;
tex->SetTextSize(0.05);
tex->SetLineWidth(2);
// for(int k=0; k<3; k++){
// for(int k=0; k<2; k++){
for(int k=2; k<3; k++){
//Which plots to make
//define input variables
vector <TFile* > filelist_kin;
// string uncert="btagup";
// if(k==1) uncert="btagdown";
string uncert="up";
if(k==1) uncert="down";
if(k==2) uncert="";
// filelist_kin.push_back(TFile::Open(("Acc_Gluino_NoLumi_pt110_8TeVxsec_BtagMap_2500GeV"+uncert+".root").c_str()));
// filelist_kin.push_back(TFile::Open("Acc_Gluino_NoLumi_pt60_pt110_8TeVxsec_BtagMap.root"));
// filelist_kin.push_back(TFile::Open(("Acc_RPV112" + uncert + ".root").c_str()));
filelist_kin.push_back(TFile::Open(("Acc_RPV" + flavor + uncert + "_Sph4.root").c_str()));
// filelist_kin.push_back(TFile::Open(("Acc_Gluino_NoLumi_pt60_pt110_8TeVxsec_BtagMap" + uncert + ".root").c_str()));
TFile f1(("RPV_" + flavor + "_accandwdithtest"+uncert+".root").c_str(), "recreate");
f1.cd();
string nameIN;
string cuts;
string prefix;
string postfix;
string folder;
string ptcut;
folder="plots_paper/";
postfix=".pdf";
/////////////Plots for each Mass separatly//////////////
/////////////----------------------------------------/////////////
for (int p = 0; p < 2; p++) {
// ptcut="112";
// if(p==1) ptcut="" + flavor;
ptcut="60";
if (p==1) {
ptcut="110";
textsiz = 0.045;
}
// string histname = string("GausWidth_vs_Mass_112") + ptcut;
string histname = string("GausWidth_vs_Mass_" + flavor) + "_" + ptcut;
cout<< histname << endl;
TGraphErrors *h_GluinoHist_Fit = (TGraphErrors*) filelist_kin[0]->Get(histname.c_str())->Clone();
histname = string("GausAcceptance_vs_Mass_" + flavor) + "_" + ptcut;
// histname = string("GausAcceptance_vs_Mass_112") + ptcut;
cout<< histname << endl;
TGraphErrors *h_GluinoHist_MCcomb = (TGraphErrors*) filelist_kin[0]->Get(histname.c_str())->Clone();
// histname = string("FullAcceptance_vs_Mass_" + flavor) + "_" + ptcut;
histname = string("GausMean_vs_Mass_" + flavor) + "_" + ptcut;
// histname = string("GausAcceptance_vs_Mass_112") + ptcut;
cout<< histname << endl;
// TGraphErrors *h_FullAccept = (TGraphErrors*) filelist_kin[0]->Get(histname.c_str())->Clone();
TGraph *h_FullAccept = (TGraph*) filelist_kin[0]->Get(histname.c_str())->Clone();
histname = string("GausMeanOffset_vs_Mass_" + flavor) + "_" + ptcut;
cout<< histname << endl;
TGraph *h_MeanOffset = (TGraph*) filelist_kin[0]->Get(histname.c_str())->Clone();
TCanvas * cGluinoFitsOpti = new TCanvas(("RPV_"+flavor +ptcut).c_str(), ("RPV_" + ptcut+"_"+cuts).c_str(), 800, 600);
//h_GluinoHist_Fit->SetFillColor(kOrange-2);
// h_GluinoHist_Fit->SetLineColor(kBlack);
string title;
string systematic = "pile-up";
// title="Gaussian Width vs. Mass for Light-ptcut RPV";
string tag = "qqb", sphericity = " Sphericity #geq 0.4";
if (flavor.compare("112") == 0)
tag = "qqq";
else if (ptcut.compare("60") == 0)
sphericity = "";
string titlepart = "Hadronic RPV #tilde{g} #rightarrow " + tag;
string titleln2 = "\\Delta = 110 GeV";
string titleln3 = "6^{th}-jet p_{T} #geq " + ptcut + " GeV";
title = titlepart;
// title = "Width for " + titlepart;
/*
if(k==0){
title="RPV gluino #bf{" + systematic + " up} m="+to_string(masses[i])+", ptcut = "+ptcut+", #Delta = 110 GeV, #bf{6^{th} Jet p_{T} = 60 GeV}";
if (i>=5 || p==0) title="RPV gluino #bf{" + systematic + " up} m="+to_string(masses[i])+", ptcut = "+ptcut+", #Delta = 110 GeV, #bf{6^{th} Jet p_{T} = 110 GeV}";
}
if(k==1){
title="RPV gluino #bf{" + systematic + " down} m="+to_string(masses[i])+", ptcut = "+ptcut+", #Delta = 110 GeV, #bf{6^{th} Jet p_{T} = 60 GeV}";
if (i>=5 || p==0) title="RPV gluino #bf{" + systematic + " down} m="+to_string(masses[i])+", ptcut = "+ptcut+", #Delta = 110 GeV, #bf{6^{th} Jet p_{T} = 110 GeV}";
}
if(k==2){
title="RPV gluino m="+to_string(masses[i])+", ptcut = "+ptcut+", #Delta = 110 GeV, #bf{6^{th} Jet p_{T} = 60 GeV}";
if (i>=5 || p==0) title="RPV gluino m="+to_string(masses[i])+", ptcut = "+ptcut+", #Delta = 110 GeV, #bf{6^{th} Jet p_{T} = 110 GeV}";
}
*/
h_GluinoHist_Fit->SetTitle(title.c_str());
float titpos = 0.2, titly = 0.89, headpos = 0.64;
if (tdrstyle == false) {
titpos -= 0.05;
titly -= 0.05;
headpos = 0.45;
}
TLatex *tex2 = new TLatex(titpos, titly, title.c_str());
tex2->SetNDC();
tex2->SetTextAlign(12); // Left-adjusted
tex2->SetTextFont(42);
tex2->SetTextSize(textsiz);
tex2->SetLineWidth(2);
TLatex *tex3 = new TLatex(titpos, titly - 0.07, titleln2.c_str());
tex3->SetNDC();
tex3->SetTextAlign(12);
tex3->SetTextFont(42);
tex3->SetTextSize(textsiz);
tex3->SetLineWidth(2);
TLatex *tex3a = new TLatex(titpos, titly - 0.14, titleln3.c_str());
tex3a->SetNDC();
tex3a->SetTextAlign(12);
tex3a->SetTextFont(42);
tex3a->SetTextSize(textsiz);
tex3a->SetLineWidth(2);
TLatex *tex4 = NULL;
if (sphericity.size() > 0) {
tex4 = new TLatex(titpos - 0.01, titly - 0.22, sphericity.c_str());
tex4->SetNDC();
tex4->SetTextAlign(12);
tex4->SetTextFont(42);
tex4->SetTextSize(textsiz);
tex4->SetLineWidth(2);
}
float legx = 0.65;
if (tdrstyle == false) {
legx = 0.56;
}
TLegend *leg = new TLegend(legx, 0.2, legx + 0.3, 0.4);
leg->SetBorderSize(1);
leg->SetTextFont(62);
leg->SetLineColor(kBlack);
leg->SetLineStyle(1);
leg->SetLineWidth(1);
leg->SetFillColor(0);
leg->SetFillStyle(1001);
//leg->SetHeader();
// leg->AddEntry(h_GluinoHist_Fit, "#splitline{Acceptance as function}{of gluino mass}","l");
h_GluinoHist_Fit->SetMarkerStyle(1);
h_GluinoHist_Fit->SetMarkerColor(kWhite);
// h_GluinoHist_Fit->SetMarkerColor(kGreen + 3);
h_GluinoHist_Fit->SetMarkerSize(0.004);
// h_GluinoHist_Fit->SetTitleSize(0.01);
TF1 *fitfunc = h_GluinoHist_Fit->GetFunction("GausWidth");
string fitnamew = "fitcopy" + ptcut;
TF1* fitfunccopy = (TF1 *) fitfunc->Clone(fitnamew.c_str());
if (fitfunc == NULL)
cout << "Can't get fit func\n";
else {
fitfunc->Delete();
fitfunccopy->SetLineWidth(3);
fitfunccopy->SetLineColor(kGreen + 3);
// fitfunc->SetLineStyle(3); // Dotted
}
float labsiz = 0.055;
h_GluinoHist_Fit->GetXaxis()->SetLabelFont(62);
h_GluinoHist_Fit->GetXaxis()->SetTitleFont(62);
h_GluinoHist_Fit->GetYaxis()->SetLabelFont(62);
h_GluinoHist_Fit->GetYaxis()->SetTitleFont(62);
h_GluinoHist_Fit->GetYaxis()->SetTitle("Gaussian width [GeV]");
float offset = 0.8;
if (tdrstyle == false) {
h_GluinoHist_Fit->GetXaxis()->SetTitleOffset(offset);
h_GluinoHist_Fit->GetYaxis()->SetTitleOffset(offset);
}
h_GluinoHist_Fit->GetXaxis()->SetTitle("Gluino mass [GeV]");
h_GluinoHist_Fit->GetXaxis()->SetTitleSize(labsiz);
h_GluinoHist_Fit->GetYaxis()->SetTitleSize(labsiz);
if (tdrstyle == false) {
float axsize = 0.04;
h_GluinoHist_Fit->GetXaxis()->SetLabelSize(axsize);
h_GluinoHist_Fit->GetYaxis()->SetLabelSize(axsize);
}
if (flavor.compare("113_223") == 0 && ptcut == "60")
h_GluinoHist_Fit->GetYaxis()->SetRangeUser(14.0, 50.0);
h_GluinoHist_Fit->Draw("APX"); // X eliminates error bars
// h_GluinoHist_Fit->Draw("AP");
// TH1 *fithist = (TH1 *) fitfunccopy->GetHistogram()->Clone(fitnamew.c_str());
TH1 *fithist = (TH1 *) fitfunccopy->GetHistogram()->Clone();
int fillcolor = kGreen + 2;
int fillstyle = 3013;
if (fithist != NULL) {
int numBins = fithist->GetNbinsX();
for (int cnt = 1; cnt <= numBins; ++cnt) {
setErr(fithist, cnt, flavor, "width");
}
fithist->SetFillColor(fillcolor);
fithist->SetFillStyle(fillstyle);
fithist->Draw("CE3SAME");
}
fitfunccopy->Draw("CSAME");
// h_GluinoHist_Fit->Draw("P"); // Draw points over fit line
leg->AddEntry(fitfunccopy, "Gaussian width", "L");
leg->AddEntry(fithist, "Uncertainty", "F");
leg->Draw();
tex->SetX(headpos);
tex->Draw();
tex2->Draw();
tex3->Draw();
tex3a->Draw();
if (tex4 != NULL)
tex4->Draw();
cGluinoFitsOpti->Write();
cGluinoFitsOpti->SaveAs((folder + "RPVwidth" +flavor + ptcut+uncert+postfix).c_str());
TCanvas * cGluinoFitsOpt2 = new TCanvas(("RPVacc_"+flavor +ptcut).c_str(), ("RPV_" + ptcut+"_"+cuts).c_str(), 800, 600);
// title="Acc. x Eff. for " + titlepart;
title= titlepart;
tex2->SetText(titpos,titly, title.c_str());
TLegend *leg2 = new TLegend(legx, 0.2, legx + 0.3, 0.4);
leg2->SetBorderSize(1);
leg2->SetTextFont(62);
leg2->SetTextSize(0.04);
leg2->SetLineColor(kBlack);
leg2->SetLineStyle(1);
leg2->SetLineWidth(1);
leg2->SetFillColor(0);
leg2->SetFillStyle(1001);
TF1 *fitfuncA = h_GluinoHist_MCcomb->GetFunction("total");
string fitname = "fitcopyA" + ptcut;
TF1* fitfunccopyA = (TF1 *) fitfuncA->Clone(fitname.c_str());
if (fitfuncA == NULL)
cout << "Can't get fit func\n";
else {
fitfuncA->Delete();
fitfunccopyA->SetLineWidth(3);
fitfunccopyA->SetLineColor(kGreen + 3);
// fitfunc->SetLineStyle(3); // Dotted
}
h_GluinoHist_MCcomb->SetMarkerStyle(1);
// h_GluinoHist_MCcomb->SetMarkerColor(kBlack);
h_GluinoHist_MCcomb->SetMarkerColor(kWhite);
h_GluinoHist_MCcomb->SetTitle(title.c_str());
h_GluinoHist_MCcomb->GetXaxis()->SetLabelFont(62);
h_GluinoHist_MCcomb->GetXaxis()->SetTitleFont(62);
h_GluinoHist_MCcomb->GetYaxis()->SetLabelFont(62);
// h_GluinoHist_MCcomb->GetYaxis()->SetLabelSize(62);
h_GluinoHist_MCcomb->GetYaxis()->SetTitleFont(62);
h_GluinoHist_MCcomb->GetYaxis()->SetTitle("Acceptance x Efficiency");
if (tdrstyle == false) {
h_GluinoHist_MCcomb->GetXaxis()->SetTitleOffset(offset);
h_GluinoHist_MCcomb->GetYaxis()->SetTitleOffset(offset + 0.35);
}
// h_GluinoHist_MCcomb->GetYaxis()->SetTitleOffset(1.4);
h_GluinoHist_MCcomb->GetXaxis()->SetTitle("Gluino mass [GeV]");
h_GluinoHist_MCcomb->GetXaxis()->SetTitleSize(labsiz);
if (tdrstyle == false)
labsiz -= 0.01;
h_GluinoHist_MCcomb->GetYaxis()->SetTitleSize(labsiz);
if (tdrstyle == false) {
float axsize = 0.035;
h_GluinoHist_MCcomb->GetXaxis()->SetLabelSize(axsize);
h_GluinoHist_MCcomb->GetYaxis()->SetLabelSize(axsize);
}
if (flavor.compare("113_223") == 0) {
float ylimit = 0.05;
if ( ptcut == "110")
ylimit = 0.022;
// gStyle->SetStatY(0.8);
h_GluinoHist_MCcomb->GetYaxis()->SetRangeUser(0.0, ylimit);
}
// h_GluinoHist_MCcomb->Draw("AL");
h_GluinoHist_MCcomb->Draw("APX");
// fitfunccopyA->Draw("C same");
// TH1 *fithistA = (TH1 *) fitfunccopyA->GetHistogram()->Clone(fitname.c_str());
TH1 *fithistA = (TH1 *) fitfunccopyA->GetHistogram()->Clone();
if (fithistA != NULL) {
int numBins = fithistA->GetNbinsX();
for (int cnt = 1; cnt <= numBins; ++cnt) {
setErr(fithistA, cnt, flavor, "acceptance");
}
fithistA->SetFillColor(fillcolor);
fithistA->SetFillStyle(fillstyle);
fithistA->Draw("CE3SAME");
}
fitfunccopyA->Draw("CSAME");
// h_GluinoHist_MCcomb->Draw("P"); // Draw points over fit line
tex->SetX(headpos);
tex->Draw();
tex2->Draw();
tex3->Draw();
tex3a->Draw();
leg2->AddEntry(fitfunccopyA, "Acc. x Eff.", "L");
leg2->AddEntry(fithistA, "Uncertainty", "F");
leg2->Draw();
if (tex4 != NULL)
tex4->Draw();
cGluinoFitsOpt2->Write();
cGluinoFitsOpt2->SaveAs((folder + "RPVacc" +flavor + ptcut+uncert+postfix).c_str());
/*
gStyle->SetStatY(0.4);
TCanvas * cGluinoFitsOpt3 = new TCanvas(("RPVfullacc_"+flavor +ptcut).c_str(), ("RPVfull_" + ptcut+"_"+cuts).c_str(), 800, 600);
title="Gaussian Mean for " + titlepart;
tex2->SetText(titpos, 0.89, title.c_str());
// title="Full Acceptance for " + titlepart;
h_FullAccept->SetMarkerStyle(1);
// h_FullAccept->SetMarkerColor(kWhite);
h_FullAccept->SetLineColor(kRed);
h_FullAccept->SetLineWidth(2.0);
// h_FullAccept->SetTitle(title.c_str());
// h_FullAccept->GetYaxis()->SetTitle("Acceptance");
h_FullAccept->GetYaxis()->SetTitle("Gaussian Mean [GeV]");
h_FullAccept->GetXaxis()->SetTitleOffset(1.3);
h_FullAccept->GetXaxis()->SetTitle("Gluino Mass [GeV]");
h_FullAccept->GetXaxis()->SetTitleSize(labsiz);
h_FullAccept->GetYaxis()->SetTitleSize(labsiz);
h_FullAccept->GetXaxis()->SetLabelSize(axsize);
h_FullAccept->GetYaxis()->SetLabelSize(axsize);
h_FullAccept->Draw("AL");
// h_FullAccept->Draw("APX");
// leg->Draw();
// h_GluinoHist_MCcomb->SetFillColor(10);
// h_GluinoHist_MCcomb->SetLineColor(kBlack);
// h_GluinoHist_MCcomb->Draw("samehist");
// f_GluinoGauss->SetLineColor(kRed);
//f_GluinoGauss->Draw("same");
//f_GluinoP4->Draw("same");
tex->SetX(titpos);
tex->Draw();
tex2->Draw();
tex3->Draw();
tex3a->Draw();
if (tex4 != NULL)
tex4->Draw();
f1.cd();
cGluinoFitsOpt3->Write();
cGluinoFitsOpt3->SaveAs((folder + "RPVmean" +flavor + ptcut+uncert+postfix).c_str());
gStyle->SetStatY(0.4);
TCanvas *cMeanOffset = new TCanvas(("RPVMeanOffset_"+flavor +ptcut).c_str(),
("RPVMeanOffset_" + ptcut + "_" + cuts).c_str(), 800, 600);
axsize = 0.035;
title="Mass Deviation for " + titlepart;
titpos = 0.35;
tex2->SetText(titpos, 0.89, title.c_str());
h_MeanOffset->SetMarkerStyle(1);
h_MeanOffset->SetLineColor(kRed);
h_MeanOffset->SetLineWidth(2.0);
h_MeanOffset->GetYaxis()->SetTitleOffset(1.3);
h_MeanOffset->GetYaxis()->SetTitle("Fractional Mass Deviation");
h_MeanOffset->GetXaxis()->SetTitleOffset(1.3);
h_MeanOffset->GetXaxis()->SetTitle("Gluino Mass [GeV]");
h_MeanOffset->GetXaxis()->SetTitleSize(labsiz);
h_MeanOffset->GetYaxis()->SetTitleSize(labsiz);
h_MeanOffset->GetXaxis()->SetLabelSize(axsize);
h_MeanOffset->GetYaxis()->SetLabelSize(axsize);
TLine *max = new TLine(0.5, 0.1, 0.5, 0.9);
max->SetLineColor(kGreen + 2);
max->SetLineWidth(3);
h_MeanOffset->Draw("AL");
max->DrawLineNDC(0.16, 0.49, 0.98, 0.49);
tex->SetX(titpos);
tex->Draw();
tex2->Draw();
tex3->SetX(titpos);
tex3->Draw();
tex3a->SetX(titpos);
tex3a->Draw();
if (tex4 != NULL) {
tex4->SetX(titpos);
tex4->Draw();
}
f1.cd();
cMeanOffset->Write();
cMeanOffset->SaveAs((folder + "RPVmassdev" +flavor + ptcut+uncert+postfix).c_str());
*/
}
}
}
void Plotter(TString plots="2btag", bool LogScale=false, TString cat = "ttbb", TString nSyst = "btagcfI"){
TString files = "../" + dirnameIn + fl;
/****************
Style
****************/
setTDRStyle();
gROOT->SetStyle("Plain");
gStyle->SetOptFit(1000);
gStyle->SetOptStat("emruo");
gStyle->SetOptStat(kFALSE);
gStyle->SetPadTickY(1);
gStyle->SetPadTickX(1);
int col_Nom = 1;
int col_Up = 2;
int col_Down = 4;
/****************
Channel
****************/
TString channel[3];
channel[0] = "mujets";
channel[1] = "ejets";
channel[2] = "lepjet";
/****************
ttbar Signal
****************/
std::vector<histos> ttbar_Nom;
std::vector<histos> ttbar_Up;
std::vector<histos> ttbar_Down;
ttbar_Nom = loadhistogramsSys(plots, files + "_ttbar_LepJetsPowhegPythia" + cat,"central");
ttbar_Up = loadhistogramsSys(plots, files + "_ttbar_LepJetsPowhegPythia" + cat, nSyst + "Up");
ttbar_Down = loadhistogramsSys(plots, files + "_ttbar_LepJetsPowhegPythia" + cat, nSyst + "Down");
setuphistograms(ttbar_Nom, col_Nom, 1);
setuphistograms(ttbar_Up, col_Up, 2);
setuphistograms(ttbar_Down, col_Down, 3);
/****************
Draw Histos
****************/
TCanvas *histocanvas;
histocanvas = new TCanvas("plots", "Plots");
for(unsigned int h=0; h<ttbar_Nom.size(); h++){
for(int ch=0; ch<3; ch++){
histocanvas->cd();
if(LogScale) histocanvas->cd()->SetLogy();
ttbar_Nom[h].hist[ch]->GetXaxis()->SetRange(ttbar_Nom[h].hist[ch]->GetXaxis()->GetFirst(),
ttbar_Nom[h].hist[ch]->GetXaxis()->GetLast());
//ttbar_Nom[h].hist[ch]->GetYaxis()->SetTitle("Events");
ttbar_Nom[h].hist[ch]->GetYaxis()->SetTitleOffset(1.2);
ttbar_Nom[h].hist[ch]->GetYaxis()->SetTitleSize(0.07);
ttbar_Nom[h].hist[ch]->GetYaxis()->SetLabelSize(0.055);
ttbar_Nom[h].hist[ch]->GetYaxis()->SetNdivisions(607);
//ttbar_Nom[h].hist[ch]->GetYaxis()->SetLabelSize(0.05);
TGaxis *hYaxis = (TGaxis*)ttbar_Nom[h].hist[ch]->GetYaxis();
//hYaxis->SetMaxDigits(3);
//ttbar_Nom[h].hist[ch]->GetXaxis()->SetLabelSize(0.0);
//ttbar_Nom[h].hist[ch]->GetXaxis()->SetTitle("");
ttbar_Nom[h].hist[ch]->GetXaxis()->SetNdivisions(509); //(402)
ttbar_Nom[h].hist[ch]->GetXaxis()->SetTitleOffset(1.1);
// Produce enough vertical space for the legend
float MaxHisto;
if(LogScale) MaxHisto = 8.0;
else MaxHisto = 1.4;
float maxh;
maxh = ttbar_Nom[h].hist[ch]->GetMaximum();
if (maxh < ttbar_Up[h].hist[ch]->GetMaximum()) maxh = ttbar_Up[h].hist[ch]->GetMaximum();
if (maxh < ttbar_Down[h].hist[ch]->GetMaximum()) maxh = ttbar_Down[h].hist[ch]->GetMaximum();
ttbar_Nom[h].hist[ch]->SetMaximum(MaxHisto*maxh);
//ttbar_Nom[h].hist[ch]->SetMinimum(0.1);
ttbar_Nom[h].hist[ch]->Draw("hist");
ttbar_Up[h].hist[ch]->Draw("histSAME");
ttbar_Down[h].hist[ch]->Draw("histSAME");
//ttbar_1[h].hist[ch]->Draw("histoSAME");
/***********************
Legends
***********************/
TLegend *leg;
float legx1=0.70;
float legy1=0.75;
float legx2=0.93;
float legy2=0.93;
leg = new TLegend(legx1,legy1,legx2,legy2);
leg->SetFillColor(0);
leg->SetLineColor(1);
leg->SetTextFont(62);
leg->SetTextSize(0.03);
leg->AddEntry((TObject*)0, cat + ": " + nSyst,"");
leg->AddEntry((TObject*)0,"","");
leg->AddEntry(ttbar_Nom[h].hist[ch],"Nom","l");
leg->AddEntry(ttbar_Up[h].hist[ch],"Up","l");
leg->AddEntry(ttbar_Down[h].hist[ch],"Down","l");
//leg->AddEntry((TObject*)0,"","");
leg->Draw("SAME");
//-------------------------------------------------------
// CMS Legend
//-------------------------------------------------------
TString htitleCMSChannel[3];
htitleCMSChannel[0] = "#mu^{#pm}+jets channel";
htitleCMSChannel[1] = "e^{#pm}+jets channel";
htitleCMSChannel[2] = "l^{#pm}+jets channel";
TLatex *titlePr;
titlePr = new TLatex(-20.,50.,"Preliminary");
titlePr->SetNDC();
titlePr->SetTextAlign(12);
titlePr->SetX(0.25);
titlePr->SetY(0.97);
titlePr->SetTextColor(2);
titlePr->SetTextFont(42);
titlePr->SetTextSize(0.05);
titlePr->SetTextSizePixels(24);
titlePr->Draw("SAME");
TLatex *title;
title = new TLatex(-20.,50.,"CMS #sqrt{s} = 13TeV, L = 36.5 fb^{-1}");
title->SetNDC();
title->SetTextAlign(12);
title->SetX(0.20);
title->SetY(0.80);
title->SetTextFont(42);
title->SetTextSize(0.03);
title->SetTextSizePixels(24);
title->Draw("SAME");
TLatex *chtitle;
chtitle = new TLatex(-20.,50.,htitleCMSChannel[ch]);
chtitle->SetNDC();
chtitle->SetTextAlign(12);
chtitle->SetX(0.20);
chtitle->SetY(0.86);
chtitle->SetTextFont(42);
chtitle->SetTextSize(0.03);
chtitle->SetTextSizePixels(24);
chtitle->Draw("SAME");
/***********************
Save Histos
***********************/
TString dirfigname_log;
if(LogScale) dirfigname_log = "_log";
else dirfigname_log = "";
TString dirfigname_pdf;
TString dirfigname_png;
dirfigname_pdf = "../" + dirnameIn + "figuresNomUpDown_" + fl + "/ttbbNomUpDown_" + cat + "/pdf" + dirfigname_log + "/";
//dirfigname_png = dirnameIn + "figuresSystComp_" + fl + "/ttbb/png" + dirfigname_log + "/";
// make a dir if it does not exist!!
gSystem->mkdir(dirfigname_pdf, kTRUE);
histocanvas->SaveAs(dirfigname_pdf + nSyst + "_" + ttbar_Nom[h].hist[ch]->GetName() + ".pdf");
//gSystem->mkdir(dirfigname_png, kTRUE);
//histocanvas->SaveAs(dirfigname_png + WJets[h].hist[ch]->GetName() + ".png");
// clear Canvas
histocanvas->Clear();
}
}
} //end Plots.C
