G03_SMCBar_WD01_extension

class beast.physicsmodel.dust.extinction_extension.G03_SMCBar_WD01_extension(*args, meta=None, name=None, **kwargs)

Bases: G03_SMCBar

dust_extinction.averages.G03_SMCBar model extended to shorter wavelengths using the dust_extinction.grain_models.WD01 SMCBar model.

Parameters:

Rv: float

R(V) = A(V)/E(B-V) = total-to-selective extinction

Raises:

InputParameterError

Input Rv values outside of defined range

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

from beast.physicsmodel.dust.extinction_extension import G03_SMCBar_WD01_extension
from dust_extinction.grain_models import WD01

fig, ax = plt.subplots()

# define the extinction model
ext_model = G03_SMCBar_WD01_extension()

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

ax.plot(1.0 / x, ext_model(x), label="G03 SMCBar WD01 ext")

dmod = WD01(modelname="SMCBar")
ax.plot(
    1.0 / x, dmod(x), label="WD01 SMCBar", linestyle="dashed", color="black"
)

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

ax.set_xscale("log")

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

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

Attributes Summary

x_range

Methods Summary

evaluate(in_x)

G03_SMCBar_WD01_extension function

Attributes Documentation

x_range = [0.3, 100.0]

Methods Documentation

evaluate(in_x)

G03_SMCBar_WD01_extension 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