Charged-particle multiplicities in proton-proton collisions at a centre-of-mass energy of 13000 GeV as a function of transverse momentum for events with the number of charged particles >=1 having transverse momentum >500 MeV and absolute(pseudorapidity) <2.5.

Charged-particle multiplicity distribution in proton-proton collisions at a centre-of-mass energy of 13000 GeV for events with the number of charged particles >=1 having transverse momentum >500 MeV and absolute(pseudorapidity) <2.5.

Average transverse momentum in proton-proton collisions at a centre-of-mass energy of 13000 GeV as a function of the number of charged particles in the event for events with the number of charged particles >=1 having transverse momentum >500 MeV and absolute(pseudorapidity) <2.5.

Extrapolated charged-particle multiplicities in proton-proton collisions at a centre-of-mass energy of 13000 GeV as a function of pseudorapidity for events with the number of charged particles >=1 having transverse momentum >500 MeV and absolute(pseudorapidity) <2.5.

Extrapolated charged-particle multiplicities in proton-proton collisions at a centre-of-mass energy of 13000 GeV as a function of transverse momentum for events with the number of charged particles >=1 having transverse momentum >500 MeV and absolute(pseudorapidity) <2.5.

Extrapolated charged-particle multiplicity distributions in proton-proton collisions at a centre-of-mass energy of 13000 GeV for events with the number of charged particles >=1 having transverse momentum >500 MeV and absolute(pseudorapidity) <2.5.

Extrapolated average transverse momentum in proton-proton collisions at a centre-of-mass energy of 13000 GeV as a function of the number of charged particles in the event for events with the number of charged particles >=1 having transverse momentum >500 MeV and absolute(pseudorapidity) <2.5.

Charged-particle multiplicities in proton-proton collisions at a centre-of-mass energy of 13000 GeV as a function of pseudorapidity for events with the number of charged particles >=1 having transverse momentum >500 MeV and absolute(pseudorapidity) <0.8.

Charged-particle multiplicities in proton-proton collisions at a centre-of-mass energy of 13000 GeV as a function of transverse momentum for events with the number of charged particles >=1 having transverse momentum >500 MeV and absolute(pseudorapidity) <0.8.

Charged-particle multiplicity distribution in proton-proton collisions at a centre-of-mass energy of 13000 GeV for events with the number of charged particles >=1 having transverse momentum >500 MeV and absolute(pseudorapidity) <0.8.

Average transverse momentum in proton-proton collisions at a centre-of-mass energy of 13000 GeV as a function of the number of charged particles in the event for events with the number of charged particles >=1 having transverse momentum >500 MeV and absolute(pseudorapidity) <0.8.

Extrapolated charged-particle multiplicities in proton-proton collisions at a centre-of-mass energy of 13000 GeV as a function of pseudorapidity for events with the number of charged particles >=1 having transverse momentum >500 MeV and absolute(pseudorapidity) <0.8.

Extrapolated charged-particle multiplicities in proton-proton collisions at a centre-of-mass energy of 13000 GeV as a function of transverse momentum for events with the number of charged particles >=1 having transverse momentum >500 MeV and absolute(pseudorapidity) <0.8.

Extrapolated charged-particle multiplicity distributions in proton-proton collisions at a centre-of-mass energy of 13000 GeV for events with the number of charged particles >=1 having transverse momentum >500 MeV and absolute(pseudorapidity) <0.8.

Extrapolated average transverse momentum in proton-proton collisions at a centre-of-mass energy of 13000 GeV as a function of the number of charged particles in the event for events with the number of charged particles >=1 having transverse momentum >500 MeV and absolute(pseudorapidity) <0.8.

Measured fiducial cross section in all $\ell^+\ell^-\ell'^+\ell'^-$ channels combined, where $\ell, \ell' = e, \mu$. The first systematic uncertainty is the combined systematic uncertainty excluding luminosity uncertainty, the second is the luminosity uncertainty.

