Merge pull request #46 from Joshiwavm/fix_notebooks_20_12

Fix notebooks 20 12
thomaswmorris · Dec 24, 2023 · 13bdc25 · 13bdc25
2 parents 49d889e + 1678717
commit 13bdc25
Show file tree

Hide file tree

Showing 3 changed files with 78 additions and 61 deletions.
diff --git a/docs/source/tutorials/Synthetic_AtLAST_observations_of_galaxies.ipynb b/docs/source/tutorials/Synthetic_AtLAST_observations_of_galaxies.ipynb
diff --git a/docs/source/tutorials/Synthetic_MUSTANG2_observations.ipynb b/docs/source/tutorials/Synthetic_MUSTANG2_observations.ipynb
diff --git a/maria/map/mappers.py b/maria/map/mappers.py
@@ -41,6 +41,7 @@ def __init__(
         filter_tods: bool = True,
         smoothing: float = 8,
         degrees: bool = True,
+        ffilter: float = 0.08,
         **kwargs,
     ):
         self.res = np.radians(res) if degrees else res
@@ -51,6 +52,7 @@ def __init__(
         self.frame = frame
         self.filter_tods = filter_tods
         self.smoothing = smoothing
+        self.ffilter = ffilter
 
         self.n_x = int(np.maximum(1, self.width / self.res))
         self.n_y = int(np.maximum(1, self.height / self.res))
@@ -132,7 +134,6 @@ def save_maps(self, filepath):
             self.header["CRVAL3"] = self.band_data[key]["nom_freq"] * 1e9
             self.header["CDELT3"] = self.band_data[key]["nom_freqwidth"] * 1e9
 
-            # save_maps[i] = self.maps[list(self.maps.keys())[i]]
             sigma_smooth = self.smoothing / 3600 / np.rad2deg(self.res) / 2.355
             save_maps[i] = self.smoothed_maps(sigma_smooth)[
                 list(self.band_data.keys())[i]
@@ -158,7 +159,7 @@ def __init__(self, **kwargs):
         self._nmtr = kwargs.get("n_modes_to_remove", 0)
 
     def _fourier_filter(self, tod_dat, tod_time):
-        ffilt = [0.08, 51.0]  # high-pass and low-pass filters, in Hz
+        ffilt = [self.ffilter, 51.0]  # high-pass and low-pass filters, in Hz
         width = 0.05
 
         n = len(tod_time)