G24_SMCBumps

class dust_extinction.averages.G24_SMCBumps(**kwargs)

Bases: BaseExtModel

Gordon et al (2024) SMC Bumps Extinction Curve

Parameters:

None

Raises:

None

Notes

From Gordon et al. (2024, ApJ, in press)

Two data points in the FUV from the data file giving the observed average were removed as they are very deviate from the FM90 parametrization. This cause the automated tests to fail.

Example showing the average curve

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

from dust_extinction.averages import G24_SMCBumps

fig, ax = plt.subplots()

# define the extinction model
ext_model = G24_SMCBumps()

# 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(x,ext_model(x),label='G24 SMC Bumps')
ax.plot(ext_model.obsdata_x, ext_model.obsdata_axav, 'ko',
        label='obsdata')

ax.set_xlabel(r'$x$ [$\mu m^{-1}$]')
ax.set_ylabel(r'$A(x)/A(V)$')

# for 2nd x-axis with lambda values
axis_xs = np.array([0.1, 0.12, 0.15, 0.2, 0.3, 0.5, 1.0])
new_ticks = 1 / axis_xs
new_ticks_labels = ["%.2f" % z for z in axis_xs]
tax = ax.twiny()
tax.set_xlim(ax.get_xlim())
tax.set_xticks(new_ticks)
tax.set_xticklabels(new_ticks_labels)
tax.set_xlabel(r"$\lambda$ [$\mu$m]")

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

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

Attributes Summary

Rv
x_range

Methods Summary

evaluate(x)

G24 SMCBumps function

Attributes Documentation

Rv = 2.55

x_range = [0.3, 10.0]

Methods Documentation

evaluate(x)

G24 SMCBumps 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