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 11 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 14 of file HcalPiecewiseLinearFunctor.cc.

References abscissae_, Exception, first, mps_fire::i, AbsHcalFunctor::isStrictlyIncreasing(), edm::second(), createJobs::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 40 of file HcalPiecewiseLinearFunctor.h.

{}

Member Function Documentation

HcalPiecewiseLinearFunctor HcalPiecewiseLinearFunctor::inverse

(

)

const

Definition at line 65 of file HcalPiecewiseLinearFunctor.cc.

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

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

References abscissae_, leftExtrapolationLinear_, 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 61 of file HcalPiecewiseLinearFunctor.cc.

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

Referenced by inverse().

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

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

overridevirtual

Implements AbsHcalFunctor.

Definition at line 35 of file HcalPiecewiseLinearFunctor.cc.

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

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

References abscissae_, leftExtrapolationLinear_, rightExtrapolationLinear_, and values_.

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

double HcalPiecewiseLinearFunctor::xmax ( ) const

overridevirtual

Reimplemented from AbsHcalFunctor.

Definition at line 89 of file HcalPiecewiseLinearFunctor.cc.

References abscissae_, and mps_fire::result.

Referenced by svgfig.XAxis::__repr__().

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

double HcalPiecewiseLinearFunctor::xmin ( ) const

overridevirtual

Reimplemented from AbsHcalFunctor.

Definition at line 82 of file HcalPiecewiseLinearFunctor.cc.

References abscissae_, and mps_fire::result.

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

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

Member Data Documentation

std::vector<double> HcalPiecewiseLinearFunctor::abscissae_

private

Definition at line 62 of file HcalPiecewiseLinearFunctor.h.

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

bool HcalPiecewiseLinearFunctor::leftExtrapolationLinear_

private

Definition at line 64 of file HcalPiecewiseLinearFunctor.h.

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

bool HcalPiecewiseLinearFunctor::rightExtrapolationLinear_

private

Definition at line 65 of file HcalPiecewiseLinearFunctor.h.

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

std::vector<double> HcalPiecewiseLinearFunctor::values_

private

Definition at line 63 of file HcalPiecewiseLinearFunctor.h.

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