$m_{\rm inv}$ dependences of $\Delta\gamma$ in 20-50$\%$ Au+Au collisions at 200 GeV.

The $\pi$ pair $\Delta\gamma$ at $m_{\rm inv} > m_{\rm inv}^{\rm low}$ for the data identified by TPC and more extended $p_T$ range identified by TPC and TOF.

The elliptic flow $v_2$ in bins of $q_2$ for the three centrality bins in the 20-50$\%$ range of Au+Au collisions at 200 GeV.

$m_{\rm inv}$ dependences of the $\Delta\gamma$ from ESE-selected event samples A (50$\%$ largest $q_{2}$) and B (50$\%$ smallest $q_{2}$) in 20-50$\%$ Au+Au collisions at 200 GeV.

$m_{\rm inv}$ dependences of the $\rm \Delta\gamma_{A} - \Delta\gamma_{B}$ in 20-50$\%$ Au+Au collisions at 200 GeV.

$\Delta\gamma_{\rm A}$ versus $\Delta\gamma_{\rm B}$ in 20-50$\%$ Au+Au collisions at 200 GeV. (The data are arranged as function of mass)

Inverse slope of mT fit of XI- and XIBAR+ produced in inelastic p+p interactions at 158 GeV/c.

The XIBAR+/Xi-ratio in inelastic p+p interactions at 158 GeV/c beam momentum.

Ratio of pT integrated yields versus rapidity of XIBAR+ and Xi- produced in inelastic p+p interactions at 158 GeV/c beam momentum.

Corrected inclusive charged-particle jet distributions in Au+Au collisions at 200 GeV for R=0.2, 0.3, and 0.4 in peripheral (60-80%) Au+Au collisions for pTlead,min = 7 GeV/c. The first uncertainty is statistical (symmetric), followed by shape uncertainty (asymmetric) and correlated uncertainty (asymmetric).

pTlead,min (7/5 GeV/c) ratios of inclusive charged-particle jet distributions in Au+Au collisions at 200 GeV for R=0.2, 0.3, and 0.4 in peripheral Au+Au collisions. The first uncertainty is statistical (symmetric), followed by shape uncertainty (asymmetric) and correlated uncertainty (asymmetric).

pTlead,min (7/5 GeV/c) ratios of inclusive charged-particle jet distributions in Au+Au collisions at 200 GeV for R=0.2, 0.3, and 0.4 in central Au+Au collisions. The first uncertainty is statistical (symmetric), followed by shape uncertainty (asymmetric) and correlated uncertainty (asymmetric).

R_CP (0-10%/60-80%) ratios of inclusive charged-particle jet distributions in Au+Au collisions at 200 GeV for R = 0.2, 0.3, and 0.4 for pTlead,min = 5 GeV/c. The first uncertainty is statistical (symmetric), followed by shape uncertainty (asymmetric) and correlated uncertainty (asymmetric). In addition, the systematic uncertainty for the T_AA normalization (30%) has to be added.

R_CP (0-10%/60-80%) ratios of inclusive charged-particle jet distributions in Au+Au collisions at 200 GeV for R = 0.2 and 0.3 for pTlead,min = 5 GeV/c. The first uncertainty is statistical (symmetric), followed by shape uncertainty (asymmetric) and correlated uncertainty (asymmetric). In the paper are shown all these three uncertainty sources (statistical, shape, correlated) summed in quadrature. The systematic uncertainty for the T_AA normalization (30%) has not been included (as for the other measurements compared with the STAR data).

R_AA for peripheral (60-80%) Au+Au collisions at 200 GeV for R = 0.2, 0.3, and 0.4 for pTlead,min = 5 GeV/c. The first uncertainty is statistical (symmetric), followed by shape uncertainty (asymmetric) and correlated uncertainty (asymmetric). In addition, the systematic uncertainty for the T_AA normalization (29%) and PYTHIA reference (22%, 20%, 18% for R = 0.2, 0.3 and 0.4) have to be added.

R_AA for central (0-10%) Au+Au collisions at 200 GeV for R = 0.2, 0.3, and 0.4 for pTlead,min = 5 GeV/c. The first uncertainty is statistical (symmetric), followed by shape uncertainty (asymmetric) and correlated uncertainty (asymmetric). In addition, the systematic uncertainty for the T_AA normalization (7%) and PYTHIA reference (22%, 20%, 18% for R = 0.2, 0.3 and 0.4) have to be added. The unbiased data points are the 4 highest for R=0.2 and 3 highest for R=0.3 and 0.4, resp.

