imap-processing 0.18.0__py3-none-any.whl → 0.19.2__py3-none-any.whl

imap_processing/idex/idex_l2b.py

@@ -17,8 +17,9 @@ Examples
     l1a_data, l1a_evt_data, l1b_evt_data = PacketParser(l0_file)
     l1b_data = idex_l1b(l1a_data)
     l1a_data = idex_l2a(l1b_data)
-    l2b_data = idex_l2b(l2a_data, [evt_data])
-    write_cdf(l2b_data)
+    l2b_and_l2c_datasets = idex_l2b(l2a_data, [evt_data])
+    write_cdf(l2b_and_l2c_datasets[0])
+    write_cdf(l2b_and_l2c_datasets[1])
 """
 
 import collections
@@ -29,9 +30,11 @@ from datetime import datetime, timedelta
 import numpy as np
 import xarray as xr
 
+from imap_processing.ena_maps.ena_maps import SkyTilingType
 from imap_processing.ena_maps.utils.spatial_utils import AzElSkyGrid
 from imap_processing.idex.idex_constants import (
     FG_TO_KG,
+    IDEX_EVENT_REFERENCE_FRAME,
     IDEX_SPACING_DEG,
     SECONDS_IN_DAY,
     IDEXEvtAcquireCodes,
@@ -81,9 +84,12 @@ LAT_BINS_EDGES = SKY_GRID.el_bin_edges
 
 def idex_l2b(
     l2a_datasets: list[xr.Dataset], evt_datasets: list[xr.Dataset]
-) -> xr.Dataset:
+) -> list[xr.Dataset]:
     """
-    Will process IDEX l2a data to create l2b data products.
+    Will process IDEX l2a data to create l2b and l2c data products.
+
+    IDEX L2B processing creates L2b and L2c at the same time because L2c needs no
+    additional dependencies and is a natural extension of L2b processing.
 
     Parameters
     ----------
