Commit 9e62ab9d authored by Guang's avatar Guang
Browse files

Set IGM absorption to 0 for X-ray wavelength

parent d278317b
......@@ -51,6 +51,7 @@ def igm_transmission(wavelength, redshift):
gamma = 0.2788 # Gamma(0.5,1) i.e., Gamma(2-beta,1) with beta = 1.5
n0 = 0.25
lambda_limit = 91.2 # Lyman limit in nm
lam_xray_max = 10**0.5 # Maximum X-ray wavelength in nm
lambda_n = np.empty(n_transitions_max)
z_n = np.empty((n_transitions_max, len(wavelength)))
......@@ -132,6 +133,9 @@ def igm_transmission(wavelength, redshift):
igm_transmission = np.exp(-tau_taun-tau_l_igm-tau_l_lls) * weight
# Reset the transmission for X-ray (E>0.4 keV) to 1
igm_transmission[np.where(wavelength <= lam_xray_max)] = 1.
return igm_transmission
......
......@@ -57,6 +57,7 @@ class Xray(SedModule):
nu_2keV = c / lam_2keV
# We define the wavelength range for the non thermal emission
# corresponding to 0.4-1200 keV
self.wave = np.logspace(-3, 0.5, 1000.)
# We compute the unobscured AGN corona X-ray emission
# normalized at 2 keV
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment