From 39ac9e6c8bde34f78e7e520b04d7f2eae90e7795 Mon Sep 17 00:00:00 2001
From: Aaron LI <aly@aaronly.me>
Date: Wed, 27 Feb 2019 20:43:58 +0800
Subject: clusters/halo: calc_acc_factor() no longer used

---
 fg21sim/configs/config.spec            |  8 -------
 fg21sim/extragalactic/clusters/halo.py | 38 ----------------------------------
 2 files changed, 46 deletions(-)

diff --git a/fg21sim/configs/config.spec b/fg21sim/configs/config.spec
index 701bc91..fd45949 100644
--- a/fg21sim/configs/config.spec
+++ b/fg21sim/configs/config.spec
@@ -456,11 +456,3 @@ stream = option("stderr", "stdout", "", default="stderr")
   f_rc = float(default=0.1)
   # The slope parameter (i.e., beta).
   beta = float(default=0.6667)
-
-  # The minimum and maximum frequencies as well as the number of frequency
-  # points for calculating the total radio emissivity, which is used to
-  # approximate the bolometric emissivity.
-  # Unit: [MHz]
-  bolo_freq_min = float(default=10)
-  bolo_freq_max = float(default=1e5)
-  bolo_freq_num = integer(default=256)
diff --git a/fg21sim/extragalactic/clusters/halo.py b/fg21sim/extragalactic/clusters/halo.py
index fbf1cb2..4835d79 100644
--- a/fg21sim/extragalactic/clusters/halo.py
+++ b/fg21sim/extragalactic/clusters/halo.py
@@ -178,9 +178,6 @@ class RadioHalo1M:
         self.injection_index = configs.getn(sec+"/injection_index")
         self.f_rc = configs.getn(sec+"/f_rc")
         self.beta = configs.getn(sec+"/beta")
-        self.bolo_freq_min = configs.getn(sec+"/bolo_freq_min")
-        self.bolo_freq_max = configs.getn(sec+"/bolo_freq_max")
-        self.bolo_freq_num = configs.getn(sec+"/bolo_freq_num")
 
         self.kT_out = configs.getn("extragalactic/clusters/kT_out")
 
@@ -577,41 +574,6 @@ class RadioHalo1M:
 
         return n_e
 
-    def calc_acc_factor(self, n_e, n_e_fiducial=None):
-        """
-        Calculate the turbulence acceleration factor, which is estimated
-        as the ratio of the bolometric emissivity between the accelerated
-        electron spectrum and the fiducial electron spectrum derived with
-        turbulent acceleration turned off.
-
-        Parameters
-        ----------
-        n_e : float 1D `~numpy.ndarray`
-            The derived (accelerated) electron spectrum.
-            Unit: [cm^-3]
-        n_e_fiducial : float 1D `~numpy.ndarray`, optional
-            The fiducial electron spectrum.
-            Unit: [cm^-3]
-
-        Returns
-        -------
-        factor : float
-            Acceleration factor of the bolometric emissivity.
-        """
-        haloem = HaloEmission(gamma=self.gamma, n_e=n_e, B=1)
-        em = haloem.calc_emissivity_bolo(freq_min=self.bolo_freq_min,
-                                         freq_max=self.bolo_freq_max,
-                                         freq_num=self.bolo_freq_num)
-
-        if n_e_fiducial is None:
-            n_e_fiducial = self.calc_electron_spectrum(fiducial=True)
-        haloem.n_e = n_e_fiducial
-        em_fiducial = haloem.calc_emissivity_bolo(
-                freq_min=self.bolo_freq_min, freq_max=self.bolo_freq_max,
-                freq_num=self.bolo_freq_num)
-
-        return em / em_fiducial
-
     def fp_injection(self, gamma, t=None):
         """
         Electron injection (rate) term for the Fokker-Planck equation.
-- 
cgit v1.2.2