Fit to the invariant mass distribution of $\mathrm{D_s^+}$ and $\mathrm{D^+}$ candidates in the $2< \pt <4$ GeV/$c$ and $80-90\%$ FT0C interval in Pb$-$Pb collisions at $\sqrt{s_\mathrm{NN}}=5.36$ 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 double-sided Crystalball functions to describe the $\mathrm{D_s^+}$ and $\mathrm{D^+}$ peaks.