Upper limits at 95% CL on $\sigma$B(pp→X→Y(bb)H) for combination as a function of m$_Y$.

Expected upper limits at 95% CL on the X → HH cross section as a function of mX from the combination of the three channels bbγγ, bbττ, and bbbb for an integrated luminosity of 3000 fb−1 with different assumptions on the systematic uncertainties.

Expected upper limits on the X → YH cross section times B(Y → bb) for mX ≤ 1000 GeV at 95 % CL as a function of m$_X$ and m$_Y$ from the combination of all considered channels projected to an integrated luminosity of 3000 fb$^{-1}$. Systematic uncertainties according to the ”S2” scenario are assumed.

Example description

Example description

Example description

Expected and observed 95% CL upper limits on $|V_\mathrm{N}|^2$ as a function of $m_\mathrm{N}$ for the mixing scenario ($r_e$, $r_\mu$, $r_\tau$) = (0.33, 0.33, 0.33) and in the Majorana scenario.

Expected and observed 95% CL upper limits on $|V_\mathrm{N}|^2$ as a function of $m_\mathrm{N}$ for the mixing scenario ($r_e$, $r_\mu$, $r_\tau$) = (0.0, 1.0, 0.0) and in the Dirac-like scenario.

Expected and observed 95% CL upper limits on $|V_\mathrm{N}|^2$ as a function of $m_\mathrm{N}$ for the mixing scenario ($r_e$, $r_\mu$, $r_\tau$) = (0.0, 0.5, 0.5) and in the Dirac-like scenario.

Expected and observed 95% CL upper limits on $|V_\mathrm{N}|^2$ as a function of $m_\mathrm{N}$ for the mixing scenario ($r_e$, $r_\mu$, $r_\tau$) = (0.5, 0.5, 0.0) and in the Dirac-like scenario.

Expected and observed 95% CL upper limits on $|V_\mathrm{N}|^2$ as a function of $m_\mathrm{N}$ for the mixing scenario ($r_e$, $r_\mu$, $r_\tau$) = (0.33, 0.33, 0.33) and in the Dirac-like scenario.

Observed 95% CL lower limits on $c\tau_\mathrm{N}$ as functions of the mixing ratios ($r_e$, $r_\mu$, $r_\tau$) for a fixed N mass of 1.0 GeV in the Majorana scenario.

Observed 95% CL lower limits on $c\tau_\mathrm{N}$ as functions of the mixing ratios ($r_e$, $r_\mu$, $r_\tau$) for a fixed N mass of 1.5 GeV in the Majorana scenario.

Observed 95% CL lower limits on $c\tau_\mathrm{N}$ as functions of the mixing ratios ($r_e$, $r_\mu$, $r_\tau$) for a fixed N mass of 2.0 GeV in the Majorana scenario.

Observed 95% CL lower limits on $c\tau_\mathrm{N}$ as functions of the mixing ratios ($r_e$, $r_\mu$, $r_\tau$) for a fixed N mass of 1.0 GeV in the Dirac-like scenario.

Observed 95% CL lower limits on $c\tau_\mathrm{N}$ as functions of the mixing ratios ($r_e$, $r_\mu$, $r_\tau$) for a fixed N mass of 1.5 GeV in the Dirac-like scenario.

Observed 95% CL lower limits on $c\tau_\mathrm{N}$ as functions of the mixing ratios ($r_e$, $r_\mu$, $r_\tau$) for a fixed N mass of 2.0 GeV in the Dirac-like scenario.

Measured natural width

The measured inclusive ratio of production cross sections

Distributions of $S_{\mathrm{ML}}$ for data, simulated background and signal events with $n_{\mathrm{track}}$ of 4. The distributions are shown for split-SUSY signals with a gluino mass of 2000 GeV and neutralino mass of 1800 GeV. Different gluino proper decay lengths are shown as $c\tau$ in the legend. All distributions are normalized to unity.

Distributions of $S_{\mathrm{ML}}$ for data, simulated background and signal events with $n_{\mathrm{track}}$ of $\geq$ 5. The distributions are shown for split-SUSY signals with a gluino mass of 2000 GeV and neutralino mass of 1900 GeV. Different gluino proper decay lengths are shown as $c\tau$ in the legend. All distributions are normalized to unity.

Distributions of $S_{\mathrm{ML}}$ for data, simulated background and signal events with $n_{\mathrm{track}}$ of $\geq$ 5. The distributions are shown for split-SUSY signals with a gluino mass of 2000 GeV and neutralino mass of 1800 GeV. Different gluino proper decay lengths are shown as $c\tau$ in the legend. All distributions are normalized to unity.

The distribution of $n_{\mathrm{track}}$ in different $S_{\mathrm{ML}}$ regions for simulated background events. Events with 0 $ < S_{\mathrm{ML}} < $ 0.2 (blue), 0.2 $ < S_{\mathrm{ML}} < $ 0.6 (red), and 0.6 $ < S_{\mathrm{ML}} < $ 1.0 (green) are compared. All distributions are normalized to unity. The similar $n_{\mathrm{track}}$ distributions demonstrate that $n_{\mathrm{track}}$ and $S_{\mathrm{ML}}$ are decorrelated.

The distribution of $d_{\mathrm{BV}}$ in $K_{\mathrm{S}}^{\mathrm{0}}$ vertices in data (black) and simulation (purple). The lower panel shows the ratio between data and simulation.

The vertex reconstruction efficiency for artificially displaced vertices in data (black) and simulation (red). In this example, the artificially displaced vertices are corrected to mimic split-SUSY signal events with gluino mass of 2000 GeV and neutralino mass of 1800 GeV.

The ML tagging efficiency for artificially displaced vertices in data (black) and simulation (red). In this example, the artificially displaced vertices are corrected to mimic split-SUSY signal events with gluino mass of 2000 GeV and neutralino mass of 1800 GeV.

Number of predicted and observed events in the control, validation, and search regions. Predictions are calculated using Eqs. (2) and (3) and fitting the data under the background-only hypothesis. Regions are organized by $S_{\mathrm{ML}}$ and $n_{\mathrm{track}}$ values, and region names corresponding with Fig. 7 are given in parentheses. The predicted number of events that pass the $S_{\mathrm{ML}}$ selection and the observed number of events that pass or fail the $S_{\mathrm{ML}}$ selection are shown in seperate rows.

The 95% CL upper limit on the product of the cross section and branching fraction squared for the split-SUSY signal model with a mass splitting of 100 GeV, shown as a function of gluino mass and $c\tau$. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the limit plot.

The 95% CL upper limit on the product of the cross section and branching fraction squared for the split-SUSY signal model with a mass splitting of 100 GeV, shown as a function of gluino mass and $c\tau$. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the limit plot.

The 95% CL upper limit on the product of the cross section and branching fraction squared for the split-SUSY model with a $c\tau$ of 10 mm, shown as a function of gluino mass and mass splitting. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the limit plot.

The 95% CL upper limit on the product of the cross section and branching fraction squared for the split-SUSY model with a $c\tau$ of 10 mm, shown as a function of gluino mass and mass splitting. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the limit plot.

The 95% CL upper limit on the product of the cross section and branching fraction squared for the GMSB SUSY signal model, shown as a function of gluino mass and $c\tau$. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the limit plot.

The 95% CL upper limit on the product of the cross section and branching fraction squared for the GMSB SUSY signal model, shown as a function of gluino mass and $c\tau$. The observed (solid black) and expected (dashed red) exclusion curves are overlaid on the limit plot.

Comparison of the observed (black points) and expected (histograms) numbers of events in nonoverlapping $m^{corr}_{\mu\mu}$ intervals in the STA-STA dimuon category, in the signal region optimized for the RPV SUSY model. Yellow and green stacked filled histograms represent mean expected background contributions from QCD and DY, respectively, while statistical uncertainties in the total expected background are shown as hatched histograms. Signal contributions expected from simulated signals indicated in the legends are shown in red and blue. Their yields are set to the corresponding median expected 95% CL exclusion limits obtained from the ensemble of both dimuon categories, scaled up as indicated in the legend to improve visibility. The last bin includes events in the histogram overflow.

Distribution of min($d_0 / \sigma_{d_0}$) for TMS-TMS dimuons with $|\Delta\Phi| < \pi/30$, for events in all mass intervals combined, for both the validation (min($d_0 / \sigma_{d_0}$) < 6) and signal (min($d_0 / \sigma_{d_0}$) > 6) regions. The number of observed events (black circles) is overlaid with the stacked histograms showing the expected numbers of QCD (yellow) and DY (green) background events. Statistical uncertainties in the total expected background are shown as hatched histograms. Signal contributions expected from simulated signals indicated in the legends are shown in red and blue. Their yields are set to the corresponding median expected 95% CL exclusion limits obtained from the ensemble of both dimuon categories, scaled up as indicated in the legend to improve visibility. Events are required to satisfy all nominal selection criteria with the exception of the $d_0 / \sigma_{d_0}$ requirement. The last bin includes events in the histogram overflow.

