G21_MWAvg

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

Bases: 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(x)

G21_MWAvg function

Attributes Documentation

Rv = 3.17
x_range = [0.03125, 1.0]

Methods Documentation

evaluate(x)[source]

G21_MWAvg function

Parameters:
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