The total coherent $\mathrm{J}/\psi$ photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$, measured in PbPb ultra-peripheral collisions at $\sqrt{s_{\mathrm{NN}}}$ = 5.02 TeV. The $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$ values used correspond to the center of each rapidity range. The theoretical uncertainties is due to the uncertainties in the photon flux.

The nuclear gluon suppression factor $R_{\mathrm{g}}^{\mathrm{Pb}}$ as a function of Bjorken $x$ extracted from the CMS measurement of the coherent $\mathrm{J}/\psi$ photoproduction in PbPb ultra-peripheral collisions at $\sqrt{s_{\mathrm{NN}}}$ = 5.02 TeV. The $x$ values are evaluated at the centers of their corresponding rapidity ranges. The theoretical uncertainties are due to the uncertainties in the photon flux and the impulse approximation model.

The nuclear gluon suppression factor $R_{\mathrm{g}}^{\mathrm{Pb}}$ as a function of Bjorken $x$ extracted from the CMS measurement of the coherent $\mathrm{J}/\psi$ photoproduction in PbPb ultra-peripheral collisions at $\sqrt{s_{\mathrm{NN}}}$ = 5.02 TeV. The $x$ values are evaluated at the centers of their corresponding rapidity ranges. The theoretical uncertainties are due to the uncertainties in the photon flux and the impulse approximation model.

The total covariance matrix of the total coherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The covariance matrix includes both the experimental and theoretical (photon flux) uncertainties. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

The photon flux values from STARLight as a function of rapidity, in different neutron multiplicity classes: 0n0n, 0nXn, and XnXn.

The experimental covariance matrix of the total coherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

The covariance matrix for the flux in the 0n0n neutron multiplicity class.

The theoretical (photon flux) covariance matrix of the total coherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

The covariance matrix for the flux in the 0nXn neutron multiplicity class.

The covariance matrix for the flux in the XnXn neutron multiplicity class.

The total covariance matrix of the total coherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The covariance matrix includes both the experimental and theoretical (photon flux) uncertainties. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

The experimental covariance matrix of the total coherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

The theoretical (photon flux) covariance matrix of the total coherent photoproduction cross section as a function of photon-nuclear center-of-mass energy per nucleon $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$. The bins are ordered as increasing in $W_{\gamma \mathrm{N}}^{\mathrm{Pb}}$.

Distribution of an input to the BDT classifier in the $3\ell$ category: The linear discrimminant $p_{T}^{miss,LD}$. The discriminant is defined by the relation $p_{T}^{miss,LD} = 0.6 p_{T}^{miss} + 0.4 H_{T}^{miss}$, where $H_{T}^{miss}$ corresponds to the magnitude of the vector $p_{T}$ sum of all electrons, muons, taus and small radius jets passing the selection critertia, and $p_{T}^{miss}$ represents the missing transverse momentum vector computed as the negative vector $p_{T}$ sum of all the particles reconstructed by the PF algorithm in the event.

Distributions of the transverse mass $m_{T}$ of the lepton that is not originating from the Z boson decay and the missing transverse momentum in the $3\ell$/WZ control region. The distributions expected for the WZ and ZZ as well as for other background processes are shown for the values of nuisance parameters obtained from the ML fit used in the signal extraction.

Distributions of the four lepton invariant mass $m_{4\ell}$ in the $4\ell$/ZZ control region. The distributions expected for the ZZ as well as for other background processes are shown for the values of nuisance parameters obtained from the ML fit used in the signal extraction.

Distributions of the transverse mass $m_{T}$ of the leading lepton and the missing transverse momentum in the $2\ell$(ss) control region. The distributions expected for the misidentified $\ell$ background as well as for other background processes is shown for the values of nuisance parameters obtained from the background-only ML fit, in which the HH signal is constrained to be zero.

Distributions of the reconstructed HH mass $m_{HH}$ in the $2\ell+2 au_{h}$ control region. The distributions expected for the misidentified $\ell/ au_{h}$ background as well as for other background processes is shown for the values of nuisance parameters obtained from the background-only ML fit, in which the HH signal is constrained to be zero.

Distribution in the output of the BDT trained for nonresonant HH production and evaluated for the benchmark scenario JHEP04 BM7 for the $2\ell$(ss) category. The SM HH signal is shown for a cross section amounting to $30$ times the value predicted in the SM. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate signal quantiles.

Distribution in the output of the BDT trained for nonresonant HH production and evaluated for the benchmark scenario JHEP04 BM7 for the $3\ell$ category. The SM HH signal is shown for a cross section amounting to $30$ times the value predicted in the SM. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate signal quantiles.

Distribution in the output of the BDT trained for nonresonant HH production and evaluated for the benchmark scenario JHEP04 BM7 for the $4\ell$ category. The SM HH signal is shown for a cross section amounting to $30$ times the value predicted in the SM. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate signal quantiles.

Distribution in the output of the BDT trained for nonresonant HH production and evaluated for the benchmark scenario JHEP04 BM7 for the $3\ell+1 au_{h}$ category. The SM HH signal is shown for a cross section amounting to $30$ times the value predicted in the SM. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate background quantiles.

Distribution in the output of the BDT trained for nonresonant HH production and evaluated for the benchmark scenario JHEP04 BM7 for the $2\ell+2 au_{h}$ category. The SM HH signal is shown for a cross section amounting to $30$ times the value predicted in the SM. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate background quantiles.

Distribution in the output of the BDT trained for nonresonant HH production and evaluated for the benchmark scenario JHEP04 BM7 for the $1\ell+3 au_{h}$ category. The SM HH signal is shown for a cross section amounting to $30$ times the value predicted in the SM. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate background quantiles.

Distribution in the output of the BDT trained for nonresonant HH production and evaluated for the benchmark scenario JHEP04 BM7 for the $4 au_{h}$ category. The SM HH signal is shown for a cross section amounting to $30$ times the value predicted in the SM. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate background quantiles.

Observed and expected $95\%$ CL upper limits on the SM HH production cross section, obtained for both individual search categories and from a simultaneous fit of all seven categories combined.

Observed and expected $95\%$ CL upper limits on the HH production cross section as a function of the Higgs boson self-coupling strength modifier $\kappa_\lambda$ for the combination of all seven categories. All Higgs boson couplings other than $\lambda$ are assumed to have the values predicted in the SM.

Observed and expected $95\%$ CL upper limits on the HH production cross section as a function of the Higgs boson self-coupling strength modifier $\kappa_\lambda$ for the seven different categories as well as their combination. All Higgs boson couplings other than $\lambda$ are assumed to have the values predicted in the SM.

Observed and expected $95\%$ CL upper limits on the HH production cross section for the twelve benchmark scenarios from doi:10.1007/JHEP04(2016)126, the additional benchmark scenario 8a from doi:10.1007/JHEP09(2018)057, the seven benchmark scenarios from doi:10.1007/JHEP03(2020)091, and for the SM.

Observed and expected $95\%$ CL upper limits on the HH production cross section for the twelve benchmark scenarios from doi:10.1007/JHEP04(2016)126, the additional benchmark scenario 8a from doi:10.1007/JHEP09(2018)057, the seven benchmark scenarios from doi:10.1007/JHEP03(2020)091, and for the SM. The limits are shown for all seven search categories as well as for their combination.

Observed and expected $95\%$ CL upper limits on the HH production cross section as a function of the effective coupling $c_{2}$ for the combination of all seven categories. All Higgs boson couplings other than $c_{2}$ are assumed to have the values predicted in the SM.

Observed and expected $95\%$ CL upper limits on the HH production cross section as a function of the effective coupling $c_{2}$ and the top Yukawa coupling modifier $\kappa_{t}$ for the combination of all seven categories. All Higgs boson couplings other than $c_{2}$ and $\kappa_{t}$ are assumed to have the values predicted in the SM. The position of the SM in the ($c_{2}-\kappa_{t}$) plane, as well as the best fit value of $(c_{2},\kappa_{t})=(1.05, 1.74)$ together with contours for the theory cross section are shown as well.

Observed and expected $95\%$ CL upper limits on the HH production cross section as a function of the Higgs boson self-coupling modifier $\kappa_{\lambda}$ and the top Yukawa coupling modifier $\kappa_{t}$ for the combination of all seven categories. All Higgs boson couplings other than $\kappa_{\lambda}$ and $\kappa_{t}$ are assumed to have the values predicted in the SM. The position of the SM in the ($\kappa_{\lambda}-\kappa_{t}$) plane, as well as the best fit value of $(\kappa_{\lambda},\kappa_{t})=(-3.54, -1.70)$ together with contours for the theory cross section are shown as well.

Observed and expected $95\%$ CL upper limits on the HH production cross section as a function of the effective coupling $c_{2}$ and the Higgs boson self-coupling modifier $\kappa_{\lambda}$ for the combination of all seven categories. All Higgs boson couplings other than $c_{2}$ and $\kappa_{\lambda}$ are assumed to have the values predicted in the SM. The position of the SM in the ($c_{2}-\kappa_{\lambda}$) plane, as well as the best fit value of $(c_{2},\kappa_{\lambda})=(-0.66, 5.33)$ together with contours for the theory cross section are shown as well.

Observed and expected $95\%$ CL upper limits on the production of new particles X of spin $0$ and mass $m_{X}$ in the range $250$-$1000$ GeV, which decay to Higgs boson pairs.

Observed and expected $95\%$ CL upper limits on the production of new particles X of spin $0$ and mass $m_{X}$ in the range $250$-$1000$ GeV, which decay to Higgs boson pairs. The Limit is shown for the seven different search categories as well as their combination.

Observed and expected $95\%$ CL upper limits on the production of new particles X of spin $2$ and mass $m_{X}$ in the range $250$-$1000$ GeV, which decay to Higgs boson pairs.

Observed and expected $95\%$ CL upper limits on the production of new particles X of spin $2$ and mass $m_{X}$ in the range $250$-$1000$ GeV, which decay to Higgs boson pairs. The Limit is shown for the seven different search categories as well as their combination.

Distribution in the output of the BDT trained for resonances of spin 2 and mass $750\,$GeV production for the $2\ell$(ss) category. The resonant HH signal is shown for a cross section amounting to $1\,$pb. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate signal quantiles.

Distribution in the output of the BDT trained for resonances of spin 2 and mass $750\,$GeV production for the $3\ell$ category. The resonant HH signal is shown for a cross section amounting to $1\,$pb. The distributions expected for the background processes are shown for the values of nuisance parameters obtained from the ML fit of the signal+background hypothesis to the data. The binning is chosen to approximate signal quantiles.

The covariance matrix containing the experimental uncertainties of the particle-level $\mathrm{t}\overline{\mathrm{t}}$ differential cross section as a function of the XCone-jet mass.

The covariance matrix containing the model uncertainties of the particle-level $\mathrm{t}\overline{\mathrm{t}}$ differential cross section as a function of the XCone-jet mass.

The normalized particle-level $\mathrm{t}\overline{\mathrm{t}}$ differential cross section in the fiducial region as a function of the XCone-jet mass.

Correlations between bins in the normalized particle-level $\mathrm{t}\overline{\mathrm{t}}$ differential cross section as a function of the XCone-jet mass.

The covariance matrix containing the statistical uncertainties of the normalized particle-level $\mathrm{t}\overline{\mathrm{t}}$ differential cross section as a function of the XCone-jet mass.

The covariance matrix containing the experimental uncertainties of the normalized particle-level $\mathrm{t}\overline{\mathrm{t}}$ differential cross section as a function of the XCone-jet mass.

The covariance matrix containing the model uncertainties of the normalized particle-level $\mathrm{t}\overline{\mathrm{t}}$ differential cross section as a function of the XCone-jet mass.

