diff --git a/cara/monte_carlo/data.py b/cara/monte_carlo/data.py
index b5e15884..9a6753b8 100644
--- a/cara/monte_carlo/data.py
+++ b/cara/monte_carlo/data.py
@@ -203,4 +203,4 @@ short_range_expiration_distributions = {
 
 
 # Derived from Fig 8 a) "stand-stand" in https://www.mdpi.com/1660-4601/17/4/1445/htm
-short_range_distances = Beta(alpha=1.4342766632654418, beta=27.49916410927064)
\ No newline at end of file
+short_range_distances = Beta(alpha=0.588715, beta=1.50214)
\ No newline at end of file
diff --git a/cara/monte_carlo/sampleable.py b/cara/monte_carlo/sampleable.py
index 5a29de81..2aeaf441 100644
--- a/cara/monte_carlo/sampleable.py
+++ b/cara/monte_carlo/sampleable.py
@@ -144,7 +144,7 @@ class Beta(SampleableDistribution):
         self.beta = beta
 
     def generate_samples(self, size: int) -> float_array_size_n:
-        return beta.rvs(a = self.alpha, b = self.beta, loc=0.5, scale=7, size=size)
+        return beta.rvs(a = self.alpha, b = self.beta, loc=0.5, scale=(2 - 0.5), size=size)
 
 
 _VectorisedFloatOrSampleable = typing.Union[