Distribution of min($d_0 / \sigma_{d_0}$) for TMS-TMS dimuons with $|\Delta\Phi| < \pi/4$, for events in all mass intervals combined, for both the validation (min($d_0 / \sigma_{d_0}$) < 6) and signal (min($d_0 / \sigma_{d_0}$) > 6) regions. The number of observed events (black circles) is overlaid with the stacked histograms showing the expected numbers of QCD (yellow) and DY (green) background events. Statistical uncertainties in the total expected background are shown as hatched histograms. Signal contributions expected from simulated signals indicated in the legends are shown in red and blue. Their yields are set to the corresponding median expected 95% CL exclusion limits obtained from the ensemble of both dimuon categories, scaled up as indicated in the legend to improve visibility. Events are required to satisfy all nominal selection criteria with the exception of the $d_0 / \sigma_{d_0}$ requirement. The last bin includes events in the histogram overflow.

Comparison of observed and expected numbers of events in bins of $m^{corr}_{\mu\mu}$ in the TMS-TMS dimuon category, in the signal regions optimized for the RPV SUSY model. The number of observed events (black circles) is overlaid with the stacked filled histograms showing the expected numbers of QCD (yellow) and DY (green) background events in bins of $m^{corr}_{\mu\mu}$ in min($d_0 / \sigma_{d_0}$) bin: 6-10. Hatched histograms show statistical uncertainties in the total expected background. Contributions expected from signal events predicted by the RPV SUSY model with the parameters indicated in the legends are shown as red and blue histograms. Their yields are set to the corresponding median expected 95% CL exclusion limits obtained from the ensemble of both dimuon categories, scaled up as indicated in the legend to improve visibility. The last bin includes events in the histogram overflow.

Comparison of observed and expected numbers of events in bins of $m^{corr}_{\mu\mu}$ in the TMS-TMS dimuon category, in the signal regions optimized for the RPV SUSY model. The number of observed events (black circles) is overlaid with the stacked filled histograms showing the expected numbers of QCD (yellow) and DY (green) background events in bins of $m^{corr}_{\mu\mu}$ in min($d_0 / \sigma_{d_0}$) bin: 10-20. Hatched histograms show statistical uncertainties in the total expected background. Contributions expected from signal events predicted by the RPV SUSY model with the parameters indicated in the legends are shown as red and blue histograms. Their yields are set to the corresponding median expected 95% CL exclusion limits obtained from the ensemble of both dimuon categories, scaled up as indicated in the legend to improve visibility. The last bin includes events in the histogram overflow.

Comparison of observed and expected numbers of events in bins of $m^{corr}_{\mu\mu}$ in the TMS-TMS dimuon category, in the signal regions optimized for the RPV SUSY model. The number of observed events (black circles) is overlaid with the stacked filled histograms showing the expected numbers of QCD (yellow) and DY (green) background events in bins of $m^{corr}_{\mu\mu}$ in min($d_0 / \sigma_{d_0}$) bin: >20. Hatched histograms show statistical uncertainties in the total expected background. Contributions expected from signal events predicted by the RPV SUSY model with the parameters indicated in the legends are shown as red and blue histograms. Their yields are set to the corresponding median expected 95% CL exclusion limits obtained from the ensemble of both dimuon categories, scaled up as indicated in the legend to improve visibility. The last bin includes events in the histogram overflow.

Comparison of observed and expected numbers of events in bins of $m_{\mu\mu}$ in the TMS-TMS dimuon category, in the signal regions optimized for the HAHM. The number of observed events (black circles) is overlaid with the stacked filled histograms showing the expected numbers of QCD (yellow) and DY (green) background events in bins of $m_{\mu\mu}$ in min($d_0 / \sigma_{d_0}$) bin: 6-10. Hatched histograms show statistical uncertainties in the total expected background. Signal contributions expected from simulated $H \rightarrow Z_DZ_D$ events with the parameters indicated in the legends are shown as red and blue histograms. Their yields are set to the corresponding median expected 95% CL exclusion limits obtained from the ensemble of both dimuon categories, scaled up as indicated in the legend to improve visibility. The last bin includes events in the histogram overflow.

Comparison of observed and expected numbers of events in bins of $m_{\mu\mu}$ in the TMS-TMS dimuon category, in the signal regions optimized for the HAHM. The number of observed events (black circles) is overlaid with the stacked filled histograms showing the expected numbers of QCD (yellow) and DY (green) background events in bins of $m_{\mu\mu}$ in min($d_0 / \sigma_{d_0}$) bin: 10-20. Hatched histograms show statistical uncertainties in the total expected background. Signal contributions expected from simulated $H \rightarrow Z_DZ_D$ events with the parameters indicated in the legends are shown as red and blue histograms. Their yields are set to the corresponding median expected 95% CL exclusion limits obtained from the ensemble of both dimuon categories, scaled up as indicated in the legend to improve visibility. The last bin includes events in the histogram overflow.

Comparison of observed and expected numbers of events in bins of $m_{\mu\mu}$ in the TMS-TMS dimuon category, in the signal regions optimized for the HAHM. The number of observed events (black circles) is overlaid with the stacked filled histograms showing the expected numbers of QCD (yellow) and DY (green) background events in bins of $m_{\mu\mu}$ in min($d_0 / \sigma_{d_0}$) bin: >20. Hatched histograms show statistical uncertainties in the total expected background. Signal contributions expected from simulated $H \rightarrow Z_DZ_D$ events with the parameters indicated in the legends are shown as red and blue histograms. Their yields are set to the corresponding median expected 95% CL exclusion limits obtained from the ensemble of both dimuon categories, scaled up as indicated in the legend to improve visibility. The last bin includes events in the histogram overflow.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m(Z_D) = 10\ GeV$, in the STA-STA and TMS-TMS dimuon categories in 2022 data and their combination.The median expected limits obtained from the STA-STA and TMS-TMS dimuon categories are shown as dashed blue and red curves, respectively; the combined median expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m(Z_D) = 20\ GeV$, in the STA-STA and TMS-TMS dimuon categories in 2022 data and their combination.The median expected limits obtained from the STA-STA and TMS-TMS dimuon categories are shown as dashed blue and red curves, respectively; the combined median expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m(Z_D) = 30\ GeV$, in the STA-STA and TMS-TMS dimuon categories in 2022 data and their combination.The median expected limits obtained from the STA-STA and TMS-TMS dimuon categories are shown as dashed blue and red curves, respectively; the combined median expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m(Z_D) = 40\ GeV$, in the STA-STA and TMS-TMS dimuon categories in 2022 data and their combination.The median expected limits obtained from the STA-STA and TMS-TMS dimuon categories are shown as dashed blue and red curves, respectively; the combined median expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m(Z_D) = 50\ GeV$, in the STA-STA and TMS-TMS dimuon categories in 2022 data and their combination.The median expected limits obtained from the STA-STA and TMS-TMS dimuon categories are shown as dashed blue and red curves, respectively; the combined median expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m(Z_D) = 60\ GeV$, in the STA-STA and TMS-TMS dimuon categories in 2022 data and their combination.The median expected limits obtained from the STA-STA and TMS-TMS dimuon categories are shown as dashed blue and red curves, respectively; the combined median expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m(Z_D) = 10\ GeV$, obtained in this analysis, the Run 2 analysis, and their combination. The observed limits in this analysis and in the Run 2 analysis are shown as blue and red curves, respectively; the median combined expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m(Z_D) = 20\ GeV$, obtained in this analysis, the Run 2 analysis, and their combination. The observed limits in this analysis and in the Run 2 analysis are shown as blue and red curves, respectively; the median combined expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m(Z_D) = 30\ GeV$, obtained in this analysis, the Run 2 analysis, and their combination. The observed limits in this analysis and in the Run 2 analysis are shown as blue and red curves, respectively; the median combined expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m(Z_D) = 40\ GeV$, obtained in this analysis, the Run 2 analysis, and their combination. The observed limits in this analysis and in the Run 2 analysis are shown as blue and red curves, respectively; the median combined expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m(Z_D) = 50\ GeV$, obtained in this analysis, the Run 2 analysis, and their combination. The observed limits in this analysis and in the Run 2 analysis are shown as blue and red curves, respectively; the median combined expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $B(H \rightarrow Z_DZ_D)$ as a function of $c\tau(Z_D)$ in the HAHM model, for $m(Z_D) = 60\ GeV$, obtained in this analysis, the Run 2 analysis, and their combination. The observed limits in this analysis and in the Run 2 analysis are shown as blue and red curves, respectively; the median combined expected limits are shown as dashed black curves; and the combined observed limits are shown as solid black curves. The green and yellow bands correspond, respectively, to the 68 and 95% quantiles for the combined expected limits.

