HcalCubicInterpolator Class Reference

#include <HcalCubicInterpolator.h>

Inheritance diagram for HcalCubicInterpolator:

Public Types
typedef std::tuple< double, double, double >	Triple

Public Member Functions
HcalCubicInterpolator	approximateInverse () const
	HcalCubicInterpolator ()
	HcalCubicInterpolator (const std::vector< Triple > &points)
double	operator() (double x) const override
double	xmax () const override
double	xmin () const override
	~HcalCubicInterpolator () override
Public Member Functions inherited from AbsHcalFunctor
bool	operator!= (const AbsHcalFunctor &r) const
bool	operator== (const AbsHcalFunctor &r) const
virtual	~AbsHcalFunctor ()

Protected Member Functions
bool	isEqual (const AbsHcalFunctor &other) const override

Private Member Functions
template<class Archive >
void	serialize (Archive &ar, unsigned)

Private Attributes
std::vector< double >	abscissae_
std::vector< double >	derivatives_
std::vector< double >	values_

Friends
class	boost::serialization::access

Additional Inherited Members
Static Protected Member Functions inherited from AbsHcalFunctor
template<class Iter >
static bool	isStrictlyDecreasing (Iter begin, Iter const end)
template<class Iter >
static bool	isStrictlyIncreasing (Iter begin, Iter const end)

Detailed Description

Definition at line 17 of file HcalCubicInterpolator.h.

Member Typedef Documentation

typedef std::tuple<double,double,double> HcalCubicInterpolator::Triple

Definition at line 21 of file HcalCubicInterpolator.h.

Constructor & Destructor Documentation

HcalCubicInterpolator::HcalCubicInterpolator

(

)

Definition at line 6 of file HcalCubicInterpolator.cc.

Referenced by approximateInverse().

{
}

HcalCubicInterpolator::HcalCubicInterpolator ( const std::vector< Triple > &  points )

explicit

Definition at line 10 of file HcalCubicInterpolator.cc.

References abscissae_, derivatives_, Exception, mps_fire::i, jetUpdater_cfi::sort, protons_cff::t, tmp, and values_.

{
    const std::size_t sz = points.size();
    if (sz)
    {
        std::vector<Triple> tmp(points);
        std::sort(tmp.begin(), tmp.end());
        abscissae_.reserve(sz);
        values_.reserve(sz);
        derivatives_.reserve(sz);
        for (std::size_t i=0; i<sz; ++i)
        {
            const Triple& t(tmp[i]);
            abscissae_.push_back(std::get<0>(t));
            values_.push_back(std::get<1>(t));
            derivatives_.push_back(std::get<2>(t));
        }
        const std::size_t szm1 = sz - 1;
        for (std::size_t i=0; i<szm1; ++i)
            if (abscissae_[i] == abscissae_[i+1])
                throw cms::Exception("In HcalCubicInterpolator constructor:"
                                            " abscissae must not coincide");
    }
}

HcalCubicInterpolator::~HcalCubicInterpolator ( )

inlineoverride

Definition at line 31 of file HcalCubicInterpolator.h.

References approximateInverse(), operator()(), x, xmax(), and xmin().

{}

Member Function Documentation

HcalCubicInterpolator HcalCubicInterpolator::approximateInverse

(

)

const

Definition at line 87 of file HcalCubicInterpolator.cc.

References abscissae_, derivatives_, Exception, HcalCubicInterpolator(), mps_fire::i, AbsHcalFunctor::isStrictlyDecreasing(), AbsHcalFunctor::isStrictlyIncreasing(), hiPixelPairStep_cff::points, and values_.

Referenced by ~HcalCubicInterpolator().

{
    const bool monotonous = isStrictlyIncreasing(values_.begin(), values_.end()) ||
                            isStrictlyDecreasing(values_.begin(), values_.end());
    if (!monotonous)
        throw cms::Exception("In HcalCubicInterpolator::inverse:"
                                    " can't invert non-monotonous functor");
    const std::size_t sz = abscissae_.size();
    std::vector<Triple> points;
    points.reserve(sz);
    for (std::size_t i=0; i<sz; ++i)
    {
        const double dydx = derivatives_[i];
        if (dydx == 0.0)
            throw cms::Exception("In HcalCubicInterpolator::inverse:"
                                        " can't invert functor with derivatives of 0");
        points.push_back(Triple(values_[i], abscissae_[i], 1.0/dydx));
    }
    return HcalCubicInterpolator(points);
}

