aboutsummaryrefslogtreecommitdiffstats
path: root/fg21sim
diff options
context:
space:
mode:
authorAaron LI <aly@aaronly.me>2017-07-23 09:51:53 +0800
committerAaron LI <aly@aaronly.me>2017-07-23 09:51:53 +0800
commitbace1e88fcf2f906a1ae17c3c10004eedb176070 (patch)
treee5ad704725ef3eba43f7825b9c78a0b93a1ab9c8 /fg21sim
parente9c6dbe6d3dafff46ae2b065ff11b9f9a5a28741 (diff)
downloadfg21sim-bace1e88fcf2f906a1ae17c3c10004eedb176070.tar.bz2
clusters/helper.py: Disable scattering on scaling relations
Signed-off-by: Aaron LI <aly@aaronly.me>
Diffstat (limited to 'fg21sim')
-rw-r--r--fg21sim/extragalactic/clusters/helper.py12
1 files changed, 8 insertions, 4 deletions
diff --git a/fg21sim/extragalactic/clusters/helper.py b/fg21sim/extragalactic/clusters/helper.py
index 5963bab..67d72de 100644
--- a/fg21sim/extragalactic/clusters/helper.py
+++ b/fg21sim/extragalactic/clusters/helper.py
@@ -84,8 +84,10 @@ def radius_halo(mass, z=0.0):
----------
Ref.[cassano2007],Fig.(11)
"""
- slope = 2.63 + np.random.normal(scale=0.5)
- intercept = 2.3 + np.random.normal(scale=0.05)
+ # slope = 2.63 + np.random.normal(scale=0.5)
+ slope = 2.63
+ # intercept = 2.3 + np.random.normal(scale=0.05)
+ intercept = 2.3
R_vir = radius_virial(mass=mass, z=z) # [kpc]
R_halo = 10 ** (slope * np.log10(R_vir) + intercept)
return R_halo
@@ -121,8 +123,10 @@ def mass_to_kT(mass, z=0.0):
The ICM mean temperature.
Unit: [keV]
"""
- A = 5.34 + np.random.normal(scale=0.22)
- alpha = 1.72 + np.random.normal(scale=0.10)
+ # A = 5.34 + np.random.normal(scale=0.22)
+ A = 5.34
+ # alpha = 1.72 + np.random.normal(scale=0.10)
+ alpha = 1.72
Ez = COSMO.E(z)
kT = 5.0 * (mass * Ez / A) ** (1/alpha)
return kT