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 $ {{{\mathrm {B}}}\to {\mathrm {D^0}}} $ nuclear modification factor $ {R_{\mathrm {AA}}} $ for PbPb collisions at ${\sqrt {\smash [b]{s_{_{\mathrm {NN}}}}}} = $ 5.02 TeV.

The nuclear modification factor of charged particles having |eta|<1 in XeXe collisions at sqrt(s_NN)=5.44 TeV. The first systematic uncertainty describes uncertainties that are not fully correlated across points, while the second systematic uncertainty is a normalization uncertainty that is fully correlated across all points.

The nuclear modification factor of charged particles having |eta|<1 in XeXe collisions at sqrt(s_NN)=5.44 TeV. The first systematic uncertainty describes uncertainties that are not fully correlated across points, while the second systematic uncertainty is a normalization uncertainty that is fully correlated across all points.

The RXePb of charged particles having |eta|<1 in XeXe collisions at sqrt(s_NN)=5.44 TeV and PbPb collisions at sqrt(s_NN)=5.02 TeV. The first systematic uncertainty describes uncertainties that are not fully correlated across points, while the second systematic uncertainty is a normalization uncertainty that is fully correlated across all points. Bins where no data point has been reported are denoted as 'empty'.

The RXePb of charged particles having |eta|<1 in XeXe collisions at sqrt(s_NN)=5.44 TeV and PbPb collisions at sqrt(s_NN)=5.02 TeV. The first systematic uncertainty describes uncertainties that are not fully correlated across points, while the second systematic uncertainty is a normalization uncertainty that is fully correlated across all points.

The RXePb of charged particles having |eta|<1 in XeXe collisions at sqrt(s_NN)=5.44 TeV and PbPb collisions at sqrt(s_NN)=5.02 TeV. The first systematic uncertainty describes uncertainties that are not fully correlated across points, while the second systematic uncertainty is a normalization uncertainty that is fully correlated across all points.

The RXePb of charged particles having |eta|<1 in XeXe collisions at sqrt(s_NN)=5.44 TeV and PbPb collisions at sqrt(s_NN)=5.02 TeV. The first systematic uncertainty describes uncertainties that are not fully correlated across points, while the second systematic uncertainty is a normalization uncertainty that is fully correlated across all points.

The average N_part, N_coll, and TAA for various centrality classes of XeXe collisions at 5.44 TeV.

Invariant transverse momentum spectra of $\omega$ production in $p$+$p$ and $d$+Au collisions at $\sqrt{s}$=200 GeV.

Invariant transverse momentum spectra of $\omega$ production in $p$+$p$ and $d$+Au collisions at $\sqrt{s}$=200 GeV.

Invariant transverse momentum spectra of $\omega$ production in Cu+Cu (left) and Au+Au (right) collisions from the $\omega \rightarrow \pi^0 \gamma$ decay channel for three centrality bins and minimum bias.

Invariant transverse momentum spectra of $\omega$ production in Cu+Cu (left) and Au+Au (right) collisions from the $\omega \rightarrow \pi^0 \gamma$ decay channel for three centrality bins and minimum bias.

Invariant transverse momentum spectra of $\omega$ production in Cu+Cu (left) and Au+Au (right) collisions from the $\omega \rightarrow \pi^0 \gamma$ decay channel for three centrality bins and minimum bias.

Invariant transverse momentum spectra of $\omega$ production in Cu+Cu (left) and Au+Au (right) collisions from the $\omega \rightarrow \pi^0 \gamma$ decay channel for three centrality bins and minimum bias.

Invariant transverse momentum spectra of $\omega$ production in Cu+Cu (left) and Au+Au (right) collisions from the $\omega \rightarrow \pi^0 \gamma$ decay channel for three centrality bins and minimum bias.

Invariant transverse momentum spectra of $\omega$ production in Cu+Cu (left) and Au+Au (right) collisions from the $\omega \rightarrow \pi^0 \gamma$ decay channel for three centrality bins and minimum bias.

Measured ${\omega}/{\pi}$ ratio as a function of $p_T$ in $p$ + $p$ collisions at $\sqrt{s}$=200 GeV.

Measured ${\omega}/{\pi}$ ratio as a function of $p_T$ in $p$ + $p$ collisions at $\sqrt{s}$=200 GeV.

Measured ${\omega}/{\pi}$ ratio as a function of $p_T$ in $p$ + $p$ collisions at $\sqrt{s}$=200 GeV.

${\omega}/{\pi}$ ratio versus transverse momentum in $d$+Au (0–88%) for $\omega \rightarrow {e}^{+}{e}^{-}$.

${\omega}/{\pi}$ ratio versus transverse momentum in $d$+Au (0–88%) for $\omega \rightarrow {\pi}^{+}{\pi}^{-}\pi^0$.

${\omega}/{\pi}$ ratio versus transverse momentum in $d$+Au (0–88%) for $\omega \rightarrow {\pi}^{0}\gamma$.

${\omega}/{\pi}$ ratio versus transverse momentum in Cu+Cu (0–94%) for $\omega \rightarrow {\pi}^{0}\gamma$.