The 95% CL upper limits on $\sigma(pp \rightarrow \tilde{q}\bar{\tilde{q}})B(\tilde{q} \rightarrow q\tilde{\chi}^{0}_{1})$ as a function of $c\tau(\tilde{\chi}^{0}_{1})$ in the RPV SUSY model, for $B(\tilde{\chi}^{0}_{1} \rightarrow \mu^{+}\mu^{-}\nu) = 0.5$ and $m(\tilde{q}) = 125\ GeV$. The observed limits for various $m(\tilde{\chi}^{0}_{1})$ indicated in the legends are shown as solid curves. The median expected limits and their 68 and 95% quantiles are shown, respectively, as dashed black curves and green and yellow bands for the case of $m(\tilde{\chi}^{0}_{1}) = 50\ GeV$ and omitted for other neutralino masses for clarity. The predicted cross section for $m(\tilde{q}) = 125\ GeV$ is 7200 pb, and falls outside the y-axis range.

The 95% CL upper limits on $\sigma(pp \rightarrow \tilde{q}\bar{\tilde{q}})B(\tilde{q} \rightarrow q\tilde{\chi}^{0}_{1})$ as a function of $c\tau(\tilde{\chi}^{0}_{1})$ in the RPV SUSY model, for $B(\tilde{\chi}^{0}_{1} \rightarrow \mu^{+}\mu^{-}\nu) = 0.5$ and $m(\tilde{q}) = 200\ GeV$. The observed limits for various $m(\tilde{\chi}^{0}_{1})$ indicated in the legends are shown as solid curves. The median expected limits and their 68 and 95% quantiles are shown, respectively, as dashed black curves and green and yellow bands for the case of $m(\tilde{\chi}^{0}_{1}) = 50\ GeV$ and omitted for other neutralino masses for clarity. The predicted cross section for $m(\tilde{q}) = 200 GeV$ is 840 pb, and falls outside the y-axis range.

The 95% CL upper limits on $\sigma(pp \rightarrow \tilde{q}\bar{\tilde{q}})B(\tilde{q} \rightarrow q\tilde{\chi}^{0}_{1})$ as a function of $c\tau(\tilde{\chi}^{0}_{1})$ in the RPV SUSY model, for $B(\tilde{\chi}^{0}_{1} \rightarrow \mu^{+}\mu^{-}\nu) = 0.5$ and $m(\tilde{q}) = 350\ GeV$. The observed limits for various $m(\tilde{\chi}^{0}_{1})$ indicated in the legends are shown as solid curves. The median expected limits and their 68 and 95% quantiles are shown, respectively, as dashed black curves and green and yellow bands for the case of $m(\tilde{\chi}^{0}_{1}) = 50\ GeV$ and omitted for other neutralino masses for clarity. The predicted cross section for $m(\tilde{q}) = 350\ GeV$ is 50 pb, and falls outside the y-axis range.

The 95% CL upper limits on $\sigma(pp \rightarrow \tilde{q}\bar{\tilde{q}})B(\tilde{q} \rightarrow q\tilde{\chi}^{0}_{1})$ as a function of $c\tau(\tilde{\chi}^{0}_{1})$ in the RPV SUSY model, for $B(\tilde{\chi}^{0}_{1} \rightarrow \mu^{+}\mu^{-}\nu) = 0.5$ and $m(\tilde{q}) = 700\ GeV$. The observed limits for various $m(\tilde{\chi}^{0}_{1})$ indicated in the legends are shown as solid curves. The median expected limits and their 68 and 95% quantiles are shown, respectively, as dashed black curves and green and yellow bands for the case of $m(\tilde{\chi}^{0}_{1}) = 50\ GeV$ and omitted for other neutralino masses for clarity. The gray horizontal line indicates the theoretical value of the squark-antisquark production cross section with the uncertainties shown as the gray shaded band.

The 95% CL upper limits on $\sigma(pp \rightarrow \tilde{q}\bar{\tilde{q}})B(\tilde{q} \rightarrow q\tilde{\chi}^{0}_{1})$ as a function of $c\tau(\tilde{\chi}^{0}_{1})$ in the RPV SUSY model, for $B(\tilde{\chi}^{0}_{1} \rightarrow \mu^{+}\mu^{-}\nu) = 0.5$ and $m(\tilde{q}) = 1150\ GeV$. The observed limits for various $m(\tilde{\chi}^{0}_{1})$ indicated in the legends are shown as solid curves. The median expected limits and their 68 and 95% quantiles are shown, respectively, as dashed black curves and green and yellow bands for the case of $m(\tilde{\chi}^{0}_{1}) = 50\ GeV$ and omitted for other neutralino masses for clarity. The gray horizontal line indicates the theoretical value of the squark-antisquark production cross section with the uncertainties shown as the gray shaded band.

The 95% CL upper limits on $\sigma(pp \rightarrow \tilde{q}\bar{\tilde{q}})B(\tilde{q} \rightarrow q\tilde{\chi}^{0}_{1})$ as a function of $c\tau(\tilde{\chi}^{0}_{1})$ in the RPV SUSY model, for $B(\tilde{\chi}^{0}_{1} \rightarrow \mu^{+}\mu^{-}\nu) = 0.5$ and $m(\tilde{q}) = 1600\ GeV$. The observed limits for various $m(\tilde{\chi}^{0}_{1})$ indicated in the legends are shown as solid curves. The median expected limits and their 68 and 95% quantiles are shown, respectively, as dashed black curves and green and yellow bands for the case of $m(\tilde{\chi}^{0}_{1}) = 50\ GeV$ and omitted for other neutralino masses for clarity. The gray horizontal line indicates the theoretical value of the squark-antisquark production cross section with the uncertainties shown as the gray shaded band.

Fractions of signal events with zero (green), one (blue), and two (red) STA muons matched to TMS muons by the STA to TMS association procedure, as a function of generated $L_{xy}$, in all HAHM signal samples combined.

Efficiencies of the Run 2 and Run 3 displaced dimuon triggers as a function of $c\tau$ for the HAHM signal events with $m(Z_D) = 50\ GeV$. The efficiency is defined as the fraction of simulated events that satisfy the requirements of the following sets of trigger paths: the Run 2 (2018) triggers (dashed black); the Run 3 (2022, L3) triggers (blue); the Run 3 (2022, L2) triggers (red); and the OR of all these triggers (Run 3 (2022), black). The lower panel shows the ratio of the overall Run 3 (2022) efficiency to the Run 2 (2018) efficiency.

Efficiencies of the Run 2 (2018) (red) and Run 3 (2022) (black) sets of displaced dimuon triggers as a function of $m(Z_D)$ for the HAHM signal events with $c\tau = 1\ cm$. The efficiency is defined as the fraction of simulated events that satisfy the detector acceptance and the requirements of the indicated set of trigger paths. The lower panel shows the ratio of the Run 3 (2022) efficiency to the Run 2 (2018) efficiency.

Efficiencies of the Run 2 (2018) (red) and Run 3 (2022) (black) sets of displaced dimuon triggers as a function of $m(Z_D)$ for the HAHM signal events with $c\tau = 10\ m$. The efficiency is defined as the fraction of simulated events that satisfy the detector acceptance and the requirements of the indicated set of trigger paths. The lower panel shows the ratio of the Run 3 (2022) efficiency to the Run 2 (2018) efficiency.

Overall selection efficiencies as a function of $c\tau(Z_D)$ for the HAHM signal with $m(Z_D) = 20\ GeV$ in different years of data taking. Efficiencies are computed as the ratios of the number of simulated signal events in which at least one dimuon candidate passes all 2016 (dashed green), 2018 (dashed red), and 2022 (solid black) trigger and offline selection criteria to the total number of simulated signal events. The lower panel shows the ratio of the 2022 efficiency to the 2018 efficiency (dashed red) and to the 2016 efficiency (dashed green).

Overall selection efficiencies as a function of $c\tau(Z_D)$ for the HAHM signal with $m(Z_D) = 50\ GeV$ in different years of data taking. Efficiencies are computed as the ratios of the number of simulated signal events in which at least one dimuon candidate passes all 2016 (dashed green), 2018 (dashed red), and 2022 (solid black) trigger and offline selection criteria to the total number of simulated signal events. The lower panel shows the ratio of the 2022 efficiency to the 2018 efficiency (dashed red) and to the 2016 efficiency (dashed green).