Relative experimental uncertainties of the particle-level $\mathrm{t}\overline{\mathrm{t}}$ differential cross section in the fiducial region as a function of the XCone-jet mass.

Relative model uncertainties of the particle-level $\mathrm{t}\overline{\mathrm{t}}$ differential cross section in the fiducial region as a function of the XCone-jet mass.

Relative experimental uncertainties of the normalized particle-level $\mathrm{t}\overline{\mathrm{t}}$ differential cross section in the fiducial region as a function of the XCone-jet mass.

Relative model uncertainties of the normalized particle-level $\mathrm{t}\overline{\mathrm{t}}$ differential cross section in the fiducial region as a function of the XCone-jet mass.

Measured and predicted fiducial cross-sections in both the lllljj and ll$\nu\nu$jj channels for the inclusive ZZjj processes. Uncertainties due to different sources are presented.

Measured and predicted fiducial cross-sections in both the lllljj and ll$\nu\nu$jj channels for the inclusive ZZjj processes. Uncertainties due to different sources are presented.

Observed and expected multivariate discriminant distribution in the $\ell\ell\ell\ell jj$ QCD CR.

Observed and expected multivariate discriminant distribution in the $\ell\ell\ell\ell jj$ QCD CR.

Observed and expected multivariate discriminant distribution in the $\ell\ell\ell\ell jj$ SR.

Observed and expected multivariate discriminant distribution in the $\ell\ell\ell\ell jj$ SR.

Observed and expected multivariate discriminant distribution in the $\ell\ell\nu\nu jj$ SR.

Observed and expected multivariate discriminant distribution in the $\ell\ell\nu\nu jj$ SR.

Ranking of the systematic uncertainties on the measured cross-section, normalised to the predicted value, in the inclusive fit to data. The impact of each nuisance parameter, $\Delta \sigma_{\text{inc}}/\sigma^{\text{pred.}}_{\text{inc}}$, is computed by comparing the nominal best-fit value of $\sigma_{\text{inc}}/\sigma^{\text{pred}}_{\text{inc}}$ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\theta$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta$ ($\pm \Delta \hat{\theta}$). The figure shows the effect of the ten most significant uncertainties.

Ranking of the systematic uncertainties on the measured cross-section, normalised to the predicted value, in the fiducial fit to data. The impact of each nuisance parameter, $\Delta \sigma_{\text{fid}}/\sigma^{\text{pred.}}_{\text{fid}}$, is computed by comparing the nominal best-fit value of $\sigma_{\text{fid}}/\sigma^{\text{pred}}_{\text{fid}}$ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\theta$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta$ ($\pm \Delta \hat{\theta}$). The figure shows the effect of the ten most significant uncertainties.

Ranking of the systematic uncertainties on the measured cross-section, normalised to the predicted value, in the fiducial fit to data. The impact of each nuisance parameter, $\Delta \sigma_{\text{fid}}/\sigma^{\text{pred.}}_{\text{fid}}$, is computed by comparing the nominal best-fit value of $\sigma_{\text{fid}}/\sigma^{\text{pred}}_{\text{fid}}$ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\theta$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta$ ($\pm \Delta \hat{\theta}$). The figure shows the effect of the ten most significant uncertainties.

Impact of different categories of systematic uncertainties on the fiducial and inclusive measurements. The quoted values are obtained by repeating the fit, fixing a set of nuisance parameters of the sources corresponding to the considered category, and subtracting in quadrature the resulting uncertainty from the total uncertainty of the nominal fit. The total uncertainty is different from the sum in quadrature of the different components due to correlations between nuisance parameters built by the fit.

Impact of different categories of systematic uncertainties on the fiducial and inclusive measurements. The quoted values are obtained by repeating the fit, fixing a set of nuisance parameters of the sources corresponding to the considered category, and subtracting in quadrature the resulting uncertainty from the total uncertainty of the nominal fit. The total uncertainty is different from the sum in quadrature of the different components due to correlations between nuisance parameters built by the fit.

Fiducial region definition

Fiducial region definition

The Bose-Einstein correlation (BEC) parameter R as a function of n_ch for HMT events using different MC generators in the calculation of R₂(Q). The uncertainties shown are statistical. The lower panel of each plot shows the ratio of the BEC parameters obtained using EPOS LHC (red circles), Pythia 8 Monash (blue squares) and Herwig++ UE-EE-5 (green triangles) compared with the parameters obtained using Pythia 8 A2. The gray band in the lower panels is the MC systematic uncertainty, obtained as explained in the text.

The Bose-Einstein correlation (BEC) parameter R as a function of k_T for MB events using different MC generators in the calculation of R₂(Q). The uncertainties shown are statistical. The lower panel of each plot shows the ratio of the BEC parameters obtained using EPOS LHC (red circles), Pythia 8 Monash (blue squares) and Herwig++ UE-EE-5 (green triangles) compared with the parameters obtained using Pythia 8 A2. The gray band in the lower panels is the MC systematic uncertainty, obtained as explained in the text.

The Bose-Einstein correlation (BEC) parameter λ as a function of k_T for MB events using different MC generators in the calculation of R₂(Q). The uncertainties shown are statistical. The lower panel of each plot shows the ratio of the BEC parameters obtained using EPOS LHC (red circles), Pythia 8 Monash (blue squares) and Herwig++ UE-EE-5 (green triangles) compared with the parameters obtained using Pythia 8 A2. The gray band in the lower panels is the MC systematic uncertainty, obtained as explained in the text.

The two-particle double-ratio correlation function, R₂(Q), for pp collisions for track p_T >100 MeV at √s=13 TeV in the multiplicity interval 71 ≤ n_ch < 80 for the minimum-bias (MB) events. The blue dashed and red solid lines show the results of the exponential and Gaussian fits, respectively. The region excluded from the fits is shown. The statistical uncertainty and the systematic uncertainty for imperfections in the data reconstruction procedure are added in quadrature.

The two-particle double-ratio correlation function, R₂(Q), for pp collisions for track p_T >100 MeV at √s=13 TeV in the multiplicity interval 231 ≤ n_ch < 300 for the high-multiplicity track (HMT) events. The blue dashed and red solid lines show the results of the exponential and Gaussian fits, respectively. The region excluded from the fits is shown. The statistical uncertainty and the systematic uncertainty for imperfections in the data reconstruction procedure are added in quadrature.

The dependence of the correlation strength, λ(m_ch), on rescaled multiplicity, m_ch, obtained from the exponential fit of the R₂(Q) correlation functions for tracks with p_T > 100  MeV and p_T > 500 MeV at √s = 13  TeV for the minimum-bias (MB) and high multiplicity track (HMT) data. The uncertainties represent the sum in quadrature of the statistical and asymmetric systematic contributions. The black and blue solid curves represent the exponential fit of λ(m_ch) for p_T >100 MeV and p_T >500 MeV, respectively.

The dependence of the correlation strength, λ(m_ch), on rescaled multiplicity, m_ch, obtained from the exponential fit of the R₂(Q) correlation functions for tracks with p_T > 100  MeV and p_T > 500 MeV at √s = 13  TeV for the minimum-bias (MB) and high multiplicity track (HMT) data. The uncertainties represent the sum in quadrature of the statistical and asymmetric systematic contributions. The black and blue solid curves represent the exponential fit of λ(m_ch) for p_T >100 MeV and p_T >500 MeV, respectively.

The dependence of the correlation strength, λ(m_ch), on rescaled multiplicity, m_ch, obtained from the exponential fit of the R₂(Q) correlation functions for tracks with p_T > 100  MeV and p_T > 500 MeV at √s = 13  TeV for the minimum-bias (MB) and high multiplicity track (HMT) data. The uncertainties represent the sum in quadrature of the statistical and asymmetric systematic contributions. The black and blue solid curves represent the exponential fit of λ(m_ch) for p_T >100 MeV and p_T >500 MeV, respectively.

The dependence of the correlation strength, λ(m_ch), on rescaled multiplicity, m_ch, obtained from the exponential fit of the R₂(Q) correlation functions for tracks with p_T > 100  MeV and p_T > 500 MeV at √s = 13  TeV for the minimum-bias (MB) and high multiplicity track (HMT) data. The uncertainties represent the sum in quadrature of the statistical and asymmetric systematic contributions. The black and blue solid curves represent the exponential fit of λ(m_ch) for p_T >100 MeV and p_T >500 MeV, respectively.

The dependence of the source radius, R(m_ch), on m_ch. The uncertainties represent the sum in quadrature of the statistical and asymmetric systematic contributions. The black and blue solid curves represent the fit of R(m_ch) for ∛m_ch < 1.2 for p_T >100 MeV and p_T >500 MeV, respectively. The black and blue dotted curves are extensions of the black and blue solid curves beyond ∛m_ch > 1.2, respectively. The black and brown dashed curves represent the saturation value of R(m_ch) for ∛m_ch > 1.45 with p_T >100 MeV and for ∛m_ch > 1.6 with p_T >500 MeV, respectively.

The dependence of the source radius, R(m_ch), on m_ch. The uncertainties represent the sum in quadrature of the statistical and asymmetric systematic contributions. The black and blue solid curves represent the fit of R(m_ch) for ∛m_ch < 1.2 for p_T >100 MeV and p_T >500 MeV, respectively. The black and blue dotted curves are extensions of the black and blue solid curves beyond ∛m_ch > 1.2, respectively. The black and brown dashed curves represent the saturation value of R(m_ch) for ∛m_ch > 1.45 with p_T >100 MeV and for ∛m_ch > 1.6 with p_T >500 MeV, respectively.

The dependence of the source radius, R(m_ch), on m_ch. The uncertainties represent the sum in quadrature of the statistical and asymmetric systematic contributions. The black and blue solid curves represent the fit of R(m_ch) for ∛m_ch < 1.2 for p_T >100 MeV and p_T >500 MeV, respectively. The black and blue dotted curves are extensions of the black and blue solid curves beyond ∛m_ch > 1.2, respectively. The black and brown dashed curves represent the saturation value of R(m_ch) for ∛m_ch > 1.45 with p_T >100 MeV and for ∛m_ch > 1.6 with p_T >500 MeV, respectively.

The dependence of the source radius, R(m_ch), on m_ch. The uncertainties represent the sum in quadrature of the statistical and asymmetric systematic contributions. The black and blue solid curves represent the fit of R(m_ch) for ∛m_ch < 1.2 for p_T >100 MeV and p_T >500 MeV, respectively. The black and blue dotted curves are extensions of the black and blue solid curves beyond ∛m_ch > 1.2, respectively. The black and brown dashed curves represent the saturation value of R(m_ch) for ∛m_ch > 1.45 with p_T >100 MeV and for ∛m_ch > 1.6 with p_T >500 MeV, respectively.

The dependence of the R(m_ch) on ∛m_ch. The uncertainties represent the sum in quadrature of the statistical and asymmetric systematic contributions. The black and blue solid curves represent the fit of R(m_ch) for ∛m_ch < 1.2 for p_T >100 MeV and p_T >500 MeV, respectively. The black and blue dotted curves are extensions of the black and blue solid curves beyond ∛m_ch > 1.2, respectively. The black and brown dashed curves represent the saturation value of R(m_ch) for ∛m_ch > 1.45 with p_T >100 MeV and for ∛m_ch > 1.6 with p_T >500 MeV, respectively