bool HcalCubicInterpolator::isEqual ( const AbsHcalFunctor &  other ) const

inlineoverrideprotectedvirtual

Implements AbsHcalFunctor.

Definition at line 42 of file HcalCubicInterpolator.h.

References abscissae_, derivatives_, trackingPlots::other, alignCSCRings::r, and values_.

    {
        const HcalCubicInterpolator& r =
            static_cast<const HcalCubicInterpolator&>(other);
        return abscissae_ == r.abscissae_ &&
               values_ == r.values_ &&
               derivatives_ == r.derivatives_;
    }

double HcalCubicInterpolator::operator() ( double  x ) const

overridevirtual

Implements AbsHcalFunctor.

Definition at line 36 of file HcalCubicInterpolator.cc.

References abscissae_, derivatives_, PVValHelper::dx, mps_fire::result, protons_cff::t, pfDeepBoostedJetPreprocessParams_cfi::upper_bound, values_, and x.

Referenced by ~HcalCubicInterpolator().

{
    double result = 0.0;
    const std::size_t sz = abscissae_.size();
    if (sz)
    {
        if (sz > 1)
        {
            const std::size_t szm1 = sz - 1;
            if (x >= abscissae_[szm1])
                result = values_[szm1] + derivatives_[szm1]*(x - abscissae_[szm1]);
            else if (x <= abscissae_[0])
                result = values_[0] + derivatives_[0]*(x - abscissae_[0]);
            else
            {
                const std::size_t cell = std::upper_bound(abscissae_.begin(), abscissae_.end(), x) -
                                         abscissae_.begin() - 1;
                const std::size_t cellp1 = cell + 1;
                const double dx = abscissae_[cellp1] - abscissae_[cell];
                const double t = (x - abscissae_[cell])/dx;
                const double onemt = 1.0 - t;
                const double h00 = onemt*onemt*(1.0 + 2.0*t);
                const double h10 = onemt*onemt*t;
                const double h01 = t*t*(3.0 - 2.0*t);
                const double h11 = t*t*onemt;
                result = h00*values_[cell] + h10*dx*derivatives_[cell] +
                         h01*values_[cellp1] - h11*dx*derivatives_[cellp1];
            }
        }
        else
            result = values_[0] + derivatives_[0]*(x - abscissae_[0]);
    }
    return result;
}

template<class Archive >

void HcalCubicInterpolator::serialize ( Archive &  ar, unsigned    )

inlineprivate

Definition at line 59 of file HcalCubicInterpolator.h.

References derivatives_.

    {
        boost::serialization::base_object<AbsHcalFunctor>(*this);
        ar & abscissae_ & values_ & derivatives_;
    }

double HcalCubicInterpolator::xmax ( ) const

overridevirtual

Reimplemented from AbsHcalFunctor.

Definition at line 79 of file HcalCubicInterpolator.cc.

References abscissae_, and mps_fire::result.

Referenced by svgfig.XAxis::__repr__(), and ~HcalCubicInterpolator().

{
    double result = 0.0;
    if (!abscissae_.empty())
        result = abscissae_.back();
    return result;
}

double HcalCubicInterpolator::xmin ( ) const

overridevirtual

Reimplemented from AbsHcalFunctor.

Definition at line 71 of file HcalCubicInterpolator.cc.

References abscissae_, and mps_fire::result.

Referenced by svgfig.XAxis::__repr__(), svgfig.Axes::__repr__(), svgfig.HGrid::__repr__(), svgfig.Grid::__repr__(), svgfig.Axes::SVG(), and ~HcalCubicInterpolator().

{
    double result = 0.0;
    if (!abscissae_.empty())
        result = abscissae_[0];
    return result;
}

Friends And Related Function Documentation

friend class boost::serialization::access

friend

Definition at line 56 of file HcalCubicInterpolator.h.

Member Data Documentation

std::vector<double> HcalCubicInterpolator::abscissae_

private

Definition at line 52 of file HcalCubicInterpolator.h.

Referenced by approximateInverse(), HcalCubicInterpolator(), isEqual(), operator()(), xmax(), and xmin().

std::vector<double> HcalCubicInterpolator::derivatives_

private

Definition at line 54 of file HcalCubicInterpolator.h.

Referenced by approximateInverse(), HcalCubicInterpolator(), isEqual(), operator()(), and serialize().

std::vector<double> HcalCubicInterpolator::values_

private

Definition at line 53 of file HcalCubicInterpolator.h.

Referenced by approximateInverse(), HcalCubicInterpolator(), isEqual(), and operator()().