Reading WUPPE Data

The WUPPE data files (see also Data Products) are written in FITS format. The fluxes for the ultraviolet spectrum are stored in the primary array. The wavelengths, Stokes parameters Q and U (in percent), and a combined error estimate for Q and U (also in percent), are stored in a binary table extension. If no polarization data are available, there is no binary table extension.

Various FITS readers may be used to read the data files. The IUEDAC software provides general FITS file reading tools written in IDL, such as ifitsrd, which can be used to read the WUPPE data files. The header, extension header, flux, wavelengths, Stokes parameters U and Q, and an error estimate are returned. The Stokes parameters must be converted to provide polarization and position angle as a function of wavelength. For instance, to read the file ab-aur_312011_2_hw.fits:

ifitsrd,'ab-aur_312011_2_hw.fits',-1,header,exthead,flux
ifitsrd,'ab-aur_312011_2_hw.fits',1,header,exthead,wave,q,u,err

The linear polarization and position angle may then be computed from Q and U. Since the data are typically noisy, the values are generally binned.

polarization = sqrt(q*q + u*u)
posangle = 0.5 * 57.2958 * atan(u/q)

For data with flux data but no polarization data, one can create the wavelength array from the FITS keywords CRVAL1 (starting wavelength in Ångstroms) and CRDELT1 (wavelength array increment).

ifitsrd,'ngc4151_811623_2_hw.fits',-1,header,exthead,flux
stpar,header,'crval1',crval1
stpar,header,'cdelt1',cdelt1
print,crval1,cdelt1           (yields 1366.00, 2.0000)
wave = fltarr(1024) + crval1
for i=1,1023 do wave(i) = wave(i-1) + cdelt1
plot,wave,flux

The error array is set to zero where data are bad quality, due to data dropouts or other problems.

The WUPPE text and Postscript files may be read or displayed using standard tools.