Robust support for semi-automated reductions of Keck/NIRSPEC data using PypeIt

Adolfo S. Carvalho,Greg Doppmann,Kyle B. Westfall,Debora Pelliccia,J. Xavier Prochaska,Joseph Hennawi,Frederick B. Davies,Max Brodheim,Feige Wang,Ryan Cooke

2024-10-26

Abstract:We present a data reduction pipeline (DRP) for Keck/NIRSPEC built as an addition to the PypeIt Python package. The DRP is capable of reducing multi-order echelle data taken both before and after the detector upgrade in 2018. As part of developing the pipeline, we implemented major improvements to the capabilities of the PypeIt package, including manual wavelength calibration for multi-order data and new output product that returns a coadded spectrum order-by-order. We also provide a procedure for correcting telluric absorption in NIRSPEC data by using the spectra of telluric standard stars taken near the time of the science spectra. At high resolutions, this is often more accurate than modeling-based approaches.

Instrumentation and Methods for Astrophysics

What problem does this paper attempt to address?

This paper aims to solve the problems encountered in Keck/NIRSPEC data processing, especially after the detector upgrade in 2018, the original data processing software could not effectively process the data of the new detector. Specifically, the paper attempts to solve the following key issues:

1. **Lack of effective data processing tools**:
- Before the detector upgrade in 2018, NIRSPEC users had three recommended software options for data processing: REDSPEC, WMKONSPEC, and NSDRP. However, after the detector upgrade, only REDSPEC was updated to support the data processing of the new detector, while NSDRP and WMKONSPEC could only be used to process the data before the upgrade. This led to the need for new data processing tools.

2. **Improve the automation and reliability of data processing**:
- In order to enable NIRSPEC users to quickly process new data and archived data, researchers developed a Python data processing pipeline (DRP) based on PypeIt. PypeIt is an open - source Python package that can automate the processing of raw observational data and calibration data from slit spectrometers. By extending the functions of PypeIt, researchers hope to provide a more powerful and reliable tool to process multi - order echelle data.

3. **Improve wavelength calibration and output format**:
- Researchers made two significant improvements to PypeIt:
- **Manual wavelength calibration**: Extended the PypeIt manual wavelength calibration procedure to enable it to process multi - order echelle data. The new procedure can process the wavelength calibration of each order one by one through the user interface and save the solutions.
- **New output format**: Introduced a new output format for high - resolution echelle data, that is, the combined spectra of each order are output in the form of an \( N\times M \) array, where \( N \) is the number of spectral orders and \( M \) is the number of columns on the detector. This format is more convenient when processing high - resolution spectra.

4. **Eliminate the influence of atmospheric absorption**:
- The paper also introduced a method to correct the atmospheric absorption in NIRSPEC data using the spectrum of a standard A0V star. The specific steps include:
1. Continuum normalization of the standard star and scientific target spectra.
2. Fit and remove the stellar atmospheric features of the standard star.
3. Correct the differences between the standard star and scientific target spectra caused by the instrument wavelength drift.
4. Uniformly scale the depth of atmospheric absorption spectral lines according to the difference in zenith distance at the time of observation.
5. Divide the corrected standard star spectrum from the scientific target spectrum to eliminate the atmospheric absorption features.

Through these improvements, researchers hope to provide a more efficient and accurate data processing tool to help NIRSPEC users better process and analyze their observational data.