From ecc35ee03c7b6d07b2a3dcaf1ce617496ebe4588 Mon Sep 17 00:00:00 2001
From: markus <marro98@gmail.com>
Date: Mon, 18 Jan 2021 16:32:04 +0100
Subject: [PATCH] add unit conversions in qr-calculation

---
 cara/montecarlo.py | 5 +++++
 1 file changed, 5 insertions(+)

diff --git a/cara/montecarlo.py b/cara/montecarlo.py
index 1d0b85a7..2dfd4a69 100644
--- a/cara/montecarlo.py
+++ b/cara/montecarlo.py
@@ -15,6 +15,11 @@ def calculate_qr(viral_load: float, emission: float, diameter: float, mask_effic
     """
     Calculates the quantum generation rate given a set of parameters.
     """
+    # Unit conversions
+    diameter *= 1e-4
+    viral_load = 10 ** viral_load
+    emission = (emission * 60) if breathing_rate is None else (emission * 1e6)
+
     volume = 4 * np.pi * (diameter / 2) ** 3 / 3
     if breathing_rate is None:
         breathing_rate = 1