3#if defined(R__HAS_VECCORE) && defined(R__HAS_VC)
36 out *= 0.3989422804014327 * inv_sigma;
66 return vecCore::Blend<::ROOT::Double_v>(
81 Double_t p10 = 2.4266795523053175e+2, q10 = 2.1505887586986120e+2, p11 = 2.1979261618294152e+1,
82 q11 = 9.1164905404514901e+1, p12 = 6.9963834886191355e+0, q12 = 1.5082797630407787e+1,
83 p13 = -3.5609843701815385e-2;
85 Double_t p20 = 3.00459261020161601e+2, q20 = 3.00459260956983293e+2, p21 = 4.51918953711872942e+2,
86 q21 = 7.90950925327898027e+2, p22 = 3.39320816734343687e+2, q22 = 9.31354094850609621e+2,
87 p23 = 1.52989285046940404e+2, q23 = 6.38980264465631167e+2, p24 = 4.31622272220567353e+1,
88 q24 = 2.77585444743987643e+2, p25 = 7.21175825088309366e+0, q25 = 7.70001529352294730e+1,
89 p26 = 5.64195517478973971e-1, q26 = 1.27827273196294235e+1, p27 = -1.36864857382716707e-7;
91 Double_t p30 = -2.99610707703542174e-3, q30 = 1.06209230528467918e-2, p31 = -4.94730910623250734e-2,
92 q31 = 1.91308926107829841e-1, p32 = -2.26956593539686930e-1, q32 = 1.05167510706793207e+0,
93 p33 = -2.78661308609647788e-1, q33 = 1.98733201817135256e+0, p34 = -2.23192459734184686e-2, q34 = 1;
101 vecCore::Mask<::ROOT::Double_v> mask3 = !(mask1 || mask2);
111 vecCore::MaskedAssign<::ROOT::Double_v>(
112 result, mask1,
v * (p10 + p11 * v2 + p12 * v4 + p13 * v6) / (q10 + q11 * v2 + q12 * v4 + v6));
113 vecCore::MaskedAssign<::ROOT::Double_v>(
116 (p20 + p21 *
v + p22 * v2 + p23 * v3 + p24 * v4 + p25 * v5 + p26 * v6 + p27 * v7) /
117 (
vecCore::math::Exp(v2) * (q20 + q21 *
v + q22 * v2 + q23 * v3 + q24 * v4 + q25 * v5 + q26 * v6 + v7)));
118 vecCore::MaskedAssign<::ROOT::Double_v>(result, mask3,
120 (
c1 + (p30 * v8 + p31 * v6 + p32 * v4 + p33 * v2 + p34) /
121 ((q30 * v8 + q31 * v6 + q32 * v4 + q33 * v2 + q34) * v2)) /
139 y * (-2.057706e-2 +
y * (2.635537e-2 +
y * (-1.647633e-2 +
y * 3.92377e-3))))))));
147 y * (3.5156229 +
y * (3.0899424 +
y * (1.2067492 +
y * (0.2659732 +
y * (3.60768e-2 +
y * 4.5813e-3)))));
154 return vecCore::Blend<::ROOT::Double_v>(ax <= 3.75, BesselI0_Split_Less(
x), BesselI0_Split_More(ax));
164 (0.39894228 +
y * (-3.988024e-2 +
166 y * (1.63801e-3 +
y * (-1.031555e-2 +
167 y * (2.282967e-2 +
y * (-2.895312e-2 +
168 y * (1.787654e-2 +
y * -4.20059e-3))))))));
169 return vecCore::Blend<::ROOT::Double_v>(
x < 0,
::ROOT::Double_v(-1.0) * result, result);
176 return x * (0.5 +
y * (0.87890594 +
177 y * (0.51498869 +
y * (0.15084934 +
y * (2.658733e-2 +
y * (3.01532e-3 +
y * 3.2411e-4))))));
184 return vecCore::Blend<::ROOT::Double_v>(ax <= 3.75, BesselI1_Split_Less(
x), BesselI1_Split_More(ax,
x));
195 1 +
y * (-0.1098628627e-2 +
y * (0.2734510407e-4 +
y * (-0.2073370639e-5 +
y * 0.2093887211e-6)));
197 -0.1562499995e-1 +
y * (0.1430488765e-3 +
y * (-0.6911147651e-5 +
y * (0.7621095161e-6 -
y * 0.934935152e-7)));
205 return (57568490574.0 +
206 y * (-13362590354.0 +
y * (651619640.7 +
y * (-11214424.18 +
y * (77392.33017 +
y * -184.9052456))))) /
207 (57568490411.0 +
y * (1029532985.0 +
y * (9494680.718 +
y * (59272.64853 +
y * (267.8532712 +
y)))));
213 return vecCore::Blend<::ROOT::Double_v>(ax < 8, BesselJ0_Split1_Less(
x), BesselJ0_Split1_More(ax));
224 1 +
y * (0.183105e-2 +
y * (-0.3516396496e-4 +
y * (0.2457520174e-5 +
y * -0.240337019e-6)));
226 0.04687499995 +
y * (-0.2002690873e-3 +
y * (0.8449199096e-5 +
y * (-0.88228987e-6 -
y * 0.105787412e-6)));
229 vecCore::MaskedAssign<::ROOT::Double_v>(result,
x < 0, -result);
238 y * (-7895059235.0 +
y * (242396853.1 +
y * (-2972611.439 +
y * (15704.48260 +
y * -30.16036606))))) /
239 (144725228442.0 +
y * (2300535178.0 +
y * (18583304.74 +
y * (99447.43394 +
y * (376.9991397 +
y)))));
245 return vecCore::Blend<::ROOT::Double_v>(ax < 8, BesselJ1_Split1_Less(
x), BesselJ1_Split1_More(ax,
x));
double beta(double x, double y)
Calculates the beta function.
Double_t Sqrt(Double_t x)
Double_t Gaus(Double_t x, Double_t mean=0, Double_t sigma=1, Bool_t norm=kFALSE)
Calculate a gaussian function with mean and sigma.
Double_t Log2(Double_t x)
Double_t BesselI1(Double_t x)
modified Bessel function K_0(x)
Double_t LaplaceDistI(Double_t x, Double_t alpha=0, Double_t beta=1)
Computes the distribution function of Laplace distribution at point x, with location parameter alpha ...
Double_t BreitWigner(Double_t x, Double_t mean=0, Double_t gamma=1)
Calculate a Breit Wigner function with mean and gamma.
Double_t LaplaceDist(Double_t x, Double_t alpha=0, Double_t beta=1)
Computes the probability density function of Laplace distribution at point x, with location parameter...
Double_t BesselJ0(Double_t x)
modified Bessel function K_1(x)
Double_t BesselJ1(Double_t x)
Bessel function J0(x) for any real x.
Double_t Freq(Double_t x)
Computation of the normal frequency function freq(x).
Double_t BesselI0(Double_t x)
integer order modified Bessel function K_n(x)