skypy
skypy copied to clipboard
Add position independent dust extinction to spectrum module
Description
In order to fully model an observed galaxy catalogue, we need to include galactic extinction to our simulation. The extinction is defined as A_lambda = -2.5 log(F_lambda/F_lambda,0)
The extinction is to be included after the redshifting of the spectrum.
The reddened magnitude in band x is then given by m_x = -2.5 log(int{R_x F 10^[-0.4A_lambda] dlambda} / int{R_x F dlambda})
, where F
is the spectrum and R_x
is the filter response.
Info
To calculate the extinction we should use the Astropy affiliated package dust_extinction
. A good example of how to use it is given here. With example 1 and the beginning of example 3, everything needed should be explained.
Input
- template spectra
Outputs
- reddened band magnitudes of template spectra