aboutsummaryrefslogtreecommitdiffstats
diff options
context:
space:
mode:
-rw-r--r--fg21sim/configs/config.spec7
-rw-r--r--fg21sim/extragalactic/clusters/halo.py5
2 files changed, 11 insertions, 1 deletions
diff --git a/fg21sim/configs/config.spec b/fg21sim/configs/config.spec
index fefb8b8..85606eb 100644
--- a/fg21sim/configs/config.spec
+++ b/fg21sim/configs/config.spec
@@ -447,3 +447,10 @@ stream = option("stderr", "stdout", "", default="stderr")
# Unit: [MHz]
fiducial_freq = float(default=150.0)
fiducial_factor = float(default=3.0, min=1)
+
+ # Parameters of the beta-model that is used to describe the gas density
+ # profile of the cluster.
+ # The fraction of the core radius to cluster's virial radius.
+ f_rc = float(default=0.1)
+ # The slope parameter (i.e., beta).
+ beta = float(default=0.8)
diff --git a/fg21sim/extragalactic/clusters/halo.py b/fg21sim/extragalactic/clusters/halo.py
index 113e434..a18aa17 100644
--- a/fg21sim/extragalactic/clusters/halo.py
+++ b/fg21sim/extragalactic/clusters/halo.py
@@ -167,6 +167,8 @@ class RadioHalo:
self.injection_index = configs.getn(comp+"/injection_index")
self.fiducial_freq = configs.getn(comp+"/fiducial_freq")
self.fiducial_factor = configs.getn(comp+"/fiducial_factor")
+ self.f_rc = configs.getn(comp+"/f_rc")
+ self.beta = configs.getn(comp+"/beta")
def _set_solver(self):
self.fpsolver = FokkerPlanckSolver(
@@ -673,7 +675,8 @@ class RadioHalo:
z = COSMO.redshift(t)
M_main = self.mass_main(t)
M_sub = self.mass_sub(t)
- return helper.calc_gas_density_profile(mass=M_main+M_sub, z=z)
+ return helper.calc_gas_density_profile(mass=M_main+M_sub, z=z,
+ f_rc=self.f_rc, beta=self.beta)
@lru_cache()
def _velocity_turb(self, t):