M3::Plotting Namespace Reference

Flexible, experiment-agnostic plotting utilities for MaCh3. More...

Classes
struct	InputFile
	Struct which wraps around the actual input file and also holds general information, and data from the file to be used by other classes. More...
class	InputManager
	This guy talks to the input files and is intended to smooth over the fact that our OA fitters use such a beautiful and diverse range of output file structures. More...
class	PlottingManager
	The main class to be used in plotting scripts. More...
class	StyleManager
	EW: provides centralized styling utilities for plots, including name prettification and style application. More...

Enumerations
enum	fileTypeEnum {   kLLH , kPostFit , kMCMC , kSigmaVar ,   kNFileTypes }
	Types of possible file that can be read. More...

Functions
TH1D	TGraphToTH1D (const TGraph &graph, const std::string &newName="", const std::string &newTitle="")
	This handy little function lets you interpret a TGraph as a TH1D. More...
std::vector< std::vector< double > >	TGraphToVector (const TGraph &graph)
	This handy little function lets you interpret a TGraph as a vector containing the same data. More...
std::vector< std::vector< double > >	TGraphToVector (const TGraph2D &graph)
	This handy little function lets you interpret a 2d TGraph as a vector containing the same data. More...
void	SetSymmetricRatioRange (const std::vector< std::unique_ptr< TH1D >> &RatioPlot)
	Set a symmetric Y-axis range around 1 for ratio plots. More...

Detailed Description

Flexible, experiment-agnostic plotting utilities for MaCh3.

Enumeration Type Documentation

fileTypeEnum

enum M3::Plotting::fileTypeEnum

Types of possible file that can be read.

Enumerator
kLLH	Log Likelihood scan.
kPostFit	Processed post fit errors.
kMCMC	MCMC chain.
kSigmaVar	Sigma variations.
kNFileTypes	Number of types of file.

Definition at line 16 of file InputManager.h.

                  {
  kLLH,         
  kPostFit,     
  kMCMC,        
  kSigmaVar,    
 
  kNFileTypes 
};

Function Documentation

SetSymmetricRatioRange()

void M3::Plotting::SetSymmetricRatioRange

(

const std::vector< std::unique_ptr< TH1D >> &

RatioPlot

)

Set a symmetric Y-axis range around 1 for ratio plots.

Definition at line 111 of file PlottingUtils.cpp.

{
  double maxz = -999;
  double minz = +999;
 
  for (int j = 0; j < static_cast<int>(RatioPlot.size()); j++) {
    if (!RatioPlot[j]) continue;
 
    for (int i = 1; i < RatioPlot[0]->GetXaxis()->GetNbins(); i++) {
      maxz = std::max(maxz, RatioPlot[j]->GetBinContent(i));
      minz = std::min(minz, RatioPlot[j]->GetBinContent(i));
    }
  }
 
  maxz = maxz * 1.001;
  minz = minz * 1.001;
 
  if (std::fabs(1 - maxz) > std::fabs(1 - minz))
    RatioPlot[0]->GetYaxis()->SetRangeUser(1 - std::fabs(1 - maxz), 1 + std::fabs(1 - maxz));
  else
    RatioPlot[0]->GetYaxis()->SetRangeUser(1 - std::fabs(1 - minz), 1 + std::fabs(1 - minz));
}

TGraphToTH1D()

TH1D M3::Plotting::TGraphToTH1D	(	const TGraph &	graph,
		const std::string &	newName,
		const std::string &	newTitle
	)

This handy little function lets you interpret a TGraph as a TH1D.

It will go through your provided graph and make a histogram binning by taking the midpoints of all the graph points x values, then fill the histogram with the graphs y values. For the first and last points it will extend the binning out of the graph bounds using the width between the outermost and second outermost points. This can be useful if e.g. you want to draw cumulative stacks of LLH scans.

Parameters

