Multi-mask least-squares deconvolution: extracting RVs using tailored masks

To push the radial velocity (RV) exoplanet detection threshold, it is crucial to find more reliable RV extraction methods. The least-squares deconvolution (LSD) technique has been used to infer the stellar magnetic flux from spectropolarimetric data for the past two decades. It relies on the assumption that stellar absorption lines are similar in shape. Although this assumption is simplistic, LSD provides a good model for intensity spectra and likewise an estimate for their Doppler shift. We present the multi-mask least-squares deconvolution (MM-LSD) RV extraction pipeline that extracts the RV from two-dimensional echelle-order spectra using LSD with multiple tailored masks after continuum normalization and telluric absorption line correction. The flexibility of LSD allows to exclude spectral lines or pixels at will, providing a means to exclude variable lines or pixels affected by instrumental problems. The MM-LSD pipeline was tested on HARPS-N data for the Sun and selected well-observed stars with 5.7 < V_mag < 12.6. For FGK-type stars with median signal-to-noise ratio above 100, the pipeline delivered RV time series with on average 12 per cent lower scatter as compared to the HARPS-N RV extraction pipeline based on the cross-correlation function technique. The MM-LSD pipeline may be used as a standalone RV code, or modified and extended to extract a proxy for the magnetic field strength.