Use sigma_T from astropy constants instead of hardcoding it, also fix…

…ing its unit.
ftsamis · Dec 17, 2016 · 136f667 · 136f667
1 parent 782a0ab
commit 136f667
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/tardis/montecarlo/base.py b/tardis/montecarlo/base.py
@@ -354,10 +354,10 @@ def from_config(cls, config):
         """
         if config.plasma.disable_electron_scattering:
             logger.warn('Disabling electron scattering - this is not physical')
-            sigma_thomson = 1e-200 / (u.cm ** 2)
+            sigma_thomson = 1e-200 * (u.cm ** 2)
         else:
             logger.debug("Electron scattering switched on")
-            sigma_thomson = 6.652486e-25 / (u.cm ** 2)
+            sigma_thomson = const.sigma_T.cgs
 
         spectrum_frequency = quantity_linspace(
             config.spectrum.stop.to('Hz', u.spectral()),