From 500130085af123c2e28ef00abf0d5bcbd64c6cd2 Mon Sep 17 00:00:00 2001 From: Alfredo Mejia-Narvaez Date: Mon, 11 Nov 2024 14:27:57 -0300 Subject: [PATCH] avoiding NaNs in SCISEN* header keywords --- python/lvmdrp/functions/fluxCalMethod.py | 10 +++++++--- 1 file changed, 7 insertions(+), 3 deletions(-) diff --git a/python/lvmdrp/functions/fluxCalMethod.py b/python/lvmdrp/functions/fluxCalMethod.py index 363775cf..0cdaf218 100644 --- a/python/lvmdrp/functions/fluxCalMethod.py +++ b/python/lvmdrp/functions/fluxCalMethod.py @@ -571,9 +571,13 @@ def fluxcal_sci_ifu_stars(in_rss, plot=True, GAIA_CACHE_DIR=None, NSCI_MAX=15): label = rss._header['CCD'] channel = label.lower() - rss.setHdrValue(f"SCISENM{label}", np.nanmean(mean_sci[1000:3000]), f"Mean scistar sensitivity in {channel}") - rss.setHdrValue(f"SCISENR{label}", np.nanmean(rms_sci[1000:3000]), f"Mean scistar sensitivity rms in {channel}") - log.info(f"Mean scistar sensitivity in {channel} : {np.nanmean(mean_sci[1000:3000])}") + mean_sci_band = np.nanmean(mean_sci[1000:3000]) + rms_sci_band = np.nanmean(rms_sci[1000:3000]) + mean_sci_band = -999.9 if np.isnan(mean_sci_band) else mean_sci_band + rms_sci_band = -999.9 if np.isnan(rms_sci_band) else rms_sci_band + rss.setHdrValue(f"SCISENM{label}", mean_sci_band, f"Mean scistar sensitivity in {channel}") + rss.setHdrValue(f"SCISENR{label}", rms_sci_band, f"Mean scistar sensitivity rms in {channel}") + log.info(f"Mean scistar sensitivity in {channel} : {mean_sci_band}") if plot: plt.ylabel("sensitivity [(ergs/s/cm^2/A) / (e-/s/A)]")