PiecewiseScalingPolynomial.h

Go to the documentation of this file.

#ifndef CondFormats_HcalObjects_PiecewiseScalingPolynomial_h_
#define CondFormats_HcalObjects_PiecewiseScalingPolynomial_h_

#include "FWCore/Utilities/interface/Exception.h"

#include "boost/serialization/vector.hpp"
#include "boost/serialization/version.hpp"

#include <vector>

class PiecewiseScalingPolynomial {
public:
  inline PiecewiseScalingPolynomial() {}

  PiecewiseScalingPolynomial(const std::vector<std::vector<double> >& coeffs, const std::vector<double>& limits);

  inline double operator()(const double x) const {
    double scale(0.0);
    if (x > 0.0) {
      const unsigned nLimits(limits_.size());
      if (nLimits) {
        const double* limits(&limits_[0]);
        unsigned which(0U);
        for (; which < nLimits; ++which)
          if (x < limits[which])
            break;
        const std::vector<double>& c(coeffs_[which]);
        const double* a = &c[0];
        for (int deg = c.size() - 1; deg >= 0; --deg) {
          scale *= x;
          scale += a[deg];
        }
      }
    }
    return scale * x;
  }

  inline bool operator==(const PiecewiseScalingPolynomial& r) const {
    return coeffs_ == r.coeffs_ && limits_ == r.limits_;
  }

  inline bool operator!=(const PiecewiseScalingPolynomial& r) const { return !(*this == r); }

private:
  bool validate() const;

  std::vector<std::vector<double> > coeffs_;
  std::vector<double> limits_;

  friend class boost::serialization::access;

  template <class Archive>
  inline void save(Archive& ar, const unsigned /* version */) const {
    if (!validate())
      throw cms::Exception("In PiecewiseScalingPolynomial::save: invalid data");
    ar & coeffs_ & limits_;
  }

  template <class Archive>
  inline void load(Archive& ar, const unsigned /* version */) {
    ar & coeffs_ & limits_;
    if (!validate())
      throw cms::Exception("In PiecewiseScalingPolynomial::load: invalid data");
  }

  BOOST_SERIALIZATION_SPLIT_MEMBER()
};

BOOST_CLASS_VERSION(PiecewiseScalingPolynomial, 1)

#endif  // CondFormats_HcalObjects_PiecewiseScalingPolynomial_h_