Constructor.
Constructor Set Vavilov parameters kappa and beta2 and define whether to calculate density fcn or distribution fcn
44 const double xp[9] = {9.29, 2.47, .89, .36, .15, .07, .03, .02, 0.0};
45 const double xq[7] = {.012, .03, .08, .26, .87, 3.83, 11.0};
47 double q, u,
x,
c1, c2, c3, c4, d1, h4, h5, h6,
q2,
x1,
d, ll, ul, xf1, xf2, rv;
61 h_[4] = 1. - beta2 * 0.42278433999999998 + 7.6 /
kappa;
64 h4 = -7.6 /
kappa - (beta2 * .57721566 + 1);
71 for (lp = 0; lp < 9; ++lp) {
76 for (lq = 0; lq < 7; ++lq) {
81 auto f2 = [h_](
double x) {
89 h_[0] =
kappa * (beta2 * .57721566 + 2.) + 9.9166128600000008;
93 h_[1] = beta2 *
kappa;
95 h_[3] = omega_ * 1.5707963250000001;
96 auto f1 = [h_](
double x) {
return h_[0] + h_[1] *
std::log(h_[2] * x) - h_[3] *
x; };
99 d =
exp(kappa * (beta2 * (.57721566 - h5) + 1.)) * .31830988654751274;
102 a_[n - 1] = omega_ * .31830988654751274;
106 for (k = 1; k <
n; ++
k) {
114 xf1 = kappa * (beta2 * c1 - c4) - x * c2;
115 xf2 = x * c1 + kappa * (c3 + beta2 * c2) +
t0_ * x;
117 d1 = q * d * omega_ *
exp(xf1);
118 a_[l - 1] = d1 *
cos(xf2);
119 b_[l - 1] = -d1 *
sin(xf2);
121 d1 = q * d *
exp(xf1) /
k;
122 a_[l - 1] = d1 *
sin(xf2);
123 b_[l - 1] = d1 *
cos(xf2);
124 a_[n - 1] += q2 *
a_[l - 1];
Sin< T >::type sin(const T &t)
int dzero(double a, double b, double &x0, double &rv, double eps, int mxf, F func)
Private version of the exponential integral.
Cos< T >::type cos(const T &t)
double f1(double x, double const *h_)
Private version of the exponential integral.
Abs< T >::type abs(const T &t)
double f2(double x, double const *h_)
void sincosint(double x, double &sint, double &cint)
const int mode_
returns the limits on the non-zero (mode=0) or normalized region (mode=1)
static const G4double kappa
double expint(double x)
Private version of the sine integral.