${\omega}/{\pi}$ ratio versus transverse momentum in Au+Au (0–92%) for $\omega \rightarrow {\pi}^{0}\gamma$.

${\omega}/{\pi}$ ratio versus transverse momentum in Au+Au (0–92%) for $\omega \rightarrow {\pi}^{0}\gamma$.

Nuclear modification factor, RdAu, measured for the ${R}_{dAu}$ in 0–88, 0–20, 20-40, and 40-60% centrality bins in $d$+Au collisions at $\sqrt{s}$=200 GeV.

Nuclear modification factor, RdAu, measured for the ${R}_{dAu}$ in 60–88% centrality bins in $d$+Au collisions at $\sqrt{s}$=200 GeV.

Nuclear modification factor, RdAu, measured for the ${R}_{dAu}$ in 0–88, 20-40, and 40–60% centrality bins in $d$+Au collisions at $\sqrt{s}$=200 GeV.

Nuclear modification factor, RdAu, measured for the ${R}_{dAu}$ in 0–20 and 60–88% centrality bins in $d$+Au collisions at $\sqrt{s}$=200 GeV.

Nuclear modification factor, RdAu, measured for the ${R}_{dAu}$ in 0–88, 0–20, 20-40, 40-60, and 60–88% centrality bins in $d$+Au collisions at $\sqrt{s}$=200 GeV.

$R_{AA}$ of the $\omega$ in Cu+Cu collisions from the $\omega \rightarrow \pi^0\gamma$ decay channel.

$R_{AA}$ of the $\omega$ in Au+Au collisions from the $\omega \rightarrow \pi^0 \gamma$ decay channel.

$R_{AA}$ of the $\omega$ in Au+Au collisions from the $\omega \rightarrow \pi^0 \gamma$ decay channel.

$R_{AA}$ of the $\omega$ in Au+Au collisions from the $\omega \rightarrow \pi^0 \gamma$ decay channel.

$R_{AA}$ of the $\omega$ in Au+Au collisions from the $\omega \rightarrow \pi^0 \gamma$ decay channel.

$R_{AA}$ of the $\omega$ in and Au+Au collisions from the $\omega \rightarrow \pi^0 \gamma$ decay channel.

($R_{AA}$) of the $\omega \rightarrow \pi^+\pi^-\pi^0$ collision for the $\omega$ meson integrated over the range $p_T$ > 7 $\mathrm{GeV}$/$c$ as a function of the number participating nucleons ($N_{part}$).

($R_{AA}$) of the $\omega \rightarrow \pi^0 \gamma$ collision for the $\omega$ meson integrated over the range $p_T$ > 7 $\mathrm{GeV}$/$c$ as a function of the number participating nucleons ($N_{part}$).

$R_{AA}$ of the $\omega \rightarrow \pi^0 \gamma$ collision for the $\omega$ meson integrated over the range $p_T$ > 7 $\mathrm{GeV}$/$c$ as a function of the number participating nucleons ($N_{part}$).

Measured $\pi^0$ $R_{AA}$ as a function of $p_T$ for the 0-10% most central collisions at $\sqrt{s_{NN}}$ = 22.4 GeV in comparison to a jet quenching calculation. The error (tot.) includes the quadratic sum of the statistical uncertainties and the point-to-point uncorrelated and correlated systematic uncertainties.

Measured $\pi^0$ $R_{AA}$ as a function of $p_T$ for the 0-10% most central collisions at $\sqrt{s_{NN}}$ = 62.4 GeV in comparison to a jet quenching calculation. The error (tot.) includes the quadratic sum of the statistical uncertainties and the point-to-point uncorrelated and correlated systematic uncertainties.

Measured $\pi^0$ $R_{AA}$ as a function of $p_T$ for the 0-10% most central collisions at $\sqrt{s_{NN}}$ = 200.0 GeV in comparison to a jet quenching calculation. The error (tot.) includes the quadratic sum of the statistical uncertainties and the point-to-point uncorrelated and correlated systematic uncertainties.

The average $R_{AA}$ in the interval 2.5 < $p_T$ < 3.5 GeV/$c$ as a function of centrality for Cu+Cu collisions at $\sqrt{s_{NN}}$ = 22.4 GeV. The error (sys.) includes the normalization and $<N_{coll}>$ uncertainties for a typical $N_{coll}$ uncertainty of 12%.

The average $R_{AA}$ in the interval 2.5 < $p_T$ < 3.5 GeV/$c$ as a function of centrality for Cu+Cu collisions at $\sqrt{s_{NN}}$ = 62.4 GeV. The error (sys.) includes the normalization and $<N_{coll}>$ uncertainties for a typical $N_{coll}$ uncertainty of 12%.

The average $R_{AA}$ in the interval 2.5 < $p_T$ < 3.5 GeV/$c$ as a function of centrality for Cu+Cu collisions at $\sqrt{s_{NN}}$ = 200 GeV. The error (sys.) includes the normalization and $<N_{coll}>$ uncertainties for a typical $N_{coll}$ uncertainty of 12%.

Neutral pion invariant yields as a function of $p_T$ measured in minimum bias and 9 centrality classes in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Neutral pion invariant yields as a function of $p_T$ measured in minimum bias and 9 centrality classes in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

$R_{AA}$ vs. $\Delta\phi$ for $\pi^0$ yields integrated over 3 < $p_T$ < 5 GeV/$c$.

$<R_{AA}>$ for $\pi^0$ yields integrated over 3 < $p_T$ < 5 GeV/$c$.

$R_{AA}$ vs. $\Delta\phi$ for $\pi^0$ yields integrated over 5 < $p_T$ < 8 GeV/$c$.

$<R_{AA}>$ for $\pi^0$ yields integrated over 5 < $p_T$ < 8 GeV/$c$.

$p_T$ spectra of charged hadrons for minimum bias collisions along with spectra for 9 centrality classes derived from the pseudo-rapidity region $|\eta|$ < 0.18. Stat. stands for statistical error, syst. stands for the systematic errors and occ. stands for occupancy error.

$p_T$ spectra of charged hadrons for minimum bias collisions along with spectra for 9 centrality classes derived from the pseudo-rapidity region $|\eta|$ < 0.18. Stat. stands for statistical error, syst. stands for the systematic errors and occ. stands for occupancy error.

Centrality dependence of $<p^{trunc}_T>$. Data points for truncated mean $p_T$ calculation for three $p_T$ ranges. The integral region is from 8 GeV/$c$ > $p_T$ > $p_{T_{min}}$ for the most peripheral bin. Since the data stop around 6.5 GeV/$c$, truncated mean $p_T$ is lower.

Ratio of charged hadron yields per nucleon-nucleon collision between central (0-10%) and peripheral (60-92%) Au + Au collisions. Stat. stands for statistical error, syst. stands for the systematic errors and occ. stands for occupancy error.

$R_{AA}$ for $(h^++h^-)/2$ and $\pi^0$ as function of $p_{T}$ for minimum bias and 9 centrality classes. Stat. stands for statistical error from Au+Au. The asymmetric uncertainty from the fit of the $\pi^0$ mass peak is denoted by fit. $\sigma_{ENBC}$ stands for error on number of binary collisions. Error from $pp$ parametrization that is not $p_T$ dependent. This error is not included in the calculation of the Upper bound of the systematic error and the Lower bound of the systematic error.

$R_{AA}$ for $(h^++h^-)/2$ and $\pi^0$ as function of $p_{T}$ for minimum bias and 9 centrality classes. Stat. stands for statistical error from Au+Au. The asymmetric uncertainty from the fit of the $\pi^0$ mass peak is denoted by fit. $\sigma_{ENBC}$ stands for error on number of binary collisions. Error from $pp$ parametrization that is not $p_T$ dependent. This error is not included in the calculation of the Upper bound of the systematic error and the Lower bound of the systematic error.

$R_{AA}$ for $(h^++h^-)/2$ and $\pi^0$ as function of $p_{T}$ for minimum bias and 9 centrality classes. Stat. stands for statistical error from Au+Au. The asymmetric uncertainty from the fit of the $\pi^0$ mass peak is denoted by fit. $\sigma_{ENBC}$ stands for error on number of binary collisions. Error from $pp$ parametrization that is not $p_T$ dependent. This error is not included in the calculation of the Upper bound of the systematic error and the Lower bound of the systematic error.

Charged hadron to $\pi^0$ ratios for minimum bias events and 9 centrality classes.

Charged hadron to $\pi^0$ ratios for minimum bias events and 9 centrality classes.

Charged hadron to $\pi^0$ ratios for minimum bias events and 9 centrality classes.

Charged hadron to $\pi^0$ ratios for minimum bias events and 9 centrality classes.

Charged hadron to $\pi^0$ ratios for minimum bias events and 9 centrality classes.

Charged hadron to $\pi^0$ ratios for minimum bias events and 9 centrality classes.

Au + Au yield integrated for $p_T$ > 4.5 (GeV/$c$) over the $N$ + $N$ yield, normalized using either $N_{coll}$ or $N_{part}$, plot as a function of $<N_{part}>$. Statistical and systematic error from Au+Au and $N$+$N$ reference includes: a) $NN$ reference error with 10.4% normalization error subtracted, this error is 11.2%, b) systematic error on Au+Au is 8.5% (correction), and c) 5% occupancy correction (it is centrality dependent). Error on $N_{coll}$+10.4% error on $N$+$N$ reference.

Au + Au yield integrated for $p_T$ > 4.5 (GeV/$c$) over the $N$ + $N$ yield, normalized using either $N_{coll}$ or $N_{part}$, plot as a function of $<N_{part}>$. Statistical and systematic error from Au+Au and $N$+$N$ reference includes: a) $NN$ reference error with 10.4% normalization error subtracted, this error is 11.2%, b) systematic error on Au+Au is 8.5% (correction), and c) 5% occupancy correction (it is centrality dependent). Error on $N_{coll}$+10.4% error on $N$+$N$ reference.

Invariant yields from the $p$+$p$ reference data. The normalization uncertainty is labeled norm, the uncertainty from the fit is labeled fit, and the uncertainty $\pi^0$ is from the $\pi^0$ ratio.