Post-fit distribution of the GNN score evaluated with $m_{H/A}$ = 400 GeV in the validation region in the 2LOS region with $\geq 8$ jets. These regions do not enter the fit. The post-fit background prediction is obtained using the post-fit nuisance parameters from the background-only fit in the control and signal regions.

Observed and expected 95% CL upper limits on the cross-section times branching ratio,$ \sigma(pp \rightarrow t\bar{t}H/A ) \times B(H/A \rightarrow t\bar{t})$, as a function of $m_{H/A}$, obtained using the 1L/2LOS final states.

Observed and expected 95% CL lower limits on tanβ as a function of the $m_{H/A}$ mass obtained using the 1L/2LOS final states, assuming a Type-II 2HDM in the alignment limit. Values of tanβ below the observed limit are excluded. The scenario where both the scalar $H$ and pseudo-scalar $A$ contribute with $m_H = m_A$ is shown.

Observed and expected 95% CL lower limits on tanβ as a function of the $m_{H/A}$ mass obtained using the 1L/2LOS final states, assuming a Type-II 2HDM in the alignment limit. Values of tanβ below the observed limit are excluded. The scenario with the contribution from scalar $H$ only is shown.

Observed and expected 95% CL lower limits on tanβ as a function of the $m_{H/A}$ mass obtained using the 1L/2LOS final states, assuming a Type-II 2HDM in the alignment limit. Values of tanβ below the observed limit are excluded. The scenario with the contribution from pseudo-scalar $A$ only is shown.

Expected and observed 95% CL upper limits on the cross-section times branching ratio,$ \sigma(pp \rightarrow t\bar{t}H/A ) \times B(H/A \rightarrow t\bar{t})$, as a function of $m_{H/A}$, obtained from the combination of the 1L/2LOS and 2LSS/ML final states.

Expected and observed 95% CL lower limits on tanβ as a function of the $m_{H/A}$ mass obtained using the 1L/2LOS and 2LSS/ML final states, assuming a Type-II 2HDM in the alignment limit. Values of tanβ below the observed limit are excluded. The scenario where both the scalar $H$ and pseudo-scalar $A$ contribute with $m_H = m_A$ is shown.

Expected and observed 95% CL lower limits on tanβ as a function of the $m_{H/A}$ mass obtained using the 1L/2LOS and 2LSS/ML final states, assuming a Type-II 2HDM in the alignment limit. Values of tanβ below the observed limit are excluded. The scenario with the contribution from scalar $H$ only is shown.

Expected and observed 95% CL lower limits on tanβ as a function of the $m_{H/A}$ mass obtained using the 1L/2LOS and 2LSS/ML final states, assuming a Type-II 2HDM in the alignment limit. Values of tanβ below the observed limit are excluded. The scenario with the contribution from pseudo-scalar $A$ only is shown.

Expected and Observed 95% CL upper limits on the $ \sigma(pp \rightarrow t\bar{t})\times B(S_8\rightarrow t\bar{t})$ production cross-section as a function of $m_{S_8}$, obtained from the combination of the 1L/2LOS and 2LSS/ML final states. The expected limits from the individual 1L/2LOS and 2LSS/ML analyses are also shown.

Post-fit distribution of the GNN scores evaluated with $m_{H/A}$ = 700 GeV in the 1L region with $\geq 10$ jets and four $b$-tagged jets. The fit is performed under the background-only hypothesis.

Post-fit distribution of the GNN scores evaluated with $m_{H/A}$ = 700 GeV in the 2LOS region with $\geq 8$ jets and $\geq 4$ $b$-tagged jets. The fit is performed under the background-only hypothesis.

Post-fit distribution of the GNN scores evaluated with $m_{H/A}=$1000 GeV,in the 1L region with $\geq 10$ jets four $b$-tagged jets. The fit is performed under the background-only hypothesis.

Post-fit distribution of the GNN scores evaluated with $m_{H/A}=$1000 GeV in the 2LOS region with $\geq 8$ jets and $\geq 4$ $b$-tagged jets. The fit is performed under the background-only hypothesis.

Post-fit distribution of the most important global features used in the GNN training, Sum of the pcb scores of the six jets with the highest scores ($\sum_{i\in[1,6]}\mathrm{pcb}_i$), in the region where the training is performed, i.e., the region with $\geq 9$ jets and $\geq3$ $b$-tagged jets in the 1L channel. The fit is performed under the background-only hypothesis using the GNN scores evaluated with $m_{H/A} = 400$ GeV.

Post-fit distribution of the most important global features used in the GNN training, Sum of the pcb scores of the six jets with the highest scores ($\sum_{i\in[1,6]}\mathrm{pcb}_i$), in the region where the training is performed i.e.the region with $\geq 7$ jets and $\geq3$ $b$-tagged jets in the 2LOS channel. The fit is performed under the background-only hypothesis using the GNN scores evaluated with $m_{H/A} = 400$ GeV.

Post-fit distribution of the most important global features used in the GNN training, $p_{\mathrm{T}}$ sum of all reconstructed leptons and jets ($\mathrm{H}_{\mathrm{T}}^{\mathrm{all}}$), in the region where the training is performed, i.e., the region with $\geq 9$ jets and $\geq3$ $b$-tagged jets in the 1L channel. The fit is performed under the background-only hypothesis using the GNN scores evaluated with $m_{H/A} = 400$ GeV.

Post-fit distribution of the most important global features used in the GNN training, $p_{\mathrm{T}}$ sum of all reconstructed leptons and jets ($\mathrm{H}_{\mathrm{T}}^{\mathrm{all}}$), in the region where the training is performed, i.e. the region with $\geq 7$ jets and $\geq3$ $b$-tagged jets in the 2LOS channel. The fit is performed under the background-only hypothesis using the GNN scores evaluated with $m_{H/A} = 400$ GeV.

