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_