From 54434d4b24104457d04af79a1a75ecf713c07bc1 Mon Sep 17 00:00:00 2001 From: Aaron LI Date: Tue, 29 Jan 2019 18:58:10 +0800 Subject: clusters/halo: Add options bolo_freq_{min,max,num} Update the default bolometric frequency band to: 10 MHz -- 100 GHz. --- fg21sim/extragalactic/clusters/halo.py | 11 +++++++++-- 1 file changed, 9 insertions(+), 2 deletions(-) (limited to 'fg21sim/extragalactic/clusters/halo.py') diff --git a/fg21sim/extragalactic/clusters/halo.py b/fg21sim/extragalactic/clusters/halo.py index e3521fe..2078ad1 100644 --- a/fg21sim/extragalactic/clusters/halo.py +++ b/fg21sim/extragalactic/clusters/halo.py @@ -166,6 +166,9 @@ class RadioHalo1M: self.injection_index = configs.getn(comp+"/injection_index") self.f_rc = configs.getn(comp+"/f_rc") self.beta = configs.getn(comp+"/beta") + self.bolo_freq_min = configs.getn(comp+"/bolo_freq_min") + self.bolo_freq_max = configs.getn(comp+"/bolo_freq_max") + self.bolo_freq_num = configs.getn(comp+"/bolo_freq_num") def _set_solver(self): self.fpsolver = FokkerPlanckSolver( @@ -510,12 +513,16 @@ class RadioHalo1M: Acceleration factor of the bolometric emissivity. """ haloem = HaloEmission(gamma=self.gamma, n_e=n_e, B=1) - em = haloem.calc_emissivity_bolo() + 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() + 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 -- cgit v1.2.2