search

bd-j / prospector

Python code for Stellar Population Inference from Spectra and SEDs
http://prospect.readthedocs.io
MIT License
157 stars 73 forks source link

Extracting emission line fluxes/ratios #316

Closed robersju closed 8 months ago

robersju commented 8 months ago

Running models with both photometry and spectroscopy of an active star-forming galaxy and am including nebular emission to fit six specific emission lines. The model itself (as far as I can tell) is properly fitting the lines and including them in the MAP spectrum, and a best-fit eline_sigma_velocity is available - this is a question about where to actually pull the emission line fluxes. Is it somewhere within the result dictionary, or is there an alternative function to call?

The nebular emission portion of the model: 

add_neb_emission: [ True] 
add_neb_continuum: [ True] 
nebemlineinspec: [False] 
marginalize_elines: [ True] 
use_eline_prior: [ True] 
lines_to_fit: ['[OII]3726' '[OII]3729' 'H beta 4861' '[OIII]4960' '[OIII]5007''H alpha 6563'] 
eline_prior_width: [0.2] 
eline_sigma: [600.]
bd-j commented 8 months ago

Hi sorry, I missed this - hope you got it sorted, but the emission line fluxes can be accessed in SpecModel()._eline_lums after a call to SpecModel().predict().