Post-fit distribution of the most important global features used in the GNN training, Number of jets, in the region where the training is performed, i.e., the region with $\geq 9$ jets and $\geq3$ $b$-tagged jets in the 1L channel. The fit is performed under the background-only hypothesis using the GNN scores evaluated with $m_{H/A} = 400$ GeV.

Post-fit distribution of the most important global features used in the GNN training, Number of jets, in the region where the training is performed, i.e., the region with $\geq 7$ jets and $\geq3$ $b$-tagged jets in the 2LOS channel. The fit is performed under the background-only hypothesis using the GNN scores evaluated with $m_{H/A} = 400$ GeV.

Transverse momentum spectrum of charged hadrons measured at midrapidity ($|y|<0.5$) in INEL>0 pp collisions at $\sqrt{s}$ = 13 TeV for different flattenicity event classes selected with the V0M estimator at forward rapidity (top figure, upper panel)

Transverse momentum spectrum of $\pi^{+} + \pi^{-}$ measured at midrapidity ($|y|<0.5$) in INEL>0 pp collisions at $\sqrt{s}$ = 13 TeV for different flattenicity event classes and HM events (0--1% V0M) in the same flattenicity event classes (bottom figure, upper panel)

Transverse momentum spectrum of $K^{+} + K^{-}$ measured at midrapidity ($|y|<0.5$) in INEL>0 pp collisions at $\sqrt{s}$ = 13 TeV for different flattenicity event classes and HM events (0--1% V0M) in the same flattenicity event classes (bottom figure, upper panel)

Transverse momentum spectrum of $p + \overline{p}$ measured at midrapidity ($|y|<0.5$) in INEL>0 pp collisions at $\sqrt{s}$ = 13 TeV for different flattenicity event classes and HM events (0--1% V0M) in the same flattenicity event classes (bottom figure, upper panel)

Transverse momentum spectrum of charged hadrons measured at midrapidity ($|y|<0.5$) in INEL>0 pp collisions at $\sqrt{s}$ = 13 TeV for different flattenicity event classes and HM events (0--1% V0M) in the same flattenicity event classes (bottom figure, upper panel)

The $p_{\rm T}$-differential $(K^{+} + K^{-})/(\pi^{+}+\pi^{-})$ particle ratio for two extremes of flattenicity event classes measured in multiplicity-integrated events in pp collisions at $\sqrt{s}$ = 13 TeV

The $p_{\rm T}$-differential $(p + \overline{p})$/($\pi^{+}+\pi^{-}$) particle ratio for two extremes of flattenicity event classes measured in multiplicity-integrated events in pp collisions at $\sqrt{s}$ = 13 TeV

The $p_{\rm T}$-differential $(K^{+} + K^{-})/(\pi^{+}+\pi^{-})$ particle ratio for two extremes of flattenicity event classes measured in the 0-1% V0M event class in pp collisions at $\sqrt{s}$ = 13 TeV

The $p_{\rm T}$-differential $(p + \overline{p})$/($\pi^{+}+\pi^{-}$) particle ratio for two extremes of flattenicity event classes measured in the 0-1% V0M event class in pp collisions at $\sqrt{s}$ = 13 TeV

Transverse momentum-integrated $(K^{+} + K^{-})/(\pi^{+}+\pi^{-})$ particle ratio as a function of the average charged-particle density calculated for different flattenicity event classes in pp collisions at $\sqrt{s}$ = 13 TeV

Transverse momentum-integrated $(p + \overline{p})$/($\pi^{+}+\pi^{-}$) particle ratio as a function of the average charged-particle density calculated for different flattenicity event classes in pp collisions at $\sqrt{s}$ = 13 TeV

Average transverse momentum of $\pi^{+}+\pi^{-}$ as a function of the average charged-particle density calculated for different flattenicity event classes in pp collisions at $\sqrt{s}$ = 13 TeV

Average transverse momentum of $K^{+} + K^{-}$ as a function of the average charged-particle density calculated for different flattenicity event classes in pp collisions at $\sqrt{s}$ = 13 TeV

Average transverse momentum of $p + \overline{p}$ as a function of the average charged-particle density calculated for different flattenicity event classes in pp collisions at $\sqrt{s}$ = 13 TeV

Distribution of the reconstructed $\tau\tau$ invariant mass ($m_{\tau\tau}$) for all events in the VBF_1 signal region for $p_{\text{T}}^{H}>200$ GeV. The observed Higgs boson signal corresponds to $(\sigma\times B)/(\sigma\times B)_{\text{SM}}\,=\,0.99$. Entries with values above the $x$-axis range are shown in the last bin of each distribution. The prediction for each sample is determined from the likelihood fit performed to measure the total $pp\rightarrow H\rightarrow\tau\tau$ cross-section.

Distribution of the reconstructed $\tau\tau$ invariant mass ($m_{\tau\tau}$) for all events in the VBF_0 signal region for $m_{jj}<1$ TeV. The observed Higgs boson signal corresponds to $(\sigma\times B)/(\sigma\times B)_{\text{SM}}\,=\,0.99$. Entries with values above the $x$-axis range are shown in the last bin of each distribution. The prediction for each sample is determined from the likelihood fit performed to measure the total $pp\rightarrow H\rightarrow \tau\tau$ cross-section.

Distribution of the reconstructed $\tau\tau$ invariant mass ($m_{\tau\tau}$) for all events in the VBF_0 signal region for $m_{jj}>1$ TeV. The observed Higgs boson signal corresponds to $(\sigma\times B)/(\sigma\times B)_{\text{SM}}\,=\,0.99$. Entries with values above the $x$-axis range are shown in the last bin of each distribution. The prediction for each sample is determined from the likelihood fit performed to measure the total $pp\rightarrow H \rightarrow \tau\tau$ cross-section.

