From a57a01fab43c6eb692fc04ccb0f0c348a1eba3b1 Mon Sep 17 00:00:00 2001 From: Aaron LI Date: Wed, 3 Jan 2018 23:30:17 +0800 Subject: clusters/halo: update radius_halo() calculation --- fg21sim/extragalactic/clusters/helper.py | 12 ++++++------ 1 file changed, 6 insertions(+), 6 deletions(-) (limited to 'fg21sim/extragalactic/clusters/helper.py') diff --git a/fg21sim/extragalactic/clusters/helper.py b/fg21sim/extragalactic/clusters/helper.py index 88c566a..39cee12 100644 --- a/fg21sim/extragalactic/clusters/helper.py +++ b/fg21sim/extragalactic/clusters/helper.py @@ -89,23 +89,23 @@ def radius_virial(mass, z=0.0): return R_vir -def radius_halo(M_main, M_sub, z=0.0, configs=CONFIGS): +def radius_halo(mass, z=0.0, configs=CONFIGS): """ - Calculate the (predicted) radius of (giant) radio halo. + Estimate the radius of (giant) radio halo. NOTE ---- The halo radius is estimated to be the same as the turbulence injection scale, i.e.: R_halo ≅ L ≅ R_vir / 3 - where R_vir the virial radius of the merged cluster. + where R_vir the virial radius of the merged (observed) cluster. Reference: [vazza2011],Sec.(3.6) Parameters ---------- - M_main, M_sub : float, `~numpy.ndarray` - Total (virial) masses of the main and sub clusters + mass : float, `~numpy.ndarray` + Cluster virial mass. Unit: [Msun] z : float, `~numpy.ndarray`, optional Redshift @@ -120,7 +120,7 @@ def radius_halo(M_main, M_sub, z=0.0, configs=CONFIGS): # Turbulence injection scale factor key = "extragalactic/halos/f_lturb" f_lturb = configs.getn(key) - R_halo = f_lturb * radius_virial(mass=M_main+M_sub, z=z) # [kpc] + R_halo = f_lturb * radius_virial(mass=mass, z=z) # [kpc] return R_halo -- cgit v1.2.2