From 6ff3b90cebff4cb6e91ec9088b837acbb269b7d8 Mon Sep 17 00:00:00 2001
From: Andrejh <andre.henriques@cern.ch>
Date: Mon, 22 Mar 2021 17:41:51 +0100
Subject: [PATCH] print values for paper

---
 cara/model_scenarios_publication.py | 36 +++++++++++++++++++++++++----
 1 file changed, 32 insertions(+), 4 deletions(-)

diff --git a/cara/model_scenarios_publication.py b/cara/model_scenarios_publication.py
index e78f2ac2..649b02c1 100644
--- a/cara/model_scenarios_publication.py
+++ b/cara/model_scenarios_publication.py
@@ -2,7 +2,7 @@ from cara import models
 from cara.montecarlo import *
 
 ######### Standard exposure models ###########
-exposure_models = [MCExposureModel(
+qR_models = [MCExposureModel(
     concentration_model=MCConcentrationModel(
         room=models.Room(volume=45),
         ventilation=models.SlidingWindow(
@@ -25,7 +25,35 @@ exposure_models = [MCExposureModel(
     exposed=models.Population(
         number=2,
         presence=models.SpecificInterval(((0, 4), (5, 9))),
-        activity=models.Activity.types['Seated'],
-        mask=models.Mask.types['Type I']
+        activity=models.Activity.types['Light activity'],
+        mask=models.Mask.types['No mask']
     )
-) for a in (1, 2, 3)]
\ No newline at end of file
+) for a in (1, 2, 3)]
+
+qR_models_shout_light = [MCExposureModel(
+    concentration_model=MCConcentrationModel(
+        room=models.Room(volume=45),
+        ventilation=models.SlidingWindow(
+            active=models.PeriodicInterval(period=120, duration=120),
+            inside_temp=models.PiecewiseConstant((0, 24), (293,)),
+            outside_temp=models.PiecewiseConstant((0, 24), (283,)),
+            window_height=1.6, opening_length=0.2,
+        ),
+        infected=MCInfectedPopulation(
+            number=1,
+            presence=models.SpecificInterval(((0, 4), (5, 9))),
+            masked=False,
+            virus=MCVirus(halflife=1.1, qID=qid),
+            expiratory_activity=3,
+            samples=2000000,
+            breathing_category=3,
+            expiratory_activity_weights=(0.7, 0.3, 0)
+        )
+    ),
+    exposed=models.Population(
+        number=2,
+        presence=models.SpecificInterval(((0, 4), (5, 9))),
+        activity=models.Activity.types['Light activity'],
+        mask=models.Mask.types['No mask']
+    )
+) for qid in (100, 500, 1000)]
\ No newline at end of file