Distribution of the reconstructed $\tau\tau$ invariant mass ($m_{\tau\tau}$) for all events in the VBF_1 signal region for $m_{jj}<1$ TeV. The observed Higgs boson signal corresponds to $(\sigma\times B)/(\sigma\times B)_{\text{SM}}\,=\,0.99$. Entries with values above the $x$-axis range are shown in the last bin of each distribution. The prediction for each sample is determined from the likelihood fit performed to measure the total $pp\rightarrow \rightarrow \tau\tau$ cross-section.

Distribution of the reconstructed $\tau\tau$ invariant mass ($m_{\tau\tau}$) for all events in the VBF_1 signal region for $m_{jj}>1$ TeV. The observed Higgs boson signal corresponds to $(\sigma\times B)/(\sigma\times B)_{\text{SM}}\,=\,0.99$. Entries with values above the $x$-axis range are shown in the last bin of each distribution. The prediction for each sample is determined from the likelihood fit performed to measure the total $pp\rightarrow H\rightarrow \tau\tau$ cross-section.

Comparison of the observed and expected event yields in the ttH categories. The window category contains events with $m_{\tau\tau} \in [100, 150]$ GeV, while events outside of this mass region are included in the sideband region. The observed Higgs boson signal corresponds to $(\sigma\times B)/(\sigma\times B)_{\text{SM}}\,=\,0.99$. The prediction for each sample is determined from the likelihood fit performed to measure the total $pp\rightarrow H\rightarrow\tau\tau$ cross-section.

The measured values of $\sigma_{H}\times B(H\rightarrow\tau\tau)$ relative to the SM expectations in the total (labelled as 'Combined') and per-production-mode fit. The total $\pm1\sigma$ uncertainty in each measurement is indicated by the error bar. The total systematic uncertainty in each measurement is also indicated.

The measured values for $\sigma_{H}\times B(H\rightarrow\tau\tau)$ relative to the SM expectations in the simplified template cross-section measurement. The total $\pm1\sigma$ uncertainty in each measurement is indicated by the error bar. The total systematic uncertainty in each measurement is also indicated.

Upper limits at 95% CL on simplified template cross-sections in the individual ttH $p_{\text{T}}^{H}$ bins relative to their SM expectation in the $H\rightarrow\tau\tau$ decay channel, derived using the CL$_\text{s}$ method. The observed upper limits, together with the expected limits both in the background-only hypothesis and in the SM hypothesis (dotted red lines) are included. In the case of the expected limits in the background-only hypothesis, one- and two-standard-deviation uncertainty bands are also shown.

Measured fiducial differential cross-sections as a function of $p_\text{T}(j_0)$. The measurements are compared with particle-level SM predictions from the POWHEG+PYTHIA8, POWHEG+HERWIG7 and MadGraph5_aMC@NLO+Pythia8 generators for the combined VBF, ggF, $VH$, and $t\bar{t}H$ Higgs boson production modes. The shaded box around each data point shows the statistical uncertainty, while the total uncertainty is indicated by the error bar. The contribution from the ggF, $VH$, and $t\bar{t}H$ production modes as predicted by POWHEG+PYTHIA8 is also shown.

Measured fiducial differential cross-sections as a function of $p_\text{T}^{H}$. The measurements are compared with particle-level SM predictions from the POWHEG+PYTHIA8, POWHEG+HERWIG7 and MadGraph5_aMC@NLO+Pythia8 generators for the combined VBF, ggF, $VH$, and $t\bar{t} H$ Higgs boson production modes. The shaded box around each data point shows the statistical uncertainty, while the total uncertainty is indicated by the error bar. The contribution from the ggF, $VH$, and $t\bar{t} H$ production modes as predicted by POWHEG+PYTHIA8 is also shown.

Measured fiducial differential cross-sections as a function of $\Delta\phi_{jj}^{\text{signed}}$. The measurements are compared with particle-level SM predictions from the POWHEG+PYTHIA8, POWHEG+HERWIG7 and MadGraph5_aMC@NLO+Pythia8 generators for the combined VBF, ggF, $VH$, and $t\bar{t}H$ Higgs boson production modes. The shaded box around each data point shows the statistical uncertainty, while the total uncertainty is indicated by the error bar. The contribution from the ggF, $VH$, and $t\bar{t} H$ production modes as predicted by POWHEG+PYTHIA8 is also shown.

Measured fiducial differential cross-sections as a function of $\Delta\phi_{jj}^{\text{signed}}$ vs $p_\text{T}^\text{H}$. The measurements are compared with particle-level SM predictions from the POWHEG+PYTHIA8, POWHEG+HERWIG7 and MadGraph5_aMC@NLO+Pythia8 generators for the combined VBF, ggF, $VH$, and $t\bar{t}H$ Higgs boson production modes. The shaded box around each data point shows the statistical uncertainty, while the total uncertainty is indicated by the error bar. The contribution from the ggF, $VH$, and $t\bar{t}H$ production modes as predicted by POWHEG+PYTHIA8 is also shown.

The unfolded data compared with the nominal SM predictions from POWHEG+PYTHIA8, as well as the predictions for several values of Wilson coefficients in the SMEFT, for $\Delta\phi_{jj}^{\text{signed}}$.

The unfolded data compared with the nominal SM predictions from POWHEG+PYTHIA8, as well as the predictions for several values of Wilson coefficients in the SMEFT, for $\Delta\phi_{jj}^{\text{signed}}$.

The unfolded data compared with the nominal SM predictions from POWHEG+PYTHIA8, as well as the predictions for several values of Wilson coefficients in the SMEFT, for $\Delta\phi_{jj}^{\text{signed}}$ vs $p_\text{T}^{H}$.

The unfolded data compared with the nominal SM predictions from POWHEG+PYTHIA8, as well as the predictions for several values of Wilson coefficients in the SMEFT, for $p_\text{T}^{H}$.

The 95% expected and observed confidence intervals for each of the six considered Wilson coefficients. Each coefficient is treated individually and both the linear and linear + quadratic models are considered when evaluating the sensitivity to the terms that enter at $\mathcal{O}(\Lambda^{-4})$. For the CP-even operators the $\Delta\phi_{jj}^{\text{signed}}$ distribution is used to extract the confidence interval. The limits are computed at a new-physics scale $\Lambda=1$ TeV.