The dependence of the R(m_ch) on ∛m_ch. The uncertainties represent the sum in quadrature of the statistical and asymmetric systematic contributions. The black and blue solid curves represent the fit of R(m_ch) for ∛m_ch < 1.2 for p_T >100 MeV and p_T >500 MeV, respectively. The black and blue dotted curves are extensions of the black and blue solid curves beyond ∛m_ch > 1.2, respectively. The black and brown dashed curves represent the saturation value of R(m_ch) for ∛m_ch > 1.45 with p_T >100 MeV and for ∛m_ch > 1.6 with p_T >500 MeV, respectively

The dependence of the R(m_ch) on ∛m_ch. The uncertainties represent the sum in quadrature of the statistical and asymmetric systematic contributions. The black and blue solid curves represent the fit of R(m_ch) for ∛m_ch < 1.2 for p_T >100 MeV and p_T >500 MeV, respectively. The black and blue dotted curves are extensions of the black and blue solid curves beyond ∛m_ch > 1.2, respectively. The black and brown dashed curves represent the saturation value of R(m_ch) for ∛m_ch > 1.45 with p_T >100 MeV and for ∛m_ch > 1.6 with p_T >500 MeV, respectively

The dependence of the R(m_ch) on ∛m_ch. The uncertainties represent the sum in quadrature of the statistical and asymmetric systematic contributions. The black and blue solid curves represent the fit of R(m_ch) for ∛m_ch < 1.2 for p_T >100 MeV and p_T >500 MeV, respectively. The black and blue dotted curves are extensions of the black and blue solid curves beyond ∛m_ch > 1.2, respectively. The black and brown dashed curves represent the saturation value of R(m_ch) for ∛m_ch > 1.45 with p_T >100 MeV and for ∛m_ch > 1.6 with p_T >500 MeV, respectively

Comparison of single-ratio two-particle correlation functions, using the unlike-charge particle (UCP) pair reference sample, for minimum-bias (MB) events, showing C₂^data(Q) (top panel) at 13 TeV (black circles) and 7 TeV (open blue circles), and the ratio of C₂^7 TeV (Q) to C₂^13 TeV (Q) (bottom panel). Comparison of C₂^data (Q) for representative multiplicity region 3.09 < m_ch ≤ 3.86. The statistical and systematic uncertainties, combined in quadrature, are presented. The systematic uncertainties include track efficiency, Coulomb correction, non-closure and multiplicity-unfolding uncertainties.

Comparison of single-ratio two-particle correlation functions, using the unlike-charge particle (UCP) pair reference sample, for minimum-bias (MB) events, showing C₂^data(Q) (top panel) at 13  TeV (black circles) and 7  TeV (open blue circles), and the ratio of C₂^7 TeV (Q) to C₂^13 TeV (Q) (bottom panel). Comparison of C₂^data (Q) for representative k_T region 400 < k_T ≤500  MeV. The statistical and systematic uncertainties, combined in quadrature, are presented. The systematic uncertainties include track efficiency, Coulomb correction, non-closure and multiplicity-unfolding uncertainties.

The k_T dependence of the correlation strength, λ(k_T), obtained from the exponential fit to the R₂(Q) correlation functions for events with multiplicity n_ch ≥ 2 and transfer momentum of tracks with p_T >100 MeV and p_T >500 MeV at √s=13 TeV for the minimum-bias (MB) and high-multiplicity track (HMT) events. The uncertainties represent the sum in quadrature of the statistical and systematic contributions. The curves represent the exponential fits to λ(k_T).

The k_T dependence of the correlation strength, λ(k_T), obtained from the exponential fit to the R₂(Q) correlation functions for events with multiplicity n_ch ≥ 2 and transfer momentum of tracks with p_T >100 MeV and p_T >500 MeV at √s=13 TeV for the minimum-bias (MB) and high-multiplicity track (HMT) events. The uncertainties represent the sum in quadrature of the statistical and systematic contributions. The curves represent the exponential fits to λ(k_T).

The k_T dependence of the correlation strength, λ(k_T), obtained from the exponential fit to the R₂(Q) correlation functions for events with multiplicity n_ch ≥ 2 and transfer momentum of tracks with p_T >100 MeV and p_T >500 MeV at √s=13 TeV for the minimum-bias (MB) and high-multiplicity track (HMT) events. The uncertainties represent the sum in quadrature of the statistical and systematic contributions. The curves represent the exponential fits to λ(k_T).

The k_T dependence of the correlation strength, λ(k_T), obtained from the exponential fit to the R₂(Q) correlation functions for events with multiplicity n_ch ≥ 2 and transfer momentum of tracks with p_T >100 MeV and p_T >500 MeV at √s=13 TeV for the minimum-bias (MB) and high-multiplicity track (HMT) events. The uncertainties represent the sum in quadrature of the statistical and systematic contributions. The curves represent the exponential fits to λ(k_T).

The k_T dependence of the source radius, R(k_T), obtained from the exponential fit to the R₂(Q) correlation functions for events with multiplicity n_ch ≥ 2 and transfer momentum of tracks with p_T >100 MeV and p_T >500 MeV at √s=13 TeV for the minimum-bias (MB) and high-multiplicity track (HMT) events. The uncertainties represent the sum in quadrature of the statistical and systematic contributions. The curves represent the exponential fits to R(k_T).

The k_T dependence of the source radius, R(k_T), obtained from the exponential fit to the R₂(Q) correlation functions for events with multiplicity n_ch ≥ 2 and transfer momentum of tracks with p_T >100 MeV and p_T >500 MeV at √s=13 TeV for the minimum-bias (MB) and high-multiplicity track (HMT) events. The uncertainties represent the sum in quadrature of the statistical and systematic contributions. The curves represent the exponential fits to R(k_T).

The k_T dependence of the source radius, R(k_T), obtained from the exponential fit to the R₂(Q) correlation functions for events with multiplicity n_ch ≥ 2 and transfer momentum of tracks with p_T >100 MeV and p_T >500 MeV at √s=13 TeV for the minimum-bias (MB) and high-multiplicity track (HMT) events. The uncertainties represent the sum in quadrature of the statistical and systematic contributions. The curves represent the exponential fits to R(k_T).

The k_T dependence of the source radius, R(k_T), obtained from the exponential fit to the R₂(Q) correlation functions for events with multiplicity n_ch ≥ 2 and transfer momentum of tracks with p_T >100 MeV and p_T >500 MeV at √s=13 TeV for the minimum-bias (MB) and high-multiplicity track (HMT) events. The uncertainties represent the sum in quadrature of the statistical and systematic contributions. The curves represent the exponential fits to R(k_T).

The two-dimensional dependence on m_ch and k_T for p_T > 100 MeV for the correlation strength, λ, obtained from the exponential fit to the R₂(Q) correlation functions using the MB sample for m_ch ≤ 3.08 and the HMT sample for m_ch > 3.08.

The two-dimensional dependence on m_ch and k_T for p_T > 100 MeV for the source radius, R, obtained from the exponential fit to the R₂(Q) correlation functions using the MB sample for m_ch ≤ 3.08 and the HMT sample for m_ch > 3.08.

The parameter λ for p_T > 100 MeV as a function of k_T in selected low m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter λ for p_T > 100 MeV as a function of k_T in selected low m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter λ for p_T > 100 MeV as a function of k_T in selected high m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter λ for p_T > 100 MeV as a function of k_T in selected high m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter λ for p_T > 100 MeV as a function of m_ch in k_T intervals between 0.1 and 0.5 GeV. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter λ for p_T > 100 MeV as a function of m_ch in k_T intervals between 0.1 and 0.5 GeV. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter λ for p_T > 100 MeV as a function of m_ch in k_T intervals between 0.5 and 1.5 GeV. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter λ for p_T > 100 MeV as a function of m_ch in k_T intervals between 0.5 and 1.5 GeV. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter R for p_T > 100 MeV as a function of k_T in selected low m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter R for p_T > 100 MeV as a function of k_T in selected low m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter R for p_T > 100 MeV as a function of k_T in selected high m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter R for p_T > 100 MeV as a function of k_T in selected high m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter R for p_T > 100 MeV as a function of m_ch in k_T intervals between 0.1 and 0.5 GeV. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter R for p_T > 100 MeV as a function of m_ch in k_T intervals between 0.1 and 0.5 GeV. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter R for p_T > 100 MeV as a function of m_ch in k_T intervals between 0.5 and 1.5 GeV. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter R for p_T > 100 MeV as a function of m_ch in k_T intervals between 0.5 and 1.5 GeV. The error bars and boxes represent the statistical and systematic contributions, respectively.

The fit parameter μ describing the dependence of the correlation strength, λ, on charged-particle scaled multiplicity, for track p_T>100 MeV and track p_T>500 MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at √s = 13 TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid (blue dashed) curve represents the exponential fit of the dependence of parameter μ on m_ch for tracks with p_T >100 MeV (p_T >500 MeV).

The fit parameter μ describing the dependence of the correlation strength, λ, on charged-particle scaled multiplicity, for track p_T>100 MeV and track p_T>500 MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at √s = 13 TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid (blue dashed) curve represents the exponential fit of the dependence of parameter μ on m_ch for tracks with p_T >100 MeV (p_T >500 MeV).

The fit parameter μ describing the dependence of the correlation strength, λ, on charged-particle scaled multiplicity, for track p_T>100 MeV and track p_T>500 MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at √s = 13 TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid (blue dashed) curve represents the exponential fit of the dependence of parameter μ on m_ch for tracks with p_T >100 MeV (p_T >500 MeV).

The fit parameter μ describing the dependence of the correlation strength, λ, on charged-particle scaled multiplicity, for track p_T>100 MeV and track p_T>500 MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at √s = 13 TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid (blue dashed) curve represents the exponential fit of the dependence of parameter μ on m_ch for tracks with p_T >100 MeV (p_T >500 MeV).

The fit parameter ν describing the dependence of the correlation strength, λ, on charged-particle scaled multiplicity, for track p_T>100 MeV and track p_T>500 MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at √s = 13 TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid (blue dashed) curve represents the exponential fit of the dependence of parameter ν on m_ch for tracks with p_T >100 MeV (p_T >500 MeV).

The fit parameter ν describing the dependence of the correlation strength, λ, on charged-particle scaled multiplicity, for track p_T>100 MeV and track p_T>500 MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at √s = 13 TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid (blue dashed) curve represents the exponential fit of the dependence of parameter ν on m_ch for tracks with p_T >100 MeV (p_T >500 MeV).

The fit parameter ν describing the dependence of the correlation strength, λ, on charged-particle scaled multiplicity, for track p_T>100 MeV and track p_T>500 MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at √s = 13 TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid (blue dashed) curve represents the exponential fit of the dependence of parameter ν on m_ch for tracks with p_T >100 MeV (p_T >500 MeV).

The fit parameter ν describing the dependence of the correlation strength, λ, on charged-particle scaled multiplicity, for track p_T>100 MeV and track p_T>500 MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at √s = 13 TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid (blue dashed) curve represents the exponential fit of the dependence of parameter ν on m_ch for tracks with p_T >100 MeV (p_T >500 MeV).

The parameter ξ describing the dependence of the source radius, R, on charged-particle scaled multiplicity, m_ch, for track p_T>100 MeV and track p_T>500 MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at √s = 13 TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid and blue dashed curves represent the saturated value of the parameter ξ for m_ch > 3.0 for tracks with p_T >100 MeV and for m_ch > 2.8 for tracks with p_T >500 MeV, respectively.

The parameter ξ describing the dependence of the source radius, R, on charged-particle scaled multiplicity, m_ch, for track p_T>100 MeV and track p_T>500 MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at √s = 13 TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid and blue dashed curves represent the saturated value of the parameter ξ for m_ch > 3.0 for tracks with p_T >100 MeV and for m_ch > 2.8 for tracks with p_T >500 MeV, respectively.