Overall selection efficiencies as a function of $c\tau(Z_D)$ for the HAHM model with $m(Z_D) = 10\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Overall selection efficiencies as a function of $c\tau(Z_D)$ for the HAHM model with $m(Z_D) = 20\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Overall selection efficiencies as a function of $c\tau(Z_D)$ for the HAHM model with $m(Z_D) = 30\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Overall selection efficiencies as a function of $c\tau(Z_D)$ for the HAHM model with $m(Z_D) = 40\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Overall selection efficiencies as a function of $c\tau(Z_D)$ for the HAHM model with $m(Z_D) = 50\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Overall selection efficiencies as a function of $c\tau(Z_D)$ for the HAHM model with $m(Z_D) = 60\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Overall selection efficiencies as a function of $c\tau(\tilde{\chi}^{0}_{1})$ for the RPV SUSY model, for events with $m(\tilde{q}) = 125\ GeV$ and $m(\tilde{\chi}^{0}_{1}) = 50\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Overall selection efficiencies as a function of $c\tau(\tilde{\chi}^{0}_{1})$ for the RPV SUSY model, for events with $m(\tilde{q}) = 200\ GeV$ and $m(\tilde{\chi}^{0}_{1}) = 50\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Overall selection efficiencies as a function of $c\tau(\tilde{\chi}^{0}_{1})$ for the RPV SUSY model, for events with $m(\tilde{q}) = 350\ GeV$ and $m(\tilde{\chi}^{0}_{1}) = 50\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Overall selection efficiencies as a function of $c\tau(\tilde{\chi}^{0}_{1})$ for the RPV SUSY model, for events with $m(\tilde{q}) = 700\ GeV$ and $m(\tilde{\chi}^{0}_{1}) = 50\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Overall selection efficiencies as a function of $c\tau(\tilde{\chi}^{0}_{1})$ for the RPV SUSY model, for events with $m(\tilde{q}) = 1150\ GeV$ and $m(\tilde{\chi}^{0}_{1}) = 50\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Overall selection efficiencies as a function of $c\tau(\tilde{\chi}^{0}_{1})$ for the RPV SUSY model, for events with $m(\tilde{q}) = 1600\ GeV$ and $m(\tilde{\chi}^{0}_{1}) = 50\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Overall selection efficiencies as a function of $c\tau(\tilde{\chi}^{0}_{1})$ for the RPV SUSY model, for events with $m(\tilde{q}) = 700\ GeV$ and $m(\tilde{\chi}^{0}_{1}) = 500\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Overall selection efficiencies as a function of $c\tau(\tilde{\chi}^{0}_{1})$ for the RPV SUSY model, for events with $m(\tilde{q}) = 1150\ GeV$ and $m(\tilde{\chi}^{0}_{1}) = 500\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Overall selection efficiencies as a function of $c\tau(\tilde{\chi}^{0}_{1})$ for the RPV SUSY model, for events with $m(\tilde{q}) = 1600\ GeV$ and $m(\tilde{\chi}^{0}_{1}) = 500\ GeV$. The plot shows efficiencies of the two dimuon categories, TMS-TMS (dashed red) and STA-STA (dashed green), as well as their combination (solid black). Each efficiency is computed as the ratio of the number of simulated signal events in which at least one dimuon candidate of a given type (or any type for the combined efficiency) passes all selection criteria (including the trigger) to the total number of simulated signal events. All efficiencies are corrected by the data-to-simulation scale factors described in the paper.

Signal efficiencies in the TMS-TMS dimuon category as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $L_{xy}^\mathrm{true}$ smaller than 20 cm in the HAHM signal model. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper.

Signal efficiencies in the STA-STA dimuon category as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $L_{xy}^\mathrm{true}$ smaller than 20 cm in the HAHM signal model. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper.

Signal efficiencies in the TMS-TMS dimuon category as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $L_{xy}^\mathrm{true}$ 20-70 cm in the HAHM signal model. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper.

Signal efficiencies in the STA-STA dimuon category as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $L_{xy}^\mathrm{true}$ 20-70 cm in the HAHM signal model. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper.

Signal efficiencies in the STA-STA dimuon category as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $L_{xy}^\mathrm{true}$ 70-500 cm in the HAHM signal model. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper.

Signal efficiencies in the TMS-TMS dimuon category as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $L_{xy}^\mathrm{true}$ smaller than 20 cm in the RPV SUSY signal model. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper.

Signal efficiencies in the STA-STA dimuon category as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $L_{xy}^\mathrm{true}$ smaller than 20 cm in the RPV SUSY signal model. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper.

Signal efficiencies in the TMS-TMS dimuon category as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $L_{xy}^\mathrm{true}$ 20-70 cm in the RPV SUSY signal model. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper.

Signal efficiencies in the STA-STA dimuon category as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $L_{xy}^\mathrm{true}$ 20-70 cm in the RPV SUSY signal model. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper.

Signal efficiencies in the STA-STA dimuon category as a function of the smaller of the two values of generated muon $p_T$ and $d_0$ in dimuons with $L_{xy}^\mathrm{true}$ 70-500 cm in the RPV SUSY signal model. The efficiency in each bin is computed as the ratio of the number of simulated signal dimuons in that bin that pass the trigger requirements and selection criteria to the total number of simulated signal dimuons in that bin and within the geometric acceptance. The geometric acceptance is defined as the generated longitudinal decay length $L_{z}$ smaller than $8\ m$ and $|\eta^\mathrm{true}|$ of both generated muons forming the dimuon smaller than 2.0. The efficiencies obtained from simulation were further corrected by the data-to-simulation scale factors described in the paper.

Unfolded E2C distributions in data compared to MC predictions.

Unfolded E2C distributions in data compared to MC predictions.

Unfolded E2C distributions in data compared to MC predictions.

Unfolded E2C distributions in data compared to MC predictions.

Unfolded E2C distributions in data compared to MC predictions.

Unfolded E3C distributions in data compared to MC predictions.

Unfolded E3C distributions in data compared to MC predictions.

Unfolded E3C distributions in data compared to MC predictions.

Unfolded E3C distributions in data compared to MC predictions.

Unfolded E3C distributions in data compared to MC predictions.

Unfolded E3C distributions in data compared to MC predictions.

Unfolded E3C distributions in data compared to MC predictions.

Unfolded E3C distributions in data compared to MC predictions.

Unfolded E3C/E2C distributions in data compared to NLO+NNLL_approx predictions. Theoretial uncertainty in each bin is handled fully correlated as shape uncertainty. For the one-sided uncertainties, we symmetrize them when performing the fit.

Unfolded E3C/E2C distributions in data compared to NLO+NNLL_approx predictions. Theoretial uncertainty in each bin is handled fully correlated as shape uncertainty. For the one-sided uncertainties, we symmetrize them when performing the fit.

Unfolded E3C/E2C distributions in data compared to NLO+NNLL_approx predictions. Theoretial uncertainty in each bin is handled fully correlated as shape uncertainty. For the one-sided uncertainties, we symmetrize them when performing the fit.

Unfolded E3C/E2C distributions in data compared to NLO+NNLL_approx predictions. Theoretial uncertainty in each bin is handled fully correlated as shape uncertainty. For the one-sided uncertainties, we symmetrize them when performing the fit.

Unfolded E3C/E2C distributions in data compared to NLO+NNLL_approx predictions. Theoretial uncertainty in each bin is handled fully correlated as shape uncertainty. For the one-sided uncertainties, we symmetrize them when performing the fit.

Unfolded E3C/E2C distributions in data compared to NLO+NNLL_approx predictions. Theoretial uncertainty in each bin is handled fully correlated as shape uncertainty. For the one-sided uncertainties, we symmetrize them when performing the fit.

Unfolded E3C/E2C distributions in data compared to NLO+NNLL_approx predictions. Theoretial uncertainty in each bin is handled fully correlated as shape uncertainty. For the one-sided uncertainties, we symmetrize them when performing the fit.

Unfolded E3C/E2C distributions in data compared to NLO+NNLL_approx predictions. Theoretial uncertainty in each bin is handled fully correlated as shape uncertainty. For the one-sided uncertainties, we symmetrize them when performing the fit.

The fitted slopes of the E3C/E2C data distributions as a function of jet pt are used to illustrate the dependency of alphas on jet pt.

Unfolded E3C/E2C distributions in data compared to MC predictions.

Unfolded E3C/E2C distributions in data compared to MC predictions.

Unfolded E3C/E2C distributions in data compared to MC predictions.

Unfolded E3C/E2C distributions in data compared to MC predictions.

Unfolded E3C/E2C distributions in data compared to MC predictions.