The 95% expected and observed confidence intervals for each of the six considered Wilson coefficients. Each coefficient is treated individually and both the linear and linear + quadratic models are considered when evaluating the sensitivity to the terms that enter at $\mathcal{O}(\Lambda^{-4})$. For the CP-odd operators the $\Delta\phi_{jj}^{\text{signed}}$ vs $p_\text{T}^{H}$ distribution is used to extract the confidence interval. The limits are computed at a new-physics scale $\Lambda=1$ TeV.

Two-dimensional observed and expected limits at 68% and 95% confidence level obtained from $c_{HW}$ and $c_{HB}$ using only the SM-dimension-6 interference in the SMEFT basis. The limits are computed at a new-physics scale $\Lambda$ = 1 TeV. In this case, the $\Delta\phi_{jj}^{signed}$ distribution is used. IMPORTANT: HEPData requires the same number of rows for each column. However, each curve of on the plot has a different number of points. For that reason, the curves are padded with data of "0.0, 0.0". The "0.0, 0.0" points should be omitted in analysis.

Two-dimensional observed and expected limits at 68% and 95% confidence level obtained from $c_{H\tilde{W}}$ and $c_{H\tilde{B}}$ using only the SM-dimension-6 interference in the SMEFT basis. The limits are computed at a new-physics scale $\Lambda$ = 1 TeV. In this purely CP-odd case, the $\Delta\phi_{jj}^{signed}$ vs $p_{T}^{H}$ distribution is used. IMPORTANT: HEPData requires the same number of rows for each column. However, each curve of on the plot has a different number of points. For that reason, the curves are padded with data of "0.0, 0.0". The "0.0, 0.0" points should be omitted in analysis.

Two-dimensional observed and expected limits at 68% and 95% confidence level obtained from $c_{HW}$ and $c_{H\tilde{W}}$ using only the SM-dimension-6 interference in the SMEFT basis. The limits are computed at a new-physics scale $\Lambda$ = 1 TeV. In this case, the $\Delta\phi_{jj}^{signed}$ distribution is used. IMPORTANT: HEPData requires the same number of rows for each column. However, each curve of on the plot has a different number of points. For that reason, the curves are padded with data of "0.0, 0.0". The "0.0, 0.0" points should be omitted in analysis.

The expected values of $\sigma_{H}\times B(H\rightarrow\tau\tau)$ relative to the SM expectations in the total (labelled as 'Combined') and per-production-mode fit. The total $\pm1\sigma$ uncertainty in each measurement is indicated by the error bar. The total systematic uncertainty in each measurement is also indicated.

The expected values for $\sigma_{H}\times B(H\rightarrow\tau\tau)$ relative to the SM expectations in the simplified template cross-section measurement. The total $\pm1\sigma$ uncertainty in each measurement is indicated by the error bar. The total systematic uncertainty in each measurement is also indicated.

The measured values for $\sigma_{H}\times B(H\rightarrow\tau\tau)$ in the simplified template cross-section measurement. The total $\pm1\sigma$ uncertainty in the measurement is indicated by the black error bars, with the individual contribution from the statistical uncertainty in blue.

The measured correlations between differential fiducial cross-sections, $\sigma$, and the normalisation factors, $\alpha$, for the $Z(\to\tau\tau) + \text{jets}$ and top backgrounds for $p_\text{T}(j_0)$. The bins denoted as $\sigma_{1,2,3,4}$ correspond to the measurements in the four bins listed elsewhere for each observable.

The measured correlations between differential fiducial cross-sections, $\sigma$, and the normalisation factors, $\alpha$, for the $Z(\to\tau\tau) + \text{jets}$ and top backgrounds for $p_\text{T}^{H}}$. The bins denoted as $\sigma_{1,2,3,4}$ correspond to the measurements in the four bins listed elsewhere for each observable.

The measured correlations between differential fiducial cross-sections, $\sigma$, and the normalisation factors, $\alpha$, for the $Z(\to\tau\tau) + \text{jets}$ and top backgrounds for $\Delta\phi_{jj}^{\text{signed}}$. The bins denoted as $\sigma_{1,2,3,4}$ correspond to the measurements in the four bins listed elsewhere for each observable.

The measured correlations between differential fiducial cross-sections, $\sigma$, and the normalisation factors, $\alpha$, for the $Z(\to\tau\tau) + \text{jets}$ and top backgrounds for $\Delta\phi_{jj}^{\text{signed}}$ vs $p_\text{T}^{H}$. The bins denoted as $\sigma_{1,2,3,4}$ correspond to the measurements in the four bins listed elsewhere for each observable.

