Time series of the RVs of the first observing season (top panel) and second observing season (bottom panel) displayed in red together with the best-fit Keplerian model (in black) and GP-fitted model used to account for stellar activity (in blue). The RVs were corrected by the systemic velocity and the RV trend. We show the 1 σ uncertainties from the GP model in shaded blue. We also show the residuals of the best model below the RVs.