Extinction is not self-consistent internally
keflavich opened this issue · comments
The Extinction module wants chi
defined in cm^2/g:
https://github.com/astrofrog/sedfitter/blob/master/sedfitter/extinction/extinction.py#L41-L43
but the function get_av
treats chi
as if it is a pure opacity:
https://github.com/astrofrog/sedfitter/blob/master/sedfitter/extinction/extinction.py#L86-L87
I guess the idea is that chi
is "tau per (cm^2 / g)", so just multiplying by -0.4 turns it to a magnitude?
In any case, I suggest the following approach for generating extinction objects, based on Karl Gordon's dust extinction:
from dust_extinction.parameter_averages import F19
# a_v to N(H) from https://arxiv.org/pdf/0903.2057.pdf
# m(particle) ~ 1.34 m(H) from memory...
guyver2009_avtocol = (2.21e21 * u.cm**-2 * (1.34*u.Da)).to(u.g/u.cm**2)
ext_wav = np.sort((np.geomspace(0.301, 8.699, 1000)/u.um).to(u.um, u.spectral()))
ext_vals = ext.evaluate(ext_wav, Rv=3.1)
extinction = Extinction()
extinction.wav = ext_wav
extinction.chi = ext_vals / guyver2009_avtocol
but I'd like to know if this looks halfway reasonable?