The parameter ξ describing the dependence of the source radius, R, on charged-particle scaled multiplicity, m_ch, for track p_T>100 MeV and track p_T>500 MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at √s = 13 TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid and blue dashed curves represent the saturated value of the parameter ξ for m_ch > 3.0 for tracks with p_T >100 MeV and for m_ch > 2.8 for tracks with p_T >500 MeV, respectively.

The parameter ξ describing the dependence of the source radius, R, on charged-particle scaled multiplicity, m_ch, for track p_T>100 MeV and track p_T>500 MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at √s = 13 TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid and blue dashed curves represent the saturated value of the parameter ξ for m_ch > 3.0 for tracks with p_T >100 MeV and for m_ch > 2.8 for tracks with p_T >500 MeV, respectively.

The parameter κ describing the dependence of the source radius, R, on charged-particle scaled multiplicity, m_ch, for track p_T>100 MeV and track p_T>500 MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at √s = 13 TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid and blue dashed curves represent the exponential fit to the parameter κ for tracks with p_T >100 MeV and for tracks with p_T >500 MeV, respectively.

The parameter κ describing the dependence of the source radius, R, on charged-particle scaled multiplicity, m_ch, for track p_T>100 MeV and track p_T>500 MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at √s = 13 TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid and blue dashed curves represent the exponential fit to the parameter κ for tracks with p_T >100 MeV and for tracks with p_T >500 MeV, respectively.

The parameter κ describing the dependence of the source radius, R, on charged-particle scaled multiplicity, m_ch, for track p_T>100 MeV and track p_T>500 MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at √s = 13 TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid and blue dashed curves represent the exponential fit to the parameter κ for tracks with p_T >100 MeV and for tracks with p_T >500 MeV, respectively.

The parameter κ describing the dependence of the source radius, R, on charged-particle scaled multiplicity, m_ch, for track p_T>100 MeV and track p_T>500 MeV in the minimum-bias (MB) and high-multiplicity track (HMT) samples at √s = 13 TeV. The error bars and boxes represent the statistical and systematic contributions, respectively. The black solid and blue dashed curves represent the exponential fit to the parameter κ for tracks with p_T >100 MeV and for tracks with p_T >500 MeV, respectively.

The two-dimensional dependence on m_ch and k_T for p_T > 500 MeV for the correlation strength, λ, obtained from the exponential fit to the R₂(Q) correlation functions using the MB sample for m_ch ≤ 3.08 and the HMT sample for m_ch > 3.08.

The two-dimensional dependence on m_ch and k_T for p_T > 500 MeV for the source radius, R, obtained from the exponential fit to the R₂(Q) correlation functions using the MB sample for m_ch ≤ 3.08 and the HMT sample for m_ch > 3.08.

The parameter λ for p_T > 500 MeV as a function of k_T in selected low m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter λ for p_T > 500 MeV as a function of k_T in selected low m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter λ for p_T > 500 MeV as a function of k_T in selected high m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter λ for p_T > 500 MeV as a function of k_T in selected high m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter λ for p_T > 500 MeV as a function of m_ch in k_T intervals between 0.5 and 1.5 GeV. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter λ for p_T > 500 MeV as a function of m_ch in k_T intervals between 0.5 and 1.5 GeV. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter R for p_T > 500 MeV as a function of k_T in selected low m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter R for p_T > 500 MeV as a function of k_T in selected low m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter R for p_T > 500 MeV as a function of k_T in selected high m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter R for p_T > 500 MeV as a function of k_T in selected high m_ch intervals. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter R for p_T > 500 MeV as a function of m_ch in k_T intervals between 0.5 and 1.5 GeV. The error bars and boxes represent the statistical and systematic contributions, respectively.

The parameter R for p_T > 500 MeV as a function of m_ch in k_T intervals between 0.5 and 1.5 GeV. The error bars and boxes represent the statistical and systematic contributions, respectively.

ATLAS and CMS results for the source radius R as a function of n_ch in pp interactions at 13 TeV. The CMS results (open circles) have been adjusted (by the CMS collaboration) to the ATLAS kinematic region∶ p_T > 100 MeV and |η|<2.5. The ATLAS uncertainties are the sum in quadrature of the statistical and asymmetric systematic uncertainties. For CMS, only the systematic uncertainties are shown since the statistical uncertainties are smaller than the marker size. The dashed blue (ATLAS) and black (CMS) lines represent the fit to ∛n_ch at low multiplicity, continued as dotted lines beyond the fit range. The solid green (ATLAS) and broken black (CMS) lines indicate the plateau level at high multiplicity.

ATLAS and CMS results for the source radius R as a function of n_ch in pp interactions at 13 TeV. The CMS results (open circles) have been adjusted (by the CMS collaboration) to the ATLAS kinematic region∶ p_T > 100 MeV and |η|<2.5. The ATLAS uncertainties are the sum in quadrature of the statistical and asymmetric systematic uncertainties. For CMS, only the systematic uncertainties are shown since the statistical uncertainties are smaller than the marker size. The dashed blue (ATLAS) and black (CMS) lines represent the fit to ∛n_ch at low multiplicity, continued as dotted lines beyond the fit range. The solid green (ATLAS) and broken black (CMS) lines indicate the plateau level at high multiplicity.

ATLAS and CMS results for the source radius R as a function of n_ch in pp interactions at 13 TeV. The CMS results (open circles) have been adjusted (by the CMS collaboration) to the ATLAS kinematic region∶ p_T > 100 MeV and |η|<2.5. The ATLAS uncertainties are the sum in quadrature of the statistical and asymmetric systematic uncertainties. For CMS, only the systematic uncertainties are shown since the statistical uncertainties are smaller than the marker size. The dashed blue (ATLAS) and black (CMS) lines represent the fit to ∛n_ch at low multiplicity, continued as dotted lines beyond the fit range. The solid green (ATLAS) and broken black (CMS) lines indicate the plateau level at high multiplicity.

ATLAS and CMS results for the source radius R as a function of ∛n_ch in pp interactions at 13 TeV. The CMS results (open circles) have been adjusted (by the CMS collaboration) to the ATLAS kinematic region∶ p_T > 100 MeV and |η|<2.5. The ATLAS uncertainties are the sum in quadrature of the statistical and asymmetric systematic uncertainties. For CMS, only the systematic uncertainties are shown since the statistical uncertainties are smaller than the marker size. The dashed blue (ATLAS) and black (CMS) lines represent the fit to ∛n_ch at low multiplicity, continued as dotted lines beyond the fit range. The solid green (ATLAS) and broken black (CMS) lines indicate the plateau level at high multiplicity.

ATLAS and CMS results for the source radius R as a function of ∛n_ch in pp interactions at 13 TeV. The CMS results (open circles) have been adjusted (by the CMS collaboration) to the ATLAS kinematic region∶ p_T > 100 MeV and |η|<2.5. The ATLAS uncertainties are the sum in quadrature of the statistical and asymmetric systematic uncertainties. For CMS, only the systematic uncertainties are shown since the statistical uncertainties are smaller than the marker size. The dashed blue (ATLAS) and black (CMS) lines represent the fit to ∛n_ch at low multiplicity, continued as dotted lines beyond the fit range. The solid green (ATLAS) and broken black (CMS) lines indicate the plateau level at high multiplicity.

ATLAS and CMS results for the source radius R as a function of ∛n_ch in pp interactions at 13 TeV. The CMS results (open circles) have been adjusted (by the CMS collaboration) to the ATLAS kinematic region∶ p_T > 100 MeV and |η|<2.5. The ATLAS uncertainties are the sum in quadrature of the statistical and asymmetric systematic uncertainties. For CMS, only the systematic uncertainties are shown since the statistical uncertainties are smaller than the marker size. The dashed blue (ATLAS) and black (CMS) lines represent the fit to ∛n_ch at low multiplicity, continued as dotted lines beyond the fit range. The solid green (ATLAS) and broken black (CMS) lines indicate the plateau level at high multiplicity.

Systematic uncertainties (in percent) in the correlation strength, λ, and source radius, R, for the exponential fit of the two-particle double-ratio correlation functions, R₂(Q), for p_T > 100 MeV at √s= 13 TeV for the MB and HMT events. The choice of MC generator gives rise to asymmetric uncertainties, denoted by uparrow and downarrow. This asymmetry propagates through to the cumulative uncertainty. The columns under ‘Uncertainty range’ show the range of systematic uncertainty from the fits in the various n_ch intervals.

The results of the fits to the dependencies of the correlation strength, λ, and source radius, R, on the average rescaled charged-particle multiplicity, m_ch, for |η| < 2.5 and both p_T > 100 MeV and p_T > 500 MeV at √s = 13 TeV for the minimum-bias (MB) and the high-multiplicity track (HMT) events. The parameters γ and δ resulting from a joint fit to the MB and HMT data are presented. The total uncertainties are shown.

The results of the fits to the dependencies of the correlation strength, λ, and source radius, R, on the pair average transverse momentum, k_T, for various functional forms and for minimum-bias (MB) and high-multiplicity track (HMT) events for p_T > 100 MeV and p_T > 500 MeV at √s = 13 TeV. The total uncertainties are shown.

The Bose-Einstein correlation (BEC) parameters λ and R as a function of n_ch and k_T using different MC generators in the calculation of R₂(Q). (a) λ versus n_ch for MB events, (b) λ versus n_ch for HMT events, (c) λ versus k_T and (d) R versus k_T for MB events. The uncertainties shown are statistical. The lower panel of each plot shows the ratio of the BEC parameters obtained using EPOS LHC (red circles), Pythia 8 Monash (blue squares) and Herwig++ UE-EE-5 (green triangles) compared with the parameters obtained using Pythia 8 A2. The gray band in the lower panels is the MC systematic uncertainty, obtained as explained in the text.

The Bose-Einstein correlation (BEC) parameters λ and R as a function of n_ch and k_T using different MC generators in the calculation of R₂(Q). (a) λ versus n_ch for MB events, (b) λ versus n_ch for HMT events, (c) λ versus k_T and (d) R versus k_T for MB events. The uncertainties shown are statistical. The lower panel of each plot shows the ratio of the BEC parameters obtained using EPOS LHC (red circles), Pythia 8 Monash (blue squares) and Herwig++ UE-EE-5 (green triangles) compared with the parameters obtained using Pythia 8 A2. The gray band in the lower panels is the MC systematic uncertainty, obtained as explained in the text.

The Bose-Einstein correlation (BEC) parameters λ and R as a function of n_ch and k_T using different MC generators in the calculation of R₂(Q). (a) λ versus n_ch for MB events, (b) λ versus n_ch for HMT events, (c) λ versus k_T and (d) R versus k_T for MB events. The uncertainties shown are statistical. The lower panel of each plot shows the ratio of the BEC parameters obtained using EPOS LHC (red circles), Pythia 8 Monash (blue squares) and Herwig++ UE-EE-5 (green triangles) compared with the parameters obtained using Pythia 8 A2. The gray band in the lower panels is the MC systematic uncertainty, obtained as explained in the text.