graph	The graph you want to convert.
newName	The new name you want to give to the histogram. If not specified, will just use the name of the graph.
newTitle	The new title you want to give to the histogram. If not specified, will just use the title of the graph.

Definition at line 10 of file PlottingUtils.cpp.

                                                                                            {
  std::string name;
  std::string title;
 
  if (newName == "")
    name = graph.GetName();
  else
    name = newName;
 
  if (newTitle == "")
    title = graph.GetTitle();
  else
    title = newTitle;
 
  int nPoints = graph.GetN();
  if(nPoints < 2){
    MACH3LOG_ERROR("Too few points in the graph.");
    throw MaCh3Exception(__FILE__,__LINE__);
  }
  std::vector<double> pointsX(nPoints);
  std::vector<double> pointsY(nPoints);
 
  // Get the points out
  Double_t x, y;
 
  for (int pointId = 0; pointId < nPoints; pointId++)
  {
    graph.GetPoint(pointId, x, y);
    pointsX[pointId] = x;
    pointsY[pointId] = y;
  }
 
  // get the bin edges
  std::vector<double> binEdges(nPoints + 1);
  binEdges[0] = pointsX[0] - (pointsX[1] - pointsX[0]) / 2.0;
  binEdges[nPoints] = pointsX[nPoints - 1] + (pointsX[nPoints - 1] - pointsX[nPoints - 2]) / 2.0;
 
  for (int pointId = 1; pointId < nPoints; pointId++)
  {
    // take the midpoint of the two surrounding points
    binEdges[pointId] = (pointsX[pointId] + pointsX[pointId - 1]) / 2.0;
  }
 
  TH1D retHist = TH1D(name.c_str(), title.c_str(), nPoints, binEdges.data());
 
  for (int binId = 0; binId < nPoints; binId++)
  {
    retHist.SetBinContent(binId + 1, pointsY[binId]);
  }
 
  return retHist;
}

TGraphToVector() [1/2]

std::vector< std::vector< double > > M3::Plotting::TGraphToVector

(

const TGraph &

graph

)

This handy little function lets you interpret a TGraph as a vector containing the same data.

Parameters

graph

The graph you want to convert.

Returns

A vector of vectors containing the data from the initial graph. The first vector is the x axis, the 2nd the y axis

Definition at line 63 of file PlottingUtils.cpp.

                                                                 {
  int nPoints = graph.GetN();
  std::vector<std::vector<double>> ret(2);
  std::vector<double> pointsX(nPoints);
  std::vector<double> pointsY(nPoints);
 
  // Get the points out
  Double_t x, y;
 
  for (int pointId = 0; pointId < nPoints; pointId++)
  {
    graph.GetPoint(pointId, x, y);
    pointsX[pointId] = x;
    pointsY[pointId] = y;
  }
 
  ret[0] = pointsX;
  ret[1] = pointsY;
 
  return ret;
}

TGraphToVector() [2/2]

std::vector< std::vector< double > > M3::Plotting::TGraphToVector

(

const TGraph2D &

graph

)

This handy little function lets you interpret a 2d TGraph as a vector containing the same data.

Parameters

graph

The graph you want to convert.

Returns

A vector of vectors containing the data from the initial graph. The first vector is the x axis, the 2nd the y axis, the 3rd is the z axis

Definition at line 86 of file PlottingUtils.cpp.

                                                                   {
  int nPoints = graph.GetN();
  std::vector<std::vector<double>> ret(3);
  std::vector<double> pointsX(nPoints);
  std::vector<double> pointsY(nPoints);
  std::vector<double> pointsZ(nPoints);
 
  // Get the points out
  Double_t x, y, z;
 
  for (int pointId = 0; pointId < nPoints; pointId++)
  {
    graph.GetPoint(pointId, x, y, z);
    pointsX[pointId] = x;
    pointsY[pointId] = y;
    pointsZ[pointId] = z;
  }
 
  ret[0] = pointsX;
  ret[1] = pointsY;
  ret[2] = pointsZ;
 
  return ret;
}