Ratio of inclusive charged-particle jet yields for R=0.2/R=0.4 for peripheral (60-80%) Au+Au collisions at 200 GeV for pTlead,min = 5 GeV/c (only unbiased region). The first uncertainty is statistical (symmetric), followed by shape uncertainty (asymmetric) and correlated uncertainty (asymmetric). Please note that in the paper, the uncertainties shown are the quadrature sums of the uncertainties listed in the table.

Ratio of inclusive charged-particle jet yields for R=0.2/R=0.4 for central (0-10%) Au+Au collisions at 200 GeV for pTlead,min = 5 GeV/c (only unbiased region). The first uncertainty is statistical (symmetric), followed by shape uncertainty (asymmetric) and correlated uncertainty (asymmetric). Please note that in the paper, the uncertainties shown are the quadrature sums of the uncertainties listed in the table.

Fit of elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\pi^{+}\pi^{-}$ photoproduction off protons with a Soeding-inspired analytic function --- statistical correlations only. Only one half of the (symmetric) matrix is stored.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\pi^{+}\pi^{-}$ photoproduction off protons, differential in the dipion mass and in bins of $W$. The tabulated cross sections are $\gamma p$ cross sections but can be converted to $ep$ cross sections using the effective photon flux $\Phi_{\gamma/e}$.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\pi^{+}\pi^{-}$ photoproduction off protons, differential in the dipion mass and in bins of $W$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction off protons, in bins of the photon-proton energy $W$. The cross section is defined as the integral of the relativistic Breit Wigner resonance in the dipion mass over the range $2m_\pi<m_{\pi\pi}<1.53$ GeV. The tabulated cross sections are $\gamma p$ cross sections but can be converted to $ep$ cross sections using the effective photon flux $\Phi_{\gamma/e}$.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction off protons in bins of the photon-proton energy $W$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Fit of elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction cross sections off protons as a function of energy. Parameters with subscript "el" and "pd" correspond to elastic and proton-dissociative cross sections, respectively.

Fit of elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction cross sections off protons as a function of energy --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Fit of elastic $\rho^0(770)$ photoproduction cross section off protons as a function of energy (various experiments)

Fit of elastic $\rho^0(770)$ photoproduction cross sections off protons as a function of energy (various experiments) --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\pi^{+}\pi^{-}$ photoproduction cross section off protons, double-differential in the dipion mass and the momentum transfer $\vert t\vert$. The tabulated cross sections are $\gamma p$ cross sections but can be converted to $ep$ cross sections using the effective photon flux $\Phi_{\gamma/e}$.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\pi^{+}\pi^{-}$ photoproduction cross section off protons, double-differential in the dipion mass and the momentum transfer $\vert t\vert$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction off protons, single-differential in the momentum transfer $\vert t\vert$. The cross section is defined as the integral of the relativistic Breit Wigner resonance in the dipion mass over the range $2m_\pi<m_{\pi\pi}<1.53$ GeV. The tabulated cross sections are $\gamma p$ cross sections but can be converted to $ep$ cross sections using the effective photon flux $\Phi_{\gamma/e}$.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction off protons, single differential in the momentum transfer $\vert t\vert$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Fit of elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction single-differential cross sections off protons as a function of momentum transfer $t$. Parameters with subscript "el" and "pd" correspond to elastic and proton-dissociative cross sections, respectively.

Fit of elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction single-differential cross sections off protons as a function of momentum transfer $t$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\pi^{+}\pi^{-}$ photoproduction double-differenial cross section off protons, as a function of the dipion mass and the momentum transfer $\vert t\vert$, in bins of the photon-proton energy $W$ The tabulated cross sections are $\gamma p$ cross sections but can be converted to $ep$ cross sections using the effective photon flux $\Phi_{\gamma/e}$.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\pi^{+}\pi^{-}$ photoproduction double-differenial cross section off protons, as a function of the dipion mass and the momentum transfer $\vert t\vert$, in bins of the photon-proton energy $W$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction off protons, single-differential in the momentum transfer $\vert t\vert$, in bins of the photon-proton energy $W$. The cross section is defined as the integral of the relativistic Breit Wigner resonance in the dipion mass over the range $2m_\pi<m_{\pi\pi}<1.53$ GeV. The tabulated cross sections are $\gamma p$ cross sections but can be converted to $ep$ cross sections using the effective photon flux $\Phi_{\gamma/e}$.

Elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction cross section off protons, single differential in the momentum transfer $\vert t\vert$ and in bins of the photon-proton energy $W$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Regge fit of elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction single-differential cross sections off protons as a function of momentum transfer $t$ and photon-proton energy $W$. Parameters with subscript "el" and "pd" correspond to elastic and proton-dissociative cross sections, respectively.

Regge fit of elastic ($m_Y=m_p$) and proton-dissociative ($1<m_Y<10$ GeV) $\rho^0(770)$ photoproduction single-differential cross sections off protons as a function of momentum transfer $t$ and photon-proton energy $W$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Fit of $b$-slopes in elastic ($m_Y=m_p$) $\rho^0(770)$ photoproduction single-differential cross sections off protons as a function of momentum transfer $t$ in bins of photon-proton energy $W$.

Fit of $b$-slopes in elastic ($m_Y=m_p$) $\rho^0(770)$ photoproduction single-differential cross sections off protons as a function of momentum transfer $t$ in bins of photon-proton energy $W$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

Fit of Pomeron trajectories in elastic ($m_Y=m_p$) $\rho^0(770)$ photoproduction single-differential cross sections off protons as a function of momentum transfer $t$ in bins of photon-proton energy $W$.

Fit of Pomeron trajectories in elastic ($m_Y=m_p$) $\rho^0(770)$ photoproduction single-differential cross sections off protons as a function of momentum transfer $t$ in bins of photon-proton energy $W$ --- statistical correlations coefficients $\rho_{ij}$ only. Only one half of the (symmetric) matrix is stored. Bins are identified by their global bin number.

$\pi^0$ spectra from figure 2a from 40-60% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\pi^0$ spectra from figure 2a from 60-80% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\eta$ spectra from figure 2b from minimum bias U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\eta$ spectra from figure 2b from 0-20% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\eta$ spectra from figure 2b from 20-40% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\eta$ spectra from figure 2b from 40-60% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\eta$ spectra from figure 2b from 60-80% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

Ratio of $\pi^0$ yield to the fit from figure 2c from minimum bias U+U collisions. Sys. uncertainties include B and C uncertainties from data points.

Ratio of $\pi^0$ yield to the fit from figure 2c from 0-20% U+U collisions. Sys. uncertainties include B and C uncertainties from data points.

Ratio of $\pi^0$ yield to the fit from figure 2c from 60-80% U+U collisions. Sys. uncertainties include B and C uncertainties from data points.

Ratio of $\eta$ yield to the fit from figure 2d from minimum bias U+U collisions. Sys. uncertainties include B and C uncertainties from data points.

Ratio of $\eta$ yield to the fit from figure 2d from 0-20% U+U collisions. Sys. uncertainties include B and C uncertainties from data points.

Ratio of $\eta$ yield to the fit from figure 2d from 40-60% U+U collisions. Sys. uncertainties include B and C uncertainties from data points.

$\eta/\pi^0$ ratio from figure 3 from minimum bias U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point.

$\eta/\pi^0$ ratio from figure 3 from 0-20% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point.

$\eta/\pi^0$ ratio from figure 3 from 20-40% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point.

$\eta/\pi^0$ ratio from figure 3 from 40-60% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point.

$\eta/\pi^0$ ratio from figure 3 from 60-80% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point.

$\eta/\pi^0$ ratio from figure 3 from minimum bias Au+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point.

$\pi^0$ $R_{AA}$ from figure 4a from 0-20% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data. Calculation is carried out for Glauber 1 set.

$\eta$ $R_{AA}$ from figure 4a from 0-20% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data. Calculation is carried out for Glauber 1 set.

$\pi^0$ $R_{AA}$ from figure 4b from 20-40% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data. Calculation is carried out for Glauber 1 set.

$\eta$ $R_{AA}$ from figure 4b from 20-40% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data. Calculation is carried out for Glauber 1 set.

$\pi^0$ $R_{AA}$ from figure 4c from 40-60% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data. Calculation is carried out for Glauber 1 set.

$\eta$ $R_{AA}$ from figure 4c from 40-60% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data. Calculation is carried out for Glauber 1 set.

$\pi^0$ $R_{AA}$ from figure 4d from minimum bias U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data. Calculation is carried out for Glauber 1 set.

$\eta$ $R_{AA}$ from figure 4d from minimum bias U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data. Calculation is carried out for Glauber 1 set.

$\pi^0$ $R_{AA}$ from figure 5a from 0-20% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data. Calculation is carried out for Glauber 2 set.

