pingu_98/cara
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

add calculate_qr

Browse source
This commit is contained in:
markus 2021-01-18 14:57:34 +01:00
parent d01ebda19e
commit d096855bec
1 changed files with 11 additions and 0 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

11
cara/montecarlo.py
View file

@ -0,0 +1,11 @@
import cara.models
import numpy as np
def calculate_qr(viral_load: float, emission_concentration: float, diameter: float, mask_efficiency: float,
copies_per_quantum: float, breathing_rate: float = 1) -> float:
"""
Calculates the quantum generation rate given a set of parameters.
"""
volume = 4 * np.pi * (diameter / 2) ** 3 / 3
return viral_load * emission_concentration * volume * (1 - mask_efficiency) * breathing_rate / copies_per_quantum