Measured unfolded differential cross-section as a function of the dilepton transverse momentum $p^{ll}_{\mathrm{T}}$. NLO predictions from Sherpa 2.2.10 and MadGraph5_aMC@NLO 2.7.3 are also shown. The uncertainty in the predictions is divided into statistical and theoretical uncertainties (scale and PDF+$\alpha_{s}$).

Measured unfolded differential cross-section as a function of the the four-body transverse momentum $p^{ll\gamma\gamma}_{\mathrm{T}}$. NLO predictions from Sherpa 2.2.10 and MadGraph5_aMC@NLO 2.7.3 are also shown. The uncertainty in the predictions is divided into statistical and theoretical uncertainties (scale and PDF+$\alpha_{s}$).

Measured unfolded differential cross-section as a function of the diphoton invariant mass $m_{\gamma\gamma}$. NLO predictions from Sherpa 2.2.10 and MadGraph5_aMC@NLO 2.7.3 are also shown. The uncertainty in the predictions is divided into statistical and theoretical uncertainties (scale and PDF+$\alpha_{s}$).

Measured unfolded differential cross-section as a function of the four-body invariant mass $m_{ll\gamma\gamma}$. NLO predictions from Sherpa 2.2.10 and MadGraph5_aMC@NLO 2.7.3 are also shown. The uncertainty in the predictions is divided into statistical and theoretical uncertainties (scale and PDF+$\alpha_{s}$).

Expected and observed $95\%$ confidence intervals for the coupling parameters $f_{T,j}/\Lambda^{4}$ of transverse dimension-8 operators. All parameter values outside of the stated range are excluded at the chosen confidence level. No unitarity constraints are applied.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,8}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,0}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,1}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,2}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,5}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,6}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,7}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

Expected and observed unitarised $95\%$ confidence intervals for the coupling parameter $f_{T,9}/\Lambda^{4}$ in the clipping energy range between 1.1 and 5 TeV. The non-unitarised limits ($E_c = \infty$) are also shown. All parameter values outside of the stated range are excluded at the chosen confidence level.

$\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

Yield ratio of non-prompt over prompt $\mathrm{D}_{s}^{+}$ for 30-50% Pb--Pb collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02$ TeV.

Yield ratio of non-prompt $\mathrm{D}_{s}^{+}$ over non-prompt $\mathrm{D}^{0}$ for pp collisions at $\sqrt{s} = 5.02$ TeV.

Yield ratio of non-prompt $\mathrm{D}_{s}^{+}$ over non-prompt $\mathrm{D}^{0}$ for 0-10% Pb--Pb collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02$ TeV.

Yield ratio of non-prompt $\mathrm{D}_{s}^{+}$ over non-prompt $\mathrm{D}^{0}$ for 30-50% Pb--Pb collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02$ TeV.

Non-prompt $\mathrm{D}_{s}^{+}$ Raa in 0-10% collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02$ TeV.

Non-prompt $\mathrm{D}_{s}^{+}$ Raa in 30-50% collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02$ TeV.

Non-prompt $\mathrm{D}_{s}^{+}$ over non-prompt $\mathrm{D}^{0}$ Raa in 0-10% collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02$ TeV.

Non-prompt $\mathrm{D}_{s}^{+}$ over non-prompt $\mathrm{D}^{0}$ Raa in 30-50% collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02$ TeV.

Non-prompt over prompt $\mathrm{D}_{s}^{+}$ Raa in 0-10% collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02$ TeV.

Non-prompt over prompt $\mathrm{D}_{s}^{+}$ Raa in 0-10% collisions at $\sqrt{s_{\mathrm{NN}}} = 5.02$ TeV.

$p_{\mathrm T}$-differential yield of $\frac{\mathrm{K^{*0}} + \overline{\mathrm{K^{*0}}}}{2}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV (Multiplicity class 10-20%).

$p_{\mathrm T}$-differential yield of $\frac{\mathrm{K^{*0}} + \overline{\mathrm{K^{*0}}}}{2}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV (Multiplicity class 20-40%).

$p_{\mathrm T}$-differential yield of $\frac{\mathrm{K^{*0}} + \overline{\mathrm{K^{*0}}}}{2}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV (Multiplicity class 40-60%).

$p_{\mathrm T}$-differential yield of $\frac{\mathrm{K^{*0}} + \overline{\mathrm{K^{*0}}}}{2}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV (Multiplicity class 60-80%).

$p_{\mathrm T}$-differential yield of $\frac{\mathrm{K^{*0}} + \overline{\mathrm{K^{*0}}}}{2}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV (Multiplicity class 80-100%).

$p_{\mathrm T}$-differential yield of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV (NSD).

$p_{\mathrm T}$-differential yield of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV (Multiplicity class 0-5%).

$p_{\mathrm T}$-differential yield of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV (Multiplicity class 5-10%).

$p_{\mathrm T}$-differential yield of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV (Multiplicity class 10-20%).

$p_{\mathrm T}$-differential yield of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV (Multiplicity class 20-40%).

$p_{\mathrm T}$-differential yield of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV (Multiplicity class 40-60%).

$p_{\mathrm T}$-differential yield of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV (Multiplicity class 60-80%).

$p_{\mathrm T}$-differential yield of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV (Multiplicity class 80-100%).

$p_{\mathrm T}$-differential $R_{\mathrm{pPb}}$ of K$^{*0}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV.

$p_{\mathrm T}$-differential $R_{\mathrm{pPb}}$ of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$8.16 TeV.

$p_{\mathrm T}$-differential $R_{\mathrm{pPb}}$ of K$^{*0}$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV.

