Fit to the invariant mass distribution of $\mathrm{D_s^+}$ and $\mathrm{D^+}$ candidates in the $8< \pt <12$ GeV/$c$ interval and the $70-100\%$ FT0M percentile in pp collisions at $\sqrt{s}=13.6$ TeV, reconstructed in the $\mathrm{D_s^+, D^+ \rightarrow \phi\pi^+ \rightarrow K^+K^-\pi^+}$ decay channel. The fit function is composed of a second-degree polynomial function to describe the combinatorial background, a template fit for the $\mathrm{D^+ \rightarrow \pi^+K^-\pi^+}$ correlated background and two gaussian functions to describe the $\mathrm{D_s^+}$ and $\mathrm{D^+}$ peaks. The correlated background is described using templates taken from MC simulations, and its normalisation is fixed to the $\mathrm{D^+}$ signal integral using the relative BR and efficiencies.