The figure shows the post-fit distribution of events as a function of $H_{\mathrm{T}}^{\ell,j} = \sum_{\ell,j} p_{T}^{\ell,j}$, scalar sum of $p_T$ for all jets and leptons in the $\ell+$jets channel, in proton-lead (p+Pb) collisions at a center-of-mass energy of $\sqrt{s_{\mathrm{NN}}} = 8.16$ TeV, with an integrated luminosity of 165 nb$^{-1}$. The data correspond to the $1\ell 1b$ $\mu$+jets channel in a pre-fit configuration. The stacked histograms represent different processes contributing to the event yield, including top quark pair production ($t\bar{t}$), single top, $W$ boson production with $b$, $c$, and light quarks, $Z$ boson production with $b$, $c$, and light quarks, diboson, and fake lepton backgrounds.

The figure shows the pre-fit distribution of events as a function of $H_{\mathrm{T}}^{\ell,j} = \sum_{\ell,j} p_{T}^{\ell,j}$, scalar sum of $p_T$ for all jets and leptons in the $\ell+$jets channel, in proton-lead (p+Pb) collisions at a center-of-mass energy of $\sqrt{s_{\mathrm{NN}}} = 8.16$ TeV, with an integrated luminosity of 165 nb$^{-1}$. The data correspond to the $1\ell 2bincl$ $e$+jets channel in a pre-fit configuration. The stacked histograms represent different processes contributing to the event yield, including top quark pair production ($t\bar{t}$), single top, $W$ boson production with $b$, $c$, and light quarks, $Z$ boson production with $b$, $c$, and light quarks, diboson, and fake lepton backgrounds.

The figure shows the post-fit distribution of events as a function of $H_{\mathrm{T}}^{\ell,j} = \sum_{\ell,j} p_{T}^{\ell,j}$, scalar sum of $p_T$ for all jets and leptons in the $\ell+$jets channel, in proton-lead (p+Pb) collisions at a center-of-mass energy of $\sqrt{s_{\mathrm{NN}}} = 8.16$ TeV, with an integrated luminosity of 165 nb$^{-1}$. The data correspond to the $1\ell 2bincl$ $e$+jets channel in a pre-fit configuration. The stacked histograms represent different processes contributing to the event yield, including top quark pair production ($t\bar{t}$), single top, $W$ boson production with $b$, $c$, and light quarks, $Z$ boson production with $b$, $c$, and light quarks, diboson, and fake lepton backgrounds.

The figure shows the pre-fit distribution of events as a function of $H_{\mathrm{T}}^{\ell,j} = \sum_{\ell,j} p_{T}^{\ell,j}$, scalar sum of $p_T$ for all jets and leptons in the $\ell+$jets channel, in proton-lead (p+Pb) collisions at a center-of-mass energy of $\sqrt{s_{\mathrm{NN}}} = 8.16$ TeV, with an integrated luminosity of 165 nb$^{-1}$. The data correspond to the $1\ell 2bincl$ $\mu$+jets channel in a pre-fit configuration. The stacked histograms represent different processes contributing to the event yield, including top quark pair production ($t\bar{t}$), single top, $W$ boson production with $b$, $c$, and light quarks, $Z$ boson production with $b$, $c$, and light quarks, diboson, and fake lepton backgrounds.

The figure shows the post-fit distribution of events as a function of $H_{\mathrm{T}}^{\ell,j} = \sum_{\ell,j} p_{T}^{\ell,j}$, scalar sum of $p_T$ for all jets and leptons in the $\ell+$jets channel, in proton-lead (p+Pb) collisions at a center-of-mass energy of $\sqrt{s_{\mathrm{NN}}} = 8.16$ TeV, with an integrated luminosity of 165 nb$^{-1}$. The data correspond to the $1\ell 2bincl$ $\mu$+jets channel in a pre-fit configuration. The stacked histograms represent different processes contributing to the event yield, including top quark pair production ($t\bar{t}$), single top, $W$ boson production with $b$, $c$, and light quarks, $Z$ boson production with $b$, $c$, and light quarks, diboson, and fake lepton backgrounds.

The figure shows the pre-fit distribution of events as a function of $H_{\mathrm{T}}^{\ell,j} = \sum_{\ell,j} p_{T}^{\ell,j}$, scalar sum of $p_T$ for all jets and leptons in the $\ell+$jets channel, in proton-lead (p+Pb) collisions at a center-of-mass energy of $\sqrt{s_{\mathrm{NN}}} = 8.16$ TeV, with an integrated luminosity of 165 nb\(^{-1}\) in the 2\(\ell\)1b signal region. The histogram represents the number of events per GeV, with the data points shown as black dots and their associated uncertainties. The stacked colored histograms indicate different processes: \(t\bar{t}\) (red), single top (orange), \(Z+b\) (dark blue), \(Z+c\) (medium blue), \(Z+\)light (light blue), diboson (gray), and fake lepton (purple). The hatched area represents the uncertainty in the predictions.

