diff options
| -rw-r--r-- | fg21sim/configs/20-extragalactic.conf.spec | 3 | ||||
| -rw-r--r-- | fg21sim/extragalactic/clusters/psformalism.py | 5 | ||||
| -rw-r--r-- | fg21sim/utils/cosmology.py | 10 |
3 files changed, 12 insertions, 6 deletions
diff --git a/fg21sim/configs/20-extragalactic.conf.spec b/fg21sim/configs/20-extragalactic.conf.spec index 3137c73..e960a5e 100644 --- a/fg21sim/configs/20-extragalactic.conf.spec +++ b/fg21sim/configs/20-extragalactic.conf.spec @@ -30,9 +30,6 @@ # native *pickle* format (i.e., .pkl) halos_dumpfile = string(default=None) - # The fraction of the dark matter mass in galaxy clusters. - f_darkmatter = float(default=0.8, min=0.5, max=1.0) - # The minimum mass for clusters when to determine the galaxy clusters # total counts and their distributions. # Unit: [Msun] diff --git a/fg21sim/extragalactic/clusters/psformalism.py b/fg21sim/extragalactic/clusters/psformalism.py index f28c0d4..508d391 100644 --- a/fg21sim/extragalactic/clusters/psformalism.py +++ b/fg21sim/extragalactic/clusters/psformalism.py @@ -42,12 +42,11 @@ class PSFormalism: """ comp = "extragalactic/clusters" self.datafile = self.configs.get_path(comp+"/ps_data") - self.f_darkmatter = self.configs.getn(comp+"/f_darkmatter") self.Mmin_cluster = self.configs.getn(comp+"/mass_min") # [Msun] @property def Mmin_halo(self): - return self.Mmin_cluster * self.f_darkmatter + return self.Mmin_cluster * COSMO.darkmatter_fraction def _load_data(self, filepath=None): """ @@ -225,7 +224,7 @@ class PSFormalism: df = pd.DataFrame(np.column_stack([z_list, M_list]), columns=["z", "mass"]) - df["mass"] /= self.f_darkmatter + df["mass"] /= COSMO.darkmatter_fraction comment = [ "cluster number counts : %d" % counts, "z : redshift", diff --git a/fg21sim/utils/cosmology.py b/fg21sim/utils/cosmology.py index d66de09..cc347e4 100644 --- a/fg21sim/utils/cosmology.py +++ b/fg21sim/utils/cosmology.py @@ -248,6 +248,16 @@ class Cosmology: """ return self.Ob0 / self.Om0 + @property + def darkmatter_fraction(self): + """ + The cosmological mean dark matter fraction (w.r.t. matter), + assumed to be *constant* regardless of redshifts! + + See also: ``self.baryon_fraction`` + """ + return 1 - self.baryon_fraction + def overdensity_virial(self, z): """ Calculate the virial overdensity, which generally used to |