LEPBandPlot: The plot for the CL bands "a la LEP". More...

#include <LEPBandPlot.h>

Inheritance diagram for LEPBandPlot:

Public Member Functions
void	draw (const char *options="") override
	Draw on canvas. More...

void	dumpToFile (const char RootFileName, const char options) override
	All the objects are written to rootfile. More...

	LEPBandPlot (const char name, const char title, const int n_points, double x_vals, double sb_vals, double b_vals, double b_rms, double *exp_vals=nullptr)
	Constructor. More...

	LEPBandPlot (const char name, const char title, const int n_points, double x_vals, double sb_vals, double b_vals, double b_up_bars1, double b_down_bars1, double b_up_bars2, double b_down_bars2, double exp_vals=nullptr)
	Constructor. More...

void	print (const char *options="") override
	Print the relevant information. More...

void	setTitle (const char *title)
	Set the title of the plot. More...

void	setXaxisTitle (const char *title)
	Set the title of the x axis. More...

	~LEPBandPlot () override
	Destructor. More...

Public Member Functions inherited from StatisticalPlot
void	dumpToImage (const char *filename)
	Write an image on disk. More...

TCanvas *	getCanvas ()
	Get the canvas. More...

bool	is_verbose ()
	get the verbosity More...

void	setCanvas (TCanvas *new_canvas)
	Set the canvas. More...

void	setVerbosity (bool verbosity)
	Set the verbosity. More...

	StatisticalPlot (const char name, const char title, bool verbosity=true)
	Constructor. More...

	~StatisticalPlot () override
	Destructor. More...

Private Attributes
TGraph *	m_b_band_graph_1sigma
	The b band 1 sigma. More...

TGraph *	m_b_band_graph_2sigma
	The b band 2 sigma. More...

TGraph *	m_b_line_graph
	The b line. More...

TGraph *	m_data_line_graph
	The data line. More...

TLegend *	m_legend
	The legend. More...

TGraph *	m_sb_line_graph
	The sb line. More...

TLine *	m_zero_line
	The line at 0. More...

Detailed Description

LEPBandPlot: The plot for the CL bands "a la LEP".

Revision: 1.4

Date: 2009/05/15 09:55:43

Author: D. Piparo (danilo.piparo<at>cern.ch), G. Schott - Universitaet Karlsruhe

This class allows to produce plots like the ones of the Tevatron HWG and originally introduced by LEP HWG. It is thought as a "enhanced" TGraph. The input to give to obtain the plot are:

The number of points
The values of the parameter on the x axis
The mean/median of the -2lnQ distribution of the sig+bkg hypothesis
The mean/median of the -2lnQ distribution of the bkg only hypothesis
The sigma (error bar size) of the -2lnQ distribution of the bkg only hypothesis
The experimental values of -2lnQ

Definition at line 41 of file LEPBandPlot.h.

Constructor & Destructor Documentation

LEPBandPlot::LEPBandPlot	(	const char *	name,
		const char *	title,
		const int	n_points,
		double *	x_vals,
		double *	sb_vals,
		double *	b_vals,
		double *	b_rms,
		double *	exp_vals = `nullptr`
	)

Constructor.

Definition at line 22 of file LEPBandPlot.cc.

