G21_MWAvg

class dust_extinction.averages.G21_MWAvg(**kwargs)[source]

Bases: dust_extinction.baseclasses.BaseExtModel

Gordon et al. (2021) Milky Way Average Extinction Curve

Parameters
None
Raises
None

Notes

From Gordon et al. (2021, ApJ, 916, 33)

Example showing the average curve

import numpy as np
import matplotlib.pyplot as plt
import astropy.units as u

from dust_extinction.averages import G21_MWAvg

fig, ax = plt.subplots()

# generate the curves and plot them
lam = np.logspace(np.log10(1.01), np.log10(31.9), num=1000)
x = (1.0/lam)/u.micron

# define the extinction model
ext_model = G21_MWAvg()

ax.plot(1.0/x,ext_model(x),label='G21_MWAvg')
ax.errorbar(1.0/ext_model.obsdata_x_irs, ext_model.obsdata_axav_irs,
            yerr=ext_model.obsdata_axav_unc_irs,
            fmt='ko', label='obsdata (IRS)')
ax.errorbar(1.0/ext_model.obsdata_x_bands, ext_model.obsdata_axav_bands,
            yerr=ext_model.obsdata_axav_unc_bands,
            fmt='g^', label='obsdata (photometry)')

ax.set_xlabel(r'$\lambda$ [$\mu m$]')
ax.set_ylabel(r'$A(x)/A(V)$')

ax.set_xscale('log')
ax.set_yscale('log')

ax.legend(loc='best')
plt.show()

(Source code, png, hires.png, pdf)

../_images/dust_extinction-averages-G21_MWAvg-1.png

Attributes Summary

Rv

x_range

Methods Summary

evaluate(in_x)

G21_MWAvg function

Attributes Documentation

Rv = 3.17
x_range = [0.03125, 1.0]

Methods Documentation

evaluate(in_x)[source]

G21_MWAvg function

Parameters
in_x: float

expects either x in units of wavelengths or frequency or assumes wavelengths in wavenumbers [1/micron]

internally wavenumbers are used

Returns
axav: np array (float)

A(x)/A(V) extinction curve [mag]

Raises
ValueError

Input x values outside of defined range