The Bose-Einstein correlation (BEC) parameters λ and R as a function of n_ch and k_T using different MC generators in the calculation of R₂(Q). (a) λ versus n_ch for MB events, (b) λ versus n_ch for HMT events, (c) λ versus k_T and (d) R versus k_T for MB events. The uncertainties shown are statistical. The lower panel of each plot shows the ratio of the BEC parameters obtained using EPOS LHC (red circles), Pythia 8 Monash (blue squares) and Herwig++ UE-EE-5 (green triangles) compared with the parameters obtained using Pythia 8 A2. The gray band in the lower panels is the MC systematic uncertainty, obtained as explained in the text.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals∶ (a) 2 < n_ch ≤ 10, (b) 11 < n_ch ≤ 20, (c) 21 < n_ch ≤ 30, (d) 31 < n_ch ≤ 40, (e) 41 < n_ch ≤ 50, (f) 51 < n_ch ≤ 60, (g) 61 < n_ch ≤ 70, (h) 71 < n_ch ≤ 80 and (i) 81 < n_ch ≤ 90. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals∶ (a) 2 < n_ch ≤ 10, (b) 11 < n_ch ≤ 20, (c) 21 < n_ch ≤ 30, (d) 31 < n_ch ≤ 40, (e) 41 < n_ch ≤ 50, (f) 51 < n_ch ≤ 60, (g) 61 < n_ch ≤ 70, (h) 71 < n_ch ≤ 80 and (i) 81 < n_ch ≤ 90. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals∶ (a) 2 < n_ch ≤ 10, (b) 11 < n_ch ≤ 20, (c) 21 < n_ch ≤ 30, (d) 31 < n_ch ≤ 40, (e) 41 < n_ch ≤ 50, (f) 51 < n_ch ≤ 60, (g) 61 < n_ch ≤ 70, (h) 71 < n_ch ≤ 80 and (i) 81 < n_ch ≤ 90. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals∶ (a) 2 < n_ch ≤ 10, (b) 11 < n_ch ≤ 20, (c) 21 < n_ch ≤ 30, (d) 31 < n_ch ≤ 40, (e) 41 < n_ch ≤ 50, (f) 51 < n_ch ≤ 60, (g) 61 < n_ch ≤ 70, (h) 71 < n_ch ≤ 80 and (i) 81 < n_ch ≤ 90. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals∶ (a) 2 < n_ch ≤ 10, (b) 11 < n_ch ≤ 20, (c) 21 < n_ch ≤ 30, (d) 31 < n_ch ≤ 40, (e) 41 < n_ch ≤ 50, (f) 51 < n_ch ≤ 60, (g) 61 < n_ch ≤ 70, (h) 71 < n_ch ≤ 80 and (i) 81 < n_ch ≤ 90. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals∶ (a) 2 < n_ch ≤ 10, (b) 11 < n_ch ≤ 20, (c) 21 < n_ch ≤ 30, (d) 31 < n_ch ≤ 40, (e) 41 < n_ch ≤ 50, (f) 51 < n_ch ≤ 60, (g) 61 < n_ch ≤ 70, (h) 71 < n_ch ≤ 80 and (i) 81 < n_ch ≤ 90. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals∶ (a) 2 < n_ch ≤ 10, (b) 11 < n_ch ≤ 20, (c) 21 < n_ch ≤ 30, (d) 31 < n_ch ≤ 40, (e) 41 < n_ch ≤ 50, (f) 51 < n_ch ≤ 60, (g) 61 < n_ch ≤ 70, (h) 71 < n_ch ≤ 80 and (i) 81 < n_ch ≤ 90. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals∶ (a) 2 < n_ch ≤ 10, (b) 11 < n_ch ≤ 20, (c) 21 < n_ch ≤ 30, (d) 31 < n_ch ≤ 40, (e) 41 < n_ch ≤ 50, (f) 51 < n_ch ≤ 60, (g) 61 < n_ch ≤ 70, (h) 71 < n_ch ≤ 80 and (i) 81 < n_ch ≤ 90. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals∶ (a) 2 < n_ch ≤ 10, (b) 11 < n_ch ≤ 20, (c) 21 < n_ch ≤ 30, (d) 31 < n_ch ≤ 40, (e) 41 < n_ch ≤ 50, (f) 51 < n_ch ≤ 60, (g) 61 < n_ch ≤ 70, (h) 71 < n_ch ≤ 80 and (i) 81 < n_ch ≤ 90. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals∶ (a) 91 < n_ch ≤ 100, (b) 101 < n_ch ≤ 125, (c) 126 < n_ch ≤ 150, (d) 151 < n_ch ≤ 200, (e) 201 < n_ch ≤ 250. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals∶ (a) 91 < n_ch ≤ 100, (b) 101 < n_ch ≤ 125, (c) 126 < n_ch ≤ 150, (d) 151 < n_ch ≤ 200, (e) 201 < n_ch ≤ 250. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals∶ (a) 91 < n_ch ≤ 100, (b) 101 < n_ch ≤ 125, (c) 126 < n_ch ≤ 150, (d) 151 < n_ch ≤ 200, (e) 201 < n_ch ≤ 250. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals∶ (a) 91 < n_ch ≤ 100, (b) 101 < n_ch ≤ 125, (c) 126 < n_ch ≤ 150, (d) 151 < n_ch ≤ 200, (e) 201 < n_ch ≤ 250. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals∶ (a) 91 < n_ch ≤ 100, (b) 101 < n_ch ≤ 125, (c) 126 < n_ch ≤ 150, (d) 151 < n_ch ≤ 200, (e) 201 < n_ch ≤ 250. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals∶ (a) 101 < n_ch ≤ 110, (b) 111 < n_ch ≤ 120, (c) 121 < n_ch ≤ 130, (d) 131 < n_ch ≤ 140, (e) 141 < n_ch ≤ 155, (f) 156 < n_ch ≤ 175, (g) 176 < n_ch ≤ 200, (h) 201 < n_ch ≤ 230 and (i) 231 < n_ch ≤ 300. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals∶ (a) 101 < n_ch ≤ 110, (b) 111 < n_ch ≤ 120, (c) 121 < n_ch ≤ 130, (d) 131 < n_ch ≤ 140, (e) 141 < n_ch ≤ 155, (f) 156 < n_ch ≤ 175, (g) 176 < n_ch ≤ 200, (h) 201 < n_ch ≤ 230 and (i) 231 < n_ch ≤ 300. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals∶ (a) 101 < n_ch ≤ 110, (b) 111 < n_ch ≤ 120, (c) 121 < n_ch ≤ 130, (d) 131 < n_ch ≤ 140, (e) 141 < n_ch ≤ 155, (f) 156 < n_ch ≤ 175, (g) 176 < n_ch ≤ 200, (h) 201 < n_ch ≤ 230 and (i) 231 < n_ch ≤ 300. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals∶ (a) 101 < n_ch ≤ 110, (b) 111 < n_ch ≤ 120, (c) 121 < n_ch ≤ 130, (d) 131 < n_ch ≤ 140, (e) 141 < n_ch ≤ 155, (f) 156 < n_ch ≤ 175, (g) 176 < n_ch ≤ 200, (h) 201 < n_ch ≤ 230 and (i) 231 < n_ch ≤ 300. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals∶ (a) 101 < n_ch ≤ 110, (b) 111 < n_ch ≤ 120, (c) 121 < n_ch ≤ 130, (d) 131 < n_ch ≤ 140, (e) 141 < n_ch ≤ 155, (f) 156 < n_ch ≤ 175, (g) 176 < n_ch ≤ 200, (h) 201 < n_ch ≤ 230 and (i) 231 < n_ch ≤ 300. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals∶ (a) 101 < n_ch ≤ 110, (b) 111 < n_ch ≤ 120, (c) 121 < n_ch ≤ 130, (d) 131 < n_ch ≤ 140, (e) 141 < n_ch ≤ 155, (f) 156 < n_ch ≤ 175, (g) 176 < n_ch ≤ 200, (h) 201 < n_ch ≤ 230 and (i) 231 < n_ch ≤ 300. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals∶ (a) 101 < n_ch ≤ 110, (b) 111 < n_ch ≤ 120, (c) 121 < n_ch ≤ 130, (d) 131 < n_ch ≤ 140, (e) 141 < n_ch ≤ 155, (f) 156 < n_ch ≤ 175, (g) 176 < n_ch ≤ 200, (h) 201 < n_ch ≤ 230 and (i) 231 < n_ch ≤ 300. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals∶ (a) 101 < n_ch ≤ 110, (b) 111 < n_ch ≤ 120, (c) 121 < n_ch ≤ 130, (d) 131 < n_ch ≤ 140, (e) 141 < n_ch ≤ 155, (f) 156 < n_ch ≤ 175, (g) 176 < n_ch ≤ 200, (h) 201 < n_ch ≤ 230 and (i) 231 < n_ch ≤ 300. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals∶ (a) 101 < n_ch ≤ 110, (b) 111 < n_ch ≤ 120, (c) 121 < n_ch ≤ 130, (d) 131 < n_ch ≤ 140, (e) 141 < n_ch ≤ 155, (f) 156 < n_ch ≤ 175, (g) 176 < n_ch ≤ 200, (h) 201 < n_ch ≤ 230 and (i) 231 < n_ch ≤ 300. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for k_T - intervals∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for k_T - intervals∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for k_T - intervals∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for k_T - intervals∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for k_T - intervals∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for k_T - intervals∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for k_T - intervals∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), for the high-multiplicity track (HMT) events using the unlike-charge particle (UCP) pairs reference sample for k_T - intervals∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals∶ (a) 2 < n_ch ≤ 9, (b) 10 < n_ch ≤ 18, (c) 19 < n_ch ≤ 27, (d) 28 < n_ch ≤ 36, (e) 37 < n_ch ≤ 45, (f) 46 < n_ch ≤ 54, (g) 55 < n_ch ≤ 63, (h) 64 < n_ch ≤ 72, (i) 73 < n_ch ≤ 81, (j) 82 < n_ch ≤ 90, (k) 91 < n_ch ≤ 113, and (l) 114 < n_ch ≤ 136. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals∶ (a) 2 < n_ch ≤ 9, (b) 10 < n_ch ≤ 18, (c) 19 < n_ch ≤ 27, (d) 28 < n_ch ≤ 36, (e) 37 < n_ch ≤ 45, (f) 46 < n_ch ≤ 54, (g) 55 < n_ch ≤ 63, (h) 64 < n_ch ≤ 72, (i) 73 < n_ch ≤ 81, (j) 82 < n_ch ≤ 90, (k) 91 < n_ch ≤ 113, and (l) 114 < n_ch ≤ 136. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals∶ (a) 2 < n_ch ≤ 9, (b) 10 < n_ch ≤ 18, (c) 19 < n_ch ≤ 27, (d) 28 < n_ch ≤ 36, (e) 37 < n_ch ≤ 45, (f) 46 < n_ch ≤ 54, (g) 55 < n_ch ≤ 63, (h) 64 < n_ch ≤ 72, (i) 73 < n_ch ≤ 81, (j) 82 < n_ch ≤ 90, (k) 91 < n_ch ≤ 113, and (l) 114 < n_ch ≤ 136. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals∶ (a) 2 < n_ch ≤ 9, (b) 10 < n_ch ≤ 18, (c) 19 < n_ch ≤ 27, (d) 28 < n_ch ≤ 36, (e) 37 < n_ch ≤ 45, (f) 46 < n_ch ≤ 54, (g) 55 < n_ch ≤ 63, (h) 64 < n_ch ≤ 72, (i) 73 < n_ch ≤ 81, (j) 82 < n_ch ≤ 90, (k) 91 < n_ch ≤ 113, and (l) 114 < n_ch ≤ 136. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals∶ (a) 2 < n_ch ≤ 9, (b) 10 < n_ch ≤ 18, (c) 19 < n_ch ≤ 27, (d) 28 < n_ch ≤ 36, (e) 37 < n_ch ≤ 45, (f) 46 < n_ch ≤ 54, (g) 55 < n_ch ≤ 63, (h) 64 < n_ch ≤ 72, (i) 73 < n_ch ≤ 81, (j) 82 < n_ch ≤ 90, (k) 91 < n_ch ≤ 113, and (l) 114 < n_ch ≤ 136. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals∶ (a) 2 < n_ch ≤ 9, (b) 10 < n_ch ≤ 18, (c) 19 < n_ch ≤ 27, (d) 28 < n_ch ≤ 36, (e) 37 < n_ch ≤ 45, (f) 46 < n_ch ≤ 54, (g) 55 < n_ch ≤ 63, (h) 64 < n_ch ≤ 72, (i) 73 < n_ch ≤ 81, (j) 82 < n_ch ≤ 90, (k) 91 < n_ch ≤ 113, and (l) 114 < n_ch ≤ 136. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals∶ (a) 2 < n_ch ≤ 9, (b) 10 < n_ch ≤ 18, (c) 19 < n_ch ≤ 27, (d) 28 < n_ch ≤ 36, (e) 37 < n_ch ≤ 45, (f) 46 < n_ch ≤ 54, (g) 55 < n_ch ≤ 63, (h) 64 < n_ch ≤ 72, (i) 73 < n_ch ≤ 81, (j) 82 < n_ch ≤ 90, (k) 91 < n_ch ≤ 113, and (l) 114 < n_ch ≤ 136. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals∶ (a) 2 < n_ch ≤ 9, (b) 10 < n_ch ≤ 18, (c) 19 < n_ch ≤ 27, (d) 28 < n_ch ≤ 36, (e) 37 < n_ch ≤ 45, (f) 46 < n_ch ≤ 54, (g) 55 < n_ch ≤ 63, (h) 64 < n_ch ≤ 72, (i) 73 < n_ch ≤ 81, (j) 82 < n_ch ≤ 90, (k) 91 < n_ch ≤ 113, and (l) 114 < n_ch ≤ 136. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals∶ (a) 2 < n_ch ≤ 9, (b) 10 < n_ch ≤ 18, (c) 19 < n_ch ≤ 27, (d) 28 < n_ch ≤ 36, (e) 37 < n_ch ≤ 45, (f) 46 < n_ch ≤ 54, (g) 55 < n_ch ≤ 63, (h) 64 < n_ch ≤ 72, (i) 73 < n_ch ≤ 81, (j) 82 < n_ch ≤ 90, (k) 91 < n_ch ≤ 113, and (l) 114 < n_ch ≤ 136. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals∶ (a) 2 < n_ch ≤ 9, (b) 10 < n_ch ≤ 18, (c) 19 < n_ch ≤ 27, (d) 28 < n_ch ≤ 36, (e) 37 < n_ch ≤ 45, (f) 46 < n_ch ≤ 54, (g) 55 < n_ch ≤ 63, (h) 64 < n_ch ≤ 72, (i) 73 < n_ch ≤ 81, (j) 82 < n_ch ≤ 90, (k) 91 < n_ch ≤ 113, and (l) 114 < n_ch ≤ 136. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals∶ (a) 2 < n_ch ≤ 9, (b) 10 < n_ch ≤ 18, (c) 19 < n_ch ≤ 27, (d) 28 < n_ch ≤ 36, (e) 37 < n_ch ≤ 45, (f) 46 < n_ch ≤ 54, (g) 55 < n_ch ≤ 63, (h) 64 < n_ch ≤ 72, (i) 73 < n_ch ≤ 81, (j) 82 < n_ch ≤ 90, (k) 91 < n_ch ≤ 113, and (l) 114 < n_ch ≤ 136. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample for n_ch - intervals∶ (a) 2 < n_ch ≤ 9, (b) 10 < n_ch ≤ 18, (c) 19 < n_ch ≤ 27, (d) 28 < n_ch ≤ 36, (e) 37 < n_ch ≤ 45, (f) 46 < n_ch ≤ 54, (g) 55 < n_ch ≤ 63, (h) 64 < n_ch ≤ 72, (i) 73 < n_ch ≤ 81, (j) 82 < n_ch ≤ 90, (k) 91 < n_ch ≤ 113, and (l) 114 < n_ch ≤ 136. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The single-ratio two-particle correlation functions, C₂^data(Q), at 7 TeV for the minimum-bias (MB) events using the unlike-charge particle (UCP) pairs reference sample k_T - intervals∶ (a) 100 < k_T ≤ 200 MeV, (b) 200 < k_T ≤ 300 MeV, (c) 300 < k_T ≤ 400 MeV, (d) 400 < k_T ≤ 500 MeV, (e) 500 < k_T ≤ 600 MeV, (f) 600 < k_T ≤ 700 MeV, (g) 700 < k_T ≤ 1000 MeV, and (h) 1000 < k_T ≤ 1500 MeV. The error bars represent the statistical uncertainties. The boxes represent the systematic uncertainties, which are the sum in quadrature of a variation of the Coulomb correction, the track reconstruction efficiency and the unfolding matrix.

