Added smoothing of probability distributions to improve robustness of results

Added Cauchy smoothing to 2D factors and Gaussian smoothing to 1D ones
to reduce significant overfitting in the high-dimensional scenario that
we approached by setting the discetization resolution to 50 for each
variable. This greatly improves the robustness of results especially in
respect with the D-statistics null distribution calculation.