|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
|
#ifndef CalibGBRMATH
#define CalibGBRMATH
#include "RVersion.h"
#if ROOT_VERSION_CODE >= ROOT_VERSION(5, 34, 00)
#include "vdt/vdtMath.h"
#else
#include "DataFormats/Math/interface/VDTMath.h"
#endif
#include <cmath>
namespace gbrmath {
inline double fast_pow(double base, double exponent) {
if (base == 0. && exponent > 0.)
return 0.;
else if (base > 0.)
return vdt::fast_exp(exponent * vdt::fast_log(base));
else
return std::nan("");
}
} // namespace gbrmath
#endif
|