aboutsummaryrefslogtreecommitdiffstats
path: root/fg21sim
diff options
context:
space:
mode:
Diffstat (limited to 'fg21sim')
-rw-r--r--fg21sim/extragalactic/clusters/helper.py14
1 files changed, 11 insertions, 3 deletions
diff --git a/fg21sim/extragalactic/clusters/helper.py b/fg21sim/extragalactic/clusters/helper.py
index 8eef6b9..4d06ceb 100644
--- a/fg21sim/extragalactic/clusters/helper.py
+++ b/fg21sim/extragalactic/clusters/helper.py
@@ -76,7 +76,7 @@ def beta_model(rho0, rc, beta):
return func
-def calc_gas_density_profile(mass, z):
+def calc_gas_density_profile(mass, z, f_rc=0.1, beta=0.8):
"""
Calculate the parameters of the β-model that is used to describe the
gas density profile.
@@ -88,6 +88,15 @@ def calc_gas_density_profile(mass, z):
Reference: [cassano2005],Sec.(4.1)
+ Parameters
+ ----------
+ f_rc : float
+ The fraction of the core radius to the virial radius.
+ Default: 0.1
+ beta : float
+ The slope parameter of the β-model.
+ Default: 0.8
+
Returns
-------
fbeta : function
@@ -95,8 +104,7 @@ def calc_gas_density_profile(mass, z):
Unit: [Msun/kpc^3]
"""
r_vir = radius_virial(mass, z) # [kpc]
- rc = 0.1 * r_vir
- beta = 0.8
+ rc = f_rc * r_vir
fint = beta_model(1, rc, beta)
v = integrate.quad(lambda r: fint(r) * r**2,
a=0, b=r_vir)[0] # [kpc^3]