Commit 82d2534a authored by Guang's avatar Guang
Browse files

calculate AGN intrinsic L_nu(2500A) for Fritz+06 model

parent c82236f2
......@@ -17,6 +17,7 @@ import numpy as np
from pcigale.data import Database
from . import SedModule
import scipy.constants as cst
class Fritz2006(SedModule):
"""Fritz et al. (2006) AGN dust torus emission
......@@ -124,6 +125,10 @@ class Fritz2006(SedModule):
sed.add_info('agn.psy', self.psy)
sed.add_info('agn.fracAGN', self.fracAGN)
# We define various constants necessary to compute the model. For
# consistency, we define the speed of light in nm s¯¹ rather than in
# m s¯¹.
c = cst.c * 1e9
# Compute the AGN luminosity
if self.fracAGN < 1.:
agn_power = luminosity * (1./(1.-self.fracAGN) - 1.)
......@@ -136,6 +141,12 @@ class Fritz2006(SedModule):
# Intrinsic (de-reddened) AGN luminosity from the central source
l_agn_intrin_agn = np.trapz(agn_power * self.fritz2006.lumin_intrin_agn,
x=self.fritz2006.wave)
# Intrinsic L_nu at rest-frame 2500A
# Calculate L_lam(2500A)
l_agn_2500A = np.interp(250, self.fritz2006.wave,
agn_power * self.fritz2006.lumin_intrin_agn)
# Convert L_lam to L_nu
l_agn_2500A *= 250**2/c
else:
raise Exception("AGN fraction is exactly 1. Behaviour "
......@@ -146,6 +157,7 @@ class Fritz2006(SedModule):
sed.add_info('agn.agn_luminosity', l_agn_agn, True)
sed.add_info('agn.luminosity', l_agn_total, True)
sed.add_info('agn.agn_intrin_luminosity', l_agn_intrin_agn, True)
sed.add_info('agn.agn_intrin_Lnu_2500A', l_agn_2500A, True)
sed.add_contribution('agn.fritz2006_therm', self.fritz2006.wave,
agn_power * self.fritz2006.lumin_therm)
......
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