change XRB gamma from 2 to 2.5 (Wu et al. 2008)

pcigale/sed_modules/xray.py

@@ -58,7 +58,6 @@ class Xray(SedModule):
         self.lam_10keV = lam_10keV
         lam_100keV = c*cst.h / (1e5*cst.eV)
         # Define frequency corresponding to 2 keV in units of Hz.
-        nu_0p5keV = c / lam_0p5keV
         nu_2keV = c / lam_2keV
 
         # We define the wavelength range for the non thermal emission
@@ -73,9 +72,9 @@ class Xray(SedModule):
         # also have two components:
         #   2.1 high-mass X-ray binaries (HMXB)
         #   2.2 low-mass X-ray binaries (LMXB)
-        # Assuming gamma=2 for both components (Lehmer et al. 2016)
+        # Assuming gamma=2.5 for both components (Wu et al. 2008)
         # normalized such that L(2-10 keV) = 1.
-        self.lumin_xrb = self.wave**-1 / np.log(5)
+        self.lumin_xrb = self.wave**-0.5 / (2*lam_2keV**0.5 - 2*lam_10keV**0.5)
 
         # We compute the unobscured AGN corona X-ray emission
         # The shape is power-law with high-E exp. cutoff