ana::FitVarsProductMarg Class Reference

See the documentation for FitVarsProductMarg. More...

#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-01-16/CAFAna/Vars/FitVarsProduct.h"

Inheritance diagram for ana::FitVarsProductMarg:

Public Member Functions
	FitVarsProductMarg (const IConstrainedFitVar &x, const IConstrainedFitVar &y)
virtual double	GetValue (const osc::IOscCalcAdjustable *osc) const
virtual void	SetValue (osc::IOscCalcAdjustable *osc, double val) const
virtual double	LowLimit () const
virtual double	HighLimit () const
double	Penalty (double val, osc::IOscCalcAdjustable *) const override
	Override the default implementation with one that pushes back towards the center of the constrained space. More...
const std::string &	ShortName () const
const std::string &	LatexName () const

Protected Member Functions
template<typename T >
T	Clamp (T val) const

Protected Attributes
const IConstrainedFitVar &	fVarX
const IConstrainedFitVar &	fVarY

Detailed Description

See the documentation for FitVarsProductMarg.

Definition at line 42 of file FitVarsProduct.h.

Constructor & Destructor Documentation

ana::FitVarsProductMarg::FitVarsProductMarg	(	const IConstrainedFitVar &	x,
		const IConstrainedFitVar &	y
	)

Definition at line 135 of file FitVarsProduct.cxx.

References ana::assert(), ana::IConstrainedFitVar::HighLimit(), and ana::IConstrainedFitVar::LowLimit().

    : IConstrainedFitVar(x.ShortName()+"*"+y.ShortName()+"_marg",
                         "Marginalization parameter for "+x.LatexName()+y.LatexName()),
      fVarX(x), fVarY(y)
  {
    assert(x.LowLimit() == 0 && x.HighLimit() == 1);
    assert(y.LowLimit() == 0 && y.HighLimit() == 1);
  }

Member Function Documentation

template<typename T >

T ana::IConstrainedFitVar::Clamp ( T  val ) const

inlineprotectedinherited

Definition at line 60 of file IFitVar.h.

References std::max(), std::min(), ana::StanExp(), ana::StanLog(), and submit_syst::x.

Referenced by ana::FitEpsMuTau::SetValue(), ana::FitDmSq32CPT::SetValue(), ana::FitVarsProduct::SetValue(), ana::FitDmSq32ScaledCPT::SetValue(), SetValue(), ana::FitDmSq32CPTHierarchy::SetValue(), FitSinSqTheta23UO::SetValue(), FitSinSqTheta23LO::SetValue(), ana::FitDmSq21CPT::SetValue(), ana::FitSinSq2Theta23CPT::SetValue(), ana::FitSinSqTheta23CPT::SetValue(), ana::FitSinSq2Theta13CPT::SetValue(), ana::FitSinSq2Theta12CPT::SetValue(), FitSinSqTheta23LowerOctant::SetValue(), FitSinSqTheta23MaxMix::SetValue(), and FitSinSqTheta23UpperOctant::SetValue().

      {
        return std::max(util::GetValAs<T>(LowLimit()), std::min(val, util::GetValAs<T>(HighLimit())));
      }

double ana::FitVarsProductMarg::GetValue ( const osc::IOscCalcAdjustable *  osc ) const

virtual

Implements ana::IFitVar.

Definition at line 146 of file FitVarsProduct.cxx.

References beta, fVarX, fVarY, and ana::IFitVar::GetValue().

  {
    // Just evaluate the formula for beta.
    double alpha, beta;
    XYToAlphaBeta(fVarX.GetValue(osc), fVarY.GetValue(osc), alpha, beta);
    return beta;
  }

virtual double ana::FitVarsProductMarg::HighLimit ( ) const

inlinevirtual

Implements ana::IConstrainedFitVar.

Definition at line 52 of file FitVarsProduct.h.

{return 1;}

const std::string& ana::IFitVar::LatexName ( ) const

inlineinherited

Definition at line 37 of file IFitVar.h.

References ana::IFitVar::fLatexName, ana::IFitVar::fShortName, and string.

Referenced by ana::FrequentistSurface::CreateHistograms(), ana::SurfaceKrige::SurfaceKrige(), and test_fitvarsproduct().

{ return fLatexName; };

virtual double ana::FitVarsProductMarg::LowLimit ( ) const

inlinevirtual

Implements ana::IConstrainedFitVar.

Definition at line 51 of file FitVarsProduct.h.

{return 0;}

double ana::IConstrainedFitVar::Penalty ( double  val, osc::IOscCalcAdjustable *    ) const

overridevirtualinherited

Override the default implementation with one that pushes back towards the center of the constrained space.

Reimplemented from ana::IFitVar.

Definition at line 8 of file IFitVar.cxx.

References hi(), ana::IConstrainedFitVar::HighLimit(), lo(), ana::IConstrainedFitVar::LowLimit(), extractScale::mean, Munits::rad, and util::sqr().

Referenced by ana::SigmaDelta::SigmaPenalty().

  {
    const auto lo = LowLimit();
    const auto 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 auto mean = (lo + hi) / 2;
    const auto 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;
  }

void ana::FitVarsProductMarg::SetValue ( osc::IOscCalcAdjustable *  osc, double  val  ) const

virtual

Implements ana::IFitVar.

Definition at line 155 of file FitVarsProduct.cxx.

References beta, ana::IConstrainedFitVar::Clamp(), fVarX, fVarY, ana::IFitVar::GetValue(), ana::IFitVar::SetValue(), submit_syst::x, and submit_syst::y.

  {
    // Switch to the other space, hold alpha constant and adjust beta, see
    // where that leaves us.
    double alpha, beta;
    XYToAlphaBeta(fVarX.GetValue(osc), fVarY.GetValue(osc), alpha, beta);
    beta = Clamp(val);

    double x, y;
    AlphaBetaToXY(alpha, beta, x, y);
    fVarX.SetValue(osc, x);
    fVarY.SetValue(osc, y);
  }

const std::string& ana::IFitVar::ShortName ( ) const

inlineinherited

Definition at line 36 of file IFitVar.h.

References ana::IFitVar::fShortName.

Referenced by ana::CovMxSurface::CovMxSurface(), ana::FrequentistSurface::FillSurface(), MakeSurface(), MakeSurfaceJoint(), plot_3flavor_withsysts(), ana::FrequentistSurface::ProgressBarTitle(), ana::MCMCSamples::SampleValue(), ana::StanFitter::transform_helper(), ana::MCMCSample::Val(), ana::MCMCSamples::VarOffset(), and ana::StanFitter::write_array().

{ return fShortName; };

Member Data Documentation

const IConstrainedFitVar& ana::FitVarsProductMarg::fVarX

protected

Definition at line 54 of file FitVarsProduct.h.

Referenced by GetValue(), and SetValue().

const IConstrainedFitVar& ana::FitVarsProductMarg::fVarY

protected

Definition at line 55 of file FitVarsProduct.h.

Referenced by GetValue(), and SetValue().

---

The documentation for this class was generated from the following files:

• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-01-16/CAFAna/Vars/FitVarsProduct.h
• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-01-16/CAFAna/Vars/FitVarsProduct.cxx