8 #include "GaudiKernel/MsgStream.h"
23 const unsigned int order) {
29 std::vector<double> res_param(
order + 3);
32 std::vector<SamplePoint> my_points(sample_points);
37 const auto [tLow ,tHigh] =
interval(sample_points);
40 double mean(0.5 * (res_param[1] + res_param[0]));
41 double length(res_param[1] - res_param[0]);
42 for (
unsigned int k = 0;
k < my_points.size();
k++) { my_points[
k].set_x1(2 * (sample_points[
k].
x1() -
mean) /
length); }
47 fitter.fit_parameters(my_points, 1, sample_points.size(), chebyshev);
48 for (
unsigned int k = 0;
k <
order + 1;
k++) { res_param[
k + 2] =
fitter.coefficients()[
k]; }
66 unsigned int nb_points(100);
67 double bin_width((
res.tUpper() -
res.tLower()) /
static_cast<double>(nb_points));
68 std::vector<double> res_param(nb_points + 2);
73 res_param[0] =
res.tLower();
75 for (
unsigned int k = 0;
k < nb_points;
k++) { res_param[
k + 2] =
res.resolution(
res.tLower() +
k *
bin_width); }
78 return rt_res_relation_look_up;