The correlation strength, λ, and source radius, R, of the exponential fits to the two-particle double-ratio correlation functions, R₂(Q), in dependence on the multiplicity, m_ch, intervals for the minimum-bias (MB) and the high-multiplicity track (HMT) events for p_T > 100 MeV at √s = 13 TeV. Statistical uncertainties for √χ²/ndf>1 are corrected by the √χ²/ndf. The total uncertainties are shown.

The correlation strength, λ, and source radius, R, of the exponential fits to the two-particle double-ratio correlation functions, R₂(Q), in dependence on the multiplicity, m_ch, intervals for the minimum-bias (MB) and the high-multiplicity track (HMT) events for p_T > 500 MeV at √s = 13 TeV. Statistical uncertainties for √χ²/ndf>1 are corrected by the √χ²/ndf. The total uncertainties are shown.

The correlation strength, λ, and source radius, R, of the exponential fits to the two-particle double-ratio correlation functions, R₂(Q), in dependence on the pair transverse momentum, k_T, intervals for the minimum-bias (MB) and the high-multiplicity track (HMT) events for p_T > 100 MeV at √s = 13 TeV. Statistical uncertainties for √χ²/ndf>1 are corrected by the √χ²/ndf. The total uncertainties are shown.

The correlation strength, λ, and source radius, R, of the exponential fits to the two-particle double-ratio correlation functions, R₂(Q), in dependence on the pair transverse momentum, k_T, intervals for the minimum-bias (MB) and the high-multiplicity track (HMT) events for p_T > 500 MeV at √s = 13 TeV. Statistical uncertainties for √χ²/ndf>1 are corrected by the √χ²/ndf. The total uncertainties are shown.

Pre-fit distribution of the reconstructed Higgs boson candidate $p_T^H$ for the dilepton $SR^{\geq 4j}_{\geq 4b}$ signal region. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations, except for the uncertainty in the $k({t\bar {t}+{\geq }1b})$ normalisation factor which is not defined pre-fit. The last bin includes the overflow.

Pre-fit distribution of the reconstructed Higgs boson candidate $p_T^H$ for the single-lepton resolved $SR^{\geq 6j}_{\geq 4b}$ signal region. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations, except for the uncertainty in the $k({t\bar {t}+{\geq }1b})$ normalisation factor which is not defined pre-fit. The last bin includes the overflow.

Pre-fit distribution of the reconstructed Higgs boson candidate $p_T^H$ for the single-lepton boosted ${{SR}_{{boosted}}}$ signal region. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations, except for the uncertainty in the $k({t\bar {t}+{\geq }1b})$ normalisation factor which is not defined pre-fit. The last bin includes the overflow.

Comparison of predicted and observed event yields in each of the control and signal regions in the dilepton channel after the fit to the data. The uncertainty band includes all uncertainties and their correlations.

Comparison of predicted and observed event yields in each of the control and signal regions in the single-lepton channels after the fit to the data. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the dilepton SRs after the inclusive fit to the data for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the dilepton SRs after the inclusive fit to the data for $120\le p_T^H<200$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the dilepton SRs after the inclusive fit to the data for $200\le p_T^H<300$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the dilepton SRs after the inclusive fit to the data for $p_{{T}}^{H}\ge 300$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton resolved SRs after the inclusive fit to the data for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton resolved SRs after the inclusive fit to the data for $120\le p_T^H<200$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton resolved SRs after the inclusive fit to the data for $200\le p_T^H<300$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton resolved SRs after the inclusive fit to the data for $300\le p_T^H<450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton resolved SRs after the inclusive fit to the data for $p_{{T}}^{H}\ge 450$ GeV (yield only). The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton boosted SRs after the inclusive fit to the data for $300\le p_T^H<450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the BDT discriminant in the single-lepton boosted SRs after the inclusive fit to the data for $p_{{T}}^{H}\ge 450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for ${\Delta R^{{avg}}_{bb}}$ after the inclusive fit to the data in the single-lepton $CR^{5j}_{{\geq}4b\ lo}$ control region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Comparison between data and prediction for ${\Delta R^{{avg}}_{bb}}$ after the inclusive fit to the data in the single-lepton $CR^{5j}_{{\geq}4b\ hi}$ control region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Post-fit yields of signal ($S$) and total background ($B$) as a function of $\log (S/B)$, compared with data. Final-discriminant bins in all dilepton and single-lepton analysis regions are combined into bins of $\log (S/B)$, with the signal normalised to the SM prediction used for the computation of $\log (S/B)$. The signal is then shown normalised to the best-fit value and the SM prediction. The lower frame reports the ratio of data to background, and this is compared with the expected ${t\bar {t}H}$-signal-plus-background yield divided by the background-only yield for the best-fit signal strength (solid red line) and the SM prediction (dashed orange line).