Cross section extrapolated to total phase space and all Z boson decays. The first systematic uncertainty is the combined systematic uncertainty excluding luminosity uncertainty, the second is the luminosity uncertainty.

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 1.0 $<p_{T}<$ 2.0 $GeV/c$ and $N_{offline}^{trk}<$ 35 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 2.0 $<p_{T}<$ 3.0 $GeV/c$ and $N_{offline}^{trk}<$ 35 bins for pp data at $\sqrt =$ 13 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 2.0 $<p_{T}<$ 3.0 $GeV/c$ and $N_{offline}^{trk}<$ 35 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 3.0 $<p_{T}<$ 4.0 $GeV/c$ and $N_{offline}^{trk}<$ 35 bins for pp data at $\sqrt =$ 13 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 3.0 $<p_{T}<$ 4.0 $GeV/c$ and $N_{offline}^{trk}<$ 35 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 0.1 $<p_{T}<$ 1.0 $GeV/c$ and 35 $<N_{offline}^{trk}<$ 80 bins for pp data at $\sqrt =$ 13 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 0.1 $<p_{T}<$ 1.0 $GeV/c$ and 35 $<N_{offline}^{trk}<$ 90 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 1.0 $<p_{T}<$ 2.0 $GeV/c$ and 35 $<N_{offline}^{trk}<$ 80 bins for pp data at $\sqrt =$ 13 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 1.0 $<p_{T}<$ 2.0 $GeV/c$ and 35 $<N_{offline}^{trk}<$ 90 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 2.0 $<p_{T}<$ 3.0 $GeV/c$ and 35 $<N_{offline}^{trk}<$ 80 bins for pp data at $\sqrt =$ 13 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 2.0 $<p_{T}<$ 3.0 $GeV/c$ and 35 $<N_{offline}^{trk}<$ 90 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 3.0 $<p_{T}<$ 4.0 $GeV/c$ and 35 $<N_{offline}^{trk}<$ 80 bins for pp data at $\sqrt =$ 13 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 3.0 $<p_{T}<$ 4.0 $GeV/c$ and 35 $<N_{offline}^{trk}<$ 90 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 0.1 $<p_{T}<$ 1.0 $GeV/c$ and 80 $<N_{offline}^{trk}<$ 105 bins for pp data at $\sqrt =$ 13 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 0.1 $<p_{T}<$ 1.0 $GeV/c$ and 90 $<N_{offline}^{trk}<$ 105 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 1.0 $<p_{T}<$ 2.0 $GeV/c$ and 80 $<N_{offline}^{trk}<$ 105 bins for pp data at $\sqrt =$ 13 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 1.0 $<p_{T}<$ 2.0 $GeV/c$ and 90 $<N_{offline}^{trk}<$ 105 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 2.0 $<p_{T}<$ 3.0 $GeV/c$ and 80 $<N_{offline}^{trk}<$ 105 bins for pp data at $\sqrt =$ 13 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 2.0 $<p_{T}<$ 3.0 $GeV/c$ and 90 $<N_{offline}^{trk}<$ 105 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 3.0 $<p_{T}<$ 4.0 $GeV/c$ and 80 $<N_{offline}^{trk}<$ 105 bins for pp data at $\sqrt =$ 13 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 3.0 $<p_{T}<$ 4.0 $GeV/c$ and 90 $<N_{offline}^{trk}<$ 105 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 0.1 $<p_{T}<$ 1.0 $GeV/c$ and $N_{offline}^{trk}>$ 105 bins for pp data at $\sqrt =$ 13 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 0.1 $<p_{T}<$ 1.0 $GeV/c$ and $N_{offline}^{trk}>$ 110 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 1.0 $<p_{T}<$ 2.0 $GeV/c$ and $N_{offline}^{trk}>$ 105 bins for pp data at $\sqrt =$ 13 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 1.0 $<p_{T}<$ 2.0 $GeV/c$ and $N_{offline}^{trk}>$ 110 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 2.0 $<p_{T}<$ 3.0 $GeV/c$ and $N_{offline}^{trk}>$ 105 bins for pp data at $\sqrt =$ 13 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 2.0 $<p_{T}<$ 3.0 $GeV/c$ and $N_{offline}^{trk}>$ 110 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 3.0 $<p_{T}<$ 4.0 $GeV/c$ and $N_{offline}^{trk}>$ 105 bins for pp data at $\sqrt =$ 13 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Correlated yield obtained with the ZYAM procedure as a function of $|\Delta\Phi|$, averaged over 2 $<|\Delta\eta|<$ 4 in for 3.0 $<p_{T}<$ 4.0 $GeV/c$ and $N_{offline}^{trk}>$ 110 bins for pp data at $\sqrt =$ 7 $TeV$. The $p_{T}$ selection applies to both particles in the pair. Only statistical uncertainties are given. The subtracted ZYAM constant is given ($C_{ZYAM}$).

