diff options
| -rw-r--r-- | fg21sim/extragalactic/clusters/main.py | 4 | ||||
| -rw-r--r-- | fg21sim/extragalactic/clusters/psformalism.py | 36 |
2 files changed, 13 insertions, 27 deletions
diff --git a/fg21sim/extragalactic/clusters/main.py b/fg21sim/extragalactic/clusters/main.py index 0251662..90aaf75 100644 --- a/fg21sim/extragalactic/clusters/main.py +++ b/fg21sim/extragalactic/clusters/main.py @@ -116,8 +116,10 @@ class GalaxyClusters: psform.calc_dndlnm() psform.write() counts = psform.calc_cluster_counts(coverage=self.sky.area) - self.catalog, self.catalog_comment = psform.sample_z_m(counts) logger.info("Simulated cluster catalog of counts %d." % counts) + z, mass, self.comments = psform.sample_z_m(counts) + self.catalog = pd.DataFrame(np.column_stack([z, mass]), + columns=["z", "mass"]) def _process_catalog(self): """ diff --git a/fg21sim/extragalactic/clusters/psformalism.py b/fg21sim/extragalactic/clusters/psformalism.py index 368349c..7c824d8 100644 --- a/fg21sim/extragalactic/clusters/psformalism.py +++ b/fg21sim/extragalactic/clusters/psformalism.py @@ -1,4 +1,4 @@ -# Copyright (c) 2017 Weitian LI <weitian@aaronly.me> +# Copyright (c) 2017-2018 Weitian LI <weitian@aaronly.me> # MIT license """ @@ -15,7 +15,6 @@ import random from functools import lru_cache import numpy as np -import pandas as pd import hmf from ...share import CONFIGS, COSMO @@ -209,10 +208,6 @@ class PSFormalism: ------- counts : int The total number of clusters within the sky patch. - - Attributes - ---------- - counts """ logger.info("Calculating the total number of clusters within " "sky patch of coverage %.1f [deg^2]" % coverage) @@ -221,13 +216,12 @@ class PSFormalism: coverage *= AUC.deg2rad**2 # [deg^2] -> [rad^2] = [sr] midx = (self.mass >= self.Mmin_halo) numgrid = self.number_grid - counts = np.sum(numgrid[:, midx]) * coverage - counts *= self.boost # number boost factor - self.counts = int(np.round(counts)) - logger.info("Total number of clusters: %d" % self.counts) - return self.counts + counts = np.sum(numgrid[:, midx]) * coverage * self.boost + counts = int(np.round(counts)) + logger.info("Total number of clusters: %d" % counts) + return counts - def sample_z_m(self, counts=None): + def sample_z_m(self, counts): """ Randomly generate the requested number of pairs of (z, M) following the halo number distribution. @@ -249,7 +243,6 @@ class PSFormalism: ---------- counts : int, optional The number of (z, mass) pairs to be sampled. - If not specified, then default to ``self.counts`` Returns ------- @@ -257,14 +250,7 @@ class PSFormalism: A Pandas data frame with 2 columns, i.e., ``z`` and ``mass``. comment : list[str] Comments to the above data frame. - - Attributes - ---------- - clusters : df - clusters_comment : comment """ - if counts is None: - counts = self.counts logger.info("Sampling (z, mass) pairs for %d clusters ..." % counts) z = self.z @@ -310,18 +296,16 @@ class PSFormalism: i += 1 logger.info("Sampled %d pairs of (z, mass) for each cluster" % counts) - df = pd.DataFrame(np.column_stack([z_list, mass_list]), - columns=["z", "mass"]) - df["mass"] /= COSMO.darkmatter_fraction + z = np.array(z_list) + mass = np.array(mass_list) / COSMO.darkmatter_fraction comment = [ "halo mass function model: %s" % self.hmf_model, "cluster minimum mass: %.2e [Msun]" % self.Mmin_cluster, "dark matter fraction: %.2f" % COSMO.darkmatter_fraction, "cluster counts: %d" % counts, "boost factor for cluster counts: %s" % self.boost, + "", "z - redshift", "mass - cluster total mass [Msun]", ] - self.clusters = df - self.clusters_comment = comment - return (df, comment) + return (z, mass, comment) |