Unfolded E3C/E2C distributions in data compared to MC predictions.

Unfolded E3C/E2C distributions in data compared to MC predictions.

Unfolded E3C/E2C distributions in data compared to MC predictions.

The chi2 scan result using eq.3 for different alphas.

The correlation matrix is composed by 10 jet pt region, each region is represented by a block in the plot. Inside each block, there are 22 xL bins same as the E2C, E3C and E3C/E2C distributions. Therefore, the x and y bins of the correlation matrix is given by, binNumber = pT_index * 22 + xL_index.

The correlation matrix is composed by 10 jet pt region, each region is represented by a block in the plot. Inside each block, there are 22 xL bins same as the E2C, E3C and E3C/E2C distributions. Therefore, the x and y bins of the correlation matrix is given by, binNumber = pT_index * 22 + xL_index.

The correlation matrix is composed by 10 jet pt region, each region is represented by a block in the plot. Inside each block, there are 22 xL bins same as the E2C, E3C and E3C/E2C distributions. Therefore, the x and y bins of the correlation matrix is given by, binNumber = pT_index * 22 + xL_index.

The energy weight of E2C as defined in eq.1 in the paper is used in the unfolding, the binning is listed in this table.

The energy weight of E3C as defined in eq.1 in the paper is used in the unfolding, the binning is listed in this table.

The elliptic flow after nonflow subtraction, $v_{2}^{sub}/4$, for $f_0(980)$ as a function of $<KE_{T}>/4$ in pPb collision at 8.16 TeV.

The elliptic flow after nonflow subtraction, $v_{2}^{sub}/2$, for $f_0(980)$ as a function of $<p_{T}>/2$ in pPb collision at 8.16 TeV.

The elliptic flow after nonflow subtraction, $v_{2}^{sub}/4$, for $f_0(980)$ as a function of $<p_{T}>/4$ in pPb collision at 8.16 TeV.

Raw primary antihelium3-to-helium3 ratio from Geant4-based MC simulations as a function of the momentum p_primary with sigma_inel(3Hebar)x1.5.

Raw primary antihelium3-to-helium3 ratio as a function of the momentum p_primary.

Raw primary antihelium3-to-helium3 ratio from Geant4-based MC simulations as a function of the momentum p_primary with default sigma_inel(3Hebar).

Raw primary antihelium3-to-helium3 ratio from Geant4-based MC simulations as a function of the momentum p_primary with sigma_inel(3Hebar)x0.5.

Raw primary antihelium3-to-helium3 ratio from Geant4-based MC simulations as a function of the momentum p_primary with sigma_inel(3Hebar)x1.5.

Inelastic interaction cross section of antihelium-3 on an averaged material element of the ALICE detector (ITS+TPC analysis).

Inelastic interaction cross section of antihelium-3 on an averaged material element of the ALICE detector (ITS+TPC+TOF analysis).

Inelastic interaction cross section of antihelium-3 on an averaged material element of the ALICE detector (ITS+TPC+TOF analysis).

Antihelium-3 cosmic-ray flux from background before the solar modulation

Antihelium-3 cosmic-ray flux from dark matter before the solar modulation

Transparency of the galaxy to antihelium-3 nuclei from background before the solar modulation

Transparency of the galaxy to antihelium-3 nuclei from dark matter before the solar modulation

Antihelium-3 cosmic-ray flux from background after the solar modulation

Antihelium-3 cosmic-ray flux from dark matter after the solar modulation

Transparency of the galaxy to antihelium-3 nuclei from background after the solar modulation

Transparency of the galaxy to antihelium-3 nuclei from dark matter after the solar modulation

Inelastic interaction cross section of antihelium-3 on proton from the fit to the ALICE data.

Inelastic interaction cross section of antihelium-3 on helium-4 from the fit to the ALICE data.

Ratio of charged-particle jet yield measured in different multiplicity classes with respect to that in MB events as a function of pT for different resolution parameters R from 0.2 to 0.7. Statistical and total systematic uncertainties are shown as vertical error bars and boxes around the data points, respectively.

Ratios of charged-particle jet spectra with R = 0.2 to that with other jet resolution parameters R from 0.3 to 0.7, shown in different V0M multiplicity classes. Statistical and systematic uncertainties are shown as vertical error bars and boxes around the data points, respectively

a) Integrated jet yields, and b) $<p_T>$ of jets with 5 $<p^{ch}_{T,jet}<$ 100 GeV/$c$ as a function of selfnormalised charged-particle multiplicity for different resolution parameters R varied from 0.2 to 0.7, with the charged-particle multiplicities provided in Ref. [67]. Statistical and systematic uncertainties are shown as vertical error bars and boxes around the data points, respectively.

Self-normalised jet yields as a function of the self-normalised charged-particle multiplicity for different resolution parameters R varied from 0.2 to 0.7 in different jet pT intervals: a) 5 $<p^{ch}_{T,jet}<$ 7 GeV/$c$, b) 9 $<p^{ch}_{T,jet}<$ 12 GeV/$c$, c) 30 $<p^{ch}_{T,jet}<$ 50 GeV/$c$, and d) 70 $<p^{ch}_{T,jet}<$ 100 GeV/$c$. The charged-particle multiplicities are taken from Ref. [67]. Statistical and systematic uncertainties are shown as vertical error bars and boxes around the data points, respectively.

Inclusive charged-particle jet cross sections in pp collisions at $\sqrt{s}$ = 13 TeV using the anti-kT algorithm for different jet resolution parameters R from 0.2 to 0.7, without UE subtraction. Statistical uncertainties are displayed as vertical error bars. The total systematic uncertainties are shown as solid boxes around the data points.

Ratio of charged-particle jet cross section for resolution parameter R = 0.2 to other radii R = X, with X ranging from 0.3 to 0.7, without UE subtraction. Data are compared with LO (PYTHIA) and NLO (POWHEG+PYTHIA8) predictions as shown in the bottom panels. The systematic uncertainties of the cross section ratios from data are indicated by solid boxes around data points in the upper panel and shaded bands around unity in the mid and lower panels. No uncertainties are shown for theoretical predictions for better visibility.

The $\mathrm{\Lambda_c^+ / D^0}$ production yield ratio in central (0--10\%) and mid-central (30--50\%) Pb--Pb collisions at $\sqrt{s_{\mathrm{NN}}}=5.02$ TeV.

The nuclear modification factor of prompt $\mathrm{\Lambda_c^+}$ as a function transverse momentum in central (0--10\%) and mid-central (30--50\%) Pb--Pb collisions at $\sqrt{s_{\mathrm{NN}}}=5.02$ TeV.

The nuclear modification factor of prompt $\mathrm{\Lambda_c^+}$ as a function transverse momentum in central (0--10\%) and mid-central (30--50\%) Pb--Pb collisions at $\sqrt{s_{\mathrm{NN}}}=5.02$ TeV.

(Anti)proton spectrum in V0M multiplicity class IV - V

(Anti)proton spectrum in V0M multiplicity class VI

(Anti)proton spectrum in V0M multiplicity class VII

(Anti)proton spectrum in V0M multiplicity class VIII

(Anti)proton spectrum in V0M multiplicity class IX

(Anti)proton spectrum in V0M multiplicity class X

(Anti)proton spectrum in V0M multiplicity class INEL > 0

(Anti)deuteron spectrum in V0M multiplicity class I

(Anti)deuteron spectrum in V0M multiplicity class II

(Anti)deuteron spectrum in V0M multiplicity class III

(Anti)deuteron spectrum in V0M multiplicity class IV - V

(Anti)deuteron spectrum in V0M multiplicity class VI

(Anti)deuteron spectrum in V0M multiplicity class VII

(Anti)deuteron spectrum in V0M multiplicity class VIII

(Anti)deuteron spectrum in V0M multiplicity class IX

(Anti)deuteron spectrum in V0M multiplicity class X

(Anti)deuteron spectrum in V0M multiplicity class INEL > 0

(Anti)helion spectrum in V0M multiplicity class I - III

(Anti)helion spectrum in V0M multiplicity class IV - X

(Anti)helion spectrum in V0M multiplicity class INEL > 0

Mean transverse momentum $<p_\mathrm{T}>$ of protons as a function of the charged particle multiplicity density at mid-rapidity

Mean transverse momentum $<p_\mathrm{T}>$ of deuterons as a function of the charged particle multiplicity density at mid-rapidity

Mean transverse momentum $<p_\mathrm{T}>$ of helions as a function of the charged particle multiplicity density at mid-rapidity

Coalescence parameter $B_2$ as a function of $p_{\mathrm{T}}$ in V0M multiplicity class I

