ProfileResult.cc

Go to the documentation of this file.

#include <atm/ProfileResult.h>
#include <utl/TabulatedFunctionErrors.h>
#include <cmath>

using namespace utl;
using namespace std;


namespace atm {

  ProfileResult::ProfileResult(const utl::TabulatedFunctionErrors& profile,
                               const Binning bX, const Binning bY) :
    fProfile(profile),
    fBinningX(bX),
    fBinningY(bY)
  { }


  double
  ProfileResult::Y(const double x)
    const
  {
    const double xx = (fBinningX == eLog) ? log(x) : x;
    return (fBinningY == eLog) ?  exp(fProfile.Y(xx)) : fProfile.Y(xx);
  }


  double
  ProfileResult::MaxX()
    const
  {
    const unsigned int n = fProfile.GetNPoints();
    const double value = n ? fProfile[n-1].X() : fProfile[0].X();
    return (fBinningX == eLog) ? exp(value) : value;
  }


  double
  ProfileResult::MinX()
    const
  {
    return (fBinningX == eLog) ? exp(fProfile[0].X()) : fProfile[0].X();
  }


  bool
  ProfileResult::IsInValidRange(const double x)
    const
  {
    return fProfile.IsInValidRange(fBinningX == eLog ? log(x) : x);
  }


  vector<double>
  ProfileResult::GetXValues()
    const
  {
    if (fBinningX == eLinear)
      return vector<double>(fProfile.XBegin(), fProfile.XEnd());
    else {
      vector<double> xv;
      xv.reserve(fProfile.GetNPoints());
      for (const auto& x : fProfile.XRange())
        xv.push_back(exp(x));
      return xv;
    }
  }


  vector<double>
  ProfileResult::GetYValues()
    const
  {
    if (fBinningY == eLinear)
      return vector<double>(fProfile.YBegin(), fProfile.YEnd());
    else {
      vector<double> yv;
      yv.reserve(fProfile.GetNPoints());
      for (const auto& y : fProfile.YRange())
        yv.push_back(exp(y));
      return yv;
    }
  }

}