From cc71b0bc70b5602dcdc9bdeffdfdc3632f92029f Mon Sep 17 00:00:00 2001 From: Aaron LI Date: Mon, 14 Aug 2017 20:49:02 +0800 Subject: clusters/halo.py: Simplify two coefficients usage Signed-off-by: Aaron LI --- fg21sim/extragalactic/clusters/halo.py | 12 +++++------- 1 file changed, 5 insertions(+), 7 deletions(-) (limited to 'fg21sim/extragalactic') diff --git a/fg21sim/extragalactic/clusters/halo.py b/fg21sim/extragalactic/clusters/halo.py index b4a504d..d5c5cf1 100644 --- a/fg21sim/extragalactic/clusters/halo.py +++ b/fg21sim/extragalactic/clusters/halo.py @@ -558,9 +558,8 @@ class RadioHalo: """ z = COSMO.redshift(t) mass = self._mass(t) - n_th = helper.density_number_thermal(mass, z) - coef = -1.20e-12 * AUC.Gyr2s # [Gyr^-1] - loss = coef * n_th * (1 + np.log(gamma/n_th) / 75) + n_th = helper.density_number_thermal(mass, z) # [cm^-3] + loss = -3.79e4 * n_th * (1 + np.log(gamma/n_th) / 75) return loss def _loss_rad(self, gamma, t): @@ -573,8 +572,7 @@ class RadioHalo: Ref.[sarazin1999],Eq.(6,7) """ z = COSMO.redshift(t) - mass = self._mass(t) - B = helper.magnetic_field(mass) - coef = -1.37e-20 * AUC.Gyr2s # [Gyr^-1] - loss = coef * gamma**2 * ((B/3.25)**2 + (1+z)**4) + mass = self._mass(t) # [Msun] + B = helper.magnetic_field(mass) # [uG] + loss = -4.32e-4 * gamma**2 * ((B/3.25)**2 + (1+z)**4) return loss -- cgit v1.2.2