Coalescence parameter $B_2$ as a function of $p_{\mathrm{T}}$ in V0M multiplicity class II

Coalescence parameter $B_2$ as a function of $p_{\mathrm{T}}$ in V0M multiplicity class III

Coalescence parameter $B_2$ as a function of $p_{\mathrm{T}}$ in V0M multiplicity class IV - V

Coalescence parameter $B_2$ as a function of $p_{\mathrm{T}}$ in V0M multiplicity class VI

Coalescence parameter $B_2$ as a function of $p_{\mathrm{T}}$ in V0M multiplicity class VII

Coalescence parameter $B_2$ as a function of $p_{\mathrm{T}}$ in V0M multiplicity class VIII

Coalescence parameter $B_2$ as a function of $p_{\mathrm{T}}$ in V0M multiplicity class IX

Coalescence parameter $B_2$ as a function of $p_{\mathrm{T}}$ in V0M multiplicity class X

Coalescence parameter $B_2$ as a function of $p_{\mathrm{T}}$ in V0M multiplicity class INEL > 0

Coalescence parameter $B_3$ as a function of $p_{\mathrm{T}}$ in V0M multiplicity class I - III

Coalescence parameter $B_3$ as a function of $p_{\mathrm{T}}$ in V0M multiplicity class IV - X

Coalescence parameter $B_3$ as a function of $p_{\mathrm{T}}$ in V0M multiplicity class INEL > 0

Coalescence parameter $B_2$ as a function of the charged particle multiplicity density at mid-rapidity

Coalescence parameter $B_3$ as a function of the charged particle multiplicity density at mid-rapidity

Deuteron over proton ratio as a function of the charged particle multiplicity density at mid-rapidity

Helion over proton ratio as a function of the charged particle multiplicity density at mid-rapidity

Transverse-momentum spectra ratios of $\mathrm{D^0}$ hadrons measured in pp collisions at $\sqrt{{s}} = 13$~TeV for different multiplicity classes selected with the $N_\mathrm{trkl}$ estimator at midrapidityto the INEL>0 event class.

Transverse-momentum spectra ratios of $\mathrm{D_s^+}$ hadrons measured in pp collisions at $\sqrt{{s}} = 13$~TeV for different multiplicity classes selected with the $N_\mathrm{trkl}$ estimator at midrapidityto the INEL>0 event class.

Transverse-momentum spectra ratios of $\mathrm{\Lambda_c^+}$ hadrons measured in pp collisions at $\sqrt{{s}} = 13$~TeV for different multiplicity classes selected with the $N_\mathrm{trkl}$ estimator at midrapidityto the INEL>0 event class.

The $\mathrm{D_s^+ / D^0}$ ratio measured in pp collisions at $\sqrt{s} = 13$~TeV for different multiplicity classes selected with the $N_\mathrm{trkl}$ estimator at midrapidity

The $\mathrm{D_s^+ / D^0}$ ratio measured in pp collisions at $\sqrt{s} = 13$~TeV for different multiplicity classes selected with the $N_\mathrm{trkl}$ estimator at midrapidity

The $\mathrm{\Lambda_c^+ / D^0}$ ratio measured in pp collisions at $\sqrt{s} = 13$~TeV for different multiplicity classes selected with the V0M estimator at forward rapidity

The $\mathrm{\Lambda_c^+ / D^0}$ ratio measured in pp collisions at $\sqrt{s} = 13$~TeV for different multiplicity classes selected with the V0M estimator at forward rapidity

The extrapolated $p_\mathrm{T}$-integrated $\mathrm{\Lambda_c^+ / D^0}$ ratio measured as a function of $\langle \mathrm{d}N_\mathrm{ch}/\mathrm{d}\eta \rangle$ in pp collisions at $\sqrt{{s}} = 13$~TeV

Transverse-momentum spectra of $\mathrm{D^0}$ hadrons measured in pp collisions at $\sqrt{{s}} = 13$~TeV for different multiplicity classes selected with the $N_\mathrm{trkl}$ estimator at midrapidity.

Transverse-momentum spectra of $\mathrm{D_s^+}$ hadrons measured in pp collisions at $\sqrt{{s}} = 13$~TeV for different multiplicity classes selected with the $N_\mathrm{trkl}$ estimator at midrapidity.

Transverse-momentum spectra of $\mathrm{\Lambda_c^+}$ hadrons measured in pp collisions at $\sqrt{{s}} = 13$~TeV for different multiplicity classes selected with the $N_\mathrm{trkl}$ estimator at midrapidity.

Transverse-momentum spectra ratios of $\mathrm{D^0}$ hadrons measured in pp collisions at $\sqrt{{s}} = 13$~TeV for different multiplicity classes selected with the $N_\mathrm{trkl}$ estimator at midrapidityto the INEL>0 event class.

Transverse-momentum spectra ratios of $\mathrm{D_s^+}$ hadrons measured in pp collisions at $\sqrt{{s}} = 13$~TeV for different multiplicity classes selected with the $N_\mathrm{trkl}$ estimator at midrapidityto the INEL>0 event class.

Transverse-momentum spectra ratios of $\mathrm{\Lambda_c^+}$ hadrons measured in pp collisions at $\sqrt{{s}} = 13$~TeV for different multiplicity classes selected with the $N_\mathrm{trkl}$ estimator at midrapidityto the INEL>0 event class.

The visible $p_\mathrm{T}$-integrated $\mathrm{D_s^+ / D^0}$ ratio measured as a function of $\langle \mathrm{d}N_\mathrm{ch}/\mathrm{d}\eta \rangle$ in pp collisions at $\sqrt{{s}} = 13$~TeV

The visible $p_\mathrm{T}$-integrated $\mathrm{\Lambda_c^+ / D^0}$ ratio measured as a function of $\langle \mathrm{d}N_\mathrm{ch}/\mathrm{d}\eta \rangle$ in pp collisions at $\sqrt{{s}} = 13$~TeV

The extrapolated $p_\mathrm{T}$-integrated $\mathrm{D_s^+ / D^0}$ ratio measured as a function of $\langle \mathrm{d}N_\mathrm{ch}/\mathrm{d}\eta \rangle$ in pp collisions at $\sqrt{{s}} = 13$~TeV

Centrality dependence of $\rho\left(v_{2}^{2}, v_{3}^{2}, [p_{\rm T}] \right)$ in Pb--Pb collisions at $\sqrt{s_{\rm NN}}$ = 5.02 TeV

$p_{\mathrm T}$-differential yield of $\frac{\mathrm{K^{*0}} + \overline{\mathrm{K^{*0}}}}{2}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV (Multiplicity class 10-20%).

$p_{\mathrm T}$-differential yield of $\frac{\mathrm{K^{*0}} + \overline{\mathrm{K^{*0}}}}{2}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV (Multiplicity class 20-40%).

$p_{\mathrm T}$-differential yield of $\frac{\mathrm{K^{*0}} + \overline{\mathrm{K^{*0}}}}{2}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV (Multiplicity class 40-60%).

$p_{\mathrm T}$-differential yield of $\frac{\mathrm{K^{*0}} + \overline{\mathrm{K^{*0}}}}{2}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV (Multiplicity class 60-80%).

$p_{\mathrm T}$-differential yield of $\frac{\mathrm{K^{*0}} + \overline{\mathrm{K^{*0}}}}{2}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV (Multiplicity class 80-100%).

$p_{\mathrm T}$-differential yield of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV (NSD).

$p_{\mathrm T}$-differential yield of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV (Multiplicity class 0-5%).

$p_{\mathrm T}$-differential yield of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV (Multiplicity class 5-10%).

$p_{\mathrm T}$-differential yield of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV (Multiplicity class 10-20%).

$p_{\mathrm T}$-differential yield of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV (Multiplicity class 20-40%).

$p_{\mathrm T}$-differential yield of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV (Multiplicity class 40-60%).

$p_{\mathrm T}$-differential yield of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV (Multiplicity class 60-80%).

$p_{\mathrm T}$-differential yield of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV (Multiplicity class 80-100%).

$p_{\mathrm T}$-differential $R_{\mathrm{pPb}}$ of K$^{*0}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV.

$p_{\mathrm T}$-differential $R_{\mathrm{pPb}}$ of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV.

$p_{\mathrm T}$-differential $R_{\mathrm{pPb}}$ of K$^{*0}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV.

$p_{\mathrm T}$-differential $R_{\mathrm{pPb}}$ of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV.

$p_{\mathrm T}$-differential $R_{\mathrm{pPb}}$ of $\Xi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV.

$p_{\mathrm T}$-differential $R_{\mathrm{pPb}}$ of $\Omega$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV.

