clusters/halo: Rewrite is_geniune() method

1 files changed, 30 insertions, 0 deletions

diff --git a/fg21sim/extragalactic/clusters/halo.py b/fg21sim/extragalactic/clusters/halo.py
index da3079b..9c25a3a 100644
--- a/fg21sim/extragalactic/clusters/halo.py
+++ b/fg21sim/extragalactic/clusters/halo.py
@@ -454,6 +454,36 @@ class RadioHalo:
 
         return n_e
 
+    def is_genuine(self, n_e):
+        """
+        Check whether the radio halo is genuine/observable by comparing the
+        emissivity to the corresponding fiducial value, which is calculated
+        from the fiducial electron spectrum derived with turbulent
+        acceleration turned off.
+
+        Parameters
+        ----------
+        n_e : float 1D `~numpy.ndarray`
+            The finally derived electron spectrum.
+            Unit: [cm^-3]
+
+        Returns
+        -------
+        genuine : bool
+            Whether the radio halo is genuine?
+        factor : float
+            Acceleration factor of the flux.
+        """
+        haloem = HaloEmission(gamma=self.gamma, n_e=n_e, B=1)
+        em = haloem.calc_emissivity(self.fiducial_freq)
+
+        ne_fiducial = self.calc_electron_spectrum(fiducial=True)
+        haloem.n_e = ne_fiducial
+        em_fiducial = haloem.calc_emissivity(self.fiducial_freq)
+
+        factor = em / em_fiducial
+        return (factor >= self.fiducial_factor, factor)
+
     def fp_injection(self, gamma, t=None):
         """
         Electron injection (rate) term for the Fokker-Planck equation.