From 9631363cf5d3c9d6c37979e0d7a79f1f4c7d849f Mon Sep 17 00:00:00 2001
From: Aaron LI <aly@aaronly.me>
Date: Wed, 26 Jul 2017 23:18:25 +0800
Subject: clusters/emission.py: Fix units/dimensions problem
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit
NOTE:
Since we just substitute ``p`` (electron momentum) for
``γ`` (electron Lorentz factor), i.e.,
N = int_p n(p) dp = int_γ n(γ) dγ ,
therefore the coefficient ``sqrt(3) e^3 B / (m c^2)`` does NOT need
modification (e.g., multiply an extra ``m c``)!
Signed-off-by: Aaron LI <aly@aaronly.me>
---
fg21sim/extragalactic/clusters/emission.py | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
(limited to 'fg21sim')
diff --git a/fg21sim/extragalactic/clusters/emission.py b/fg21sim/extragalactic/clusters/emission.py
index 9e097d8..2454644 100644
--- a/fg21sim/extragalactic/clusters/emission.py
+++ b/fg21sim/extragalactic/clusters/emission.py
@@ -185,7 +185,8 @@ class SynchrotronEmission:
# Integrate over energy ``gamma`` in logarithmic grid
j_nu = integrate.simps(s1d*self.gamma, np.log(self.gamma))
- coef = np.sqrt(3) * AC.e**3 * self.B / AC.c
+ B_gauss = self.B * 1e-6 # [uG] -> [G]
+ coef = np.sqrt(3) * AC.e**3 * B_gauss / AU.mec2
j_nu *= coef
return j_nu
--
cgit v1.2.2