void drawBeamSpotPdfs (TDirectory* directory, const char* coord,
const char* fname)
{
TH3* obsHisto = directory->Get("PVobs");
TH3* estHisto = directory->Get("PVest");
if ( obsHisto==0 || estHisto==0 ) return;
std::string fullName("PV");
fullName += coord;
if ( fname ) fullName += fname;
else fullName += directory->GetName();
TCanvas* c = new TCanvas(fullName.c_str(),fullName.c_str());
c->SetLogy(true);
TH1* obsHisto1D = obsHisto->Project3D(coord);
TH1* estHisto1D = estHisto->Project3D(coord);
TAxis* xaxis = obsHisto1D->GetXaxis();
string atitle("Primary vertex ");
atitle += coord;
atitle += " [cm]";
xaxis->SetTitle(atitle.c_str());
if ( c->GetLogy() ) obsHisto1D->SetMinimum(0.5);
obsHisto1D->SetMarkerStyle(21);
// obsHisto1D->SetLineColor(2);
// obsHisto1D->SetMarkerColor(2);
obsHisto1D->Draw("E0");
estHisto1D->SetLineColor(2);
estHisto1D->SetLineWidth(2);
estHisto1D->Draw("hist same c");
if ( strcmp(coord,"x") || strcmp(coord,"y") ) xaxis->SetNdivisions(508);
string epsName = fullName + ".eps";
c->SaveAs(epsName.c_str());
string pngName = fullName + ".png";
c->SaveAs(pngName.c_str());
}
//~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
//
// drawFigure6
//
// logy = 0
// logy = 1
//
//~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
void drawFigure6(Int_t logy = 0)
{
gInterpreter->ExecuteMacro("WZPaperStyle.C");
gSystem->mkdir("pdf", kTRUE);
gSystem->mkdir("png", kTRUE);
// Read the input file
//----------------------------------------------------------------------------
TString name = "WZ_PtZ_plot_allCh_largeATGC";
TFile* file = new TFile("rootfiles/" + name + ".root", "read");
TCanvas* c1 = (TCanvas*)file->Get("c1_allCh");
TH1F* data = (TH1F*)c1->FindObject("histo_data_3e");
TH1F* aTGC_dk = (TH1F*)c1->FindObject("histo_aTGC_dk_3e");
TH1F* aTGC_lam = (TH1F*)c1->FindObject("histo_aTGC_lam_3e");
TH1F* aTGC_dg = (TH1F*)c1->FindObject("histo_aTGC_dg_3e");
TH1F* WZ = (TH1F*)c1->FindObject("histo_SM_3e");
TH1F* fakes = (TH1F*)c1->FindObject("fake_0");
TH1F* ZZ = (TH1F*)c1->FindObject("total_bkg_rebined_ZZ_0");
TH1F* Zgamma = (TH1F*)c1->FindObject("total_bkg_rebined_Zgamma_0");
TH1F* WV = (TH1F*)c1->FindObject("total_bkg_rebined_WV_0");
TH1F* VVV = (TH1F*)c1->FindObject("total_bkg_rebined_VVV_0");
WZ->SetFillColor(kOrange-2);
WZ->SetLineColor(kOrange-2);
Zgamma->SetFillColor(kRed+1);
Zgamma->SetLineColor(kRed+1);
ZZ->SetFillColor(kRed+1);
ZZ->SetLineColor(kRed+1);
fakes->SetFillColor(kGray+1);
fakes->SetLineColor(kGray+1);
WV->SetFillColor(kRed+1);
WV->SetLineColor(kRed+1);
VVV->SetFillColor(kRed+1);
VVV->SetLineColor(kRed+1);
// Draw
//----------------------------------------------------------------------------
TCanvas* canvas = new TCanvas("canvas", "canvas");
canvas->SetLogy(logy);
data->Draw("ep");
// Axis labels
//----------------------------------------------------------------------------
TAxis* xaxis = data->GetXaxis();
TAxis* yaxis = data->GetYaxis();
xaxis->SetLabelFont ( 42);
xaxis->SetLabelOffset(0.01);
xaxis->SetLabelSize (0.05);
xaxis->SetNdivisions ( 505);
xaxis->SetTitleFont ( 42);
xaxis->SetTitleOffset( 1.2);
xaxis->SetTitleSize (0.05);
yaxis->SetLabelFont ( 42);
yaxis->SetLabelOffset(0.01);
yaxis->SetLabelSize (0.05);
yaxis->SetNdivisions ( 505);
yaxis->SetTitleFont ( 42);
yaxis->SetTitleOffset( 1.6);
yaxis->SetTitleSize (0.05);
xaxis->SetTitle("p_{T}^{Z} (GeV)");
yaxis->SetTitle(Form("Events / %.0f GeV", data->GetBinWidth(0)));
// Adjust scale
//----------------------------------------------------------------------------
Float_t theMax = GetMaximumIncludingErrors(data);
Float_t theMaxMC = GetMaximumIncludingErrors(aTGC_dk);
if (theMaxMC > theMax) theMax = theMaxMC;
if (canvas->GetLogy()) {
data->SetMaximum(15 * theMax);
data->SetMinimum(1);
} else {
data->SetMaximum(1.2 * theMax);
}
// Legend
//----------------------------------------------------------------------------
Double_t x0;
Double_t y0;
if (logy)
{
x0 = 0.630;
y0 = 0.765;
DrawTLegend(x0 - 0.37, y0 + 2.*(_yoffset+0.001), data, " Data", "ep");
DrawTLegend(x0 - 0.37, y0 + 1.*(_yoffset+0.001), aTGC_dk, " WZ aTGC (#Delta#kappa^{Z} = 0.6)", "l");
DrawTLegend(x0 - 0.37, y0, aTGC_dg, " WZ aTGC (#Deltag^{Z}_{1} = -0.06)", "l");
DrawTLegend(x0 - 0.37, y0 - 1.*(_yoffset+0.001), aTGC_lam, " WZ aTGC (#lambda = 0.04)", "l");
DrawTLegend(x0, y0 + 2.*(_yoffset+0.001), WZ, " WZ", "f");
DrawTLegend(x0, y0 + 1.*(_yoffset+0.001), fakes, " Non-prompt leptons", "f");
DrawTLegend(x0, y0, ZZ, " MC background", "f");
}
else
{
x0 = 0.570;
y0 = 0.755;
DrawTLegend(x0, y0 + 2.*(_yoffset+0.001), data, " Data", "ep");
DrawTLegend(x0, y0 + 1.*(_yoffset+0.001), aTGC_dk, " WZ aTGC (#Delta#kappa^{Z} = 0.6)", "l");
DrawTLegend(x0, y0, aTGC_dg, " WZ aTGC (#Deltag^{Z}_{1} = -0.06)", "l");
DrawTLegend(x0, y0 - 1.*(_yoffset+0.001), aTGC_lam, " WZ aTGC (#lambda = 0.04)", "l");
DrawTLegend(x0, y0 - 2.*(_yoffset+0.001), WZ, " WZ", "f");
DrawTLegend(x0, y0 - 3.*(_yoffset+0.001), fakes, " Non-prompt leptons", "f");
DrawTLegend(x0, y0 - 4.*(_yoffset+0.001), ZZ, " MC background", "f");
}
// Finish it
//----------------------------------------------------------------------------
data->SetTitle("");
if (logy)
{
DrawTLatex(_cmsTextFont, 0.190, 0.94, 0.055, 11, "CMS");
// DrawTLatex(_extraTextFont, 0.315, 0.94, 0.030, 11, "Preliminary");
}
else
{
DrawTLatex(_cmsTextFont, 0.215, 0.891, 0.055, 13, "CMS");
// DrawTLatex(_extraTextFont, 0.215, 0.837, 0.030, 13, "Preliminary");
}
DrawTLatex(_lumiTextFont, 0.940, 0.94, 0.040, 31, "19.6 fb^{-1} (8 TeV)");
WZ ->Draw("hist,same");
aTGC_dk ->Draw("hist,same");
aTGC_lam->Draw("hist,same");
aTGC_dg ->Draw("hist,same");
fakes ->Draw("hist,same");
ZZ ->Draw("hist,same");
Zgamma ->Draw("hist,same");
WV ->Draw("hist,same");
VVV ->Draw("hist,same");
data ->Draw("ep,same");
canvas->GetFrame()->DrawClone();
canvas->RedrawAxis();
canvas->Update();
TString cname = name;
if (logy) cname += "_log_range";
canvas->SaveAs("pdf/" + cname + ".pdf");
canvas->SaveAs("png/" + cname + ".png");
}