diff --git a/imap_processing/lo/l2/lo_l2.py b/imap_processing/lo/l2/lo_l2.py index 8785d5485..50a91ed9a 100644 --- a/imap_processing/lo/l2/lo_l2.py +++ b/imap_processing/lo/l2/lo_l2.py @@ -94,6 +94,10 @@ def lo_l2( logger.info("Step 3: Converting to dataset and adding geometric factors") dataset = sky_map.to_dataset() + + if cg_correction: + dataset["energy_sc"] /= dataset["exposure_factor"] + dataset = add_geometric_factors(dataset, map_descriptor.species) logger.info("Step 4: Calculating rates and intensities") @@ -394,6 +398,8 @@ def process_single_pset( pset_processed, energy_values_ev ) + pset_processed["energy_sc"] *= pset_processed["exposure_factor"] + # Always calculate ram-mask to identify ram/anti-ram bins pset_processed = calculate_ram_mask(pset_processed) @@ -869,8 +875,8 @@ def calculate_all_rates_and_intensities( logger.info("Interpolating map intensities to helio-frame energies") dataset = interpolate_map_flux_to_helio_frame( dataset, - dataset["energy"], - dataset["energy"], + dataset["energy"] * 1000, + dataset["energy"] * 1000, ["ena_intensity", "bg_intensity"], )