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

[Abridged] Detecting cosmic ray hits (cosmics) in fiber-fed 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. The 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. We developed a novel algorithm, PyCosmic, which combines the edge-detection algorithm of L.A.Cosmic with a point-spread function convolution scheme. We generated mock data to compute the efficiency of different algorithms for a wide range of characteristic fibre-fed IFS datasets using the PMAS and VIMOS IFS instruments as representative cases. 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.

[1]  P. Weilbacher,et al.  P3D: a general data-reduction tool for fiber-fed integral-field spectrographs , 2010, 1002.4406.

[2]  K. Horne,et al.  AN OPTIMAL EXTRACTION ALGORITHM FOR CCD SPECTROSCOPY. , 1986 .

[3]  Christopher J. Conselice,et al.  The Relationship between Stellar Light Distributions of Galaxies and Their Formation Histories , 2003 .

[4]  James E. Rhoads Cosmic‐Ray Rejection by Linear Filtering of Single Images , 2000 .

[5]  Kevin A. Pimbblet,et al.  Evaluation of Cosmic Ray Rejection Algorithms on Single-Shot Exposures , 2005, Publications of the Astronomical Society of Australia.

[6]  L. Shamir,et al.  A fuzzy logic‐based algorithm for cosmic‐ray hit rejection from single images , 2005 .

[7]  P. Dokkum Cosmic-Ray Rejection by Laplacian Edge Detection , 2001, astro-ph/0108003.

[8]  H Germany,et al.  PMAS: The Potsdam Multi‐Aperture Spectrophotometer. I. Design, Manufacture, and Performance , 2005, astro-ph/0502581.

[9]  A. Hopkins,et al.  The Sydney-AAO Multi-object Integral field spectrograph , 2011, 1112.3367.

[10]  I. A. Steele,et al.  A fully automated data reduction pipeline for the FRODOSpec integral field spectrograph , 2011, 1112.2574.

[11]  W. Pych A Fast Algorithm for Cosmic‐Ray Removal from Single Images , 2003, astro-ph/0311290.

[12]  David Dunlap Observatory A Fast Algorithm for Cosmic Rays Removal from Single Images , 2003 .

[13]  Yan-mei Chen,et al.  ACTIVE GALACTIC NUCLEI WITH DOUBLE-PEAKED NARROW LINES: ARE THEY DUAL ACTIVE GALACTIC NUCLEI? , 2009, 0910.0580.

[14]  Zhongfu Ye,et al.  Cosmic-Ray Detection Based on Gray-Scale Morphology of Spectroscopic CCD Images , 2009, Publications of the Astronomical Society of Australia.

[15]  Martin M. Roth,et al.  Spatially resolved spectroscopy of planetary nebulae and their halos - I. Five galactic disk objects , 2008, 0802.3813.

[16]  M. Birchall,et al.  Optimal Extraction of Fibre Optic Spectroscopy , 2009, Publications of the Astronomical Society of Australia.

[17]  Andreas Kelz,et al.  PMAS: The Potsdam Multi-Aperture Spectrophotometer. II. The Wide Integral Field Unit PPak , 2006 .

[18]  A. Quirrenbach,et al.  CALIFA, the Calar Alto Legacy Integral Field Area survey : I. Survey presentation , 2011, 1111.0962.

[19]  Richard L. White,et al.  DECISION TREES FOR AUTOMATED IDENTIFICATION OF COSMIC-RAY HITS IN HUBBLE SPACE TELESCOPE IMAGES , 1995 .