The performance of jet energy resolution (JER) for jets with |y| < 2.1 evaluated as a function of pT_truth in different centrality bins. Simulated hard scatter events were overlaid onto events from a dedicated sample of minimum-bias Xe+Xe data. The fit parameters are listed in a sperate table (Extras 1)

The relative magnitude of systematic uncertainties for per-pair normalized xJ distributions in 0-10% Xe+Xe centrality

The relative magnitude of systematic uncertainties for absolutely normalized xJ distributions in 0-10% Xe+Xe centrality

The relative magnitude of systematic uncertainties for rho distributions for leading jets in 0-10% Xe+Xe centrality

Per-pair normalized xJ distribution evaluated in four centrality intervals and given pT1 interval.

Per-pair normalized xJ distribution evaluated in four centrality intervals and given pT1 interval.

Per-pair normalized xJ distribution evaluated in four centrality intervals and given pT1 interval.

Absolutely normalized xJ distribution evaluated in four centrality intervals and given pT1 interval.

Absolutely normalized xJ distribution evaluated in four centrality intervals and given pT1 interval.

Absolutely normalized xJ distribution evaluated in four centrality intervals and given pT1 interval.

Per-pair normalized xJ distribution in Xe+Xe collisions.

Per-pair normalized xJ distribution in Pb+Pb collisions.

Per-pair normalized xJ distribution in Xe+Xe collisions.

Per-pair normalized xJ distribution in Pb+Pb collisions.

Per-pair normalized xJ distribution in Xe+Xe collisions.

Per-pair normalized xJ distribution in Pb+Pb collisions.

Absolutely normalized xJ distribution in Xe+Xe collisions.

Absolutely normalized xJ distribution in Xe+Xe collisions.

Absolutely normalized xJ distribution in Xe+Xe collisions.

Absolutely normalized xJ distribution in Xe+Xe collisions.

Absolutely normalized xJ distribution in Xe+Xe collisions.

Absolutely normalized xJ distribution in Xe+Xe collisions.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of leading jet pT in the same centrality intervals.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of leading jet pT in the same centrality intervals.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of leading jet pT in the same centrality intervals.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of subleading jet pT in the same centrality intervals.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of subleading jet pT in the same centrality intervals.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of subleading jet pT in the same centrality intervals.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of leading jet pT in the same SUM ETFCal intervals (selecting equivalent event activity)

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of leading jet pT in the same SUM ETFCal intervals (selecting equivalent event activity)

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of leading jet pT in the same SUM ETFCal intervals (selecting equivalent event activity)

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of subleading jet pT in the same SUM ETFCal intervals (selecting equivalent event activity)

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of subleading jet pT in the same SUM ETFCal intervals (selecting equivalent event activity)

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of subleading jet pT in the same SUM ETFCal intervals (selecting equivalent event activity)

Parameter a,b, and c from JER fits in Figure 1b.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Probability of finding at least one TAR jet, where the p_T-leading TAR jet passes the m_Wcand and D₂^&beta;=1 requirements, as a function of m_s. The probability is determined in a sample of signal events with m_Z'=2100 GeV, with the preselections applied.

Probability of finding a W_had candidate reconstructed as a pair of R=0.4 PFlow jets, as a function of m_s. The probability is determined in a sample of signal events with m_Z'=500 GeV, with the preselections applied that do not pass the requirements of the merged category.

Probability of finding a W_had candidate reconstructed as a pair of R=0.4 PFlow jets, as a function of m_s. The probability is determined in a sample of signal events with m_Z'=1000 GeV, with the preselections applied that do not pass the requirements of the merged category.

Probability of finding a W_had candidate reconstructed as a pair of R=0.4 PFlow jets, as a function of m_s. The probability is determined in a sample of signal events with m_Z'=1700 GeV, with the preselections applied that do not pass the requirements of the merged category.

Probability of finding a W_had candidate reconstructed as a pair of R=0.4 PFlow jets, as a function of m_s. The probability is determined in a sample of signal events with m_Z'=2100 GeV, with the preselections applied that do not pass the requirements of the merged category.

Observed exclusion contour at 95&#37; C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_&chi;=1, m_&chi;=200 GeV, and sin&theta;=0.01.

