separated noise and beam effect

Author: samueldmcdermott

deepcmbsim/noise.py

@@ -26,6 +26,10 @@ def max_multipole(beamfwhm_arcmin, additional_factor=3):
     """
     return 180 * 60 * additional_factor / beamfwhm_arcmin
 
+def beam_effect(beamfwhm_arcmin, lmax):
+    ells = np.arange(lmax + 1)
+    return np.exp(ells * (ells + 1) * (beamfwhm_arcmin * arcmin_to_rad)**2 / (8 * np.log(2)))
+
 def detector_white_noise(noise_uKarcmin, beamfwhm_arcmin, lmax, TT=True, units_uK = True):
     """
     describes white (no angular scale) noise from a detector
@@ -54,4 +58,4 @@ def detector_white_noise(noise_uKarcmin, beamfwhm_arcmin, lmax, TT=True, units_u
     noise_uKarcmin = noise_uKarcmin if TT else noise_uKarcmin * np.sqrt(2)
     noise_uKarcmin = noise_uKarcmin if units_uK else noise_uKarcmin / 2.72548e6
     ells = np.arange(lmax+1)
-    return (noise_uKarcmin * arcmin_to_rad) ** 2 * np.exp(ells * (ells + 1) * (beamfwhm_arcmin * arcmin_to_rad) ** 2 / (8 * np.log(2)))
+    return (noise_uKarcmin * arcmin_to_rad) ** 2