def luminosity_for_filter(self, filter, unit="W/micron"):
"""
This function ...
:param filter:
:param unit:
:return:
"""
#luminosity = filter.integrate(self.sed["Wavelength"], self.sed["Luminosity"])
luminosity = filter.convolve(self.sed.wavelengths(unit="micron", asarray=True), self.sed.luminosities(unit="W/micron", asarray=True)) # also in W/micron
luminosity = luminosity * Unit("W/micron")
# Return the luminosity in the desired unit
return luminosity.to(unit)
# -----------------------------------------------------------------
评论列表
文章目录