If we consider the planet orbits to be circular, \(\mathbf {v}_{\rm orb}\)  from Equation REF  could then provide a complete model of the RV data as the sum
of shift (planet)-driven velocity variations and the shape (stellar activity)-driven velocity variations [1]}. Treating this as a least-square problem enabled us to solve for the unknowns, the amplitudes, and phases of \(\mathbf {F}\) . The vector \({\mathbf {\theta }_{\rm orb}}\)  can then be evaluated by minimizing the \(\chi ^2\) . Similar Implementation has been employed by [2]} as well. The detailed algorithm is given in Appendix REF .
