diff options
Diffstat (limited to 'fg21sim')
| -rw-r--r-- | fg21sim/configs/config.spec | 8 | ||||
| -rw-r--r-- | fg21sim/extragalactic/clusters/halo.py | 38 |
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. |