Even though HARPS is probably the most precise astronomical spectrometer ever built,
the Cross-Correlation Function (CCF) data analysis method that has been commonly used
to analyze this data, is suboptimal in the sense that it does not exploit the full Doppler
information in the stellar spectrum (Pepe et al. 2002). For this reason, instead of using the
CCF RVs provided by the ESO archive, we use a least-square template-matching method
to derive new RV measurements. Thanks to the instrumental stability and the excellent
wavelength calibration provided by the HARPS-ESO data reduction software, the model
required to match each observation to a high signal-to-noise-ratio (SNR) template only needs
to include a Doppler offset and a multiplicative polynomial to correct for the flux variability
across each echelle order. The template is obtained by co-adding all the spectra after a
preliminary RV measurement is obtained using the highest SNR observation. The leastsquares
matching technique has been used on HARPS data before. An example are the RV
measurements on GJ 1214 (V=14.57) used to derive the mass of the transiting super-Earth
reported in Charbonneau et al. (2009). The performance and description of our software
tool, HARPS-TERRA (Template Enhanced Radial velocity Re-analysis Application) on a
representative sample of stars can be found in Anglada-Escud´e & Butler (2012).