$\pi^0$ $R_{AA}$ from figure 5a from 0-5% Au+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\pi^0$ $R_{AA}$ from figure 5a from 40-60% U+U collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data. Calculation is carried out for Glauber 2 set.

$\pi^0$ $R_{AA}$ from figure 5a from 40-50% Au+Au collisions. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\pi^0$ integrated $R_{AA}$ for $p_T>5$ GeV/$c$ in U+U collisions as a function of centrality from figure 6a. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data. Calculation is carried out for Glauber 1 set.

$\eta$ integrated $R_{AA}$ for $p_T>5$ GeV/$c$ in U+U collisions as a function of centrality from figure 6a. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data. Calculation is carried out for Glauber 1 set.

$\pi^0$ integrated $R_{AA}$ for $p_T>5$ GeV/$c$ in Au+Au collisions as a function of centrality from figure 6a. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data.

$\eta$ integrated $R_{AA}$ for $p_T>5$ GeV/$c$ in U+U collisions as a function of centrality from figure 6b. Type A uncertainties are uncorrelated point-to-point. Type B uncertainties are correlated point-to-point. Type C uncertainties affect the scale of the data. Calculation is carried out for Glauber 2 set.

Number of participant nucleons of in dep U$+$U collisions

Number of participant nucleons of in dep U$+$U collisions

Number of participant nucleons of in dep Au$+$Au collisions

Dielectron excess spectra as in pp collisions at $\sqrt{s}$ = 13 TeV as a function of $m_{\rm ee}$ after subtraction of the hadronicdecay cocktail. Electrons are measured within $|\eta_{\rm e}| < 0.8$ and $p_{\rm T,e} > 0.075$ GeV/$c$.

Dielectron excess spectra as in pp collisions at $\sqrt{s}$ = 13 TeV in the invariant mass interval 0.15<$m_{\rm ee}$<0.6 as a function of $p_{\rm T,ee}$ after subtraction of the hadronicdecay cocktail. Electrons are measured within $|\eta_{\rm e}| < 0.8$ and $p_{\rm T,e} > 0.075$ GeV/$c$.

Integrated away-side $\gamma_{dir}$-h per-trigger $I_{AA}$ and $I_{dA}$ of Au+Au, d+Au, and p+p, as a function of $\xi$.

$I_{AA}$ vs $\xi$ for direct photon $p_{T}^{\gamma}$ of 5-7 GeV/$c$, 7-9 GeV/$c$, and 9-12 GeV/$c$.

$I_{AA}$ vs $\xi$ for direct photon $p_{T}^{\gamma}$ of 5-7 GeV/$c$, 7-9 GeV/$c$, and 9-12 GeV/$c$.

$I_{AA}$ vs $\xi$ for direct photon $p_{T}^{\gamma}$ of 5-7 GeV/$c$, 7-9 GeV/$c$, and 9-12 GeV/$c$.

$I_{AA}$ as a function of $\xi$ for direct photon $p_{T}^{\gamma}$ of 5-7, 7-9, and 9-12 GeV/$c$. Three away-side integration ranges are chosen to calculate the per-trigger yield and the corresponding $I_{AA}$, $|\Delta\phi-\pi|<\pi/2$, $|\Delta\phi-\pi|<\pi/3$, and $|\Delta\phi-\pi|<\pi/6$.

$I_{AA}$ as a function of $\xi$ for direct photon $p_{T}^{\gamma}$ of 5-7, 7-9, and 9-12 GeV/$c$. Three away-side integration ranges are chosen to calculate the per-trigger yield and the corresponding $I_{AA}$, $|\Delta\phi-\pi|<\pi/2$, $|\Delta\phi-\pi|<\pi/3$, and $|\Delta\phi-\pi|<\pi/6$.

$I_{AA}$ as a function of $\xi$ for direct photon $p_{T}^{\gamma}$ of 5-7, 7-9, and 9-12 GeV/$c$. Three away-side integration ranges are chosen to calculate the per-trigger yield and the corresponding $I_{AA}$, $|\Delta\phi-\pi|<\pi/2$, $|\Delta\phi-\pi|<\pi/3$, and $|\Delta\phi-\pi|<\pi/6$.

Ratios of $I_{AA}$ as a function of direct photon $p_T$ for three different away-side integration ranges.

Measured $I_{AA}$ for direct photon $p_T$ of 5-7, 7-9, and 9-12 GeV/$c$, as a function of $\xi$, are compared with theoretical model calculations.