Associated yield for the near side of the correlation function averaged over 2 $<|\Delta\eta|<$ 4 and integrated over the region $|\Delta\Phi| < \Delta\Phi_{ZYAM}$ for pp data at $\sqrt{s} =$ 13 $TeV$. The associated yield as a function of $p_{T}$ for events with $N^{offline}_{trk} \geq$ 105. The $p_{T}$ value for each $p_{T}$ bin is the average $p_{T}$ value.

Associated yield for the near side of the correlation function averaged over 2 $<|\Delta\eta|<$ 4 for pp data at $\sqrt{s} =$ 7 $TeV$. The associated yield as a function of $p_{T}$ for events with $N^{offline}_{trk} \geq$ 110. The $p_{T}$ value for each $p_{T}$ bin is the average $p_{T}$ value.

Associated yield for the near side of the correlation function averaged over 2 $<|\Delta\eta|<$ 4 and integrated over the region $|\Delta\Phi| < \Delta\Phi_{ZYAM}$ for pp data at $\sqrt{s} =$ 13 $TeV$. The associated yield as a function of $N_{trk}^{offline}$ for events with 1.0 $< p_{T} <$ 2.0 GeV/c. The $N_{trk}^{offline}$ value for each $N_{trk}^{offline}$ bin is the average $N_{trk}^{offline}$ value.

Associated yield for the near side of the correlation function averaged over 2 $<|\Delta\eta|<$ 4 for pp data at $\sqrt{s} =$ 7 $TeV$. The associated yield as a function of $N_{trk}^{offline}$ for events with 1.0 $< p_{T} <$ 2.0 GeV/c. The $N_{trk}^{offline}$ value for each $N_{trk}^{offline}$ bin is the average $N_{trk}^{offline}$ value.

Associated yield for the near side of the correlation function averaged over 2 $<|\Delta\eta|<$ 4 pPb data at $\sqrt{s} =$ 5.02 $TeV$. The associated yield as a function of $N_{trk}^{offline}$ for events with 1.0 $< p_{T} <$ 2.0 GeV/c. The $N_{trk}^{offline}$ value for each $N_{trk}^{offline}$ bin is the average $N_{trk}^{offline}$ value.

Associated yield for the near side of the correlation function averaged over 2 $<|\Delta\eta|<$ 4 PbPb data at $\sqrt{s} =$ 2.76 $TeV$. The associated yield as a function of $N_{trk}^{offline}$ for events with 1.0 $< p_{T} <$ 2.0 GeV/c. The $N_{trk}^{offline}$ value for each $N_{trk}^{offline}$ bin is the average $N_{trk}^{offline}$ value.

Differential production cross-sections for prompt $D^{*+} + D^{*-}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic.

The ratios of differential production cross-sections, $R_{13/7}$, for prompt $D^{0} + \bar{D}^{0}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic.

The ratios of differential production cross-sections, $R_{13/7}$, for prompt $D^{+} + D^{-}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic.

The ratios of differential production cross-sections, $R_{13/7}$, for prompt $D_{s}^{+} + D_{s}^{-}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic.

The ratios of differential production cross-sections, $R_{13/7}$, for prompt $D^{*+} + D^{*-}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic.

The ratios of differential production cross-section-times-branching fraction measurements for prompt $D^{+}$ and $D^{0}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic. All values are given in percent.

The ratios of differential production cross-section-times-branching fraction measurements for prompt $D_{s}^{+}$ and $D^{0}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic. All values are given in percent.

The ratios of differential production cross-section-times-branching fraction measurements for prompt $D^{*+}$ and $D^{0}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic. All values are given in percent.

The ratios of differential production cross-section-times-branching fraction measurements for prompt $D_{s}^{+}$ and $D^{+}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic. All values are given in percent.

The ratios of differential production cross-section-times-branching fraction measurements for prompt $D^{*+}$ and $D^{+}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic. All values are given in percent.

The ratios of differential production cross-section-times-branching fraction measurements for prompt $D_{s}^{+}$ and $D^{*+}$ mesons in bins of $(p_{\mathrm{T}}, y)$. The first uncertainty is statistical, and the second is the total systematic. All values are given in percent.

Ratios of transverse-momentum distributions of charged particles in three intervals of multiplicities to the respective one for inclusive (INEL>0) collisions. The spectra were normalised by the integral prior to division.

