NParameters.cxx

#include <vector>
#include <TMath.h>
#include "NLogger.h"
 
#include "NParameters.h"
 
ClassImp(Ndmspc::NParameters);
 
namespace Ndmspc {
NParameters::NParameters()
    : TNamed("parameters", "Parameters"), fHisto(nullptr), fNames({})
{
}
 
NParameters::NParameters(std::vector<std::string> parNames, const char * name, const char * title)
    : TNamed(name, title), fNames(parNames)
{
 
  if (parNames.empty()) {
    NLogTrace("NParameters::NParameters: No parameter names provided, creating empty parameters histogram.");
    return;
  }
  fHisto = new TH1D("ParametersHisto", "ParametersHisto", parNames.size(), 0, parNames.size());
  // set parameter names as labels
  for (size_t i = 0; i < parNames.size(); i++) {
    fHisto->GetXaxis()->SetBinLabel(i + 1, parNames[i].c_str());
  }
 
  fHisto->Sumw2(kFALSE); // Disable sum of squares of weights for error calculation
}
NParameters::~NParameters()
{
  delete fHisto;
}
 
void NParameters::Print(Option_t * /*option*/) const
{
  fHisto->Print("all");
}
 
bool NParameters::SetParameter(int bin, Double_t value, Double_t error)
{
  if (bin < 1 || bin > fHisto->GetNbinsX()) {
    return false;
  }
  fHisto->SetBinContent(bin, value);
  fHisto->SetBinError(bin, error);
  return true;
}
 
bool NParameters::SetParameter(const char * parName, Double_t value, Double_t error)
{
  int bin = fHisto->GetXaxis()->FindBin(parName);
  if (bin < 1 || bin > fHisto->GetNbinsX()) {
    NLogError("NParameters::SetParameter: Parameter name '%s' not found !!!", parName);
    return false;
  }
  fHisto->SetBinContent(bin, value);
  fHisto->SetBinError(bin, error);
  return true;
}
 
Double_t NParameters::GetParameter(int bin) const
{
  if (bin < 1 || bin > fHisto->GetNbinsX()) {
    NLogError("NParameters::GetParameter: Parameter index '%d' out of range !!!", bin);
    return TMath::QuietNaN();
  }
  return fHisto->GetBinContent(bin);
}
 
Double_t NParameters::GetParameter(const char * parName) const
{
  int bin = fHisto->GetXaxis()->FindBin(parName);
  if (bin < 1 || bin > fHisto->GetNbinsX()) {
    NLogError("NParameters::GetParameter: Parameter name '%s' not found !!!", parName);
    return TMath::QuietNaN();
  }
  return fHisto->GetBinContent(bin);
}
Double_t NParameters::GetParameterError(int bin) const
{
  if (bin < 1 || bin > fHisto->GetNbinsX()) {
    NLogError("NParameters::GetParError: Parameter index '%d' out of range !!!", bin);
    return TMath::QuietNaN();
  }
  return fHisto->GetBinError(bin);
}
Double_t NParameters::GetParameterError(const char * parName) const
{
  int bin = fHisto->GetXaxis()->FindBin(parName);
  if (bin < 1 || bin > fHisto->GetNbinsX()) {
    NLogError("NParameters::GetParError: Parameter name '%s' not found !!!", parName);
    return TMath::QuietNaN();
  }
  return fHisto->GetBinError(bin);
}
 
} // namespace Ndmspc