utils/cosmology.py: Add property "baryon_fraction"

XXX: assumed to be independent of redshifts!

Signed-off-by: Aaron LI <aly@aaronly.me>

1 files changed, 9 insertions, 0 deletions

diff --git a/fg21sim/utils/cosmology.py b/fg21sim/utils/cosmology.py
index 9da3d80..e4979b1 100644
--- a/fg21sim/utils/cosmology.py
+++ b/fg21sim/utils/cosmology.py
@@ -199,6 +199,15 @@ class Cosmology:
         """
         return self.Om0 * (1+z)**3 / self.E(z)**2
 
+    @property
+    def baryon_fraction(self):
+        """
+        The cosmological mean baryon fraction (w.r.t. matter).
+
+        XXX: assumed to be *constant* regardless of redshifts!
+        """
+        return self.Ob0 / self.Om0
+
     def overdensity_virial(self, z):
         """
         Calculate the virial overdensity, which generally used to