Commit 67550872 authored by Guang's avatar Guang Committed by Guang
Browse files

only apply E(B-V) to type 1

parent a5182b96
...@@ -138,9 +138,14 @@ class Fritz2006(SedModule): ...@@ -138,9 +138,14 @@ class Fritz2006(SedModule):
# The extinction factor, flux_1/flux_0 # The extinction factor, flux_1/flux_0
ext_fac = 10**(A_lam/-2.5) ext_fac = 10**(A_lam/-2.5)
# Calculate the new AGN SED shape after extinction # Calculate the new AGN SED shape after extinction
# The direct and scattered components (line-of-sight) are extincted if self.psy > (90-self.opening_angle):
lumin_agn_new = self.fritz2006.lumin_agn * ext_fac # The direct and scattered components (line-of-sight) are extincted for type-1 AGN
lumin_scatt_new = self.fritz2006.lumin_scatt * ext_fac lumin_agn_new = self.fritz2006.lumin_agn * ext_fac
lumin_scatt_new = self.fritz2006.lumin_scatt * ext_fac
else:
# Keep the direct and scatter components for type-2
lumin_agn_new = self.fritz2006.lumin_agn
lumin_scatt_new = self.fritz2006.lumin_scatt
# Calculate the total extincted luminosity averaged over all directions # Calculate the total extincted luminosity averaged over all directions
# note that self.opening_angle is different from Fritz's definiting! # note that self.opening_angle is different from Fritz's definiting!
l_ext = np.trapz(self.fritz2006.lumin_intrin_agn * (1-ext_fac), l_ext = np.trapz(self.fritz2006.lumin_intrin_agn * (1-ext_fac),
......
Supports Markdown
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