G24_SMCBumps¶
- class dust_extinction.averages.G24_SMCBumps(**kwargs)[source]¶
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
Methods Summary
evaluate
(in_x)G24 SMCBumps function
Attributes Documentation
- Rv = 2.55¶
- x_range = [0.3, 10.0]¶
Methods Documentation
- evaluate(in_x)[source]¶
G24 SMCBumps 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