XMM-Newton Science Analysis System


ssclib (ssclib-4.33) [xmmsas_20190531_1155-18.0.0]

Integrated Gaussian probability

This function returns the probability that a gaussian-distributed variable $y$ will depart from the mean $y_{\rm {mean}}$ by greater than $abs(y-y_{\rm {mean}})$. For values obeying a gaussian distribution of standard deviation $\sigma$, the probability $P$ of a given $y$ value (or greater) occurring by chance is

$$P(y) = 1 - erf \big[ abs(y - y_{\rm {mean}})/\sigma\sqrt{2} \big].$$

  function integratedGaussProb(testY, meanY, sigma) result(probability)
    real(single), intent(in) :: testY, meanY, sigma
    real(single) :: probability
  end function integratedGaussProb