Measured $I_{AA}$ for direct photon $p_T$ of 5-7, 7-9, and 9-12 GeV/$c$, as a function of $\xi$, are compared with theoretical model calculations.

Measured $I_{AA}$ for direct photon $p_T$ of 5-7, 7-9, and 9-12 GeV/$c$, as a function of $\xi$, are compared with theoretical model calculations.

$b\bar{b}$ invariant yield measured via B meson decay to like-sign dimuons as a function of pair $p_{T}$. The statistical and systematic uncertainties vary point-to-point and are listed for each measured value. An additional global systematic uncertainty is provided in each column heading, which applies to all data points per column.

$\lambda_{\phi}$ measured in $J/\psi$ transverse momentum bins of 0.0 < $p_T$ < 3.0 GeV/$c$ and 3.0 < $p_T$ < 10.0 GeV/$c$ overlaid with NRQCD predictions in the Helicity and Collins-Soper frames.

$\lambda_{\theta \phi}$ measured in $J/\psi$ transverse momentum bins of 0.0 < $p_T$ < 3.0 GeV/$c$ and 3.0 < $p_T$ < 10.0 GeV/$c$ overlaid with NRQCD predictions in the Helicity and Collins-Soper frames.

$\lambda_{\theta \phi}$ measured in $J/\psi$ transverse momentum bins of 0.0 < $p_T$ < 3.0 GeV/$c$ and 3.0 < $p_T$ < 10.0 GeV/$c$ overlaid with NRQCD predictions in the Helicity and Collins-Soper frames.

$\bar{\lambda}$ measured in $J/\psi$ transverse momentum bins of 0.0 < $p_T$ < 3.0 GeV/$c$ and 3.0 < $p_T$ < 10.0 GeV/$c$ overlaid with CSM and NRQCD predictions in the Helicity frame.

$\bar{\lambda}$ measured in $J/\psi$ transverse momentum bins of 0.0 < $p_T$ < 3.0 GeV/$c$ and 3.0 < $p_T$ < 10.0 GeV/$c$ overlaid with CSM and NRQCD predictions in the Helicity frame.

Invariant $F$ measured in $J/\psi$ transverse momentum bins of 0.0 < $p_T$ < 3.0 GeV/$c$ and 3.0 < $p_T$ < 10.0 GeV/$c$ overlaid with CSM and NRQCD predictions in the Helicity frame.

Invariant $F$ measured in $J/\psi$ transverse momentum bins of 0.0 < $p_T$ < 3.0 GeV/$c$ and 3.0 < $p_T$ < 10.0 GeV/$c$ overlaid with CSM and NRQCD predictions in the Helicity frame.

Event yields for a simulated $ZH(\rightarrow aa\rightarrow (b\bar{b})(b\bar{b}))$ sample with $m_a = 17.5\,\text{GeV}$. The signal sample is produced with cross section equals to the standard model $pp\to ZH$, i.e. $0.88\,\text{pb}$. Cut 0 corresponds to the initial number of events. Cut 1 requires the single lepton trigger. Cut 2 requires 2 identified leptons. Cut 3 requires the Z-boson mass window. Cut 4 requires 2 reconstructed $a\to b\bar{b}$ candidates. Cut 5a requires 2 identified $a\to b\bar{b}$ candidates in the 1HPC1LPC region. Cut 6a requires the 2 $a\to b\bar{b}$ candidates in the 1HPC1LPC region to be inside the Higgs mass window. Cut 5b requires 2 identified $a\to b\bar{b}$ candidates in the 2HPC region. Cut 6b requires the 2 $a\to b\bar{b}$ candidates in the 2HPC region to be inside the Higgs mass window.

Background yield table for Z+jets, $t\bar{t}$, and rare sources. Observed data yield. Signal $ZH(\rightarrow aa\rightarrow (b\bar{b})(b\bar{b}))$ yield with $m_a = 20\,\text{GeV}$. The signal sample is produced with cross section equals to the standard model $pp\to ZH$, i.e. $0.88\,\text{pb}$. The table includes the yields in two signal regions with leptons consistent with an on-shell Z-boson decay, one with 2 $a\to b\bar{b}$ candidates in the 2HPC region and one with 2 $a\to b\bar{b}$ candidates in the 1HPC1LPC region. The table also includes the yields in four control regions, one with leptons consistent with an on-shell Z-boson decay and 2 $a\to b\bar{b}$ candidates in the Low Purity Category (LPC), and three others where the leptons are not consistent an on-shell Z-boson decay.