The figure shows the post-fit distribution of events as a function of $H_{\mathrm{T}}^{\ell,j} = \sum_{\ell,j} p_{T}^{\ell,j}$, scalar sum of $p_T$ for all jets and leptons in the $\ell+$jets channel, in proton-lead (p+Pb) collisions at a center-of-mass energy of $\sqrt{s_{\mathrm{NN}}} = 8.16$ TeV, with an integrated luminosity of 165 nb\(^{-1}\) in the 2\(\ell\)1b signal region. The histogram represents the number of events per GeV, with the data points shown as black dots and their associated uncertainties. The stacked colored histograms indicate different processes: \(t\bar{t}\) (red), single top (orange), \(Z+b\) (dark blue), \(Z+c\) (medium blue), \(Z+\)light (light blue), diboson (gray), and fake lepton (purple). The hatched area represents the uncertainty in the predictions.

The figure shows the pre-fit distribution of events as a function of $H_{\mathrm{T}}^{\ell,j} = \sum_{\ell,j} p_{T}^{\ell,j}$, scalar sum of $p_T$ for all jets and leptons in the $\ell+$jets channel, in proton-lead (p+Pb) collisions at a center-of-mass energy of $\sqrt{s_{\mathrm{NN}}} = 8.16$ TeV, with an integrated luminosity of 165 nb\(^{-1}\) in the 2\(\ell\)2bincl signal region. The histogram represents the number of events per GeV, with the data points shown as black dots and their associated uncertainties. The stacked colored histograms indicate different processes: \(t\bar{t}\) (red), single top (orange), \(Z+b\) (dark blue), \(Z+c\) (medium blue), \(Z+\)light (light blue), diboson (gray), and fake lepton (purple). The hatched area represents the uncertainty in the predictions.

The figure shows the post-fit distribution of events as a function of $H_{\mathrm{T}}^{\ell,j} = \sum_{\ell,j} p_{T}^{\ell,j}$, scalar sum of $p_T$ for all jets and leptons in the $\ell+$jets channel, in proton-lead (p+Pb) collisions at a center-of-mass energy of $\sqrt{s_{\mathrm{NN}}} = 8.16$ TeV, with an integrated luminosity of 165 nb\(^{-1}\) in the 2\(\ell\)2bincl signal region. The histogram represents the number of events per GeV, with the data points shown as black dots and their associated uncertainties. The stacked colored histograms indicate different processes: \(t\bar{t}\) (red), single top (orange), \(Z+b\) (dark blue), \(Z+c\) (medium blue), \(Z+\)light (light blue), diboson (gray), and fake lepton (purple). The hatched area represents the uncertainty in the predictions.

The correlation matrix of the fit parameters for the six combined regions, including the signal regions (electron+jets: 1-lepton 1-bjet and 1-lepton 2-bjet inclusive, muon+jets: 1-lepton 1-bjet and 1-lepton 2-bjet inclusive, dilepton: 2-lepton 1-bjet and 2-lepton 2-bjet inclusive) is shown.

