From d7a604de1a74d17b5d3f2e560e57482d840f51bd Mon Sep 17 00:00:00 2001 From: markus Date: Fri, 26 Feb 2021 12:31:06 +0100 Subject: [PATCH] stop accumulating qD during breaks --- cara/montecarlo.py | 7 ++++++- 1 file changed, 6 insertions(+), 1 deletion(-) diff --git a/cara/montecarlo.py b/cara/montecarlo.py index 94b634a6..58543492 100644 --- a/cara/montecarlo.py +++ b/cara/montecarlo.py @@ -1063,8 +1063,13 @@ def compare_concentration_curves(exp_models: typing.List[MCExposureModel], label ax.set_ylim(ax.get_ylim()[0], ax.get_ylim()[1] * 1.2) factors = [0.6 * model.exposed.activity.inhalation_rate * (1 - model.exposed.mask.η_inhale) for model in exp_models] + present_indexes = np.array([exp_models[0].exposed.person_present(t) for t in times]) + modified_concentrations = [np.array(c) for c in concentrations] + for mc in modified_concentrations: + mc[~present_indexes] = 0 + qds = [[np.trapz(c[:i + 1], times[:i + 1]) * factor for i in range(len(times))] - for c, factor in zip(concentrations, factors)] + for c, factor in zip(modified_concentrations, factors)] plt.suptitle("") plt.xlabel("Exposure time ($h$)")