Invariant-mass distributions $M(\mathrm{K}\pi\pi(\mathrm{K}^0_{\rm S}) - M(\mathrm{K}\pi\pi)$ of ${\rm D}_{\rm s1}^+$ (top panels) and ${\rm D}_{\rm s2}^+$ (bottom panels) candidates and charge conjugates in the $2 < p_{\rm T} < 24$ GeV/$c$ interval, for the INEL>0 and HMV0 (0--0.1% V0M) samples. The blue solid lines show the total fit functions described in the text, and the red dashed lines represent the combinatorial background. The raw-yield ($S$) values are reported together with their statistical uncertainties, as well as the estimated significance of the signal (Signif.). The bottom sub-panels of each panel display the differences between the data and the fit, normalized by the experimental uncertainties, to show the agreement between the data and the fit.