PyWiFeS / pipeline

The Python data reduction pipeline for WiFeS
6 stars 26 forks source link

Extraction and splice of the spectra in NS and Classical mode #46

Closed felipeji closed 4 months ago

felipeji commented 4 months ago

New routines added to extract and splice the spectra from the generated data cubes in nod and shuffle and classical observing modes. The cubes' extraction is the current last step of the pipeline. Thus, these processes are added at the end of the pipeline. To do that, the main changes and additions are as follows:

  1. New file src/pywifes/extract_spec.py with new functions that will perform the extraction of both blue and red spectra from the data cubes using auto_extract(). To do that, we first detect the 3 brightest sources in the field of view for each of the blue and red data cubes excluding 3 lines/rows of pixels in the edges of the images to avoid issues in the detector (param border_width). Then, define circular apertures to extract the spectra of the sources. After that, we subtract the sky measured from an outer circular annulus aperture to the source in the aperture, but in the classical mode only (in nod and shuffle data, the sky is already subtracted earlier in the pipeline). Finally, we extract the spectra in each arm, blue and red, using auto_extract(). The function plot_apertures() will be a data product plot so the users have information of the apertures around the detected sources.

  2. New file src/pywifes/splice_spec.py with functions from previous scripts available within the team to splice the extracted spectra, that is, to join the blue and red spectra in one single spectrum. The main function is splice().

  3. New routines added in the reduction_scripts/reduce_data.py file to execute the above functions at the end of the pipeline.