The unfolded dijet balance distribution, $(1/N_{dijet})(dN_{dijet}/dx_{j})$, as function of $x_{j}$ for the $10-60$, $60-120$, $120-185$, $185-250$ and $250-400$ multiplicity ranges with leading and subleading jets at forward and midrapidity regions, respectively.

The unfolded dijet balance distribution, $(1/N_{dijet})(dN_{dijet}/dx_{j})$, as function of $x_{j}$ for the $10-60$, $60-120$, $120-185$, $185-250$ and $250-400$ multiplicity ranges with leading and subleading jets at backward and midrapidity regions, respectively.

The ratio of unfolded dijet balance distribution to the lower multiplicity range, $10-60$, as function of $x_{j}$ for the $60-120$, $120-185$, $185-250$ and $250-400$ ranges with both jets at the midrapidity regions.

The ratio of unfolded dijet balance distribution to the lower multiplicity range, $10-60$, as function of $x_{j}$ for the $60-120$, $120-185$, $185-250$ and $250-400$ ranges with leading and subleading jets at midrapidity and forward regions, respectively.

The ratio of unfolded dijet balance distribution to the lower multiplicity range, $10-60$, as function of $x_{j}$ for the $60-120$, $120-185$, $185-250$ and $250-400$ ranges with leading and subleading jets at midrapidity and backward regions, respectively.

The ratio of unfolded dijet balance distribution to the lower multiplicity range, $10-60$, as function of $x_{j}$ for the $60-120$, $120-185$, $185-250$ and $250-400$ ranges with leading and subleading jets at forward and midrapidity regions, respectively.

The ratio of unfolded dijet balance distribution to the lower multiplicity range, $10-60$, as function of $x_{j}$ for the $60-120$, $120-185$, $185-250$ and $250-400$ ranges with leading and subleading jets at backward and midrapidity regions, respectively.

Mean xj ratio, $\langle x_{j} \rangle$, in the range $10-60$ as unfolded for different multiplicities, for different leading and subleading jet combinations for the midrapidity, forward, and backward regions

The reconstructed (raw) dijet balance distribution, $(1/N_{dijet})(dN_{dijet}/dx_{j})$, as function of $x_{j}$ for the $10-60$, $60-120$, $120-185$, $185-250$ and $250-400$ multiplicity ranges with both jets at the midrapidity regions.

The reconstructed (raw) dijet balance distribution, $(1/N_{dijet})(dN_{dijet}/dx_{j})$, as function of $x_{j}$ for the $10-60$, $60-120$, $120-185$, $185-250$ and $250-400$ multiplicity ranges with leading and subleading jets at midrapidity and forward regions, respectively.

The reconstructed (raw) dijet balance distribution, $(1/N_{dijet})(dN_{dijet}/dx_{j})$, as function of $x_{j}$ for the $10-60$, $60-120$, $120-185$, $185-250$ and $250-400$ multiplicity ranges with leading and subleading jets at midrapidity and backward regions, respectively.

The reconstructed (raw) dijet balance distribution, $(1/N_{dijet})(dN_{dijet}/dx_{j})$, as function of $x_{j}$ for the $10-60$, $60-120$, $120-185$, $185-250$ and $250-400$ multiplicity ranges with leading and subleading jets at forward and midrapidity regions, respectively.

The reconstructed (raw) dijet balance distribution, $(1/N_{dijet})(dN_{dijet}/dx_{j})$, as function of $x_{j}$ for the $10-60$, $60-120$, $120-185$, $185-250$ and $250-400$ multiplicity ranges with leading and subleading jets at backward and midrapidity regions, respectively.

The ratio of reconstructed (raw) dijet balance distribution to the lower multiplicity range, $10-60$, as function of $x_{j}$ for the $60-120$, $120-185$, $185-250$ and $250-400$ ranges with both jets at the midrapidity regions.

The ratio of reconstructed (raw) dijet balance distribution to the lower multiplicity range, $10-60$, as function of $x_{j}$ for the $60-120$, $120-185$, $185-250$ and $250-400$ ranges with leading and subleading jets at midrapidity and forward regions, respectively.

The ratio of reconstructed (raw) dijet balance distribution to the lower multiplicity range, $10-60$, as function of $x_{j}$ for the $60-120$, $120-185$, $185-250$ and $250-400$ ranges with leading and subleading jets at midrapidity and backward regions, respectively.

The ratio of reconstructed (raw) dijet balance distribution to the lower multiplicity range, $10-60$, as function of $x_{j}$ for the $60-120$, $120-185$, $185-250$ and $250-400$ ranges with leading and subleading jets at forward and midrapidity regions, respectively.

The ratio of reconstructed (raw) dijet balance distribution to the lower multiplicity range, $10-60$, as function of $x_{j}$ for the $60-120$, $120-185$, $185-250$ and $250-400$ ranges with leading and subleading jets at backward and midrapidity regions, respectively.

Mean xj ratio, $\langle x_{j} \rangle$, in the range $10-60$ as reconstructed (raw) for different multiplicities, for different leading and subleading jet combinations for the midrapidity, forward, and backward regions

Normalized $\sigma_{\psi(2S)}/\sigma_{J/\psi}$ in $6.5<p_T<30.0\,GeV$ and $ 1 < y_{CM} < 1.935$ as functions of normalized $\text{N}^{{\text{corr.}}}_{\text{track}}$

Normalized $\sigma_{\psi(2S)}/\sigma_{J/\psi}$ in $6.5<p_T<30.0\,GeV$ and $ -2.865 < y_{CM} < 1.935$ as functions of normalized $\text{N}^{{\text{corr.}}}_{\text{track}}$

Normalized $\sigma_{\psi(2S)}/\sigma_{J/\psi}$ in $3.0<p_T<6.5\,GeV$ and $ -2.865 < y_{CM} < -2 $ as functions of normalized $\text{N}^{{\text{corr.}}}_{\text{track}}$

Normalized $\sigma_{\psi(2S)}/\sigma_{J/\psi}$ in $3.0<p_T<6.5\,GeV$ and $ 1 < y_{CM} < 1.935 $ as functions of normalized $\text{N}^{{\text{corr.}}}_{\text{track}}$

The second-order Fourier sine coefficients of $\Lambda$ polarization along the beam direction as functions of $N_\mathrm{trk}^\mathrm{offline}$ in pPb collisions at 8.16 TeV extracted with HF based event planes. The HF-based event plane is reconstructed only with the HF detectors toward the Pb going direction.

The second-order Fourier sine coefficients of $\bar{\Lambda}$ polarization along the beam direction as functions of $N_\mathrm{trk}^\mathrm{offline}$ in pPb collisions at 8.16 TeV extracted with HF based event planes. The HF-based event plane is reconstructed only with the HF detectors toward the Pb going direction.

The second-order Fourier sine coefficients of $\Lambda+\bar{\Lambda}$ polarization along the beam direction as functions of $N_\mathrm{trk}^\mathrm{offline}$ in pPb collisions at 8.16 TeV extracted with HF based event planes. The HF-based event plane is reconstructed only with the HF detectors toward the Pb going direction.

The elliptic anisotropy $v_2\{4\}$ as a function of $\langle N_{trk}^{offline} \rangle$ with 4-subevent method compared between pPb collisions at 8.16 TeV and PbPb collisions at 5.02 TeV with $p_{T}^{POI}>6$ GeV.

The elliptic anisotropy $v_2\{4\}$ as a function of $\langle N_{trk}^{corrected} \rangle$ with 4-subevent method compared between pPb collisions at 8.16 TeV and PbPb collisions at 5.02 TeV with $p_{T}^{POI}>6$ GeV.

The differential cumulant $d_2\{4\}$ as a function of $p_T$ in generator level HIJING for $60 \le N_{trk}^{gen} <120$ are compared to the experimental pPb results with $60 \le N_{trk}^{offline} <120$ and $185 \le N_{trk}^{offline} <250$.

Measured inclusive fiducial H->ZZ->4l cross section as a function of the center-of-mass energy.

Postfit reconstructed distributions of the Higgs boson transverse momentum.

Postfit reconstructed distributions of the Higgs boson rapidity.

Higgs fiducial cross section in bins of pT of the Higgs boson.

Higgs fiducial cross section in bins of absolute value of rapidity of the Higgs boson.

Correlation matrix of fit in different final states

Correlation matrix of fit in pT(4l) bins

Correlation matrix of fit in abs(y(4l)) bins

Tuned intrinsic kT parameters ShowerHandler:IntrinsicPtGaussian in Herwig with the underlying-event tune CH2 at nucleon-nucleon center-of-mass energy from 38.8 GeV to 13 TeV.

Tuned intrinsic kT parameters ShowerHandler:IntrinsicPtGaussian in Herwig with the underlying-event tune CH3 at nucleon-nucleon center-of-mass energy from 38.8 GeV to 13 TeV.

Tuned intrinsic kT parameters BeamRemnants:PrimordialkThard in Pythia with the underlying-event tune CP5 and ISR starting scale SpaceShower:pT0Ref = 1 GeV at nucleon-nucleon center-of-mass energy from 38.8 GeV to 13 TeV.

