diff options
| author | Aaron LI <aly@aaronly.me> | 2017-07-23 09:51:53 +0800 |
|---|---|---|
| committer | Aaron LI <aly@aaronly.me> | 2017-07-23 09:51:53 +0800 |
| commit | bace1e88fcf2f906a1ae17c3c10004eedb176070 (patch) | |
| tree | e5ad704725ef3eba43f7825b9c78a0b93a1ab9c8 /fg21sim/extragalactic | |
| parent | e9c6dbe6d3dafff46ae2b065ff11b9f9a5a28741 (diff) | |
| download | fg21sim-bace1e88fcf2f906a1ae17c3c10004eedb176070.tar.bz2 | |
clusters/helper.py: Disable scattering on scaling relations
Signed-off-by: Aaron LI <aly@aaronly.me>
Diffstat (limited to 'fg21sim/extragalactic')
| -rw-r--r-- | fg21sim/extragalactic/clusters/helper.py | 12 |
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 |