From 52020a3296fb0d425f15593469fcc0829b9fcad9 Mon Sep 17 00:00:00 2001 From: Jon Niehof Date: Wed, 14 Jan 2026 11:53:59 -0700 Subject: [PATCH] Harrison L3 - JNIE - IMAP Lo: Metadata updates for background-subtracted ISN --- .../cdf/config/imap_lo_l3_variable_attrs.yaml | 60 ++++++------------- imap_l3_processing/lo/lo_processor.py | 6 +- 2 files changed, 22 insertions(+), 44 deletions(-) diff --git a/imap_l3_processing/cdf/config/imap_lo_l3_variable_attrs.yaml b/imap_l3_processing/cdf/config/imap_lo_l3_variable_attrs.yaml index 67110cd8..f01484e4 100644 --- a/imap_l3_processing/cdf/config/imap_lo_l3_variable_attrs.yaml +++ b/imap_l3_processing/cdf/config/imap_lo_l3_variable_attrs.yaml @@ -609,29 +609,6 @@ bg_rates_sys_err: LABL_PTR_2: longitude_label LABL_PTR_3: latitude_label SCALETYP: linear -bg_subtracted_stat_err: - NAME: bg_subtracted_stat_err - DATA_TYPE: CDF_REAL4 - CATDESC: The statisical err of the background subtraction? - DEPEND_0: epoch - DEPEND_1: energy - DEPEND_2: longitude - DEPEND_3: latitude - VAR_TYPE: data - RECORD_VARYING: RV - DISPLAY_TYPE: map_image - VARIABLE_PURPOSE: PRIMARY_VAR,SUMMARY - FIELDNAM: BG subtracted stat err - FORMAT: F12.1 - LABLAXIS: BG Statistical Error - UNITS: cm -2 s -1 sr -1 keV -1 - VALIDMIN: 0.0 - VALIDMAX: 1.00E+9 - FILLVAL: -1.00E+31 - LABL_PTR_1: energy_label - LABL_PTR_2: longitude_label - LABL_PTR_3: latitude_label - SCALETYP: linear ena_count_rate: NAME: ena_count_rate DATA_TYPE: CDF_REAL4 @@ -680,10 +657,10 @@ ena_count_rate_stat_uncert: LABL_PTR_2: longitude_label LABL_PTR_3: latitude_label SCALETYP: linear -isn_bg_rate_subtracted: - NAME: isn_bg_rate_subtracted +isn_rate_bg_subtracted: + NAME: isn_rate_bg_subtracted DATA_TYPE: CDF_REAL4 - CATDESC: ISN Background rate. + CATDESC: ISN rate, background subtracted. DEPEND_0: epoch DEPEND_1: energy DEPEND_2: longitude @@ -692,23 +669,24 @@ isn_bg_rate_subtracted: RECORD_VARYING: RV DISPLAY_TYPE: map_image VARIABLE_PURPOSE: PRIMARY_VAR,SUMMARY - FIELDNAM: Counts + FIELDNAM: ISN Rate bg sub FORMAT: F12.1 - LABLAXIS: ISN BG Rate + LABLAXIS: ISN Rate UNITS: count s-1 - VALIDMIN: 0.0 + VALIDMIN: -1.00E+9 VALIDMAX: 1.00E+9 FILLVAL: -1.00E+31 LABL_PTR_1: energy_label LABL_PTR_2: longitude_label LABL_PTR_3: latitude_label - DELTA_PLUS_VAR: isn_bg_rate_subtracted_stat_err - DELTA_MINUS_VAR: isn_bg_rate_subtracted_stat_err + DELTA_PLUS_VAR: isn_rate_bg_subtracted_stat_uncert + DELTA_MINUS_VAR: isn_rate_bg_subtracted_stat_uncert SCALETYP: linear -isn_bg_rate_subtracted_stat_err: - NAME: isn_bg_rate_subtracted_stat_err + VAR_NOTE: Negative values indicate poor counting statistics and should not be ignored; average over time. +isn_rate_bg_subtracted_stat_uncert: + NAME: isn_rate_bg_subtracted_stat_uncert DATA_TYPE: CDF_REAL4 - CATDESC: ISN Background rate subtracted err. + CATDESC: ISN rate, background subtracted statistical uncertainty. DEPEND_0: epoch DEPEND_1: energy DEPEND_2: longitude @@ -717,9 +695,9 @@ isn_bg_rate_subtracted_stat_err: RECORD_VARYING: RV DISPLAY_TYPE: map_image VARIABLE_PURPOSE: PRIMARY_VAR,SUMMARY - FIELDNAM: ISN Background rate subtracted error + FIELDNAM: ISN rate bg sub stat unc FORMAT: F12.1 - LABLAXIS: ISN BG Subtracted Rate + LABLAXIS: Statistical Unc. UNITS: count s-1 VALIDMIN: 0.0 VALIDMAX: 1.00E+9 @@ -728,10 +706,10 @@ isn_bg_rate_subtracted_stat_err: LABL_PTR_2: longitude_label LABL_PTR_3: latitude_label SCALETYP: linear -isn_bg_rate_subtracted_sys_uncert: - NAME: isn_bg_rate_subtracted_sys_uncert +isn_rate_bg_subtracted_sys_err: + NAME: isn_rate_bg_subtracted_sys_err DATA_TYPE: CDF_REAL4 - CATDESC: ISN Background rate systematic uncertainty + CATDESC: ISN rate, background subtracted non-statistical error. DEPEND_0: epoch DEPEND_1: energy DEPEND_2: longitude @@ -740,9 +718,9 @@ isn_bg_rate_subtracted_sys_uncert: RECORD_VARYING: RV DISPLAY_TYPE: map_image VARIABLE_PURPOSE: PRIMARY_VAR,SUMMARY - FIELDNAM: ISN Background subtracted rate systematic uncertainty + FIELDNAM: ISN rate bg sub sys err FORMAT: F12.1 - LABLAXIS: ISN BG Subtracted Rate Systematic Uncertainty + LABLAXIS: Non-statistical Error UNITS: count s-1 VALIDMIN: 0.0 VALIDMAX: 1.00E+9 diff --git a/imap_l3_processing/lo/lo_processor.py b/imap_l3_processing/lo/lo_processor.py index 9459b2aa..634b7a19 100644 --- a/imap_l3_processing/lo/lo_processor.py +++ b/imap_l3_processing/lo/lo_processor.py @@ -117,10 +117,10 @@ def isn_background_subtraction(isn_rate_data: ISNRateData) -> ISNBackgroundSubtr ena_count_rate_sys_uncert=np.zeros_like(isn_rate_data.ena_count_rate), latitude=isn_rate_data.latitude, longitude=isn_rate_data.longitude, - isn_bg_rate_subtracted=isn_rate_background_subtracted, - isn_bg_rate_subtracted_stat_err=np.sqrt(np.square(isn_rate_data.ena_count_rate_stat_uncert) + np.square( + isn_rate_bg_subtracted=isn_rate_background_subtracted, + isn_rate_bg_subtracted_stat_unc=np.sqrt(np.square(isn_rate_data.ena_count_rate_stat_uncert) + np.square( isn_rate_data.bg_rates_stat_uncert)), - isn_bg_rate_subtracted_sys_uncert=isn_rate_data.bg_rates_sys_err + isn_rate_bg_subtracted_sys_err=isn_rate_data.bg_rates_sys_err ) return ISNBackgroundSubtractedMapData(isn_rate_map_data=map_data)