diff options
| author | Aaron LI <aly@aaronly.me> | 2017-08-14 20:49:02 +0800 |
|---|---|---|
| committer | Aaron LI <aly@aaronly.me> | 2017-08-14 20:49:02 +0800 |
| commit | cc71b0bc70b5602dcdc9bdeffdfdc3632f92029f (patch) | |
| tree | 1d6f754618fae8ad9c8e6341594ae47180619f60 /fg21sim | |
| parent | 1127396c18feed61c13e6c76bef48fca1123c33b (diff) | |
| download | fg21sim-cc71b0bc70b5602dcdc9bdeffdfdc3632f92029f.tar.bz2 | |
clusters/halo.py: Simplify two coefficients usage
Signed-off-by: Aaron LI <aly@aaronly.me>
Diffstat (limited to 'fg21sim')
| -rw-r--r-- | fg21sim/extragalactic/clusters/halo.py | 12 |
1 files changed, 5 insertions, 7 deletions
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 |