//------------------------------------------------------------------------------
// SetAxis
//------------------------------------------------------------------------------
void SetAxis(TH1* hist,
TString xtitle,
TString ytitle,
Float_t size,
Float_t offset)
{
TAxis* xaxis = (TAxis*)hist->GetXaxis();
TAxis* yaxis = (TAxis*)hist->GetYaxis();
xaxis->SetLabelFont(42);
yaxis->SetLabelFont(42);
xaxis->SetTitleFont(42);
yaxis->SetTitleFont(42);
xaxis->SetLabelOffset(0.025);
yaxis->SetLabelOffset(0.025);
xaxis->SetTitleOffset(1.4);
yaxis->SetTitleOffset(offset);
xaxis->SetLabelSize(size);
yaxis->SetLabelSize(size);
xaxis->SetTitleSize(size);
yaxis->SetTitleSize(size);
xaxis->SetTitle(xtitle);
yaxis->SetTitle(ytitle);
xaxis->SetNdivisions(505);
yaxis->SetNdivisions(505);
yaxis->CenterTitle();
gPad->GetFrame()->DrawClone();
gPad->RedrawAxis();
}
void showGraph(double canvasSizeX, double canvasSizeY,
TGraph* graph,
bool useLogScaleX, double xMin, double xMax, const std::string& xAxisTitle, double xAxisOffset,
bool useLogScaleY, double yMin, double yMax, const std::string& yAxisTitle, double yAxisOffset,
const std::string& outputFileName)
{
TCanvas* canvas = new TCanvas("canvas", "canvas", canvasSizeX, canvasSizeY);
canvas->SetFillColor(10);
canvas->SetBorderSize(2);
canvas->SetTopMargin(0.05);
canvas->SetLeftMargin(0.19);
canvas->SetBottomMargin(0.19);
canvas->SetRightMargin(0.05);
canvas->SetLogx(useLogScaleX);
canvas->SetLogy(useLogScaleY);
TH1* dummyHistogram = new TH1D("dummyHistogram", "dummyHistogram", 10, xMin, xMax);
dummyHistogram->SetTitle("");
dummyHistogram->SetStats(false);
dummyHistogram->SetMinimum(yMin);
dummyHistogram->SetMaximum(yMax);
dummyHistogram->Draw("axis");
TAxis* xAxis = dummyHistogram->GetXaxis();
xAxis->SetTitle(xAxisTitle.data());
xAxis->SetTitleOffset(xAxisOffset);
xAxis->SetTitleSize(0.065);
xAxis->SetLabelSize(0.055);
xAxis->SetLabelOffset(0.01);
xAxis->SetTickLength(0.055);
xAxis->SetNdivisions(505);
TAxis* yAxis = dummyHistogram->GetYaxis();
yAxis->SetTitle(yAxisTitle.data());
yAxis->SetTitleOffset(yAxisOffset);
yAxis->SetTitleSize(0.070);
yAxis->SetLabelSize(0.055);
yAxis->SetLabelOffset(0.01);
yAxis->SetTickLength(0.055);
yAxis->SetNdivisions(505);
graph->SetMarkerColor(1);
graph->SetLineColor(1);
graph->Draw("p");
canvas->Update();
size_t idx = outputFileName.find_last_of('.');
std::string outputFileName_plot = std::string(outputFileName, 0, idx);
if ( useLogScaleY ) outputFileName_plot.append("_log");
else outputFileName_plot.append("_linear");
if ( idx != std::string::npos ) canvas->Print(std::string(outputFileName_plot).append(std::string(outputFileName, idx)).data());
canvas->Print(std::string(outputFileName_plot).append(".png").data());
//canvas->Print(std::string(outputFileName_plot).append(".pdf").data());
//canvas->Print(std::string(outputFileName_plot).append(".root").data());
delete dummyHistogram;
delete canvas;
}
void showHistograms(double canvasSizeX, double canvasSizeY,
TH1* histogram,
double xMin, double xMax, const std::string& xAxisTitle, double xAxisOffset,
bool useLogScale, double yMin, double yMax, const std::string& yAxisTitle, double yAxisOffset,
const std::string& outputFileName)
{
TCanvas* canvas = new TCanvas("canvas", "canvas", canvasSizeX, canvasSizeY);
canvas->SetFillColor(10);
canvas->SetBorderSize(2);
canvas->SetLeftMargin(0.15);
canvas->SetBottomMargin(0.15);
canvas->SetLogy(useLogScale);
histogram->SetTitle("");
histogram->SetStats(true);
histogram->SetMinimum(yMin);
histogram->SetMaximum(yMax);
histogram->SetLineColor(1);
histogram->SetLineWidth(2);
histogram->SetMarkerColor(1);
histogram->SetMarkerStyle(20);
histogram->SetMarkerSize(1.5);
histogram->Draw("hist");
TAxis* xAxis = histogram->GetXaxis();
xAxis->SetRangeUser(xMin, xMax);
xAxis->SetTitle(xAxisTitle.data());
xAxis->SetTitleSize(0.060);
xAxis->SetTitleOffset(xAxisOffset);
xAxis->SetLabelSize(0.050);
xAxis->SetNdivisions(505);
TAxis* yAxis = histogram->GetYaxis();
yAxis->SetTitle(yAxisTitle.data());
yAxis->SetTitleSize(0.060);
yAxis->SetTitleOffset(yAxisOffset);
yAxis->SetLabelSize(0.050);
yAxis->SetNdivisions(505);
canvas->Update();
size_t idx = outputFileName.find_last_of('.');
std::string outputFileName_plot = std::string(outputFileName, 0, idx);
if ( useLogScale ) outputFileName_plot.append("_log");
else outputFileName_plot.append("_linear");
if ( idx != std::string::npos ) canvas->Print(std::string(outputFileName_plot).append(std::string(outputFileName, idx)).data());
//canvas->Print(std::string(outputFileName_plot).append(".png").data());
canvas->Print(std::string(outputFileName_plot).append(".pdf").data());
//canvas->Print(std::string(outputFileName_plot).append(".root").data());
delete canvas;
}
void LEPStyle::setHist()
{
_hist = new TH1F( "bidon", "bidon", 100, _xmin, _xmax );
_ax = _hist->GetXaxis();
_ay = _hist->GetYaxis();
_ax->SetTitle(_title);
_ax->SetNdivisions(_ndivx);
if( _logX )
{
_ax->SetMoreLogLabels();
}
_ax->SetLabelSize(_scale*0.060);
// _ax->SetLabelOffset(_scale*0.007);
_ax->SetTitleSize(_scale*0.070);
_ax->SetTitleOffset(0.9);
_ax->CenterTitle( _centerTitle );
_ay->SetNdivisions(0);
}
void SetAxisLabels(TH1& hist, char* xtitle, char* ytitle="",
double xoffset=1.1, double yoffset=1.4) {
TAxis* x = hist.GetXaxis();
TAxis* y = hist.GetYaxis();
x->SetTitle(xtitle);
x->SetTitleSize(0.06);
x->SetLabelSize(0.05);
x->SetTitleOffset(xoffset);
x->SetNdivisions(505);
y->SetTitle(ytitle);
y->SetTitleSize(0.06);
y->SetLabelSize(0.05);
y->SetTitleOffset(yoffset);
y->SetNoExponent();
hist.SetLineWidth(2);
hist.SetMarkerStyle(20);
std::stringstream str;
str << "Events / " << (int) lumi << " pb^{-1} ";
std::string defYtitle = str.str();
if(ytitle=="") y->SetTitle( defYtitle.c_str() );
}
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());
}
void makePlot(const std::string& inputFilePath, const std::string& canvasName, const std::string& sample, int massPoint, const std::string& channel, double k,
const std::string& inputFileName, const std::string& outputFilePath, const std::string& outputFileName)
{
std::string inputFileName_full = Form("%s%s", inputFilePath.data(), inputFileName.data());
TFile* inputFile = new TFile(inputFileName_full.data());
if ( !inputFile ) {
std::cerr << "Failed to open input file = " << inputFileName_full << " !!" << std::endl;
assert(0);
}
inputFile->ls();
TCanvas* canvas = dynamic_cast<TCanvas*>(inputFile->Get(canvasName.data()));
if ( !canvas ) {
std::cerr << "Failed to load canvas = " << canvasName << " !!" << std::endl;
assert(0);
}
int idxPad = -1;
if ( massPoint == 90 ) idxPad = 1;
if ( massPoint == 125 ) idxPad = 2;
if ( massPoint == 200 ) idxPad = 3;
if ( massPoint == 300 ) idxPad = 4;
if ( massPoint == 500 ) idxPad = 5;
if ( massPoint == 800 ) idxPad = 6;
if ( !(idxPad >= 1 && idxPad <= 6) ) {
std::cerr << "Invalid sample = " << sample << " !!" << std::endl;
assert(0);
}
TVirtualPad* pad = canvas->GetPad(idxPad);
std::cout << "pad = " << pad << ": ClassName = " << pad->ClassName() << std::endl;
TCanvas* canvas_new = new TCanvas("canvas_new", "canvas_new", 900, 800);
canvas_new->SetFillColor(10);
canvas_new->SetBorderSize(2);
canvas_new->SetTopMargin(0.065);
canvas_new->SetLeftMargin(0.17);
canvas_new->SetBottomMargin(0.165);
canvas_new->SetRightMargin(0.015);
canvas_new->SetLogx(true);
canvas_new->SetLogy(true);
canvas_new->Draw();
canvas_new->cd();
//TList* pad_primitives = canvas->GetListOfPrimitives();
TList* pad_primitives = pad->GetListOfPrimitives();
TH1* histogramCA = 0;
TH1* histogramSVfit = 0;
TH1* histogramSVfitMEMkEq0 = 0;
TH1* histogramSVfitMEMkNeq0 = 0;
TIter pad_nextObj(pad_primitives);
while ( TObject* obj = pad_nextObj() ) {
std::string objName = "";
if ( dynamic_cast<TNamed*>(obj) ) objName = (dynamic_cast<TNamed*>(obj))->GetName();
std::cout << "obj = " << obj << ": name = " << objName << ", type = " << obj->ClassName() << std::endl;
TH1* tmpHistogram = dynamic_cast<TH1*>(obj);
if ( tmpHistogram ) {
std::cout << "tmpHistogram:"
<< " fillColor = " << tmpHistogram->GetFillColor() << ", fillStyle = " << tmpHistogram->GetFillStyle() << ","
<< " lineColor = " << tmpHistogram->GetLineColor() << ", lineStyle = " << tmpHistogram->GetLineStyle() << ", lineWidth = " << tmpHistogram->GetLineWidth() << ","
<< " markerColor = " << tmpHistogram->GetMarkerColor() << ", markerStyle = " << tmpHistogram->GetMarkerStyle() << ", markerSize = " << tmpHistogram->GetMarkerSize() << ","
<< " integral = " << tmpHistogram->Integral() << std::endl;
std::cout << "(mean = " << tmpHistogram->GetMean() << ", rms = " << tmpHistogram->GetRMS() << ": rms/mean = " << (tmpHistogram->GetRMS()/tmpHistogram->GetMean()) << ")" << std::endl;
if ( tmpHistogram->GetLineColor() == 416 ) histogramCA = tmpHistogram;
if ( tmpHistogram->GetLineColor() == 600 ) histogramSVfit = tmpHistogram;
if ( tmpHistogram->GetLineColor() == 616 ) histogramSVfitMEMkEq0 = tmpHistogram;
if ( tmpHistogram->GetLineColor() == 632 ) histogramSVfitMEMkNeq0 = tmpHistogram;
}
}
if ( !(histogramCA && histogramSVfit && histogramSVfitMEMkEq0 && histogramSVfitMEMkNeq0) ) {
std::cerr << "Failed to load histograms !!" << std::endl;
assert(0);
}
//gStyle->SetLineStyleString(2,"40 10 10 10 10 10 10 10");
//gStyle->SetLineStyleString(3,"25 15");
//gStyle->SetLineStyleString(4,"60 25");
//int colors[4] = { kBlack, kGreen - 6, kBlue - 7, kMagenta - 7 };
int colors[4] = { 28, kGreen - 6, kBlue - 7, kBlack };
//int lineStyles[4] = { 2, 3, 4, 1 };
int lineStyles[4] = { 7, 1, 1, 1 };
//int lineWidths[4] = { 3, 3, 4, 3 };
int lineWidths[4] = { 3, 3, 1, 1 };
int markerStyles[4] = { 20, 25, 21, 24 };
int markerSizes[4] = { 2, 2, 2, 2 };
histogramCA->SetFillColor(0);
histogramCA->SetFillStyle(0);
histogramCA->SetLineColor(colors[0]);
histogramCA->SetLineStyle(lineStyles[0]);
histogramCA->SetLineWidth(lineWidths[0]);
histogramCA->SetMarkerColor(colors[0]);
histogramCA->SetMarkerStyle(markerStyles[0]);
histogramCA->SetMarkerSize(markerSizes[0]);
histogramSVfit->SetFillColor(0);
histogramSVfit->SetFillStyle(0);
histogramSVfit->SetLineColor(colors[1]);
histogramSVfit->SetLineStyle(lineStyles[1]);
histogramSVfit->SetLineWidth(lineWidths[1]);
histogramSVfit->SetMarkerColor(colors[1]);
histogramSVfit->SetMarkerStyle(markerStyles[1]);
histogramSVfit->SetMarkerSize(markerSizes[1]);
histogramSVfitMEMkEq0->SetFillColor(0);
histogramSVfitMEMkEq0->SetFillStyle(0);
histogramSVfitMEMkEq0->SetLineColor(colors[2]);
histogramSVfitMEMkEq0->SetLineStyle(lineStyles[2]);
histogramSVfitMEMkEq0->SetLineWidth(lineWidths[2]);
histogramSVfitMEMkEq0->SetMarkerColor(colors[2]);
histogramSVfitMEMkEq0->SetMarkerStyle(markerStyles[2]);
histogramSVfitMEMkEq0->SetMarkerSize(markerSizes[2]);
// CV: fix pathological bins at high mass for which dN/dm increases
int numBins = histogramSVfitMEMkEq0->GetNbinsX();
for ( int idxBin = 1; idxBin <= numBins; ++idxBin ) {
double binCenter = histogramSVfitMEMkEq0->GetBinCenter(idxBin);
if ( (channel == "#tau_{h}#tau_{h}" && massPoint == 500 && binCenter > 1500.) ||
(channel == "#tau_{h}#tau_{h}" && massPoint == 800 && binCenter > 2000.) ||
(channel == "#mu#tau_{h}" && massPoint == 500 && binCenter > 1500.) ||
(channel == "#mu#tau_{h}" && massPoint == 800 && binCenter > 2500.) ) {
histogramSVfitMEMkEq0->SetBinContent(idxBin, 0.);
}
}
histogramSVfitMEMkNeq0->SetFillColor(0);
histogramSVfitMEMkNeq0->SetFillStyle(0);
histogramSVfitMEMkNeq0->SetLineColor(colors[3]);
histogramSVfitMEMkNeq0->SetLineStyle(lineStyles[3]);
histogramSVfitMEMkNeq0->SetLineWidth(lineWidths[3]);
histogramSVfitMEMkNeq0->SetMarkerColor(colors[3]);
histogramSVfitMEMkNeq0->SetMarkerStyle(markerStyles[3]);
histogramSVfitMEMkNeq0->SetMarkerSize(markerSizes[3]);
TAxis* xAxis = histogramCA->GetXaxis();
xAxis->SetTitle("m_{#tau#tau} [GeV]");
xAxis->SetTitleOffset(1.15);
xAxis->SetTitleSize(0.070);
xAxis->SetTitleFont(42);
xAxis->SetLabelOffset(0.010);
xAxis->SetLabelSize(0.055);
xAxis->SetLabelFont(42);
xAxis->SetTickLength(0.040);
xAxis->SetNdivisions(510);
//double xMin = 20.;
//double xMax = xAxis->GetXmax();
//xAxis->SetRangeUser(xMin, xMax);
TAxis* yAxis = histogramCA->GetYaxis();
yAxis->SetTitle("dN/dm_{#tau#tau} [1/GeV]");
yAxis->SetTitleOffset(1.20);
yAxis->SetTitleSize(0.070);
yAxis->SetTitleFont(42);
yAxis->SetLabelOffset(0.010);
yAxis->SetLabelSize(0.055);
yAxis->SetLabelFont(42);
yAxis->SetTickLength(0.040);
yAxis->SetNdivisions(505);
double massPoint_double = 0.;
if ( massPoint == 90 ) massPoint_double = 91.2;
else massPoint_double = massPoint;
double dLog = (TMath::Log(5.*massPoint_double) - TMath::Log(50.))/25.; // xMin = 50, xMax = 5*massPoint, numBins = 25
double binWidth = TMath::Exp(TMath::Log(massPoint_double) + 0.5*dLog) - TMath::Exp(TMath::Log(massPoint_double) - 0.5*dLog);
double sf_binWidth = 1./binWidth;
std::cout << "massPoint = " << massPoint << ": sf_binWidth = " << sf_binWidth << std::endl;
histogramCA->SetTitle("");
histogramCA->SetStats(false);
histogramCA->SetMaximum(sf_binWidth*0.79);
histogramCA->SetMinimum(sf_binWidth*1.1e-4);
histogramCA->Draw("hist");
histogramSVfit->Draw("histsame");
//histogramSVfitMEMkEq0->Draw("histsame");
histogramSVfitMEMkEq0->Draw("epsame");
//histogramSVfitMEMkNeq0->Draw("histsame");
histogramSVfitMEMkNeq0->Draw("epsame");
histogramCA->Draw("axissame");
//TPaveText* label_sample = new TPaveText(0.21, 0.86, 0.46, 0.94, "NDC");
TPaveText* label_sample = new TPaveText(0.1700, 0.9475, 0.4600, 1.0375, "NDC");
label_sample->SetFillStyle(0);
label_sample->SetBorderSize(0);
label_sample->AddText(sample.data());
label_sample->SetTextFont(42);
label_sample->SetTextSize(0.055);
label_sample->SetTextColor(1);
label_sample->SetTextAlign(13);
label_sample->Draw();
//TLegend* legend_new = new TLegend(0.225, 0.52, 0.41, 0.82, NULL, "brNDC");
TLegend* legend_new = new TLegend(0.30, 0.30, 0.80, 0.80, NULL, "brNDC");
legend_new->SetFillColor(10);
legend_new->SetFillStyle(0);
legend_new->SetBorderSize(0);
legend_new->SetTextFont(42);
legend_new->SetTextSize(0.055);
legend_new->SetTextColor(1);
legend_new->SetMargin(0.20);
legend_new->AddEntry(histogramCA, "CA", "l");
legend_new->AddEntry(histogramSVfit, "SVfit", "l");
//legend_new->AddEntry(histogramSVfitMEMkEq0, "SVfitMEM (k=0)", "l");
legend_new->AddEntry(histogramSVfitMEMkEq0, "SVfitMEM (k=0)", "p");
//legend_new->AddEntry(histogramSVfitMEMkNeq0, Form("SVfitMEM(k=%1.0f)", k), "l");
legend_new->AddEntry(histogramSVfitMEMkNeq0, Form("SVfitMEM (k=%1.0f)", k), "p");
//legend_new->Draw();
double label_channel_y0;
if ( channel == "e#mu" ) label_channel_y0 = 0.9275;
else if ( channel == "#mu#tau_{h}" ) label_channel_y0 = 0.9400;
else if ( channel == "#tau_{h}#tau_{h}" ) label_channel_y0 = 0.9350;
else {
std::cerr << "Invalid channel = " << channel << " !!" << std::endl;
assert(0);
}
TPaveText* label_channel = new TPaveText(0.895, label_channel_y0, 0.975, label_channel_y0 + 0.055, "NDC");
label_channel->SetFillStyle(0);
label_channel->SetBorderSize(0);
label_channel->AddText(channel.data());
label_channel->SetTextFont(62);
label_channel->SetTextSize(0.055);
label_channel->SetTextColor(1);
label_channel->SetTextAlign(31);
label_channel->Draw();
canvas_new->Update();
std::string outputFileName_full = Form("%s%s", outputFilePath.data(), outputFileName.data());
size_t idx = outputFileName_full.find_last_of('.');
std::string outputFileName_plot = std::string(outputFileName_full, 0, idx);
canvas_new->Print(std::string(outputFileName_plot).append(".pdf").data());
canvas_new->Print(std::string(outputFileName_plot).append(".root").data());
std::string channel_string;
if ( channel == "e#mu" ) channel_string = "emu";
else if ( channel == "#mu#tau_{h}" ) channel_string = "muhad";
else if ( channel == "#tau_{h}#tau_{h}" ) channel_string = "hadhad";
else {
std::cerr << "Invalid channel = " << channel << " !!" << std::endl;
assert(0);
}
std::string outputFileName_legend = Form("makeSVfitMEM_PerformancePlots_legend_%s.pdf", channel_string.data());
makePlot_legend(legend_new, outputFilePath, outputFileName_legend);
delete label_sample;
delete legend_new;
delete label_channel;
delete canvas_new;
delete inputFile;
}
//~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
//
// drawFigure7
//
// parameter = "dg0"
// parameter = "lam0"
// parameter = "dk0"
//
//~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
void drawFigure7(TString parameter = "dk0")
{
gInterpreter->ExecuteMacro("WZPaperStyle.C");
gSystem->mkdir("pdf", kTRUE);
gSystem->mkdir("png", kTRUE);
// Individual settings
//----------------------------------------------------------------------------
if (parameter.Contains("dg0"))
{
xtitle = "#Delta#kappa^{Z}";
ytitle = "#lambda";
xmin = -0.5;
xmax = 0.5;
ymin = -0.039;
ymax = 0.053;
}
else if (parameter.Contains("lam0"))
{
xtitle = "#Delta#kappa^{Z}";
ytitle = "#Deltag^{Z}_{1}";
xmin = -0.53;
xmax = 0.53;
ymin = -0.05;
ymax = 0.09;
}
else if (parameter.Contains("dk0"))
{
xtitle = "#Deltag^{Z}_{1}";
ytitle = "#lambda";
xmin = -0.059;
xmax = 0.059;
ymin = -0.037;
ymax = 0.052;
}
// Read the input file
//----------------------------------------------------------------------------
TFile* file = new TFile("rootfiles/contours_" + parameter + "_2Dpol2.root", "read");
TGraph* cont_exp_68 = (TGraph*)file->Get("cont_exp_68");
TGraph* cont_exp_95 = (TGraph*)file->Get("cont_exp_95");
TGraph* cont_exp_99 = (TGraph*)file->Get("cont_exp_99");
TGraph* cont_obs_95 = (TGraph*)file->Get("cont_obs_95");
TGraph* bestFit = (TGraph*)file->Get("bestFit");
bestFit->SetMarkerSize(1.5);
bestFit->SetMarkerStyle(34);
TGraph* g_large = new TGraph(2);
g_large->SetPoint(0, -1.0, 1.0);
g_large->SetPoint(1, 1.0, 1.0);
// Draw
//----------------------------------------------------------------------------
TCanvas* canvas = new TCanvas("canvas", "canvas");
TMultiGraph* mg = new TMultiGraph();
mg->Add(cont_exp_68);
mg->Add(cont_exp_95);
mg->Add(cont_exp_99);
mg->Add(cont_obs_95);
mg->Add(bestFit, "p");
mg->Add(g_large);
mg->Draw("ac");
mg->SetMinimum(ymin);
mg->SetMaximum(ymax);
// Axis labels
//----------------------------------------------------------------------------
TAxis* xaxis = mg->GetXaxis();
TAxis* yaxis = mg->GetYaxis();
xaxis->SetLabelFont ( 42);
xaxis->SetLabelOffset( 0.01);
xaxis->SetLabelSize ( 0.05);
xaxis->SetNdivisions ( 505);
xaxis->SetTitle (xtitle);
xaxis->SetTitleFont ( 42);
xaxis->SetTitleOffset( 1.2);
xaxis->SetTitleSize ( 0.05);
xaxis->SetLimits(xmin, xmax);
yaxis->SetLabelFont ( 42);
yaxis->SetLabelOffset( 0.01);
yaxis->SetLabelSize ( 0.05);
yaxis->SetNdivisions ( 505);
yaxis->SetTitle (ytitle);
yaxis->SetTitleFont ( 42);
yaxis->SetTitleOffset( 1.6);
yaxis->SetTitleSize ( 0.05);
canvas->Modified();
// Legend
//----------------------------------------------------------------------------
Double_t x0 = 0.580;
Double_t y0 = 0.755;
DrawTLegend(x0 - 0.36, y0 + 2.*(_yoffset+0.001), (TH1F*)cont_exp_68, " Expected 68% CL", "l");
DrawTLegend(x0 - 0.36, y0 + 1.*(_yoffset+0.001), (TH1F*)cont_exp_95, " Expected 95% CL", "l");
DrawTLegend(x0, y0 + 2.*(_yoffset+0.001), (TH1F*)cont_exp_99, " Expected 99% CL", "l");
DrawTLegend(x0, y0 + 1.*(_yoffset+0.001), (TH1F*)cont_obs_95, " Observed 95% CL", "l");
DrawTLegend(x0, y0, (TH1F*)bestFit, " Best fit", "p");
// Finish it
//----------------------------------------------------------------------------
DrawTLatex(_cmsTextFont, 0.190, 0.94, 0.055, 11, "CMS");
// DrawTLatex(_extraTextFont, 0.315, 0.94, 0.030, 11, "Preliminary");
DrawTLatex(_lumiTextFont, 0.940, 0.94, 0.040, 31, "19.6 fb^{-1} (8 TeV)");
canvas->SaveAs("pdf/lZ_dkg_2dlimit_" + parameter + "_2Dpol2_deltaNLL.pdf");
canvas->SaveAs("png/lZ_dkg_2dlimit_" + parameter + "_2Dpol2_deltaNLL.png");
}
void makeTable( int num )
{
// TString rootfilename;
// if(num<9) rootfilename = Form("root_electroneff_%d_eta_et.root", num);
// else if(num<17) rootfilename = Form("root_electroneff_10%d_eta_et.root", num-8);
// else if(num<25) rootfilename = Form("root_electroneff_20%d_eta_et.root", num-16);
// else if(num<33) rootfilename = Form("root_electroneff_30%d_eta_et.root", num-24);
// else rootfilename = Form("root_electroneff_40%d_eta_et.root", num-32);
TString rootfilename;
if(num==0) rootfilename = "root_electroneff_WenuOff_eta_et.root";
if(num==1) rootfilename = "root_electroneff_WenuOn_eta_et.root";
if(num==2) rootfilename = "root_electroneff_ZeeOff_eta_et.root";
if(num==3) rootfilename = "root_electroneff_ZeeOn_eta_et.root";
TString pname = "";
if(num==0) pname = "WenuOff-";
if(num==1) pname = "WenuOn-";
if(num==2) pname = "ZeeOff-";
if(num==3) pname = "ZeeOn-";
// if(num==1) pname = "WenuSelection-A1";
// if(num==2) pname = "WenuSelection-A2";
// if(num==3) pname = "WenuSelection-A3";
// if(num==4) pname = "WenuSelection-B1";
// if(num==5) pname = "WenuSelection-B2";
// if(num==6) pname = "WenuSelection-B3";
// if(num==7) pname = "ZeeSelection-A";
// if(num==8) pname = "ZeeSelection-B";
// if(num==9) pname = "WenuSelection-A1-noId-";
// if(num==10) pname = "WenuSelection-A2-noId-";
// if(num==11) pname = "WenuSelection-A3-noId-";
// if(num==12) pname = "WenuSelection-B1-noId-";
// if(num==13) pname = "WenuSelection-B2-noId-";
// if(num==14) pname = "WenuSelection-B3-noId-";
// if(num==15) pname = "ZeeSelection-A-noId-";
// if(num==16) pname = "ZeeSelection-B-noId-";
// if(num==17) pname = "WenuSelection-A1-noTrackIso-";
// if(num==18) pname = "WenuSelection-A2-noTrackIso-";
// if(num==19) pname = "WenuSelection-A3-noTrackIso-";
// if(num==20) pname = "WenuSelection-B1-noTrackIso-";
// if(num==21) pname = "WenuSelection-B2-noTrackIso-";
// if(num==22) pname = "WenuSelection-B3-noTrackIso-";
// if(num==23) pname = "ZeeSelection-A-noTrackIso-";
// if(num==24) pname = "ZeeSelection-B-noTrackIso-";
// if(num==25) pname = "WenuSelection-A1-noIso-";
// if(num==26) pname = "WenuSelection-A2-noIso-";
// if(num==27) pname = "WenuSelection-A3-noIso-";
// if(num==28) pname = "WenuSelection-B1-noIso-";
// if(num==29) pname = "WenuSelection-B2-noIso-";
// if(num==30) pname = "WenuSelection-B3-noIso-";
// if(num==31) pname = "ZeeSelection-A-noIso-";
// if(num==32) pname = "ZeeSelection-B-noIso-";
// if(num==33) pname = "WenuSelection-A1-noSigmaEtaEta-";
// if(num==34) pname = "WenuSelection-A2-noSigmaEtaEta-";
// if(num==35) pname = "WenuSelection-A3-noSigmaEtaEta-";
// if(num==36) pname = "WenuSelection-B1-noSigmaEtaEta-";
// if(num==37) pname = "WenuSelection-B2-noSigmaEtaEta-";
// if(num==38) pname = "WenuSelection-B3-noSigmaEtaEta-";
// if(num==39) pname = "ZeeSelection-A-noSigmaEtaEta-";
// if(num==40) pname = "ZeeSelection-B-noSigmaEtaEta-";
TString label = TString("electroneff-") + pname;
gROOT->ProcessLine(".L tdrstyleNew.C");
setTDRStyle();
tdrStyle->SetPadLeftMargin(0.16);
tdrStyle->SetPadRightMargin(0.1);
tdrStyle->SetPadTopMargin(0.08);
tdrStyle->SetLegendBorderSize(0);
TFile* myhistos = new TFile(rootfilename);
TH1F* h0 = (TH1F*) myhistos->Get("sbs_eff_EtDet");
TH1F* h1 = (TH1F*) myhistos->Get("sbs_eff_EtaDet");
TH2F* h2 = (TH2F*) myhistos->Get("sbs_eff_EtDet_EtaDet");
// TH1F *base = (TH1F*)f.Get("base");
// TH1F *destination = (TH1F*)f.Get("destination");
// TGraphAsymmErrors *g1 = new TGraphAsymmErrors();
// g1->BayesDivide(destination, base, "");
// g1->GetYaxis()->SetRangeUser(0.5, 1.05);
h0->SetMarkerStyle(20);
h1->SetMarkerStyle(20);
ScaleTo10pb( *h0 );
ScaleTo10pb( *h1 );
// TH1F* h01 = (TH1F*) myhistos->Get("fit_eff_Et");
// TH1F* h11 = (TH1F*) myhistos->Get("fit_eff_Eta");
// h01->SetMarkerColor(4);
// h01->SetLineColor(4);
// h11->SetMarkerColor(4);
// h11->SetLineColor(4);
// ScaleTo10pb( *h01 );
// ScaleTo10pb( *h11 );
TH1F* h02 = (TH1F*) myhistos->Get("truth_eff_EtDet");
TH1F* h12 = (TH1F*) myhistos->Get("truth_eff_EtaDet");
h02->SetMarkerColor(2);
h02->SetLineColor(2);
h02->SetMarkerStyle(24);
h12->SetMarkerColor(2);
h12->SetLineColor(2);
h12->SetMarkerStyle(24);
ScaleTo10pb( *h02 );
ScaleTo10pb( *h12 );
TAxis* xaxis = h0->GetXaxis();
TAxis* yaxis = h0->GetYaxis();
xaxis->SetTitle("E_{T} (GeV)");
xaxis->SetTitleSize(0.04);
xaxis->SetTitleOffset(1.3);
xaxis->SetNdivisions(505);
yaxis->SetTitle("Efficiency ");
yaxis->SetTitleSize(0.04);
yaxis->SetTitleOffset(1.6);
h0->SetMinimum(0.5);
h0->SetMaximum(1.1);
TAxis* xaxis = h1->GetXaxis();
TAxis* yaxis = h1->GetYaxis();
xaxis->SetTitle("#eta");
xaxis->SetTitleSize(0.04);
xaxis->SetTitleOffset(1.3);
xaxis->SetNdivisions(505);
yaxis->SetTitle("Efficiency ");
yaxis->SetTitleSize(0.04);
yaxis->SetTitleOffset(1.6);
h1->SetMinimum(0.5);
h1->SetMaximum(1.1);
TAxis* xaxis = h2->GetXaxis();
TAxis* yaxis = h2->GetYaxis();
xaxis->SetTitle("E_{T} (GeV)");
xaxis->SetTitleSize(0.04);
xaxis->SetTitleOffset(1.3);
xaxis->SetNdivisions(505);
yaxis->SetTitle("#eta ");
yaxis->SetTitleSize(0.04);
yaxis->SetTitleOffset(1.7);
double l1 = 0.2, l2 =0.3;
double p1 = 0.5;
if(num==0) p1 = 0.82;
if(num==2) p1 = 0.4;
TCanvas* can0 = new TCanvas("can0", "", 500, 500);
gStyle->SetOptStat(0);
h0->Draw("e");
// h01->Draw("esame");
h02->Draw("esame");
leg_hist = new TLegend(0.3, l1, 0.8, l2);
leg_hist->AddEntry( h0,"Sideband subtracted","P");
// leg_hist->AddEntry( h01,"fit","le");
leg_hist->AddEntry( h02,"Monte Carlo truth","P");
leg_hist->SetFillColor(0);
leg_hist->SetMargin(0.15);
leg_hist->Draw();
TLatex* CMS = new TLatex();
CMS->SetTextAlign(12);
CMS->SetTextSize(0.04);
CMS->SetNDC();
CMS->DrawLatex(0.25, p1, "CMS Preliminary, #intLdt = 10 pb^{-1}");
can0->Update();
can0->SaveAs( label+TString("et.eps") );
can0->SaveAs( label+TString("et.gif") );
can0->SaveAs( label+TString("et.root") );
can0->SaveAs( label+TString("et.pdf") );
p1 = 0.5;
if(num==0 || num==2) p1 = 0.82;
TCanvas* can1 = new TCanvas("can1", "", 500, 500);
gStyle->SetOptStat(0);
h1->Draw("e");
// h11->Draw("esame");
h12->Draw("esame");
leg_hist = new TLegend(0.3,l1,0.8,l2);
leg_hist->AddEntry( h1,"Sideband subtracted","P");
// leg_hist->AddEntry( h11,"fit","le");
leg_hist->AddEntry( h12,"Monte Carlo truth","P");
leg_hist->SetFillColor(0);
leg_hist->SetMargin(0.15);
leg_hist->Draw();
TLatex* CMS2 = new TLatex();
CMS2->SetTextAlign(12);
CMS2->SetTextSize(0.04);
CMS2->SetNDC();
CMS2->DrawLatex(0.25,p1,"CMS Preliminary, #intLdt = 10 pb^{-1}");
can1->Update();
can1->SaveAs( label+TString("eta.eps") );
can1->SaveAs( label+TString("eta.gif") );
can1->SaveAs( label+TString("eta.root") );
can1->SaveAs( label+TString("eta.pdf") );
tdrStyle->SetPadRightMargin(0.2);
TCanvas* can2 = new TCanvas("can2", "", 500, 500);
gStyle->SetOptStat(0);
h2->Draw("colz");
can2->SaveAs( label+TString("et_eta.eps") );
can2->SaveAs( label+TString("et_eta.gif") );
can2->SaveAs( label+TString("et_eta.root") );
can2->SaveAs( label+TString("et_eta.pdf") );
delete can0;
delete can1;
delete can2;
delete myhistos;
}
/*
* Main Function
*/
void ratio() {
gStyle->SetFrameLineWidth(1);
//Number of Measurements
const int NUM = 2;
//Measurements and uncertainties
// mean, -stat, +stat, -syst, +syst
double m[NUM][5] = {
0.98, 0.10, 0.10, 0.14, 0.14,
//1.33, 0.32, 0.32, 0.22, 0.22
1.26, 0.37, 0.37, 0.46, 0.46
};
//Theory and uncertainties
// mean, -uncert, +uncert
double t[NUM][3] = {
1.0, 0.114, 0.114,
1.0, 0.111, 0.111
};
// label text, sub-label text
// Note: TString does not work, b/c one cannot pass an array of TStrings
// as an argument to a function
char label[NUM][2][100] = {
"Z#gamma#gamma", "",
//"W#gamma#gamma", ""
"W#gamma#gamma (#mu)", ""
};
// format:
// # color, bgColor, fontSytle, linewidth, markerStyle
int aux[NUM][5] = {
1, 10, 42, 2, 20,
1, 10, 42, 2, 20
};
// determning the x size of the plot
double lowX = LOW_X;
double uppX = UPP_X;
TH2F* lft = new TH2F("lft", "", 50, lowX, uppX, 1, 0.0, 1.0);
// height = NUM*unitHeight + 2*spacers + 1*bottomMargin + 0.5*topMargin
const double unitHeight = 50.0; // even number
const double height = (double(NUM)+2.0)*unitHeight+100.0+30.0;
// how much to step each time to cover the vertical range of the histo in
// exactly NUM+2 steps
const double vstep = unitHeight/(height-100.0-30.0);
const double width = 800.0;
printf("Canvas: width=%d, height=%d\n",
TMath::Nint(width), TMath::Nint(height));
printf("Y-step = %6.4f\n", vstep);
// Set canvas and margins
TCanvas* canvas = new TCanvas("canvas", "canvas", 200, 0,
TMath::Nint(width), TMath::Nint(height));
canvas->SetFillColor(10);
canvas->SetRightMargin(20.0/width);
canvas->SetLeftMargin(20.0/width);
canvas->SetBottomMargin(56.0/height);
canvas->SetTopMargin(30.0/height);
canvas->Draw();
canvas->cd();
//printf("TopMargin : %6.4f\n", canvas->GetTopMargin());
//printf("BottomMargin: %6.4f\n", canvas->GetBottomMargin());
TAxis* xaxis = lft->GetXaxis();
TAxis* yaxis = lft->GetYaxis();
xaxis->CenterTitle(kTRUE);
xaxis->SetTitleSize(0.07);
xaxis->SetTitleFont(62);
xaxis->SetTitleOffset(1.1);
xaxis->SetNdivisions(6,5,0);
xaxis->SetLabelOffset(0.01);
xaxis->SetLabelSize(0.05);
xaxis->SetLabelFont(42);
yaxis->SetLabelSize(0.0);
yaxis->SetNdivisions(-1);
lft->SetXTitle("Cross Section Ratio #sigma_{Exp} / #sigma_{Theory}");
lft->SetYTitle("");
lft->SetStats(kFALSE);
lft->SetTitle("");
lft->Draw();
// Draw Theory Bands
for (int i=0; i!=NUM; ++i) {
drawTheory(i, t[i][0], t[i][1], t[i][2], vstep);
}
canvas->RedrawAxis();
gPad->RedrawAxis();
// Draw Measurements
for (int i=0; i!=NUM; ++i) {
drawMeasurement(i, m[i], label[i], aux[i], vstep, lft, canvas);
}
//Set CMS Preliminary marker
int iPeriod=2; // 8 TeV
int iPos = 0; //
CMS_lumi( canvas, iPeriod, iPos );
// Print
//canvas->Print("Ratio_CrossSec_Exp_Theory.pdf");
canvas->Draw();
}
void draw() {
gROOT->Reset();
// first make a canvas and a 2D histogram for the axes
TCanvas* canvas = new TCanvas("canvas", "canvas", 10, 10, 500, 500);
canvas->SetFillColor(0);
canvas->SetBorderMode(0);
canvas->SetFrameBorderMode(0); // need this to turn off red hist frame!
gROOT->SetStyle("Plain");
canvas->UseCurrentStyle();
gPad->SetLeftMargin(0.15);
gPad->SetRightMargin(0.05);
gPad->SetTopMargin(0.07);
gPad->SetBottomMargin(0.17);
gStyle->SetOptStat(0);
gStyle->SetTitleBorderSize(0);
gStyle->SetTitleSize(0.04);
gStyle->SetTextFont(42);
gStyle->SetTextSize(0.04);
gStyle->SetTitleFont(42, "hxy"); // for histogram and axis title
gStyle->SetLabelFont(42, "xyz"); // for axis labels (values)
gStyle->SetTitleOffset(0.8, "h"); // what does this do?
gStyle->SetTitleX(0.15);
gStyle->SetTitleY(0.99);
gROOT->ForceStyle();
// can make histogram or alternatively use the histograms automatically
// connected to the TF1 or TGraph objects
double xMin = 0.1;
double xMax = 100.;
// double yMin = 0.1;
// double yMax = 10.;
double yMin = 0.;
double yMax = 8.;
TH2F* axhist = new TH2F("axhist", "title", 10, xMin, xMax, 10, yMin, yMax);
axhist->SetTitle("");
axhist->SetXTitle("b");
axhist->SetYTitle("med[Z|s]");
gPad->SetLogx(1);
gPad->SetLogy(0);
double u[20];
double x[20][500];
// Read in data from file and insert in TTree
TString fileName;
// cout << "Enter file name: ";
// cin >> fileName;
fileName = "medsig_s5_rel_bi.txt";
ifstream inFile;
inFile.open(fileName);
if (inFile.fail()) {
cout << "Couldn't open file!" << endl;
exit(1);
}
bool readLine = true;
int lineNum = 0;
int ncol;
while ( readLine ){
TString line;
stringstream ss;
line.ReadLine(inFile);
readLine = inFile.good();
if ( readLine ) {
TString firstChar = line(0,1);
bool useLine = firstChar != "#";
if ( useLine ){
int i = 0;
stringstream ss;
ss << line; // put whole line into ss
TString token;
bool getToken = true;
while ( getToken ) {
ss >> token; // extracts one token
if ( token.Length() > 0 ) {
u[i] = token.Atof();
i++;
}
else {
getToken = false;
}
} // getToken
ncol = i;
for (int i=0; i<ncol; i++){
x[i][lineNum] = u[i];
}
lineNum++;
} // useLine
} // readLine
} // readLine
int n = lineNum;
inFile.close();
// for (int i=0; i<n; i++){
// cout << i << " " << x[0][i] << " " << x[3][i] << " "
// << x[4][i] << endl;
// }
TGraph* tg1 = new TGraph(n, x[0], x[1]);
TGraph* tg2 = new TGraph(n, x[0], x[2]);
TGraph* tg3 = new TGraph(n, x[0], x[3]);
TGraph* tg4 = new TGraph(n, x[0], x[4]);
TGraph* tg5 = new TGraph(n, x[0], x[5]);
TGraph* tg6 = new TGraph(n, x[0], x[6]);
TGraph* tg7 = new TGraph(n, x[0], x[7]);
TGraph* tg8 = new TGraph(n, x[0], x[8]);
TGraph* tg9 = new TGraph(n, x[0], x[9]);
TGraph* tg10 = new TGraph(n, x[0], x[10]);
TGraph* tg11 = new TGraph(n, x[0], x[11]);
TGraph* tg12 = new TGraph(n, x[0], x[12]);
TGraph* tg13 = new TGraph(n, x[0], x[13]);
TGraph* tg14 = new TGraph(n, x[0], x[14]);
TGraph* tg15 = new TGraph(n, x[0], x[15]);
TAxis* xa = axhist->GetXaxis();
TAxis* ya = axhist->GetYaxis();
xa->SetTitleOffset(1.2); // factor multiplies default offset
ya->SetTitleOffset(1.1);
xa->SetLabelOffset(0.005);
ya->SetLabelOffset(0.005);
xa->SetTickLength(0.015); // default = 0.03
ya->SetTickLength(0.015); // default = 0.03
xa->SetTitleSize(0.05);
ya->SetTitleSize(0.05);
// gPad->SetLogx(1);
// xa->SetLimits(90., 700.);
xa->SetNdivisions(-5); // negative value should force number of divisions?
ya->SetNdivisions(-4);
xa->SetLabelSize(0.05);
ya->SetLabelSize(0.05);
// Draw axes and then add stuff
// kDot=1, kPlus, kStar, kCircle=4, kMultiply=5,
// kFullDotSmall=6, kFullDotMedium=7, kFullDotLarge=8,
// kFullCircle=20, kFullSquare=21, kFullTriangleUp=22,
// kFullTriangleDown=23, kOpenCircle=24, kOpenSquare=25,
// kOpenTriangleUp=26, kOpenDiamond=27, kOpenCross=28,
// kFullStar=29, kOpenStar=30
axhist->Draw();
tg1->SetLineColor(kRed);
tg1->SetLineWidth(2);
tg1->SetLineStyle(2);
tg1->SetMarkerColor(kRed);
tg1->SetMarkerSize(0.8);
tg1->SetMarkerStyle(20);
tg1->Draw("L,same"); // or P for points
tg2->SetLineColor(kRed);
tg2->SetLineWidth(2);
tg2->SetLineStyle(2);
tg2->SetMarkerColor(kRed);
tg2->SetMarkerSize(0.8);
tg2->SetMarkerStyle(20);
tg2->Draw("L,same"); // or P for points
tg3->SetLineColor(kRed);
tg3->SetLineWidth(2);
tg3->SetLineStyle(2);
tg3->SetMarkerColor(kRed);
tg3->SetMarkerSize(0.8);
tg3->SetMarkerStyle(20);
// tg3->Draw("L,same"); // or P for points
tg4->SetLineColor(kBlue);
tg4->SetLineWidth(2);
tg4->SetLineStyle(1);
tg4->SetMarkerColor(kBlue);
tg4->SetMarkerSize(0.8);
tg4->SetMarkerStyle(20);
tg4->Draw("L,same"); // or P for points
tg5->SetLineColor(kBlue);
tg5->SetLineWidth(2);
tg5->SetLineStyle(1);
tg5->SetMarkerColor(kBlue);
tg5->SetMarkerSize(0.8);
tg5->SetMarkerStyle(20);
tg5->Draw("L,same"); // or P for points
tg6->SetLineColor(kBlue);
tg6->SetLineWidth(2);
tg6->SetLineStyle(1);
tg6->SetMarkerColor(kBlue);
tg6->SetMarkerSize(0.8);
tg6->SetMarkerStyle(20);
// tg6->Draw("L,same"); // or P for points
tg7->SetLineColor(kRed);
tg7->SetLineWidth(2);
tg7->SetLineStyle(1);
tg7->SetMarkerColor(kRed);
tg7->SetMarkerSize(0.8);
tg7->SetMarkerStyle(21);
// tg7->Draw("P,same"); // or P for points
tg8->SetLineColor(kRed);
tg8->SetLineWidth(2);
tg8->SetLineStyle(1);
tg8->SetMarkerColor(kRed);
tg8->SetMarkerSize(0.8);
tg8->SetMarkerStyle(21);
// tg8->Draw("P,same"); // or P for points
tg9->SetLineColor(kRed);
tg9->SetLineWidth(2);
tg9->SetLineStyle(2);
tg9->SetMarkerColor(kRed);
tg9->SetMarkerSize(0.8);
tg9->SetMarkerStyle(21);
// tg9->Draw("P,same"); // or P for points
tg10->SetLineColor(kBlack);
tg10->SetLineWidth(2);
tg10->SetLineStyle(2);
tg10->SetMarkerColor(kBlack);
tg10->SetMarkerSize(0.8);
tg10->SetMarkerStyle(20);
tg10->Draw("P,same"); // or P for points
tg11->SetLineColor(kBlack);
tg11->SetLineWidth(2);
tg11->SetLineStyle(2);
tg11->SetMarkerColor(kBlack);
tg11->SetMarkerSize(0.8);
tg11->SetMarkerStyle(20);
tg11->Draw("P,same"); // or P for points
tg12->SetLineColor(kBlack);
tg12->SetLineWidth(2);
tg12->SetLineStyle(2);
tg12->SetMarkerColor(kBlack);
tg12->SetMarkerSize(0.8);
tg12->SetMarkerStyle(20);
// tg12->Draw("P,same"); // or P for points
tg13->SetLineColor(kBlack);
tg13->SetLineWidth(2);
tg13->SetLineStyle(3);
tg13->SetMarkerColor(kBlack);
tg13->SetMarkerSize(0.8);
tg13->SetMarkerStyle(20);
// tg13->Draw("L,same"); // or P for points
tg14->SetLineColor(kBlack);
tg14->SetLineWidth(2);
tg14->SetLineStyle(3);
tg14->SetMarkerColor(kBlack);
tg14->SetMarkerSize(0.8);
tg14->SetMarkerStyle(20);
// tg14->Draw("L,same"); // or P for points
tg15->SetLineColor(kBlack);
tg15->SetLineWidth(2);
tg15->SetLineStyle(3);
tg15->SetMarkerColor(kBlack);
tg15->SetMarkerSize(0.8);
tg15->SetMarkerStyle(20);
// tg15->Draw("L,same"); // or P for points
TLegend* leg = new TLegend(0.53, 0.48, 0.95, .73); // x1, y1, x2, y2
leg->SetTextSize(0.05);
leg->SetTextFont(42);
leg->SetBorderSize(0);
leg->SetFillColor(0);
leg->AddEntry(tg1, " s / #sqrt{b + #sigma_{b}^{2}}", "l");
leg->AddEntry(tg4, " Z_{A}", "l");
// leg->AddEntry(tg7, " #sqrt{q0}, MC median", "p");
leg->AddEntry(tg10, " Monte Carlo", "p");
// leg->AddEntry(tg13, " Z_{bi}" , "l");
// leg->AddEntry(tg4, "s = 0.03", "l");
leg->Draw();
TLatex* tl = new TLatex();
tl->SetTextAlign(11);
tl->SetTextSize(0.05);
tl->SetTextFont(42);
tl->SetNDC();
// tl->DrawLatex(.76, 0.33, "s = 2");
// tl->DrawLatex(.76, 0.465, "s = 5");
// tl->DrawLatex(.76, 0.6, "s = 10");
// tl->DrawLatex(.76, 0.73, "s = 20");
// tl->DrawLatex(.5, .77, "#sigma_{b}/b = 0.2, 0.5, 1");
tl->DrawLatex(.5, .77, "#sigma_{b}/b = 0.2, 0.5");
tl->DrawLatex(.5, .85, "s = 5");
// Fix idiotic problem with frame
TLine* tli = new TLine();
tli->SetLineStyle(1);
tli->SetLineWidth(1);
tli->DrawLine(xMin, yMin, xMax, yMin);
tli->DrawLine(xMax, yMin, xMax, yMax);
tli->DrawLine(xMin, yMax, xMax, yMax);
tli->DrawLine(xMin, yMin, xMin, yMax);
TPostScript psfile("medsig_s5_rel_bi.eps", 113); // 113 makes eps
canvas->Draw();
psfile.Close();
// canvas->Print("plot.gif", "gif");
}
void plot_ZJetBalance( int index, std::string JetResponseFilename,
Float_t genMean[], Float_t genSigma[],
Float_t recoMean[], Float_t recoSigma[]) {
////////////// Defining the L3JetCorrector ///////////////////////
double p[6];
p[0] = 10.0;
p[1] = 800.0;
// // with MEAN response
// p[2] = 0.953124;
// p[3] = 4.87151;
// p[4] = 2.83723;
// p[5] = 2.91468;
// // with MPV response
// p[2] = 0.978755;
// p[3] = 2.1759;
// p[4] = 2.25031;
// p[5] = 0.0;
// // with TruncMEAN 2 sigma response
// p[2] = 1.00299;
// p[3] = 3.83695;
// p[4] = 2.87351;
// p[5] = 1.6071;
// // with TruncMEAN 1.5 sigma response
p[2] = 1.00299;
p[3] = 3.83695;
p[4] = 2.87351;
p[5] = 1.6071;
// // with TruncMEAN 1.5 sigma iterative 3 times
// p[2] = 1.16793;
// p[3] = 8.8575;
// p[4] = 5.01816;
// p[5] = 5.03132;
// // with TruncMEAN 1.5 sigma iterative n times
// p[2] = 1.23392;
// p[3] = 9.33578;
// p[4] = 5.60468;
// p[5] = 5.52937;
// // with TruncMEAN 1 sigma response
// p[2] = 0.690152;
// p[3] = 2.97111;
// p[4] = 1.70586;
// p[5] = 0.72181;
// // with MPV: when p[5] is allowed to go negative
// p[2] = 0.893141;
// p[3] = 0.00355417;
// p[4] = 0.0123656;
// p[5] = -1.00176;
// // from dijet MC truth
// p[2] = 0.996998;
// p[3] = 4.39412;
// p[4] = 2.96134;
// p[5] = 1.69966;
if(storeResponsHistogramsInRootFile == true) {
TFile respHistFile(JetResponseFilename.c_str(),"RECREATE");
}
TH1F* responseHistGen[nZPtBins];
TH1F* responseHistReco[nZPtBins];
TH1F* genJetpt[nZPtBins];
TH1F* recoZpt[nZPtBins];
TH1F* caloJetpt[nZPtBins];
TH1F* Zpt[nZPtBins];
for(int i=0; i<nZPtBins; i++) {
responseHistGen[i] = new TH1F(TString("responseHistGen_")+ptBin[i],
"",40, 0.0, 2.0);
TAxis* responseHistGenx = responseHistGen[i]->GetXaxis();
TAxis* responseHistGeny = responseHistGen[i]->GetYaxis();
responseHistGenx->SetTitle("p_{T}^{jet} / p_{T}^{Z} ");
responseHistGeny->SetTitle("Events / 0.05");
responseHistGeny->SetTitleOffset(1.4);
responseHistGenx->SetNdivisions(505);
responseHistGeny->SetNdivisions(505);
responseHistReco[i] = new TH1F(TString("responseHistReco_")+ptBin[i],
"",40,0.0,2.0);
responseHistReco[i]->SetLineColor(2);
responseHistReco[i]->SetMarkerColor(2);
genJetpt[i] = new TH1F(TString("genJetpt_")+ptBin[i],
"", 140, 0, 800);
caloJetpt[i] = new TH1F(TString("caloJetpt_")+ptBin[i],
"", 140, 0, 800);
recoZpt[i] = new TH1F(TString("recoZpt_")+ptBin[i],
"", 140, 0, 800);
responseHistReco[i]->Sumw2();
responseHistGen[i]->Sumw2();
genJetpt[i]->Sumw2();
caloJetpt[i]->Sumw2();
recoZpt[i]->Sumw2();
}
for (Long64_t entry =0; entry < mychain->GetEntries(); entry++) {
mychain->GetEntry(entry);
if(entry%100000==0) std::cout<<"**** Event # "<< entry <<std::endl;
// Fill generator level quantities
int leadGenIndex=-1, secondGenIndex=-1;
FindLeadIndex(JetGenPt[index], JetGenEta[index], JetGenPhi[index],
eMinusEtaGen, eMinusPhiGen,
ePlusEtaGen, ePlusPhiGen, leadGenIndex, secondGenIndex);
if(leadGenIndex !=-1 && secondGenIndex !=-1 &&
(ePlusPtGen>20.0) && (eMinusPtGen>20.0) &&
((fabs(ePlusEtaGen)<1.4442) ||
(fabs(ePlusEtaGen)>1.560 && fabs(ePlusEtaGen)<2.5)) &&
((fabs(eMinusEtaGen)<1.4442) ||
(fabs(eMinusEtaGen)>1.560 && fabs(eMinusEtaGen)<2.5)) &&
(fabs(JetGenEta[index][leadGenIndex])<1.3) ) {
float leadGenJetPt = JetGenPt[index][leadGenIndex];
float secondGenJetPt = JetGenPt[index][secondGenIndex];
double wt = GetWeight( mychain->GetFile()->GetName() );
double ptRatioGen = (double) leadGenJetPt/ (double) Z_PtGen;
if((double) secondGenJetPt/(double) Z_PtGen < 0.1) {
for(int bin=0; bin<nZPtBins; bin++) { //begin Z pT bin loop
if( (Z_PtGen > theZPt[bin]) && (Z_PtGen < theZPt[bin+1]) ) {
responseHistGen[bin]->Fill( ptRatioGen, wt );
genJetpt[bin]->Fill(leadGenJetPt, wt );
} // end if loop
} // end Z pT loop
}
}
// Fill reco level quantities
int leadRecoIndex=-1, secondRecoIndex=-1;
FindLeadIndex(JetRecoPt[index], JetRecoEta[index], JetRecoPhi[index],
eMinusEta, eMinusPhi,
ePlusEta, ePlusPhi, leadRecoIndex, secondRecoIndex);
if(leadRecoIndex !=-1 && secondRecoIndex !=-1) {
float leadRecoJetPt;
float secondRecoJetPt;
float leadRecoJetEta;
double dPhiReco;
double L3scale = 1.0;
double apt = JetRecoPt[index][leadRecoIndex];
if( apt< p[0] ) apt = p[0];
if( apt> p[1] ) apt = p[1];
double log10pt = log10(apt);
double result = p[2]+p[3]/(pow(log10pt,p[4])+p[5]);
if(usingCorrectedCaloJetPt) L3scale = result;
leadRecoJetPt = L3scale * JetRecoPt[index][leadRecoIndex];
secondRecoJetPt = JetRecoPt[index][secondRecoIndex];
leadRecoJetEta = JetRecoEta[index][leadRecoIndex];
dPhiReco = dPhi(JetRecoPhi[index][leadRecoIndex], Z_Phi);
bool pass = BoolCutResult( Z_Pt, mZee, eMinusPt, eMinusEta, ePlusPt,
ePlusEta, eMinus_trackiso, ePlus_trackiso,
eMinus_ecaliso, ePlus_ecaliso,
eMinus_hcaliso, ePlus_hcaliso,
iseMinusLoose, isePlusLoose,
leadRecoJetPt, leadRecoJetEta,
secondRecoJetPt, dPhiReco);
if( ! pass ) continue;
double wt = GetWeight( mychain->GetFile()->GetName() );
double ptRatioReco = (double) leadRecoJetPt/ (double) Z_Pt;
for(int bin=0; bin<nZPtBins; bin++) { //begin Z pT bin loop
if( (Z_Pt > theZPt[bin]) && (Z_Pt < theZPt[bin+1]) ) {
responseHistReco[bin]->Fill( ptRatioReco, wt );
caloJetpt[bin]->Fill(leadRecoJetPt, wt );
recoZpt[bin]->Fill( Z_Pt, wt );
} // end if loop
} // end Z pT loop
}
} // end TTree loop
// Fill the mean and error vectors
for(int i=0; i<nZPtBins; i++) {
if(ScaleToLuminosity) {
for(int j=0; j<40; j++) {
double err = sqrt( responseHistReco[i]->GetBinContent(j));
responseHistReco[i]->SetBinError(j, err);
err = sqrt( responseHistGen[i]->GetBinContent(j));
responseHistGen[i]->SetBinError(j, err);
}
}
genMean[i] = (Float_t) responseHistGen[i]->GetMean(1);
genSigma[i] = (Float_t) responseHistGen[i]->GetMean(11);
recoMean[i] = (Float_t) responseHistReco[i]->GetMean(1);
recoSigma[i] = (Float_t) responseHistReco[i]->GetMean(11);
}
// plot the pT balance response histograms
if(makeplot_ZptBalance==true) {
for(int i=0; i<nZPtBins; i++) {
TString plotname = Form("ptBalance-allCuts_%d_%d", (int) theZPt[i],
(int) theZPt[i+1] );
PlotOnCanvas( *responseHistGen[i], *responseHistReco[i], plotname);
}
}
// Now write all the histograms in a ROOT file
if(storeResponsHistogramsInRootFile == true) {
respHistFile.cd();
for(int i=0; i<nZPtBins; i++) {
responseHistGen[i]->Write();
responseHistReco[i]->Write();
genJetpt[i]->Write();
caloJetpt[i]->Write();
recoZpt[i]->Write();
}
respHistFile.Close();
}
// clean up the memory
delete [] responseHistGen;
delete [] responseHistReco;
delete [] genJetpt;
delete [] caloJetpt;
delete [] recoZpt;
}
void makePlot(const std::string& inputFilePath, const std::string& canvasName, const std::string& histogram, const std::string& channel, const std::string& xAxisTitle, const std::string& yAxisTitle,
const std::string& inputFileName, const std::string& outputFilePath, const std::string& outputFileName)
{
std::string inputFileName_full = Form("%s%s", inputFilePath.data(), inputFileName.data());
TFile* inputFile = new TFile(inputFileName_full.data());
if ( !inputFile ) {
std::cerr << "Failed to open input file = " << inputFileName_full << " !!" << std::endl;
assert(0);
}
inputFile->ls();
TCanvas* canvas = dynamic_cast<TCanvas*>(inputFile->Get(canvasName.data()));
if ( !canvas ) {
std::cerr << "Failed to load canvas = " << canvasName << " !!" << std::endl;
assert(0);
}
int idxPad = -1;
if ( histogram == "mVis" ) idxPad = 1;
if ( histogram == "mTauTau" ) idxPad = 2;
if ( !(idxPad >= 1 && idxPad <= 2) ) {
std::cerr << "Invalid histogram = " << histogram << " !!" << std::endl;
assert(0);
}
TVirtualPad* pad = canvas->GetPad(idxPad);
std::cout << "pad = " << pad << ": ClassName = " << pad->ClassName() << std::endl;
TCanvas* canvas_new = new TCanvas("canvas_new", "canvas_new", 900, 850);
canvas_new->SetFillColor(10);
canvas_new->SetBorderSize(2);
canvas_new->SetTopMargin(0.065);
canvas_new->SetLeftMargin(0.17);
canvas_new->SetBottomMargin(0.155);
canvas_new->SetRightMargin(0.045);
canvas_new->SetLogx(false);
canvas_new->SetLogy(false);
canvas_new->Draw();
canvas_new->cd();
//TList* pad_primitives = canvas->GetListOfPrimitives();
TList* pad_primitives = pad->GetListOfPrimitives();
TH1* histogramDYJets = 0;
TH1* histogramHiggs125 = 0;
TH1* histogramHiggs200 = 0;
TH1* histogramHiggs300 = 0;
TIter pad_nextObj(pad_primitives);
while ( TObject* obj = pad_nextObj() ) {
std::string objName = "";
if ( dynamic_cast<TNamed*>(obj) ) objName = (dynamic_cast<TNamed*>(obj))->GetName();
std::cout << "obj = " << obj << ": name = " << objName << ", type = " << obj->ClassName() << std::endl;
TH1* tmpHistogram = dynamic_cast<TH1*>(obj);
if ( tmpHistogram ) {
std::cout << "tmpHistogram:"
<< " fillColor = " << tmpHistogram->GetFillColor() << ", fillStyle = " << tmpHistogram->GetFillStyle() << ","
<< " lineColor = " << tmpHistogram->GetLineColor() << ", lineStyle = " << tmpHistogram->GetLineStyle() << ", lineWidth = " << tmpHistogram->GetLineWidth() << ","
<< " markerColor = " << tmpHistogram->GetMarkerColor() << ", markerStyle = " << tmpHistogram->GetMarkerStyle() << ", markerSize = " << tmpHistogram->GetMarkerSize() << ","
<< " integral = " << tmpHistogram->Integral() << std::endl;
if ( tmpHistogram->GetFillColor() == 0 ) histogramDYJets = tmpHistogram;
if ( tmpHistogram->GetFillColor() == 632 ) histogramHiggs125 = tmpHistogram;
if ( tmpHistogram->GetFillColor() == 616 ) histogramHiggs200 = tmpHistogram;
if ( tmpHistogram->GetFillColor() == 600 ) histogramHiggs300 = tmpHistogram;
}
}
if ( !(histogramDYJets && histogramHiggs125 && histogramHiggs200 && histogramHiggs300) ) {
std::cerr << "Failed to load histograms !!" << std::endl;
assert(0);
}
int lineColors[4] = { 1, 2, 6, 4 };
int lineStyles[4] = { 1, 1, 1, 1 };
int lineWidths[4] = { 2, 2, 2, 2 };
int fillColors[4] = { 10, 2, 6, 4 };
int fillStyles[4] = { 0, 3002, 3004, 3005 };
histogramDYJets->SetFillColor(fillColors[0]);
histogramDYJets->SetFillStyle(fillStyles[0]);
histogramDYJets->SetLineColor(lineColors[0]);
histogramDYJets->SetLineStyle(lineStyles[0]);
histogramDYJets->SetLineWidth(lineWidths[0]);
histogramHiggs125->SetFillColor(fillColors[1]);
histogramHiggs125->SetFillStyle(fillStyles[1]);
histogramHiggs125->SetLineColor(lineColors[1]);
histogramHiggs125->SetLineStyle(lineStyles[1]);
histogramHiggs125->SetLineWidth(lineWidths[1]);
histogramHiggs200->SetFillColor(fillColors[2]);
histogramHiggs200->SetFillStyle(fillStyles[2]);
histogramHiggs200->SetLineColor(lineColors[2]);
histogramHiggs200->SetLineStyle(lineStyles[2]);
histogramHiggs200->SetLineWidth(lineWidths[2]);
histogramHiggs300->SetFillColor(fillColors[3]);
histogramHiggs300->SetFillStyle(fillStyles[3]);
histogramHiggs300->SetLineColor(lineColors[3]);
histogramHiggs300->SetLineStyle(lineStyles[3]);
histogramHiggs300->SetLineWidth(lineWidths[3]);
TAxis* xAxis = histogramHiggs300->GetXaxis();
if ( histogram == "mVis" ) xAxis->SetRangeUser(0,350);
else xAxis->SetRangeUser(0,450);
xAxis->SetTitle(xAxisTitle.data());
xAxis->SetTitleOffset(1.00);
xAxis->SetTitleSize(0.070);
xAxis->SetTitleFont(42);
xAxis->SetLabelOffset(0.010);
xAxis->SetLabelSize(0.050);
xAxis->SetLabelFont(42);
xAxis->SetTickLength(0.040);
xAxis->SetNdivisions(505);
//double xMin = 20.;
//double xMax = xAxis->GetXmax();
//xAxis->SetRangeUser(xMin, xMax);
TAxis* yAxis = histogramHiggs300->GetYaxis();
yAxis->SetTitle(yAxisTitle.data());
yAxis->SetTitleOffset(1.20);
yAxis->SetTitleSize(0.070);
yAxis->SetTitleFont(42);
yAxis->SetLabelOffset(0.010);
yAxis->SetLabelSize(0.055);
yAxis->SetLabelFont(42);
yAxis->SetTickLength(0.040);
yAxis->SetNdivisions(505);
histogramHiggs300->SetTitle("");
histogramHiggs300->SetStats(false);
histogramHiggs300->SetMaximum(1.2*histogramDYJets->GetMaximum());
histogramHiggs300->SetMinimum(0.);
histogramHiggs300->Draw("hist");
histogramHiggs200->Draw("histsame");
histogramHiggs125->Draw("histsame");
histogramDYJets->Draw("histsame");
histogramHiggs300->Draw("axissame");
TLegend* legend_new = new TLegend(0.50, 0.62, 0.87, 0.92, NULL, "brNDC");
legend_new->SetFillColor(10);
legend_new->SetFillStyle(0);
legend_new->SetBorderSize(0);
legend_new->SetTextFont(42);
legend_new->SetTextSize(0.055);
legend_new->SetTextColor(1);
legend_new->SetMargin(0.20);
legend_new->AddEntry(histogramDYJets, "Z/#gamma* #rightarrow #tau#tau", "f");
legend_new->AddEntry(histogramHiggs125, "H(125 GeV) #rightarrow #tau#tau", "f");
legend_new->AddEntry(histogramHiggs200, "H(200 GeV) #rightarrow #tau#tau", "f");
legend_new->AddEntry(histogramHiggs300, "H(300 GeV) #rightarrow #tau#tau", "f");
legend_new->Draw();
double label_channel_y0;
if ( channel == "e#mu" ) label_channel_y0 = 0.9275;
else if ( channel == "#mu#tau_{h}" ) label_channel_y0 = 0.9400;
else if ( channel == "#tau_{h}#tau_{h}" ) label_channel_y0 = 0.9350;
else {
std::cerr << "Invalid channel = " << channel << " !!" << std::endl;
assert(0);
}
TPaveText* label_channel = new TPaveText(0.855, label_channel_y0, 0.945, label_channel_y0 + 0.055, "NDC");
label_channel->SetFillStyle(0);
label_channel->SetBorderSize(0);
label_channel->AddText(channel.data());
label_channel->SetTextFont(62);
label_channel->SetTextSize(0.055);
label_channel->SetTextColor(1);
label_channel->SetTextAlign(31);
label_channel->Draw();
canvas_new->Update();
std::string outputFileName_full = Form("%s%s", outputFilePath.data(), outputFileName.data());
size_t idx = outputFileName_full.find_last_of('.');
std::string outputFileName_plot = std::string(outputFileName_full, 0, idx);
canvas_new->Print(std::string(outputFileName_plot).append(".pdf").data());
//canvas_new->Print(std::string(outputFileName_plot).append(".root").data());
delete legend_new;
delete label_channel;
delete canvas_new;
delete inputFile;
}
//------------------------------------------------------------------------------
// DrawZPeak
//------------------------------------------------------------------------------
void DrawZPeak(TString energy)
{
if (energy.Contains("7TeV"))
{
_lumiText = "4.9 fb^{-1} (7 TeV)";
}
else
{
_lumiText = "19.6 fb^{-1} (8 TeV)";
}
// Read the input file
//----------------------------------------------------------------------------
TString name = "invMass2Lep_";
TFile* file = new TFile("rootfiles/" + name + energy + ".root", "read");
TH1F* data;
TH1F* WZ;
TH1F* fakes;
TH1F* ZZ;
TH1F* Zgamma;
TH1F* WV;
TH1F* VVV;
TH1F* allmc;
if (energy.Contains("7TeV"))
{
data = (TH1F*)file->Get("hZMass_Sel_datahist_COMB");
WZ = (TH1F*)file->Get("hZMass_Sel_wz_COMB");
fakes = (TH1F*)file->Get("hZMass_Sel_datadriven_COMB");
ZZ = (TH1F*)file->Get("hZMass_Sel_zz_COMB");
Zgamma = (TH1F*)file->Get("hZMass_Sel_zg_COMB");
WV = (TH1F*)file->Get("hZMass_Sel_wz_COMB");
VVV = (TH1F*)file->Get("hZMass_Sel_wz_COMB");
allmc = (TH1F*)file->Get("all_estimates_with_error_COMB");
}
else if (energy.Contains("8TeV"))
{
data = (TH1F*)file->Get("h_data");
WZ = (TH1F*)file->Get("h_WZ");
fakes = (TH1F*)file->Get("h_Fakes");
ZZ = (TH1F*)file->Get("h_ZZ");
Zgamma = (TH1F*)file->Get("h_ZGamma");
WV = (TH1F*)file->Get("h_WV");
VVV = (TH1F*)file->Get("h_VVV");
allmc = (TH1F*)file->Get("h_All");
}
WZ->SetFillColor(kOrange-2);
WZ->SetLineColor(kOrange-2);
Zgamma->SetFillColor(kRed+1); // kRed+2
Zgamma->SetLineColor(kRed+1); // kRed+2
ZZ->SetFillColor(kRed+1);
ZZ->SetLineColor(kRed+1);
fakes->SetFillColor(kGray+1);
fakes->SetLineColor(kGray+1);
data->SetMarkerStyle(kFullCircle);
allmc->SetFillColor (kBlack);
allmc->SetFillStyle (3345);
allmc->SetLineColor (kWhite);
allmc->SetLineWidth (0);
allmc->SetMarkerColor(kOrange-2);
allmc->SetMarkerSize (0);
THStack* hs = new THStack();
if (energy.Contains("8TeV"))
{
WV->SetFillColor(kRed+1); // kAzure
WV->SetLineColor(kRed+1); // kAzure
VVV->SetFillColor(kRed+1); // kBlack
VVV->SetLineColor(kRed+1); // kBlack
hs->Add(VVV);
hs->Add(WV);
}
hs->Add(Zgamma);
hs->Add(ZZ);
hs->Add(fakes);
hs->Add(WZ);
// Draw
//----------------------------------------------------------------------------
TCanvas* canvas = new TCanvas(energy, energy);
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.3);
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->SetRangeUser(68, 112);
xaxis->SetTitle("m_{#font[12]{ll}} (GeV)");
yaxis->SetTitle(Form("Events / %.0f GeV", data->GetBinWidth(0)));
// Adjust scale
//----------------------------------------------------------------------------
Float_t theMax = GetMaximumIncludingErrors(data);
Float_t theMaxMC = GetMaximumIncludingErrors(allmc);
if (theMaxMC > theMax) theMax = theMaxMC;
data->SetMaximum(1.15 * theMax);
// Legend
//----------------------------------------------------------------------------
Double_t x0 = 0.635;
Double_t y0 = 0.770;
DrawTLegend(x0, y0 + 2.*(_yoffset+0.001), data, " Data", "ep");
DrawTLegend(x0, y0 + 1.*(_yoffset+0.001), WZ, " WZ", "f");
DrawTLegend(x0, y0, fakes, " Non-prompt leptons", "f");
DrawTLegend(x0, y0 - 1.*(_yoffset+0.001), ZZ, " MC background", "f");
DrawTLegend(x0, y0 - 2.*(_yoffset+0.001), allmc, " stat. #oplus syst.", "f");
// Finish it
//----------------------------------------------------------------------------
data->SetTitle("");
DrawTLatex(_cmsTextFont, 0.215, 0.880, 0.055, 13, "CMS");
// DrawTLatex(_extraTextFont, 0.215, 0.826, 0.030, 13, "Preliminary");
DrawTLatex(_lumiTextFont, 0.940, 0.940, 0.040, 31, _lumiText);
hs ->Draw("hist,same");
allmc->Draw("e2,same");
data ->Draw("ep,same");
canvas->GetFrame()->DrawClone();
canvas->RedrawAxis();
canvas->Update();
canvas->SaveAs("pdf/" + name + energy + ".pdf");
canvas->SaveAs("png/" + name + energy + ".png");
}
void plot3(TString infile = "fp-d", TString pltmd = "cos") {
// CHECK FOR RIGHT INPUT ////////////////////////////////////////////////
string strpltmd = pltmd, filename = infile, strfile = infile;
if( (strpltmd.compare("cos") != 0 ) &&
(strpltmd.compare("sin") != 0 ) &&
(strpltmd.compare("tan") != 0 ) &&
(strpltmd.compare("mmp") != 0 ) ) {error(4);};
// GLOBAL VARIABLES ////////////////////////////////////////////////////
Int_t file, point, color, style;
Float_t fits2b, fittph, tphold, fitsph, fitcph, fitx, fitxmin, fitxmax = -1.0;
Float_t xVal, yVal;
Float_t xMin = 100000, xMax = -1.0, yMin = 100000, yMax = -1.0;
Float_t MZ, MZmin, MW, MWmin, Mmin = 100000;
Float_t Cz1, Cz2, Cz3, Cw1, Cw2, Cw3, Cw4, C1, C2;
Float_t phiMin, phiMax, cphmin, cphmax, sphmin, sphmax;
// CUSTOMIZE PLOT ///////////////////////////////////////////////////////
gROOT->Reset();
gROOT->SetStyle("Plain");
gStyle->SetTitleBorderSize(0);
gStyle->SetPalette(1);
TCanvas *MyC = new TCanvas("MyC","Plot of the GAPP fit results",200,10,700,500);
Float_t mmlegxmin, mmlegxmax, mmlegymin, mmlegymax;
Float_t s2blegxmin, s2blegymin, s2blegxmax, s2blegymax;
Float_t lblxmin, lblxmax, lblymin, lblymax;
string plottitle = infile; //"Model: " + infile + " | Plot: ";
sToUpper(plottitle);
string xtitle, ytitle, NPleg, SMleg, display;
NPleg = "#font[52]{M_{H}^{(NP)}, #bar{m}_{t}^{(NP)}}";
SMleg = "#font[52]{M_{H}^{(SM)}, #bar{m}_{t}^{(SM)}}";
if (strpltmd.compare("tan") == 0) {
// plottitle += "#font[42]{tan^{2}(#tilde{#phi}) over }#font[52]{#tilde{x}}#font[42]{.}";
xtitle = "#font[52]{#tilde{x}}";
ytitle = "#font[42]{tan^{2}(#tilde{#phi})}";
display = "C";
mmlegxmin = 0.15;
mmlegxmax = 0.30;
mmlegymin = 0.75;
mmlegymax = 0.85;
s2blegxmin = 0.15;
s2blegxmax = 0.40;
s2blegymin = 0.30;
s2blegymax = 0.50;
lblxmin = 0.88;
lblxmax = 0.88;
lblymin = 0.60;
lblymax = 0.65;
} else if (strpltmd.compare("cos") == 0) {
// plottitle += "#font[42]{cos(#tilde{#phi}) over }#font[52]{#tilde{x}}#font[42]{.}";
xtitle = "#font[52]{#tilde{x}}";
ytitle = "#font[42]{cos(#tilde{#phi})}";
display = "C";
mmlegxmin = 0.65;
mmlegxmax = 0.80;
mmlegymin = 0.65;
mmlegymax = 0.75;
s2blegxmin = 0.60;
s2blegxmax = 0.85;
s2blegymin = 0.30;
s2blegymax = 0.50;
lblxmin = 0.35;
lblxmax = 0.50;
lblymin = 0.45;
lblymax = 0.50;
} else if (strpltmd.compare("sin") == 0) {
// plottitle += "#font[42]{sin(#tilde{#phi}) over }#font[52]{#tilde{x}}#font[42]{.}";
xtitle = "#font[52]{#tilde{x}}";
ytitle = "#font[42]{sin(#tilde{#phi})}";
display = "C";
mmlegxmin = 0.15;
mmlegxmax = 0.30;
mmlegymin = 0.75;
mmlegymax = 0.85;
s2blegxmin = 0.60;
s2blegxmax = 0.85;
s2blegymin = 0.30;
s2blegymax = 0.50;
lblxmin = 0.65;
lblxmax = 0.80;
lblymin = 0.60;
lblymax = 0.65;
} else if (strpltmd.compare("mmp") == 0) {
// plottitle += "#font[42]{Masses of the new heavy gauge bosons.}";
xtitle = "#font[52]{M_{Z'}}#font[42]{ (TeV)}";
ytitle = "#font[52]{M_{W'}}#font[42]{ (TeV)}";
display = "C";
mmlegxmin = 0.20;
mmlegxmax = 0.35;
mmlegymin = 0.45;
mmlegymax = 0.55;
s2blegxmin = 0.20;
s2blegxmax = 0.45;
s2blegymin = 0.60;
s2blegymax = 0.80;
lblxmin = 0.65;
lblxmax = 0.80;
lblymin = 0.30;
lblymax = 0.35;
};
// PREPARE BOSON MASSES AND PHI BOUNDS //////////////////////////////////
string mdl(filename,0,2);
if ( (mdl.compare("lr") == 0) ||
(mdl.compare("lp") == 0) ||
(mdl.compare("hp") == 0) ||
(mdl.compare("fp") == 0) ) {
phiMin = 5.600; phiMax = 84.400;
string Higgs(filename,3,1);
if (Higgs.compare("d") == 0) {
Cz1 = 11.95349795785275;
Cz2 = 30.63269990028513;
Cz3 = 42.58619785813789;
Cw1 = 21.29309892906894;
Cw2 = 9.339600971216193;
Cw3 = 30.63269990028513;
Cw4 = 42.58619785813789;
}
else if (Higgs.compare("t") == 0) {
Cz1 = 5.976748978926375;
Cz2 = 30.63269990028513;
Cz3 = 85.17239571627579;
Cw1 = 15.05649464522066;
Cw2 = 3.302047590161717;
Cw3 = 21.66058982554409;
Cw4 = 60.22597858088265;
}
}
else if ( (mdl.compare("uu") == 0) ||
(mdl.compare("nu") == 0) ) {
phiMin = 10.179, phiMax = 79.821;
C1 = 94.0397928463607;
C2 = 77.1253849720165;
} else {error(6);}
cphmin = cos(TMath::Pi()*phiMin/180.0)*cos(TMath::Pi()*phiMin/180.0);
cphmax = cos(TMath::Pi()*phiMax/180.0)*cos(TMath::Pi()*phiMax/180.0);
sphmin = sin(TMath::Pi()*phiMin/180.0)*sin(TMath::Pi()*phiMin/180.0);
sphmax = sin(TMath::Pi()*phiMax/180.0)*sin(TMath::Pi()*phiMax/180.0);
// LOOP OVER ROOT FILES ////////////////////////////////////////////////
for(file=0; file<=1; file++) {
if(file==0) string epsfile = filename + "_" + strpltmd + ".eps";
if(file==1) string filename = filename + "_sm";
string rootname = filename + ".root";
TFile *rootfile = TFile::Open(rootname.c_str());
if(rootfile == NULL) error(1);
TTree *tree = (TTree*)rootfile->Get(filename.c_str());
if(tree == NULL) error(2);
TBranch *fits2bbranch = (TBranch*)tree->GetBranch("fits2b");
TBranch *fittphbranch = (TBranch*)tree->GetBranch("fittph");
TBranch *fitxbranch = (TBranch*)tree->GetBranch("fitx");
if( (fits2bbranch == NULL) ||
(fittphbranch == NULL) ||
(fitxbranch == NULL) ) error(3);
tree->SetBranchAddress("fits2b",&fits2b);
tree->SetBranchAddress("fittph",&fittph);
tree->SetBranchAddress("fitx", &fitx);
// GET ARRAYS ///////////////////////////////////////////////////////////
Int_t Npoints = (Int_t)tree->GetEntries();
Int_t tphStep = 0;
Float_t tphMax = -1.0;
for(point=0; point<Npoints; point++) {
tree->GetEntry(point);
if( fittph > tphMax ) {tphStep++; tphMax = fittph;}
};
const int tphSteps = tphStep;
Float_t xArray[tphSteps], yArray[tphSteps], zArray[tphSteps];
tphStep = -1, tphold = -1.0, fitxmin = 100000;
for(point=0; point<Npoints; point++) {
tree->GetEntry(point);
if(fittph > tphold) {tphStep++; fitxmin = 100000;}
fitsph = fittph / (1.0 + fittph);
fitcph = 1.0 - fitsph;
if (strpltmd.compare("tan") == 0) {
xVal = fitx;
yVal = fittph;
} else if (strpltmd.compare("cos") == 0) {
xVal = fitx;
yVal = sqrt(fitcph);
} else if (strpltmd.compare("sin") == 0) {
xVal = fitx;
yVal = sqrt(fitsph);
} else if (strpltmd.compare("mmp") == 0) {
if (fitsph != 0.0) {
if ( (mdl.compare("lr") == 0) ||
(mdl.compare("lp") == 0) ||
(mdl.compare("hp") == 0) ||
(mdl.compare("fp") == 0) ) {
MMI(Cz1,Cz2,Cz3,Cw1,Cw2,Cw3,Cw4,fitx,fitsph,fitcph,fits2b,xVal,yVal);
} else if ( (mdl.compare("uu") == 0) ||
(mdl.compare("nu") == 0) ) {
MMII(C1,C2,fitx,fitsph,fitcph,xVal,yVal);
}
}
}
if( (strpltmd.compare("mmp") == 0) && (tphStep==1) ) {
xArray[0] = xArray[1];
yArray[0] = yArray[1];
zArray[0] = zArray[1];
}
if(fitx>fitxmax) fitxmax = fitx;
if(fitx<fitxmin) {
xArray[tphStep] = xVal;
yArray[tphStep] = yVal;
zArray[tphStep] = fits2b;
fitxmin = fitx;
}
tphold = fittph;
}
if(file==0) TGraph *NPplot = new TGraph(tphSteps,xArray,yArray);
if(file==1) TGraph *SMplot = new TGraph(tphSteps,xArray,yArray);
TMarker *NPmrk[tphSteps], *SMmrk[tphSteps];
MZmin = 100000, MWmin = 100000;
for(tphStep=0; tphStep<tphSteps; tphStep++){
marker(zArray[tphStep],color,style);
if(file==0) { NPmrk[tphStep] = new TMarker(xArray[tphStep],yArray[tphStep],style);
NPmrk[tphStep]->SetMarkerSize(0.8);
NPmrk[tphStep]->SetMarkerColor(color);}
if(file==1) { SMmrk[tphStep] = new TMarker(xArray[tphStep],yArray[tphStep],style);
SMmrk[tphStep]->SetMarkerSize(0.8);
SMmrk[tphStep]->SetMarkerColor(color);}
if( (strpltmd.compare("mmp") != 0) || (tphStep !=0 )) {
if (xArray[tphStep] < xMin) xMin = xArray[tphStep];
if (xArray[tphStep] > xMax) xMax = xArray[tphStep];
if (yArray[tphStep] < yMin) yMin = yArray[tphStep];
if (yArray[tphStep] > yMax) yMax = yArray[tphStep];
}
if( ((strfile.compare("uu-d") != 0) || (strfile.compare("nu-d") != 0)) && (strpltmd.compare("cos") == 0) ) {
fitx = xArray[tphStep];
fitcph = yArray[tphStep]*yArray[tphStep];
fitsph = 1.0 - fitcph;
if( (cphmax < fitcph) && (fitcph < cphmin) ) {
MMI(Cz1,Cz2,Cz3,Cw1,Cw2,Cw3,Cw4,fitx,fitsph,fitcph,fits2b,MZ,MW);
if(MZ < MZmin) MZmin = MZ;
if(MW < MWmin) MWmin = MW;
}
}
if( ((strfile.compare("uu-d") == 0) || (strfile.compare("nu-d") == 0)) && (strpltmd.compare("sin") == 0) ) {
fitx = xArray[tphStep];
fitsph = yArray[tphStep]*yArray[tphStep];
fitcph = 1.0 - fitsph;
if( (sphmin < fitsph) && (fitsph<sphmax) ) {
MMII(C1,C2,fitx,fitsph,fitcph,MZ,MW);
if(MZ < Mmin) { Mmin = MZ; cout << MZ << "\t" << sqrt(fitsph) << endl;}
}
}
}
cout << "(" << file << ") " << "MZmin: " << MZmin << endl;
cout << "(" << file << ") " << "MWmin: " << MWmin << endl;
}
// CREATE PLOTS /////////////////////////////////////////////////////////
NPplot->SetLineStyle(2);
NPplot->SetMarkerStyle(20);
NPplot->SetMarkerSize(0.4);
SMplot->SetMarkerStyle(20);
SMplot->SetMarkerSize(0.4);
if(strpltmd.compare("cos") == 0) {xMin = 0.0; yMin = 0.0; yMax = 1.0;}
if(strpltmd.compare("sin") == 0) {yMin = 0.0; yMax = 1.0;}
if(strpltmd.compare("mmp") == 0) {xMin = 0.0; xMax = 5.8/1.1; yMin = 0.0; yMax = 5;}
TH1F* frame = MyC->DrawFrame(0.9*xMin,0.9*yMin,1.1*xMax,1.0*yMax);
frame->SetTitle(plottitle.c_str());
TAxis *xaxis = frame->GetXaxis();
TAxis *yaxis = frame->GetYaxis();
xaxis->SetTitle(xtitle.c_str());
xaxis->CenterTitle();
xaxis->SetTitleOffset(1.);
xaxis->SetDecimals();
xaxis->SetLabelSize(0.03);
xaxis->SetLabelOffset(0.01);
yaxis->SetNdivisions(505);
yaxis->SetTitle(ytitle.c_str());
yaxis->CenterTitle();
yaxis->SetTitleOffset(1.2);
yaxis->SetDecimals();
yaxis->SetLabelSize(0.03);
yaxis->SetLabelOffset(0.01);
TLegend *mmleg = new TLegend(mmlegxmin,mmlegymin,mmlegxmax,mmlegymax);
mmleg->AddEntry(NPplot,NPleg.c_str(),"l");
mmleg->AddEntry(SMplot,SMleg.c_str(),"l");
mmleg->SetTextSize(0.025);
mmleg->SetFillStyle(0);
if( (strfile.compare("uu-d") != 0) && (strfile.compare("nu-d") != 0) ) {
for(tphStep=0; tphStep<tphSteps; tphStep++){NPmrk[tphStep]->Draw(); SMmrk[tphStep]->Draw();}
}
Float_t xdummy[1] = {0.0}, ydummy[1] = {0.0};
TGraph *circle = new TGraph(1,xdummy,ydummy);
circle->SetMarkerStyle(24);
circle->SetMarkerColor(kGreen+1);
circle->SetMarkerSize(0.8);
TGraph *square = new TGraph(1,xdummy,ydummy);
square->SetMarkerStyle(25);
square->SetMarkerColor(kCyan+1);
square->SetMarkerSize(0.8);
TGraph *triangle = new TGraph(1,xdummy,ydummy);
triangle->SetMarkerStyle(26);
triangle->SetMarkerColor(kBlue+1);
triangle->SetMarkerSize(0.8);
TGraph *diamond = new TGraph(1,xdummy,ydummy);
diamond->SetMarkerStyle(27);
diamond->SetMarkerColor(kMagenta+1);
diamond->SetMarkerSize(0.8);
TLegend *s2bleg = new TLegend(s2blegxmin,s2blegymin,s2blegxmax,s2blegymax);
s2bleg->AddEntry(circle,"#font[42]{0.00 < sin^{2}(2#tilde{#beta}) #leq 0.25}","p");
s2bleg->AddEntry(square,"#font[42]{0.25 < sin^{2}(2#tilde{#beta}) #leq 0.50}","p");
s2bleg->AddEntry(triangle,"#font[42]{0.50 < sin^{2}(2#tilde{#beta}) #leq 0.75}","p");
s2bleg->AddEntry(diamond,"#font[42]{0.75 < sin^{2}(2#tilde{#beta}) #leq 1.00}","p");
s2bleg->SetTextSize(0.025);
s2bleg->SetFillStyle(0);
NPplot->Draw(display.c_str());
SMplot->Draw(display.c_str());
mmleg->Draw();
if( (strfile.compare("uu-d") != 0) && (strfile.compare("nu-d") != 0) ) s2bleg->Draw();
// BOUNDS ON PHI //////////////////////////////////////////////////////
Int_t i;
const int iSteps = 100;
fitxmin = 1.0, fitxmax *= 1.5;
Float_t deltax = (fitxmax-fitxmin)/iSteps;
Float_t phixmin0[iSteps], phixmax0[iSteps], phixmin1[iSteps], phixmax1[iSteps];
Float_t phiymin0[iSteps], phiymax0[iSteps], phiymin1[iSteps], phiymax1[iSteps];
if ( (strpltmd.compare("tan") == 0) ||
(strpltmd.compare("cos") == 0) ||
(strpltmd.compare("sin") == 0) ) {
for(i=0; i<100; i++) {
fitx = fitxmin + i*deltax;
phixmin0[i] = fitx;
if (strpltmd.compare("tan") == 0) { phiymin0[i] = sphmin / cphmin; phiymax0[i] = sphmax / cphmax; }
if (strpltmd.compare("cos") == 0) { phiymin0[i] = sqrt(cphmin); phiymax0[i] = sqrt(cphmax); }
if (strpltmd.compare("sin") == 0) { phiymin0[i] = sqrt(sphmin); phiymax0[i] = sqrt(sphmax); }
}
TGraph *phiMin0 = new TGraph(iSteps,phixmin0,phiymin0);
TGraph *phiMax0 = new TGraph(iSteps,phixmin0,phiymax0);
} else if (strpltmd.compare("mmp") == 0) {
if ( (mdl.compare("lr") == 0) ||
(mdl.compare("lp") == 0) ||
(mdl.compare("hp") == 0) ||
(mdl.compare("fp") == 0) ) {
for(i=0; i<100; i++) {
fitx = fitxmin + i*deltax;
MMI(Cz1,Cz2,Cz3,Cw1,Cw2,Cw3,Cw4,fitx,sphmin,cphmin,0.0,phixmin0[i],phiymin0[i]);
MMI(Cz1,Cz2,Cz3,Cw1,Cw2,Cw3,Cw4,fitx,sphmin,cphmin,1.0,phixmin1[i],phiymin1[i]);
MMI(Cz1,Cz2,Cz3,Cw1,Cw2,Cw3,Cw4,fitx,sphmax,cphmax,0.0,phixmax0[i],phiymax0[i]);
MMI(Cz1,Cz2,Cz3,Cw1,Cw2,Cw3,Cw4,fitx,sphmax,cphmax,1.0,phixmax1[i],phiymax1[i]);
}
TGraph *phiMin0 = new TGraph(iSteps,phixmin0,phiymin0);
TGraph *phiMin1 = new TGraph(iSteps,phixmin1,phiymin1);
TGraph *phiMax0 = new TGraph(iSteps,phixmax0,phiymax0);
TGraph *phiMax1 = new TGraph(iSteps,phixmax1,phiymax1);
phiMin1->SetLineStyle(7);
phiMin1->SetMarkerStyle(22);
phiMin1->SetMarkerSize(1.0);
phiMax1->SetLineStyle(7);
phiMax1->SetMarkerStyle(22);
phiMax1->SetMarkerSize(1.0);
// phiMin1->Draw("C");
// phiMax1->Draw("C");
} else if ( (mdl.compare("uu") == 0) ||
(mdl.compare("nu") == 0) ) {
for(i=0; i<100; i++) {
fitx = fitxmin + i*deltax;
MMII(C1,C2,fitx,sphmin,cphmin,phixmin0[i],phiymin0[i]);
MMII(C1,C2,fitx,sphmax,cphmax,phixmax0[i],phiymax0[i]);
}
TGraph *phiMin0 = new TGraph(iSteps,phixmin0,phiymin0);
TGraph *phiMax0 = new TGraph(iSteps,phixmax0,phiymax0);
}
}
phiMin0->SetLineStyle(3);
phiMin0->SetMarkerStyle(20);
phiMin0->SetMarkerSize(0.4);
phiMax0->SetLineStyle(3);
phiMax0->SetMarkerStyle(20);
phiMax0->SetMarkerSize(0.4);
phiMin0->Draw("C");
phiMax0->Draw("C");
// LABEL ALLOWED REGION ///////////////////////////////////////////////
TPaveText *allowed = new TPaveText(lblxmin,lblymin,lblxmax,lblymax,"NDC");
TText *text = allowed->AddText("#font[42]{allowed (95% CL)}");
allowed->SetTextSize(0.04);
if (strpltmd.compare("tan") == 0) text->SetTextAngle(270);
allowed->SetFillStyle(0);
allowed->SetLineColor(0);
allowed->SetBorderSize(1);
allowed->Draw();
// SAVE GRAPHIC ///////////////////////////////////////////////////////
MyC->Print(epsfile.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");
}