@@ -94,16 +100,18 @@ def idex_l2b(
 
     Returns
     -------
-    l2b_dataset : xarray.Dataset
-        The``xarray`` dataset containing the science data and supporting metadata.
+    list[xarray.Dataset]
+        The``xarray`` datasets containing the l2b and l2c science data and supporting
+        metadata.
     """
     logger.info(
-        f"Running IDEX L2B processing on datasets: "
-        f"{[ds.attrs['Logical_source'] for ds in l2a_datasets]}"
+        "Running IDEX L2B and L2C processing on L2a datasets. NOTE: L2C datasets are "
+        "processed at the same time as L2B datasets because L2C needs no additional "
+        "dependencies."
     )
-
     # create the attribute manager for this data level
-    idex_attrs = get_idex_attrs("l2b")
+    idex_l2b_attrs = get_idex_attrs("l2b")
+    idex_l2c_attrs = get_idex_attrs("l2c")
     evt_dataset = xr.concat(evt_datasets, dim="epoch")
 
     # Concat all the l2a datasets together
@@ -133,125 +141,138 @@ def idex_l2b(
         daily_on_percentage,
     )
     # Create l2b Dataset
-    charge_bins = np.arange(len(CHARGE_BIN_EDGES) - 1)
-    mass_bins = np.arange(len(CHARGE_BIN_EDGES) - 1)
-    spin_phase_bins = np.arange(len(SPIN_PHASE_BIN_EDGES) - 1)
+    charge_bin_means = np.sqrt(CHARGE_BIN_EDGES[:-1] * CHARGE_BIN_EDGES[1:])
+    mass_bin_means = np.sqrt(MASS_BIN_EDGES[:-1] * MASS_BIN_EDGES[1:])
+    spin_phase_means = (SPIN_PHASE_BIN_EDGES[:-1] + SPIN_PHASE_BIN_EDGES[1:]) / 2
+
+    # Define xarrays that are shared between l2b and l2c
     epoch = xr.DataArray(
         name="epoch",
         data=daily_epoch,
         dims="epoch",
-        attrs=idex_attrs.get_variable_attributes("epoch", check_schema=False),
+        attrs=idex_l2b_attrs.get_variable_attributes("epoch", check_schema=False),
     )
-    data_vars = {
+
+    common_vars = {
         "impact_day_of_year": xr.DataArray(
             name="impact_day_of_year",
             data=epoch_doy_unique,
             dims="epoch",
-            attrs=idex_attrs.get_variable_attributes("impact_day_of_year"),
-        ),
-        "rate_calculation_quality_flags": xr.DataArray(
-            name="rate_calculation_quality_flags",
-            data=rate_quality_flags,
-            dims="epoch",
-            attrs=idex_attrs.get_variable_attributes("rate_calculation_quality_flags"),
+            attrs=idex_l2b_attrs.get_variable_attributes("impact_day_of_year"),
         ),
         "charge_labels": xr.DataArray(
             name="impact_charge_labels",
-            data=charge_bins.astype(str),
+            data=charge_bin_means.astype(str),
             dims="impact_charge",
-            attrs=idex_attrs.get_variable_attributes(
+            attrs=idex_l2b_attrs.get_variable_attributes(
                 "charge_labels", check_schema=False
             ),
         ),
-        "spin_phase_labels": xr.DataArray(
-            name="spin_phase_labels",
-            data=spin_phase_bins.astype(str),
-            dims="spin_phase",
-            attrs=idex_attrs.get_variable_attributes(
-                "spin_phase_labels", check_schema=False
-            ),
-        ),
         "mass_labels": xr.DataArray(
             name="mass_labels",
-            data=mass_bins.astype(str),
+            data=mass_bin_means.astype(str),
             dims="mass",
-            attrs=idex_attrs.get_variable_attributes("mass_labels", check_schema=False),
-        ),
-        "rectangular_lon_pixel_label": xr.DataArray(
-            name="rectangular_lon_pixel_label",
-            data=SKY_GRID.az_bin_midpoints.astype(str),
-            dims="rectangular_lon_pixel",
-            attrs=idex_attrs.get_variable_attributes(
-                "rectangular_lon_pixel_label", check_schema=False
-            ),
-        ),
-        "rectangular_lat_pixel_label": xr.DataArray(
-            name="rectangular_lat_pixel_label",
-            data=SKY_GRID.el_bin_midpoints.astype(str),
-            dims="rectangular_lat_pixel",
-            attrs=idex_attrs.get_variable_attributes(
-                "rectangular_lat_pixel_label", check_schema=False
+            attrs=idex_l2b_attrs.get_variable_attributes(
+                "mass_labels", check_schema=False
             ),
         ),
         "impact_charge": xr.DataArray(
             name="impact_charge",
-            data=charge_bins,
+            data=charge_bin_means,
             dims="impact_charge",
-            attrs=idex_attrs.get_variable_attributes(
+            attrs=idex_l2b_attrs.get_variable_attributes(
                 "impact_charge", check_schema=False
             ),
         ),
         "mass": xr.DataArray(
             name="mass",
-            data=mass_bins,
+            data=mass_bin_means,
             dims="mass",
-            attrs=idex_attrs.get_variable_attributes("mass", check_schema=False),
+            attrs=idex_l2b_attrs.get_variable_attributes("mass", check_schema=False),
         ),
+    }
+    l2b_vars = common_vars | {
         "spin_phase": xr.DataArray(
             name="spin_phase",
-            data=spin_phase_bins,
+            data=spin_phase_means,
             dims="spin_phase",
-            attrs=idex_attrs.get_variable_attributes("spin_phase", check_schema=False),
+            attrs=idex_l2b_attrs.get_variable_attributes(
+                "spin_phase", check_schema=False
+            ),
         ),
-        "rectangular_lon_pixel": xr.DataArray(
-            name="rectangular_lon_pixel",
-            data=SKY_GRID.az_bin_midpoints,
-            dims="rectangular_lon_pixel",
-            attrs=idex_attrs.get_variable_attributes(
-                "rectangular_lon_pixel", check_schema=False
+        "spin_phase_labels": xr.DataArray(
+            name="spin_phase_labels",
+            data=spin_phase_means.astype(str),
+            dims="spin_phase",
+            attrs=idex_l2b_attrs.get_variable_attributes(
+                "spin_phase_labels", check_schema=False
             ),
         ),
-        "rectangular_lat_pixel": xr.DataArray(
-            name="rectangular_lat_pixel",
-            data=SKY_GRID.el_bin_midpoints,
-            dims="rectangular_lat_pixel",
-            attrs=idex_attrs.get_variable_attributes(
-                "rectangular_lat_pixel", check_schema=False
+        "rate_calculation_quality_flags": xr.DataArray(
+            name="rate_calculation_quality_flags",
+            data=rate_quality_flags,
+            dims="epoch",
+            attrs=idex_l2b_attrs.get_variable_attributes(
+                "rate_calculation_quality_flags"
             ),
         ),
         "counts_by_charge": xr.DataArray(
             name="counts_by_charge",
             data=counts_by_charge.astype(np.int64),
             dims=("epoch", "impact_charge", "spin_phase"),
-            attrs=idex_attrs.get_variable_attributes("counts_by_charge"),
+            attrs=idex_l2b_attrs.get_variable_attributes("counts_by_charge"),
         ),
         "counts_by_mass": xr.DataArray(
             name="counts_by_mass",
             data=counts_by_mass.astype(np.int64),
             dims=("epoch", "mass", "spin_phase"),
-            attrs=idex_attrs.get_variable_attributes("counts_by_mass"),
+            attrs=idex_l2b_attrs.get_variable_attributes("counts_by_mass"),
         ),
         "rate_by_charge": xr.DataArray(
             name="rate_by_charge",
             data=rate_by_charge,
             dims=("epoch", "impact_charge", "spin_phase"),
-            attrs=idex_attrs.get_variable_attributes("rate_by_charge"),
+            attrs=idex_l2b_attrs.get_variable_attributes("rate_by_charge"),
         ),
         "rate_by_mass": xr.DataArray(
             name="rate_by_mass",
             data=rate_by_mass,
             dims=("epoch", "mass", "spin_phase"),
-            attrs=idex_attrs.get_variable_attributes("rate_by_mass"),
+            attrs=idex_l2b_attrs.get_variable_attributes("rate_by_mass"),
+        ),
+    }
+    l2c_vars = common_vars | {
+        "rectangular_lon_pixel_label": xr.DataArray(
+            name="rectangular_lon_pixel_label",
+            data=SKY_GRID.az_bin_midpoints.astype(str),
+            dims="rectangular_lon_pixel",
+            attrs=idex_l2c_attrs.get_variable_attributes(
+                "rectangular_lon_pixel_label", check_schema=False
+            ),
+        ),
+        "rectangular_lat_pixel_label": xr.DataArray(
+            name="rectangular_lat_pixel_label",
+            data=SKY_GRID.el_bin_midpoints.astype(str),
+            dims="rectangular_lat_pixel",
+            attrs=idex_l2c_attrs.get_variable_attributes(
+                "rectangular_lat_pixel_label", check_schema=False
+            ),
+        ),
+        "rectangular_lon_pixel": xr.DataArray(
+            name="rectangular_lon_pixel",
+            data=SKY_GRID.az_bin_midpoints,
+            dims="rectangular_lon_pixel",
+            attrs=idex_l2c_attrs.get_variable_attributes(
+                "rectangular_lon_pixel", check_schema=False
+            ),
+        ),
+        "rectangular_lat_pixel": xr.DataArray(
+            name="rectangular_lat_pixel",
+            data=SKY_GRID.el_bin_midpoints,
+            dims="rectangular_lat_pixel",
+            attrs=idex_l2c_attrs.get_variable_attributes(
+                "rectangular_lat_pixel", check_schema=False
+            ),
         ),
         "counts_by_charge_map": xr.DataArray(
             name="counts_by_charge_map",
@@ -262,7 +283,7 @@ def idex_l2b(
                 "rectangular_lon_pixel",
                 "rectangular_lat_pixel",
             ),
-            attrs=idex_attrs.get_variable_attributes("counts_by_charge_map"),
+            attrs=idex_l2c_attrs.get_variable_attributes("counts_by_charge_map"),
         ),
         "counts_by_mass_map": xr.DataArray(
             name="counts_by_mass_map",
@@ -273,7 +294,7 @@ def idex_l2b(
                 "rectangular_lon_pixel",
                 "rectangular_lat_pixel",
             ),
-            attrs=idex_attrs.get_variable_attributes("counts_by_mass_map"),
+            attrs=idex_l2c_attrs.get_variable_attributes("counts_by_mass_map"),
         ),
         "rate_by_charge_map": xr.DataArray(
             name="rate_by_charge_map",
@@ -284,7 +305,7 @@ def idex_l2b(
                 "rectangular_lon_pixel",
                 "rectangular_lat_pixel",
             ),
-            attrs=idex_attrs.get_variable_attributes("rate_by_charge_map"),
+            attrs=idex_l2c_attrs.get_variable_attributes("rate_by_charge_map"),
         ),
         "rate_by_mass_map": xr.DataArray(
             name="rate_by_mass_map",
@@ -295,18 +316,31 @@ def idex_l2b(
                 "rectangular_lon_pixel",
                 "rectangular_lat_pixel",
             ),
-            attrs=idex_attrs.get_variable_attributes("rate_by_mass_map"),
+            attrs=idex_l2c_attrs.get_variable_attributes("rate_by_mass_map"),
         ),
     }
+
     l2b_dataset = xr.Dataset(
         coords={"epoch": epoch},
-        data_vars=data_vars,
-        attrs=idex_attrs.get_global_attributes("imap_idex_l2b_sci"),
+        data_vars=l2b_vars,
+        attrs=idex_l2b_attrs.get_global_attributes("imap_idex_l2b_sci"),
+    )
+    l2c_dataset = xr.Dataset(
+        coords={"epoch": epoch},
+        data_vars=l2c_vars,
     )
+    # Add map attributes
+    map_attrs = {
+        "sky_tiling_type": SkyTilingType.RECTANGULAR.value,
+        "Spacing_degrees": str(IDEX_SPACING_DEG),
+        "Spice_reference_frame": IDEX_EVENT_REFERENCE_FRAME.name,
+    } | idex_l2c_attrs.get_global_attributes("imap_idex_l2c_sci-rectangular")
 
-    logger.info("IDEX L2B science data processing completed.")
+    l2c_dataset.attrs.update(map_attrs)
 
-    return l2b_dataset
+    logger.info("IDEX L2B and L2C science data processing completed.")
+
+    return [l2b_dataset, l2c_dataset]
 
 
 def compute_counts_by_charge_and_mass(
@@ -329,37 +363,11 @@ def compute_counts_by_charge_and_mass(
         dataset, Two 4D arrays containing counts by charge or mass, and by lon and lat
         for each dataset, and a 1D array of daily epoch values.
     """
-    # Initialize arrays to hold counts.
-    # There should be 4 spin phase bins, 10 charge bins, and 10 mass bins.
-    # The first bin for charge and mass is for values below the first bin edge.
-    counts_by_charge = np.zeros(
-        (
-            len(epoch_doy_unique),
-            len(CHARGE_BIN_EDGES) - 1,
-            len(SPIN_PHASE_BIN_EDGES) - 1,
-        ),
-    )
-    counts_by_mass = np.zeros(
-        (len(epoch_doy_unique), len(MASS_BIN_EDGES) - 1, len(SPIN_PHASE_BIN_EDGES) - 1),
-    )
-    # Initialize arrays to hold count maps. Each map is a 3 or 4D array with shape
-    # (epoch, 10 [charge or mass], 60 [longitude bins], 30 [latitude bins]).
-    counts_by_charge_map = np.zeros(
-        (
-            len(epoch_doy_unique),
-            len(CHARGE_BIN_EDGES) - 1,
-            len(LON_BINS_EDGES) - 1,
-            len(LAT_BINS_EDGES) - 1,
-        ),
-    )
-    counts_by_mass_map = np.zeros(
-        (
-            len(epoch_doy_unique),
-            len(MASS_BIN_EDGES) - 1,
-            len(LON_BINS_EDGES) - 1,
-            len(LAT_BINS_EDGES) - 1,
-        ),
-    )
+    # Initialize lists to hold counts.
+    counts_by_charge = []
+    counts_by_mass = []
+    counts_by_charge_map = []
+    counts_by_mass_map = []
     daily_epoch = np.zeros(len(epoch_doy_unique), dtype=np.float64)
     for i in range(len(epoch_doy_unique)):
         doy = epoch_doy_unique[i]
@@ -379,43 +387,44 @@ def compute_counts_by_charge_and_mass(
         latitude = l2a_dataset["latitude"].data[current_day_indices]
         # Convert units
         mass_vals = FG_TO_KG * np.atleast_1d(mass_vals)
-        # Bin masses
-        binned_mass = np.asarray(np.digitize(mass_vals, bins=MASS_BIN_EDGES))
-        # Bin charges
-        binned_charge = np.asarray(np.digitize(charge_vals, bins=CHARGE_BIN_EDGES))
         # Bin spin phases
         binned_spin_phase = bin_spin_phases(spin_phase_angles)
-        # Bin longitude and latitude into the rectangular grid.
-        binned_longitude = np.asarray(np.digitize(longitude, bins=LON_BINS_EDGES))
+        # Clip arrays to ensure that the values are within the valid range of bins.
         # Latitude should be binned with the right edge included. 90 is a valid latitude
-        binned_latitude = np.asarray(np.digitize(latitude, bins=LAT_BINS_EDGES))
-        # Clip latitude value above the right edge to be in the last bin
-        binned_latitude = np.clip(binned_latitude, 1, len(LAT_BINS_EDGES) - 1)
-        # If the values in the array are beyond the bounds of bins, 0 or len(bins) it is
-        # returned as such. In this case, the desired result is to place the values
-        # beyond the first or last bin into the first or last bin, respectively.
-        binned_charge = np.clip(binned_charge, 1, len(CHARGE_BIN_EDGES) - 1)
-        binned_mass = np.clip(binned_mass, 1, len(MASS_BIN_EDGES) - 1)
-
-        # Count dust events for each spin phase, mass bin, charge bin, and bin into
-        # a rectangular grid
-        for mass_bin, charge_bin, spin_phase_bin, lon_bin, lat_bin in zip(
-            binned_mass,
-            binned_charge,
-            binned_spin_phase,
-            binned_longitude,
-            binned_latitude,
-        ):
-            counts_by_mass[i, mass_bin - 1, spin_phase_bin] += 1
-            counts_by_charge[i, charge_bin - 1, spin_phase_bin] += 1
-            counts_by_mass_map[i, mass_bin - 1, lon_bin - 1, lat_bin - 1] += 1
-            counts_by_charge_map[i, charge_bin - 1, lon_bin - 1, lat_bin - 1] += 1
+        latitude = np.clip(latitude, -90, 90)
+        mass_vals = np.clip(mass_vals, MASS_BIN_EDGES[0], MASS_BIN_EDGES[-1])
+        charge_vals = np.clip(charge_vals, CHARGE_BIN_EDGES[0], CHARGE_BIN_EDGES[-1])
+
+        counts_by_mass.append(
+            np.histogramdd(
+                np.column_stack([mass_vals, binned_spin_phase]),
+                bins=[MASS_BIN_EDGES, np.arange(5)],
+            )[0]
+        )
+        counts_by_charge.append(
+            np.histogramdd(
+                np.column_stack([charge_vals, binned_spin_phase]),
+                bins=[CHARGE_BIN_EDGES, np.arange(5)],
+            )[0]
+        )
+        counts_by_mass_map.append(
+            np.histogramdd(
+                np.column_stack([mass_vals, longitude, latitude]),
+                bins=[MASS_BIN_EDGES, LON_BINS_EDGES, LAT_BINS_EDGES],
+            )[0]
+        )
+        counts_by_charge_map.append(
+            np.histogramdd(
+                np.column_stack([charge_vals, longitude, latitude]),
+                bins=[CHARGE_BIN_EDGES, LON_BINS_EDGES, LAT_BINS_EDGES],
+            )[0]
+        )
 
     return (
-        counts_by_charge,
-        counts_by_mass,
-        counts_by_charge_map,
-        counts_by_mass_map,
+        np.stack(counts_by_charge),
+        np.stack(counts_by_mass),
+        np.stack(counts_by_charge_map),
+        np.stack(counts_by_mass_map),
         daily_epoch,
     )
 
@@ -598,7 +607,7 @@ def get_science_acquisition_timestamps(
     epochs = evt_dataset["epoch"][sc_indices].data
     # Now the state change values and check if it is either a science
     # acquisition start or science acquisition stop event.
-    for v1, v2, epoch in zip(val1, val2, epochs):
+    for v1, v2, epoch in zip(val1, val2, epochs, strict=False):
         # An "acquire" start will have val1=ACQSETUP and val2=ACQ
         # An "acquire" stop will have val1=ACQ and val2=CHILL
         if (v1, v2) == (IDEXEvtAcquireCodes.ACQSETUP, IDEXEvtAcquireCodes.ACQ):
@@ -637,6 +646,12 @@ def get_science_acquisition_on_percentage(evt_dataset: xr.Dataset) -> dict:
     """
     # Get science acquisition start and stop times
     evt_logs, evt_time, evt_values = get_science_acquisition_timestamps(evt_dataset)
+    if len(evt_time) == 0:
+        logger.warning(
+            "No science acquisition events found in event dataset. Returning empty "
+            "uptime percentages. All rate variables will be set to -1."
+        )
+        return {}
     # Track total and 'on' durations per day
     daily_totals: collections.defaultdict = defaultdict(timedelta)
     daily_on: collections.defaultdict = defaultdict(timedelta)

imap_processing/idex/idex_utils.py

@@ -1,7 +1,5 @@
 """Contains helper functions to support IDEX processing."""
 
-from typing import Optional
-
 import xarray as xr
 
 from imap_processing.cdf.imap_cdf_manager import ImapCdfAttributes
@@ -31,7 +29,7 @@ def setup_dataset(
     dataset: xr.Dataset,
     match_strings: list,
     idex_attrs: ImapCdfAttributes,
-    data_vars: Optional[dict] = None,
+    data_vars: dict | None = None,
 ) -> xr.Dataset:
     """
     Initialize a dataset and copy over any dataArrays.

imap_processing/lo/ancillary_data/imap_lo_hydrogen-geometric-factor_v001.csv

@@ -0,0 +1,75 @@
+incident_E-Step,Observed_E-Step,Cntr_E,Cntr_E_unc,GF_Dbl_all,GF_Dbl_all_unc,GF_Trpl_all,GF_Trpl_all_unc,GF_Dbl_H,GF_Dbl_H_unc,GF_Trpl_H,GF_Trpl_H_unc
+Hi_Res,,[keV],[keV],[cm^2sr keV/keV],[cm^2sr keV/keV],[cm^2sr keV/keV],[cm^2sr keV/keV],[cm^2sr keV/keV],[cm^2sr keV/keV],[cm^2sr keV/keV],[cm^2sr keV/keV]
+1,1,0.015,0.00135,5.35E-05,4.82E-06,2.20E-05,1.87E-06,5.32E-05,4.26E-06,2.19E-05,1.86E-06
+2,1,0.015,0.00135,3.10E-05,2.79E-06,1.25E-05,1.06E-06,3.09E-05,2.47E-06,1.27E-05,1.08E-06
+2,2,0.029,0.00261,1.04E-04,9.36E-06,4.27E-05,3.63E-06,1.03E-04,8.26E-06,4.24E-05,3.61E-06
+3,1,0.015,0.00135,4.13E-05,3.72E-06,1.68E-05,1.43E-06,4.16E-05,3.33E-06,1.71E-05,1.45E-06
+3,2,0.029,0.00261,1.28E-04,1.15E-05,5.02E-05,4.26E-06,1.30E-04,1.04E-05,5.35E-05,4.55E-06
+3,3,0.055,0.00495,1.59E-04,1.43E-05,6.54E-05,5.56E-06,1.59E-04,1.27E-05,6.51E-05,5.54E-06
+4,1,0.015,0.00135,9.85E-06,8.86E-07,3.73E-06,3.17E-07,7.97E-06,6.38E-07,3.28E-06,2.79E-07
+4,2,0.029,0.00261,1.71E-05,1.54E-06,6.54E-06,5.56E-07,1.68E-05,1.34E-06,6.89E-06,5.86E-07
+4,3,0.055,0.00495,6.72E-05,6.05E-06,2.58E-05,2.19E-06,6.92E-05,5.53E-06,2.84E-05,2.42E-06
+4,4,0.11,0.0099,1.79E-04,1.61E-05,7.36E-05,6.25E-06,1.78E-04,1.42E-05,7.30E-05,6.20E-06
+5,1,0.015,0.00135,1.62E-05,1.46E-06,5.05E-06,4.29E-07,7.86E-06,6.29E-07,3.23E-06,2.75E-07
+5,2,0.029,0.00261,1.57E-05,1.41E-06,5.16E-06,4.39E-07,9.72E-06,7.77E-07,3.99E-06,3.39E-07
+5,3,0.055,0.00495,1.45E-05,1.31E-06,4.98E-06,4.24E-07,1.49E-05,1.20E-06,6.14E-06,5.22E-07
+5,4,0.11,0.0099,5.81E-05,5.23E-06,2.15E-05,1.83E-06,6.00E-05,4.80E-06,2.47E-05,2.10E-06
+5,5,0.209,0.01881,1.77E-04,1.59E-05,7.26E-05,6.17E-06,1.76E-04,1.41E-05,7.22E-05,6.14E-06
+6,1,0.015,0.00135,1.34E-05,1.21E-06,4.33E-06,3.68E-07,9.40E-06,7.52E-07,3.86E-06,3.28E-07
+6,2,0.029,0.00261,1.61E-05,1.45E-06,4.64E-06,3.94E-07,1.11E-05,8.92E-07,4.58E-06,3.89E-07
+6,3,0.055,0.00495,1.74E-05,1.57E-06,5.02E-06,4.27E-07,1.42E-05,1.14E-06,5.84E-06,4.96E-07
+6,4,0.11,0.0099,1.45E-05,1.30E-06,4.28E-06,3.64E-07,1.56E-05,1.25E-06,6.40E-06,5.44E-07
+6,5,0.209,0.01881,5.81E-05,5.23E-06,2.06E-05,1.75E-06,6.01E-05,4.81E-06,2.47E-05,2.10E-06
+6,6,0.439,0.03951,2.06E-04,1.85E-05,8.46E-05,7.19E-06,2.06E-04,1.64E-05,8.45E-05,7.18E-06
+7,1,0.015,0.00135,5.32E-06,4.79E-07,1.51E-06,1.29E-07,5.05E-06,4.04E-07,2.08E-06,1.76E-07
+7,2,0.029,0.00261,1.61E-05,1.45E-06,5.15E-06,4.38E-07,1.39E-05,1.11E-06,5.70E-06,4.85E-07
+7,3,0.055,0.00495,1.96E-05,1.77E-06,6.10E-06,5.19E-07,1.42E-05,1.14E-06,5.84E-06,4.96E-07
+7,4,0.11,0.0099,2.08E-05,1.87E-06,6.43E-06,5.46E-07,1.65E-05,1.32E-06,6.78E-06,5.77E-07
+7,5,0.209,0.01881,2.52E-05,2.27E-06,8.69E-06,7.39E-07,2.61E-05,2.09E-06,1.07E-05,9.11E-07
+7,6,0.439,0.03951,1.93E-04,1.74E-05,6.94E-05,5.90E-06,1.98E-04,1.59E-05,8.15E-05,6.93E-06
+7,7,0.872,0.07848,3.82E-04,3.44E-05,1.57E-04,1.33E-05,3.82E-04,3.05E-05,1.57E-04,1.33E-05
+8,1,0.015,0.00135,5.28E-05,4.75E-06,2.25E-05,1.91E-06,4.10E-05,3.28E-06,1.68E-05,1.43E-06
+8,2,0.029,0.00261,6.14E-06,5.53E-07,1.96E-06,1.66E-07,5.25E-06,4.20E-07,2.16E-06,1.83E-07
+8,3,0.055,0.00495,1.13E-05,1.02E-06,3.39E-06,2.88E-07,9.37E-06,7.50E-07,3.85E-06,3.27E-07
+8,4,0.11,0.0099,1.16E-05,1.04E-06,3.83E-06,3.25E-07,8.69E-06,6.95E-07,3.57E-06,3.03E-07
+8,5,0.209,0.01881,1.09E-05,9.78E-07,3.55E-06,3.02E-07,8.41E-06,6.73E-07,3.46E-06,2.94E-07
+8,6,0.439,0.03951,3.29E-05,2.96E-06,1.28E-05,1.08E-06,3.28E-05,2.62E-06,1.35E-05,1.15E-06
+8,7,0.872,0.07848,3.95E-04,3.56E-05,1.44E-04,1.23E-05,4.04E-04,3.23E-05,1.66E-04,1.41E-05
+8,8,1.821,0.16389,5.41E-04,4.87E-05,2.22E-04,1.89E-05,5.41E-04,4.33E-05,2.22E-04,1.89E-05
+Hi_Thr,,,,,,,,,,,
+1,1,0.016,0.00144,8.92E-05,8.03E-06,3.67E-05,3.12E-06,8.87E-05,7.10E-06,3.65E-05,3.10E-06
+2,1,0.016,0.00144,5.16E-05,4.65E-06,2.08E-05,1.76E-06,5.15E-05,4.12E-06,2.12E-05,1.80E-06
+2,2,0.03,0.0027,1.73E-04,1.56E-05,7.12E-05,6.06E-06,1.72E-04,1.38E-05,7.07E-05,6.01E-06
+3,1,0.016,0.00144,6.88E-05,6.19E-06,2.80E-05,2.38E-06,6.94E-05,5.55E-06,2.85E-05,2.42E-06
+3,2,0.03,0.0027,2.13E-04,1.91E-05,8.36E-05,7.11E-06,2.17E-04,1.74E-05,8.92E-05,7.58E-06
+3,3,0.056,0.00504,2.65E-04,2.39E-05,1.09E-04,9.26E-06,2.64E-04,2.11E-05,1.09E-04,9.23E-06
+4,1,0.016,0.00144,1.64E-05,1.48E-06,6.22E-06,5.29E-07,1.33E-05,1.06E-06,5.46E-06,4.64E-07
+4,2,0.03,0.0027,2.84E-05,2.56E-06,1.09E-05,9.27E-07,2.80E-05,2.24E-06,1.15E-05,9.77E-07
+4,3,0.056,0.00504,1.12E-04,1.01E-05,4.30E-05,3.66E-06,1.15E-04,9.22E-06,4.74E-05,4.03E-06
+4,4,0.111,0.00999,2.98E-04,2.68E-05,1.23E-04,1.04E-05,2.96E-04,2.37E-05,1.22E-04,1.03E-05
+5,1,0.016,0.00144,2.70E-05,2.43E-06,8.42E-06,7.15E-07,1.31E-05,1.05E-06,5.38E-06,4.58E-07
+5,2,0.03,0.0027,2.61E-05,2.35E-06,8.61E-06,7.31E-07,1.62E-05,1.30E-06,6.65E-06,5.66E-07
+5,3,0.056,0.00504,2.42E-05,2.18E-06,8.30E-06,7.06E-07,2.49E-05,1.99E-06,1.02E-05,8.70E-07
+5,4,0.111,0.00999,9.69E-05,8.72E-06,3.58E-05,3.04E-06,1.00E-04,8.00E-06,4.11E-05,3.49E-06
+5,5,0.21,0.0189,2.94E-04,2.65E-05,1.21E-04,1.03E-05,2.93E-04,2.34E-05,1.20E-04,1.02E-05
+6,1,0.016,0.00144,2.24E-05,2.01E-06,7.22E-06,6.14E-07,1.57E-05,1.25E-06,6.44E-06,5.47E-07
+6,2,0.03,0.0027,2.68E-05,2.41E-06,7.73E-06,6.57E-07,1.86E-05,1.49E-06,7.64E-06,6.49E-07
+6,3,0.056,0.00504,2.91E-05,2.61E-06,8.37E-06,7.12E-07,2.37E-05,1.89E-06,9.73E-06,8.27E-07
+6,4,0.111,0.00999,2.41E-05,2.17E-06,7.13E-06,6.06E-07,2.60E-05,2.08E-06,1.07E-05,9.07E-07
+6,5,0.21,0.0189,9.68E-05,8.71E-06,3.44E-05,2.92E-06,1.00E-04,8.02E-06,4.12E-05,3.50E-06
+6,6,0.44,0.0396,3.43E-04,3.09E-05,1.41E-04,1.20E-05,3.43E-04,2.74E-05,1.41E-04,1.20E-05
+7,1,0.016,0.00144,8.87E-06,7.98E-07,2.52E-06,2.14E-07,8.42E-06,6.73E-07,3.46E-06,2.94E-07
+7,2,0.03,0.0027,2.68E-05,2.41E-06,8.59E-06,7.30E-07,2.31E-05,1.85E-06,9.50E-06,8.08E-07
+7,3,0.056,0.00504,3.27E-05,2.95E-06,1.02E-05,8.64E-07,2.37E-05,1.89E-06,9.73E-06,8.27E-07
+7,4,0.111,0.00999,3.47E-05,3.12E-06,1.07E-05,9.10E-07,2.75E-05,2.20E-06,1.13E-05,9.61E-07
+7,5,0.21,0.0189,4.20E-05,3.78E-06,1.45E-05,1.23E-06,4.35E-05,3.48E-06,1.79E-05,1.52E-06
+7,6,0.44,0.0396,3.21E-04,2.89E-05,1.16E-04,9.84E-06,3.31E-04,2.65E-05,1.36E-04,1.16E-05
+7,7,0.873,0.07857,6.36E-04,5.73E-05,2.61E-04,2.22E-05,6.36E-04,5.09E-05,2.61E-04,2.22E-05
+8,1,0.016,0.00144,8.79E-05,7.91E-06,3.75E-05,3.19E-06,6.83E-05,5.46E-06,2.81E-05,2.39E-06
+8,2,0.03,0.0027,1.02E-05,9.21E-07,3.26E-06,2.77E-07,8.75E-06,7.00E-07,3.60E-06,3.06E-07
+8,3,0.056,0.00504,1.88E-05,1.70E-06,5.65E-06,4.80E-07,1.56E-05,1.25E-06,6.42E-06,5.46E-07
+8,4,0.111,0.00999,1.93E-05,1.73E-06,6.38E-06,5.42E-07,1.45E-05,1.16E-06,5.95E-06,5.06E-07
+8,5,0.21,0.0189,1.81E-05,1.63E-06,5.92E-06,5.03E-07,1.40E-05,1.12E-06,5.76E-06,4.90E-07
+8,6,0.44,0.0396,5.48E-05,4.93E-06,2.13E-05,1.81E-06,5.46E-05,4.37E-06,2.25E-05,1.91E-06
+8,7,0.873,0.07857,6.59E-04,5.93E-05,2.40E-04,2.04E-05,6.74E-04,5.39E-05,2.77E-04,2.35E-05
+8,8,1.822,0.16398,9.01E-04,8.11E-05,3.70E-04,3.15E-05,9.01E-04,7.21E-05,3.70E-04,3.15E-05

imap_processing/lo/ancillary_data/imap_lo_oxygen-geometric-factor_v001.csv

@@ -0,0 +1,75 @@
+incident_E-Step,Observed_E-Step,Cntr_E,Cntr_E_unc,GF_Dbl_all,GF_Dbl_all_unc,GF_Trpl_all,GF_Trpl_all_unc,GF_Dbl_O,GF_Dbl_O_unc,GF_Trpl_O,GF_Trpl_O_unc
+Hi_Res,,[keV],[keV],[cm^2sr keV/keV],[cm^2sr keV/keV],[cm^2sr keV/keV],[cm^2sr keV/keV],[cm^2sr keV/keV],[cm^2sr keV/keV],[cm^2sr keV/keV],[cm^2sr keV/keV]
+1,1,0.016,0.00144,3.8700E-04,3.4400E-04,8.9100E-05,8.0200E-05,4.2600E-04,3.4100E-04,1.0600E-04,9.1500E-05
+2,1,0.016,0.00144,2.8700E-04,2.5600E-04,6.6100E-05,5.9500E-05,2.7700E-04,2.2100E-04,6.9200E-05,5.9500E-05
+2,2,0.032,0.00288,4.4800E-04,3.9900E-04,9.6700E-05,8.7000E-05,4.7800E-04,3.8200E-04,1.2000E-04,1.0300E-04
+3,1,0.016,0.00144,5.2900E-04,4.7100E-04,1.3200E-04,1.1900E-04,4.1800E-04,3.3400E-04,1.0400E-04,8.9800E-05
+3,2,0.032,0.00288,5.9300E-04,5.2800E-04,1.4100E-04,1.2700E-04,6.3000E-04,5.0400E-04,1.5800E-04,1.3600E-04
+3,3,0.065,0.00585,6.9900E-04,6.2200E-04,1.6600E-04,1.5000E-04,7.6400E-04,6.1100E-04,1.9100E-04,1.6400E-04
+4,1,0.016,0.00144,6.6400E-04,5.9100E-04,1.7700E-04,1.5900E-04,4.1600E-04,3.3300E-04,1.0400E-04,8.9400E-05
+4,2,0.032,0.00288,1.1500E-03,1.0200E-03,2.6900E-04,2.4200E-04,1.0900E-03,8.7300E-04,2.7300E-04,2.3500E-04
+4,3,0.065,0.00585,2.3600E-03,2.1000E-03,5.1700E-04,4.6500E-04,2.5200E-03,2.0100E-03,6.2900E-04,5.4100E-04
+4,4,0.135,0.01215,8.1500E-04,7.2600E-04,1.9100E-04,1.7200E-04,8.7400E-04,7.0000E-04,2.1900E-04,1.8800E-04
+5,1,0.016,0.00144,2.8100E-04,2.5000E-04,8.8400E-05,7.9600E-05,1.2500E-04,9.9900E-05,3.1200E-05,2.6800E-05
+5,2,0.032,0.00288,3.8100E-04,3.4000E-04,1.0800E-04,9.7300E-05,2.0800E-04,1.6600E-04,5.1900E-05,4.4600E-05
+5,3,0.065,0.00585,5.8000E-04,5.1600E-04,1.3800E-04,1.2400E-04,5.7100E-04,4.5700E-04,1.4300E-04,1.2300E-04
+5,4,0.135,0.01215,5.4100E-04,4.8100E-04,1.2100E-04,1.0800E-04,5.6900E-04,4.5500E-04,1.4200E-04,1.2200E-04
+5,5,0.279,0.02511,8.4100E-04,7.4900E-04,2.0500E-04,1.8500E-04,8.6900E-04,6.9500E-04,2.1700E-04,1.8700E-04
+6,1,0.016,0.00144,2.2100E-04,1.9700E-04,7.1000E-05,6.3900E-05,4.7600E-05,3.8100E-05,1.1900E-05,1.0200E-05
+6,2,0.032,0.00288,3.1900E-04,2.8400E-04,9.7800E-05,8.8000E-05,1.3800E-04,1.1100E-04,3.4600E-05,2.9800E-05
+6,3,0.065,0.00585,4.4900E-04,4.0000E-04,1.2400E-04,1.1100E-04,1.7000E-04,1.3600E-04,4.2500E-05,3.6500E-05
+6,4,0.135,0.01215,6.7700E-04,6.0300E-04,1.5400E-04,1.3800E-04,5.1000E-04,4.0800E-04,1.2700E-04,1.1000E-04
+6,5,0.279,0.02511,1.2800E-03,1.1400E-03,2.9100E-04,2.6200E-04,1.1700E-03,9.3300E-04,2.9100E-04,2.5100E-04
+6,6,0.601,0.05409,8.3600E-04,7.4400E-04,1.9800E-04,1.7800E-04,8.2700E-04,6.6100E-04,2.0700E-04,1.7800E-04
+7,1,0.016,0.00144,1.8500E-04,1.6500E-04,5.9600E-05,5.3700E-05,5.6700E-05,4.5400E-05,1.4200E-05,1.2200E-05
+7,2,0.032,0.00288,2.3400E-04,2.0800E-04,7.3700E-05,6.6300E-05,9.2100E-05,7.3700E-05,2.3000E-05,1.9800E-05
+7,3,0.065,0.00585,3.2000E-04,2.8500E-04,9.2000E-05,8.2800E-05,1.0600E-04,8.5000E-05,2.6600E-05,2.2800E-05
+7,4,0.135,0.01215,4.4900E-04,4.0000E-04,1.2000E-04,1.0800E-04,1.7800E-04,1.4200E-04,4.4400E-05,3.8200E-05
+7,5,0.279,0.02511,7.1800E-04,6.3900E-04,1.7600E-04,1.5800E-04,6.1600E-04,4.9300E-04,1.5400E-04,1.3200E-04
+7,6,0.601,0.05409,1.0300E-03,9.1800E-04,2.4500E-04,2.2000E-04,9.9900E-04,7.9900E-04,2.5000E-04,2.1500E-04
+7,7,1.206,0.10854,9.6000E-04,8.5400E-04,2.4400E-04,2.1900E-04,9.6100E-04,7.6900E-04,2.4000E-04,2.0700E-04
+8,1,0.016,0.00144,2.8900E-04,2.5800E-04,9.1800E-05,8.2700E-05,4.4800E-05,3.5800E-05,1.1200E-05,9.6300E-06
+8,2,0.032,0.00288,3.4200E-04,3.0500E-04,1.0500E-04,9.4300E-05,6.8300E-05,5.4600E-05,1.7100E-05,1.4700E-05
+8,3,0.065,0.00585,1.9300E-04,1.7200E-04,5.6900E-05,5.1200E-05,4.9700E-05,3.9800E-05,1.2400E-05,1.0700E-05
+8,4,0.135,0.01215,2.7000E-04,2.4000E-04,7.2300E-05,6.5000E-05,6.1300E-05,4.9000E-05,1.5300E-05,1.3200E-05
+8,5,0.279,0.02511,3.7600E-04,3.3500E-04,9.7400E-05,8.7700E-05,1.5500E-04,1.2400E-04,3.8800E-05,3.3300E-05
+8,6,0.601,0.05409,7.3600E-04,6.5500E-04,1.8400E-04,1.6600E-04,6.0000E-04,4.8000E-04,1.5000E-04,1.2900E-04
+8,7,1.206,0.10854,1.0600E-03,9.4100E-04,2.5600E-04,2.3100E-04,9.7300E-04,7.7900E-04,2.4300E-04,2.0900E-04
+8,8,2.361,0.21249,9.8400E-04,8.7500E-04,2.5500E-04,2.2900E-04,9.3600E-04,7.4900E-04,2.3400E-04,2.0100E-04
+Hi_Thr,,,,,,,,,,,
+1,1,0.017,0.00153,6.4500E-04,5.7400E-04,1.4900E-04,1.3400E-04,7.1000E-04,5.6800E-04,1.7700E-04,1.5300E-04
+2,1,0.017,0.00153,4.7900E-04,4.2600E-04,1.1000E-04,9.9200E-05,4.6100E-04,3.6900E-04,1.1500E-04,9.9200E-05
+2,2,0.033,0.00297,7.4600E-04,6.6400E-04,1.6100E-04,1.4500E-04,7.9700E-04,6.3700E-04,1.9900E-04,1.7100E-04
+3,1,0.017,0.00153,8.8200E-04,7.8500E-04,2.2000E-04,1.9800E-04,6.9600E-04,5.5700E-04,1.7400E-04,1.5000E-04
+3,2,0.033,0.00297,9.8800E-04,8.7900E-04,2.3500E-04,2.1100E-04,1.0500E-03,8.4100E-04,2.6300E-04,2.2600E-04
+3,3,0.066,0.00594,1.1700E-03,1.0400E-03,2.7700E-04,2.4900E-04,1.2700E-03,1.0200E-03,3.1800E-04,2.7400E-04
+4,1,0.017,0.00153,1.1100E-03,9.8600E-04,2.9500E-04,2.6600E-04,6.9300E-04,5.5400E-04,1.7300E-04,1.4900E-04
+4,2,0.033,0.00297,1.9200E-03,1.7100E-03,4.4800E-04,4.0300E-04,1.8200E-03,1.4600E-03,4.5500E-04,3.9100E-04
+4,3,0.066,0.00594,3.9300E-03,3.4900E-03,8.6100E-04,7.7500E-04,4.1900E-03,3.3500E-03,1.0500E-03,9.0100E-04
+4,4,0.136,0.01224,1.3600E-03,1.2100E-03,3.1800E-04,2.8600E-04,1.4600E-03,1.1700E-03,3.6400E-04,3.1300E-04
+5,1,0.017,0.00153,4.6800E-04,4.1700E-04,1.4700E-04,1.3300E-04,2.0800E-04,1.6700E-04,5.2000E-05,4.4700E-05
+5,2,0.033,0.00297,6.3600E-04,5.6600E-04,1.8000E-04,1.6200E-04,3.4600E-04,2.7700E-04,8.6500E-05,7.4400E-05
+5,3,0.066,0.00594,9.6600E-04,8.6000E-04,2.3000E-04,2.0700E-04,9.5100E-04,7.6100E-04,2.3800E-04,2.0400E-04
+5,4,0.136,0.01224,9.0100E-04,8.0200E-04,2.0100E-04,1.8100E-04,9.4800E-04,7.5800E-04,2.3700E-04,2.0400E-04
+5,5,0.28,0.0252,1.4000E-03,1.2500E-03,3.4200E-04,3.0800E-04,1.4500E-03,1.1600E-03,3.6200E-04,3.1100E-04
+6,1,0.017,0.00153,3.6800E-04,3.2800E-04,1.1800E-04,1.0700E-04,7.9400E-05,6.3500E-05,1.9800E-05,1.7100E-05
+6,2,0.033,0.00297,5.3200E-04,4.7300E-04,1.6300E-04,1.4700E-04,2.3100E-04,1.8500E-04,5.7700E-05,4.9600E-05
+6,3,0.066,0.00594,7.4800E-04,6.6600E-04,2.0600E-04,1.8500E-04,2.8300E-04,2.2600E-04,7.0800E-05,6.0900E-05
+6,4,0.136,0.01224,1.1300E-03,1.0000E-03,2.5600E-04,2.3000E-04,8.4900E-04,6.8000E-04,2.1200E-04,1.8300E-04
+6,5,0.28,0.0252,2.1400E-03,1.9100E-03,4.8500E-04,4.3700E-04,1.9400E-03,1.5500E-03,4.8600E-04,4.1800E-04
+6,6,0.602,0.05418,1.3900E-03,1.2400E-03,3.3000E-04,2.9700E-04,1.3800E-03,1.1000E-03,3.4400E-04,2.9600E-04
+7,1,0.017,0.00153,3.0900E-04,2.7500E-04,9.9400E-05,8.9500E-05,9.4500E-05,7.5600E-05,2.3600E-05,2.0300E-05
+7,2,0.033,0.00297,3.9000E-04,3.4700E-04,1.2300E-04,1.1100E-04,1.5400E-04,1.2300E-04,3.8400E-05,3.3000E-05
+7,3,0.066,0.00594,5.3400E-04,4.7500E-04,1.5300E-04,1.3800E-04,1.7700E-04,1.4200E-04,4.4300E-05,3.8100E-05
+7,4,0.136,0.01224,7.4900E-04,6.6700E-04,2.0000E-04,1.8000E-04,2.9600E-04,2.3700E-04,7.4000E-05,6.3700E-05
+7,5,0.28,0.0252,1.2000E-03,1.0700E-03,2.9300E-04,2.6400E-04,1.0300E-03,8.2100E-04,2.5700E-04,2.2100E-04
+7,6,0.602,0.05418,1.7200E-03,1.5300E-03,4.0800E-04,3.6700E-04,1.6700E-03,1.3300E-03,4.1600E-04,3.5800E-04
+7,7,1.207,0.10863,1.6000E-03,1.4200E-03,4.0600E-04,3.6500E-04,1.6000E-03,1.2800E-03,4.0000E-04,3.4400E-04
+8,1,0.017,0.00153,4.8200E-04,4.2900E-04,1.5300E-04,1.3800E-04,7.4700E-05,5.9700E-05,1.8700E-05,1.6100E-05
+8,2,0.033,0.00297,5.7000E-04,5.0800E-04,1.7500E-04,1.5700E-04,1.1400E-04,9.1100E-05,2.8500E-05,2.4500E-05
+8,3,0.066,0.00594,3.2100E-04,2.8600E-04,9.4800E-05,8.5300E-05,8.2900E-05,6.6300E-05,2.0700E-05,1.7800E-05
+8,4,0.136,0.01224,4.5000E-04,4.0100E-04,1.2000E-04,1.0800E-04,1.0200E-04,8.1700E-05,2.5500E-05,2.2000E-05
+8,5,0.28,0.0252,6.2700E-04,5.5800E-04,1.6200E-04,1.4600E-04,2.5800E-04,2.0700E-04,6.4600E-05,5.5600E-05
+8,6,0.602,0.05418,1.2300E-03,1.0900E-03,3.0700E-04,2.7600E-04,1.0000E-03,8.0000E-04,2.5000E-04,2.1500E-04
+8,7,1.207,0.10863,1.7600E-03,1.5700E-03,4.2700E-04,3.8400E-04,1.6200E-03,1.3000E-03,4.0600E-04,3.4900E-04
+8,8,2.362,0.21258,1.6400E-03,1.4600E-03,4.2500E-04,3.8200E-04,1.5600E-03,1.2500E-03,3.9000E-04,3.3500E-04