$p_{\mathrm T}$-differential $R_{\mathrm{pPb}}$ of $\phi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV.

$p_{\mathrm T}$-differential $R_{\mathrm{pPb}}$ of $\Xi$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV.

$p_{\mathrm T}$-differential $R_{\mathrm{pPb}}$ of $\Omega$ in p-Pb collisions at $\sqrt{s_{\mathrm{NN}}}~=~$5.02 TeV.

Differential cross sections for the EW W$\gamma$jj production. Given that the ranges of some variables extend to infinity, the last bins accommodate all the events up to infinity as marked by the bin label.

Differential cross sections for the EW W$\gamma$jj production. Given that the ranges of some variables extend to infinity, the last bins accommodate all the events up to infinity as marked by the bin label.

Differential cross sections for the EW W$\gamma$jj production. Given that the ranges of some variables extend to infinity, the last bins accommodate all the events up to infinity as marked by the bin label.

Differential cross sections for the EW W$\gamma$jj production. Given that the ranges of some variables extend to infinity, the last bins accommodate all the events up to infinity as marked by the bin label.

Differential cross sections for the EW W$\gamma$jj production. Given that the ranges of some variables extend to infinity, the last bins accommodate all the events up to infinity as marked by the bin label.

Differential cross sections for the EW + QCD W$\gamma$jj production. Given that the ranges of some variables extend to infinity, the last bins accommodate all the events up to infinity as marked by the bin label.

Differential cross sections for the EW +QCD W$\gamma$jj production. Given that the ranges of some variables extend to infinity, the last bins accommodate all the events up to infinity as marked by the bin label.

Differential cross sections for the EW + QCD W$\gamma$jj production. Given that the ranges of some variables extend to infinity, the last bins accommodate all the events up to infinity as marked by the bin label.

Differential cross sections for the EW + QCD W$\gamma$jj production. Given that the ranges of some variables extend to infinity, the last bins accommodate all the events up to infinity as marked by the bin label.

Differential cross sections for the EW + QCD W$\gamma$jj production. Given that the ranges of some variables extend to infinity, the last bins accommodate all the events up to infinity as marked by the bin label.

Differential cross sections for the EW + QCD W$\gamma$jj production. Given that the ranges of some variables extend to infinity, the last bins accommodate all the events up to infinity as marked by the bin label.

The $m_{W\gamma}$ distribution for muon events (electron events are also used and the muon events are shown here as an example.) satisfying the aQGC region selection and used to set constraints on the anomalous gauge coupling parameters (left). The gray line represents a nonzero $f_{M,2}/Λ^4$ setting. Events with $m_{W\gamma} > 1500 GeV$ are included in the last bin.

Exclusion limits at the 95% CL for each aQGC coefficient, derived from the mW$\gamma$ distribution, assuming all other coefficients are set to zero. Unitarity bounds corresponding to each operator are also listed. All coupling parameter limits are in TeV^{−4}, while Ubound values are in TeV.

Results summary.

Measured and predicted fiducial cross-sections in both the lllljj and ll$\nu\nu$jj channels for the inclusive ZZjj processes. Uncertainties due to different sources are presented.

Measured and predicted fiducial cross-sections in both the lllljj and ll$\nu\nu$jj channels for the inclusive ZZjj processes. Uncertainties due to different sources are presented.

Observed and expected multivariate discriminant distribution in the $\ell\ell\ell\ell jj$ QCD CR.

Observed and expected multivariate discriminant distribution in the $\ell\ell\ell\ell jj$ QCD CR.

Observed and expected multivariate discriminant distribution in the $\ell\ell\ell\ell jj$ SR.

Observed and expected multivariate discriminant distribution in the $\ell\ell\ell\ell jj$ SR.

Observed and expected multivariate discriminant distribution in the $\ell\ell\nu\nu jj$ SR.

Observed and expected multivariate discriminant distribution in the $\ell\ell\nu\nu jj$ SR.

Charged hadron $dN_{ch}/d\eta$ as a function of pseudorapidity in high-multiplicity 0%-5% central $p$+Al collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Charged hadron $dN_{ch}/d\eta$ as a function of pseudorapidity in various mutiplicity classes of $^3$He+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Charged hadron $dN_{ch}/d\eta$ as a function of pseudorapidity in various mutiplicity classes of $d$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Charged hadron $dN_{ch}/d\eta$ as a function of pseudorapidity in various mutiplicity classes of $p$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Charged hadron $dN_{ch}/d\eta$ as a function of pseudorapidity in various mutiplicity classes of $p$+Al collisions at $\sqrt{s_{NN}}$ = 200 GeV.

Midrapidity charged hadron $dN_{ch}/d\eta$ per participating quark pair ($N_{qp}$/2) as a function of the number of participating quarks ($N_{qp}$).

Midrapidity charged hadron $dN_{ch}/d\eta$ per participating quark pair ($N_{qp}$/2) as a function of the number of participating quarks ($N_{qp}$).

Midrapidity charged hadron $dN_{ch}/d\eta$ per participating quark pair ($N_{qp}$/2) as a function of the number of participating quarks ($N_{qp}$).

Midrapidity charged hadron $dN_{ch}/d\eta$ per participating quark pair ($N_{qp}$/2) as a function of the number of participating quarks ($N_{qp}$).

Midrapidity charged hadron $dN_{ch}/d\eta$ per participating quark pair ($N_{qp}$/2) as a function of the number of participating quarks ($N_{qp}$).

Elliptic flow $v_2$ as a function of pseudorapidity in high-multiplicity 0%-5% central $p$+Al, $p$+Au, $d$+Au, and $^3$He+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV.