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 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.

A fast and portable Re-Implementation of Piskunov and Valenti's Optimal-Extraction Algorithm with improved Cosmic-Ray Removal and Optimal Sky Subtraction

A Probabilistic Method of Background Removal for High Energy Astrophysics Data

X-Sifter: detecting transients in X-ray data using the optimal Poisson matched filter

Non-negative matrix factorization approach to sky subtraction for optical spectroscopy

Sky Subtraction in an Era of Low Surface Brightness Astronomy

Proper image subtraction - optimal transient detection, photometry and hypothesis testing

BIPP: An efficient HPC implementation of the Bluebild algorithm for radio astronomy

Flat-relative optimal extraction. A quick and efficient algorithm for stabilised spectrographs

A background-estimation technique for the detection of extended gamma-ray structures with IACTs

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

cuFFS: A GPU-accelerated code for Fast Faraday Rotation Measure Synthesis

Image reduction pipeline for the detection of variable sources in highly crowded fields

CLEAN algorithm implementation comparisons between popular software packages

PulsarX: a new pulsar searching package -I. A high performance folding program for pulsar surveys

Astrophysical foreground cleanup using non-local means

Detection and Removal of Artifacts in Astronomical Images

PUMA: The Positional Update and Matching Algorithm

Image Subtraction in Fourier Space

PURIFYing real radio interferometric observations

Optimal photometry of point sources: Joint source flux and background determination on array detectors -- from theory to practical implementation