Invariant mass distribution of J/$\psi$ candidates in pp collisions at $\sqrt{s}=$ 13 TeV for 8 $< p_{\rm T}< $ 15 GeV/$c$ (left) and 30 $< p_{\rm T}< $ 40 GeV/$c$ (right). The gray open markers depict the distribution for e$^+$e$^-$ pairs where at least one track could be matched to an EMCal cluster. It is scaled by $1/2$ and $1/4$, respectively, for the different $p_{\rm T}$ intervals, to enhance the visibility. The black closed markers represent the distribution after applying the EMCal PID selections on at least one of the J/$\psi$ decay products. The combined signal and background fit is shown by the blue line, while the polynomial background and pure signal fit are shown as gray and red lines respectively [75].