scipy.stats.truncnorm uses standardized bounds: a = (lower - loc) / scale.
Using a=0 set the lower bound at loc (the mean distance), cutting off the
entire lower half of the Gaussian. Fix: a = -loc/scale so the lower bound
is at distance=0 (physical constraint), preserving the full prior shape.