|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
|
#ifndef CondFormats_HcalObjects_HcalCubicInterpolator_h
#define CondFormats_HcalObjects_HcalCubicInterpolator_h
#include <vector>
#include <tuple>
#include "CondFormats/HcalObjects/interface/AbsHcalFunctor.h"
#include "boost/serialization/access.hpp"
#include "boost/serialization/version.hpp"
#include "boost/serialization/vector.hpp"
//
// Cubic Hermite spline interpolator in 1-d. See, for example,
// http://en.wikipedia.org/wiki/Cubic_Hermite_spline
//
class HcalCubicInterpolator : public AbsHcalFunctor {
public:
// Order: abscissa, desired value, derivative
typedef std::tuple<double, double, double> Triple;
// Dummy constructor, to be used for deserialization only
HcalCubicInterpolator();
// Normal constructor from the set of interpolated points.
// The points will be sorted internally, so they can be
// given in arbitrary order.
explicit HcalCubicInterpolator(const std::vector<Triple>& points);
inline ~HcalCubicInterpolator() override {}
double operator()(double x) const override;
double xmin() const override;
double xmax() const override;
// Cubic approximation to the inverse curve (note, not the real
// solution of the direct cubic equation). Use at your own risc.
HcalCubicInterpolator approximateInverse() const;
protected:
inline bool isEqual(const AbsHcalFunctor& other) const override {
const HcalCubicInterpolator& r = static_cast<const HcalCubicInterpolator&>(other);
return abscissae_ == r.abscissae_ && values_ == r.values_ && derivatives_ == r.derivatives_;
}
private:
std::vector<double> abscissae_;
std::vector<double> values_;
std::vector<double> derivatives_;
friend class boost::serialization::access;
template <class Archive>
inline void serialize(Archive& ar, unsigned /* version */) {
boost::serialization::base_object<AbsHcalFunctor>(*this);
ar & abscissae_ & values_ & derivatives_;
}
};
BOOST_CLASS_VERSION(HcalCubicInterpolator, 1)
BOOST_CLASS_EXPORT_KEY(HcalCubicInterpolator)
#endif // CondFormats_HcalObjects_HcalCubicInterpolator_h
|