Expected exclusion contour at 95&#37; C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_&chi;=1, m_&chi;=200 GeV, and sin&theta;=0.01.

Expected+1&sigma; exclusion contour at 95&#37; C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_&chi;=1, m_&chi;=200 GeV, and sin&theta;=0.01.

Expected-1&sigma; exclusion contour at 95&#37; C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_&chi;=1, m_&chi;=200 GeV, and sin&theta;=0.01.

Expected+2&sigma; exclusion contour at 95&#37; C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_&chi;=1, m_&chi;=200 GeV, and sin&theta;=0.01.

Expected-2&sigma; exclusion contour at 95&#37; C.L. for the dark Higgs model in the (m_Z', m_s) plane for g_q=0.25, g_&chi;=1, m_&chi;=200 GeV, and sin&theta;=0.01.

Observed upper limits at 95&#37; C.L. on &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) for m_Z'=0.5 TeV as a function of m_s. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) process for m_Z'=0.5 TeV, shown in dashed blue.

Observed upper limits at 95&#37; C.L. on &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) for m_Z'=1 TeV as a function of m_s. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) process for m_Z'=1 TeV, shown in dashed blue.

Observed upper limits at 95&#37; C.L. on &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) for m_Z'=1.7 TeV as a function of m_s. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) process for m_Z'=1.7 TeV, shown in dashed blue.

Observed upper limits at 95&#37; C.L. on &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) for m_Z'=2.1 TeV as a function of m_s. The expected limits, varied up and down by one and two standard deviations, are shown as green and yellow bands, respectively. The observed and expected limits are compared to the theoretical LO cross section for the &sigma;(pp &rarr; s &chi;&chi;) &times; B(s &rarr; W^&pm; W^&mp;) process for m_Z'=2.1 TeV, shown in dashed blue.

Data overlaid on SM background yields stacked in each SR and CR category after the fit to data ('Post-fit'). The yields in the SR are broken down into their contributions to the individual bins. The maximum-likelihood estimators are set to the conditional values of the CR-only fit, and propagated to SR and CRs.

Dominant sources of uncertainty for three dark Higgs scenarios after the fit to data. The uncertainties are quantified in terms of their contribution to the fitted signal uncertainty that is expressed relative to the theory prediction. Three representative dark Higgs signal scenarios with g_q=0.25, g_&chi;=1.0, sin&theta;=0.01 and m_&chi;=200 GeV are considered, which are indicated using the (m_Z', m_s) format in units of GeV in the table columns.

Cumulative efficiencies in the merged category for three representative dark Higgs signal scenarios with g_q=0.25, g_&chi;=1.0, sin&theta;=0.01, m_Z' = 1 TeV, and m_&chi;=200 GeV considering s&rarr;W(&ell;&nu;)W(qq) decays only.

Cumulative efficiencies in the resolved category for three representative dark Higgs signal scenarios with g_q=0.25, g_&chi;=1.0, sin&theta;=0.01, m_Z' = 1 TeV, and m_&chi;=200 GeV considering s&rarr;W(&ell;&nu;)W(qq) decays only.

Theoretical cross section for &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) for each of the dark Higgs signal points at m_Z′ ={300, 350, 400, 500, 750} GeV, with g_q = 0.25, g<sub>&chi; = 1.0, sin&theta; = 0.01, m_Z′ = 1 TeV , and m_&chi; = 200 GeV. Also shown are experimentally excluded cross sections of &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) (Obs.) together with the expectations (Exp.) varied up and down by one standard deviation (&plusmn;1&sigma;).

Theoretical cross section for &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) for each of the dark Higgs signal points at m_Z′ ={1000, 1700} GeV, with g_q = 0.25, g<sub>&chi; = 1.0, sin&theta; = 0.01, m_Z′ = 1 TeV , and m_&chi; = 200 GeV. Also shown are experimentally excluded cross sections of &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) (Obs.) together with the expectations (Exp.) varied up and down by one standard deviation (&plusmn;1&sigma;).

