Accurate evaluation of an integral involving the product of two Bessel functions and a Gaussian