From 086bcdaabbcb714cc40af8bd70b1a6e4fe5c0af7 Mon Sep 17 00:00:00 2001 From: Kris Thielemans Date: Tue, 24 Sep 2024 14:22:38 +0100 Subject: [PATCH] use SensitivityModel for multiplication to cope with TOF --- SIRF_data_preparation/data_QC.py | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) diff --git a/SIRF_data_preparation/data_QC.py b/SIRF_data_preparation/data_QC.py index 9176da6..784e88d 100644 --- a/SIRF_data_preparation/data_QC.py +++ b/SIRF_data_preparation/data_QC.py @@ -178,7 +178,9 @@ def main(argv=None): acquired_data = STIR.AcquisitionData(os.path.join(srcdir, 'prompts.hs')) additive_term = STIR.AcquisitionData(os.path.join(srcdir, 'additive_term.hs')) mult_factors = STIR.AcquisitionData(os.path.join(srcdir, 'mult_factors.hs')) - background = additive_term * mult_factors + sensitivity_factors = STIR.AcquisitionSensitivityModel(mult_factors) + sensitivity_factors.set_up(acquired_data) + background = sensitivity_factors.forward(additive_term) plot_sinogram_profile(acquired_data, background, srcdir=srcdir) OSEM_image = plot_image_if_exists(os.path.join(srcdir, 'OSEM_image'), **slices)