diff options
| author | Aaron LI <aly@aaronly.me> | 2018-12-22 23:34:39 +0800 |
|---|---|---|
| committer | Aaron LI <aly@aaronly.me> | 2018-12-22 23:34:39 +0800 |
| commit | 1947fffba78c7b287594851c99cbd0dff4d22bfa (patch) | |
| tree | 217b19710f920e5ab979621c06f9e96c1928f478 | |
| parent | 1703eb5e04f9e35a56e5a764a8b08be1350d5326 (diff) | |
| download | fg21sim-1947fffba78c7b287594851c99cbd0dff4d22bfa.tar.bz2 | |
clusters/halo: Update is_genuine() to return flux factor also
| -rw-r--r-- | fg21sim/extragalactic/clusters/halo.py | 6 |
1 files changed, 5 insertions, 1 deletions
diff --git a/fg21sim/extragalactic/clusters/halo.py b/fg21sim/extragalactic/clusters/halo.py index 5580c14..e1e5e2e 100644 --- a/fg21sim/extragalactic/clusters/halo.py +++ b/fg21sim/extragalactic/clusters/halo.py @@ -520,6 +520,9 @@ class RadioHalo: 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=self.B_obs, radius=self.radius, @@ -530,7 +533,8 @@ class RadioHalo: haloem.n_e = ne_fiducial flux_fiducial = haloem.calc_flux(self.fiducial_freq) - return flux >= flux_fiducial * self.fiducial_factor + factor = flux / flux_fiducial + return (flux >= self.fiducial_factor, factor) def fp_injection(self, gamma, t=None): """ |