Differential cross-sections as a function of $p_\perp$ integrated over $y$ for prompt $J/\psi$ mesons. The first uncertainties are statistical and the second (third) are uncorrelated (correlated) systematic uncertainties amongst bins.

Differential cross-sections as a function of $p_\perp$ integrated over $y$ for $J/\psi$-from-$b$ mesons. The first uncertainties are statistical and the second (third) are uncorrelated (correlated) systematic uncertainties amongst bins.

Differential cross-sections as a function of $y$ integrated over $p_\perp$ for prompt $J/\psi$ mesons. The first uncertainties are statistical and the second (third) are the correlated (uncorrelated) systematic uncertainties.

Differential cross-sections as a function of $y$ integrated over $p_\perp$ for $J/\psi$-from-$b$ mesons. The first uncertainties are statistical and the second (third) are the correlated (uncorrelated) systematic uncertainties.

Ratios of differential cross-sections between measurements at $\sqrt{s} = 13$ TeV and $\sqrt{s} = 8$ TeV as a function of $p_\perp$ in bins of $y$ for prompt $J/\psi$ mesons. The systematic errors are negligible.

Ratios of differential cross-sections between measurements at $\sqrt{s} = 13$ TeV and $\sqrt{s} = 8$ TeV as a function of $p_\perp$ in bins of $y$ for $J/\psi$-from-$b$ mesons. The systematic errors are negligible.

Ratios of differential cross-sections between measurements at $\sqrt{s} = 13$ TeV and $\sqrt{s} = 8$ TeV as a function of $y$ integrated over $p_\perp$ for prompt $J/\psi$ mesons. The systematic errors are negligible.

Ratios of differential cross-sections between measurements at $\sqrt{s} = 13$ TeV and $\sqrt{s} = 8$ TeV as a function of $y$ integrated over $p_\perp$ for $J/\psi$-from-$b$ mesons. The systematic errors are negligible.

Ratios of differential cross-sections between measurements at $\sqrt{s} = 13$ TeV and $\sqrt{s} = 8$ TeV as a function of $p_\perp$ integrated over $y$ for prompt $J/\psi$ mesons.

Ratios of differential cross-sections between measurements at $\sqrt{s} = 13$ TeV and $\sqrt{s} = 8$ TeV as a function of $p_\perp$ integrated over $y$ for $J/\psi$-from-$b$ mesons. The systematic errors are negligible.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

No description provided.

Summary of the impact of different categories of systematic uncertainties on the background estimate in BDT bins, after the background-only fit. The individual uncertainties can be correlated and do not necessarily add in quadrature to equal the total background uncertainty.

Observed and expected exclusion limits at 95% CL on BR(H$\to\gamma\gamma_d$) as function of the $\gamma_d$ mass.

Background, data and signal yields in bins of BDT, in the SR and in the VV$\gamma$ CR, before the background-only fit. The expectations for the signal are shown for the massless dark photon and for dark photon mass values of 20 GeV and 40 GeV, assuming BR(H$\to\gamma\gamma_d$) = 5%. Uncertainties are statistical for data, while for backgrounds include statistical and systematic sources.

Number of events after each cut in SR for 139fb$^{-1}$ for ZH, H$\to\gamma\gamma_d$ considering BR=5% and massless dark photon. Generator-level filters are applied on the total events. Only the statistical uncertainty is shown.

Distribution of the BDT classifier response in data and for the expected SM background after the background-only fit, in the ee channel. The expectations for signal are also shown for the massless dark photon and for dark photon mass values of 20 GeV and 40 GeV, assuming BR(H$\to\gamma\gamma_d$) = 5%. Uncertainties shown are statistical for data, while for backgrounds include statistical and systematic sources.

Distribution of the BDT classifier response for data and for the expected SM background after the background-only fit, in the $\mu\mu$ channel. The expectations for signal are also shown for the massless dark photon and for dark photon mass values of 20 GeV and 40 GeV, assuming BR(H$\to\gamma\gamma_d$) = 5%. Uncertainties shown are statistical for data, while for backgrounds include statistical and systematic sources.

Distribution of jet girth of CMS data, SM multijet MC events, and various signal samples. Events are required to pass the trigger requirements and also have 4 jets with $p_\mathrm{T}>100\mathrm{GeV}$. The distribution from various processes have their total count normalized to 1.