$p_\mathrm{T}$-differential production yield of prompt $\mathrm{D_s^+}$ mesons in the 30-50% centrality class in Pb-Pb collisions at $\sqrt{{s_\mathrm{NN}}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y|<0.5$. Branching ratio of $\mathrm{D_s^+\to \phi\pi^+\to K^-K^+\pi^+}$: 0.0224

$p_\mathrm{T}$-differential nuclear modification factor $R_\mathrm{AA}$ of prompt $\mathrm{D_s^+}$ mesons in the 0-10% centrality class in Pb-Pb collisions at $\sqrt{{s_\mathrm{NN}}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y|<0.5$.

$p_\mathrm{T}$-differential nuclear modification factor $R_\mathrm{AA}$ of prompt $\mathrm{D_s^+}$ mesons in the 0-10% centrality class in Pb-Pb collisions at $\sqrt{{s_\mathrm{NN}}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y|<0.5$.

$p_\mathrm{T}$-differential nuclear modification factor $R_\mathrm{AA}$ of prompt $\mathrm{D_s^+}$ mesons in the 30-50% centrality class in Pb-Pb collisions at $\sqrt{{s_\mathrm{NN}}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y|<0.5$.

$p_\mathrm{T}$-differential nuclear modification factor $R_\mathrm{AA}$ of prompt $\mathrm{D_s^+}$ mesons in the 30-50% centrality class in Pb-Pb collisions at $\sqrt{{s_\mathrm{NN}}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y|<0.5$.

$p_\mathrm{{T}}$-differential production ratios of prompt $\mathrm{D}_s^+/\mathrm{D}^0$ mesons in the 0-10% centrality class in Pb-Pb collisions at $\sqrt{{s_\mathrm{NN}}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y|<0.5$. Branching ratio of $\mathrm{D_s^+\to \phi\pi^+\to K^-K^+\pi^+}$: 0.0224 Branching ratio of $\mathrm{D^0 \to K^-\pi^+}$: 0.0395

$p_\mathrm{{T}}$-differential production ratios of prompt $\mathrm{D}_s^+/\mathrm{D}^0$ mesons in the 0-10% centrality class in Pb-Pb collisions at $\sqrt{{s_\mathrm{NN}}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y|<0.5$. Branching ratio of $\mathrm{D_s^+\to \phi\pi^+\to K^-K^+\pi^+}$: 0.0224 Branching ratio of $\mathrm{D^0 \to K^-\pi^+}$: 0.0395

$p_\mathrm{{T}}$-differential production ratios of prompt $\mathrm{D}_s^+/\mathrm{D}^0$ mesons in the 30-50% centrality class in Pb-Pb collisions at $\sqrt{{s_\mathrm{NN}}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y|<0.5$. Branching ratio of $\mathrm{D_s^+\to \phi\pi^+\to K^-K^+\pi^+}$: 0.0224 Branching ratio of $\mathrm{D^0 \to K^-\pi^+}$: 0.0395

$p_\mathrm{{T}}$-differential production ratios of prompt $\mathrm{D}_s^+/\mathrm{D}^0$ mesons in the 30-50% centrality class in Pb-Pb collisions at $\sqrt{{s_\mathrm{NN}}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $|y|<0.5$. Branching ratio of $\mathrm{D_s^+\to \phi\pi^+\to K^-K^+\pi^+}$: 0.0224 Branching ratio of $\mathrm{D^0 \to K^-\pi^+}$: 0.0395

$p_\mathrm{{T}}$-differential prompt $\mathrm{D}_s^+/\mathrm{D}^0$ ratios in Pb-Pb collisions in the 0-10% centrality class at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$ divided by those in pp collisions at the same energy.

$p_\mathrm{{T}}$-differential prompt $\mathrm{D}_s^+/\mathrm{D}^0$ ratios in Pb-Pb collisions in the 0-10% centrality class at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$ divided by those in pp collisions at the same energy.

$p_\mathrm{{T}}$-differential prompt $\mathrm{D}_s^+/\mathrm{D}^0$ ratios in Pb-Pb collisions in the 30-50% centrality class at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$ divided by those in pp collisions at the same energy.

$p_\mathrm{{T}}$-differential prompt $\mathrm{D}_s^+/\mathrm{D}^0$ ratios in Pb-Pb collisions in the 30-50% centrality class at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$ divided by those in pp collisions at the same energy.

$p_\mathrm{T}$-differential elliptic flow coefficient $v_2$ of prompt $\mathrm{D_s^+}$ mesons in the 30-50% centrality class in Pb-Pb collisions at $\sqrt{{s_\mathrm{NN}}}=5.02~\mathrm{{TeV}}$.

$p_\mathrm{T}$-differential elliptic flow coefficient $v_2$ of prompt $\mathrm{D_s^+}$ mesons in the 30-50% centrality class in Pb-Pb collisions at $\sqrt{{s_\mathrm{NN}}}=5.02~\mathrm{{TeV}}$.

Comparison of the azimuthal-correlation distributions of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) with $3 < p_{\rm T} < 5$ GeV/$c$ and charged particles with $1 < p_{\rm T} < 2$ GeV/$c$, in pp collisions at $\sqrt{s} =$ 5 TeV and 13 TeV, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) with $8 < p_{\rm T} < 16$ GeV/$c$ and charged particles with $1 < p_{\rm T} < 2$ GeV/$c$, in pp collisions at $\sqrt{s} =$ 5 TeV and 13 TeV, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) with $16 < p_{\rm T} < 24$ GeV/$c$ and charged particles with $1 < p_{\rm T} < 2$ GeV/$c$, in pp collisions at $\sqrt{s} =$ 5 TeV and 13 TeV, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the near-side associated peak yields measured in pp collisions at $\sqrt{s} = $5, 7, and 13 TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $0.3 < p_{\rm T} < 1$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the near-side associated peak yields measured in pp collisions at $\sqrt{s}=$ 5 TeV and 13 TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $1 < p_{\rm T} < 2$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the near-side associated peak width measured in pp collisions at $\sqrt{s} = $5, 7, and 13 TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $0.3 < p_{\rm T} < 1$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the near-side associated peak width measured in pp collisions at $\sqrt{s}=$ 5 TeV and 13 TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $1 < p_{\rm T} < 2$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the away-side associated peak yields measured in pp collisions at $\sqrt{s}=$ 5 TeV and 13 TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $0.3 < p_{\rm T} < 1$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the away-side associated peak yields measured in pp collisions at $\sqrt{s}=$ 5 TeV and 13 TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $1 < p_{\rm T} < 2$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the away-side associated peak width measured in pp collisions at $\sqrt{s}=$ 5 TeV and 13 TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $0.3 < p_{\rm T} < 1$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the away-side associated peak width measured in pp collisions at $\sqrt{s}=$ 5 TeV and 13 TeV, for correlations of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $1 < p_{\rm T} < 2$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the azimuthal-correlation distributions of of D$^{0}$ mesons with $3 < p_{\rm T} < 5$ GeV/$c$ and charged particles with $p_{\rm T} > 0.3$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D$^{0}$ mesons with $5 < p_{\rm T} < 8$ GeV/$c$ and charged particles with $p_{\rm T} > 0.3$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D$^{0}$ mesons with $8 < p_{\rm T} < 16$ GeV/$c$ and charged particles with $p_{\rm T} > 0.3$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D$^{0}$ mesons with $16 < p_{\rm T} < 24$ GeV/$c$ and charged particles with $p_{\rm T} > 0.3$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D$^{0}$ mesons with $3 < p_{\rm T} < 5$ GeV/$c$ and charged particles with $0.3 < p_{\rm T} < 1$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D$^{0}$ mesons with $5 < p_{\rm T} < 8$ GeV/$c$ and charged particles with $0.3 < p_{\rm T} < 1$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D$^{0}$ mesons with $8 < p_{\rm T} < 16$ GeV/$c$ and charged particles with $0.3 < p_{\rm T} < 1$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D$^{0}$ mesons with $16 < p_{\rm T} < 24$ GeV/$c$ and charged particles with $0.3 < p_{\rm T} < 1$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D$^{0}$ mesons with $3 < p_{\rm T} < 5$ GeV/$c$ and charged particles with $1 < p_{\rm T} < 2$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D$^{0}$ mesons with $5 < p_{\rm T} < 8$ GeV/$c$ and charged particles with $1 < p_{\rm T} < 2$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D$^{0}$ mesons with $8 < p_{\rm T} < 16$ GeV/$c$ and charged particles with $1 < p_{\rm T} < 2$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D$^{0}$ mesons with $16 < p_{\rm T} < 24$ GeV/$c$ and charged particles with $1 < p_{\rm T} < 2$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D$^{0}$ mesons with $3 < p_{\rm T} < 5$ GeV/$c$ and charged particles with $2 < p_{\rm T} < 3$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D$^{0}$ mesons with $5 < p_{\rm T} < 8$ GeV/$c$ and charged particles with $2 < p_{\rm T} < 3$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D$^{0}$ mesons with $8 < p_{\rm T} < 16$ GeV/$c$ and charged particles with $2 < p_{\rm T} < 3$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the azimuthal-correlation distributions of D$^{0}$ mesons with $16 < p_{\rm T} < 24$ GeV/$c$ and charged particles with $2 < p_{\rm T} < 3$ GeV/$c$, in pp collisions at $\sqrt{s} = $13 TeV in four different V0 multiplicity classes, after baseline subtraction. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$. The azimuthal-correlation distributions are reported in the range $0 < \Delta\varphi < \pi$.

Comparison of the near-side associated peak yields among four different V0 multiplicity classes measured in pp collisions at $\sqrt{s} = $13 TeV, for correlations of D$^{0}$ mesons and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $p_{\rm T} > 0.3 $ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the near-side associated peak yields among four different V0 multiplicity classes measured in pp collisions at $\sqrt{s} = $13 TeV, for correlations of D$^{0}$ mesons and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $0.3 < p_{\rm T} < 1 $ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the near-side associated peak yields among four different V0 multiplicity classes measured in pp collisions at $\sqrt{s} = $13 TeV, for correlations of D$^{0}$ mesons and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $1 < p_{\rm T} < 2 $ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the near-side associated peak yields among four different V0 multiplicity classes measured in pp collisions at $\sqrt{s} = $13 TeV, for correlations of D$^{0}$ mesons and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $2 < p_{\rm T} < 3 $ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the near-side associated peak width among four different V0 multiplicity classes measured in pp collisions at $\sqrt{s} = $13 TeV, for correlations of D$^{0}$ mesons and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $ p_{\rm T} > 0.3 $ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the near-side associated peak width among four different V0 multiplicity classes measured in pp collisions at $\sqrt{s} = $13 TeV, for correlations of D$^{0}$ mesons and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $0.3 < p_{\rm T} < 1 $ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the near-side associated peak width among four different V0 multiplicity classes measured in pp collisions at $\sqrt{s} = $13 TeV, for correlations of D$^{0}$ mesons and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $1 < p_{\rm T} < 2 $ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the near-side associated peak width among four different V0 multiplicity classes measured in pp collisions at $\sqrt{s} = $13 TeV, for correlations of D$^{0}$ mesons and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $2 < p_{\rm T} < 3 $ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the away-side associated peak yields among four different V0 multiplicity classes measured in pp collisions at $\sqrt{s} = $13 TeV, for correlations of D$^{0}$ mesons and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $ p_{\rm T} > 0.3 $ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Comparison of the away-side associated peak width among four different V0 multiplicity classes measured in pp collisions at $\sqrt{s} = $13 TeV, for correlations of D$^{0}$ mesons and charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $ p_{\rm T} > 0.3 $ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Near-side associated peak yield of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $ p_{\rm T} > 0.3 $ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Near-side associated peak yield of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $0.3 < p_{\rm T} < 1$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Near-side associated peak yield of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $1 < p_{\rm T} < 2$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Near-side associated peak yield of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $2 < p_{\rm T} < 3$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Near-side associated peak width of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $ p_{\rm T} > 0.3 $ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Near-side associated peak width of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $0.3 < p_{\rm T} < 1$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Near-side associated peak width of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $1 < p_{\rm T} < 2$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Near-side associated peak width of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $2 < p_{\rm T} < 3$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Away-side associated peak yield of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $ p_{\rm T} > 0.3 $ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Away-side associated peak yield of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $0.3 < p_{\rm T} < 1$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Away-side associated peak yield of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $1 < p_{\rm T} < 2$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Away-side associated peak yield of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $2 < p_{\rm T} < 3$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Away-side associated peak width of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $ p_{\rm T} > 0.3 $ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Away-side associated peak width of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $0.3 < p_{\rm T} < 1$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Away-side associated peak width of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $1 < p_{\rm T} < 2$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Away-side associated peak width of D mesons (average of D$^{0}$, D$^{+}$, and D$^{*+}$) correlated with charged particles, as a function of the D-meson $p_{\rm T}$, for associated particle $2 < p_{\rm T} < 3$ GeV/$c$. Rapidity range for the D mesons is $|y^{\rm D}_{\rm cms}| < 0.5$. Correlations are integrated for $|\Delta\eta|=|\eta_{\rm ch}-\eta_{\rm D}| < 1$.

Transverse momentum distributions d$N$/(d$y$d$p_{\mathrm{T}}$) of prompt D$^{+}$ (middle) meson in the 30-50% (diamonds) centrality classes in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV in |y|<0.5. The uncertainties on the BRs are quoted separately and the horizontal bars represent bin widths. The points in the figure of the paper are scaled for visibility by a factor 10$^{-1}$. Branching ratio of D$^+$ --> K$\pi\pi$ = 0.0938

Transverse momentum distributions d$N$/(d$y$d$p_{\mathrm{T}}$) of prompt D*$^{+}$ (right) meson in the 0-10% (crosses) centrality classes in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV in |y|<0.5. The uncertainties on the BRs are quoted separately and the horizontal bars represent bin widths. Branching ratio of D*$^+$ --> D$^0\pi$ = 0.677

Transverse momentum distributions d$N$/(d$y$d$p_{\mathrm{T}}$) of prompt D*$^{+}$ (right) meson in the 30-50% (diamonds) centrality classes in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV in |y|<0.5. The uncertainties on the BRs are quoted separately and the horizontal bars represent bin widths. The points in the figure of the paper are scaled for visibility by a factor 10$^{-1}$. Branching ratio of D*$^+$ --> D$^0\pi$ = 0.677

Ratio D$^{+}$/D$^{0}$ in central Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV, in |y|<0.5 and as a function of $p_\mathrm{T}$. Statistical (bars) and systematic (boxes) uncertainties are shown. Branching ratio of D$^0$ --> K$\pi$ = 0.0395 Branching ratio of D$^+$ --> K$\pi\pi$ = 0.0938

Ratio D$^{+}$/D$^{0}$ in semicentral Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV, in |y|<0.5 and as a function of $p_\mathrm{T}$. Statistical (bars) and systematic (boxes) uncertainties are shown. Branching ratio of D$^0$ --> K$\pi$ = 0.0395 Branching ratio of D$^+$ --> K$\pi\pi$ = 0.0938

Ratio D*$^{+}$/D$^{0}$ in central Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV, in |y|<0.5 and as a function of $p_\mathrm{T}$. Statistical (bars) and systematic (boxes) uncertainties are shown. Branching ratio of D$^0$ --> K$\pi$ = 0.0395 Branching ratio of D*$^+$ --> D$^0\pi$ = 0.677

Ratio D*$^{+}$/D$^{0}$ in semicentral Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV, in |y|<0.5 and as a function of $p_\mathrm{T}$. Statistical (bars) and systematic (boxes) uncertainties are shown. Branching ratio of D$^0$ --> K$\pi$ = 0.0395 Branching ratio of D*$^+$ --> D$^0\pi$ = 0.677

Nuclear modification factor of D$^{0}$ meson in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV, in the 0-10% centrality class, in |y|<0.5 and as a function of $p_\mathrm{T}$.

Nuclear modification factor of D$^{+}$ meson in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV, in the 0-10% centrality class, in |y|<0.5 and as a function of $p_\mathrm{T}$.

Nuclear modification factor of D*$^{+}$ meson in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV, in the 0-10% centrality class, in |y|<0.5 and as a function of $p_\mathrm{T}$.

Nuclear modification factor of D$^{0}$ meson in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV, in the 30-50% centrality class, in |y|<0.5 and as a function of $p_\mathrm{T}$.

Nuclear modification factor of D$^{+}$ meson in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV, in the 30-50% centrality class, in |y|<0.5 and as a function of $p_\mathrm{T}$.

Nuclear modification factor of D*$^{+}$ meson in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV, in the 30-50% centrality class, in |y|<0.5 and as a function of $p_\mathrm{T}$.

Average nuclear modification factor of D$^{0}$, D$^{+}$ and D*$^{+}$ meson in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV, in the 0-10% centrality class, in |y|<0.5 and as a function of $p_\mathrm{T}$.

Average nuclear modification factor of D$^{0}$, D$^{+}$ and D*$^{+}$ meson in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV, in the 30-50% centrality class, in |y|<0.5 and as a function of $p_\mathrm{T}$.

$p_\mathrm{T}$-integrated nuclear modification factor of D$^{0}$ meson in central (0-10%) and semicentral (30-50%) Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02$ TeV, in |y|<0.5.