IFitVar.cxx

Go to the documentation of this file.

#include "sbnana/CAFAna/Core/IFitVar.h"

#include "sbnana/CAFAna/Core/MathUtil.h"

#include <cmath>

namespace ana
{
  //----------------------------------------------------------------------
  double IConstrainedFitVar::Penalty(double val,
                                     osc::IOscCalcAdjustable*) const
  {
    const double lo = LowLimit();
    const double hi = HighLimit();

    if(val >= lo && val <= hi) return 0;

    // Try to direct fit back towards centre of the space. Engineer penalty to
    // be zero at the limits.
    const double mean = (lo+hi)/2;
    const double rad = (hi-lo)/2;
    return util::sqr((val-mean)/rad)-1;


    //    if(val < lo) return util::sqr(lo-val);
    //    if(val > hi) return util::sqr(val-hi);
    //    return 0;
  }

  //----------------------------------------------------------------------
  double IConstrainedFitVar::Clamp(double val) const
  {
    return std::max(LowLimit(), std::min(val, HighLimit()));
  }
} // namespace