Distribution of uGNN output score of CMS data, SM multijet MC events, and various signal samples. Events are required to pass the trigger requirements and also have 4 jets with $p_\mathrm{T}>100\mathrm{GeV}$. The distribution from various processes have their total count normalized to 1.

Distribution of aGNN output score of CMS data, SM multijet MC events, and various signal samples. Events are required to pass the trigger requirements and also have 4 jets with $p_\mathrm{T}>100\mathrm{GeV}$. The distribution from various processes have their total count normalized to 1.

Template fit of the DeepJet discriminator used to determine the b jet fraction of the non-EMJ tagged jets in data events that pass the 1-EMJ selection of the event selection criteria defined in the paper.

The EJ tagger misidentification probability for b quark jets and light jets as a function of jet $p_\mathrm{T}$ for the taggers u-tag 1, u-tag 2, and u-tag 3, as defined in the paper, evaluated using data and generator-level flavor information from simulated samples in events containing a high-$p_\mathrm{T}$ photon.

The EJ tagger misidentification probability for b quark jets and light jets as a function of jet $p_\mathrm{T}$ for the taggers uGNN tag 1, uGNN tag 2, and uGNN tag 3, as defined in the paper, evaluated using data and generator-level flavor information from simulated samples in events containing a high-$p_\mathrm{T}$ photon.

The EJ tagger misidentification probability for b quark jets and light jets as a function of jet $p_\mathrm{T}$ for the taggers u-tag 4 and u-tag 5, as defined in the paper, evaluated using data and generator-level flavor information from simulated samples in events containing a high-$p_\mathrm{T}$ photon.

The expected and observed 95% CL upper limits on the production cross section for various signal models in the unflavored scenario model-generic and GNN EMJ tagging methods.

The expected and observed 95% CL upper limits on the production cross section for various signal models in the flavor-aligned scenario model-generic and GNN EMJ tagging methods.

The observed yield of events in data satisfying the validation selection criteria with at least two jets passing the corresponding validation tag, and the estimation based on the misidentification rate calculated using validation events with exactly one jet passing the validation tagger scaled by the factor given in Equation (4) of the paper.

The estimated number of events from the background prediction based on control samples in data and the observed event yields.

Distribution of event $H_\mathrm{T}$ of CMS data, SM multijet MC events, and various signal samples. Events are required to pass the trigger requirements and also have 4 jets with $p_\mathrm{T}>100\mathrm{GeV}$. The distribution from various processes have their total count normalized to 1.

Distribution of leading jet $p_\mathrm{T}$ of CMS data, SM multijet MC events, and various signal samples. Events are required to pass the trigger requirements and also have 4 jets with $p_\mathrm{T}>100\mathrm{GeV}$. The distribution from various processes have their total count normalized to 1.

Distribution of second leading jet $p_\mathrm{T}$ of CMS data, SM multijet MC events, and various signal samples. Events are required to pass the trigger requirements and also have 4 jets with $p_\mathrm{T}>100\mathrm{GeV}$. The distribution from various processes have their total count normalized to 1.

Distribution of third leading jet $p_\mathrm{T}$ of CMS data, SM multijet MC events, and various signal samples. Events are required to pass the trigger requirements and also have 4 jets with $p_\mathrm{T}>100\mathrm{GeV}$. The distribution from various processes have their total count normalized to 1.

Distribution of forth leading jet $p_\mathrm{T}$ of CMS data, SM multijet MC events, and various signal samples. Events are required to pass the trigger requirements and also have 4 jets with $p_\mathrm{T}>100\mathrm{GeV}$. The distribution from various processes have their total count normalized to 1.

Distribution of Track $d_{xy}$ of CMS data, SM multijet MC events, and various signal samples. Events are required to pass the trigger requirements and also have 4 jets with $p_\mathrm{T}>100\mathrm{GeV}$. The distribution from various processes have their total count normalized to 1.

Distribution of track $D_{N}$ of CMS data, SM multijet MC events, and various signal samples. Events are required to pass the trigger requirements and also have 4 jets with $p_\mathrm{T}>100\mathrm{GeV}$. The distribution from various processes have their total count normalized to 1.

Distribution of track $d_{xy}$ of CMS data, SM multijet MC events, and various signal samples. Events are required to pass the trigger requirements and also have 4 jets with $p_\mathrm{T}>100\mathrm{GeV}$. The distribution from various processes have their total count normalized to 1.