Theoretical cross section for &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) for each of the dark Higgs signal points at m_Z′ ={2100, 2500, 2900, 3300} GeV, with g_q = 0.25, g<sub>&chi; = 1.0, sin&theta; = 0.01, m_Z′ = 1 TeV , and m_&chi; = 200 GeV. Also shown are experimentally excluded cross sections of &sigma;(pp &rarr; s&chi;&chi;) &times; B(s &rarr; W^&plusmn;W^&#8723;) (Obs.) together with the expectations (Exp.) varied up and down by one standard deviation (&plusmn;1&sigma;).

The unfolded differential charge asymmetry as a function of the transverse momentum of the top pair system. The measured values are given with statistical and systematic uncertainties. The SM theory predictions calculated at NNLO in QCD and NLO in EW theory are listed. The scale uncertainty is obtained by varying renormalisation and factorisation scales independently by a factor of 2 or 0.5 around $\mu_0$ to calculate the maximum and minimum value of the asymmetry, respectively. The nominal value $\mu_0$ is chosen as $H_T/4$. The variations in which one scale is multiplied by 2 while the other scale is divided by 2 are excluded. Finally, the scale and MC integration uncertainties are added in quadrature.

The unfolded differential charge asymmetry as a function of the longitudinal boost of the top pair system. The measured values are given with statistical and systematic uncertainties. The SM theory predictions calculated at NNLO in QCD and NLO in EW theory are listed. The scale uncertainty is obtained by varying renormalisation and factorisation scales independently by a factor of 2 or 0.5 around $\mu_0$ to calculate the maximum and minimum value of the asymmetry, respectively. The nominal value $\mu_0$ is chosen as $H_T/4$. The variations in which one scale is multiplied by 2 while the other scale is divided by 2 are excluded. Finally, the scale and MC integration uncertainties are added in quadrature.

The unfolded inclusive leptonic asymmetry. The unfolded $A_C^{\ell\bar{\ell}}$ is obtained in the reduced phase-space defined by the requirement $|\Delta |\eta_{\ell\bar{\ell}}||<2.5$. The measured values are given with statistical and systematic uncertainties. The SM theory predictions calculated at NLO in QCD and NLO in EW theory are listed. The theory uncertainty is obtained by varying both scales by a factor of 0.5 or 2.0 to calculate the minimum and maximum value of the asymmetry, respectively.

The unfolded differential leptonic asymmetry as a function of the invariant mass of the di-lepton pair. The unfolded $A_C^{\ell\bar{\ell}}$ is obtained in the reduced phase-space defined by the requirement $|\Delta |\eta_{\ell\bar{\ell}}||<2.5$. The measured values are given with statistical and systematic uncertainties. The SM theory predictions calculated at NLO in QCD and NLO in EW theory are listed. The theory uncertainty is obtained by varying both scales by a factor of 0.5 or 2.0 to calculate the minimum and maximum value of the asymmetry, respectively.

The unfolded differential leptonic asymmetry as a function of the transverse momentum of the di-lepton pair. The unfolded $A_C^{\ell\bar{\ell}}$ is obtained in the reduced phase-space defined by the requirement $|\Delta |\eta_{\ell\bar{\ell}}||<2.5$. The measured values are given with statistical and systematic uncertainties. The SM theory predictions calculated at NLO in QCD and NLO in EW theory are listed. The theory uncertainty is obtained by varying both scales by a factor of 0.5 or 2.0 to calculate the minimum and maximum value of the asymmetry, respectively.

The unfolded differential leptonic asymmetry as a function of the longitudinal boost of the di-lepton pair. The unfolded $A_C^{\ell\bar{\ell}}$ is obtained in the reduced phase-space defined by the requirement $|\Delta |\eta_{\ell\bar{\ell}}||<2.5$. The measured values are given with statistical and systematic uncertainties. The SM theory predictions calculated at NLO in QCD and NLO in EW theory are listed. The theory uncertainty is obtained by varying both scales by a factor of 0.5 or 2.0 to calculate the minimum and maximum value of the asymmetry, respectively.

Individual 68% and 95% CL bounds on the relevant Wilson coefficients of the SM Effective Field Theory in units of $\text{TeV}^{-2}$. The bounds are derived from the $A_C^{t\bar{t}}$ inclusive measurement. The experimental uncertainties are accounted for, in the form of the complete covariance matrix that keeps track of correlations between bins for the differential measurement. The theory uncertainty from the NNLO QCD + NLO EW calculation is included by explicitly varying the renormalization and factorization scales, or the parton density functions, in the calculation and registering the variations in the intervals.

Individual 68% and 95% CL bounds on the relevant Wilson coefficients of the SM Effective Field Theory in units of $\text{TeV}^{-2}$. The bounds are derived from the $A_C^{t\bar{t}}$ vs $m_{t\bar{t}}$ measurement. The experimental uncertainties are accounted for, in the form of the complete covariance matrix that keeps track of correlations between bins for the differential measurement. The theory uncertainty from the NNLO QCD + NLO EW calculation is included by explicitly varying the renormalization and factorization scales, or the parton density functions, in the calculation and registering the variations in the intervals.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ inclusive measurement. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $\beta_{z,t\bar{t}}$ measurement for $\beta_{z,t\bar{t}}$ $\in$ [0,0.3]. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $\beta_{z,t\bar{t}}$ measurement for $\beta_{z,t\bar{t}}$ $\in$ [0.3,0.6]. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $\beta_{z,t\bar{t}}$ measurement for $\beta_{z,t\bar{t}}$ $\in$ [0.6,0.8]. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $\beta_{z,t\bar{t}}$ measurement for $\beta_{z,t\bar{t}}$ $\in$ [0.8,1]. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $m_{t\bar{t}}$ measurement for $m_{t\bar{t}}$ < 500 GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $m_{t\bar{t}}$ measurement for $m_{t\bar{t}}$ $\in$ [500,750] GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $m_{t\bar{t}}$ measurement for $m_{t\bar{t}}$ $\in$ [750,1000] GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $m_{t\bar{t}}$ measurement for $m_{t\bar{t}}$ $\in$ [1000,1500] GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $m_{t\bar{t}}$ measurement for $m_{t\bar{t}}$ > 1500 GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $p_{T,t\bar{t}}$ measurement for $p_{T,t\bar{t}}$ < 30 GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $p_{T,t\bar{t}}$ measurement for $p_{T,t\bar{t}}$ $\in$ [30,120] GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{t\bar{t}}$ vs $p_{T,t\bar{t}}$ measurement for $p_{T,t\bar{t}}$ > 120 GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ inclusive measurement. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $\beta_{z,\ell\bar{\ell}}$ measurement for $\beta_{z,\ell\bar{\ell}}$ $\in$[0,0.3]. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $\beta_{z,\ell\bar{\ell}}$ measurement for $\beta_{z,\ell\bar{\ell}}$ $\in$[0.3,0.6]. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $\beta_{z,\ell\bar{\ell}}$ measurement for $\beta_{z,\ell\bar{\ell}}$ $\in$[0.6,0.8]. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $\beta_{z,\ell\bar{\ell}}$ measurement for $\beta_{z,\ell\bar{\ell}}$ $\in$[0.8,1]. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $m_{\ell\bar{\ell}}$ measurement for $m_{\ell\bar{\ell}}$ < 200 GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $m_{\ell\bar{\ell}}$ measurement for $m_{\ell\bar{\ell}}$ $\in$ [200,300] GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $m_{\ell\bar{\ell}}$ measurement for $m_{\ell\bar{\ell}}$ $\in$ [300,400] GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $m_{\ell\bar{\ell}}$ measurement for $m_{\ell\bar{\ell}}$ > 400 GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $p_{T,\ell\bar{\ell}}$ measurement for $p_{T,\ell\bar{\ell}}$ < 20 GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $p_{T,\ell\bar{\ell}}$ measurement for $p_{T,\ell\bar{\ell}}$ $\in$ [20, 70] GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Ranking of the systematic uncertainties with marginalisation for the $A_C^{\ell\ell}$ vs $p_{T,\ell\bar{\ell}}$ measurement for $p_{T,\ell\bar{\ell}}$ > 70 GeV. The effect on unfolded $A_C$ for down and up variation of the systematic uncertainty is shown, respectively. The pulls and constraints of the ranked NPs are obtained from data.

Post-marginalisation correlation coefficients $\rho_{ij}$ of nuisance parameters for the $A_C^{t\bar{t}}$ inclusive measurement. Only $|\rho_{ij}| > 0.05$ values are included.

Post-marginalisation correlation coefficients $\rho_{ij}$ of nuisance parameters for the $A_C^{t\bar{t}}$ vs $m_{t\bar{t}}$ measurement. Only $|\rho_{ij}| > 0.05$ values are included.

Post-marginalisation correlation coefficients $\rho_{ij}$ of nuisance parameters for the $A_C^{t\bar{t}}$ vs $p_{T,t\bar{t}}$ measurement. Only $|\rho_{ij}| > 0.05$ values are included.

Post-marginalisation correlation coefficients $\rho_{ij}$ of nuisance parameters for the $A_C^{t\bar{t}}$ vs $\beta_{z,t\bar{t}}$ measurement. Only $|\rho_{ij}| > 0.05$ values are included.

Post-marginalisation correlation coefficients $\rho_{ij}$ of nuisance parameters for the $A_C^{\ell\ell}$ inclusive measurement. Only $|\rho_{ij}| > 0.05$ values are included.

Post-marginalisation correlation coefficients $\rho_{ij}$ of nuisance parameters for the $A_C^{\ell\ell}$ vs $m_{\ell\bar{\ell}}$ measurement. Only $|\rho_{ij}| > 0.05$ values are included.

Post-marginalisation correlation coefficients $\rho_{ij}$ of nuisance parameters for the $A_C^{\ell\ell}$ vs $p_{T,\ell\bar{\ell}}$ measurement. Only $|\rho_{ij}| > 0.05$ values are included.

Post-marginalisation correlation coefficients $\rho_{ij}$ of nuisance parameters for the $A_C^{\ell\ell}$ vs $\beta_{z,\ell\bar{\ell}}$ measurement. Only $|\rho_{ij}| > 0.05$ values are included.

Covariance matrix for the $A_C^{t\bar{t}}$ vs $m_{t\bar{t}}$ measurement.

Covariance matrix for the $A_C^{t\bar{t}}$ vs $p_{T,t\bar{t}}$ measurement.

Covariance matrix for the $A_C^{t\bar{t}}$ vs $\beta_{z,t\bar{t}}$ measurement.

Covariance matrix for the $A_C^{\ell\ell}$ vs $m_{\ell\bar{\ell}}$ measurement.

Covariance matrix for the $A_C^{\ell\ell}$ vs $p_{T,\ell\bar{\ell}}$ measurement.

Covariance matrix for the $A_C^{\ell\ell}$ vs $\beta_{z,\ell\bar{\ell}}$ measurement.

Exponent according to the Eq. 5 as a function of transverse momenta extracted from p+Au/p+Al and $^{3}$HeAu/p+Au collision systems. The uncertainties from the $N_{coll}$ calculations and from the overall normalization of p+p cancel in these ratios

Nuclear modification factors in p+Al, p+Au, d+Au, and $^{3}$He+Au in five centrality bins and for inelastic collisions at $\sqrt{s}$ = 200 GeV. The right boxes are the uncertainties of the $N_{coll}$ collisions from the Glauber model, while the left box represents the overall normalization uncertainty from p+p collisions

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations. Centrality vs. the number of collisions is tabulated here.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations. Centrality vs. Average $R_{xA}$ is tabulated here.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations. Centrality vs. the number of collisions is tabulated here.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations. Centrality vs. Average $R_{xA}$ is tabulated here.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations. Centrality vs. the number of collisions per projectile participant is tabulated here.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations. Centrality vs. Average $R_{xA}$ is tabulated here.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations. Centrality vs. the number of collisions per projectile participant is tabulated here.

Average $R_{xA}$ in two different $p_T$ regions versus the number of collisions (panels a,b) and the number of collisions per projectile participant (panels c,d). The tilted error bars represent the anti-correlated uncertainty on the y and x-axis due to the Ncoll calculations. Centrality vs. Average $R_{xA}$ is tabulated here.

$R_{xA}$ for inelastic collisions compared to three different nuclear PDF calculations and their uncertainties. The data points include the statistical and systematical uncertainties. Relevant data has been tabulated in Figure 7.

The upper panel (a) shows the <$R_{xA}$> above $p_T$ =8 Gev/c as a function $N_{coll}$/$N_{proj}$. The data are compared to predictions from [49] for the consequences of high-$x$ nucleon size fluctuations. The lower panels show the $R_{xA}$ as a function of $p_T$ for most central on left (b) and most peripheral on right panel (c) collisions. Panel (a) has been tabulated in Figure 10 and Panel (b) and (c) in Figure 9.

The upper panel (a) shows the <$R_{xA}$> above $p_T$ =8 Gev/c as a function $N_{coll}$/$N_{proj}$. The data are compared to predictions from [49] for the consequences of high-$x$ nucleon size fluctuations. The lower panels show the $R_{xA}$ as a function of $p_T$ for most central on left (b) and most peripheral on right panel (c) collisions. Panel (a) has been tabulated in Figure 10 and Panel (b) and (c) in Figure 9. Centrality vs. the number of collisions per projectile participant in panel (a) is tabulated here.

The upper panel (a) shows the <$R_{xA}$> above $p_T$ =8 Gev/c as a function $N_{coll}$/$N_{proj}$. The data are compared to predictions from [49] for the consequences of high-$x$ nucleon size fluctuations. The lower panels show the $R_{xA}$ as a function of $p_T$ for most central on left (b) and most peripheral on right panel (c) collisions. Panel (a) has been tabulated in Figure 10 and Panel (b) and (c) in Figure 9. Centrality vs. Average $R_{xA}$ in panel (a) is tabulated here.

The upper panel (a) shows the <$R_{xA}$> above $p_T$ =8 Gev/c as a function $N_{coll}$/$N_{proj}$. The data are compared to predictions from [49] for the consequences of high-$x$ nucleon size fluctuations. The lower panels show the $R_{xA}$ as a function of $p_T$ for most central on left (b) and most peripheral on right panel (c) collisions. Panel (a) has been tabulated in Figure 10 and Panel (b) and (c) in Figure 9. Panel(b) and panel(c) are tabulated here.

Integrated yields for 1--2 GeV/c in panel (a) and 2--3 GeV/c in panel (b) as a function of charged particle multiplicity density at mid rapidity.

Integrated yields for 1--2 GeV/c in panel (a) and 2--3 GeV/c in panel (b) as a function of charged particle multiplicity density at mid rapidity. Centrality vs. charged particle multiplicity density at mid rapidity is tabulated here.

Integrated yields for 1--2 GeV/c in panel (a) and 2--3 GeV/c in panel (b) as a function of charged particle multiplicity density at mid rapidity. Centrality vs. integrated yield is tabulated here.

Integrated yields for 1--2 GeV/c in panel (a) and 2--3 GeV/c in panel (b) as a function of charged particle multiplicity density at mid rapidity.

Integrated yields for 1--2 GeV/c in panel (a) and 2--3 GeV/c in panel (b) as a function of charged particle multiplicity density at mid rapidity. Centrality vs. charged particle multiplicity density at mid rapidity is tabulated here.

Integrated yields for 1--2 GeV/c in panel (a) and 2--3 GeV/c in panel (b) as a function of charged particle multiplicity density at mid rapidity. Centrality vs. integrated yield is tabulated here.

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.

Direct photon spectra(Nucl. Part. Phys. 23, A1 (1997)) normalized by $(dN_{ch}/d\eta)^{1.25}$ for in p+p at $\sqrt{s_{NN}}$= 63 GeV.

Direct photon spectra(Sov. J. Nucl. Phys. 51, 836 (1990)) normalized by $(dN_{ch}/d\eta)^{1.25}$ for in p+p at $\sqrt{s_{NN}}$= 63 GeV.

Direct photon spectra normalized by $(dN_{ch}/d\eta)^{1.25}$ for in Au+Au at $\sqrt{s_{NN}}$= 62.4 GeV.

Direct photon spectra normalized by $(dN_{ch}/d\eta)^{1.25}$ for in Au+Au at $\sqrt{s_{NN}}$= 39.0 GeV.

Direct photon spectra(Physical Review C91, 064904 (2015)) normalized by $(dN_{ch}/d\eta)^{1.25}$ for in Au+Au at $\sqrt{s_{NN}}$= 200 GeV.

Direct photon spectra(Physical Review Letters 104, 132301 (2010)) normalized by $(dN_{ch}/d\eta)^{1.25}$ for in Au+Au at $\sqrt{s_{NN}}$= 200 GeV.

Direct photon spectra(Physical Review Letters 109, 152302 (2012)) normalized by $(dN_{ch}/d\eta)^{1.25}$ for in Au+Au at $\sqrt{s_{NN}}$= 200 GeV.

Direct photon spectra(Physical Review C98, 054902 (2018)) normalized by $(dN_{ch}/d\eta)^{1.25}$ for in Cu+Cu at $\sqrt{s_{NN}}$= 200 GeV.

Direct photon spectra(Physics Letters B754 235 (2016)) normalized by $(dN_{ch}/d\eta)^{1.25}$ for in Pb+Pb at $\sqrt{s_{NN}}$= 2760 GeV.

Number of binary collisions N_{coll} vs. Charged particle multiplicity in Au+Au at various cm energies.

Charged particle multiplicity vs. centrality in Au+Au at various cm energies.

Number of binary collisions N_{coll} vs. centrality in Au+Au at various cm energies.

Integrated direct photon yield vs. Charged particle multiplicity in various systems at various cm energies.

Charged particle multiplicity vs. centrality in various systems at various cm energies.

Integrated direct photon yield vs. centrality in various systems at various cm energies.

Integrated direct photon yield vs. Charged particle multiplicity in various systems at various cm energies.

Charged particle multiplicity vs. centrality in various systems at various cm energies.

Integrated direct photon yield vs. centrality in various systems at various cm energies.

Integrated per-trigger yields as a function of associated particle angle relative to $\psi_2$ event plane $\phi^a - \psi_2$ integrated over the near side $|\Delta\phi|<\pi/4$.

Integrated per-trigger yields as a function of associated particle angle relative to $\psi_2$ event plane $\phi^a - \psi_2$ integrated over the away side $|\Delta\phi|<\pi/4$.

Integrated per-trigger yields as a function of associated particle angle relative to $\psi_2$ event plane for near-side $|\Delta\phi|<\pi/4$ and for away-side $|\Delta\phi|<\pi/4$.

Integrated per-trigger yields as a function of associated particle angle relative to $\psi_2$ event plane for near-side $|\Delta\phi|<\pi/4$ and for away-side $|\Delta\phi|<\pi/4$.

Per-trigger yields $Y(\Delta\phi)$ of dihadrons pairs measured in Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200~{\rm GeV}$ after subtracting the underlying event model with several $p_T$ selections. Systematic uncertainties due to track matching are given.

Per-trigger yields $Y(\Delta\phi)$ of dihadrons pairs measured in Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200~{\rm GeV}$ after subtracting the underlying event model with several $p_T$ selections. Systematic uncertainties due to track matching are given.

$\Psi_{2}$-dependent $C(\Delta\phi)$ for $2<p_{T}^{trig}<4$ GeV/c and $1<p_{T}^{assoc}<2$ GeV/c.

$\Psi_{2}$-dependent $C(\Delta\phi)$ for $2<p_{T}^{trig}<4$ GeV/c and $2<p_{T}^{assoc}<4$ GeV/c.

$\Psi_{2}$-dependent $C(\Delta\phi)$ for $4<p_{T}^{trig}<10$ GeV/c and $2<p_{T}^{assoc}<4$ GeV/c.

$\Psi_{3}$-dependent $C(\Delta\phi)$ for $2<p_{T}^{trig}<4$ GeV/c and $1<p_{T}^{assoc}<2$ GeV/c.

$\Psi_{3}$-dependent $C(\Delta\phi)$ for $2<p_{T}^{trig}<4$ GeV/c and $2<p_{T}^{assoc}<4$ GeV/c.

$\Psi_{3}$-dependent $C(\Delta\phi)$ for $4<p_{T}^{trig}<10$ GeV/c and $1<p_{T}^{assoc}<2$ GeV/c.

$\Psi_{2}$-dependent $Y(\Delta\phi)$ for $2<p_{T}^{trig}<4$ GeV/c and $1<p_{T}^{assoc}<2$ GeV/c.

$\Psi_{2}$-dependent $C(\Delta\phi)$ for $2<p_{T}^{trig}<4$ GeV/c and $2<p_{T}^{assoc}<4$ GeV/c.

$\Psi_{2}$-dependent $C(\Delta\phi)$ for $4<p_{T}^{trig}<10$ GeV/c and $2<p_{T}^{assoc}<4$ GeV/c.

$\Psi_{3}$-dependent $C(\Delta\phi)$ for $2<p_{T}^{trig}<4$ GeV/c and $1<p_{T}^{assoc}<2$ GeV/c.

$\Psi_{3}$-dependent $C(\Delta\phi)$ for $2<p_{T}^{trig}<4$ GeV/c and $2<p_{T}^{assoc}<4$ GeV/c.

$\Psi_{3}$-dependent $C(\Delta\phi)$ for $4<p_{T}^{trig}<10$ GeV/c and $2<p_{T}^{assoc}<4$ GeV/c.

$\Psi_2$-dependent C for 2 $< p_T^t <$ 4 and 1 $< p_T^a <$ 2 GeV/c.

$\Psi_2$-dependent C for 2 $< p_T^t <$ 4 and 2 $< p_T^a <$ 4 GeV/c.

$\Psi_2$-dependent C for 4 $< p_T^t <$ 10 and 2 $< p_T^a <$ 4 GeV/c.

$\Psi_2$-dependent Y for 2 $< p_T^t <$ 4 and 1 $< p_T^a <$ 2 GeV/c.

$\Psi_2$-dependent Y for 2 $< p_T^t <$ 4 and 2 $< p_T^a <$ 4 GeV/c.

$\Psi_2$-dependent Y for 4 $< p_T^t <$ 10 and 2 $< p_T^a <$ 4 GeV/c.

$\Psi_3$-dependent C for 2 $< p_T^t <$ 4 and 1 $< p_T^a <$ 2 GeV/c.

$\Psi_3$-dependent C for 2 $< p_T^t <$ 4 and 2 $< p_T^a <$ 4 GeV/c.

$\Psi_3$-dependent C for 4 $< p_T^t <$ 10 and 2 $< p_T^a <$ 4 GeV/c.

$\Psi_3$-dependent Y for 2 $< p_T^t <$ 4 and 1 $< p_T^a <$ 2 GeV/c.

$\Psi_3$-dependent Y for 2 $< p_T^t <$ 4 and 2 $< p_T^a <$ 4 GeV/c.

$\Psi_3$-dependent Y for 4 $< p_T^t <$ 10 and 2 $< p_T^a <$ 4 GeV/c.

Fourier fit coefficients for CNT-MPCN ($d$-going) correlations, as a function of collision system and $\pi^0$ $p_T$: (a) the negative of the dipole coefficient, $-c_1$; (b) the quadrupole coefficient $c_2$; (c) the ratio ${-c_2}/{c_1}$.

Fourier fit coefficients for CNT-MPCN ($d$-going) correlations, as a function of collision system and $\pi^0$ $p_T$: Fractional systematic uncertainty on the quadrupole coefficient $c_2$ for $d$+Au.

Fourier fit coefficients for CNT-MPCN ($d$-going) correlations, as a function of collision system and $\pi^0$ $p_T$: Fractional systematic uncertainty on the quadrupole coefficient $c_2$ for $d$+Au and $p$+$p$.

System centrality dependence of correlation coefficients $-c_1$, $c_2$, and the ratio ${-c_2}/{c_1}$ as a function of $N_{part}$.

System centrality dependence of correlation coefficients $-c_1$, $c_2$, and the ratio ${-c_2}/{c_1}$ as a function of $N_{part}$.

Summary of $c\bar{c}$ and $b\bar{b}$ cross sections measured in $d$+Au collisions expressed as nucleon-nucleon equivalent cross section by dividing the $d$+Au cross section by the average number of binary nucleon-nucleon collisions $N_{coll}$ = 7.6 $\pm$ 0.4.