Left and middle: additional plots to Fig. 13 restricted to bound orbits recovered only by our MCMC fit, showing i) posterior distributions of semi-major axis and orbital period (left); ii) bidimensional density maps, involving the semi-major axis a versus the eccentricity e and the inclination i. The color scale appearing to the right of the plots also applies to all similar plots in Fig. 13. Right: the same bidimensional density maps as shown in Fig. 14, showing orbit solutions recovered by our LSMC fit with semi-major axes smaller than 10′′.