Distribution of jet $\overline{\tau}_{2/1}$ of CMS data, SM multijet MC events, and various signal samples. Events are required to pass the trigger requirements and also have 4 jets with $p_\mathrm{T}>100\mathrm{GeV}$. The distribution from various processes have their total count normalized to 1.

Distribution of $|\Delta\eta|$ between jet and tracks of CMS data, SM multijet MC events, and various signal samples. Events are required to pass the trigger requirements and also have 4 jets with $p_\mathrm{T}>100\mathrm{GeV}$. The distribution from various processes have their total count normalized to 1.

Distribution of $|\Delta\phi|$ between jet and tracks of CMS data, SM multijet MC events, and various signal samples. Events are required to pass the trigger requirements and also have 4 jets with $p_\mathrm{T}>100\mathrm{GeV}$. The distribution from various processes have their total count normalized to 1.

Distribution of track $d_{xy}$ of CMS data, SM multijet MC events, and various signal samples. Events are required to pass the trigger requirements and also have 4 jets with $p_\mathrm{T}>100\mathrm{GeV}$. The distribution from various processes have their total count normalized to 1.

Distribution of track $p_\mathrm{T}$ of CMS data, SM multijet MC events, and various signal samples. Events are required to pass the trigger requirements and also have 4 jets with $p_\mathrm{T}>100\mathrm{GeV}$. The distribution from various processes have their total count normalized to 1.

Distribution of track $d_{z}$ of CMS data, SM multijet MC events, and various signal samples. Events are required to pass the trigger requirements and also have 4 jets with $p_\mathrm{T}>100\mathrm{GeV}$. The distribution from various processes have their total count normalized to 1.

Distribution of track $p_\mathrm{T}$ fraction of CMS data, SM multijet MC events, and various signal samples. Events are required to pass the trigger requirements and also have 4 jets with $p_\mathrm{T}>100\mathrm{GeV}$. The distribution from various processes have their total count normalized to 1.

The EJ tagger misidentification probability for b quark jets and light jets as a function of jet $p_\mathrm{T}$ for the taggers a-tag 1, a-tag 2, a-tag 3, and a-tag 4, as defined in the paper, evaluated using data and generator-level flavor information from simulated samples in events containing a high-$p_\mathrm{T}$ photon.

The EJ tagger misidentification probability for b quark jets and light jets as a function of jet $p_\mathrm{T}$ for the taggers aGNN tag 1, aGNN tag 2, and aGNN tag 3, as defined in the paper, evaluated using data and generator-level flavor information from simulated samples in events containing a high-$p_\mathrm{T}$ photon.

Signal acceptance of signal models in the unflavored scenario with using the designated cut sets with the model agnostic and GNN EJ tagging methods.

Signal acceptance of signal models in the flavor-aligned scenario with using the designated cut sets with the model agnostic and GNN EJ tagging methods.

Tagging efficiency and its uncertainty of the taggers used to select signal jets in the unflavored scenario, as a function of the $p_\mathrm{T}$ and $|\eta|$ of the signal jet.

Tagging efficiency and its uncertainty of the taggers used to select signal jets in the unflavored scenario, as a function of the $p_\mathrm{T}$ and $|\eta|$ of the signal jet.

Tagging efficiency and its uncertainty of the taggers used to select signal jets in the unflavored scenario, as a function of the $p_\mathrm{T}$ and $|\eta|$ of the signal jet.

Tagging efficiency and its uncertainty of the taggers used to select signal jets in the unflavored scenario, as a function of the $p_\mathrm{T}$ and $|\eta|$ of the signal jet.

Tagging efficiency and its uncertainty of the taggers used to select signal jets in the unflavored scenario, as a function of the $p_\mathrm{T}$ and $|\eta|$ of the signal jet.

Tagging efficiency and its uncertainty of the taggers used to select signal jets in the unflavored scenario, as a function of the $p_\mathrm{T}$ and $|\eta|$ of the signal jet.

Tagging efficiency and its uncertainty of the taggers used to select signal jets in the unflavored scenario, as a function of the $p_\mathrm{T}$ and $|\eta|$ of the signal jet.

Tagging efficiency and its uncertainty of the taggers used to select signal jets in the unflavored scenario, as a function of the $p_\mathrm{T}$ and $|\eta|$ of the signal jet.

