From 793ec38289ef2f79a3606ad3973b88df9e44190e Mon Sep 17 00:00:00 2001
From: "andrea.giannetti" <andrea.giannetti@inaf.it>
Date: Tue, 16 Jul 2024 11:07:49 -0600
Subject: [PATCH] Revised prior for non-detections in density inference.
---
etl/prs/prepare_inference_input.py | 5 ++++-
1 file changed, 4 insertions(+), 1 deletion(-)
diff --git a/etl/prs/prepare_inference_input.py b/etl/prs/prepare_inference_input.py
index 345a5e0..d37e7a6 100644
--- a/etl/prs/prepare_inference_input.py
+++ b/etl/prs/prepare_inference_input.py
@@ -40,7 +40,10 @@ for line_id in ['86', '87', '88', '257', '256', '381', '380']:
intensities_dict[line_id] = {}
uncertainties_dict[line_id] = {}
for row in data_dict:
- intensities_dict[line_id][str(row['source_name'])] = row[f'area_{line_id}']
+ if row[f'area_{line_id}'] >= 3 * row[f'area_unc_{line_id}']:
+ intensities_dict[line_id][str(row['source_name'])] = row[f'area_{line_id}']
+ else:
+ intensities_dict[line_id][str(row['source_name'])] = row[f'area_unc_{line_id}']
uncertainties_dict[line_id][str(row['source_name'])] = row[f'area_unc_{line_id}']
config_dict['measured_integrated_intensities'] = intensities_dict
config_dict['integrated_intensities_uncertainties'] = uncertainties_dict
--
GitLab