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
virtual double	operator() (double x) const override
virtual double	xmax () const override
virtual double	xmin () const override
virtual	~HcalPiecewiseLinearFunctor ()
Public Member Functions inherited from AbsHcalFunctor
bool	operator!= (const AbsHcalFunctor &r) const
bool	operator== (const AbsHcalFunctor &r) const
virtual	~AbsHcalFunctor ()

Protected Member Functions
virtual 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(), 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");
    }
}

virtual HcalPiecewiseLinearFunctor::~HcalPiecewiseLinearFunctor ( )

inlinevirtual

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);
}

virtual 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_, 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()().