Tagging efficiency and its uncertainty of the taggers used to select signal jets in the unflavored scenario, as a function of the $p_\mathrm{T}$ and $|\eta|$ of the signal jet.

Tagging efficiency and its uncertainty of the taggers used to select signal jets in the unflavored scenario, as a function of the $p_\mathrm{T}$ and $|\eta|$ of the signal jet.

Tagging efficiency and its uncertainty of the taggers used to select signal jets in the flavor-aligned scenario, as a function of the $p_\mathrm{T}$ and $|\eta|$ of the signal jet.

Tagging efficiency and its uncertainty of the taggers used to select signal jets in the flavor-aligned scenario, as a function of the $p_\mathrm{T}$ and $|\eta|$ of the signal jet.

Tagging efficiency and its uncertainty of the taggers used to select signal jets in the flavor-aligned scenario, as a function of the $p_\mathrm{T}$ and $|\eta|$ of the signal jet.

Tagging efficiency and its uncertainty of the taggers used to select signal jets in the flavor-aligned scenario, as a function of the $p_\mathrm{T}$ and $|\eta|$ of the signal jet.

Tagging efficiency and its uncertainty of the taggers used to select signal jets in the flavor-aligned scenario, as a function of the $p_\mathrm{T}$ and $|\eta|$ of the signal jet.

Tagging efficiency and its uncertainty of the taggers used to select signal jets in the flavor-aligned scenario, as a function of the $p_\mathrm{T}$ and $|\eta|$ of the signal jet.

Tagging efficiency and its uncertainty of the taggers used to select signal jets in the flavor-aligned scenario, as a function of the $p_\mathrm{T}$ and $|\eta|$ of the signal jet.

Tagging efficiency and its uncertainty of the taggers used to select signal jets in the flavor-aligned scenario, as a function of the $p_\mathrm{T}$ and $|\eta|$ of the signal jet.

Tagging efficiency and its uncertainty of the taggers used to select signal jets in the flavor-aligned scenario, as a function of the $p_\mathrm{T}$ and $|\eta|$ of the signal jet.

Distribution of jet-associated tracks of jets in SM multijet processes that pass the basic kinematic selection, the uGNN tag 1 tagging criteria, and aGNN tag 1 EMJ tagging criteria as a function of $\Delta R$ between the jet and track of interest and $\mathrm{sign}(d_{xy})\cdot\ln\left( 1+\left\lvert\frac{d_{xy}}{1\mathrm{cm}}\right\rvert\right)$, with $d_{xy}$ being the transverse impact parameter of the track of interest.

Example cut flow used for selecting signals events using the u-set 4 selection criteria. Both background and MC events has been normalized to unity.

Example cut flow used for selecting signals events using the u-set 1 selection criteria. Both background and MC events has been normalized to unity.

Example cut flow used for selecting signals events using the u-set 2 selection criteria. Both background and MC events has been normalized to unity.

Example cut flow used for selecting signals events using the u-set 3 selection criteria. Both background and MC events has been normalized to unity.

Example cut flow used for selecting signals events using the u-set 5 selection criteria. Both background and MC events has been normalized to unity.

Example cut flow used for selecting signals events using the uGNN set 3 selection criteria. Both background and MC events has been normalized to unity.

Example cut flow used for selecting signals events using the uGNN set 1 selection criteria. Both background and MC events has been normalized to unity.

Example cut flow used for selecting signals events using the uGNN set 2 selection criteria. Both background and MC events has been normalized to unity.

Example cut flow used for selecting signals events using the a-set 3 selection criteria. Both background and MC events has been normalized to unity.

Example cut flow used for selecting signals events using the a-set 1 selection criteria. Both background and MC events has been normalized to unity.

Example cut flow used for selecting signals events using the a-set 2 selection criteria. Both background and MC events has been normalized to unity.

Example cut flow used for selecting signals events using the a-set 4 selection criteria. Both background and MC events has been normalized to unity.

Example cut flow used for selecting signals events using the a-set 5 selection criteria. Both background and MC events has been normalized to unity.

Example cut flow used for selecting signals events using the aGNN set 3 selection criteria. Both background and MC events has been normalized to unity.

Example cut flow used for selecting signals events using the aGNN set 1 selection criteria. Both background and MC events has been normalized to unity.

Example cut flow used for selecting signals events using the aGNN set 2 selection criteria. Both background and MC events has been normalized to unity.