HcalPiecewiseLinearFunctor Class Reference

#include <HcalPiecewiseLinearFunctor.h>

Inheritance diagram for HcalPiecewiseLinearFunctor:

Public Member Functions
	HcalPiecewiseLinearFunctor ()
	HcalPiecewiseLinearFunctor (const std::vector< std::pair< double, double > > &points, bool leftExtrapolationLinear, bool rightExtrapolationLinear)
HcalPiecewiseLinearFunctor	inverse () const
bool	isStrictlyMonotonous () const
double	operator() (double x) const override
double	xmax () const override
double	xmin () const override
	~HcalPiecewiseLinearFunctor () 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_
bool	leftExtrapolationLinear_
bool	rightExtrapolationLinear_
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 HcalPiecewiseLinearFunctor.h.

Constructor & Destructor Documentation

HcalPiecewiseLinearFunctor::HcalPiecewiseLinearFunctor

(

)

Definition at line 14 of file HcalPiecewiseLinearFunctor.cc.

Referenced by inverse().

    : leftExtrapolationLinear_(false), rightExtrapolationLinear_(false)
{
}

HcalPiecewiseLinearFunctor::HcalPiecewiseLinearFunctor	(	const std::vector< std::pair< double, double > > &	points,
		bool	leftExtrapolationLinear,
		bool	rightExtrapolationLinear
	)

Definition at line 19 of file HcalPiecewiseLinearFunctor.cc.

References abscissae_, Exception, plotBeamSpotDB::first, mps_fire::i, AbsHcalFunctor::isStrictlyIncreasing(), hiPixelPairStep_cff::points, edm::second(), jetUpdater_cfi::sort, tmp, and values_.

    : leftExtrapolationLinear_(leftExtrapolationLinear),
      rightExtrapolationLinear_(rightExtrapolationLinear)
{
    const std::size_t sz = points.size();
    if (sz)
    {
        std::vector<std::pair<double, double> > tmp(points);
        std::sort(tmp.begin(), tmp.end());
        abscissae_.reserve(sz);
        values_.reserve(sz);
        for (std::size_t i=0; i<sz; ++i)
        {
            abscissae_.push_back(tmp[i].first);
            values_.push_back(tmp[i].second);
        }
        if (!isStrictlyIncreasing(abscissae_.begin(), abscissae_.end()))
            throw cms::Exception("In HcalPiecewiseLinearFunctor constructor:"
                                        " abscissae must not coincide");
    }
}

HcalPiecewiseLinearFunctor::~HcalPiecewiseLinearFunctor ( )

inlineoverride

Definition at line 41 of file HcalPiecewiseLinearFunctor.h.

References inverse(), isStrictlyMonotonous(), operator()(), x, xmax(), and xmin().

{}

Member Function Documentation

HcalPiecewiseLinearFunctor HcalPiecewiseLinearFunctor::inverse

(

)

const

Definition at line 89 of file HcalPiecewiseLinearFunctor.cc.

References abscissae_, Exception, HcalPiecewiseLinearFunctor(), mps_fire::i, isStrictlyMonotonous(), checklumidiff::l, leftExtrapolationLinear_, hiPixelPairStep_cff::points, alignCSCRings::r, rightExtrapolationLinear_, std::swap(), and values_.

Referenced by ~HcalPiecewiseLinearFunctor().

{
    if (!isStrictlyMonotonous())
        throw cms::Exception("In HcalPiecewiseLinearFunctor::inverse:"
                                    " can't invert non-monotonous functor");
    const std::size_t sz = abscissae_.size();
    std::vector<std::pair<double, double> > points;
    points.reserve(sz);
    for (std::size_t i=0; i<sz; ++i)
        points.push_back(std::make_pair(values_[i], abscissae_[i]));
    bool l = leftExtrapolationLinear_;
    bool r = rightExtrapolationLinear_;
    if (values_[0] > values_[sz - 1])
        std::swap(l, r);
    return HcalPiecewiseLinearFunctor(points, l, r);
}

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

inlineoverrideprotectedvirtual

Implements AbsHcalFunctor.

Definition at line 55 of file HcalPiecewiseLinearFunctor.h.

References abscissae_, leftExtrapolationLinear_, trackingPlots::other, alignCSCRings::r, rightExtrapolationLinear_, and values_.

    {
        const HcalPiecewiseLinearFunctor& r =
            static_cast<const HcalPiecewiseLinearFunctor&>(other);
        return abscissae_ == r.abscissae_ &&
               values_ == r.values_ &&
               leftExtrapolationLinear_ == r.leftExtrapolationLinear_ &&
               rightExtrapolationLinear_ == r.rightExtrapolationLinear_;
    }

bool HcalPiecewiseLinearFunctor::isStrictlyMonotonous

(

)

const

Definition at line 83 of file HcalPiecewiseLinearFunctor.cc.

References AbsHcalFunctor::isStrictlyDecreasing(), AbsHcalFunctor::isStrictlyIncreasing(), and values_.

Referenced by inverse(), and ~HcalPiecewiseLinearFunctor().

{
    return isStrictlyIncreasing(values_.begin(), values_.end()) ||
           isStrictlyDecreasing(values_.begin(), values_.end());
}

double HcalPiecewiseLinearFunctor::operator() ( double  x ) const

overridevirtual

Implements AbsHcalFunctor.

Definition at line 44 of file HcalPiecewiseLinearFunctor.cc.

References abscissae_, interpolateSimple(), leftExtrapolationLinear_, mps_fire::result, rightExtrapolationLinear_, pfDeepBoostedJetPreprocessParams_cfi::upper_bound, values_, and x.

Referenced by ~HcalPiecewiseLinearFunctor().

{
    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])
            {
                if (rightExtrapolationLinear_)
                    result = interpolateSimple(abscissae_[sz-2], abscissae_[szm1],
                                               values_[sz-2], values_[szm1], x);
                else
                    result = values_[szm1];
            }
            else if (x <= abscissae_[0])
            {
                if (leftExtrapolationLinear_)
                    result = interpolateSimple(abscissae_[0], abscissae_[1],
                                               values_[0], values_[1], x);
                else
                    result = values_[0];
            }
            else
            {
                const std::size_t cell = std::upper_bound(abscissae_.begin(), abscissae_.end(), x) -
                                         abscissae_.begin() - 1;
                result = interpolateSimple(abscissae_[cell], abscissae_[cell+1],
                                           values_[cell], values_[cell+1], x);
            }
        }
        else
            result = values_[0];
    }
    return result;
}

template<class Archive >

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

inlineprivate

Definition at line 74 of file HcalPiecewiseLinearFunctor.h.

References rightExtrapolationLinear_.

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

double HcalPiecewiseLinearFunctor::xmax ( ) const

overridevirtual

Reimplemented from AbsHcalFunctor.

Definition at line 114 of file HcalPiecewiseLinearFunctor.cc.

References abscissae_, and mps_fire::result.

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

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

double HcalPiecewiseLinearFunctor::xmin ( ) const

overridevirtual

Reimplemented from AbsHcalFunctor.

Definition at line 106 of file HcalPiecewiseLinearFunctor.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 ~HcalPiecewiseLinearFunctor().

{
    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 71 of file HcalPiecewiseLinearFunctor.h.

Member Data Documentation

std::vector<double> HcalPiecewiseLinearFunctor::abscissae_

private

Definition at line 66 of file HcalPiecewiseLinearFunctor.h.

Referenced by HcalPiecewiseLinearFunctor(), inverse(), isEqual(), operator()(), xmax(), and xmin().

bool HcalPiecewiseLinearFunctor::leftExtrapolationLinear_

private

Definition at line 68 of file HcalPiecewiseLinearFunctor.h.

Referenced by inverse(), isEqual(), and operator()().

bool HcalPiecewiseLinearFunctor::rightExtrapolationLinear_

private

Definition at line 69 of file HcalPiecewiseLinearFunctor.h.

Referenced by inverse(), isEqual(), operator()(), and serialize().

std::vector<double> HcalPiecewiseLinearFunctor::values_

private

Definition at line 67 of file HcalPiecewiseLinearFunctor.h.

Referenced by HcalPiecewiseLinearFunctor(), inverse(), isEqual(), isStrictlyMonotonous(), and operator()().