Commit 6a57af0d authored by Guang's avatar Guang
Browse files

change the cut-off energy from 300keV to 100keV

parent c3b1e2c6
...@@ -54,7 +54,7 @@ class Xray(SedModule): ...@@ -54,7 +54,7 @@ class Xray(SedModule):
# Define wavelenght corresponding to some energy in units of nm. # Define wavelenght corresponding to some energy in units of nm.
lam_0p5keV = c*cst.h / (5e2*cst.eV) lam_0p5keV = c*cst.h / (5e2*cst.eV)
lam_2keV = c*cst.h / (2e3*cst.eV) lam_2keV = c*cst.h / (2e3*cst.eV)
lam_100keV = c*cst.h / (1e5*cst.eV) lam_300keV = c*cst.h / (3e5*cst.eV)
self.lam_2keV = lam_2keV self.lam_2keV = lam_2keV
self.lam_10keV = c*cst.h / (1e4*cst.eV) self.lam_10keV = c*cst.h / (1e4*cst.eV)
# Define frequency corresponding to 2 keV in units of Hz. # Define frequency corresponding to 2 keV in units of Hz.
...@@ -78,11 +78,12 @@ class Xray(SedModule): ...@@ -78,11 +78,12 @@ class Xray(SedModule):
# We compute the unobscured AGN corona X-ray emission # We compute the unobscured AGN corona X-ray emission
# The shape is power-law with high-E exp. cutoff # The shape is power-law with high-E exp. cutoff
# with cut-off E=300 keV (Ueda+2014; Aird+2015)
self.lumin_corona = self.wave**(self.gam - 3.) * \ self.lumin_corona = self.wave**(self.gam - 3.) * \
np.exp(-lam_100keV/self.wave) np.exp(-lam_300keV/self.wave)
# Normaliz the SED at 2 keV # Normaliz the SED at 2 keV
self.lumin_corona /= lam_2keV**(self.gam - 3.) * \ self.lumin_corona /= lam_2keV**(self.gam - 3.) * \
np.exp(-lam_100keV/lam_2keV) np.exp(-lam_300keV/lam_2keV)
def process(self, sed): def process(self, sed):
"""Add the X-ray contribution. """Add the X-ray contribution.
......
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