Tuned intrinsic kT parameters ShowerHandler:IntrinsicPtGaussian in Herwig with the underlying-event tune CH3 and ISR starting scale SudakovCommon:pTmin=0.7 GeV at nucleon-nucleon center-of-mass energy from 38.8 GeV to 13 TeV.

Tuned intrinsic kT parameters BeamRemnants:PrimordialkThard in Pythia with the underlying-event tune CP5 at proton-proton center-of-mass energy 38.8 GeV versus the dilepton mass range of the process.

Tuned intrinsic kT parameters BeamRemnants:PrimordialkThard in Pythia with the underlying-event tune CP5 at proton-proton center-of-mass energy 8000 GeV versus the dilepton mass range of the process.

Tuned intrinsic kT parameters BeamRemnants:PrimordialkThard in Pythia with the underlying-event tune CP5 at proton-nucleon center-of-mass energy 8160 GeV versus the dilepton mass range of the process.

Tuned intrinsic kT parameters BeamRemnants:PrimordialkThard in Pythia with the underlying-event tune CP5 at proton-nucleon center-of-mass energy 13000 GeV versus the dilepton mass range of the process.

Tuned intrinsic kT parameters ShowerHandler:IntrinsicPtGaussian in Herwig with the underlying-event tune CH2 at proton-proton center-of-mass energy 38.8 GeV versus the dilepton mass range of the process.

Tuned intrinsic kT parameters ShowerHandler:IntrinsicPtGaussian in Herwig with the underlying-event tune CH2 at proton-proton center-of-mass energy 8000 GeV versus the dilepton mass range of the process.

Tuned intrinsic kT parameters ShowerHandler:IntrinsicPtGaussian in Herwig with the underlying-event tune CH2 at proton-nucleon center-of-mass energy 8160 GeV versus the dilepton mass range of the process.

Tuned intrinsic kT parameters ShowerHandler:IntrinsicPtGaussian in Herwig with the underlying-event tune CH2 at proton-proton center-of-mass energy 13000 GeV versus the dilepton mass range of the process.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-proton center-of-mass energy 38.8 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-proton center-of-mass energy 62 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-proton center-of-mass energy 200 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-anti-proton center-of-mass energy 1800 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-anti-proton center-of-mass energy 1960 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-proton center-of-mass energy 2760 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-proton center-of-mass energy 8000 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-nucleon center-of-mass energy 8160 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-proton center-of-mass energy 13000 GeV when fitting to the CMS measurement.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters among various generator settings at proton-proton center-of-mass energy 13000 GeV when fitting to the LHCb measurement.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-proton center-of-mass energy 38.8 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-proton center-of-mass energy 62 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-proton center-of-mass energy 200 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-anti-proton center-of-mass energy 1800 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-anti-proton center-of-mass energy 1960 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-proton center-of-mass energy 2760 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-proton center-of-mass energy 8000 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-nucleon center-of-mass energy 8160 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-proton center-of-mass energy 13000 GeV when fitting to the CMS measurement.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Pythia at proton-proton center-of-mass energy 13000 GeV when fitting to the LHCb measurement.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-proton center-of-mass energy 38.8 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-proton center-of-mass energy 62 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-proton center-of-mass energy 200 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-anti-proton center-of-mass energy 1800 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-anti-proton center-of-mass energy 1960 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-proton center-of-mass energy 2760 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-proton center-of-mass energy 8000 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-nucleon center-of-mass energy 8160 GeV.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-proton center-of-mass energy 13000 GeV when fitting to the CMS measurement.

Covariance of the data uncertainty contribution to the tuned intrinsic kT parameters between the two generator settings of Herwig at proton-proton center-of-mass energy 13000 GeV when fitting to the LHCb measurement.

The $R_{\mathrm{HL}}$ for $\mathrm{B}^+$ in the full $p_{\mathrm{T}}$ range and as a function of the multiplicity density. The three different points in the plot represent B+ events with multiplicities: $60 \leq$ $N_{\mathrm{ch}}$ $< 85$, $85 \leq$ $N_{\mathrm{ch}}$ $< 110$, and $110\leq$ $N_{\mathrm{ch}}$ $\leq 250$, respectively. The error bars correspond to the statistical uncertainty, and the boxes represent the sum in quadrature of systematic uncertainties.

Z differential cross section in $p_{\mathrm{T}}$ bins divided into classes of multiplicity. Vertical bars represent total uncertainties.

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.

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

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

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