diff --git a/cara/mc-output_IGH_paper.py b/cara/mc-output_IGH_paper.py index 1c798f91..79befe96 100644 --- a/cara/mc-output_IGH_paper.py +++ b/cara/mc-output_IGH_paper.py @@ -1,179 +1,223 @@ from cara.montecarlo import * from cara.model_scenarios_IGH_paper import * -##Fig 1a -#compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0],classroom_model_IGH_no_mask_windows_open_breaks[0], +# #Fig 1c +# compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0],classroom_model_IGH_no_mask_windows_open_breaks[0], # classroom_model_IGH_no_mask_windows_open_alltimes[0]], # labels=['Windows closed', 'Window open during breaks', 'Window open at all times'], # colors=['tomato', 'lightskyblue', 'limegreen', '#1f77b4', 'seagreen', 'lightskyblue', 'deepskyblue'], -# title='No mask - Spring/Summer period' +# title='' # ) -#print_qd_info(classroom_model_IGH_no_mask_windows_closed[0]) -#print_qd_info(classroom_model_IGH_no_mask_windows_open_breaks[0]) -#print_qd_info(classroom_model_IGH_no_mask_windows_open_alltimes[0]) +# print("\n############ Fig 1c ###########\n") +# print_qd_info(classroom_model_IGH_no_mask_windows_closed[0]) +# print_qd_info(classroom_model_IGH_no_mask_windows_open_breaks[0]) +# print_qd_info(classroom_model_IGH_no_mask_windows_open_alltimes[0]) # -##Fig 1b -#compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0],classroom_model_IGH_no_mask_2windows_open_breaks[0], -# classroom_model_IGH_no_mask_2windows_open_alltimes[0]], -# labels=['Windows closed', '2 windows open during breaks', '2 windows open at all times'], -# colors=['tomato', 'lightskyblue', 'limegreen', '#1f77b4', 'seagreen', 'lightskyblue', 'deepskyblue'], -# title='No mask - Spring/Summer period' -# ) -#print_qd_info(classroom_model_IGH_no_mask_windows_closed[0]) -#print_qd_info(classroom_model_IGH_no_mask_2windows_open_breaks[0]) -#print_qd_info(classroom_model_IGH_no_mask_2windows_open_alltimes[0]) - -##Fig 1c -#compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0],classroom_model_IGH_no_mask_windows_open_breaks_winter[0], +# #Fig 1d +# compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed_winter[0],classroom_model_IGH_no_mask_windows_open_breaks_winter[0], # classroom_model_IGH_no_mask_windows_open_alltimes_winter[0], classroom_model_IGH_no_mask_windows_fully_open_alltimes_winter[0]], # labels=['Windows closed', 'Window fully open during breaks', 'Window slightly open at all times', 'Window fully open at all times'], # colors=['tomato', 'lightskyblue', '#1f77b4', 'limegreen', 'seagreen', 'lightskyblue', 'deepskyblue'], -# title='No mask - Winter period' +# title='' # ) -#print_qd_info(classroom_model_IGH_no_mask_windows_closed[0]) -#print_qd_info(classroom_model_IGH_no_mask_windows_open_breaks_winter[0]) -#print_qd_info(classroom_model_IGH_no_mask_windows_open_alltimes_winter[0]) -#print_qd_info(classroom_model_IGH_no_mask_windows_fully_open_alltimes_winter[0]) +# print("\n############ Fig 1d ###########\n") +# print_qd_info(classroom_model_IGH_no_mask_windows_closed_winter[0]) +# print_qd_info(classroom_model_IGH_no_mask_windows_open_breaks_winter[0]) +# print_qd_info(classroom_model_IGH_no_mask_windows_open_alltimes_winter[0]) +# print_qd_info(classroom_model_IGH_no_mask_windows_fully_open_alltimes_winter[0]) # - -##Fig 1d -#compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0],classroom_model_IGH_no_mask_2windows_open_breaks_winter[0], +# #Fig 1e +# compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0],classroom_model_IGH_no_mask_2windows_open_breaks[0], +# classroom_model_IGH_no_mask_2windows_open_alltimes[0]], +# labels=['Windows closed', '2 windows open during breaks', '2 windows open at all times'], +# colors=['tomato', 'lightskyblue', 'limegreen', '#1f77b4', 'seagreen', 'lightskyblue', 'deepskyblue'], +# title='' +# ) +# print("\n############ Fig 1e ###########\n") +# print_qd_info(classroom_model_IGH_no_mask_windows_closed[0]) +# print_qd_info(classroom_model_IGH_no_mask_2windows_open_breaks[0]) +# print_qd_info(classroom_model_IGH_no_mask_2windows_open_alltimes[0]) +# +# #Fig 1f +# compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed_winter[0],classroom_model_IGH_no_mask_2windows_open_breaks_winter[0], # classroom_model_IGH_no_mask_2windows_open_alltimes_winter[0], classroom_model_IGH_no_mask_2windows_fully_open_alltimes_winter[0]], # labels=['Windows closed', '2 windows fully open during breaks', '2 windows slighty open at all times', '2 windows fully open at all times'], # colors=['tomato', 'lightskyblue', '#1f77b4', 'limegreen', 'seagreen', 'lightskyblue', 'deepskyblue'], -# title='No mask - Winter period' +# title='' # ) -#print_qd_info(classroom_model_IGH_no_mask_windows_closed[0]) -#print_qd_info(classroom_model_IGH_no_mask_2windows_open_breaks_winter[0]) -#print_qd_info(classroom_model_IGH_no_mask_2windows_open_alltimes_winter[0]) -#print_qd_info(classroom_model_IGH_no_mask_2windows_fully_open_alltimes_winter[0]) +# print("\n############ Fig 1f ###########\n") +# print_qd_info(classroom_model_IGH_no_mask_windows_closed_winter[0]) +# print_qd_info(classroom_model_IGH_no_mask_2windows_open_breaks_winter[0]) +# print_qd_info(classroom_model_IGH_no_mask_2windows_open_alltimes_winter[0]) +# print_qd_info(classroom_model_IGH_no_mask_2windows_fully_open_alltimes_winter[0]) # -##Fig 1e -#compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0],classroom_model_IGH_no_mask_6windows_open_breaks[0], +# #Fig 1g +# compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0],classroom_model_IGH_no_mask_6windows_open_breaks[0], # classroom_model_IGH_no_mask_6windows_open_alltimes[0]], # labels=['Windows closed', '6 windows open during breaks', '6 windows open at all times'], # colors=['tomato', 'lightskyblue', 'limegreen', '#1f77b4', 'seagreen', 'lightskyblue', 'deepskyblue'], -# title='No mask - Spring/Summer period' +# title='' # ) -#print_qd_info(classroom_model_IGH_no_mask_windows_closed[0]) -#print_qd_info(classroom_model_IGH_no_mask_6windows_open_breaks[0]) -#print_qd_info(classroom_model_IGH_no_mask_6windows_open_alltimes[0]) +# print("\n############ Fig 1g ###########\n") +# print_qd_info(classroom_model_IGH_no_mask_windows_closed[0]) +# print_qd_info(classroom_model_IGH_no_mask_6windows_open_breaks[0]) +# print_qd_info(classroom_model_IGH_no_mask_6windows_open_alltimes[0]) # -##Fig 1f -#compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0],classroom_model_IGH_no_mask_6windows_open_breaks_winter[0], +# #Fig 1h +# compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed_winter[0],classroom_model_IGH_no_mask_6windows_open_breaks_winter[0], # classroom_model_IGH_no_mask_6windows_open_alltimes_winter[0], classroom_model_IGH_no_mask_6windows_fully_open_alltimes_winter[0]], # labels=['Windows closed', '6 windows fully open during breaks', '6 windows slighty open at all times', '6 windows fully open at all times'], # colors=['tomato', 'lightskyblue', '#1f77b4', 'limegreen', 'seagreen', 'lightskyblue', 'deepskyblue'], -# title='No mask - Winter period' +# title='' # ) -#print_qd_info(classroom_model_IGH_no_mask_windows_closed[0]) -#print_qd_info(classroom_model_IGH_no_mask_6windows_open_breaks_winter[0]) -#print_qd_info(classroom_model_IGH_no_mask_6windows_open_alltimes_winter[0]) -#rint_qd_info(classroom_model_IGH_no_mask_6windows_fully_open_alltimes_winter[0]) +# print("\n############ Fig 1h ###########\n") +# print_qd_info(classroom_model_IGH_no_mask_windows_closed_winter[0]) +# print_qd_info(classroom_model_IGH_no_mask_6windows_open_breaks_winter[0]) +# print_qd_info(classroom_model_IGH_no_mask_6windows_open_alltimes_winter[0]) +# print_qd_info(classroom_model_IGH_no_mask_6windows_fully_open_alltimes_winter[0]) # -##Fig 2 -#compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0],classroom_model_IGH_no_mask_windows_closed_1HEPA[0], +# #Fig 1i +# compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0], +# classroom_model_IGH_no_mask_6windows_open_breaks_endOfClass[0]], +# labels=['Windows closed', +# '6 windows open during breaks + end of classes'], +# colors=['tomato', '#1f77b4'], +# title='' +# ) +# print("\n############ Fig 1i ###########\n") +# print_qd_info(classroom_model_IGH_no_mask_windows_closed[0]) +# print_qd_info(classroom_model_IGH_no_mask_6windows_open_breaks_endOfClass[0]) +# +# #Fig 1j +# compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed_winter[0], +# classroom_model_IGH_no_mask_6windows_open_breaks_endOfClass_winter[0]], +# labels=['Windows closed', +# '6 windows open during breaks + end of classes'], +# colors=['tomato', '#1f77b4'], +# title='' +# ) +# print("\n############ Fig 1j ###########\n") +# print_qd_info(classroom_model_IGH_no_mask_windows_closed_winter[0]) +# print_qd_info(classroom_model_IGH_no_mask_6windows_open_breaks_endOfClass_winter[0]) +# +# +# +# #Fig 2a +# compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0],classroom_model_IGH_no_mask_windows_closed_1HEPA[0], # classroom_model_IGH_no_mask_windows_closed_2HEPA[0]], # labels=['Windows closed', 'Windows closed + 1 HEPA filter (2.5 ACH)', 'Windows closed + 2 HEPA filter (5 ACH)'], # colors=['tomato', 'lightskyblue', 'limegreen', '#1f77b4', 'seagreen', 'lightskyblue', 'deepskyblue'], -# title='No mask' +# title='' # ) -#print_qd_info(classroom_model_IGH_no_mask_windows_closed[0]) -#print_qd_info(classroom_model_IGH_no_mask_windows_closed_1HEPA[0]) -#print_qd_info(classroom_model_IGH_no_mask_windows_closed_2HEPA[0]) +# print("\n############ Fig 2a ###########\n") +# print_qd_info(classroom_model_IGH_no_mask_windows_closed[0]) +# print_qd_info(classroom_model_IGH_no_mask_windows_closed_1HEPA[0]) +# print_qd_info(classroom_model_IGH_no_mask_windows_closed_2HEPA[0]) # -##Fig 3 -#compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0],classroom_model_IGH_with_mask_windows_closed[0]], +# #Fig 2b +# compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed_winter[0],classroom_model_IGH_no_mask_windows_closed_1HEPA_winter[0], +# classroom_model_IGH_no_mask_windows_closed_2HEPA_winter[0]], +# labels=['Windows closed', 'Windows closed + 1 HEPA filter (2.5 ACH)', 'Windows closed + 2 HEPA filter (5 ACH)'], +# colors=['tomato', 'lightskyblue', 'limegreen', '#1f77b4', 'seagreen', 'lightskyblue', 'deepskyblue'], +# title='' +# ) +# print("\n############ Fig 2b ###########\n") +# print_qd_info(classroom_model_IGH_no_mask_windows_closed_winter[0]) +# print_qd_info(classroom_model_IGH_no_mask_windows_closed_1HEPA_winter[0]) +# print_qd_info(classroom_model_IGH_no_mask_windows_closed_2HEPA_winter[0]) +# +# #Fig 2c +# compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0],classroom_model_IGH_with_mask_windows_closed[0]], # labels=['Windows closed', 'Windows closed + surgical type masks'], # colors=['tomato', 'lightskyblue', 'limegreen', '#1f77b4', 'seagreen', 'lightskyblue', 'deepskyblue'], -# title='Use of mask' +# title='' # ) -#print_qd_info(classroom_model_IGH_no_mask_windows_closed[0]) -#print_qd_info(classroom_model_IGH_with_mask_windows_closed[0]) +# print("\n############ Fig 2c ###########\n") +# print_qd_info(classroom_model_IGH_no_mask_windows_closed[0]) +# print_qd_info(classroom_model_IGH_with_mask_windows_closed[0]) # -##Fig 4a -#compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0],classroom_model_IGH_no_mask_2windows_open_breaks_winter_bis[0], +# #Fig 2d +# compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed_winter[0],classroom_model_IGH_with_mask_windows_closed_winter[0]], +# labels=['Windows closed', 'Windows closed + surgical type masks'], +# colors=['tomato', 'lightskyblue', 'limegreen', '#1f77b4', 'seagreen', 'lightskyblue', 'deepskyblue'], +# title='' +# ) +# print("\n############ Fig 2d ###########\n") +# print_qd_info(classroom_model_IGH_no_mask_windows_closed_winter[0]) +# print_qd_info(classroom_model_IGH_with_mask_windows_closed_winter[0]) +# +# +# #Fig 3b +# compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed_winter[0],classroom_model_IGH_no_mask_2windows_open_alltimes_winter_bis[0], +# classroom_model_IGH_with_mask_2windows_open_alltimes_winter[0]], +# labels=['Windows closed', '2 windows open at all times', '2 windows open at all times + surgical masks'], +# colors=['tomato', 'lightskyblue', 'limegreen', '#1f77b4', 'seagreen', 'lightskyblue', 'deepskyblue'], +# title='' +# ) +# print("\n############ Fig 3b ###########\n") +# print_qd_info(classroom_model_IGH_no_mask_windows_closed_winter[0]) +# print_qd_info(classroom_model_IGH_no_mask_2windows_open_alltimes_winter_bis[0]) +# print_qd_info(classroom_model_IGH_with_mask_2windows_open_alltimes_winter[0]) +# +# #Fig 3c +# compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed_winter[0],classroom_model_IGH_with_mask_windows_closed_1HEPA[0], +# classroom_model_IGH_with_mask_windows_closed_2HEPA[0]], +# labels=['No masks', '1 HEPA filter (2.5 ACH) + surgical masks', '2 HEPA filter (5 ACH) + surgical masks'], +# colors=['tomato', 'lightskyblue', 'limegreen', '#1f77b4', 'seagreen', 'lightskyblue', 'deepskyblue'], +# title='' +# ) +# print("\n############ Fig 3c ###########\n") +# print_qd_info(classroom_model_IGH_no_mask_windows_closed[0]) +# print_qd_info(classroom_model_IGH_with_mask_windows_closed_1HEPA[0]) +# print_qd_info(classroom_model_IGH_with_mask_windows_closed_2HEPA[0]) +# +# #Fig 3d +# compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed_winter[0],classroom_model_IGH_with_mask_2windows_open_alltimes_winter_1HEPA[0], +# classroom_model_IGH_with_mask_2windows_open_alltimes_winter_2HEPA[0]], +# labels=['Windows closed', '2 windows open at all times + 1 HEPA filter (2.5 ACH) + surgical masks', '2 windows open at all times + 2 HEPA filter (5 ACH) + surgical masks'], +# colors=['tomato', 'lightskyblue', 'limegreen', '#1f77b4', 'seagreen', 'lightskyblue', 'deepskyblue'], +# title='' +# ) +# print("\n############ Fig 3d ###########\n") +# print_qd_info(classroom_model_IGH_no_mask_windows_closed_winter[0]) +# print_qd_info(classroom_model_IGH_with_mask_2windows_open_alltimes_winter_1HEPA[0]) +# print_qd_info(classroom_model_IGH_with_mask_2windows_open_alltimes_winter_2HEPA[0]) +# +# +# #Hourly breaks +# print("\n############ Hourly breaks - summer ###########\n") +# print_qd_info(classroom_model_IGH_no_mask_windows_open_breaks_endOfClass[0]) +# print_qd_info(classroom_model_IGH_no_mask_2windows_open_breaks_endOfClass[0]) +# print_qd_info(classroom_model_IGH_no_mask_6windows_open_breaks_endOfClass[0]) +# +# print("\n############ Hourly breaks - winter ###########\n") +# print_qd_info(classroom_model_IGH_no_mask_windows_open_breaks_endOfClass_winter[0]) +# print_qd_info(classroom_model_IGH_no_mask_2windows_open_breaks_endOfClass_winter[0]) +# print_qd_info(classroom_model_IGH_no_mask_6windows_open_breaks_endOfClass_winter[0]) + + + +#>>>>>>>>>>> Aux >>>>>>>> + + +#print_qd_info(classroom_model_IGH_no_mask_windows_closed[0]) +#print_qd_info(classroom_model_IGH_no_mask_6windows_open_breaks_endOfClass[0]) +#print_qd_info(classroom_model_IGH_no_mask_2windows_open_breaks_endOfClass_winter[0]) +#print_qd_info(classroom_model_IGH_no_mask_windows_open_breaks_endOfClass_winter[0]) + +# compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0],classroom_model_IGH_no_mask_2windows_open_breaks_winter_bis[0], # classroom_model_IGH_with_mask_2windows_open_breaks_winter[0]], # labels=['Windows closed', '2 windows open during breaks', '2 windows open during breaks + surgical masks'], # colors=['tomato', 'lightskyblue', 'limegreen', '#1f77b4', 'seagreen', 'lightskyblue', 'deepskyblue'], # title='Combination of measures (winter period) ' # ) -# -##Fig 4b -#compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0],classroom_model_IGH_with_mask_2windows_open_breaks_winter_1HEPA[0], + +# compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0],classroom_model_IGH_with_mask_2windows_open_breaks_winter_1HEPA[0], # classroom_model_IGH_with_mask_2windows_open_breaks_winter_2HEPA[0]], # labels=['Windows closed', '2 windows open during breaks + 1 HEPA filter (2.5 ACH) + surgical masks', '2 windows open during breaks + 2 HEPA filter (5 ACH) + surgical masks'], # colors=['tomato', 'lightskyblue', 'limegreen', '#1f77b4', 'seagreen', 'lightskyblue', 'deepskyblue'], # title='Combination of measures (winter period) ' # ) -# -##Fig 4c -#compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0],classroom_model_IGH_no_mask_2windows_open_alltimes_winter_bis[0], -# classroom_model_IGH_with_mask_2windows_open_alltimes_winter[0]], -# labels=['Windows closed', '2 windows open at all times', '2 windows open at all times + surgical masks'], -# colors=['tomato', 'lightskyblue', 'limegreen', '#1f77b4', 'seagreen', 'lightskyblue', 'deepskyblue'], -# title='Combination of measures (winter period) ' -# ) -#print_qd_info(classroom_model_IGH_no_mask_windows_closed[0]) -#print_qd_info(classroom_model_IGH_no_mask_2windows_open_alltimes_winter_bis[0]) -#print_qd_info(classroom_model_IGH_with_mask_2windows_open_alltimes_winter[0]) -# -##Fig 4d -#compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0],classroom_model_IGH_with_mask_2windows_open_alltimes_winter_1HEPA[0], -# classroom_model_IGH_with_mask_2windows_open_alltimes_winter_2HEPA[0]], -# labels=['Windows closed', '2 windows open at all times + 1 HEPA filter (2.5 ACH) + surgical masks', '2 windows open at all times + 2 HEPA filter (5 ACH) + surgical masks'], -# colors=['tomato', 'lightskyblue', 'limegreen', '#1f77b4', 'seagreen', 'lightskyblue', 'deepskyblue'], -# title='Combination of measures (winter period) ' -# ) -#print_qd_info(classroom_model_IGH_no_mask_windows_closed[0]) -#print_qd_info(classroom_model_IGH_with_mask_2windows_open_alltimes_winter_1HEPA[0]) -#print_qd_info(classroom_model_IGH_with_mask_2windows_open_alltimes_winter_2HEPA[0]) -# - -#compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0], -# classroom_model_IGH_no_mask_6windows_open_breaks_endOfClass[0]], -# labels=['Windows closed', -# '6 windows open during breaks + end of classes'], -# colors=['tomato', '#1f77b4'], -# title='Without mask - Spring/Summer period' -# ) -# -#print_qd_info(classroom_model_IGH_no_mask_windows_closed[0]) -##print_qd_info(classroom_model_IGH_no_mask_6windows_open_alltimes[0]) -#print_qd_info(classroom_model_IGH_no_mask_6windows_open_breaks_endOfClass[0]) - -#compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0], -# classroom_model_IGH_no_mask_6windows_open_breaks_endOfClass_winter[0]], -# labels=['Windows closed', -# '6 windows open during breaks + end of classes'], -# colors=['tomato', '#1f77b4'], -# title='Without mask - Winter period' -# ) -# -#print_qd_info(classroom_model_IGH_no_mask_windows_closed[0]) -#print_qd_info(classroom_model_IGH_no_mask_6windows_open_alltimes_winter[0]) -#print_qd_info(classroom_model_IGH_no_mask_6windows_open_breaks_endOfClass_winter[0]) - -#print_qd_info(classroom_model_IGH_no_mask_2windows_open_breaks_endOfClass_winter[0]) -#print_qd_info(classroom_model_IGH_no_mask_windows_open_breaks_endOfClass_winter[0]) -#print_qd_info(classroom_model_IGH_no_mask_windows_closed[0]) - -#print_qd_info(classroom_model_IGH_no_mask_windows_closed[0]) -#print_qd_info(classroom_model_IGH_no_mask_6windows_open_breaks_endOfClass[0]) -#print_qd_info(classroom_model_IGH_no_mask_2windows_open_breaks_endOfClass_winter[0]) -#print_qd_info(classroom_model_IGH_no_mask_windows_open_breaks_endOfClass_winter[0]) - -compare_concentration_curves_virus_IGH_paper([classroom_model_IGH_no_mask_windows_closed[0],classroom_model_IGH_with_mask_windows_closed_1HEPA[0], - classroom_model_IGH_with_mask_windows_closed_2HEPA[0]], - labels=['No masks', '1 HEPA filter (2.5 ACH) + surgical masks', '2 HEPA filter (5 ACH) + surgical masks'], - colors=['tomato', 'lightskyblue', 'limegreen', '#1f77b4', 'seagreen', 'lightskyblue', 'deepskyblue'], - title='Combination of measures (windows closed) ' - ) - -print_qd_info(classroom_model_IGH_no_mask_windows_closed[0]) -print_qd_info(classroom_model_IGH_with_mask_windows_closed_1HEPA[0]) -print_qd_info(classroom_model_IGH_with_mask_windows_closed_2HEPA[0]) - #compare_concentration_curves([classroom_model_IGH_no_mask_windows_closed[1],classroom_model_IGH_no_mask_2windows_open_breaks[1],classroom_model_IGH_no_mask_2windows_open_breaks_endOfClass[1], # classroom_model_IGH_no_mask_2windows_open_alltimes[1], classroom_model_IGH_no_mask_windows_closed_1HEPA[1], diff --git a/cara/model_scenarios_IGH_paper.py b/cara/model_scenarios_IGH_paper.py index 84be2561..5acc598b 100644 --- a/cara/model_scenarios_IGH_paper.py +++ b/cara/model_scenarios_IGH_paper.py @@ -3,7 +3,7 @@ from cara.montecarlo import * #IGH paper -# Baseline +# Baseline scenarios classroom_model_IGH_no_mask_windows_closed = [MCExposureModel( concentration_model=MCConcentrationModel( room=models.Room(volume=160), @@ -54,7 +54,7 @@ classroom_model_IGH_no_mask_windows_closed_winter = [MCExposureModel( ) )for qid in (100, 60)] -#Fig 1a +#Fig 1c classroom_model_IGH_no_mask_windows_open_breaks = [MCExposureModel( concentration_model=MCConcentrationModel( room=models.Room(volume=160), @@ -109,7 +109,89 @@ classroom_model_IGH_no_mask_windows_open_alltimes = [MCExposureModel( ) )for qid in (100, 60)] -#Fig 1b +#Fig 1d +classroom_model_IGH_no_mask_windows_open_breaks_winter = [MCExposureModel( + concentration_model=MCConcentrationModel( + room=models.Room(volume=160), + ventilation=models.SlidingWindow( + active=models.SpecificInterval(((1.5, 2), (3.5, 4.5), (6, 6.5))), + inside_temp=models.PiecewiseConstant((0, 24), (295,)), + outside_temp=models.PiecewiseConstant((0, 24), (278,)), + window_height=1.6, opening_length=0.6, + ), + infected=MCInfectedPopulation( + number=1, + presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), + masked=False, + virus=MCVirus(halflife=6.43, qID=qid), + expiratory_activity=2, + samples=200000, + breathing_category=3, + ) + ), + exposed=models.Population( + number=19, + presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), + activity=models.Activity.types['Seated'], + mask=models.Mask.types['No mask'] + ) +)for qid in (100, 60)] + +classroom_model_IGH_no_mask_windows_open_alltimes_winter = [MCExposureModel( + concentration_model=MCConcentrationModel( + room=models.Room(volume=160), + ventilation=models.SlidingWindow( + active=models.PeriodicInterval(period=120, duration=120), + inside_temp=models.PiecewiseConstant((0, 24), (295,)), + outside_temp=models.PiecewiseConstant((0, 24), (278,)), + window_height=1.6, opening_length=0.2, + ), + infected=MCInfectedPopulation( + number=1, + presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), + masked=False, + virus=MCVirus(halflife=6.43, qID=qid), + expiratory_activity=2, + samples=200000, + breathing_category=3, + ) + ), + exposed=models.Population( + number=19, + presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), + activity=models.Activity.types['Seated'], + mask=models.Mask.types['No mask'] + ) +)for qid in (100, 60)] + +classroom_model_IGH_no_mask_windows_fully_open_alltimes_winter = [MCExposureModel( + concentration_model=MCConcentrationModel( + room=models.Room(volume=160), + ventilation=models.SlidingWindow( + active=models.PeriodicInterval(period=120, duration=120), + inside_temp=models.PiecewiseConstant((0, 24), (295,)), + outside_temp=models.PiecewiseConstant((0, 24), (278,)), + window_height=1.6, opening_length=0.6, + ), + infected=MCInfectedPopulation( + number=1, + presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), + masked=False, + virus=MCVirus(halflife=6.43, qID=qid), + expiratory_activity=2, + samples=200000, + breathing_category=3, + ) + ), + exposed=models.Population( + number=19, + presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), + activity=models.Activity.types['Seated'], + mask=models.Mask.types['No mask'] + ) +)for qid in (100, 60)] + +#Fig 1e classroom_model_IGH_no_mask_2windows_open_breaks = [MCExposureModel( concentration_model=MCConcentrationModel( room=models.Room(volume=160), @@ -164,89 +246,7 @@ classroom_model_IGH_no_mask_2windows_open_alltimes = [MCExposureModel( ) )for qid in (100, 60)] -#Fig 1c -classroom_model_IGH_no_mask_windows_open_breaks_winter = [MCExposureModel( - concentration_model=MCConcentrationModel( - room=models.Room(volume=160), - ventilation=models.SlidingWindow( - active=models.SpecificInterval(((1.5, 2), (3.5, 4.5), (6, 6.5))), - inside_temp=models.PiecewiseConstant((0, 24), (295,)), - outside_temp=models.PiecewiseConstant((0, 24), (278,)), - window_height=1.6, opening_length=0.6, - ), - infected=MCInfectedPopulation( - number=1, - presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), - masked=False, - virus=MCVirus(halflife=1.1, qID=qid), - expiratory_activity=2, - samples=200000, - breathing_category=3, - ) - ), - exposed=models.Population( - number=19, - presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), - activity=models.Activity.types['Seated'], - mask=models.Mask.types['No mask'] - ) -)for qid in (100, 60)] - -classroom_model_IGH_no_mask_windows_open_alltimes_winter = [MCExposureModel( - concentration_model=MCConcentrationModel( - room=models.Room(volume=160), - ventilation=models.SlidingWindow( - active=models.PeriodicInterval(period=120, duration=120), - inside_temp=models.PiecewiseConstant((0, 24), (295,)), - outside_temp=models.PiecewiseConstant((0, 24), (278,)), - window_height=1.6, opening_length=0.2, - ), - infected=MCInfectedPopulation( - number=1, - presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), - masked=False, - virus=MCVirus(halflife=1.1, qID=qid), - expiratory_activity=2, - samples=200000, - breathing_category=3, - ) - ), - exposed=models.Population( - number=19, - presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), - activity=models.Activity.types['Seated'], - mask=models.Mask.types['No mask'] - ) -)for qid in (100, 60)] - -classroom_model_IGH_no_mask_windows_fully_open_alltimes_winter = [MCExposureModel( - concentration_model=MCConcentrationModel( - room=models.Room(volume=160), - ventilation=models.SlidingWindow( - active=models.PeriodicInterval(period=120, duration=120), - inside_temp=models.PiecewiseConstant((0, 24), (295,)), - outside_temp=models.PiecewiseConstant((0, 24), (278,)), - window_height=1.6, opening_length=0.6, - ), - infected=MCInfectedPopulation( - number=1, - presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), - masked=False, - virus=MCVirus(halflife=1.1, qID=qid), - expiratory_activity=2, - samples=200000, - breathing_category=3, - ) - ), - exposed=models.Population( - number=19, - presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), - activity=models.Activity.types['Seated'], - mask=models.Mask.types['No mask'] - ) -)for qid in (100, 60)] - -#Fig 1d +#Fig 1f classroom_model_IGH_no_mask_2windows_open_breaks_winter = [MCExposureModel( concentration_model=MCConcentrationModel( room=models.Room(volume=160), @@ -260,7 +260,7 @@ classroom_model_IGH_no_mask_2windows_open_breaks_winter = [MCExposureModel( number=1, presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), masked=False, - virus=MCVirus(halflife=1.1, qID=qid), + virus=MCVirus(halflife=6.43, qID=qid), expiratory_activity=2, samples=200000, breathing_category=3, @@ -287,7 +287,7 @@ classroom_model_IGH_no_mask_2windows_open_alltimes_winter = [MCExposureModel( number=1, presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), masked=False, - virus=MCVirus(halflife=1.1, qID=qid), + virus=MCVirus(halflife=6.43, qID=qid), expiratory_activity=2, samples=200000, breathing_category=3, @@ -314,7 +314,7 @@ classroom_model_IGH_no_mask_2windows_fully_open_alltimes_winter = [MCExposureMod number=1, presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), masked=False, - virus=MCVirus(halflife=1.1, qID=qid), + virus=MCVirus(halflife=6.43, qID=qid), expiratory_activity=2, samples=200000, breathing_category=3, @@ -328,7 +328,7 @@ classroom_model_IGH_no_mask_2windows_fully_open_alltimes_winter = [MCExposureMod ) )for qid in (100, 60)] -#Fig 1e +#Fig 1g classroom_model_IGH_no_mask_6windows_open_breaks = [MCExposureModel( concentration_model=MCConcentrationModel( room=models.Room(volume=160), @@ -383,7 +383,7 @@ classroom_model_IGH_no_mask_6windows_open_alltimes = [MCExposureModel( ) )for qid in (100, 60)] -#Fig 1f +#Fig 1h classroom_model_IGH_no_mask_6windows_open_breaks_winter = [MCExposureModel( concentration_model=MCConcentrationModel( room=models.Room(volume=160), @@ -397,7 +397,7 @@ classroom_model_IGH_no_mask_6windows_open_breaks_winter = [MCExposureModel( number=1, presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), masked=False, - virus=MCVirus(halflife=1.1, qID=qid), + virus=MCVirus(halflife=6.43, qID=qid), expiratory_activity=2, samples=200000, breathing_category=3, @@ -424,7 +424,7 @@ classroom_model_IGH_no_mask_6windows_open_alltimes_winter = [MCExposureModel( number=1, presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), masked=False, - virus=MCVirus(halflife=1.1, qID=qid), + virus=MCVirus(halflife=6.43, qID=qid), expiratory_activity=2, samples=200000, breathing_category=3, @@ -447,6 +447,34 @@ classroom_model_IGH_no_mask_6windows_fully_open_alltimes_winter = [MCExposureMod outside_temp=models.PiecewiseConstant((0, 24), (278,)), window_height=1.6, opening_length=6*0.6, ), + infected=MCInfectedPopulation( + number=1, + presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), + masked=False, + virus=MCVirus(halflife=6.43, qID=qid), + expiratory_activity=2, + samples=200000, + breathing_category=3, + ) + ), + exposed=models.Population( + number=19, + presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), + activity=models.Activity.types['Seated'], + mask=models.Mask.types['No mask'] + ) +)for qid in (100, 60)] + +#Fig 1i +classroom_model_IGH_no_mask_6windows_open_breaks_endOfClass = [MCExposureModel( + concentration_model=MCConcentrationModel( + room=models.Room(volume=160), + ventilation=models.SlidingWindow( + active=models.SpecificInterval(((0.75, 0.92), (1.5, 2), (2.75, 2.92), (3.5, 4.5), (5.25, 5.42), (6, 6.5), (7.25, 7.42))), + inside_temp=models.PiecewiseConstant((0, 24), (295,)), + outside_temp=models.PiecewiseConstant((0, 24), (291,)), + window_height=1.6, opening_length=6*0.6, + ), infected=MCInfectedPopulation( number=1, presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), @@ -465,8 +493,36 @@ classroom_model_IGH_no_mask_6windows_fully_open_alltimes_winter = [MCExposureMod ) )for qid in (100, 60)] +#Fig 1j +classroom_model_IGH_no_mask_6windows_open_breaks_endOfClass_winter = [MCExposureModel( + concentration_model=MCConcentrationModel( + room=models.Room(volume=160), + ventilation=models.SlidingWindow( + active=models.SpecificInterval(((0.75, 0.92), (1.5, 2), (2.75, 2.92), (3.5, 4.5), (5.25, 5.42), (6, 6.5), (7.25, 7.42))), + inside_temp=models.PiecewiseConstant((0, 24), (295,)), + outside_temp=models.PiecewiseConstant((0, 24), (278,)), + window_height=1.6, opening_length=6*0.6, + ), + infected=MCInfectedPopulation( + number=1, + presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), + masked=False, + virus=MCVirus(halflife=6.43, qID=qid), + expiratory_activity=2, + samples=200000, + breathing_category=3, + ) + ), + exposed=models.Population( + number=19, + presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), + activity=models.Activity.types['Seated'], + mask=models.Mask.types['No mask'] + ) +)for qid in (100, 60)] -#Fig 2 + +#Fig 2a classroom_model_IGH_no_mask_windows_closed_1HEPA = [MCExposureModel( concentration_model=MCConcentrationModel( room=models.Room(volume=160), @@ -529,7 +585,70 @@ classroom_model_IGH_no_mask_windows_closed_2HEPA = [MCExposureModel( ) )for qid in (100, 60)] -#Fig 3 +#Fig 2b +classroom_model_IGH_no_mask_windows_closed_1HEPA_winter = [MCExposureModel( + concentration_model=MCConcentrationModel( + room=models.Room(volume=160), + ventilation=models.MultipleVentilation( + ventilations=( + models.AirChange( + active=models.PeriodicInterval(period=120, duration=120), + air_exch=0.25, + ), + models.HEPAFilter(active=models.PeriodicInterval(period=120, duration=120), + q_air_mech=430) + ) + ), + infected=MCInfectedPopulation( + number=1, + presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), + masked=False, + virus=MCVirus(halflife=6.43, qID=qid), + expiratory_activity=2, + samples=200000, + breathing_category=3, + ) + ), + exposed=models.Population( + number=19, + presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), + activity=models.Activity.types['Seated'], + mask=models.Mask.types['No mask'] + ) +)for qid in (100, 60)] + +classroom_model_IGH_no_mask_windows_closed_2HEPA_winter = [MCExposureModel( + concentration_model=MCConcentrationModel( + room=models.Room(volume=160), + ventilation=models.MultipleVentilation( + ventilations=( + models.AirChange( + active=models.PeriodicInterval(period=120, duration=120), + air_exch=0.25, + ), + models.HEPAFilter(active=models.PeriodicInterval(period=120, duration=120), + q_air_mech=860) + ) + ), + infected=MCInfectedPopulation( + number=1, + presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), + masked=False, + virus=MCVirus(halflife=6.43, qID=qid), + expiratory_activity=2, + samples=200000, + breathing_category=3, + ) + ), + exposed=models.Population( + number=19, + presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), + activity=models.Activity.types['Seated'], + mask=models.Mask.types['No mask'] + ) +)for qid in (100, 60)] + +#Fig 2c classroom_model_IGH_with_mask_windows_closed = [MCExposureModel( concentration_model=MCConcentrationModel( room=models.Room(volume=160), @@ -555,48 +674,19 @@ classroom_model_IGH_with_mask_windows_closed = [MCExposureModel( ) )for qid in (100, 60)] -#Fig 4a -classroom_model_IGH_no_mask_2windows_open_breaks_winter_bis = [MCExposureModel( +#Fig 2d +classroom_model_IGH_with_mask_windows_closed_winter = [MCExposureModel( concentration_model=MCConcentrationModel( room=models.Room(volume=160), - ventilation=models.SlidingWindow( - active=models.SpecificInterval(((1.5, 2), (3.5, 4.5), (6, 6.5))), - inside_temp=models.PiecewiseConstant((0, 24), (295,)), - outside_temp=models.PiecewiseConstant((0, 24), (278,)), - window_height=1.6, opening_length=2*0.6, - ), - infected=MCInfectedPopulation( - number=1, - presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), - masked=False, - virus=MCVirus(halflife=1.1, qID=qid), - expiratory_activity=2, - samples=200000, - breathing_category=3, - ) - ), - exposed=models.Population( - number=19, - presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), - activity=models.Activity.types['Seated'], - mask=models.Mask.types['No mask'] - ) -)for qid in (100, 60)] - -classroom_model_IGH_with_mask_2windows_open_breaks_winter = [MCExposureModel( - concentration_model=MCConcentrationModel( - room=models.Room(volume=160), - ventilation=models.SlidingWindow( - active=models.SpecificInterval(((1.5, 2), (3.5, 4.5), (6, 6.5))), - inside_temp=models.PiecewiseConstant((0, 24), (295,)), - outside_temp=models.PiecewiseConstant((0, 24), (278,)), - window_height=1.6, opening_length=2 * 0.6, + ventilation=models.AirChange( + active=models.PeriodicInterval(period=120, duration=120), + air_exch=0.25, ), infected=MCInfectedPopulation( number=1, presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), masked=True, - virus=MCVirus(halflife=1.1, qID=qid), + virus=MCVirus(halflife=6.43, qID=qid), expiratory_activity=2, samples=200000, breathing_category=3, @@ -610,74 +700,7 @@ classroom_model_IGH_with_mask_2windows_open_breaks_winter = [MCExposureModel( ) )for qid in (100, 60)] -#Fig 4b -classroom_model_IGH_with_mask_2windows_open_breaks_winter_1HEPA = [MCExposureModel( - concentration_model=MCConcentrationModel( - room=models.Room(volume=160), - ventilation=models.MultipleVentilation( - ventilations=( - models.SlidingWindow( - active=models.SpecificInterval(((1.5, 2), (3.5, 4.5), (6, 6.5))), - inside_temp=models.PiecewiseConstant((0, 24), (295,)), - outside_temp=models.PiecewiseConstant((0, 24), (278,)), - window_height=1.6, opening_length=2 * 0.6, - ), - models.HEPAFilter(active=models.PeriodicInterval(period=120, duration=120), - q_air_mech=430) - ) - ), - infected=MCInfectedPopulation( - number=1, - presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), - masked=True, - virus=MCVirus(halflife=1.1, qID=qid), - expiratory_activity=2, - samples=200000, - breathing_category=3, - ) - ), - exposed=models.Population( - number=19, - presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), - activity=models.Activity.types['Seated'], - mask=models.Mask.types['Type I'] - ) -)for qid in (100, 60)] - -classroom_model_IGH_with_mask_2windows_open_breaks_winter_2HEPA = [MCExposureModel( - concentration_model=MCConcentrationModel( - room=models.Room(volume=160), - ventilation=models.MultipleVentilation( - ventilations=( - models.SlidingWindow( - active=models.SpecificInterval(((1.5, 2), (3.5, 4.5), (6, 6.5))), - inside_temp=models.PiecewiseConstant((0, 24), (295,)), - outside_temp=models.PiecewiseConstant((0, 24), (278,)), - window_height=1.6, opening_length=2 * 0.6, - ), - models.HEPAFilter(active=models.PeriodicInterval(period=120, duration=120), - q_air_mech=2*430) - ) - ), - infected=MCInfectedPopulation( - number=1, - presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), - masked=True, - virus=MCVirus(halflife=1.1, qID=qid), - expiratory_activity=2, - samples=200000, - breathing_category=3, - ) - ), - exposed=models.Population( - number=19, - presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), - activity=models.Activity.types['Seated'], - mask=models.Mask.types['Type I'] - ) -)for qid in (100, 60)] - -#Fig 4c +#Fig 3b classroom_model_IGH_no_mask_2windows_open_alltimes_winter_bis = [MCExposureModel( concentration_model=MCConcentrationModel( room=models.Room(volume=160), @@ -691,7 +714,7 @@ classroom_model_IGH_no_mask_2windows_open_alltimes_winter_bis = [MCExposureModel number=1, presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), masked=False, - virus=MCVirus(halflife=1.1, qID=qid), + virus=MCVirus(halflife=6.43, qID=qid), expiratory_activity=2, samples=200000, breathing_category=3, @@ -718,294 +741,7 @@ classroom_model_IGH_with_mask_2windows_open_alltimes_winter = [MCExposureModel( number=1, presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), masked=True, - virus=MCVirus(halflife=1.1, qID=qid), - expiratory_activity=2, - samples=200000, - breathing_category=3, - ) - ), - exposed=models.Population( - number=19, - presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), - activity=models.Activity.types['Seated'], - mask=models.Mask.types['Type I'] - ) -)for qid in (100, 60)] - -#Fig 4d -classroom_model_IGH_with_mask_2windows_open_alltimes_winter_1HEPA = [MCExposureModel( - concentration_model=MCConcentrationModel( - room=models.Room(volume=160), - ventilation=models.MultipleVentilation( - ventilations=( - models.SlidingWindow( - active=models.PeriodicInterval(period=120, duration=120), - inside_temp=models.PiecewiseConstant((0, 24), (295,)), - outside_temp=models.PiecewiseConstant((0, 24), (278,)), - window_height=1.6, opening_length=2 * 0.2, - ), - models.HEPAFilter(active=models.PeriodicInterval(period=120, duration=120), - q_air_mech=430) - ) - ), - infected=MCInfectedPopulation( - number=1, - presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), - masked=True, - virus=MCVirus(halflife=1.1, qID=qid), - expiratory_activity=2, - samples=200000, - breathing_category=3, - ) - ), - exposed=models.Population( - number=19, - presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), - activity=models.Activity.types['Seated'], - mask=models.Mask.types['Type I'] - ) -)for qid in (100, 60)] - -classroom_model_IGH_with_mask_2windows_open_alltimes_winter_2HEPA = [MCExposureModel( - concentration_model=MCConcentrationModel( - room=models.Room(volume=160), - ventilation=models.MultipleVentilation( - ventilations=( - models.SlidingWindow( - active=models.PeriodicInterval(period=120, duration=120), - inside_temp=models.PiecewiseConstant((0, 24), (295,)), - outside_temp=models.PiecewiseConstant((0, 24), (278,)), - window_height=1.6, opening_length=2 * 0.2, - ), - models.HEPAFilter(active=models.PeriodicInterval(period=120, duration=120), - q_air_mech=2*430) - ) - ), - infected=MCInfectedPopulation( - number=1, - presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), - masked=True, - virus=MCVirus(halflife=1.1, qID=qid), - expiratory_activity=2, - samples=200000, - breathing_category=3, - ) - ), - exposed=models.Population( - number=19, - presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), - activity=models.Activity.types['Seated'], - mask=models.Mask.types['Type I'] - ) -)for qid in (100, 60)] - -###### - - -classroom_model_IGH_no_mask_windows_open_breaks_endOfClass = [MCExposureModel( - concentration_model=MCConcentrationModel( - room=models.Room(volume=160), - ventilation=models.SlidingWindow( - active=models.SpecificInterval(((0.75, 0.92), (1.5, 2), (2.75, 2.92), (3.5, 4.5), (5.25, 5.42), (6, 6.5), (7.25, 7.42))), - inside_temp=models.PiecewiseConstant((0, 24), (295,)), - outside_temp=models.PiecewiseConstant((0, 24), (291,)), - window_height=1.6, opening_length=0.6, - ), - infected=MCInfectedPopulation( - number=1, - presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), - masked=False, - virus=MCVirus(halflife=1.1, qID=qid), - expiratory_activity=2, - samples=200000, - breathing_category=3, - ) - ), - exposed=models.Population( - number=19, - presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), - activity=models.Activity.types['Seated'], - mask=models.Mask.types['No mask'] - ) -)for qid in (100, 60)] - -classroom_model_IGH_no_mask_2windows_open_breaks_endOfClass = [MCExposureModel( - concentration_model=MCConcentrationModel( - room=models.Room(volume=160), - ventilation=models.SlidingWindow( - active=models.SpecificInterval(((0.75, 0.92), (1.5, 2), (2.75, 2.92), (3.5, 4.5), (5.25, 5.42), (6, 6.5), (7.25, 7.42))), - inside_temp=models.PiecewiseConstant((0, 24), (295,)), - outside_temp=models.PiecewiseConstant((0, 24), (291,)), - window_height=1.6, opening_length=2*0.6, - ), - infected=MCInfectedPopulation( - number=1, - presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), - masked=False, - virus=MCVirus(halflife=1.1, qID=qid), - expiratory_activity=2, - samples=200000, - breathing_category=3, - ) - ), - exposed=models.Population( - number=19, - presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), - activity=models.Activity.types['Seated'], - mask=models.Mask.types['No mask'] - ) -)for qid in (100, 60)] - -classroom_model_IGH_no_mask_6windows_open_breaks_endOfClass = [MCExposureModel( - concentration_model=MCConcentrationModel( - room=models.Room(volume=160), - ventilation=models.SlidingWindow( - active=models.SpecificInterval(((0.75, 0.92), (1.5, 2), (2.75, 2.92), (3.5, 4.5), (5.25, 5.42), (6, 6.5), (7.25, 7.42))), - inside_temp=models.PiecewiseConstant((0, 24), (295,)), - outside_temp=models.PiecewiseConstant((0, 24), (291,)), - window_height=1.6, opening_length=6*0.6, - ), - infected=MCInfectedPopulation( - number=1, - presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), - masked=False, - virus=MCVirus(halflife=1.1, qID=qid), - expiratory_activity=2, - samples=200000, - breathing_category=3, - ) - ), - exposed=models.Population( - number=19, - presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), - activity=models.Activity.types['Seated'], - mask=models.Mask.types['No mask'] - ) -)for qid in (100, 60)] - -classroom_model_IGH_no_mask_windows_open_breaks_endOfClass_winter = [MCExposureModel( - concentration_model=MCConcentrationModel( - room=models.Room(volume=160), - ventilation=models.SlidingWindow( - active=models.SpecificInterval(((0.75, 0.92), (1.5, 2), (2.75, 2.92), (3.5, 4.5), (5.25, 5.42), (6, 6.5), (7.25, 7.42))), - inside_temp=models.PiecewiseConstant((0, 24), (295,)), - outside_temp=models.PiecewiseConstant((0, 24), (278,)), - window_height=1.6, opening_length=0.6, - ), - infected=MCInfectedPopulation( - number=1, - presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), - masked=False, - virus=MCVirus(halflife=1.1, qID=qid), - expiratory_activity=2, - samples=200000, - breathing_category=3, - ) - ), - exposed=models.Population( - number=19, - presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), - activity=models.Activity.types['Seated'], - mask=models.Mask.types['No mask'] - ) -)for qid in (100, 60)] - -classroom_model_IGH_no_mask_2windows_open_breaks_endOfClass_winter = [MCExposureModel( - concentration_model=MCConcentrationModel( - room=models.Room(volume=160), - ventilation=models.SlidingWindow( - active=models.SpecificInterval(((0.75, 0.92), (1.5, 2), (2.75, 2.92), (3.5, 4.5), (5.25, 5.42), (6, 6.5), (7.25, 7.42))), - inside_temp=models.PiecewiseConstant((0, 24), (295,)), - outside_temp=models.PiecewiseConstant((0, 24), (278,)), - window_height=1.6, opening_length=2*0.6, - ), - infected=MCInfectedPopulation( - number=1, - presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), - masked=False, - virus=MCVirus(halflife=1.1, qID=qid), - expiratory_activity=2, - samples=200000, - breathing_category=3, - ) - ), - exposed=models.Population( - number=19, - presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), - activity=models.Activity.types['Seated'], - mask=models.Mask.types['No mask'] - ) -)for qid in (100, 60)] - -classroom_model_IGH_no_mask_6windows_open_breaks_endOfClass_winter = [MCExposureModel( - concentration_model=MCConcentrationModel( - room=models.Room(volume=160), - ventilation=models.SlidingWindow( - active=models.SpecificInterval(((0.75, 0.92), (1.5, 2), (2.75, 2.92), (3.5, 4.5), (5.25, 5.42), (6, 6.5), (7.25, 7.42))), - inside_temp=models.PiecewiseConstant((0, 24), (295,)), - outside_temp=models.PiecewiseConstant((0, 24), (278,)), - window_height=1.6, opening_length=6*0.6, - ), - infected=MCInfectedPopulation( - number=1, - presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), - masked=False, - virus=MCVirus(halflife=1.1, qID=qid), - expiratory_activity=2, - samples=200000, - breathing_category=3, - ) - ), - exposed=models.Population( - number=19, - presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), - activity=models.Activity.types['Seated'], - mask=models.Mask.types['No mask'] - ) -)for qid in (100, 60)] - - -classroom_model_IGH_with_mask_windows_open_breaks_endOfClass = [MCExposureModel( - concentration_model=MCConcentrationModel( - room=models.Room(volume=160), - ventilation=models.SlidingWindow( - active=models.SpecificInterval(((0.75, 0.92), (1.5, 2), (2.75, 2.92), (3.5, 4.5), (5.25, 5.42), (6, 6.5), (7.25, 7.42))), - inside_temp=models.PiecewiseConstant((0, 24), (295,)), - outside_temp=models.PiecewiseConstant((0, 24), (291,)), - window_height=1.6, opening_length=0.6, - ), - infected=MCInfectedPopulation( - number=1, - presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), - masked=True, - virus=MCVirus(halflife=1.1, qID=qid), - expiratory_activity=2, - samples=200000, - breathing_category=3, - ) - ), - exposed=models.Population( - number=19, - presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), - activity=models.Activity.types['Seated'], - mask=models.Mask.types['Type I'] - ) -)for qid in (100, 60)] - -classroom_model_IGH_with_mask_windows_open_alltimes = [MCExposureModel( - concentration_model=MCConcentrationModel( - room=models.Room(volume=160), - ventilation=models.SlidingWindow( - active=models.PeriodicInterval(period=120, duration=120), - inside_temp=models.PiecewiseConstant((0, 24), (295,)), - outside_temp=models.PiecewiseConstant((0, 24), (291,)), - window_height=1.6, opening_length=0.6, - ), - infected=MCInfectedPopulation( - number=1, - presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), - masked=True, - virus=MCVirus(halflife=1.1, qID=qid), + virus=MCVirus(halflife=6.43, qID=qid), expiratory_activity=2, samples=200000, breathing_category=3, @@ -1019,6 +755,7 @@ classroom_model_IGH_with_mask_windows_open_alltimes = [MCExposureModel( ) )for qid in (100, 60)] +#Fig 3c classroom_model_IGH_with_mask_windows_closed_1HEPA = [MCExposureModel( concentration_model=MCConcentrationModel( room=models.Room(volume=160), @@ -1079,4 +816,361 @@ classroom_model_IGH_with_mask_windows_closed_2HEPA = [MCExposureModel( activity=models.Activity.types['Seated'], mask=models.Mask.types['Type I'] ) -)for qid in (100, 60)] \ No newline at end of file +)for qid in (100, 60)] + +#Fig 3d +classroom_model_IGH_with_mask_2windows_open_alltimes_winter_1HEPA = [MCExposureModel( + concentration_model=MCConcentrationModel( + room=models.Room(volume=160), + ventilation=models.MultipleVentilation( + ventilations=( + models.SlidingWindow( + active=models.PeriodicInterval(period=120, duration=120), + inside_temp=models.PiecewiseConstant((0, 24), (295,)), + outside_temp=models.PiecewiseConstant((0, 24), (278,)), + window_height=1.6, opening_length=2 * 0.2, + ), + models.HEPAFilter(active=models.PeriodicInterval(period=120, duration=120), + q_air_mech=430) + ) + ), + infected=MCInfectedPopulation( + number=1, + presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), + masked=True, + virus=MCVirus(halflife=6.43, qID=qid), + expiratory_activity=2, + samples=200000, + breathing_category=3, + ) + ), + exposed=models.Population( + number=19, + presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), + activity=models.Activity.types['Seated'], + mask=models.Mask.types['Type I'] + ) +)for qid in (100, 60)] + +classroom_model_IGH_with_mask_2windows_open_alltimes_winter_2HEPA = [MCExposureModel( + concentration_model=MCConcentrationModel( + room=models.Room(volume=160), + ventilation=models.MultipleVentilation( + ventilations=( + models.SlidingWindow( + active=models.PeriodicInterval(period=120, duration=120), + inside_temp=models.PiecewiseConstant((0, 24), (295,)), + outside_temp=models.PiecewiseConstant((0, 24), (278,)), + window_height=1.6, opening_length=2 * 0.2, + ), + models.HEPAFilter(active=models.PeriodicInterval(period=120, duration=120), + q_air_mech=2*430) + ) + ), + infected=MCInfectedPopulation( + number=1, + presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), + masked=True, + virus=MCVirus(halflife=6.43, qID=qid), + expiratory_activity=2, + samples=200000, + breathing_category=3, + ) + ), + exposed=models.Population( + number=19, + presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), + activity=models.Activity.types['Seated'], + mask=models.Mask.types['Type I'] + ) +)for qid in (100, 60)] + + +###### + +classroom_model_IGH_no_mask_windows_open_breaks_endOfClass = [MCExposureModel( + concentration_model=MCConcentrationModel( + room=models.Room(volume=160), + ventilation=models.SlidingWindow( + active=models.SpecificInterval(((0.75, 0.92), (1.5, 2), (2.75, 2.92), (3.5, 4.5), (5.25, 5.42), (6, 6.5), (7.25, 7.42))), + inside_temp=models.PiecewiseConstant((0, 24), (295,)), + outside_temp=models.PiecewiseConstant((0, 24), (291,)), + window_height=1.6, opening_length=0.6, + ), + infected=MCInfectedPopulation( + number=1, + presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), + masked=False, + virus=MCVirus(halflife=1.1, qID=qid), + expiratory_activity=2, + samples=200000, + breathing_category=3, + ) + ), + exposed=models.Population( + number=19, + presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), + activity=models.Activity.types['Seated'], + mask=models.Mask.types['No mask'] + ) +)for qid in (100, 60)] + +classroom_model_IGH_no_mask_2windows_open_breaks_endOfClass = [MCExposureModel( + concentration_model=MCConcentrationModel( + room=models.Room(volume=160), + ventilation=models.SlidingWindow( + active=models.SpecificInterval(((0.75, 0.92), (1.5, 2), (2.75, 2.92), (3.5, 4.5), (5.25, 5.42), (6, 6.5), (7.25, 7.42))), + inside_temp=models.PiecewiseConstant((0, 24), (295,)), + outside_temp=models.PiecewiseConstant((0, 24), (291,)), + window_height=1.6, opening_length=2*0.6, + ), + infected=MCInfectedPopulation( + number=1, + presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), + masked=False, + virus=MCVirus(halflife=1.1, qID=qid), + expiratory_activity=2, + samples=200000, + breathing_category=3, + ) + ), + exposed=models.Population( + number=19, + presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), + activity=models.Activity.types['Seated'], + mask=models.Mask.types['No mask'] + ) +)for qid in (100, 60)] + +classroom_model_IGH_no_mask_windows_open_breaks_endOfClass_winter = [MCExposureModel( + concentration_model=MCConcentrationModel( + room=models.Room(volume=160), + ventilation=models.SlidingWindow( + active=models.SpecificInterval(((0.75, 0.92), (1.5, 2), (2.75, 2.92), (3.5, 4.5), (5.25, 5.42), (6, 6.5), (7.25, 7.42))), + inside_temp=models.PiecewiseConstant((0, 24), (295,)), + outside_temp=models.PiecewiseConstant((0, 24), (278,)), + window_height=1.6, opening_length=0.6, + ), + infected=MCInfectedPopulation( + number=1, + presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), + masked=False, + virus=MCVirus(halflife=6.43, qID=qid), + expiratory_activity=2, + samples=200000, + breathing_category=3, + ) + ), + exposed=models.Population( + number=19, + presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), + activity=models.Activity.types['Seated'], + mask=models.Mask.types['No mask'] + ) +)for qid in (100, 60)] + +classroom_model_IGH_no_mask_2windows_open_breaks_endOfClass_winter = [MCExposureModel( + concentration_model=MCConcentrationModel( + room=models.Room(volume=160), + ventilation=models.SlidingWindow( + active=models.SpecificInterval(((0.75, 0.92), (1.5, 2), (2.75, 2.92), (3.5, 4.5), (5.25, 5.42), (6, 6.5), (7.25, 7.42))), + inside_temp=models.PiecewiseConstant((0, 24), (295,)), + outside_temp=models.PiecewiseConstant((0, 24), (278,)), + window_height=1.6, opening_length=2*0.6, + ), + infected=MCInfectedPopulation( + number=1, + presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), + masked=False, + virus=MCVirus(halflife=6.43, qID=qid), + expiratory_activity=2, + samples=200000, + breathing_category=3, + ) + ), + exposed=models.Population( + number=19, + presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), + activity=models.Activity.types['Seated'], + mask=models.Mask.types['No mask'] + ) +)for qid in (100, 60)] + + + +classroom_model_IGH_no_mask_2windows_open_breaks_winter_bis = [MCExposureModel( + concentration_model=MCConcentrationModel( + room=models.Room(volume=160), + ventilation=models.SlidingWindow( + active=models.SpecificInterval(((1.5, 2), (3.5, 4.5), (6, 6.5))), + inside_temp=models.PiecewiseConstant((0, 24), (295,)), + outside_temp=models.PiecewiseConstant((0, 24), (278,)), + window_height=1.6, opening_length=2*0.6, + ), + infected=MCInfectedPopulation( + number=1, + presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), + masked=False, + virus=MCVirus(halflife=6.43, qID=qid), + expiratory_activity=2, + samples=200000, + breathing_category=3, + ) + ), + exposed=models.Population( + number=19, + presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), + activity=models.Activity.types['Seated'], + mask=models.Mask.types['No mask'] + ) +)for qid in (100, 60)] + +classroom_model_IGH_with_mask_2windows_open_breaks_winter = [MCExposureModel( + concentration_model=MCConcentrationModel( + room=models.Room(volume=160), + ventilation=models.SlidingWindow( + active=models.SpecificInterval(((1.5, 2), (3.5, 4.5), (6, 6.5))), + inside_temp=models.PiecewiseConstant((0, 24), (295,)), + outside_temp=models.PiecewiseConstant((0, 24), (278,)), + window_height=1.6, opening_length=2 * 0.6, + ), + infected=MCInfectedPopulation( + number=1, + presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), + masked=True, + virus=MCVirus(halflife=6.43, qID=qid), + expiratory_activity=2, + samples=200000, + breathing_category=3, + ) + ), + exposed=models.Population( + number=19, + presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), + activity=models.Activity.types['Seated'], + mask=models.Mask.types['Type I'] + ) +)for qid in (100, 60)] + +classroom_model_IGH_with_mask_2windows_open_breaks_winter_1HEPA = [MCExposureModel( + concentration_model=MCConcentrationModel( + room=models.Room(volume=160), + ventilation=models.MultipleVentilation( + ventilations=( + models.SlidingWindow( + active=models.SpecificInterval(((1.5, 2), (3.5, 4.5), (6, 6.5))), + inside_temp=models.PiecewiseConstant((0, 24), (295,)), + outside_temp=models.PiecewiseConstant((0, 24), (278,)), + window_height=1.6, opening_length=2 * 0.6, + ), + models.HEPAFilter(active=models.PeriodicInterval(period=120, duration=120), + q_air_mech=430) + ) + ), + infected=MCInfectedPopulation( + number=1, + presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), + masked=True, + virus=MCVirus(halflife=6.43, qID=qid), + expiratory_activity=2, + samples=200000, + breathing_category=3, + ) + ), + exposed=models.Population( + number=19, + presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), + activity=models.Activity.types['Seated'], + mask=models.Mask.types['Type I'] + ) +)for qid in (100, 60)] + +classroom_model_IGH_with_mask_2windows_open_breaks_winter_2HEPA = [MCExposureModel( + concentration_model=MCConcentrationModel( + room=models.Room(volume=160), + ventilation=models.MultipleVentilation( + ventilations=( + models.SlidingWindow( + active=models.SpecificInterval(((1.5, 2), (3.5, 4.5), (6, 6.5))), + inside_temp=models.PiecewiseConstant((0, 24), (295,)), + outside_temp=models.PiecewiseConstant((0, 24), (278,)), + window_height=1.6, opening_length=2 * 0.6, + ), + models.HEPAFilter(active=models.PeriodicInterval(period=120, duration=120), + q_air_mech=2*430) + ) + ), + infected=MCInfectedPopulation( + number=1, + presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), + masked=True, + virus=MCVirus(halflife=6.43, qID=qid), + expiratory_activity=2, + samples=200000, + breathing_category=3, + ) + ), + exposed=models.Population( + number=19, + presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), + activity=models.Activity.types['Seated'], + mask=models.Mask.types['Type I'] + ) +)for qid in (100, 60)] + +classroom_model_IGH_with_mask_windows_open_breaks_endOfClass = [MCExposureModel( + concentration_model=MCConcentrationModel( + room=models.Room(volume=160), + ventilation=models.SlidingWindow( + active=models.SpecificInterval(((0.75, 0.92), (1.5, 2), (2.75, 2.92), (3.5, 4.5), (5.25, 5.42), (6, 6.5), (7.25, 7.42))), + inside_temp=models.PiecewiseConstant((0, 24), (295,)), + outside_temp=models.PiecewiseConstant((0, 24), (291,)), + window_height=1.6, opening_length=0.6, + ), + infected=MCInfectedPopulation( + number=1, + presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), + masked=True, + virus=MCVirus(halflife=1.1, qID=qid), + expiratory_activity=2, + samples=200000, + breathing_category=3, + ) + ), + exposed=models.Population( + number=19, + presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), + activity=models.Activity.types['Seated'], + mask=models.Mask.types['Type I'] + ) +)for qid in (100, 60)] + +classroom_model_IGH_with_mask_windows_open_alltimes = [MCExposureModel( + concentration_model=MCConcentrationModel( + room=models.Room(volume=160), + ventilation=models.SlidingWindow( + active=models.PeriodicInterval(period=120, duration=120), + inside_temp=models.PiecewiseConstant((0, 24), (295,)), + outside_temp=models.PiecewiseConstant((0, 24), (291,)), + window_height=1.6, opening_length=0.6, + ), + infected=MCInfectedPopulation( + number=1, + presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), + masked=True, + virus=MCVirus(halflife=1.1, qID=qid), + expiratory_activity=2, + samples=200000, + breathing_category=3, + ) + ), + exposed=models.Population( + number=19, + presence=models.SpecificInterval(((0, 1.5), (2, 3.5), (4.5, 6), (6.5, 8))), + activity=models.Activity.types['Seated'], + mask=models.Mask.types['Type I'] + ) +)for qid in (100, 60)] + + + diff --git a/cara/montecarlo.py b/cara/montecarlo.py index abd222c9..4151b21f 100644 --- a/cara/montecarlo.py +++ b/cara/montecarlo.py @@ -1202,12 +1202,13 @@ def compare_concentration_curves_virus_IGH_paper(exp_models: typing.List[MCExpos times = np.arange(start, stop, TIME_STEP) concentrations = [[np.mean(model.concentration_model.concentration(t))*exp_models[0].concentration_model.virus.qID for t in times] for model in exp_models] - fig, ax = plt.subplots() + fig, ax = plt.subplots(figsize=(8.5,6)) for c, label, color in zip(concentrations, labels, colors): ax.plot(times, c, label=label, color=color) ax.legend(loc='upper left') - ax.set_ylim(ax.get_ylim()[0], ax.get_ylim()[1] * 1.2) + #ax.set_ylim(ax.get_ylim()[0], ax.get_ylim()[1] * 1.2) + ax.set_ylim(ax.get_ylim()[0], 230) ax.spines["right"].set_visible(False) factors = [0.6 * model.exposed.activity.inhalation_rate * (1 - model.exposed.mask.η_inhale) for model in exp_models] @@ -1230,7 +1231,8 @@ def compare_concentration_curves_virus_IGH_paper(exp_models: typing.List[MCExpos ax1.spines["right"].set_linestyle("--") ax1.spines["right"].set_linestyle((0,(1,5))) ax1.set_ylabel('Mean cumulative dose\n(virion)', fontsize=14) - ax1.set_ylim(ax1.get_ylim()[0], ax1.get_ylim()[1] * 1.2) + #ax1.set_ylim(ax1.get_ylim()[0], ax1.get_ylim()[1] * 1.2) + ax1.set_ylim(ax.get_ylim()[0], 300) #ax2 = ax.twinx() #ax2.spines["right"].set_position(("axes", 1.15))