aboutsummaryrefslogtreecommitdiffstats
path: root/fg21sim
diff options
context:
space:
mode:
authorAaron LI <aly@aaronly.me>2017-07-26 23:18:25 +0800
committerAaron LI <aly@aaronly.me>2017-07-26 23:18:25 +0800
commit9631363cf5d3c9d6c37979e0d7a79f1f4c7d849f (patch)
tree75a5e401a27e0877af1f515f589a105b7e09e024 /fg21sim
parentfe714b495697cad9caa3362e9222fa7c910056de (diff)
downloadfg21sim-9631363cf5d3c9d6c37979e0d7a79f1f4c7d849f.tar.bz2
clusters/emission.py: Fix units/dimensions problem
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>
Diffstat (limited to 'fg21sim')
-rw-r--r--fg21sim/extragalactic/clusters/emission.py3
1 files changed, 2 insertions, 1 deletions
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