Adjust uncertainty for initial "frozen" maps

imap_processing/hi/hi_l2.py

@@ -37,8 +37,8 @@
 FULL_EXPOSURE_TIME_AVERAGE_SET = {"bg_rate", "bg_rate_sys_err", "obs_date", "energy_sc"}
 
 # Calibration systematic uncertainty as a fraction of intensity.
-# This represents the current calibration uncertainty (22%).
-CALIBRATION_UNCERTAINTY_FRACTION = 0.22
+# This represents the current calibration uncertainty (40%).
+CALIBRATION_UNCERTAINTY_FRACTION = 0.40
 
 
 # =============================================================================

imap_processing/hi/utils.py

@@ -141,7 +141,7 @@ class HiConstants:
     # ESAs 7, 8, 9 get an extra 0.0025/s uncertainty to account for possible
     # unidentified additional background in these ESA steps.
     UPPER_ESA_EXTRA_BACKGROUND_UNC = xr.DataArray(
-        np.array([0.0025, 0.0025, 0.0025]),
+        np.array([0.0025, 0.0025, 0.0055]),
         dims=["esa_energy_step"],
         coords={"esa_energy_step": np.array([7, 8, 9], dtype=int)},
     )

imap_processing/tests/hi/test_hi_l1c.py

@@ -1248,21 +1248,29 @@ def test_pset_backgrounds_esa_7_8_9_extra_uncertainty(mock_compute_background_co
 
     # With 0 counts, Poisson uncertainty is 0, so the only uncertainties are:
     # - EXCESS_BACKGROUND_COUNT_RATE_UNC = 0.001 (for all ESAs)
-    # - UPPER_ESA_EXTRA_BACKGROUND_UNC = 0.0025 (for ESAs 7,8,9 only)
+    # - UPPER_ESA_EXTRA_BACKGROUND_UNC = 0.0025 for ESAs 7 and 8, 0.0055 for
+    #   ESA 9
 
     # Expected uncertainty for ESA 1 (low ESA, no extra uncertainty):
     # sqrt(0 + 0.001^2 + 0) = 0.001
     expected_unc_esa1 = np.sqrt(
         0**2 + HiConstants.EXCESS_BACKGROUND_COUNT_RATE_UNC**2 + 0**2
     )
 
-    # Expected uncertainty for ESAs 7 and 9 (high ESAs, with extra uncertainty):
+    # Expected uncertainty for ESAs 7 and 8 (high ESAs, with extra uncertainty):
     # sqrt(0 + 0.001^2 + 0.0025^2) = sqrt(0.000001 + 0.00000625) ≈ 0.002693
-    expected_unc_esa7_9 = np.sqrt(
+    expected_unc_esa7 = np.sqrt(
         0**2
         + HiConstants.EXCESS_BACKGROUND_COUNT_RATE_UNC**2
         + 0.0025**2  # UPPER_ESA_EXTRA_BACKGROUND_UNC
     )
+    # Expected uncertainty for ESA 9:
+    # sqrt(0 + 0.001^2 + 0.0055^2) = sqrt(0.000001 + 0.00000625) ≈ 0.002693
+    expected_unc_esa9 = np.sqrt(
+        0**2
+        + HiConstants.EXCESS_BACKGROUND_COUNT_RATE_UNC**2
+        + 0.0055**2  # UPPER_ESA_EXTRA_BACKGROUND_UNC
+    )
 
     # ESA 1 is at index 0, ESA 7 at index 1, ESA 9 at index 2 in the output
     # Use isel to select by position
@@ -1278,15 +1286,15 @@ def test_pset_backgrounds_esa_7_8_9_extra_uncertainty(mock_compute_background_co
     unc_esa7 = unc_result.isel(esa_energy_step=1).values
     np.testing.assert_allclose(
         unc_esa7,
-        expected_unc_esa7_9,
+        expected_unc_esa7,
         rtol=1e-6,
         err_msg="ESA 7 uncertainty should include UPPER_ESA_EXTRA_BACKGROUND_UNC",
     )
 
     unc_esa9 = unc_result.isel(esa_energy_step=2).values
     np.testing.assert_allclose(
         unc_esa9,
-        expected_unc_esa7_9,
+        expected_unc_esa9,
         rtol=1e-6,
         err_msg="ESA 9 uncertainty should include UPPER_ESA_EXTRA_BACKGROUND_UNC",
     )