Skip to content

Code to create the figures for the SpeX NIR dust extinction paper.

License

Notifications You must be signed in to change notification settings

mdecleir/spex_nir_extinction

Repository files navigation

SpeX NIR Extinction Paper

This repository contains code to measure NIR extinction curves for a sample of stars, based on their SpeX spectra. The scripts in the plotting folder produce the figures in the SpeX NIR dust extinction paper by Decleir et al. (submitted to ApJ).

This work is still in development. Use at your own risk.

Contributors

Marjorie Decleir

License

This project is Copyright (c) Marjorie Decleir and licensed under the terms of the BSD 3-Clause license (see the LICENSE file for more information).

Dependencies

This project uses other repositories:

  • dust_extinction provides a conversion function (AxAvToExv), average extinction curves and grain models from the literature.
  • measure_extinction contains the routines to read the data and calculate extinction curves.

The data files (including final spectra and extinction curves) related to this work, can be found on Zenodo.

Use

  • To create, scale and plot SpeX spectra of all stars in the sample, run:

    python prepare_all_spex_spectra.py
    
  • To plot all spectra (Figs. 2-4), and create IR color-color plots (Fig. 6), run:

    python plotting/plot_spex_spec.py
    
  • To calculate, fit and plot extinction curves for all the reddened stars in the sample, run:

    python calc_fit_plot_ext.py
    
  • To plot all extinction curves (Figs. 5, 7, 10, 11, 12), run:

    python plotting/plot_spex_ext.py
    
  • To create tables and plot results (Figs. 8-9, Tables 2-3), run:

    python plotting/plot_results.py
    
  • To measure, fit and plot the R(V) dependence (Figs. 13-15, Table 4), run:

    python RV_dependence.py