NuDockServerBase.cpp File Reference

#include "NuDockServerBase.h"
#include "Samples/SampleHandlerBase.h"

Include dependency graph for NuDockServerBase.cpp:

Go to the source code of this file.

Functions
void	TH1toResponse (const TH1 *mc_hist, const std::string &sample_title, nlohmann::json &response)
void	TH2toResponse (const TH2 *mc_hist, const std::string &sample_title, nlohmann::json &response)

Function Documentation

TH1toResponse()

void TH1toResponse	(	const TH1 *	mc_hist,
		const std::string &	sample_title,
		nlohmann::json &	response
	)

Definition at line 149 of file NuDockServerBase.cpp.

                                                                                              {
// ***************************************************************************
  auto hist1d = dynamic_cast<const TH1D*>(mc_hist);
  if (!hist1d) {
    MACH3LOG_ERROR("{}: input histogram is not TH1D", __func__);
    throw MaCh3Exception(__FILE__,__LINE__);
  }
  const TAxis *ax = mc_hist->GetXaxis();
  std::string xtitle = ax->GetTitle();
  int nxbins = ax->GetNbins();
  std::vector<double> xbins(nxbins+1);
  std::vector<double> binvals(nxbins);
 
  xbins[0] = ax->GetBinLowEdge(1);
  for (int ix = 1; ix <= nxbins; ++ix) {
    xbins[ix] = ax->GetBinUpEdge(ix);
    binvals[ix-1] = mc_hist->GetBinContent(ix);
  }
  response["axis_titles"][sample_title] = {xtitle};
  response["bin_edges"][sample_title] = {xbins};
  response["bin_values"][sample_title] = binvals;
}

TH2toResponse()

void TH2toResponse	(	const TH2 *	mc_hist,
		const std::string &	sample_title,
		nlohmann::json &	response
	)

Definition at line 173 of file NuDockServerBase.cpp.

                                                                                              {
// ***************************************************************************
  auto hist2d = dynamic_cast<const TH2D*>(mc_hist);
  if (!hist2d) {
    MACH3LOG_ERROR("{}: input histogram is not TH2D", __func__);
    throw MaCh3Exception(__FILE__,__LINE__);
  }
  const TAxis *x_axis = mc_hist->GetXaxis();
  std::string xtitle = x_axis->GetTitle();
  int nxbins = x_axis->GetNbins();
  std::vector<double> xbins(nxbins+1);
 
  const TAxis *y_axis = mc_hist->GetYaxis();
  std::string ytitle = y_axis->GetTitle();
  int nybins = y_axis->GetNbins();
  std::vector<double> ybins(nybins+1);
 
  std::vector<std::vector<double>> binvals(nxbins, std::vector<double>(nybins));
 
  xbins[0] = x_axis->GetBinLowEdge(1);
  ybins[0] = y_axis->GetBinLowEdge(1);
  for (int ix = 1; ix <= nxbins; ++ix) {
    for (int iy = 1; iy <= nybins; ++iy) {
      xbins[ix] = x_axis->GetBinUpEdge(ix);
      ybins[iy] = y_axis->GetBinUpEdge(iy);
      binvals[ix-1][iy-1] = mc_hist->GetBinContent(ix, iy);
    }
  }
  response["axis_titles"][sample_title] = {xtitle, ytitle};
  response["bin_edges"][sample_title] = {xbins, ybins};
  response["bin_values"][sample_title] = binvals;
}