Volume 545, September 2012
|Number of page(s)||10|
|Section||Numerical methods and codes|
|Published online||20 September 2012|
PyCosmic: a robust method to detect cosmics in CALIFA and other fiber-fed integral-field spectroscopy datasets⋆
1 Leibniz-Institut für Astrophysik Potsdam (AIP), An der Sternwarte 16, 14482 Potsdam, Germany
2 Instituto de Astrofísica de Andalucía (CSIC), Camino Bajo de Huétor s/n Aptdo. 3004, 18080 Granada, Spain
3 Centro Astronómico Hispano Alemán de Calar Alto (CSIC-MPIA), 4004 Almería, Spain
Received: 25 July 2012
Accepted: 8 August 2012
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 point-spread 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.
Key words: techniques: image processing / instrumentation: miscellaneous
PyCosmic is freely available as a Python-based stand-alone program at http://pycosmic.sf.net for download.
© ESO, 2012
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.