Non-prompt $\mathrm{D}^0$ and $\mathrm{D}^+$ average $p_\mathrm{{T}}$-differential $R_{\mathrm{pPb}}$ in p--Pb collisions at $\sqrt{{s_\mathrm{NN}}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $-0.96 < y_{\mathrm{cms}} < 0.04$.

Non-prompt $\Lambda_{c}^{+}$ $p_\mathrm{{T}}$-differential $R_{\mathrm{pPb}}$ in p--Pb collisions at $\sqrt{{s_\mathrm{NN}}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $-0.96 < y_{\mathrm{cms}} < 0.04$.

Non-prompt $\mathrm{D}^0$ $p_\mathrm{{T}}$-integrated $R_{\mathrm{pPb}}$ in p--Pb collisions at $\sqrt{{s_\mathrm{NN}}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $-0.96 < y_{\mathrm{cms}} < 0.04$.

Non-prompt $\mathrm{D}^+$ $p_\mathrm{{T}}$-integrated $R_{\mathrm{pPb}}$ p--Pb collisions at $\sqrt{{s_\mathrm{NN}}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $-0.96 < y_{\mathrm{cms}} < 0.04$.

Non-prompt $\mathrm{D}^+$ / non-prompt $\mathrm{D}^0$ yield ratio in p--Pb collisions at $\sqrt{{s_\mathrm{NN}}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $-0.96 < y_{\mathrm{cms}} < 0.04$.

Non-prompt $\Lambda_{c}^{+}$ / non-prompt $\mathrm{D}^0$ yield ratio in p--Pb collisions at $\sqrt{{s_\mathrm{NN}}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $-0.96 < y_{\mathrm{cms}} < 0.04$.

Ratio of the $p_\mathrm{{T}}$-integrated ($2 < p_\mathrm{{T}} < 24$ $\mathrm{{GeV}}/c$) $R_{\mathrm{pPb}}$ of non-prompt $\Lambda_{c}^{+}$ and non-prompt $\mathrm{D}^0$ in p--Pb collisions at $\sqrt{{s_\mathrm{NN}}}=5.02~\mathrm{{TeV}}$ in the rapidity interval $-0.96 < y_{\mathrm{cms}} < 0.04$.

Transverse momentum spectra of antideuterons measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as shown in Fig. 1 (right panel). Rapidity interval 0.3 to 0.4.

Transverse momentum spectra of antideuterons measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as shown in Fig. 1 (right panel). Rapidity interval 0.4 to 0.5.

Transverse momentum spectra of antideuterons measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as shown in Fig. 1 (right panel). Rapidity interval 0.5 to 0.6.

Transverse momentum spectra of antideuterons measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as shown in Fig. 1 (right panel). Rapidity interval 0.6 to 0.7.

Transverse momentum spectra of antiprotons measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as shown in Fig. 1 (left panel). Rapidity interval 0 to 0.1.

Transverse momentum spectra of antiprotons measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as shown in Fig. 1 (left panel). Rapidity interval 0.1 to 0.2.

Transverse momentum spectra of antiprotons measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as shown in Fig. 1 (left panel). Rapidity interval 0.2 to 0.3.

Transverse momentum spectra of antiprotons measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as shown in Fig. 1 (left panel). Rapidity interval 0.3 to 0.4.

Transverse momentum spectra of antiprotons measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as shown in Fig. 1 (left panel). Rapidity interval 0.4 to 0.5.

Transverse momentum spectra of antiprotons measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as shown in Fig. 1 (left panel). Rapidity interval 0.5 to 0.6.

Transverse momentum spectra of antiprotons measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as shown in Fig. 1 (left panel). Rapidity interval 0.6 to 0.7.

Integrated yields of antideuterons measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as a function of absolute value of rapidity. The figure is shown in Fig. 2 (right panel).

Integrated yields of antiprotons measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as a function of absolute value of rapidity. The figure is shown in Fig. 2 (left panel).

coalescence parameter B2 measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as a function of transverse momentum per nucleon. The figure is shown in Fig. 3, rapidity interval 0.0 to 0.1, in absolute value.

coalescence parameter B2 measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as a function of transverse momentum per nucleon. The figure is shown in Fig. 3, rapidity interval 0.1 to 0.2, in absolute value.

coalescence parameter B2 measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as a function of transverse momentum per nucleon. The figure is shown in Fig. 3, rapidity interval 0.2 to 0.3, in absolute value.

coalescence parameter B2 measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as a function of transverse momentum per nucleon. The figure is shown in Fig. 3, rapidity interval 0.3 to 0.4, in absolute value.

coalescence parameter B2 measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as a function of transverse momentum per nucleon. The figure is shown in Fig. 3, rapidity interval 0.4 to 0.5, in absolute value.

coalescence parameter B2 measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as a function of transverse momentum per nucleon. The figure is shown in Fig. 3, rapidity interval 0.5 to 0.6, in absolute value.

coalescence parameter B2 measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as a function of transverse momentum per nucleon. The figure is shown in Fig. 3, rapidity interval 0.6 to 0.7, in absolute value.

coalescence parameter B2 measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as a function of rapidity, for fixed intervals of transverse momentum per nucleon. The data is shown in Fig. 4, in the PT per nucleon interval from 0.9 to 1.0.

coalescence parameter B2 measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as a function of rapidity, for fixed intervals of transverse momentum per nucleon. The data is shown in Fig. 4, in the PT per nucleon interval from 1.0 to 1.1.

coalescence parameter B2 measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as a function of rapidity, for fixed intervals of transverse momentum per nucleon. The data is shown in Fig. 4, in the PT per nucleon interval from 1.1 to 1.3.

coalescence parameter B2 measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as a function of rapidity, for fixed intervals of transverse momentum per nucleon. The data is shown in Fig. 4, in the PT per nucleon interval from 1.3 to 1.5.

coalescence parameter B2 measured in pp collisions at centre-of-mass per nucleon-nucleon energy of 13 TeV, as a function of rapidity, for fixed intervals of transverse momentum per nucleon. The data is shown in Fig. 4, in the PT per nucleon interval from 1.5 to 1.7.

Total covariance matrix for the unfolded Lund Jet Plane in the top jet selection.

Statistical covariance matrix for the unfolded Lund Jet Plane in the top jet selection.

Unfolded Lund Jet Plane in the W jet selection.

Unfolded Lund Jet Plane in the W jet selection, $2.11<\ln(1/z)<2.46$.

Unfolded Lund Jet Plane in the W jet selection, $1.60<\ln(R/\Delta R)<2.00$.

Total covariance matrix for the unfolded Lund Jet Plane in the W jet selection.

Statistical covariance matrix for the unfolded Lund Jet Plane in the W jet selection.

95&percnt; CL expected and observed exclusion limits on the branching ratio of the Higgs boson to dark photons as a function of the &gamma;_d mass for the HAHM model.

95&percnt; CL expected and observed exclusion limits on the branching ratio of the Higgs boson to dark photons as a function of the &gamma;_d mass for the FRVZ model.

90&percnt; CL exclusion limit contours on the branching ratio of the Higgs boson to dark photons as functions of the &gamma;_d mass and kinetic mixing parameter &epsilon; for the HAHM model. These exclusions assume Higgs boson branching fractions to dark photons within a 0.01&percnt; and 10&percnt; range.

90&percnt; CL exclusion limit contours on the branching ratio of the Higgs boson to dark photons as functions of the &gamma;_d mass and kinetic mixing parameter &epsilon; for the FRVZ model. These exclusions assume Higgs boson branching fractions to dark photons within a 0.01&percnt; and 10&percnt; range.

The reconstruction efficiency for &mu;LJ objects produced in the decay of a &gamma;_d into a muon pair. The reconstruction efficiency is shown for different dark photon mass hypotheses generated according to the HAHM signal model. The efficiency is shown as a function of the dark photon p_T.

The reconstruction efficiency for &mu;LJ objects produced in the decay of a &gamma;_d into a muon pair. The reconstruction efficiency is shown for different dark photon mass hypotheses generated according to the FRVZ signal model. The efficiency is shown as a function of the dark photon p_T.

The reconstruction efficiency for &mu;LJ objects produced in the decay of a &gamma;_d into a muon pair. The reconstruction efficiency is shown for different dark photon mass hypotheses generated according to the HAHM signal model. The efficiency is shown as a function of the &Delta;R opening between the decay products of a dark photon.

The reconstruction efficiency for &mu;LJ objects produced in the decay of a &gamma;_d into a muon pair. The reconstruction efficiency is shown for different dark photon mass hypotheses generated according to the FRVZ signal model. The efficiency is shown as a function of the &Delta;R opening between the decay products of a dark photon.

The reconstruction efficiency for &mu;LJ objects produced in the decay of a &gamma;_d into a muon pair. The reconstruction efficiency is shown for a dark photon mass of 0.4 GeV, generated according to the HAHM and FRVZ signal models. The efficiency is shown as a function of the &gamma;_d p_T.

The reconstruction efficiency for &mu;LJ objects produced in the decay of a &gamma;_d into a muon pair. The reconstruction efficiency is shown for a dark photon mass of 0.4 GeV, generated according to the HAHM and FRVZ signal models. The efficiency is shown as a function of the &Delta;R opening between the decay products of a dark photon.

The reconstruction efficiency for eLJ objects produced in the decay of a &gamma;_d into an electron pair. The reconstruction efficiency is shown for different dark photon mass hypotheses generated according to the HAHM signal model. The efficiency is shown as a function of the dark photon p_T.

The reconstruction efficiency for eLJ objects produced in the decay of a &gamma;_d into an electron pair. The reconstruction efficiency is shown for different dark photon mass hypotheses generated according to the FRVZ signal model. The efficiency is shown as a function of the dark photon p_T.

The reconstruction efficiency for eLJ objects produced in the decay of a &gamma;_d into a pair of electrons. The reconstruction efficiency is shown for different dark photon mass hypotheses generated according to the HAHM signal model. The efficiency is shown as a function of the &Delta;R opening between the decay products of a dark photon. The decrease in efficiency, visible at 0.1, is due to the transition between the region of the phase space where eLJs are identified from merged EM showers with two tracks, to the one where eLJs are obtained by the clustering of two resolved electrons.

The reconstruction efficiency for eLJ objects produced in the decay of a &gamma;_d into a pair of electrons. The reconstruction efficiency is shown for different dark photon mass hypotheses generated according to the FRVZ signal model. The efficiency is shown as a function of the &Delta;R opening between the decay products of a dark photon. The decrease in efficiency, visible at 0.1, is due to the transition between the region of the phase space where eLJs are identified from merged EM showers with two tracks, to the one where eLJs are obtained by the clustering of two resolved electrons.

The reconstruction efficiency for eLJ objects produced in the decay of a &gamma;_d into a pair of electrons. The reconstruction efficiency is shown for a dark photon mass of 0.1 GeV, generated according to the HAHM and FRVZ signal models. The efficiency is shown as a function of the &gamma;_d p_T.

The reconstruction efficiency for eLJ objects produced in the decay of a &gamma;_d into a pair of electrons. The reconstruction efficiency is shown for a dark photon mass of 0.1 GeV, generated according to the HAHM and FRVZ signal models. The efficiency is shown as a function of the &Delta;R opening between the decay products of a dark photon.

95&percnt; CL expected and observed exclusion limits on the branching ratio of the Higgs boson to dark photons as a function of the &gamma;_d mass for the HAHM model, obtained from the &mu;LJ&ndash;&mu;LJ signal region.

95&percnt; CL expected and observed exclusion limits on the branching ratio of the Higgs boson to dark photons as a function of the &gamma;_d mass for the FRVZ model, obtained from the &mu;LJ&ndash;&mu;LJ signal region.

95&percnt; CL expected and observed exclusion limits on the branching ratio of the Higgs boson to dark photons as a function of the &gamma;_d mass for the HAHM model, obtained from the &mu;LJ&ndash;eLJ signal region.

95&percnt; CL expected and observed exclusion limits on the branching ratio of the Higgs boson to dark photons as a function of the &gamma;_d mass for the FRVZ model, obtained from the &mu;LJ&ndash;eLJ signal region.

Acceptance times efficiencies (A&times;&varepsilon;) as a function of the &gamma;_d mass for the HAHM process for the three signal regions of the analysis, eLJ-eLJ, eLJ-&mu;LJ and &mu;LJ-&mu;LJ. The decrease of the A&times;&varepsilon; is expected for increasing masses of the &gamma;_d, as its decay products become resolved and are no longer reconstructed as LJs. The local minimum around m_{&gamma;_d}=2 GeV observed for the eLJ-&mu;LJ distribution is expected, as the reconstruction efficiency for eLJs with merged EM showers decreases before the reconstruction of eLJ with two resolved electrons becomes efficient.

Acceptance times efficiencies (A&times;&varepsilon;) as a function of the &gamma;_d mass for the FRVZ process for the three signal regions of the analysis, eLJ-eLJ, eLJ-&mu;LJ and &mu;LJ-&mu;LJ. The decrease of the A&times;&varepsilon; is expected for increasing masses of the &gamma;_d, as its decay products become resolved and are no longer reconstructed as LJs. The local minimum around m_{&gamma;_d}=2 GeV observed for the eLJ-&mu;LJ distribution is expected, as the reconstruction efficiency for eLJs with merged EM showers decreases before the reconstruction of eLJ with two resolved electrons becomes efficient.

The transverse-mass dependence of the correlation-strength parameter $\lambda$ in 30-40% centrality bin obtained from Lévy fits with Eq. (9).

The transverse-mass dependence of the correlation-strength parameter $\lambda$ in 40-50% centrality bin obtained from Lévy fits with Eq. (9).

The transverse-mass dependence of the correlation-strength parameter $\lambda$ in 50-60% centrality bin obtained from Lévy fits with Eq. (9).

The transverse-mass dependence of the Lévy-scale parameter $R$ in 0-10% centrality bin obtained from Lévy fits with Eq. (9).

The transverse-mass dependence of the Lévy-scale parameter $R$ in 10-20% centrality bin obtained from Lévy fits with Eq. (9).

The transverse-mass dependence of the Lévy-scale parameter $R$ in 20-30% centrality bin obtained from Lévy fits with Eq. (9).

The transverse-mass dependence of the Lévy-scale parameter $R$ in 30-40% centrality bin obtained from Lévy fits with Eq. (9).

The transverse-mass dependence of the Lévy-scale parameter $R$ in 40-50% centrality bin obtained from Lévy fits with Eq. (9).

The transverse-mass dependence of the Lévy-scale parameter $R$ in 50-60% centrality bin obtained from Lévy fits with Eq. (9).

The transverse-mass dependence of the Lévy-index of stability parameter $\alpha$ in 0-10% centrality bin obtained from Lévy fits with Eq. (9).

The transverse-mass dependence of the Lévy-index of stability parameter $\alpha$ in 10-20% centrality bin obtained from Lévy fits with Eq. (9).

The transverse-mass dependence of the Lévy-index of stability parameter $\alpha$ in 20-30% centrality bin obtained from Lévy fits with Eq. (9).

The transverse-mass dependence of the Lévy-index of stability parameter $\alpha$ in 30-40% centrality bin obtained from Lévy fits with Eq. (9).

The transverse-mass dependence of the Lévy-index of stability parameter $\alpha$ in 40-50% centrality bin obtained from Lévy fits with Eq. (9).

The transverse-mass dependence of the Lévy-index of stability parameter $\alpha$ in 50-60% centrality bin obtained from Lévy fits with Eq. (9).

The transverse-mass dependence of the normalized correlation-strength parameter $\lambda/\lambda_{\rm max}$ for 0-10% centrality intervals. For each centrality bin the data are fitted with the Gaussian function of Eq. (15). This parameterization can describe the data and is discussed in detail in Section VI B. The fit parameters with the statistic and systematic uncertainties are shown in Fig. 7.

The transverse-mass dependence of the normalized correlation-strength parameter $\lambda/\lambda_{\rm max}$ for 10-20% centrality intervals. For each centrality bin the data are fitted with the Gaussian function of Eq. (15). This parameterization can describe the data and is discussed in detail in Section VI B. The fit parameters with the statistic and systematic uncertainties are shown in Fig. 7.

The transverse-mass dependence of the normalized correlation-strength parameter $\lambda/\lambda_{\rm max}$ for 20-30% centrality intervals. For each centrality bin the data are fitted with the Gaussian function of Eq. (15). This parameterization can describe the data and is discussed in detail in Section VI B. The fit parameters with the statistic and systematic uncertainties are shown in Fig. 7.

The transverse-mass dependence of the normalized correlation-strength parameter $\lambda/\lambda_{\rm max}$ for 30-40% centrality intervals. For each centrality bin the data are fitted with the Gaussian function of Eq. (15). This parameterization can describe the data and is discussed in detail in Section VI B. The fit parameters with the statistic and systematic uncertainties are shown in Fig. 7.

The transverse-mass dependence of the normalized correlation-strength parameter $\lambda/\lambda_{\rm max}$ for 40-50% centrality interval. For each centrality bin the data are fitted with the Gaussian function of Eq. (15). This parameterization can describe the data and is discussed in detail in Section VI B. The fit parameters with the statistic and systematic uncertainties are shown in Fig. 7.

The transverse-mass dependence of the normalized correlation-strength parameter $\lambda/\lambda_{\rm max}$ for 50-60% centrality interval. For each centrality bin the data are fitted with the Gaussian function of Eq. (15). This parameterization can describe the data and is discussed in detail in Section VI B. The fit parameters with the statistic and systematic uncertainties are shown in Fig. 7.

The transverse-mass dependence of the inverse square of the Lévy-scale parameter $R$, shown in 0-10% centrality bin. The values and uncertainties of the A and B linear fit parameters with the statistic and systematic uncertainties are shown in Fig. 8.

The transverse-mass dependence of the inverse square of the Lévy-scale parameter $R$, shown in 10-20% centrality bin. The values and uncertainties of the A and B linear fit parameters with the statistic and systematic uncertainties are shown in Fig. 8.

The transverse-mass dependence of the inverse square of the Lévy-scale parameter $R$, shown in 20-30% centrality bin. The values and uncertainties of the A and B linear fit parameters with the statistic and systematic uncertainties are shown in Fig. 8.

The transverse-mass dependence of the inverse square of the Lévy-scale parameter $R$, shown in 30-40% centrality bin. The values and uncertainties of the A and B linear fit parameters with the statistic and systematic uncertainties are shown in Fig. 8.

The transverse-mass dependence of the inverse square of the Lévy-scale parameter $R$, shown in 40-50% centrality bin. The values and uncertainties of the A and B linear fit parameters with the statistic and systematic uncertainties are shown in Fig. 8.

The transverse-mass dependence of the inverse square of the Lévy-scale parameter $R$, shown in 50-60% centrality bin. The values and uncertainties of the A and B linear fit parameters with the statistic and systematic uncertainties are shown in Fig. 8.

The transverse-mass dependence of the $\widehat{R}$ parameter, shown in 0-10% centrality bin. The values and uncertainties of the linear fit parameters $\widehat{A}$ and $\widehat{B}$ with the statistic and systematic uncertainties are shown in Fig. 9

The transverse-mass dependence of the $\widehat{R}$ parameter, shown in 10-20% centrality bin. The values and uncertainties of the linear fit parameters $\widehat{A}$ and $\widehat{B}$ with the statistic and systematic uncertainties are shown in Fig. 9

The transverse-mass dependence of the $\widehat{R}$ parameter, shown in 20-30% centrality bin. The values and uncertainties of the linear fit parameters $\widehat{A}$ and $\widehat{B}$ with the statistic and systematic uncertainties are shown in Fig. 9

The transverse-mass dependence of the $\widehat{R}$ parameter, shown in 30-40% centrality bin. The values and uncertainties of the linear fit parameters $\widehat{A}$ and $\widehat{B}$ with the statistic and systematic uncertainties are shown in Fig. 9

The transverse-mass dependence of the $\widehat{R}$ parameter, shown in 40-50% centrality bin. The values and uncertainties of the linear fit parameters $\widehat{A}$ and $\widehat{B}$ with the statistic and systematic uncertainties are shown in Fig. 9

The transverse-mass dependence of the $\widehat{R}$ parameter, shown in 50-60% centrality bin. The values and uncertainties of the linear fit parameters $\widehat{A}$ and $\widehat{B}$ with the statistic and systematic uncertainties are shown in Fig. 9

The $H$ parameter that characterize the magnitude of the $\lambda/\lambda_{max}(m_T)$ and is defined in Eq. (15)., shown as functions of $N_{part}$.

The $\sigma$ parameter that characterize the width of the $\lambda/\lambda_{max}(m_T)$ and is defined in Eq. (15)., shown as functions of $N_{part}$.

The parameter $A$ of the affine linear fit to the inverse square of the Lévy-scale parameter that are defined in Eq.(12) as function of $N_{part}$

The parameter $B$ of the affine linear fit to the inverse square of the Lévy-scale parameter that are defined in Eq.(12) as function of $N_{part}$

The parameter $\widehat{A}$ of the affine linear fit to the inverse of the $\widehat{R}$ parameter that are defined in Eq.(14) as function of $N_{part}$.

The parameter $\widehat{B}$ of the affine linear fit to the inverse of the $\widehat{R}$ parameter that are defined in Eq.(14) as function of $N_{part}$.

The centrality dependence of $\alpha_0$, the $m_T$-averaged value of $\alpha$ as a function of $N_{part}$.

The Lévy scale parameter $R$ as a function of $N^{1/3}_{part}$ at $m_T=0.326$ GeV.

The Lévy scale parameter $R$ as a function of $N^{1/3}_{part}$ at $m_T=0.438$ GeV.

The Lévy scale parameter $R$ as a function of $N^{1/3}_{part}$ at $m_T=0.553$ GeV.

The Lévy scale parameter $R$ as a function of $N^{1/3}_{part}$ at $m_T=0.670$ GeV.

The Lévy scale parameter $R$ as a function of $N^{1/3}_{part}$ at $m_T=0.787$ GeV.

Monte Carlo calcuations with no mass drop assumed in 0-10%.

Monte Carlo calcuations with no mass drop assumed in 10-20%.

Monte Carlo calcuations with no mass drop assumed in 20-30%.

Monte Carlo calcuations with no mass drop assumed in 30-40%.

Monte Carlo calcuations with no mass drop assumed in 40-50%.

Monte Carlo calcuations with no mass drop assumed in 50-60%.

Monte Carlo calcuations with the best mass drop assumed in 0-10%.

Monte Carlo calcuations with the best mass drop assumed in 10-20%.

Monte Carlo calcuations with the best mass drop assumed in 20-30%.

Monte Carlo calcuations with the best mass drop assumed in 30-40%.

Monte Carlo calcuations with the best mass drop assumed in 40-50%.

Monte Carlo calcuations with the best mass drop assumed in 50-60%.

The best values of $B^{-1}_{\eta\prime}$ spectrum of the $\eta\prime$ condensate in the six centrality classes.

The centrality dependence of the best values of the in-medium mass of the $m^{*}_{\eta\prime}$.

The CL maps of the comparison of the MC to data in 0-10% (see Fig. 11).

The CL maps of the comparison of the MC to data in 10-20% (see Fig. 11).

The CL maps of the comparison of the MC to data in 20-30% (see Fig. 11).

The CL maps of the comparison of the MC to data in 30-40% (see Fig. 11).

The CL maps of the comparison of the MC to data in 40-50% (see Fig. 11).

The CL maps of the comparison of the MC to data in 50-60% (see Fig. 11).

The isolated-photon cross section ratio, 13 TeV to 7 TeV, from pQCD NLO calculations with JETPHOX. The calculations were obtained by choosing factorisation, normalisation, and fragmentation scales. The parton distribution function used in the calculations is NNPDF4.0 for the $\sqrt{s}=$13 TeV measurement and CT14 for the $\sqrt{s}=$7 TeV measurement. The fragmentation function is BFG II for both measurements.

Differential cross section of isolated photons as a function of $x_\mathrm{T}^{\gamma}$ measured in pp collisions at $\sqrt{s}=$13 TeV.