Sample & Hold: Estimate jitter delay of an ADC
I think the jitter to SNR formula is based on a gaussian jitter.

You use "rand" which yields uniformly distributed numbers. You should use randn() instead which yields numbers distributed according to a gaussian distribution.