The impact of systematic uncertainties on the fitted signal-strength parameter $\hat{\mu}$ for the combined (six signal regions ($e+$jets: $1\ell1b$ and $1\ell2b\mathrm{incl}$, $\mu+$jets: $1\ell1b$ and $1\ell2b\mathrm{incl}$, dilepton: $2\ell1b$ and $2\ell2b\mathrm{incl}$) fit of all channels. Only the 15 most significant systematic uncertainties are shown and listed in decreasing order of their impact on $\mu$ on the $y$-axis. The empty (filled) blue/cyan boxes correspond to the pre-fit (post-fit) impact on $\mu$, referring to the upper $x$-axis. The impact of each systematic uncertainty, $\Delta \mu$, is calculated by comparing the nominal best-fit value of $\mu$ with the result of the fit when fixing the corresponding nuisance parameter $\theta$ to its best-fit value $\hat{\theta}$ shifted by its pre-fit (post-fit) uncertainties $\hat{\theta} \pm \Delta \theta(\hat{\theta} \pm \Delta \hat{\theta})$. The black points, which refer to the lower $x$-axis, show the pulls of the fitted nuisance parameters, i.e., the deviations of the fitted parameters $\hat{\theta}$ from their nominal values $\theta_0$, normalized to their nominal uncertainties $\Delta \theta$. The black lines show the post-fit uncertainties of the nuisance parameters, relative to their nominal uncertainties, which are indicated by the dashed lines.

Breakdown of relative systematic uncertainties in the measured ${t\bar{t}}$ cross-section. The quoted uncertainties are obtained by repeating the fit with a group of nuisance parameters fixed to their fitted values and subtracting in quadrature the resulting total uncertainty from the uncertainty of the complete fit. However, the total uncertainty is not the quadratic sum of the grouped impacts, as this approach neglects the correlation among the different groups.

The signal strength $\mu_{t\bar{t}}$, defined as the ratio of the observed signal for the combined six signal regions ($e+$jets: $1\ell1b$ and $1\ell2b\mathrm{incl}$, $\mu+$jets: $1\ell1b$ and $1\ell2b\mathrm{incl}$, dilepton: $2\ell1b$ and $2\ell2b\mathrm{incl}$.) to the SM expectation with no nPDF effects included, is measured using a binned profile-likelihood method. The parameter $\mu_{t\bar{t}}$ is determined by the fit to the $H_{T}^{\ell,j}$ data distributions in the six signal regions, where the $H_{T}^{\ell,j}$ variable is defined as the scalar sum of the transverse momenta of the leptons and HI jets. Over all 9% uncertainties is observed

The figure shows the measurement of the top quark pair production cross-section, $\sigma_{t\bar{t}}$, in proton-lead (p+Pb) collisions at a center-of-mass energy of $\sqrt{s_{\mathrm{NN}}} = 8.16$ TeV, based on 165 nb$^{-1}$ of data. The central measured value is shown alongside theoretical predictions using various parton distribution functions (PDFs): MCFM TUJU21, MCFM nNNPDF30, MCFM nCTEQ15HQ, and MCFM EPPS21. The green and yellow bands represent the total and statistical uncertainties of the data, respectively. Additional measurements for comparison include the CMS result at 8.16 TeV for p+Pb collisions and an extrapolated ATLAS+CMS result for pp collisions at 8 TeV. Relevant references are indicated: PRL 119 (2017) 242001 and JHEP 07 (2023) 213.

The ATLAS measurement of \( R_{p\text{Pb}} \) as a function of the nuclear modification factor in proton-lead (\( p\) + Pb) collisions at a center-of-mass energy per nucleon pair \( \sqrt{s_{NN}} = 8.16 \) TeV, with an integrated luminosity of 165 nb\(^{-1}\). The black points represent theoretical predictions from various MCFM parton distribution functions (PDFs): TUJU21, nNNPDF30, nCTEQ15HQ, and EPPS21. The green band indicates the total uncertainty of the data, while the yellow band represents the statistical uncertainty. The dashed line at \( R_{p\text{Pb}} = 1 \) shows the expected value in the absence of nuclear modifications.

Distribution of $v_{3}$ from MBT events plotted as a function of the A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of $v_{3}$ from $p_{T}^{jet}>75$ GeV events plotted as a function of the A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of $v_{3}$ from $p_{T}^{jet}>100$ GeV events plotted as a function of the A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of $v_{2}$ from MBT events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of $v_{2}$ from $p_{T}^{jet}>75$ GeV events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of $v_{2}$ from $p_{T}^{jet}>100$ GeV events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of $v_{2}$ from MBT events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of $v_{2}$ from $p_{T}^{jet}>75$ GeV events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of $v_{2}$ from $p_{T}^{jet}>100$ GeV events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of $v_{2}$ from MBT events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (9-100 GeV).

Distribution of $v_{2}$ from $p_{T}^{jet}>75$ GeV events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (9-100 GeV).

Distribution of $v_{2}$ from $p_{T}^{jet}>100$ GeV events plotted as a function of the event centrality and A-particle $p_\mathrm{T}$ in (9-100 GeV).

Distribution of relative UE-UE pair fractions from MBT events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative HS-HS pair fractions from MBT events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative HS-UE pair fractions from MBT events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative UE-HS pair fractions from MBT events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative UE-UE pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative HS-HS pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative HS-UE pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative UE-HS pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative UE-UE pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative HS-HS pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative HS-UE pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative UE-HS pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of A-particle $p_\mathrm{T}$ for 0-5% centrality.

Distribution of relative UE-UE pair fractions from MBT events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of relative HS-UE pair fractions from MBT events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of relative UE-UE pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of relative HS-UE pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of relative UE-UE pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of relative HS-UE pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (0.5-2 GeV).

Distribution of relative UE-UE pair fractions from MBT events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of relative HS-UE pair fractions from MBT events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of relative UE-UE pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of relative HS-UE pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of relative UE-UE pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of relative HS-UE pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (2-9 GeV).

Distribution of relative UE-UE pair fractions from MBT events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (9-100 GeV).

Distribution of relative HS-UE pair fractions from MBT events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (9-100 GeV).

Distribution of relative UE-UE pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (9-100 GeV).

Distribution of relative HS-UE pair fractions from $p_{T}^{jet}>75$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (9-100 GeV).

Distribution of relative UE-UE pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (9-100 GeV).

Distribution of relative HS-UE pair fractions from $p_{T}^{jet}>100$ GeV events plotted as a function of centrality A-particle $p_\mathrm{T}$ in (9-100 GeV).

Extrapolation factors AWW -- used to extrapolate from the fiducial particle level to the total phase space.

Cross sections and ratios (13/8 TeV) -- experimental and PDF uncertainties are treated as uncorrelated, scale uncertainties are correlated between the different centre-of-mass energies.