References mps_fire::i, m_b_band_graph_1sigma, m_b_band_graph_2sigma, m_b_line_graph, m_data_line_graph, m_legend, m_sb_line_graph, and m_zero_line.

     : StatisticalPlot(name, title, false) {
   // bkg hypothesis line
   m_b_line_graph = new TGraph(n_points, x_vals, b_vals);
   m_b_line_graph->SetLineWidth(2);
   m_b_line_graph->SetLineStyle(2);
   m_b_line_graph->SetFillColor(kWhite);
 
   // bkg hypothesis band 1 sigma
   m_b_band_graph_1sigma = new TGraphErrors(n_points, x_vals, b_vals, nullptr, b_rms);
   m_b_band_graph_1sigma->SetFillColor(kGreen);
   m_b_band_graph_1sigma->SetLineColor(kGreen);
   m_b_band_graph_1sigma->SetMarkerColor(kGreen);
 
   // Make the band 2 times wider:
   double* b_2rms = new double[n_points];
   for (int i = 0; i < n_points; ++i)
     b_2rms[i] = 2 * b_rms[i];
 
   // bkg hypothesis band 2 sigma
   m_b_band_graph_2sigma = new TGraphErrors(n_points, x_vals, b_vals, nullptr, b_2rms);
   m_b_band_graph_2sigma->SetFillColor(kYellow);
   m_b_band_graph_2sigma->SetFillColor(kYellow);
   m_b_band_graph_2sigma->SetLineColor(kYellow);
   m_b_band_graph_2sigma->SetMarkerColor(kYellow);
   m_b_band_graph_2sigma->GetYaxis()->SetTitle("-2lnQ");
 
   // sig+bkg hypothesis line
   m_sb_line_graph = new TGraph(n_points, x_vals, sb_vals);
   m_sb_line_graph->SetLineWidth(2);
   m_sb_line_graph->SetLineStyle(4);
   m_sb_line_graph->SetLineColor(kRed);
   m_sb_line_graph->SetFillColor(kWhite);
 
   // The points of the data
   if (exp_vals != nullptr) {
     m_data_line_graph = new TGraph(n_points, x_vals, exp_vals);
     m_data_line_graph->SetLineWidth(2);
     m_data_line_graph->SetFillColor(kWhite);
   } else
     m_data_line_graph = nullptr;
 
   // Line for 0
   m_zero_line = new TLine(m_b_line_graph->GetXaxis()->GetXmin(), 0, m_b_line_graph->GetXaxis()->GetXmax(), 0);
 
   // The legend
 
   m_legend = new TLegend(0.75, 0.78, 0.98, 0.98);
   m_legend->SetName("Confidence Levels");
   m_legend->AddEntry(m_b_band_graph_1sigma, "-2lnQ #pm 1#sigma");
   m_legend->AddEntry(m_b_band_graph_2sigma, "-2lnQ #pm 2#sigma");
   m_legend->AddEntry(m_b_line_graph, "-2lnQ_{B}");
   m_legend->AddEntry(m_sb_line_graph, "-2lnQ_{SB}");
   if (m_data_line_graph != nullptr)
     m_legend->AddEntry(m_data_line_graph, "-2lnQ_{Obs}");
 
   m_legend->SetFillColor(0);
 
   delete[] b_2rms;
 }

LEPBandPlot::LEPBandPlot	(	const char *	name,
		const char *	title,
		const int	n_points,
		double *	x_vals,
		double *	sb_vals,
		double *	b_vals,
		double *	b_up_bars1,
		double *	b_down_bars1,
		double *	b_up_bars2,
		double *	b_down_bars2,
		double *	exp_vals = `nullptr`
	)

Constructor.

Definition at line 92 of file LEPBandPlot.cc.

References m_b_band_graph_1sigma, m_b_band_graph_2sigma, m_b_line_graph, m_data_line_graph, m_legend, m_sb_line_graph, and m_zero_line.

     : StatisticalPlot(name, title, false) {
   // bkg hypothesis line
   m_b_line_graph = new TGraph(n_points, x_vals, b_vals);
   m_b_line_graph->SetLineWidth(2);
   m_b_line_graph->SetLineStyle(2);
   m_b_line_graph->SetFillColor(kWhite);
 
   // bkg hypothesis band 1 sigma
   m_b_band_graph_1sigma = new TGraphAsymmErrors(n_points, x_vals, b_vals, nullptr, nullptr, b_down_bars1, b_up_bars1);
   m_b_band_graph_1sigma->SetFillColor(kGreen);
   m_b_band_graph_1sigma->SetLineColor(kGreen);
   m_b_band_graph_1sigma->SetMarkerColor(kGreen);
 
   // bkg hypothesis band 2 sigma
   m_b_band_graph_2sigma = new TGraphAsymmErrors(n_points, x_vals, b_vals, nullptr, nullptr, b_down_bars2, b_up_bars2);
   m_b_band_graph_2sigma->SetFillColor(kYellow);
   m_b_band_graph_2sigma->SetFillColor(kYellow);
   m_b_band_graph_2sigma->SetLineColor(kYellow);
   m_b_band_graph_2sigma->SetMarkerColor(kYellow);
   m_b_band_graph_2sigma->GetYaxis()->SetTitle("-2lnQ");
 
   // sig+bkg hypothesis line
   m_sb_line_graph = new TGraph(n_points, x_vals, sb_vals);
   m_sb_line_graph->SetLineWidth(2);
   m_sb_line_graph->SetLineStyle(4);
   m_sb_line_graph->SetLineColor(kRed);
   m_sb_line_graph->SetFillColor(kWhite);
 
   // The points of the data
   if (exp_vals != nullptr) {
     m_data_line_graph = new TGraph(n_points, x_vals, exp_vals);
     m_data_line_graph->SetLineWidth(2);
     m_data_line_graph->SetFillColor(kWhite);
   } else
     m_data_line_graph = nullptr;
 
   // Line for 0
   m_zero_line = new TLine(m_b_line_graph->GetXaxis()->GetXmin(), 0, m_b_line_graph->GetXaxis()->GetXmax(), 0);
 
   // The legend
 
   m_legend = new TLegend(0.75, 0.78, 0.98, 0.98);
   m_legend->SetName("Confidence Levels");
   m_legend->AddEntry(m_b_band_graph_1sigma, "-2lnQ #pm 1#sigma");
   m_legend->AddEntry(m_b_band_graph_2sigma, "-2lnQ #pm 2#sigma");
   m_legend->AddEntry(m_b_line_graph, "-2lnQ_{B}");
   m_legend->AddEntry(m_sb_line_graph, "-2lnQ_{SB}");
   if (m_data_line_graph != nullptr)
     m_legend->AddEntry(m_data_line_graph, "-2lnQ_{Obs}");
 
   m_legend->SetFillColor(0);
 }