Comparison between data and prediction for the reconstruction BDT score for the Higgs boson candidate identified using Higgs boson information, after the inclusive fit to the data in the dilepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the average $\Delta \eta $ between $b$-tagged jets, after the inclusive fit to the data in the dilepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the likelihood discriminant, after the inclusive fit to the data in the single-lepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the average $\Delta R$ for all possible combinations of $b$-tagged jet pairs, after the inclusive fit to the data in the single-lepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the DNN $P(H)$ output for the Higgs boson candidate after the inclusive fit to the data in the single-lepton boosted channel for $300\le p_T^H<450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the DNN $P(H)$ output for the Higgs boson candidate after the inclusive fit to the data in the single-lepton boosted channel for $p_{{T}}^{H}\ge 450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Post-fit distribution of the reconstructed Higgs boson candidate mass for the dilepton $SR^{\geq 4j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Post-fit distribution of the reconstructed Higgs boson candidate mass for the single-lepton resolved $SR^{\geq 6j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Post-fit distribution of the reconstructed Higgs boson candidate mass for the single-lepton boosted ${{SR}_{{boosted}}}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Fitted values of the ${t\bar {t}H}$ signal strength parameter in the individual channels and in the inclusive signal-strength measurement.

Ranking of the 20 nuisance parameters with the largest post-fit impact on $\mu $ in the fit. Nuisance parameters corresponding to statistical uncertainties in the simulated event samples are not included. The empty blue rectangles correspond to the pre-fit impact on $\mu $ and the filled blue ones to the post-fit impact on $\mu $, both referring to the upper scale. The impact of each nuisance parameter, $\Delta \mu $, is computed by comparing the nominal best-fit value of $\mu $ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\hat{\theta }$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta $ ($\pm \Delta \hat{\theta }$). The black points show the pulls of the nuisance parameters relative to their nominal values, $\theta _0$. These pulls and their relative post-fit errors, $\Delta \hat{\theta }/\Delta \theta $, refer to the lower scale. The `ljets' (`dilep') label refers to the single-lepton (dilepton) channel.

Pre-fit distribution of the number of jets in the dilepton $SR^{\geq 4j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the Standard Model expectation. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations, except the uncertainty in the $k({t\bar {t}+{\geq }1b})$ normalisation factor that is not defined pre-fit.

Pre-fit distribution of the number of jets in the single-lepton resolved $SR^{\geq 6j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the Standard Model expectation. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations, except the uncertainty in the $k({t\bar {t}+{\geq }1b})$ normalisation factor that is not defined pre-fit.

Pre-fit distribution of the number of jets in the single-lepton boosted ${{SR}_{{boosted}}}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the Standard Model expectation. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations, except the uncertainty in the $k({t\bar {t}+{\geq }1b})$ normalisation factor that is not defined pre-fit.

Post-fit distribution of the number of jets in the dilepton $SR^{\geq 4j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Post-fit distribution of the number of jets in the single-lepton resolved $SR^{\geq 6j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Post-fit distribution of the number of jets in the single-lepton boosted ${{SR}_{{boosted}}}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Post-fit distribution of the reconstructed Higgs boson candidate $p_T^H$ for the dilepton $SR^{\geq 4j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The last bin includes the overflow.

Post-fit distribution of the reconstructed Higgs boson candidate $p_T^H$ for the single-lepton resolved $SR^{\geq 6j}_{\geq 4b}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The last bin includes the overflow.

Post-fit distribution of the reconstructed Higgs boson candidate $p_T^H$ for the single-lepton boosted ${{SR}_{{boosted}}}$ signal region. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The last bin includes the overflow.

Signal-strength measurements in the individual STXS ${\hat{p}_{{T}}^{H}}$ bins, as well as the inclusive signal strength.

95% CL simplified template cross-section upper limits in the individual STXS ${\hat{p}_{{T}}^{H}}$ bins, as well as the inclusive limit. The observed limits are shown (solid black lines), together with the expected limits both in the background-only hypothesis (dotted black lines) and in the SM hypothesis (dotted red lines). In the case of the expected limits in the background-only hypothesis, one- and two-standard-deviation uncertainty bands are also shown. The hatched uncertainty bands correspond to the theory uncertainty in the fiducial cross-section prediction in each bin.

The ratios $S/B$ (black solid line, referring to the vertical axis on the left) and $S/\sqrt{B}$ (red dashed line, referring to the vertical axis on the right) for each category in the inclusive analysis in the dilepton channel (left) and in the single-lepton channels (right), where $S$ ($B$) is the number of selected signal (background) events predicted by the simulation and normalised to a luminosity of 139 fb$^{-1}$ .

Comparison between data and prediction for the $\Delta R$ between the Higgs candidate and the ${t\bar {t}}$ candidate system, after the inclusive fit to the data in the dilepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the number of $b$-tagged jet pairs with an invariant mass within 30 GeV of 125 GeV, after the inclusive fit to the data in the dilepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the reconstruction BDT score for the Higgs boson candidate identified using Higgs boson information, after the inclusive fit to the data in the single-lepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the $\Delta R$ between the two highest ${p_{{T}}}$ $b$-tagged jets, after the inclusive fit to the data in the single-lepton resolved channel for $0\le p_T^H<120$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations.

Comparison between data and prediction for the sum of $b$-tagging discriminants of jets from Higgs, hadronic top and leptonic top candidates, after the inclusive fit to the data in the single-lepton boosted channel for $300\le p_T^H<450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Comparison between data and prediction for the sum of $b$-tagging discriminants of jets from Higgs, hadronic top and leptonic top candidates, after the inclusive fit to the data in the single-lepton boosted channel for $p_{{T}}^{H}\ge 450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Comparison between data and prediction for the hadronic top candidate invariant mass, after the inclusive fit to the data in the single-lepton boosted channel for $300\le p_T^H<450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Comparison between data and prediction for the hadronic top candidate invariant mass, after the inclusive fit to the data in the single-lepton boosted channel for $p_{{T}}^{H}\ge 450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Comparison between data and prediction for the fraction of the sum of $b$-tagging discriminants of all jets not associated to the Higgs or hadronic top candidates, after the inclusive fit to the data in the single-lepton boosted channel for $300\le p_T^H<450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Comparison between data and prediction for the fraction of the sum of $b$-tagging discriminants of all jets not associated to the Higgs or hadronic top candidates, after the inclusive fit to the data in the single-lepton boosted channel for $p_{{T}}^{H}\ge 450$ GeV. The ${t\bar {t}H}$ signal yield (solid red) is normalised to the fitted $\mu $ value from the inclusive fit. The dashed line shows the ${t\bar {t}H}$ signal distribution normalised to the total background prediction. The uncertainty band includes all uncertainties and their correlations. The first (last) bin includes the underflow (overflow).

Ranking of the 20 nuisance parameters with the largest post-fit impact on $\mu $ in the STXS fit for $0\le {\hat{p}_{{T}}^{H}}<120$ GeV. Nuisance parameters corresponding to statistical uncertainties in the simulated event samples are not included. The empty blue rectangles correspond to the pre-fit impact on $\mu $ and the filled blue ones to the post-fit impact on $\mu $, both referring to the upper scale. The impact of each nuisance parameter, $\Delta \mu $, is computed by comparing the nominal best-fit value of $\mu $ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\hat{\theta }$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta $ ($\pm \Delta \hat{\theta }$). The black points show the pulls of the nuisance parameters relative to their nominal values, $\theta _0$. These pulls and their relative post-fit errors, $\Delta \hat{\theta }/\Delta \theta $, refer to the lower scale. For experimental uncertainties that are decomposed into several independent sources, NP X corresponds to the X$^{th}$ nuisance parameter, ordered by their impact on $\mu $. The `ljets' (`dilep') label refers to the single-lepton (dilepton) channel.

Ranking of the 20 nuisance parameters with the largest post-fit impact on $\mu $ in the STXS fit for $120\le {\hat{p}_{{T}}^{H}}<200$ GeV. Nuisance parameters corresponding to statistical uncertainties in the simulated event samples are not included. The empty blue rectangles correspond to the pre-fit impact on $\mu $ and the filled blue ones to the post-fit impact on $\mu $, both referring to the upper scale. The impact of each nuisance parameter, $\Delta \mu $, is computed by comparing the nominal best-fit value of $\mu $ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\hat{\theta }$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta $ ($\pm \Delta \hat{\theta }$). The black points show the pulls of the nuisance parameters relative to their nominal values, $\theta _0$. These pulls and their relative post-fit errors, $\Delta \hat{\theta }/\Delta \theta $, refer to the lower scale. For experimental uncertainties that are decomposed into several independent sources, NP X corresponds to the X$^{th}$ nuisance parameter, ordered by their impact on $\mu $. The `ljets' (`dilep') label refers to the single-lepton (dilepton) channel.

Ranking of the 20 nuisance parameters with the largest post-fit impact on $\mu $ in the STXS fit for $200\le {\hat{p}_{{T}}^{H}}<300$ GeV. Nuisance parameters corresponding to statistical uncertainties in the simulated event samples are not included. The empty blue rectangles correspond to the pre-fit impact on $\mu $ and the filled blue ones to the post-fit impact on $\mu $, both referring to the upper scale. The impact of each nuisance parameter, $\Delta \mu $, is computed by comparing the nominal best-fit value of $\mu $ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\hat{\theta }$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta $ ($\pm \Delta \hat{\theta }$). The black points show the pulls of the nuisance parameters relative to their nominal values, $\theta _0$. These pulls and their relative post-fit errors, $\Delta \hat{\theta }/\Delta \theta $, refer to the lower scale. For experimental uncertainties that are decomposed into several independent sources, NP X corresponds to the X$^{th}$ nuisance parameter, ordered by their impact on $\mu $. The `ljets' (`dilep') label refers to the single-lepton (dilepton) channel.

Ranking of the 20 nuisance parameters with the largest post-fit impact on $\mu $ in the STXS fit for $300\le {\hat{p}_{{T}}^{H}}<450$ GeV. Nuisance parameters corresponding to statistical uncertainties in the simulated event samples are not included. The empty blue rectangles correspond to the pre-fit impact on $\mu $ and the filled blue ones to the post-fit impact on $\mu $, both referring to the upper scale. The impact of each nuisance parameter, $\Delta \mu $, is computed by comparing the nominal best-fit value of $\mu $ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\hat{\theta }$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta $ ($\pm \Delta \hat{\theta }$). The black points show the pulls of the nuisance parameters relative to their nominal values, $\theta _0$. These pulls and their relative post-fit errors, $\Delta \hat{\theta }/\Delta \theta $, refer to the lower scale. For experimental uncertainties that are decomposed into several independent sources, NP X corresponds to the X$^{th}$ nuisance parameter, ordered by their impact on $\mu $. The `ljets' (`dilep') label refers to the single-lepton (dilepton) channel.

Ranking of the 20 nuisance parameters with the largest post-fit impact on $\mu $ in the STXS fit for ${\hat{p}_{{T}}^{H}}\ge 450$ GeV. Nuisance parameters corresponding to statistical uncertainties in the simulated event samples are not included. The empty blue rectangles correspond to the pre-fit impact on $\mu $ and the filled blue ones to the post-fit impact on $\mu $, both referring to the upper scale. The impact of each nuisance parameter, $\Delta \mu $, is computed by comparing the nominal best-fit value of $\mu $ with the result of the fit when fixing the considered nuisance parameter to its best-fit value, $\hat{\theta }$, shifted by its pre-fit (post-fit) uncertainties $\pm \Delta \theta $ ($\pm \Delta \hat{\theta }$). The black points show the pulls of the nuisance parameters relative to their nominal values, $\theta _0$. These pulls and their relative post-fit errors, $\Delta \hat{\theta }/\Delta \theta $, refer to the lower scale. For experimental uncertainties that are decomposed into several independent sources, NP X corresponds to the X$^{th}$ nuisance parameter, ordered by their impact on $\mu $. The `ljets' (`dilep') label refers to the single-lepton (dilepton) channel.

95% confidence level upper limits on signal-strength measurements in the individual STXS ${\hat{p}_{{T}}^{H}}$ bins, as well as the inclusive signal-strength limit, after the fit used to extract multiple signal-strength parameters. The observed limits are shown (solid black lines), together with the expected limits both in the background-only hypothesis (dotted black lines) and in the SM hypothesis (dotted red lines). In the case of the expected limits in the background-only hypothesis, one- and two-standard-deviation uncertainty bands are also shown.

Post-fit correlation matrix (in percentages) between the $\mu $ values obtained in the STXS bins.

Performance of the Higgs boson reconstruction algorithms. For each row of `truth' ${\hat{p}_{{T}}^{H}}$, the matrix shows (in percentages) the fraction of Higgs boson candidates which are truth-matched to ${b\bar {b}}$ decays, with reconstructed $p_T^H$ in the various bins of the dilepton (left), single lepton resolved (middle) and boosted (right) channels.

Pre-fit event yields in the dilepton signal regions and control regions. All uncertainties are included except the $k({t\bar {t}+{\geq }1b})$ uncertainty that is not defined pre-fit. For the ${t\bar {t}H}$ signal, the pre-fit yield values correspond to the theoretical prediction and corresponding uncertainties. `Other sources' refers to s-channel, t-channel, $tW$, $tWZ$, $tZq$, $Z+$ jets and diboson events.

Post-fit event yields in the dilepton signal regions and control regions, after the inclusive fit in all channels. All uncertainties are included, taking into account correlations. For the ${t\bar {t}H}$ signal, the post-fit yield and uncertainties correspond to those in the inclusive signal-strength measurement. `Other sources' refers to s-channel, t-channel, $tW$, $tWZ$, $tZq$, $Z+$ jets and diboson events.

Pre-fit event yields in the single-lepton resolved and boosted signal regions and control regions. All uncertainties are included except the $k({t\bar {t}+{\geq }1b})$ uncertainty that is not defined pre-fit. For the ${t\bar {t}H}$ signal, the pre-fit yield values correspond to the theoretical prediction and corresponding uncertainties. `Other top sources' refers to s-channel, t-channel, $tWZ$ and $tZq$ events.

Post-fit event yields in the single-lepton resolved and boosted signal regions and control regions, after the inclusive fit in all channels. All uncertainties are included, taking into account correlations. For the ${t\bar {t}H}$ signal, the post-fit yield and uncertainties correspond to those in the inclusive signal-strength measurement. `Other top sources' refers to s-channel, t-channel, $tWZ$ and $tZq$ events.

Breakdown of the contributions to the uncertainties in $\mu$. The contributions from the different sources of uncertainty are evaluated after the fit. The $\Delta \mu $ values are obtained by repeating the fit after having fixed a certain set of nuisance parameters corresponding to a group of systematic uncertainties, and then evaluating $(\Delta \mu)^2$ by subtracting the resulting squared uncertainty of $\mu $ from its squared uncertainty found in the full fit. The same procedure is followed when quoting the effect of the ${t\bar {t}+{\geq }1b}$ normalisation. The total uncertainty is different from the sum in quadrature of the different components due to correlations between nuisance parameters existing in the fit.

Fraction (in percentages) of signal events, after SR and CR selections, originating from $b\bar {b}$, $WW$ and other remaining Higgs boson decay modes in the dilepton channel.

Fraction (in percentages) of signal events, after SR and CR selections, originating from $b\bar {b}$, $WW$ and other remaining Higgs boson decay modes in the single-lepton channels.

Predicted SM ${t\bar {t}H}$ cross-section in each of the five STXS ${\hat{p}_{{T}}^{H}}$ bins and signal acceptance times efficiency (including all event selection criteria) in each STXS bin as well as for the inclusive ${\hat{p}_{{T}}^{H}}$ range.

Number of expected signal events before the fit, after each selection requirement applied to enter the dilepton channel $SR^{\geq 4j}_{\geq 4b}$ region. All ${t\bar {t}H}$ signal events are included, regardless of the $H$ or ${t\bar {t}H}$ decay mode. All object corrections are applied, except for the initial number of events which is calculated using the NLO QCD+EW theoretical prediction. All quoted numbers are rounded to unity. More details on the selection criteria can be found in the text.

Number of expected signal events before the fit, after each selection requirement applied to enter the single-lepton channel resolved $SR^{\geq 6j}_{\geq 4b}$ region. All ${t\bar {t}H}$ signal events are included, regardless of the $H$ or ${t\bar {t}H}$ decay mode. All object corrections are applied, except for the initial number of events which is calculated using the NLO QCD+EW theoretical prediction. All quoted numbers are rounded to unity. More details on the selection criteria can be found in the text.

Number of expected signal events before the fit, after each selection requirement applied to enter the single-lepton channel boosted $SR_{boosted}$ region. All ${t\bar {t}H}$ signal events are included, regardless of the $H$ or ${t\bar {t}H}$ decay mode. All object corrections are applied, except for the initial number of events which is calculated using the NLO QCD+EW theoretical prediction. All quoted numbers are rounded to unity. More details on the selection criteria can be found in the text.

The measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window, in bins of the dilepton transverse momentum. The values are normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $106 \le M_{ll} < 170$ GeV mass window, in bins of the dilepton transverse momentum. The values are normalized by the bin width.

The measured differential cross section in the $106 \le M_{ll} < 170$ GeV mass window, in bins of the dilepton transverse momentum. The values are normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $170 \le M_{ll} < 350$ GeV mass window, in bins of the dilepton transverse momentum. The values are normalized by the bin width.

The measured differential cross section in the $170 \le M_{ll} < 350$ GeV mass window, in bins of the dilepton transverse momentum. The values are normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $350 \le M_{ll} < 1000$ GeV mass window, in bins of the dilepton transverse momentum. The values are normalized by the bin width.

The measured differential cross section in the $350 \le M_{ll} < 1000$ GeV mass window, in bins of the dilepton transverse momentum. The values are normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $50 \le M_{ll} < 76$ GeV mass window, in bins of the dilepton transverse momentum. At least one jet is required. The values are normalized by the bin width.

The measured differential cross section in the $50 \le M_{ll} < 76$ GeV mass window, in bins of the dilepton transverse momentum. At least one jet is required. The values are normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window, in bins of the dilepton transverse momentum. At least one jet is required. The values are normalized by the bin width.

The measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window, in bins of the dilepton transverse momentum. At least one jet is required. The values are normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $106 \le M_{ll} < 170$ GeV mass window, in bins of the dilepton transverse momentum. At least one jet is required. The values are normalized by the bin width.

The measured differential cross section in the $106 \le M_{ll} < 170$ GeV mass window, in bins of the dilepton transverse momentum. At least one jet is required. The values are normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $170 \le M_{ll} < 350$ GeV mass window, in bins of the dilepton transverse momentum. At least one jet is required. The values are normalized by the bin width.

The measured differential cross section in the $170 \le M_{ll} < 350$ GeV mass window, in bins of the dilepton transverse momentum. At least one jet is required. The values are normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $50 \le M_{ll} < 76$ GeV mass window, in bins of the $\varphi^*$ variable. The values are normalized by the bin width.

The measured differential cross section in the $50 \le M_{ll} < 76$ GeV mass window, in bins of the $\varphi^*$ variable. The values are normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window, in bins of the $\varphi^*$ variable. The values are normalized by the bin width.

The measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window, in bins of the $\varphi^*$ variable. The values are normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $106 \le M_{ll} < 170$ GeV mass window, in bins of the $\varphi^*$ variable. The values are normalized by the bin width.

The measured differential cross section in the $106 \le M_{ll} < 170$ GeV mass window, in bins of the $\varphi^*$ variable. The values are normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $170 \le M_{ll} < 350$ GeV mass window, in bins of the $\varphi^*$ variable. The values are normalized by the bin width.

The measured differential cross section in the $170 \le M_{ll} < 350$ GeV mass window, in bins of the $\varphi^*$ variable. The values are normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $350 \le M_{ll} < 1000$ GeV mass window, in bins of the $\varphi^*$ variable. The values are normalized by the bin width.

The measured differential cross section in the $350 \le M_{ll} < 1000$ GeV mass window, in bins of the $\varphi^*$ variable. The values are normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $50 \le M_{ll} < 76$ GeV mass window, in bins of the dilepton transverse momentum, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. The values are not normalized by the bin width.

The measured differential cross section in the $50 \le M_{ll} < 76$ GeV mass window, in bins of the dilepton transverse momentum, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. The values are not normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $106 \le M_{ll} < 170$ GeV mass window, in bins of the dilepton transverse momentum, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. The values are not normalized by the bin width.

The measured differential cross section in the $106 \le M_{ll} < 170$ GeV mass window, in bins of the dilepton transverse momentum, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. The values are not normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $170 \le M_{ll} < 350$ GeV mass window, in bins of the dilepton transverse momentum, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. The values are not normalized by the bin width.

The measured differential cross section in the $170 \le M_{ll} < 350$ GeV mass window, in bins of the dilepton transverse momentum, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. The values are not normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $350 \le M_{ll} < 1000$ GeV mass window, in bins of the dilepton transverse momentum, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. The values are not normalized by the bin width.

The measured differential cross section in the $350 \le M_{ll} < 1000$ GeV mass window, in bins of the dilepton transverse momentum, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. The values are not normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $50 \le M_{ll} < 76$ GeV mass window, in bins of the dilepton transverse momentum, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. At least one jet is required. The values are not normalized by the bin width.

The measured differential cross section in the $50 \le M_{ll} < 76$ GeV mass window, in bins of the dilepton transverse momentum, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. At least one jet is required. The values are not normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $106 \le M_{ll} < 170$ GeV mass window, in bins of the dilepton transverse momentum, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. At least one jet is required. The values are not normalized by the bin width.

The measured differential cross section in the $106 \le M_{ll} < 170$ GeV mass window, in bins of the dilepton transverse momentum, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. At least one jet is required. The values are not normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $170 \le M_{ll} < 350$ GeV mass window, in bins of the dilepton transverse momentum, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. At least one jet is required. The values are not normalized by the bin width.

The measured differential cross section in the $170 \le M_{ll} < 350$ GeV mass window, in bins of the dilepton transverse momentum, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. At least one jet is required. The values are not normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $50 \le M_{ll} < 76$ GeV mass window, in bins of the $\varphi^*$ variable, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. At least one jet is required. The values are not normalized by the bin width.

The measured differential cross section in the $50 \le M_{ll} < 76$ GeV mass window, in bins of the $\varphi^*$ variable, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. At least one jet is required. The values are not normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $106 \le M_{ll} < 170$ GeV mass window, in bins of the $\varphi^*$ variable, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. At least one jet is required. The values are not normalized by the bin width.

The measured differential cross section in the $106 \le M_{ll} < 170$ GeV mass window, in bins of the $\varphi^*$ variable, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. At least one jet is required. The values are not normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $170 \le M_{ll} < 350$ GeV mass window, in bins of the $\varphi^*$ variable, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. At least one jet is required. The values are not normalized by the bin width.

The measured differential cross section in the $170 \le M_{ll} < 350$ GeV mass window, in bins of the $\varphi^*$ variable, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. At least one jet is required. The values are not normalized by the bin width. This entry contains the covariance matrix of the results.

The measured differential cross section in the $350 \le M_{ll} < 1000$ GeV mass window, in bins of the $\varphi^*$ variable, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. At least one jet is required. The values are not normalized by the bin width.

The measured differential cross section in the $350 \le M_{ll} < 1000$ GeV mass window, in bins of the $\varphi^*$ variable, divided by the measured differential cross section in the $76 \le M_{ll} < 106$ GeV mass window. At least one jet is required. The values are not normalized by the bin width. This entry contains the covariance matrix of the results.

Response matrix for pT ll mass 50-76 (electron channel)

Response matrix for pT ll mass 50-76 (muon channel)

Response matrix for pT ll mass 76-106 (electron channel)

Response matrix for pT ll mass 76-106 (muon channel)

Response matrix for pT ll mass 106-170 (electron channel)

Response matrix for pT ll mass 106-170 (muon channel)

Response matrix for pT ll mass 170-350 (electron channel)

Response matrix for pT ll mass 170-350 (muon channel)

Response matrix for pT ll mass 350-1000 (electron channel)

Response matrix for pT ll mass 50-76 jet (electron channel)

Response matrix for pT ll mass 50-76 jet (muon channel)

Response matrix for pT ll mass 76-106 jet (electron channel)

Response matrix for pT ll mass 76-106 jet (muon channel)

Response matrix for pT ll mass 106-170 jet (electron channel)

Response matrix for pT ll mass 106-170 jet (muon channel)

Response matrix for pT ll mass 170-350 jet (electron channel)

Response matrix for pT ll mass 170-350 jet (muon channel)

Response matrix for phistar mass 50-76 (electron channel)

Response matrix for phistar mass 50-76 (muon channel)

Response matrix for phistar mass 76-106 (electron channel)

Response matrix for phistar mass 76-106 (muon channel)

Response matrix for phistar mass 106-170 (electron channel)

Response matrix for phistar mass 106-170 (muon channel)

Response matrix for phistar mass 170-350 (electron channel)

Response matrix for phistar mass 170-350 (muon channel)

Response matrix for phistar mass 350-1000 (electron channel)

Response matrix for phistar mass 350-1000 (muon channel)