Skip to content
Snippets Groups Projects
Commit 793ec382 authored by Andrea Giannetti's avatar Andrea Giannetti
Browse files

Revised prior for non-detections in density inference.

parent 65fb2825
No related branches found
No related tags found
No related merge requests found
Show whitespace changes
Inline Side-by-side
etl/prs/prepare_inference_input.py
+ 4
− 1
View file @ 793ec382
...@@ -40,7 +40,10 @@ for line_id in ['86', '87', '88', '257', '256', '381', '380']: ...@@ -40,7 +40,10 @@ for line_id in ['86', '87', '88', '257', '256', '381', '380']:
intensities_dict[line_id] = {} intensities_dict[line_id] = {}
uncertainties_dict[line_id] = {} uncertainties_dict[line_id] = {}
for row in data_dict: for row in data_dict:
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}'] 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}'] uncertainties_dict[line_id][str(row['source_name'])] = row[f'area_unc_{line_id}']
config_dict['measured_integrated_intensities'] = intensities_dict config_dict['measured_integrated_intensities'] = intensities_dict
config_dict['integrated_intensities_uncertainties'] = uncertainties_dict config_dict['integrated_intensities_uncertainties'] = uncertainties_dict
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment