linestacker: a spectral line stacking tool for interferometric data

Abstract

ABSTRACT linestacker is a new open access and open source tool for stacking of spectral lines in interferometric data. linestacker is an ensemble of casa tasks, and can stack both 3D cubes or already extracted spectra. The algorithm is tested on increasingly complex simulated data sets, mimicking Atacama Large Millimeter/submillimeter Array, and Karl G. Jansky Very Large Array observations of [C ii] and CO(3–2) emission lines, from z ∼ 7 and z ∼ 4 galaxies, respectively. We find that the algorithm is very robust, successfully retrieving the input parameters of the stacked lines in all cases with an accuracy ≳90 per cent. However, we distinguish some specific situations showcasing the intrinsic limitations of the method. Mainly that high uncertainties on the redshifts (Δz > 0.01) can lead to poor signal-to-noise ratio improvement, due to lines being stacked on shifted central frequencies. Additionally, we give an extensive description of the embedded statistical tools included in linestacker: mainly bootstrapping, rebinning, and subsampling. Velocity rebinning is applied on the data before stacking and proves necessary when studying line profiles, in order to avoid artificial spectral features in the stack. Subsampling is useful to sort the stacked sources, allowing to find a subsample maximizing the searched parameters, while bootstrapping allows to detect inhomogeneities in the stacked sample. linestacker is a useful tool for extracting the most from spectral observations of various types.