PyCosmic: a robust method to detect cosmics in CALIFA and other fiber-fed integral-field spectroscopy datasets

Abstract

Context. Detecting cosmic ray hits (cosmics) in fiber-fed integral-field spectroscopy (IFS) data of single exposures is a challenging task because of the complex signal recorded by IFS instruments. Existing detection algorithms are commonly found to be unreliable in the case of IFS data, and the optimal parameter settings are usually unknown a priori for a given dataset. Aims. The Calar Alto legacy integral field area (CALIFA) survey generates hundreds of IFS datasets for which a reliable and robust detection algorithm for cosmics is required as an important part of the fully automatic CALIFA data reduction pipeline. Such a new algorithm needs to be tested against the performance of the commonly used algorithms L.A.Cosmic and DCR. General recommendations for the usage and optimal parameter settings of each algorithm have not yet been systematically studied for fiber-fed IFS datasets to guide users in their choice. Methods. We developed a novel algorithm, PyCosmic, which combines the edge-detection algorithm of L.A.Cosmic with a pointspread function convolution scheme. We generated mock data to compute the efficiency of different algorithms for a wide range of characteristic fiber-fed IFS datasets using the Potsdam Multi-Aperture Spectrophotometer (PMAS) and the VIsible MultiObject Spectrograph (VIMOS) IFS instruments as representative cases. Results. PyCosmic is the only algorithm that achieves an acceptable detection performance for CALIFA data. We find that PyCosmic is the most robust tool with a detection rate of 90% and a false detection rate 5% for any of the tested IFS data. It has one less free parameter than the L.A.Cosmic algorithm. Only for strongly undersampled IFS data does L.A.Cosmic exceed the performance of PyCosmic by a few per cent. DCR never reaches the efficiency of the other two algorithms and should only be used if computational speed is a concern. Thus, PyCosmic appears to be the most versatile cosmics detection algorithm for IFS data. It is implemented in the new CALIFA data reduction pipeline as well as in recent versions of the multi-instrument IFS pipeline P3D. Although PyCosmic has been optimized for IFS data, we have also successfully applied it to longslit data and anticipate that good results will be achieved with imaging data.

Unraveling the Complex Structure of Agn-Driven Outflows. I. Kinematics and Sizes

“…4 First, we subtracted the CCD bias from all frames by using the standard bias frame and removed cosmic rays by using the PyCosmic routine (Husemann et al 2012), which shows robust results for cosmic-ray removal in fiber-fed IFU data. Second, we obtained the extraction solution based on the flat-field images from the afternoon calibration and corrected for the fiber-to-fiber variation by > erg s −1 and z 0.1 < from Woo et al (2016).…”

CALIFA, the Calar Alto Legacy Integral Field Area survey

“…Cosmic rays were removed and masked according to Husemann et al (2012); ii) relative offsets in the tracing due to flexure were estimated with respect to the continuum and arclamp calibration frames and the wavelength solution was applied to each individual science frame; iii) the stray-light map was reconstructed using the gaps between the fiber traces 3 and subtracted from the calibration and science exposures; iv) an optimal extraction algorithm (Horne 1986) was used to extract the spectra based on measurements of the positions and widths derived from the continuum lamp. The optimal extraction relies on a good characterization of the shape of the spectra along the cross-dispersion axis.…”

The Limited Impact of Outflows: Integral-field Spectroscopy of 20 Local AGNs

“…of ∼0.02Å. Then we removed cosmic rays on each science frame by using the PyCosmic routine, developed for a robust cosmic-ray removal of fiber-fed IFU data (Husemann et al 2012). After the cosmic-ray removal, we extracted the object and sky spectra by using the tracing solutions.…”