aboutsummaryrefslogtreecommitdiffstats
path: root/fg21sim
diff options
context:
space:
mode:
Diffstat (limited to 'fg21sim')
-rw-r--r--fg21sim/configs/config.spec8
-rw-r--r--fg21sim/extragalactic/clusters/halo.py38
2 files changed, 0 insertions, 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.