LEPBandPlot::~LEPBandPlot ( )

override

Destructor.

Definition at line 158 of file LEPBandPlot.cc.

References m_b_band_graph_1sigma, m_b_band_graph_2sigma, m_b_line_graph, m_data_line_graph, m_legend, m_sb_line_graph, and m_zero_line.

                           {
   delete m_b_line_graph;
   delete m_sb_line_graph;
 
   delete m_b_band_graph_1sigma;
   delete m_b_band_graph_2sigma;
 
   if (m_data_line_graph != nullptr)
     delete m_data_line_graph;
 
   delete m_zero_line;
 
   delete m_legend;
 }

Member Function Documentation

void LEPBandPlot::draw ( const char * options = "" )

overridevirtual

Draw on canvas.

Implements StatisticalPlot.

Definition at line 191 of file LEPBandPlot.cc.

References StatisticalPlot::getCanvas(), m_b_band_graph_1sigma, m_b_band_graph_2sigma, m_b_line_graph, m_data_line_graph, m_legend, m_sb_line_graph, m_zero_line, runTheMatrix::opt, and StatisticalPlot::setCanvas().

                                           {
   setCanvas(new TCanvas(GetName(), GetTitle()));
   getCanvas()->cd();
 
   getCanvas()->SetGridx();
   getCanvas()->SetGridy();
 
   TString opt(options);
   // Bands
   if (opt.Contains("4") == 0) {
     m_b_band_graph_2sigma->Draw("A3");
     m_b_band_graph_1sigma->Draw("3");
   } else {
     m_b_band_graph_2sigma->Draw("A4");
     m_b_band_graph_1sigma->Draw("4");
   }
 
   // Lines
   if (opt.Contains("4") == 0) {
     m_b_line_graph->Draw("L");
     m_sb_line_graph->Draw("L");
   } else {
     m_b_line_graph->Draw("C");
     m_sb_line_graph->Draw("C");
   }
 
   if (m_data_line_graph != nullptr)
     m_data_line_graph->Draw("L");
 
   m_zero_line->Draw("Same");
 
   // Legend
   m_legend->Draw("Same");
 }

void LEPBandPlot::dumpToFile	(	const char *	RootFileName,
		const char *	options
	)

overridevirtual

All the objects are written to rootfile.

Implements StatisticalPlot.

Definition at line 228 of file LEPBandPlot.cc.

References m_b_band_graph_1sigma, m_b_band_graph_2sigma, m_b_line_graph, m_data_line_graph, m_legend, m_sb_line_graph, and indexGen::ofile.

                                                                           {
   TFile ofile(RootFileName, options);
   ofile.cd();
 
   // Bands
   m_b_band_graph_2sigma->Write("bkg_band_2sigma");
   m_b_band_graph_1sigma->Draw("bkg_band_1sigma");
 
   // Lines
   m_b_line_graph->Draw("bkg_line");
   m_sb_line_graph->Draw("sigbkg_line");
   if (m_data_line_graph)
     m_data_line_graph->Draw("observed_line");
 
   // Legend
   m_legend->Draw("IamTheLegend");
 
   ofile.Close();
 }