A tool for quickly selecting skylines in a sky-spectrum and then performing a wavelength dependent fit.
- Python >= 3.6
- Usual data-analysis modules (numpy, scipy, matplotlib)
- mh.spectra -- another repo of mine for working with astrophysical spectra
To start without any spectrum loaded
fit_skylines
or
fit_skylines my_spectrum.dat
to read from a text file containing a spectrum.
If you have already written some lines to disk and want to reload them
fit_skylines --readlines
Sometimes you don't have a sky-spectrum nor an arc-spectrum to use, just an observed spectrum. Fortunately the sky spectrum is often embedded in the flux errors, though strictly it is the flux variances we need to get the correct line widths. For this we have the option
fit_skylines --usevar
Other options are available for the skyline box width and polynomial order.
Description of all options
fit_skylines -h/--help