From e3fb2c4c3fc3f616295ee57bb03a7c739b5cf73c Mon Sep 17 00:00:00 2001
From: Aaron LI <aly@aaronly.me>
Date: Tue, 24 Oct 2017 16:26:53 +0800
Subject: clusters/halo: Add method "_magnetic_field()"

---
 fg21sim/extragalactic/clusters/halo.py | 25 +++++++++++++++++++++++--
 1 file changed, 23 insertions(+), 2 deletions(-)

(limited to 'fg21sim/extragalactic/clusters')

diff --git a/fg21sim/extragalactic/clusters/halo.py b/fg21sim/extragalactic/clusters/halo.py
index 1efaf8f..575decd 100644
--- a/fg21sim/extragalactic/clusters/halo.py
+++ b/fg21sim/extragalactic/clusters/halo.py
@@ -746,6 +746,28 @@ class RadioHalo:
         mass = rate * (t - t_merger) + self.M_main
         return mass
 
+    def _magnetic_field(self, t):
+        """
+        Calculate the mean magnetic field strength of the main cluster mass
+        at the given (cosmic) time.
+
+        Parameters
+        ----------
+        t : float
+            The (cosmic) time/age.
+            Unit: [Gyr]
+
+        Returns
+        -------
+        B : float
+            The mean magnetic field strength of the main cluster.
+            Unit: [uG]
+        """
+        z = COSMO.redshift(t)
+        mass = self._mass(t)  # [Msun]
+        B = helper.magnetic_field(mass=mass, z=z)  # [uG]
+        return B
+
     def _loss_ion(self, gamma, t):
         """
         Energy loss through ionization and Coulomb collisions.
@@ -783,8 +805,7 @@ class RadioHalo:
         ----------
         Ref.[sarazin1999],Eq.(6,7)
         """
+        B = self._magnetic_field(t)  # [uG]
         z = COSMO.redshift(t)
-        mass = self._mass(t)  # [Msun]
-        B = helper.magnetic_field(mass=mass, z=z)  # [uG]
         loss = -4.32e-4 * gamma**2 * ((B/3.25)**2 + (1+z)**4)
         return loss
-- 
cgit v1.2.2