Menu:


Update:

16th January 2024







  back 

With this function, one can read the standard output fits file of the TIGRE/HEROS automatic data reduction pipeline. The notation of the main output fits file is Sci_...sp_ech_main.fits. A secondary output is a small fits file. The notation is Sci_...sp_ech_small.fits. This file include only the blaze normalised spectrum.
In this function, the barycentric velocity correction is performed automatically if one don't set barrcorr=no. The barycentric velosity correction is estimated during the reduction and listed in the image header. Furthermore, one can also commit a radial velocity and the wavelength is shifted by this value.
The syntax is:

function readheros(file [, radvel=Value, barycorr=no, small=yes])

Input:
          file: Name of the fits file
Optional:
          radvel: Radial velocity [km/s]
          barycorr: If one set no than the radial velocity and barycentric correction are not performed
          small: One use small=yes if one read the small fits file
Returns:
For a spectrum from 1 exposure and from more than 1 where the single images are co-add:
          result[0]: Header of the image
          result[1]: Header of the spectrum
          result[2]: Array of the wavelength of each order
          result[3]: Array of the blaze normalised spectrum of each order
          result[4]: Array of the blaze function (flat field spectrum) of each order
          result[5]: Array of the signal to noise of each order
          result[6]: Array of the normalised spectrum of each order
          result[7]: Array of the merged wavelength
          result[8]: Array of the merged blaze normalised spectrum
          result[9]: Array of the error of the merged blaze normelised spectrum
          result[10]: Array of the merged normalised spectrum
          result[11]: Array of the error of the merged normalised spectrum
          result[12]: Array of the global response curve
          result[13]: Array of the correction of the small scale variations in the spectrum

For a spectrum from more than 1 exposure and the single spectra are available:
          result[0]: Header of the spectrum
          result[1]: Header of the image
          result[2]: Array of the wavelength of each order
          result[3]: Array of the blaze normalised spectrum of each order
          result[4]: Array of the blaze function (flat field spectrum) of each order
          result[5]: Array of the signal to noise of each order
          result[6]: Array of the normalised spectrum of each order
          result[7]: Array of the spectrum of each exposure and each order
          result[8]: Array of the merged wavelength
          result[9]: Array of the merged blaze normalised spectrum
          result[10]: Array of the error of the merged blaze normelised spectrum
          result[11]: Array of the merged normalised spectrum
          result[12]: Array of the error of the merged normalised spectrum
          result[13]: Array of the global response curve
          result[14]: Array of the correction of the small scale variations in the spectrum

For the small fits file:
          result[0]: Header of the spectrum
          result[1]: Array of the wavelength
          result[2]: Array of the spectrum

Example: