From 57e07da3790ccc6c75123745786cc58cbbf876a2 Mon Sep 17 00:00:00 2001 From: markus Date: Fri, 19 Feb 2021 19:01:37 +0100 Subject: [PATCH] add classroom_model --- cara/montecarlo.py | 29 ++++++++++++++++++++++++++++- 1 file changed, 28 insertions(+), 1 deletion(-) diff --git a/cara/montecarlo.py b/cara/montecarlo.py index f68c5f67..121b45df 100644 --- a/cara/montecarlo.py +++ b/cara/montecarlo.py @@ -963,7 +963,34 @@ exposure_models = [MCExposureModel( ) ) for qid in (100, 60)] -print(np.mean(large_population_baselines[0].infection_probability())) +classroom_model = MCExposureModel( + concentration_model=MCConcentrationModel( + room=models.Room(volume=160), + ventilation=models.SlidingWindow( + active=models.PeriodicInterval(period=120, duration=10), + inside_temp=models.PiecewiseConstant((0, 24), (293,)), + outside_temp=models.PiecewiseConstant((0, 24), (283,)), + window_height=1.6, opening_length=0.6, + ), + infected=MCInfectedPopulation( + number=1, + presence=models.SpecificInterval(((0, 2), (2.5, 4), (5, 7), (7.5, 9))), + masked=False, + virus=MCVirus(halflife=1.1, qID=60), + expiratory_activity=2, + samples=200000, + breathing_category=3, + ) + ), + exposed=models.Population( + number=19, + presence=models.SpecificInterval(((0, 2), (2.5, 4), (5, 7), (7.5, 9))), + activity=models.Activity.types['Seated'], + mask=models.Mask.types['No mask'] + ) +) + +plot_concentration_curve(classroom_model) #generate_cdf_curves_vs_qr(masked=False)