diff --git a/cara/monte_carlo/data.py b/cara/monte_carlo/data.py
index b3724524..759e8c06 100644
--- a/cara/monte_carlo/data.py
+++ b/cara/monte_carlo/data.py
@@ -38,28 +38,30 @@ symptomatic_vl_frequencies = LogCustomKernel(
     kernel_bandwidth=0.1
 )
 
+# From https://doi.org/10.1093/cid/ciaa1579
+infectious_virus_distribution = Uniform(0.15, 0.45)
 
 # From CERN-OPEN-2021-04 and refererences therein
 virus_distributions = {
     'SARS_CoV_2': mc.SARSCoV2(
                 viral_load_in_sputum=symptomatic_vl_frequencies,
                 infectious_dose=100,
-                viable_to_RNA=Uniform(0.15, 0.45),
+                viable_to_RNA=infectious_virus_distribution,
                 ),
     'SARS_CoV_2_B117': mc.SARSCoV2(
                 viral_load_in_sputum=symptomatic_vl_frequencies,
                 infectious_dose=60,
-                viable_to_RNA=Uniform(0.15, 0.45),
+                viable_to_RNA=infectious_virus_distribution,
                 ),
     'SARS_CoV_2_P1': mc.SARSCoV2(
                 viral_load_in_sputum=symptomatic_vl_frequencies,
                 infectious_dose=100/2.25,
-                viable_to_RNA=Uniform(0.15, 0.45),
+                viable_to_RNA=infectious_virus_distribution,
                 ),
     'SARS_CoV_2_B16172': mc.SARSCoV2(
                 viral_load_in_sputum=symptomatic_vl_frequencies,
                 infectious_dose=60/1.6,
-                viable_to_RNA=Uniform(0.15, 0.45),
+                viable_to_RNA=infectious_virus_distribution,
                 ),
 }
 
@@ -73,4 +75,4 @@ mask_distributions = {
     'FFP2': mc.Mask(Uniform(0.83, 0.91)),
 }
 
-infectious_virus_distribution = Uniform(0.15, 0.45)
+