From beed5c6db0e1cef7f8faa03a459c447f279fd138 Mon Sep 17 00:00:00 2001
From: Lorenzo Moneta <Lorenzo.Moneta@cern.ch>
Date: Fri, 1 Mar 2019 14:32:07 +0100
Subject: [PATCH] apply fix in VavilovFast distribution as suggested by L.
Grzanks (ROOT-9899)
---
math/mathmore/src/VavilovFast.cxx | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/math/mathmore/src/VavilovFast.cxx b/math/mathmore/src/VavilovFast.cxx
index 45721cb8c70..1d8200b8b6a 100644
--- a/math/mathmore/src/VavilovFast.cxx
+++ b/math/mathmore/src/VavilovFast.cxx
@@ -382,7 +382,7 @@ double VavilovFast::Pdf (double x) const
double s = 1 + fHC[7]*h[9];
for (k=2; k<=6; k++)
s += fHC[k]*h[k+1];
- if (s>0) v = fHC[8]*std::exp(-0.5*xx*xx);
+ if (s>0) v = fHC[8]*s*std::exp(-0.5*xx*xx);
}
else if (fItype == 2) {
double xx = x*x;
@@ -502,7 +502,7 @@ double VavilovFast::Quantile (double z) const {
for (int k = 2; k <= 6; ++k) {
y += fHC[k]*h[k+1];
}
- if (y > 0) fu = fHC[8]*std::exp(-0.5*x*x);
+ if (y > 0) fu = fHC[8]*y*std::exp(-0.5*x*x);
}
else if (fItype == 2) {
double x = rlam*rlam;
--
GitLab