diff options
| author | Aaron LI <aly@aaronly.me> | 2017-07-26 23:18:25 +0800 |
|---|---|---|
| committer | Aaron LI <aly@aaronly.me> | 2017-07-26 23:18:25 +0800 |
| commit | 9631363cf5d3c9d6c37979e0d7a79f1f4c7d849f (patch) | |
| tree | 75a5e401a27e0877af1f515f589a105b7e09e024 | |
| parent | fe714b495697cad9caa3362e9222fa7c910056de (diff) | |
| download | fg21sim-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>
| -rw-r--r-- | fg21sim/extragalactic/clusters/emission.py | 3 |
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 |