aboutsummaryrefslogtreecommitdiffstats
path: root/fg21sim/extragalactic
diff options
context:
space:
mode:
authorAaron LI <aly@aaronly.me>2019-01-17 15:16:12 +0800
committerAaron LI <aly@aaronly.me>2019-01-17 15:16:12 +0800
commit576bdca6372710bdd52107e57e944d843bd89a11 (patch)
tree790087506b46d11333afc998b67ae226ce4ec21b /fg21sim/extragalactic
parentab00ba4ff010e11d40e292186d5bc15a52a9b90a (diff)
downloadfg21sim-576bdca6372710bdd52107e57e944d843bd89a11.tar.bz2
clusters/main: Account for baryon fraction in _simulate_catalog()
Diffstat (limited to 'fg21sim/extragalactic')
-rw-r--r--fg21sim/extragalactic/clusters/main.py14
-rw-r--r--fg21sim/extragalactic/clusters/psformalism.py3
2 files changed, 11 insertions, 6 deletions
diff --git a/fg21sim/extragalactic/clusters/main.py b/fg21sim/extragalactic/clusters/main.py
index 2f27c07..52e179c 100644
--- a/fg21sim/extragalactic/clusters/main.py
+++ b/fg21sim/extragalactic/clusters/main.py
@@ -115,9 +115,17 @@ class GalaxyClusters:
psform.write()
counts = psform.calc_cluster_counts(coverage=self.sky.area)
z, mass, self.comments = psform.sample_z_m(counts)
- self.catalog = []
- for z_, m_ in zip(z, mass):
- self.catalog.append(OrderedDict([("z", z_), ("mass", m_)]))
+ dm_frac = 1 - COSMO.baryon_fraction
+ self.catalog = [OrderedDict([("z", z_),
+ ("mass_dm", m_),
+ ("mass", m_ / dm_frac)])
+ for z_, m_ in zip(z, mass)]
+ self.comments += [
+ "",
+ "z - redshift",
+ "mass_dm - [Msun] dark matter halo mass",
+ "mass - [Msun] cluster total mass",
+ ]
logger.info("Simulated a catalog of %d clusters" % counts)
def _process_catalog(self):
diff --git a/fg21sim/extragalactic/clusters/psformalism.py b/fg21sim/extragalactic/clusters/psformalism.py
index dff69a5..d5c9fa5 100644
--- a/fg21sim/extragalactic/clusters/psformalism.py
+++ b/fg21sim/extragalactic/clusters/psformalism.py
@@ -298,8 +298,5 @@ class PSFormalism:
"cluster minimum mass: %.2e [Msun]" % self.Mmin,
"cluster counts: %d" % counts,
"boost factor for cluster counts: %s" % self.boost,
- "",
- "z - redshift",
- "mass - [Msun] cluster total mass",
]
return (z, mass, comment)