RRP89_MWGC

class dust_extinction.averages.RRP89_MWGC(*args, meta=None, name=None, **kwargs)[source]

Bases: BaseExtModel

Reike, Rieke, & Paul (1989) MW Average Extinction Curve

Parameters:
None
Raises:
None

Notes

From Rieke, Rieke, & Paul (1989, ApJ, 336, 752)

Example showing the average curve

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

from dust_extinction.averages import RRP89_MWGC

fig, ax = plt.subplots()

# define the extinction model
ext_model = RRP89_MWGC()

# generate the curves and plot them
x = np.arange(1.0/ext_model.x_range[1], 1.0/ext_model.x_range[0], 0.1) * u.micron

ax.plot(x,ext_model(x),label='RRP89_MWGC')
ax.plot(1.0/ext_model.obsdata_x, ext_model.obsdata_axav, 'ko',
        label='obsdata')

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

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

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

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

Attributes Summary

Rv

obsdata_axav

obsdata_elvebv

obsdata_tolerance

obsdata_x

x_range

Methods Summary

evaluate(in_x)

RRP89 MWGC function

Attributes Documentation

Rv = 3.09
obsdata_axav = array([0.48220065, 0.2815534 , 0.17475728, 0.11197411, 0.06796117,        0.03236246, 0.02588997, 0.11650485, 0.07119741, 0.08090615,        0.03559871])
obsdata_elvebv = array([-1.6  , -2.22 , -2.55 , -2.744, -2.88 , -2.99 , -3.01 , -2.73 ,        -2.87 , -2.84 , -2.98 ])
obsdata_tolerance = 1e-06
obsdata_x = array([1.11111111, 0.8       , 0.625     , 0.45454545, 0.28571429,        0.20833333, 0.125     , 0.10526316, 0.09433962, 0.09090909,        0.07692308])
x_range = [0.07692307692307693, 1.1111111111111112]

Methods Documentation

evaluate(in_x)[source]

RRP89 MWGC 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