Comparison between runs at fixed cosmology obtained using different samplers, carried out using a reference galaxy power spectrum (Model 3 HOD sample at z = 0.9). The left panel shows the posterior distribution of the nuisance parameters with $b_{{\mathcal{G}}_2}$ fixed to the local Lagrangian relation (Eq. 35), while in the right panel all the parameters are left free to vary. Different colours correspond to different samplers, as listed in the legend.