Fully corrected assorted $\pi^{0}$-hadron conditional pair distributions for $d$+Au collisions centrality 0-88% and $p$+$p$ collisions. The trigger $\pi^0$s are within 5 < $p_{T,trig}$ < 10 GeV/$c$ and are correlated with hadrons with $p_{T,assoc}$ of 1.5-2 GeV/$c$, 2-3GeV/$c$, 3-4 GeV/$c$, and 4-5 GeV/$c$.

Near- and far-side widths as a function of $p_{T, assoc}$ for charged hadron azimuthal correlations from minimumbias $d$+Au collisions.

Near-side width, far-side width as a function of $p_{T,assoc}$ for a charged pion and neutral pion triggers from the $p_{T,trig}$ range of 5-10 GeV/$c$ in minimum bias $d$+Au collisions.

Near-side width, far-side width as a function of $p_{T,assoc}$ for a charged pion and neutral pion triggers from the $p_{T,trig}$ range of 5-10 GeV/$c$ in minimum bias $d$+Au collisions.

Near- and far-side conditional yields as a function of $p_{T, assoc}$ for charged hadron triggers (2.5−4 GeV/$c$) and associated charged hadrons from $d$+Au collisions.

Near- and far-side conditional yields as a function of $p_{T, assoc}$ for charged hadron triggers (4−6 GeV/$c$) and associated charged hadrons from $d$+Au collisions.

Near- and far-side conditional yields as a function of $p_{T, assoc}$ for charged pion triggers (5−10 GeV/$c$) and associated charged hadrons from minimum-bias $d$+Au collisions.

Near- and far-side conditional yields as a function of $p_{T, assoc}$ for neutral pion triggers (5−10 GeV/$c$) and associated charged hadrons from minimum-bias $d$+Au collisions.

Near- and far-side widths and conditional yields as a function of $N_{coll}$ for charged hadron triggers (2.5−4 GeV/$c$) and associated charged hadrons (1–2.5 GeV/$c$) from $d$+Au collisions.

Near- and far-side widths and conditional yields as a function of centrality for neutral pion triggers (5−10 GeV/$c$) and associated charged hadrons (2–3 GeV/$c$) from $d$+Au collisions.

The extracted $\sqrt{<j^2_T>}$ for minimum-bias $d$+Au collisions from all four dihadron correlations. The trigger $p_T$ ranges are 3 - 5 GeV/$c$, 5 - 10 GeV/$c$ and 5 - 10 GeV/$c$ for $h^{\pm}$ - $h^{\pm}$, $\pi^0$ - $h^{\pm}$ and $\pi^{\pm}$ - $h^{\pm}$ assorted-$p_T$ correlations.

The extracted $\sqrt{<j^2_T>}$ for minimum-bias $d$+Au collisions from all four dihadron correlations. The trigger $p_T$ ranges are 3 - 5 GeV/$c$, 5 - 10 GeV/$c$ and 5 - 10 GeV/$c$ for $h^{\pm}$ - $h^{\pm}$, $\pi^0$ - $h^{\pm}$ and $\pi^{\pm}$ - $h^{\pm}$ assorted-$p_T$ correlations.

The extracted $\sqrt{<j^2_T>}$ for minimum-bias $d$+Au collisions from all four dihadron correlations. The trigger $p_T$ ranges are 3 - 5 GeV/$c$, 5 - 10 GeV/$c$ and 5 - 10 GeV/$c$ for $h^{\pm}$ - $h^{\pm}$, $\pi^0$ - $h^{\pm}$ and $\pi^{\pm}$ - $h^{\pm}$ assorted-$p_T$ correlations.

The extracted $\sqrt{<j^2_T>}$ for minimum-bias $d$+Au collisions from all four dihadron correlations. The trigger $p_T$ ranges are 3 - 5 GeV/$c$, 5 - 10 GeV/$c$ and 5 - 10 GeV/$c$ for $h^{\pm}$ - $h^{\pm}$, $\pi^0$ - $h^{\pm}$ and $\pi^{\pm}$ - $h^{\pm}$ assorted-$p_T$ correlations.

$<sin^2(\phi_{jj})>$ for minimum bias $d$+Au collisions as function of associated particle $p_T$ for $\pi^{\pm}$ - $h^{\pm}$ and $\pi^0$ - $h^{\pm}$ correlations where the trigger particles have a $p_T$ between 5 - 10 GeV/$c$.

$<sin^2(\phi_{jj})>$ for minimum bias $d$+Au collisions as function of associated particle $p_T$ for $\pi^{\pm}$ - $h^{\pm}$ and $\pi^0$ - $h^{\pm}$ correlations where the trigger particles have a $p_T$ between 5 - 10 GeV/$c$.

$<sin^2(\phi_{jj})>$ for minimum bias $d$+Au collisions as function of trigger particle $p_T$ for $\pi^{\pm}$ - $h^{\pm}$ and $\pi^0$ - $h^{\pm}$ correlations.

$<sin^2(\phi_{jj})>$ for minimum bias $d$+Au collisions as function of trigger particle $p_T$ for $\pi^{\pm}$ - $h^{\pm}$ and $\pi^0$ - $h^{\pm}$ correlations.

Near-side and far-side $p_{out}$ distributions for minimum bias $d$+Au collisions obtained from $\pi^{\pm}$ - $h^{\pm}$ correlation. The trigger is 5 < $p_{T,trig}$ < 10 GeV/$c$, the associated particle is 0.5 < $p_{T,assoc}$ < 5.0 GeV/$c$.

Conditional yield as a function of $x_E$ for near-side and far-side for $\pi^{\pm}$ - $h^{\pm}$ and $\pi^0$ - $h^{\pm}$ from minimum-bias $d$+Au collisions.

Conditional yield as a function of $x_E$ for near-side and far-side for $\pi^{\pm}$ - $h^{\pm}$ and $\pi^0$ - $h^{\pm}$ from minimum-bias $d$+Au collisions.

Conditional yield as a function of $x_E$ for near-side and far-side correlations for $\pi^{\pm}$ - $h^{\pm}$ correlations from minimum-bias $d$+Au collisions. The trigger pions are 5 < $p_{T,trig}$ < 6 GeV/$c$.

Conditional yield as a function of $x_E$ for near-side and far-side correlation for $\pi^{\pm}$ - $h^{\pm}$ correlation for several different trigger $p_T$s for minimum-bias $d$+Au collisions.

Far-side conditional yield as a function of $p_{T,trig}$ for different ranges of $x_E$ for $\pi^{\pm}$ - $h^{\pm}$ correlation for minimum-bias $d$+Au collisions.

Far-side conditional yield as a function of $p_{T,trig}$ for different ranges of $x_E$ for $\pi^{\pm}$ - $h^{\pm}$ correlation for minimum-bias $d$+Au collisions.

Far-side conditional yield as a function of $p_{T,trig}$ for different ranges of $x_E$ for $\pi^{\pm}$ - $h^{\pm}$ correlation for minimum-bias $d$+Au collisions.

Far-side conditional yield as a function of $p_{T,trig}$ for different ranges of $x_E$ for $\pi^{\pm}$ - $h^{\pm}$ correlation for minimum-bias $d$+Au collisions.

The comparison of the $\sqrt{<j^2_T>}$ values between $d$+Au and $p$+$p$ for the $\pi^{\pm}$ - $h^{\pm}$ correlations and $\pi^0$ - $h^{\pm}$ correlations. The trigger pion range is 5 - 10 GeV/$c$.

The comparison of the $\sqrt{<j^2_T>}$ values between $d$+Au and $p$+$p$ for the $\pi^{\pm}$ - $h^{\pm}$ correlations and $\pi^0$ - $h^{\pm}$ correlations. The trigger pion range is 5 - 10 GeV/$c$.

Quadrature difference between minimum-bias $d$+Au and $p$+$p$ $<sin^2(\phi_{jj})>$ values.

Quadrature difference between minimum-bias $d$+Au and $p$+$p$ $<sin^2(\phi_{jj})>$ values.

The comparison of the $p_{out}$ distribution at the near-side and far-side between central $d$+Au collisions and $p$+$p$ collisions. Results are obtained for $\pi^{\pm}$ - $h^{\pm}$ correlations with the associated hadron range 0.5 < $p_{T,assoc}$ < 5 GeV/$c$ and trigger pion range of 5 < $p_{T,trig}$ < 10 GeV/$c$.

The comparison of the $x_E$ distribution from $\pi^{\pm}$ - $h^{\pm}$ correlation at the near-side and far-side between minimum bias $d$+Au collisions and $p$+$p$ collisions. The trigger $\pi^{\pm}$ are from 5 - 10 GeV/$c$.

Far-side conditional yield as a function of $p_{T,trig}$ for different ranges of $x_E$ from $p$+$p$ collisions, triggers are $\pi^{\pm}$ from 5 - 10 GeV/$c$.

Far-side conditional yield as a function of $p_{T,trig}$ for different ranges of $x_E$ from $p$+$p$ collisions, triggers are $\pi^{\pm}$ from 5 - 10 GeV/$c$.

Far-side conditional yield as a function of $p_{T,trig}$ for different ranges of $x_E$ from $p$+$p$ collisions, triggers are $\pi^{\pm}$ from 5 - 10 GeV/$c$.

The fitted fractional change in $dN/d_{x_E}$ per unit $p_{T,trig}$ of the far-side conditional yield for different ranges of $x_E$.

Centrality dependent near and far $CY(p_T)$ for $\pi^{\pm}$ - $h^{\pm}$ correlations and $\pi^0$ - $h^{\pm}$ correlations.

Centrality dependent near and far $CY(p_T)$ for $\pi^{\pm}$ - $h^{\pm}$ correlations and $\pi^0$ - $h^{\pm}$ correlations.

Centrality dependence to the ratio of near and far $CY(p_T)$ for $d$+Au to $p$+$p$.

Centrality dependence to the ratio of near and far $CY(p_T)$ for $d$+Au to $p$+$p$.

Near- and far-side widths as a function of $p_{T, trig}$ for charged pion triggers and associated charged hadrons (2–4.5 Gev/$c$) from minimum-bias $d$+Au collisions.

Near- and far-side widths as a function of $p_{T, trig}$ for neutral pion triggers and associated charged hadrons (2–4.5 Gev/$c$) from minimum-bias $d$+Au collisions.

Corrected $m(HH)$ distribution in the resolved $4b$ signal region (dots), after the fit under the background-only hypothesis. The error bars on the $4b$ points represent the Poisson uncertainties corresponding to their event yields. The background model (teal histogram) is shown with its total post-fit uncertainty (gray band). The final bin includes overflow. Representative spin-0 signal hypotheses (dashed, dotted, and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. The bottom panel shows the difference between the $4b$ distribution and the background model, relative to the background model. No significant excess of data relative to the SM background is observed.

Corrected $m(HH)$ distribution in the resolved $4b$ signal region (dots), after the fit under the background-only hypothesis. The error bars on the $4b$ points represent the Poisson uncertainties corresponding to their event yields. The background model (teal histogram) is shown with its total post-fit uncertainty (gray band). The final bin includes overflow. Representative spin-2 signal hypotheses (dashed, dotted, and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. The bottom panel shows the difference between the $4b$ distribution and the background model, relative to the background model. No significant excess of data relative to the SM background is observed.

Cumulative signal acceptance times efficiency as a function of the resonance mass for various selection steps in the boosted channel. The steps up to the $b$-tag categorization are shown for the spin-0.

Cumulative signal acceptance times efficiency as a function of the resonance mass for various selection steps in the boosted channel. The steps up to the $b$-tag categorization are shown for the spin-2.

Cumulative signal acceptance times efficiency as a function of the resonance mass for various selection steps in the boosted channel. The efficiencies of the three b-tag categories are shown for the spin-0 scenario; this efficiency is obtained after the other selection steps including the SR definition. The signal efficiency in the 4b region has a maximum around 1.5 TeV. Above that value the track-jets start to merge together, and for the highest resonance masses the 2b category becomes the most efficient.

Cumulative signal acceptance times efficiency as a function of the resonance mass for various selection steps in the boosted channel. The efficiencies of the three b-tag categories are shown for the spin-2 scenario; this efficiency is obtained after the other selection steps including the SR definition. The signal efficiency in the 4b region has a maximum around 1.5 TeV. Above that value the track-jets start to merge together, and for the highest resonance masses the 2b category becomes the most efficient.

Comparison of the background model (stacked histograms) with data (dots) in the $2b$ validation region. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background uncertainty (gray band) is computed by adding all individual components in quadrature and is not allowed to extend below zero.

Comparison of the background model (stacked histograms) with data (dots) in the $3b$ validation region. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background uncertainty (gray band) is computed by adding all individual components in quadrature and is not allowed to extend below zero.

Comparison of the background model (stacked histograms) with data (dots) in the $4b$ validation region. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background uncertainty (gray band) is computed by adding all individual components in quadrature and is not allowed to extend below zero.

The $m(HH)$ distributions in the boosted $2b$ signal regions (dots), after the fit under the background-only hypothesis. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background model (stacked histogram) is shown with its total post-fit uncertainty (gray band). The uncertainty bands are defined using an ensemble of curves constructed by sampling a multivariate Gaussian probability density function built from the covariance matrix of the fit. Representative spin-0 signal hypotheses (dashed and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. No significant excess of data relative to the SM background is observed.

The $m(HH)$ distributions in the boosted $2b$ signal regions (dots), after the fit under the background-only hypothesis. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background model (stacked histogram) is shown with its total post-fit uncertainty (gray band). The uncertainty bands are defined using an ensemble of curves constructed by sampling a multivariate Gaussian probability density function built from the covariance matrix of the fit. Representative spin-2 signal hypotheses (dashed and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. No significant excess of data relative to the SM background is observed.

The $m(HH)$ distributions in the boosted $3b$ signal regions (dots), after the fit under the background-only hypothesis. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background model (stacked histogram) is shown with its total post-fit uncertainty (gray band). The uncertainty bands are defined using an ensemble of curves constructed by sampling a multivariate Gaussian probability density function built from the covariance matrix of the fit. Representative spin-0 signal hypotheses (dashed and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. No significant excess of data relative to the SM background is observed.

The $m(HH)$ distributions in the boosted $3b$ signal regions (dots), after the fit under the background-only hypothesis. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background model (stacked histogram) is shown with its total post-fit uncertainty (gray band). The uncertainty bands are defined using an ensemble of curves constructed by sampling a multivariate Gaussian probability density function built from the covariance matrix of the fit. Representative spin-2 signal hypotheses (dashed and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. No significant excess of data relative to the SM background is observed.

The $m(HH)$ distributions in the boosted $4b$ signal regions (dots), after the fit under the background-only hypothesis. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background model (stacked histogram) is shown with its total post-fit uncertainty (gray band). The uncertainty bands are defined using an ensemble of curves constructed by sampling a multivariate Gaussian probability density function built from the covariance matrix of the fit. Representative spin-0 signal hypotheses (dashed and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. No significant excess of data relative to the SM background is observed.

The $m(HH)$ distributions in the boosted $4b$ signal regions (dots), after the fit under the background-only hypothesis. The error bars on the data points represent the Poisson uncertainties corresponding to their event yields. The background model (stacked histogram) is shown with its total post-fit uncertainty (gray band). The uncertainty bands are defined using an ensemble of curves constructed by sampling a multivariate Gaussian probability density function built from the covariance matrix of the fit. Representative spin-2 signal hypotheses (dashed and dashed-dotted lines) are overlaid, normalized to the overall expected limits on their cross-sections. No significant excess of data relative to the SM background is observed.

Expected (dashed black lines) and observed (solid black lines) 95% CL upper limits on the cross-section of resonant $HH$ production in the spin-0 signal models. The $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty ranges for the expected limits (colored bands) are shown. Expected limits using each of the resolved and boosted channels individually (dashed colored lines) are shown. The nominal $H\rightarrow b\bar{b}$ branching ratio is taken as 0.582.

Expected (dashed black lines) and observed (solid black lines) 95% CL upper limits on the cross-section of resonant $HH$ production in the spin-2 signal models. The $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty ranges for the expected limits (colored bands) are shown. Expected limits using each of the resolved and boosted channels individually (dashed colored lines) are shown. The theoretical prediction for the bulk RS model with $k/\bar{M}_{\text{Pl}} = 1$ (solid red line) is shown; the decrease below 350 GeV is due to a sharp reduction in the $G^{*}_{\text{KK}} \rightarrow HH$ branching ratio. The nominal $H\rightarrow b\bar{b}$ branching ratio is taken as 0.582.

Two-dimensional distribution of the invariant mass $m_{DV}$ and the track multiplicity in the High-pT jet SR for observed data events

Two-dimensional distribution of the invariant mass $m_{DV}$ and the track multiplicity in the High-pT jet SR for expected signal events in the strong gluino pair pair production model with m(gluino)=1.8 TeV, m(chi0)=0.2 TeV, tau(chi0)=0.1 ns

Two-dimensional distribution of the invariant mass $m_{DV}$ and the track multiplicity in the Trackless jet SR for observed data events

Two-dimensional distribution of the invariant mass $m_{DV}$ and the track multiplicity in the Trackless jet SR for expected signal events in the electroweak pair production model

Expected exclusion limits at 95% CL on the lifetime and mass of the neutralino in electroweakino pair production models

Expected (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in electroweakino pair production models

Expected (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in electroweakino pair production models

Observed exclusion limits at 95% CL on the lifetime and mass of the neutralino in electroweakino pair production models

Observed (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in electroweakino pair production models

Observed (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in electroweakino pair production models

Expected exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.4 TeV

Expected (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.4 TeV

Expected (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.4 TeV

Observed exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.4 TeV

Observed (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.4 TeV

Observed (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.4 TeV

Exclusion limits at 95% CL on the production cross section in the electroweak pair production model.

Exclusion limits at 95% CL on the production cross section in the strong gluino pair production models and m(gluino)=2.4 TeV

Expected exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.0 TeV

Expected (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.0 TeV

Expected (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.0 TeV

Observed exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.0 TeV

Observed (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.0 TeV

Observed (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.0 TeV

Expected exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.2 TeV

Expected (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.2 TeV

Expected (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.2 TeV

Observed exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.2 TeV

Observed (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.2 TeV

Observed (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.2 TeV

Expected exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=50 GeV

Expected (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=50 GeV

Expected (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=50 GeV

Observed exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=50 GeV

Observed (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=50 GeV

Observed (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=50 GeV

Expected exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=450 GeV

Expected (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=450 GeV

Expected (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=450 GeV

Observed exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=450 GeV

Observed (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=450 GeV

Observed (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=450 GeV

Expected exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.01 ns

Expected (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.01 ns

Expected (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.01 ns

Observed exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.01 ns

Observed (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.01 ns

Observed (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.01 ns

Expected exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.1 ns

Expected (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.1 ns

Expected (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.1 ns

Observed exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.1 ns

Observed (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.1 ns

Observed (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.1 ns

Expected exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=1 ns

Expected (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=1 ns

Expected (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=1 ns

Observed exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=1 ns

Observed (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=1 ns

Observed (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=1 ns

Expected exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=10 ns

Expected (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=10 ns

Expected (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=10 ns

Observed exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=10 ns

Observed (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=10 ns

Observed (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=10 ns

Exclusion limits at 95% CL on the production cross section in the strong gluino pair production models and m($ ilde{\chi}^0_1$)=1.25 TeV

Acceptance cutflow for the High-pT SR for representative points in the strong gluino pair production model. See additional resources for more information.

Acceptance cutflow for the Trackless SR for representative points in the electroweak pair production model. See additional resources for more information.

Acceptance cutflow for the Trackless SR for representative points in the electroweak pair production model with heavy-flavor quarks final state. See additional resources for more information.

Acceptance cutflow for the High-pT SR for representative points in the electroweak pair production model with heavy-flavor quarks final state. See additional resources for more information.

Reinterpretation Material: Event-level Efficiency for HighPt SR selections, R < 1150 mm

Reinterpretation Material: Event-level Efficiency for HighPt SR selections, R [1150, 3870] mm

Reinterpretation Material: Event-level Efficiency for HighPt SR selections, R > 3870 mm

Reinterpretation Material: Event-level Efficiency for Trackless SR selections, R < 1150 mm

Reinterpretation Material: Event-level Efficiency for Trackless SR selections, R [1150, 3870] mm

Reinterpretation Material: Event-level Efficiency for Trackless SR selections, R > 3870 mm

Reinterpretation Material: Vertex-level Efficiency for R < 22 mm

Reinterpretation Material: Vertex-level Efficiency for R [22, 25] mm

Reinterpretation Material: Vertex-level Efficiency for R [25, 29] mm

Reinterpretation Material: Vertex-level Efficiency for R [29, 38] mm

Reinterpretation Material: Vertex-level Efficiency for R [38, 46] mm

Reinterpretation Material: Vertex-level Efficiency for R [46, 73] mm

Reinterpretation Material: Vertex-level Efficiency for R [73, 84] mm

Reinterpretation Material: Vertex-level Efficiency for R [84, 111] mm

Reinterpretation Material: Vertex-level Efficiency for R [111, 120] mm

Reinterpretation Material: Vertex-level Efficiency for R [120, 145] mm

Reinterpretation Material: Vertex-level Efficiency for R [145, 180] mm

Reinterpretation Material: Vertex-level Efficiency for R [180, 300] mm

Cutflow (acceptance x efficiency) for the High-pT SR for representative points in the strong gluino pair production model. See additional resources for more information.

Cutflow (acceptance x efficiency) for the Trackless SR for representative points in the electroweak pair production model. See additional resources for more information.

Cutflow (acceptance x efficiency) for the Trackless SR for representative points in the electroweak pair production model with heavy-flavor quarks. See additional resources for more information.

Cutflow (acceptance x efficiency) for the High-pT SR for representative points in the electroweak pair production model with heavy-flavor quarks. See additional resources for more information.

Expected (dashed line) value of $-2\ln\Lambda$ as a function of $\kappa_{2V}$, when all other coupling modifiers are fixed to their SM predictions.

Likelihood observed contours at 68% CL (solid line) and 95% CL (dashed line) in the $\kappa_{\lambda}$, $\kappa_{2V}$ parameter space, when all other coupling modifiers are fixed to their SM predictions. The best-fit value, denoted by a cross in the plot, corresponds to the very first entry in the table.

Likelihood expected contours at 68% CL (teal shaded region) and 95% CL (yellow shaded region) in the $\kappa_{\lambda}$, $\kappa_{2V}$ parameter space, when all other coupling modifiers are fixed to their SM predictions.In the plot, the SM prediction is indicated by the star.

Likelihood observed contours at 68% CL (solid line) and 95% CL (dashed line) in the $c_{gghh}$ versus $c_{hhh}$ HEFT parameter space, with the remaining coefficient fixed to its SM value. The best-fit value, denoted by a cross in the plot, corresponds to the very first entry in the table.

Likelihood expected contours at 68% CL (teal shaded region) 95% CL (yellow shaded region) in the $c_{gghh}$ versus $c_{hhh}$ HEFT parameter space, with the remaining coefficient fixed to its SM value. In the plot, the SM prediction is indicated by the star.

Likelihood observed contours at 68% CL (solid line) and 95% CL (dashed line) in the $c_{tthh}$ versus $c_{hhh}$ HEFT parameter space, with the remaining coefficient fixed to its SM value. The best-fit value, denoted by a cross in the plot, corresponds to the very first entry in the table.

Likelihood exptected contours at 68% CL (teal shaded region) and 95% CL (yellow shaded region) in the $c_{tthh}$ versus $c_{hhh}$ HEFT parameter space, with the remaining coefficient fixed to its SM value. In the plot, the SM prediction is indicated by the star.

The observed (filled circles) and expected (hollow circles) 95% CL upper limits on the $HH$ ggF production cross-section in the Standard Model and for seven HEFT benchmark points. The teal and yellow colored bands indicate the $\pm 1\sigma$ and $\pm 2\sigma$ variations on the expected limit due to statistical and systematic uncertainties. The predicted cross-sections of each of the models under consideration are shown by the red crosses. The contribution from VBF production to the total $HH$ production cross-section is neglected.

Likelihood observed contours at 68% CL (solid line) and 95% CL (dashed line) in the $c_{{H}\boxed{}}$ versus $c_{H}$ SMEFT parameter space, with the remaining coefficient fixed to its SM value. The best-fit value, denoted by a cross in the plot, corresponds to the very first entry in the table.

Likelihood expected contours at 68% CL (teal shaded region) 95% CL (yellow shaded region) in the $c_{{H}\boxed{}}$ versus $c_{H}$ SMEFT parameter space, with the remaining coefficient fixed to its SM value. In the plot, the SM prediction is indicated by the star.

The acceptance times efficiency for the signal ggF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category. The other coupling modifiers affecting ggF $HH$ production are fixed to their SM predictions. The bands indicate the simulation statistical uncertainty.

The acceptance times efficiency for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category. When varying $\kappa_{\lambda}$, the other coupling modifiers affecting the VBF $HH$ production mode are set to their SM predictions. The bands indicate the simulation statistical uncertainty.

The acceptance times efficiency for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{2V}$ in each analysis category. When varying $\kappa_{2V}$, the other coupling modifiers affecting the VBF $HH$ production mode are set to their SM predictions. The bands indicate the simulation statistical uncertainty.

Expected yields of the signal $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ after applying the analysis preselection and in each analysis category after the final selection. The $HH$ yield is obtained considering both the $ggF$ and $VBF$ contributions. The other coupling modifiers affecting $HH$ production are fixed to their SM predictions. The bands indicate the simulation statistical uncertainty.

The acceptance times efficiency for the signal ggF $HH$ process in each analysis category as a function of the $c_{hhh}$ HEFT coefficients. The dashed lines denote values that are excluded at 95% CL. The bottom panels show the efficiency of the sum of all analysis categories.

The acceptance times efficiency for the signal ggF $HH$ process in each analysis category as a function of the $c_{tthh}$ HEFT coefficients. The dashed lines denote values that are excluded at 95% CL. The bottom panels show the efficiency of the sum of all analysis categories.

The acceptance times efficiency for the signal ggF $HH$ process in each analysis category as a function of the $c_{gghh}$ HEFT coefficients. The dashed lines denote values that are excluded at 95% CL. The bottom panels show the efficiency of the sum of all analysis categories.

The acceptance times efficiency for the signal ggF $HH$ process in each analysis category as a function of the $c_{H}$ SMEFT coefficients. The dashed lines denote values that are excluded at 95% CL. The bottom panels show the efficiency of the sum of all analysis categories.

The acceptance times efficiency for the signal ggF $HH$ process in each analysis category as a function of the $c_{{H}\boxed{}}$ SMEFT coefficients. The dashed lines denote values that are excluded at 95% CL. The bottom panels show the efficiency of the sum of all analysis categories.

The yield of the signal ggF $HH$ process in each analysis category as a function of the $c_{hhh}$ HEFT coefficients. The dashed lines denote values that are excluded at 95% CL. The bottom panels show the efficiency of the sum of all analysis categories.

The yield of the signal ggF $HH$ process in each analysis category as a function of the $c_{tthh}$ HEFT coefficients. The dashed lines denote values that are excluded at 95% CL. The bottom panels show the efficiency of the sum of all analysis categories.

The yield of the signal ggF $HH$ process in each analysis category as a function of the $c_{gghh}$ HEFT coefficients. The dashed lines denote values that are excluded at 95% CL. The bottom panels show the efficiency of the sum of all analysis categories.

The yield of the signal ggF $HH$ process in each analysis category as a function of the $c_{H}$ SMEFT coefficients. The dashed lines denote values that are excluded at 95% CL. The bottom panels show the efficiency of the sum of all analysis categories.

The yield of the signal ggF $HH$ process in each analysis category as a function of the $c_{{H}\boxed{}}$ SMEFT coefficients. The dashed lines denote values that are excluded at 95% CL. The bottom panels show the efficiency of the sum of all analysis categories.

Distributions of the diphoton invariant mass for events in data (black dots with error bars) compared to the sum of the expected signal and backgrounds (histograms) in the low mass 1 region.

Distributions of the diphoton invariant mass for events in data (black dots with error bars) compared to the sum of the expected signal and backgrounds (histograms) in the low mass 2 region.

Distributions of the diphoton invariant mass for events in data (black dots with error bars) compared to the sum of the expected signal and backgrounds (histograms) in the low mass 3 region.

Distributions of the diphoton invariant mass for events in data (black dots with error bars) compared to the sum of the expected signal and backgrounds (histograms) in the low mass 4 region.

Distributions of the diphoton invariant mass for events in data (black dots with error bars) compared to the sum of the expected signal and backgrounds (histograms) in the high mass 1 region.

Distributions of the diphoton invariant mass for events in data (black dots with error bars) compared to the sum of the expected signal and backgrounds (histograms) in the high mass 2 region.

Distributions of the diphoton invariant mass for events in data (black dots with error bars) compared to the sum of the expected signal and backgrounds (histograms) in the high mass 3 region.

Observed (solid line) value of $-2\ln\Lambda$ as a function of the $c_{hhh}$ HEFT coefficients when all other coupling modifiers are fixed to the SM value

Expected (dashed line) value of $-2\ln\Lambda$ as a function of the $c_{hhh}$ HEFT coefficients when all other coupling modifiers are fixed to the SM value

Observed (solid line) value of $-2\ln\Lambda$ as a function of the $c_{gghh}$ HEFT coefficients when all other coupling modifiers are fixed to the SM value

Expected (dashed line) value of $-2\ln\Lambda$ as a function of the $c_{gghh}$ HEFT coefficients when all other coupling modifiers are fixed to the SM value

Observed (solid line) value of $-2\ln\Lambda$ as a function of the $c_{tthh}$ HEFT coefficients when all other coupling modifiers are fixed to the SM value

Expected (dashed line) value of $-2\ln\Lambda$ as a function of the $c_{tthh}$ HEFT coefficients when all other coupling modifiers are fixed to the SM value

Observed (solid line) value of $-2\ln\Lambda$ as a function of the $c_{H}$ SMEFT coefficients when all other coupling modifiers are fixed to the SM value.

Expected (dashed line) value of $-2\ln\Lambda$ as a function of the $c_{H}$ SMEFT coefficients when all other coupling modifiers are fixed to the SM value.

Observed (solid line) value of $-2\ln\Lambda$ as a function of the $c_{{H}\boxed{}}$ SMEFT coefficients when all other coupling modifiers are fixed to the SM value.

Expected (dashed line) value of $-2\ln\Lambda$ as a function of the $c_{{H}\boxed{}}$ SMEFT coefficients when all other coupling modifiers are fixed to the SM value.

Summary of the expected (blue) and observed (orange) one-dimensional constraints at 68% (solid error bars) and 95% (dashed error bars) CL on the coefficients of selected HEFT (top) and SMEFT (bottom) operators describing Higgs boson interactions affecting $HH$ production through gluon--gluon fusion. The results are obtained from one dimensional scans of the profile log-likelihood assuming that all other Wilson coefficients are fixed to their SM values. The contribution from VBF $HH$ production is neglected.

Centrality dependence of (a) $v_2${2} and (b) $v_2${4}. (a) The red points indicate no pseudorapidity gap whereas the magenta points indicate a pseudorapidity gap of |$\Delta\eta$| > 2.0. (b) The black points indicate $v_2${4} with no pseudorapidity gap, the blue points indicate a two-subevent method with |$\Delta\eta$| > 2.0 but where some short-range pairs are allowed, and the red points indicate a two-subevent method with |$\Delta\eta$| > 2.0 where no short-range pairs are allowed.

Centrality dependence of (a) $v_2${2} and (b) $v_2${4}. (a) The red points indicate no pseudorapidity gap whereas the magenta points indicate a pseudorapidity gap of |$\Delta\eta$| > 2.0. (b) The black points indicate $v_2${4} with no pseudorapidity gap, the blue points indicate a two-subevent method with |$\Delta\eta$| > 2.0 but where some short-range pairs are allowed, and the red points indicate a two-subevent method with |$\Delta\eta$| > 2.0 where no short-range pairs are allowed.

Multi-particle $v_2$ as a function of centrality in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The magenta open diamonds indicate the $v_2${SP }, the blue open squares indicate $v_2${4}, the black open circles indicate $v_2${6}, and the green filled diamonds indicate $v_2${8}.

Multi-particle $v_2$ as a function of centrality in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The magenta open diamonds indicate the $v_2${SP }, the blue open squares indicate $v_2${4}, the black open circles indicate $v_2${6}, and the green filled diamonds indicate $v_2${8}.

Multi-particle $v_2$ as a function of centrality in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The magenta open diamonds indicate the $v_2${SP }, the blue open squares indicate $v_2${4}, the black open circles indicate $v_2${6}, and the green filled diamonds indicate $v_2${8}.

Multi-particle $v_2$ as a function of centrality in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The magenta open diamonds indicate the $v_2${SP }, the blue open squares indicate $v_2${4}, the black open circles indicate $v_2${6}, and the green filled diamonds indicate $v_2${8}.

Cumulant method estimate of $\sigma_{v_{2}}$ / $\langle v_{2} \rangle$ as a function of centrality in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The data are shown as black open squares. The same calculation as done in data is done in ampt, shown as a solid green line. Calculations of $\sigma_{\epsilon_{2}}$ /$ \langle\epsilon_{2}\rangle$ performed in the Monte Carlo Glauber model are shown as blue lines. The solid blue line is the Monte Carlo Glauber calculation done using the same estimate as the data, the dashed blue line is the direct calculation of the moments of the MC Glauber $\epsilon_{2}$ distribution.

Centrality dependence of $v_3${2, |$\Delta\eta$| > 2} in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV, shown as magenta diamonds. The systematic uncertainty is indicated as a shaded magenta band. Also shown as black squares are $\sqrt{\langle v^{2}_{3}\rangle}$ as determined from the folding analysis, which is shown in the next section.

Centrality dependence of $v_3${2, |$\Delta\eta$| > 2} in Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV, shown as magenta diamonds. The systematic uncertainty is indicated as a shaded magenta band. Also shown as black squares are $\sqrt{\langle v^{2}_{3}\rangle}$ as determined from the folding analysis, which is shown in the next section.

Centrality dependence of $c_3${4} for Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. (a) Calculations using both arms: $c_3${4} (black circles), $c_3${4}$_{ab|ab}$ (blue diamonds), $c_3${4}$_{aa|bb}$ (red squares), and comparison to STAR [36] (black stars). (b) Comparison of $c_3${4} determined using both arms (open symbols) and a single arm (closed symbols). Note that the open black circles are the same in (a) and (b).

Centrality dependence of $c_3${4} for Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. (a) Calculations using both arms: $c_3${4} (black circles), $c_3${4}$_{ab|ab}$ (blue diamonds), $c_3${4}$_{aa|bb}$ (red squares), and comparison to STAR [36] (black stars). (b) Comparison of $c_3${4} determined using both arms (open symbols) and a single arm (closed symbols). Note that the open black circles are the same in (a) and (b).

Centrality dependence of $c_3${4} for Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. (a) Calculations using both arms: $c_3${4} (black circles), $c_3${4}$_{ab|ab}$ (blue diamonds), $c_3${4}$_{aa|bb}$ (red squares), and comparison to STAR [36] (black stars). (b) Comparison of $c_3${4} determined using both arms (open symbols) and a single arm (closed symbols). Note that the open black circles are the same in (a) and (b).

Folding results for (a) $v_2$, (b) $\sigma_{v_{2}}$ , and (c) $\sigma_{v_{2}}$ / $\langle v_{2}\rangle$. The black lines above and below the points indicate the systematic uncertainties. The red (green) boxes indicate the statistical uncertainties at the 68.27% (95.45%) confidence level. In the case of $\langle v_{2}\rangle$, the statistical uncertainties at the 68.27% confidence level are too small to be seen, and the uncertainties at the 95.45% confidence level are visible but noticeably smaller than the marker size. Shown as blue squares are the same quantities as determined using the cumulant based calculation—these points are slightly offset in the x-coordinate to improve visibility.

Folding results for (a) $v_2$, (b) $\sigma_{v_{2}}$ , and (c) $\sigma_{v_{2}}$ / $\langle v_{2}\rangle$. The black lines above and below the points indicate the systematic uncertainties. The red (green) boxes indicate the statistical uncertainties at the 68.27% (95.45%) confidence level. In the case of $\langle v_{2}\rangle$, the statistical uncertainties at the 68.27% confidence level are too small to be seen, and the uncertainties at the 95.45% confidence level are visible but noticeably smaller than the marker size. Shown as blue squares are the same quantities as determined using the cumulant based calculation—these points are slightly offset in the x-coordinate to improve visibility.

Folding results for (a) $v_2$, (b) $\sigma_{v_{2}}$ , and (c) $\sigma_{v_{2}}$ / $\langle v_{2}\rangle$. The black lines above and below the points indicate the systematic uncertainties. The red (green) boxes indicate the statistical uncertainties at the 68.27% (95.45%) confidence level. In the case of $\langle v_{2}\rangle$, the statistical uncertainties at the 68.27% confidence level are too small to be seen, and the uncertainties at the 95.45% confidence level are visible but noticeably smaller than the marker size. Shown as blue squares are the same quantities as determined using the cumulant based calculation—these points are slightly offset in the x-coordinate to improve visibility.

Folding results for (a) $\langle v_{3}\rangle$, (b) $\sigma_{v_{3}}$ , and (c) $\sigma_{v_{3}}$ /$\langle v_{3}\rangle$. The black lines above and below the points indicate the systematic uncertainties. The red (green) boxes indicate the statistical uncertainties at the 68.27% (95.45%) confidence level. The $\sigma_{v_{3}}$ /$\langle v_{3}\rangle$ values are all ≈ 0.52, the apparent limiting value of this quantity for the Bessel-Gaussian distribution.

Folding results for (a) $\langle v_{3}\rangle$, (b) $\sigma_{v_{3}}$ , and (c) $\sigma_{v_{3}}$ /$\langle v_{3}\rangle$. The black lines above and below the points indicate the systematic uncertainties. The red (green) boxes indicate the statistical uncertainties at the 68.27% (95.45%) confidence level. The $\sigma_{v_{3}}$ /$\langle v_{3}\rangle$ values are all ≈ 0.52, the apparent limiting value of this quantity for the Bessel-Gaussian distribution.

Folding results for (a) $\langle v_{3}\rangle$, (b) $\sigma_{v_{3}}$ , and (c) $\sigma_{v_{3}}$ /$\langle v_{3}\rangle$. The black lines above and below the points indicate the systematic uncertainties. The red (green) boxes indicate the statistical uncertainties at the 68.27% (95.45%) confidence level. The $\sigma_{v_{3}}$ /$\langle v_{3}\rangle$ values are all ≈ 0.52, the apparent limiting value of this quantity for the Bessel-Gaussian distribution.

$v_2$ for inclusive charged hadrons in Au+Au at $\sqrt{s_{NN}}$ = 200 GeV.

$v_2$ for inclusive charged hadrons in Au+Au at $\sqrt{s_{NN}}$ = 200 GeV.

$v_2$ for inclusive charged hadrons in Au+Au at $\sqrt{s_{NN}}$ = 200 GeV.

$v_2$ for inclusive charged hadrons in Au+Au at $\sqrt{s_{NN}}$ = 200 GeV.

$v_2$ for inclusive charged hadrons in Au+Au at $\sqrt{s_{NN}}$ = 200 GeV.

$v_2$ for inclusive charged hadrons in Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV.

$v_2$ for inclusive charged hadrons in Au+Au at $\sqrt{s_{NN}}$ = 62.4 GeV.

$v_2$ for inclusive charged hadrons in Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV compared with 200 GeV.

$v_2$ for inclusive charged hadrons in Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV compared with 200 GeV.

$v_2$ for inclusive charged hadrons in Cu+Cu at $\sqrt{s_{NN}}$ = 62.4 GeV compared with 200 GeV.

Comparison of integrated $v_2$ at $\sqrt{s_{NN}}$ = 62.4 and 200 GeV. Au+Au reaction.

Comparison of integrated $v_2$ at $\sqrt{s_{NN}}$ = 62.4 and 200 GeV. Au+Au reaction.

Comparison of integrated $v_2$ at $\sqrt{s_{NN}}$ = 62.4 and 200 GeV. Cu+Cu reaction.

Comparison of integrated $v_2$ at $\sqrt{s_{NN}}$ = 62.4 and 200 GeV. Cu+Cu reaction.

The comparison of integrated $v_2$ as a function of centrality.

The comparison of the normalized $v_2$/$\epsilon$ vs. centrality.

The comparison of integrated $v_2$ as a function of centrality.

The comparison of the normalized $v_2$/$\epsilon$ vs. centrality.

The comparison of integrated $v_2$ as a function of centrality.

The comparison of the normalized $v_2$/$\epsilon$ vs. centrality.

The comparison of integrated $v_2$ as a function of centrality.

The comparison of the normalized $v_2$/$\epsilon$ vs. centrality.

Comparison of $v_2$ ($p_T$) at 200 GeV for two example systems with different collision size.

Comparison of $v_2$ ($p_T$) at 200 GeV for two example systems with different collision size.

Comparison of $v_2$ ($p_T$) at 200 GeV for two example systems with different collision size.

$v_2$ vs. $p_T$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 62.4 and 200 GeV for centralities given.

$v_2$ vs. $p_T$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 62.4 and 200 GeV for centralities given.

$v_2$ vs. $p_T$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 62.4 and 200 GeV for centralities given.

$v_2$ vs. $p_T$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 62.4 and 200 GeV for centralities given.

$v_2$ vs. $p_T$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 62.4 and 200 GeV for centralities given.

$v_2$ vs. $p_T$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 62.4 and 200 GeV for centralities given.

$v_2$ vs. $p_T$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 62.4 and 200 GeV for centralities given.

$v_2$ vs. $p_T$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 62.4 and 200 GeV for centralities given.

$v_2$ vs. $p_T$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 62.4 and 200 GeV for centralities given.

$v_2$ vs. $p_T$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 62.4 and 200 GeV for centralities given.

Comparison of $v_2$ between $\sqrt{s_{NN}}$ = 62.4 GeV and 200 GeV for $\pi$/$K$/$p$ emitted from central Au+Au collisions.

Comparison of $v_2$ between $\sqrt{s_{NN}}$ = 62.4 GeV and 200 GeV for $\pi$/$K$/$p$ emitted from central Au+Au collisions.

Comparison of $v_2$ between $\sqrt{s_{NN}}$ = 62.4 GeV and 200 GeV for $\pi$/$K$/$p$ emitted from central Au+Au collisions.

Comparison of $v_2$ between $\sqrt{s_{NN}}$ = 62.4 GeV and 200 GeV for $\pi$/$K$/$p$ emitted from central Au+Au collisions.

Comparison of $v_2$ between $\sqrt{s_{NN}}$ = 62.4 GeV and 200 GeV for $\pi$/$K$/$p$ emitted from central Au+Au collisions.

Comparison of $v_2$ between $\sqrt{s_{NN}}$ = 62.4 GeV and 200 GeV for $\pi$/$K$/$p$ emitted from central Au+Au collisions.

Comparison of $v_2$ between $\sqrt{s_{NN}}$ = 62.4 GeV and 200 GeV for $\pi$/$K$/$p$ emitted from central Au+Au collisions.

Comparison of the $v_2$ of particles, antiparticles, for a minimum bias sample at 200 GeV and central 62.4 GeV Au+Au collisions.

Comparison of the $v_2$ of particles, antiparticles, for a minimum bias sample at 200 GeV and central 62.4 GeV Au+Au collisions.

Comparison of the $v_2$ of particles, antiparticles, for a minimum bias sample at 200 GeV and central 62.4 GeV Au+Au collisions.

Comparison of the $v_2$ of particles, antiparticles, for a minimum bias sample at 200 GeV and central 62.4 GeV Au+Au collisions.

Comparison of the $v_2$ of particles, antiparticles, for a minimum bias sample at 200 GeV and central 62.4 GeV Au+Au collisions.

Comparison of the $v_2$ of particles, antiparticles, for a minimum bias sample at 200 GeV and central 62.4 GeV Au+Au collisions.

The ratio $v_2$/$n_q$ vs. $p_T$/$n_q$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 200 GeV collisions.

The ratio $v_2$/$n_q$ vs. $p_T$/$n_q$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 200 GeV collisions.

The ratio $v_2$/$n_q$ vs. $p_T$/$n_q$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 200 GeV collisions.

The ratio $v_2$/$n_q$ vs. $p_T$/$n_q$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 200 GeV collisions.

The ratio $v_2$/$n_q$ vs. $p_T$/$n_q$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 200 GeV collisions.

The ratio $v_2$/$n_q$ vs. $p_T$/$n_q$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 200 GeV collisions.

The ratio $v_2$/$n_q$ vs. $p_T$/$n_q$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 200 GeV collisions.

The ratio $v_2$/$n_q$ vs. $p_T$/$n_q$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 200 GeV collisions.

The ratio $v_2$/$n_q$ vs. $p_T$/$n_q$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 200 GeV collisions.

The ratio $v_2$/$n_q$ vs. $p_T$/$n_q$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 200 GeV collisions.

The ratio $v_2$/$n_q$ vs. ${KE}_T$/$n_q$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 200 GeV collisions.

The ratio $v_2$/$n_q$ vs. ${KE}_T$/$n_q$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 200 GeV collisions.

The ratio $v_2$/$n_q$ vs. ${KE}_T$/$n_q$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 200 GeV collisions.

The ratio $v_2$/$n_q$ vs. ${KE}_T$/$n_q$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 200 GeV collisions.

The ratio $v_2$/$n_q$ vs. ${KE}_T$/$n_q$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 200 GeV collisions.

The ratio $v_2$/$n_q$ vs. ${KE}_T$/$n_q$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 200 GeV collisions.

The ratio $v_2$/$n_q$ vs. ${KE}_T$/$n_q$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 200 GeV collisions.

The ratio $v_2$/$n_q$ vs. ${KE}_T$/$n_q$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 200 GeV collisions.

The ratio $v_2$/$n_q$ vs. ${KE}_T$/$n_q$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 200 GeV collisions.

The ratio $v_2$/$n_q$ vs. ${KE}_T$/$n_q$ for $\pi$/$K$/$p$ emitted from Au+Au at 62.4 and 200 GeV and Cu+Cu at 200 GeV collisions.

The ratio $v_2$/$n_q$ vs. ${KE}_T$/$n_q$ for the indicated hadrons emitted from central Au+Au collisions at 62.4 GeV.

The ratio $v_2$/$n_q$ vs. ${KE}_T$/$n_q$ for the indicated hadrons emitted from central Au+Au collisions at 62.4 GeV.

The ratio $v_2$/$n_q$ vs. ${KE}_T$/$n_q$ for the indicated hadrons emitted from central Au+Au collisions at 62.4 GeV.

$v_2$ vs. $p_T$ and $v_2$/($\epsilon * N^{1/3}_{part} * n_q$) vs. ${KE}_T$/$n_q$ for $\pi$/$K$/$p$ in Au+Au at 200 GeV, in Au+Au at 62.4 GeV, and in Cu+Cu at 200 GeV. The values of $v_2$ and $p_T$ in Au+Au at 200 GeV, in Au+Au at 62.4 GeV, and in Cu+Cu at 200 GeV are the same for as figure 14, and the values of $v_2$, $n_q$, and $KE_T$ in Au+Au at 200 GeV, in Au+Au at 62.4 GeV, and in Cu+Cu at 200 GeV are the same for as figure 18.

Expected exclusion limits at 95% CL from Fig 7(a) for $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

One $\sigma$ band of expected exclusion limits at 95% CL from Fig 7(a) for $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Observed exclusion limits at 95% CL from Fig 7(c) for $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Positive one $\sigma$ observed exclusion limits at 95% CL from Fig 7(c) for $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Negative one $\sigma$ observed exclusion limits at 95% CL from Fig 7(c) for $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Expected exclusion limits at 95% CL from Fig 7(c) for $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

One $\sigma$ band of expected exclusion limits at 95% CL from Fig 7(c) for $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Observed exclusion limits at 95% CL from Fig 7(f) for $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Positive one $\sigma$ observed exclusion limits at 95% CL from Fig 7(f) for $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Negative one $\sigma$ observed exclusion limits at 95% CL from Fig 7(f) for $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Expected exclusion limits at 95% CL from Fig 7(f) for $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

One $\sigma$ band of expected exclusion limits at 95% CL from Fig 7(f) for $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Observed exclusion limits at 95% CL from Fig 7(e) for direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Positive one $\sigma$ observed exclusion limits at 95% CL from Fig 7(e) for direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Negative one $\sigma$ observed exclusion limits at 95% CL from Fig 7(e) for direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Expected exclusion limits at 95% CL from Fig 7(e) for direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

One $\sigma$ band of expected exclusion limits at 95% CL from Fig 7(e) for direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Observed exclusion limits at 95% CL from Fig 7(b) for $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Positive one $\sigma$ observed exclusion limits at 95% CL from Fig 7(b) for $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Negative one $\sigma$ observed exclusion limits at 95% CL from Fig 7(b) for $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Expected exclusion limits at 95% CL from Fig 7(b) for $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

One $\sigma$ band of expected exclusion limits at 95% CL from Fig 7(b) for $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Observed exclusion limits at 95% CL from Fig 7(d) for $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Positive one $\sigma$ observed exclusion limits at 95% CL from Fig 7(d) for $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Negative one $\sigma$ observed exclusion limits at 95% CL from Fig 7(d) for $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Expected exclusion limits at 95% CL from Fig 7(d) for $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

One $\sigma$ band of expected exclusion limits at 95% CL from Fig 7(d) for $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

N-1 distribution for $m_{\mathrm{eff}}$of observed data and expected background in SRGGWZ-H.

N-1 distribution for $E_{\mathrm{T}}^{\mathrm{miss}}$of observed data and expected background in SRGGSlep-M.

N-1 distribution for $\sum{p_{\mathrm{T}}^\mathrm{jet}}$of observed data and expected background in SRUDD-ge2b.

N-1 distribution for $m_{\mathrm{eff}}$of observed data and expected background in SRLQD.

N-1 distribution for $m_{\mathrm{eff}}$of observed data and expected background in SRSSWZ-H.

N-1 distribution for $m_{\mathrm{eff}}$of observed data and expected background in SRSSSlep-H(loose).

Signal acceptance for SRGGWZ-H signal region from Fig 10(c) in a SUSY scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRGGWZ-H signal region from Fig 15(c) in a SUSY scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRGGWZ-M signal region from Fig 10(b) in a SUSY scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRGGWZ-M signal region from Fig 15(b) in a SUSY scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRGGWZ-L signal region from Fig 10(a) in a SUSY scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRGGWZ-L signal region from Fig 15(a) in a SUSY scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRGGSlep-L signal region from Fig 12(a) in a SUSY scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRGGSlep-L signal region from Fig 17(a) in a SUSY scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRGGSlep-M signal region from Fig 12(b) in a SUSY scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRGGSlep-M signal region from Fig 17(b) in a SUSY scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRGGSlep-H signal region from Fig 12(c) in a SUSY scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRGGSlep-H signal region from Fig 17(c) in a SUSY scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRUDD-1b signal region from Fig 14(b) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal efficiency for SRUDD-1b signal region from Fig 19(b) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal acceptance for SRUDD-2b signal region from Fig 14(c) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal efficiency for SRUDD-2b signal region from Fig 19(c) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal acceptance for SRUDD-ge2b signal region from Fig 14(d) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal efficiency for SRUDD-ge2b signal region from Fig 19(d) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal acceptance for SRUDD-ge3b signal region from Fig 14(e) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal efficiency for SRUDD-ge3b signal region from Fig 19(e) in a SUSY scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Signal acceptance for SRLQD signal region from Fig 14(a) in a SUSY scenario where direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Signal efficiency for SRLQD signal region from Fig 19(a) in a SUSY scenario where direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Signal acceptance for SRSSWZ-L signal region from Fig 11(a) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSWZ-L signal region from Fig 16(a) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSWZ-ML signal region from Fig 11(b) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSWZ-ML signal region from Fig 16(b) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSWZ-MH signal region from Fig 11(c) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSWZ-MH signal region from Fig 16(c) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSWZ-H signal region from Fig 11(d) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSWZ-H signal region from Fig 16(d) in a SUSY scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSSlep-H signal region from Fig 13(d) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSSlep-H signal region from Fig 18(d) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSSlep-MH signal region from Fig 13(c) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSSlep-MH signal region from Fig 18(c) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSSlep-L signal region from Fig 13(a) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSSlep-L signal region from Fig 18(a) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSSlep-ML signal region from Fig 13(b) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSSlep-ML signal region from Fig 18(b) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal acceptance for SRSSSlep-H(loose) signal region from Fig 13(e) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Signal efficiency for SRSSSlep-H(loose) signal region from Fig 18(e) in a SUSY scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRGGWZ-H in a susy scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1400 GeV, $m(\tilde{\chi_{1}^{0}})$ = 1000 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRGGWZ-M in a susy scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1400 GeV, $m(\tilde{\chi_{1}^{0}})$ = 1000 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRGGWZ-L in a susy scenario where $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1400 GeV, $m(\tilde{\chi_{1}^{0}})$ = 1000 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRGGSlep-L in a susy scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 2000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 500 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRGGSlep-M in a susy scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 2000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 500 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRGGSlep-H in a susy scenario where $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 2000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 500 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRUDD-1b in a susy scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1600 GeV, $m(\tilde{t})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRUDD-2b in a susy scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1600 GeV, $m(\tilde{t})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRUDD-ge2b in a susy scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1600 GeV, $m(\tilde{t})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRUDD-ge3b in a susy scenario where $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 1600 GeV, $m(\tilde{t})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRLQD in a susy scenario where direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$. The masses of the superpartners involved in the process are set to $m(\tilde{g})$ = 2200 GeV, $m(\tilde{\chi_{1}^{0}})$ = 1870 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSWZ-L in a susy scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 800 GeV, $m(\tilde{\chi_{1}^{0}})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSWZ-ML in a susy scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 800 GeV, $m(\tilde{\chi_{1}^{0}})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSWZ-MH in a susy scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 800 GeV, $m(\tilde{\chi_{1}^{0}})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSWZ-H in a susy scenario where $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 800 GeV, $m(\tilde{\chi_{1}^{0}})$ = 600 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSSlep-H in a susy scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 1000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 800 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSSlep-MH in a susy scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 1000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 800 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSSlep-L in a susy scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 1000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 800 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSSlep-ML in a susy scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 1000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 800 GeV. Only statistical uncertainties are shown.

Number of signal events expected for 139 fb$^{-1}$ at different stages of the event selection for the signal region SRSSSlep-H(loose) in a susy scenario where $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$. The masses of the superpartners involved in the process are set to $m(\tilde{q})$ = 1000 GeV, $m(\tilde{\chi_{1}^{0}})$ = 800 GeV. Only statistical uncertainties are shown.

Cross-section upper limits at 95% CL from Fig1(a) for $\tilde{g}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Cross-section upper limits at 95% CL from Fig1(c) for $\tilde{g}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

Cross-section upper limits at 95% CL from Fig1(f) for $\tilde{g}$ decays into anti-top and $\tilde{t}$ and $\tilde{t}$ decays via a non-zero RPV coupling $\lambda''$

Cross-section upper limits at 95% CL from Fig1(e) for direct $\tilde{\chi_{1}^{0}}$ decay into SM leptons and quarks via a non-zero RPV coupling $\lambda'$

Cross-section upper limits at 95% CL from Fig1(b) for $\tilde{q}$ decays into SM gauge bosons and $\tilde{\chi}^{0}_{1}$

Cross-section upper limits at 95% CL from Fig1(d) for $\tilde{q}$ decays into sleptons and subsequently to SM leptons and $\tilde{\chi}^{0}_{1}$

The expected minus one standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $ee$ channel for the case with mass thresholds.

The expected plus two standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $ee$ channel for the case with mass thresholds.

The expected minus two standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $ee$ channel for the case with mass thresholds.

The observed exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case with mass thresholds.

The median expected exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case with mass thresholds.

The expected plus one standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case with mass thresholds.

The expected minus one standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case with mass thresholds.

The expected plus two standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case with mass thresholds.

The expected minus two standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case with mass thresholds.

The observed exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $ee$ channel for the case without mass thresholds.

The median expected exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $ee$ channel for the case without mass thresholds.

The expected plus one standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $ee$ channel for the case without mass thresholds.

The expected minus one standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $ee$ channel for the case without mass thresholds.

The expected plus two standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $ee$ channel for the case without mass thresholds.

The expected minus two standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $ee$ channel for the case without mass thresholds.

The observed exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case without mass thresholds.

The median expected exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case without mass thresholds.

The expected plus one standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case without mass thresholds.

The expected minus one standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case without mass thresholds.

The expected plus two standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case without mass thresholds.

The expected minus two standard deviation exclusion limit at 95% CL for the clockwork gravity model projected in the $k–M_{5}$ parameter space for the $\gamma\gamma$ channel for the case without mass thresholds.

The minimum $\chi^2$ and the parameters $T_{fo}$ and $\beta^{max}_T$ for each of the five centrality selections. The best fit parameters are determined by averaging all parameter pairs within the 1$\sigma$ contour. The errors correspond to the standard deviation of the parameter pairs within the 1$\sigma$ $\chi^2$ contour.

Fit parameters for each particle species using Equations $\frac{dN}{dy}$ = $C \cdot$ ($N_{part}$)$^{\alpha_{part}}$ and $\frac{dN}{dy}$ = $C^{\prime} \cdot$ ($N_{coll}$)$^{\alpha_{coll}}$.

Values of the parameters $n_{pp}$ and $x$ from fitting Equation $R \equiv \frac{dN/dy}{N_{part}}$ = (1-$x$) $\cdot n_{pp} \frac{1}{2}$ + $x \cdot n_{pp} \frac{N_{coll}}{N_{part}}$ = $n_{pp}$ [$\frac{1}{2}$ + $x$($\frac{N_{coll}}{N_{part}}$ - $\frac{1}{2})$]to the observed $dN$/$dy$ per $N_{part}$.

Invariant yields for $\pi^{\pm}$ measured in minimum-bias events at midrapidity and normalized to one rapidity unit.

Invariant yields for $K^{\pm}$ measured in minimum-bias events at midrapidity and normalized to one rapidity unit.

Invariant yields for (anti)$p$ measured in minimum-bias events at midrapidity and normalized to one rapidity unit.

Pion invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

Pion invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

Kaon invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

Kaon invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

(Anti)proton invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

(Anti)proton invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

(Anti)proton invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

(Anti)proton invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

(Anti)proton invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

(Anti)proton invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

(Anti)proton invariant yields in each event centrality and $p_T$ bin measured at midrapidity, normalized to one rapidity unit.

Minimum-bias invariant yields for $\pi^{\pm}$ in equal $p_T$ bins.

Minimum-bias invariant yields for $K^{\pm}$ in equal $p_T$ bins.

Minimum-bias invariant yields for (anti)$p$ in equal $p_T$ bins.

Pion invariant yields in each event centrality normalized to one rapidity unit at midrapidity.

Pion invariant yields in each event centrality normalized to one rapidity unit at midrapidity.

Kaon invariant yields in each event centrality normalized to one rapidity unit at midrapidity.

(Anti)proton invariant yields in each event centrality normalized to one rapidity unit at midrapidity.

(Anti)proton invariant yields in each event centrality normalized to one rapidity unit at midrapidity.

The integrated mean $p_T$ for pions, kaons, and (anti)protons produced in the five different classes of event centrality.

The expansion parameters $T_{fo}$ and $\beta^{max}_T$ as a function of the number of participants.

<p>Charged hadron spectra for centrality selected d+Au collisions.</p>

<p>Charged hadron spectra for centrality selected d+Au collisions.</p>

<p>Charged hadron spectra for centrality selected N+Au collisions. All centralities.</p>

<p>Charged hadron spectra for centrality selected N+Au collisions. Centrality A (Central).</p>

<p>Charged hadron spectra for centrality selected N+Au collisions. Centrality B.</p>

<p>Charged hadron spectra for centrality selected N+Au collisions. Centrality C.</p>

<p>Charged hadron spectra for centrality selected N+Au collisions. Centrality D (Peripheral).</p>

<p>$R_{dAu}$ as a function of $p_T$.</p>

<p>$R_{dAu}$ as a function of $p_T$.</p>

<p>$R_{dAu}$ as a function of $p_T$.</p>

<p>$R_{dAu}$ as a function of $p_T$.</p>

<p>$R_{NAu}$ as a function of $p_T$. Centrality A (Peripheral).</p>

<p>$R_{NAu}$ as a function of $p_T$. Centrality B.</p>

<p>$R_{dAu}$ as a function of $p_T$. Centrality C.</p>

<p>$R_{dAu}$ as a function of $p_T$. Centrality D (Central).</p>

<p>$R_{dAu}$ values averaged in three momentum ranges as functions of $\nu−1$. Horizontal bars show the uncertainty in the value of $\nu$ for each centrality class.</p>

<p>$R_{NAu}$ values averaged in three momentum ranges as functions of $\nu−1$. Horizontal bars show the uncertainty in the value of $\nu$ for each centrality class.</p>

<p>Charged hadron spectra for centrality selected d+Au collisions.</p>

Predicted and observed yields as a function of $m_{4\ell}$ in the VBS-Suppressed region. Overflow events are included in the last bin of the distribution.

Relative systematic uncertainties in the unfolded differential cross section for 4ljj production in the VBS-Enhanced region as a function of m_{jj}

Relative systematic uncertainties in the unfolded differential cross section for 4ljj production in the VBS-Enhanced region as a function of m_{4\ell}

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $m_{jj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $m_{jj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $m_{4\ell}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $m_{4\ell}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $ p_{T,4l}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $ p_{T,4l}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $\Delta\phi_{jj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $\Delta\phi_{jj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $|\Delta y_{jj}|$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $|\Delta y_{jj}|$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $cos(\theta^{*}_{12})$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $cos(\theta^{*}_{12})$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $cos(\theta^{*}_{34})$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $cos(\theta^{*}_{34})$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $p_{T,jj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $p_{T,jj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $p_{T,4ljj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $p_{T,4ljj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $S_{T,4ljj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $S_{T,4ljj}$. Overflow events are included in the last bin of the distribution.

Statistical correlation between different bins of the unfolded cross-section for all the measured observables in the VBS-Enhanced region.

Statistical correlation between different bins of the unfolded cross-section for all the measured observables in the VBS-Suppressed region.

The second-order Fourier coefficients $v_2 (p_T)$ for charge-combined identified hadrons $\pi^{\pm}$, $K^{\pm}$, $ p$, and $\bar{p}$ measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV for the centrality classes marked in each panel.

The second-order Fourier coefficients $v_2 (p_T)$ for charge-combined identified hadrons $\pi^{\pm}$, $K^{\pm}$, $ p$, and $\bar{p}$ measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV for the centrality classes marked in each panel.

The second-order Fourier coefficients $v_2 (p_T)$ for charge-combined identified hadrons $\pi^{\pm}$, $K^{\pm}$, $ p$, and $\bar{p}$ measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV for the centrality classes marked in each panel.

The second-order Fourier coefficients $v_2 (p_T)$ for charge-combined identified hadrons $\pi^{\pm}$, $K^{\pm}$, $ p$, and $\bar{p}$ measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV for the centrality classes marked in each panel.

The second-order Fourier coefficients $v_2 (p_T)$ for charge-combined identified hadrons $\pi^{\pm}$, $K^{\pm}$, $ p$, and $\bar{p}$ measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV for the centrality classes marked in each panel.

The third-order Fourier coefficients $v_3(p_T)$ for charge-combined identified hadrons $\pi^{\pm}$, $K^{\pm}$, $p$, and $\bar{p}$ measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV for 0%-30% centrality.

The third-order Fourier coefficients $v_3(p_T)$ for charge-combined identified hadrons $\pi^{\pm}$, $K^{\pm}$, $p$, and $\bar{p}$ measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV for 0%-30% centrality.

The third-order Fourier coefficients $v_3(p_T)$ for charge-combined identified hadrons $\pi^{\pm}$, $K^{\pm}$, $p$, and $\bar{p}$ measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV for 0%-30% centrality.

The first-order Fourier coefficients $v_1(p_T)$ for charge-combined $\pi^{\pm}$, $K^{\pm}$, $ p$, and $\bar{p}$ identified hadrons measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV for 10%-50% centrality.

The first-order Fourier coefficients $v_1(p_T)$ for charge-combined $\pi^{\pm}$, $K^{\pm}$, $ p$, and $\bar{p}$ identified hadrons measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV for 10%-50% centrality.

Scaled second-order Fourier coefficients $v_2(p_T)/\epsilon_2$ for charged hadrons measured at midrapidity in Cu + Cu collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/\epsilon_2$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/\epsilon_2$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/\epsilon_2$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/\epsilon_2$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/\epsilon_2$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/\epsilon_2$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/\epsilon_2$ for charged hadrons measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/(\epsilon_2 N^{1/3}_{part})$ for charged hadrons measured at midrapidity in Cu + Cu collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/(\epsilon_2 N^{1/3}_{part})$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/(\epsilon_2 N^{1/3}_{part})$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/(\epsilon_2 N^{1/3}_{part})$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/(\epsilon_2 N^{1/3}_{part})$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/(\epsilon_2 N^{1/3}_{part})$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/(\epsilon_2 N^{1/3}_{part})$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled second-order Fourier coefficients $v_2(p_T)/(\epsilon_2 N^{1/3}_{part})$ for charged hadrons measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled third-order Fourier coefficients $v_3(p_T)/\epsilon_3$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled third-order Fourier coefficients $v_3(p_T)/\epsilon_3$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled third-order Fourier coefficients $v_3(p_T)/\epsilon_3$ for charged hadrons measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled third-order Fourier coefficients $\frac{v_3}{($_3 N_{part}^{1/3})}$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled third-order Fourier coefficients $\frac{v_3}{(e_3 N_{part}^{1/3})}$ for charged hadrons measured at midrapidity in Au + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

Scaled third-order Fourier coefficients $\frac{v_3}{(e_3 N_{part}^{1/3})}$ for charged hadrons measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV.

The second-order Fourier coefficients $v_2$($p_T$ ) for charged hadrons measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV in comparison with hydrodynamics calculations for the centrality classes marked in each panel.

The third-order Fourier coefficients $v_3(p_T)$ for charged hadrons measured at midrapidity in Cu + Au collisions at $\sqrt{S_{NN}}$ = 200 GeV in comparison with hydrodynamics calculations for the centrality classes marked in each panel.

The ratio of $\psi(2S)$ to $J/\psi(1S)$ cross sections, as measured in the dimuon decay channel in pp collisions at 200 GeV over the rapidity interval $1.2<|y|<2.2$, as a function of transverse momentum $p_{T}$.

The double ratios $[\sigma_{\psi(2S)}/\sigma_{\psi(1S)}]_{pA}/[\sigma_{\psi(2S)}/\sigma_{\psi(1S)}]_{pp}$. The global systematic uncertainty is $\pm 14\%$.

The charged-particle multiplicity $dN_{ch}/d\eta$ determined from AMPT simulations and the nuclear overlap $S_T$ from Glauber simulations for the data shown in Fig. 5. Here we assume $dN/d\eta$ = $\frac{3}{2}dN_{ch}/d\eta$.

The charged-particle multiplicity $dN_{ch}/d\eta$ determined from AMPT simulations and the nuclear overlap $S_T$ from Glauber simulations for the data shown in Fig. 5. Here we assume $dN/d\eta$ = $\frac{3}{2}dN_{ch}/d\eta$.

The charged-particle multiplicity $dN_{ch}/d\eta$ determined from AMPT simulations and the nuclear overlap $S_T$ from Glauber simulations for the data shown in Fig. 5. Here we assume $dN/d\eta$ = $\frac{3}{2}dN_{ch}/d\eta$.

Data table for $\sigma_{\psi(2S)}/\sigma_{\psi(1S)}$ ratios, measured in the dimuon decay channel. Here "fwd" refers to the forward rapidity interval $1.2<y<2.2$, "bwd" refers to the backward rapidity interval $-2.2<y<-1.2$, and "comb" refers to the rapidity interval $1.2<|y|<2.2$.

Distribution of $m_{\mathrm{JJ}}$ in the SR (black points) compared to the background-only expectation, where the shape is derived from the CR and the normalisation is fitted in the SR.

Observed and expected limits at 95% CL on the cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model A (see the text for details). The darker and lighter shaded bands around the expected limits represent the $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty range, respectively.

The predicted theory cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model A (see the text for details). In the theory model, a universal coupling of the $Z^\prime$ to each SM quark $g_q$ is assumed and the dark quarks are all represented by one species with an effective coupling $g_{q_d}$ to the $Z^\prime$. For model A, the predicted cross-section is shown for $g_q=0.05$ and $g_{q_d} = 0.2$.

Observed and expected limits at 95% CL on the cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model B (see the text for details). The darker and lighter shaded bands around the expected limits represent the $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty range, respectively.

The predicted theory cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model B (see the text for details). In the theory model, a universal coupling of the $Z^\prime$ to each SM quark $g_q$ is assumed and the dark quarks are all represented by one species with an effective coupling $g_{q_d}$ to the $Z^\prime$. For model B, the predicted cross-section is shown for $g_q=0.15$ and $g_{q_d} = 0.5$.

Observed and expected limits at 95% CL on the cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model C (see the text for details). The darker and lighter shaded bands around the expected limits represent the $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty range, respectively.

The predicted theory cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model C (see the text for details). In the theory model, a universal coupling of the $Z^\prime$ to each SM quark $g_q$ is assumed and the dark quarks are all represented by one species with an effective coupling $g_{q_d}$ to the $Z^\prime$. For model C, the predicted cross-section is shown for $g_q=0.15$ and $g_{q_d}= 0.5$.

Observed and expected limits at 95% CL on the cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model D (see the text for details). The darker and lighter shaded bands around the expected limits represent the $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty range, respectively.

The predicted theory cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model D (see the text for details). In the theory model, a universal coupling of the $Z^\prime$ to each SM quark $g_q$ is assumed and the dark quarks are all represented by one species with an effective coupling $g_{q_d}$ to the $Z^\prime$. For model D, the predicted cross-section is shown for $g_q=0.05$ and $g_{q_d} = 0.2$.

Relative efficiency (in %) of each stage of the analysis selection with respect to the previous one for models A through D and $m_{Z^\prime}=2.5$ TeV. Approximately 90000 events were generated per signal mass hypothesis.

Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64822.9 $\pm$ 254.6, NPowHeg truth =634276.

Unfolded dR distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64273.2 $\pm$ 253.5, NPowHeg truth =338671.

Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 65362.4 $\pm$ 255.7 , NPowHeg truth =634214.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63486.8 $\pm$ 252.0, NPowHeg truth =333348.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64432.6 $\pm$ 253.8, NPowHeg truth =624059.

Unfolded $M(l^{+}\gamma)$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57095.9 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =76327.

Unfolded $M(l^{-}\gamma)$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57130 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =763273.

Unfolded $M(l^{+}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58310.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58320.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Unfolded dR distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57705.4 $\pm$ 254 , NPowHeg truth =301655, N Sherpa truth =763259.

Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58869.2 $\pm$ 256, NPowHeg truth =564660, NS herpa truth =813995.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57358.6 $\pm$ 254 , NPowHeg truth =298662, N Sherpa truth =756609.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58226.4 $\pm$ 254, NPowHeg truth =558949, N Sherpa truth =806988.

Unfolded $M(l^{+}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63717.4 $\pm$ 252.4, NPowHeg truth =338714.

Unfolded $M(l^{-}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63855.8 $\pm$ 252.7 , NPowHeg truth =338708.

Unfolded $M(l^{+}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64809.8 $\pm$ 254.6, NPowHeg truth =634285.

Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64822.9 $\pm$ 254.6, NPowHeg truth =634276.

Unfolded dR distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64273.2 $\pm$ 253.5, NPowHeg truth =338671.

Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 65362.4 $\pm$ 255.7 , NPowHeg truth =634214.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63486.8 $\pm$ 252.0, NPowHeg truth =333348.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64432.6 $\pm$ 253.8, NPowHeg truth =624059.

Unfolded $M(l^{+}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57095.9 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =76327.

Unfolded $M(l^{-}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57130 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =763273.

Unfolded $M(l^{+}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58310.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58320.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Unfolded dR distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57705.4 $\pm$ 254 , NPowHeg truth =301655, N Sherpa truth =763259.

Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58869.2 $\pm$ 256, NPowHeg truth =564660, NS herpa truth =813995.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57358.6 $\pm$ 254 , NPowHeg truth =298662, N Sherpa truth =756609.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58226.4 $\pm$ 254, NPowHeg truth =558949, N Sherpa truth =806988.

Truth $P_{T}^{Z}$ distribution after reweghting to data.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63717.4 $\pm$ 252.4, NPowHeg truth =338714.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63855.8 $\pm$ 252.7 , NPowHeg truth =338708.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64809.8 $\pm$ 254.6, NPowHeg truth =634285.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64822.9 $\pm$ 254.6, NPowHeg truth =634276.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64273.2 $\pm$ 253.5, NPowHeg truth =338671.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 65362.4 $\pm$ 255.7 , NPowHeg truth =634214.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63486.8 $\pm$ 252.0, NPowHeg truth =333348.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64432.6 $\pm$ 253.8, NPowHeg truth =624059.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57095.9 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =76327.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57130 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =763273.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58310.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58320.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57705.4 $\pm$ 254 , NPowHeg truth =301655, N Sherpa truth =763259.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58869.2 $\pm$ 256, NPowHeg truth =564660, NS herpa truth =813995.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57358.6 $\pm$ 254 , NPowHeg truth =298662, N Sherpa truth =756609.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58226.4 $\pm$ 254, NPowHeg truth =558949, N Sherpa truth =806988.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63717.4 $\pm$ 252.4, NPowHeg truth =338714.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63855.8 $\pm$ 252.7 , NPowHeg truth =338708.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64809.8 $\pm$ 254.6, NPowHeg truth =634285.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64822.9 $\pm$ 254.6, NPowHeg truth =634276.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64273.2 $\pm$ 253.5, NPowHeg truth =338671.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 65362.4 $\pm$ 255.7 , NPowHeg truth =634214.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63486.8 $\pm$ 252.0, NPowHeg truth =333348.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64432.6 $\pm$ 253.8, NPowHeg truth =624059.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57095.9 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =76327.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57130 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =763273.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58310.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58320.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57705.4 $\pm$ 254 , NPowHeg truth =301655, N Sherpa truth =763259.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58869.2 $\pm$ 256, NPowHeg truth =564660, NS herpa truth =813995.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57358.6 $\pm$ 254 , NPowHeg truth =298662, N Sherpa truth =756609.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58226.4 $\pm$ 254, NPowHeg truth =558949, N Sherpa truth =806988.

Unfolded $M(l^{+}\gamma\gamma)$ distribution for $Z \to ll\gamma\gamma$ process with dressed leptons.

Unfolded $M(l^{+}\gamma\gamma)$ distribution for $Z \to ll\gamma\gamma$ process with bare leptons.

Unfolded dR distribution for $Z \to ll\gamma\gamma$ process with dressed leptons for leading photon.

Unfolded dR distribution for $Z \to ll\gamma\gamma$ process with bare leptons for leading photon.

Unfolded dR distribution for $Z \to ll\gamma\gamma$ process with dressed leptons for second photon.

Unfolded dR distribution for $Z \to ll\gamma\gamma$ process with bare leptons for second photon.

Unfolded dRgg distribution for $Z \to ll\gamma\gamma$ process with dressed leptons.

Unfolded dRgg distribution for $Z \to ll\gamma\gamma$ process with bare leptons.

Unfolded $P_{T}^{\gamma 1}$ distribution for $Z \to ll\gamma\gamma$ process with dressed leptons for leading photon.

Unfolded $P_{T}^{\gamma 1}$ distribution for $Z \to ll\gamma\gamma$ process with bare leptons for leading photon.

Unfolded $P_{T}^{\gamma 2}$ distribution for $Z \to ll\gamma\gamma$ process with dressed leptons for second photon.

Unfolded $P_{T}^{\gamma 2}$ distribution for $Z \to ll\gamma\gamma$ process with bare leptons for second photon.

Measured $p_T$ cross sections in full-lepton phase space for 1.2 < |y| < 1.6.

Measured $p_T$ cross sections in full-lepton phase space for 1.6 < |y| < 2.0.

Measured $p_T$ cross sections in full-lepton phase space for 2.0 < |y| < 2.4.

Measured $p_T$ cross sections in full-lepton phase space for 2.4 < |y| < 2.8.

Measured $p_T$ cross sections in full-lepton phase space for 2.8 < |y| < 3.2.

Measured cross sections in full-lepton phase space as a function of |y|.

Normalised measured $p_T$ cross sections in full-lepton phase space for |y| < 1.6.

Efficiency corrected cumulants of net-charge distributions as a function of $\langle N_{part} \rangle$ from Au+Au collisions at different collision energies.

$\langle N_{part} \rangle$ dependence of efficiency corrected $\mu / \sigma^2$ of net-charge distributions for Au+Au collisions at different collision energies.

$\langle N_{part} \rangle$ dependence of efficiency corrected $S \sigma$ of net-charge distributions for Au+Au collisions at different collision energies.

$\langle N_{part} \rangle$ dependence of efficiency corrected $\kappa \sigma^2$ of net-charge distributions for Au+Au collisions at different collision energies.

$\langle N_{part} \rangle$ dependence of efficiency corrected $S \sigma^3 / \mu$ of net-charge distributions for Au+Au collisions at different collision energies.

The energy dependence of efficiency corrected $\mu / \sigma^2$, $S \sigma$, $\kappa \sigma^2$, and $S \sigma^3 / \mu$ of netcharge distributions for central (0%–5%) Au+Au collisions.

The energy dependence of the chemical freeze-out parameter $\mu_B$.

Away-side jet widths from a Gaussian fit by $h^{\pm}$ partner momentum for various $\pi^0$ trigger momenta in Au+Au collisions.

Away-side $I_{AA}$ for the entire away-side $|\Delta \phi - \pi| < \pi /2$ selection vs. $h^{\pm}$ partner momentum for various $\pi^0$ trigger momenta. A point-to-point uncorrelated 6% normalization uncertainty (mainly due to efficiency corrections) applies to all measurements.

Away-side $I_{AA}$ for the entire away-side $|\Delta \phi - \pi| < \pi /2$ selection vs. $h^{\pm}$ partner momentum for various $\pi^0$ trigger momenta. A point-to-point uncorrelated 6% normalization uncertainty (mainly due to efficiency corrections) applies to all measurements.

Away-side $I_{AA}$ for a narrow “head” $|\Delta \phi - \pi| < \pi /6$ selection vs. $h^{\pm}$ partner momentum for various $\pi^0$ trigger momenta. A point-to-point uncorrelated 6% normalization uncertainty (mainly due to efficiency corrections) applies to all measurements.

Away-side $I_{AA}$ for a narrow “head” $|\Delta \phi - \pi| < \pi /6$ selection vs. $h^{\pm}$ partner momentum for various $\pi^0$ trigger momenta. A point-to-point uncorrelated 6% normalization uncertainty (mainly due to efficiency corrections) applies to all measurements.

Supplemental near-side jet widths from a Gaussian fit by $h^{\pm}$ partner momentum for various $\pi^0$ trigger momenta in $p+p$ collisions.

Supplemental near-side jet widths from a Gaussian fit by $h^{\pm}$ partner momentum for various $\pi^0$ trigger momenta in Au+Au collisions.

Supplemental near-side jet widths from a Gaussian fit by $h^{\pm}$ partner momentum for various $\pi^0$ trigger momenta in Au+Au collisions.

Supplemental away-side $I_{AA}$ for $|\Delta \phi - \pi| < \pi /3$ selection vs. $h^{\pm}$ partner momentum for various $\pi^0$ trigger momenta. A point-to-point uncorrelated 6% normalization uncertainty (mainly due to efficiency corrections) applies to all measurements.

Supplemental away-side $I_{AA}$ for $|\Delta \phi - \pi| < \pi /3$ selection vs. $h^{\pm}$ partner momentum for various $\pi^0$ trigger momenta. A point-to-point uncorrelated 6% normalization uncertainty (mainly due to efficiency corrections) applies to all measurements.

$p+p$ yield vs $p_T^a$ supplemental tables.

Au+Au yield vs $p_T^a$ supplemental tables.

Au+Au yield vs $p_T^a$ supplemental tables.

$z_{g}$ for Matched Trigger jets in pp$+$AuAu Data anti-kT R$=$0.4

$z_{g}$ for HardCore Recoil jets in AuAu Data anti-kT R$=$0.4

$z_{g}$ for HardCore Recoil jets in pp$+$AuAu Data anti-kT R$=$0.4

$z_{g}$ for Matched Recoil jets in AuAu Data anti-kT R$=$0.4

$z_{g}$ for Matched Recoil jets in pp$+$AuAu Data anti-kT R$=$0.4

$R_{g}$ for HardCore Trigger jets in AuAu Data anti-kT R$=$0.4

$R_{g}$ for HardCore Trigger jets in pp$+$AuAu Data anti-kT R$=$0.4

$R_{g}$ for Matched Trigger jets in AuAu Data anti-kT R$=$0.4

$R_{g}$ for Matched Trigger jets in pp$+$AuAu Data anti-kT R$=$0.4

$R_{g}$ for HardCore Recoil jets in AuAu Data anti-kT R$=$0.4

$R_{g}$ for HardCore Recoil jets in pp$+$AuAu Data anti-kT R$=$0.4

$R_{g}$ for Matched Recoil jets in AuAu Data anti-kT R$=$0.4

$R_{g}$ for Matched Recoil jets in pp$+$AuAu Data anti-kT R$=$0.4

Fully corrected two subjet $z_{sj}$ for $15 < p_{\rm{T, jet}} < 20$ GeV/c, R$=$0.4 anti-kT and R$=$0.1 subjets

Fully corrected two subjet $z_{sj}$ for $20 < p_{\rm{T, jet}} < 25$ GeV/c, R$=$0.4 anti-kT and R$=$0.1 subjets

Fully corrected two subjet $z_{sj}$ for $25 < p_{\rm{T, jet}} < 30$ GeV/c, R$=$0.4 anti-kT and R$=$0.1 subjets

Fully corrected two subjet $z_{sj}$ for $30 < p_{\rm{T, jet}} < 40$ GeV/c, R$=$0.4 anti-kT and R$=$0.1 subjets

Fully corrected two subjet $\theta_{sj}$ for $15 < p_{\rm{T, jet}} < 20$ GeV/c, R$=$0.4 anti-kT and R$=$0.1 subjets

Fully corrected two subjet $\theta_{sj}$ for $20 < p_{\rm{T, jet}} < 25$ GeV/c, R$=$0.4 anti-kT and R$=$0.1 subjets

Fully corrected two subjet $\theta_{sj}$ for $25 < p_{\rm{T, jet}} < 30$ GeV/c, R$=$0.4 anti-kT and R$=$0.1 subjets

Fully corrected two subjet $\theta_{sj}$ for $30 < p_{\rm{T, jet}} < 40$ GeV/c, R$=$0.4 anti-kT and R$=$0.1 subjets

$z_{sj}$ for HardCore Trigger jets in AuAu Data anti-kT R$=$0.4

$z_{sj}$ for HardCore Trigger jets in pp$+$AuAu Data anti-kT R$=$0.4

$z_{sj}$ for Matched Trigger jets in AuAu Data anti-kT R$=$0.4

$z_{sj}$ for Matched Trigger jets in pp$+$AuAu Data anti-kT R$=$0.4

$z_{sj}$ for HardCore Recoil jets in AuAu Data anti-kT R$=$0.4

$z_{sj}$ for HardCore Recoil jets in pp$+$AuAu Data anti-kT R$=$0.4

$z_{sj}$ for Matched Recoil jets in AuAu Data anti-kT R$=$0.4

$z_{sj}$ for Matched Recoil jets in pp$+$AuAu Data anti-kT R$=$0.4

$\theta_{sj}$ for HardCore Trigger jets in AuAu Data anti-kT R$=$0.4

$\theta_{sj}$ for HardCore Trigger jets in pp$+$AuAu Data anti-kT R$=$0.4

$\theta_{sj}$ for Matched Trigger jets in AuAu Data anti-kT R$=$0.4

$\theta_{sj}$ for Matched Trigger jets in pp$+$AuAu Data anti-kT R$=$0.4

$\theta_{sj}$ for HardCore Recoil jets in AuAu Data anti-kT R$=$0.4

$\theta_{sj}$ for HardCore Recoil jets in pp$+$AuAu Data anti-kT R$=$0.4

$\theta_{sj}$ for Matched Recoil jets in AuAu Data anti-kT R$=$0.4

$\theta_{sj}$ for Matched Recoil jets in pp$+$AuAu Data anti-kT R$=$0.4

$A_{J}$ for HardCore R=0.4 anti-kT jets with $\theta_{sj} > 0.1$ in Au$+$Au collisions

$A_{J}$ for HardCore R=0.4 anti-kT jets with $\theta_{sj} > 0.1$ in pp$+$AuAu reference

$A_{J}$ for HardCore R=0.4 anti-kT jets with $0.1 < \theta_{sj} < 0.2$ in Au$+$Au collisions

$A_{J}$ for HardCore R=0.4 anti-kT jets with $0.1 < \theta_{sj} < 0.2$ in pp$+$AuAu reference

$A_{J}$ for HardCore R=0.4 anti-kT jets with $0.2 < \theta_{sj} < 0.3$ in Au$+$Au collisions

$A_{J}$ for HardCore R=0.4 anti-kT jets with $0.2 < \theta_{sj} < 0.3$ in pp$+$AuAu reference

$A_{J}$ for Matched R=0.4 anti-kT jets with $\theta_{sj} > 0.1$ in Au$+$Au collisions

$A_{J}$ for Matched R=0.4 anti-kT jets with $\theta_{sj} > 0.1$ in pp$+$AuAu reference

$A_{J}$ for Matched R=0.4 anti-kT jets with $0.1 < \theta_{sj} < 0.2$ in Au$+$Au collisions

$A_{J}$ for Matched R=0.4 anti-kT jets with $0.1 < \theta_{sj} < 0.2$ in pp$+$AuAu reference

$A_{J}$ for Matched R=0.4 anti-kT jets with $0.2 < \theta_{sj} < 0.3$ in Au$+$Au collisions

$A_{J}$ for Matched R=0.4 anti-kT jets with $0.2 < \theta_{sj} < 0.3$ in pp$+$AuAu reference

tracking efficiency in pp collisions as a function of track momentum

tracking efficiency in AuAu collisions as a function of track momentum

absolutely normalized dijet yields scaled by 1/<TAA> in 20-40% central PbPb collisions

absolutely normalized dijet yields scaled by 1/<TAA> in 40-60% central PbPb collisions

absolutely normalized dijet yields scaled by 1/<TAA> in 60-80% central PbPb collisions

self normalized dijets from pp collisions

self normalized dijet distributions in 0-10% central PbPb collisions

self normalized dijet distributions in 10-20% central PbPb collisions

self normalized dijet distributions in 20-40% central PbPb collisions

self normalized dijet distributions in 40-60% central PbPb collisions

self normalized dijet distributions in 60-80% central PbPb collisions

leading jet RAA^pair in 0-10% central PbPb collisions

subleading jet RAA^pair in 0-10% central PbPb collisions

leading jet RAA^pair in 10-20% central PbPb collisions

subleading jet RAA^pair in 10-20% central PbPb collisions

leading jet RAA^pair in 20-40% central PbPb collisions

subleading jet RAA^pair in 20-40% central PbPb collisions

leading jet RAA^pair in 40-60% central PbPb collisions

subleading jet RAA^pair in 40-60% central PbPb collisions

leading jet RAA^pair in 60-80% central PbPb collisions

subleading jet RAA^pair in 60-80% central PbPb collisions

ratio of subleading jet RAA^pair to leading jet RAA^pair in PbPb collisions

\Delta v_1 between proton and anti-proton vs rapidity in Au+Au 200 GeV 50-80% centrality.

\Delta v_1 between proton and anti-proton vs rapidity in Zr+Zr and Ru+Ru (combined) 200 GeV 50-80% centrality.

\Delta v_1 between proton and anti-proton vs rapidity in Au+Au 27 GeV 50-80% centrality.

\Delta dv1/dy as a function of centrality in Au+Au 200 GeV.

\Delta dv1/dy as a function of centrality in Ru+Ru and Zr+Zr 200 GeV

\Delta dv1/dy as a function of centrality in Au+Au 27 GeV.

Total efficiencies for the selection of signal events after categorization and application of the tighter photon $E_{\mathrm{T}}^{\gamma}$ selection used to optimize the signal significance spin-1 $q\bar{q'}\to X^{\pm} \to W^{\pm}\gamma$. In addition to the total efficiency, contributions to the signal selection from each of the separate event categories are shown. The efficiencies calculated from MC samples with $W/Z$ hadronic decays are shown as the points on each curve. The line presents interpolated results.

Total efficiencies for the selection of signal events after categorization and application of the tighter photon $E_{\mathrm{T}}^{\gamma}$ selection used to optimize the signal significance spin-2 $q\bar{q}\to X^0 \to Z\gamma$. In addition to the total efficiency, contributions to the signal selection from each of the separate event categories are shown. The efficiencies calculated from MC samples with $W/Z$ hadronic decays are shown as the points on each curve. The line presents interpolated results.

Total efficiencies for the selection of signal events after categorization and application of the tighter photon $E_{\mathrm{T}}^{\gamma}$ selection used to optimize the signal significance spin-2 $gg\to X^0 \to Z\gamma$. In addition to the total efficiency, contributions to the signal selection from each of the separate event categories are shown. The efficiencies calculated from MC samples with $W/Z$ hadronic decays are shown as the points on each curve. The line presents interpolated results.

The $m_{J\gamma}$ distributions of data events selected for the spin-1 $q\bar{q'}\to X^{\pm} \to W^{\pm}\gamma$ search in the D2 category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-1 $q\bar{q'}\to X^{\pm} \to W^{\pm}\gamma$ search in the WMASS category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-0 $gg \to X^{0} \to Z\gamma$ search in the BTAG category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-0 $gg \to X^{0} \to Z\gamma$ search in the D2 category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-0 $gg \to X^{0} \to Z\gamma$ search in the ZMASS category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-2 $gg \to X^{0} \to Z\gamma$ search in the BTAG category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-2 $gg \to X^{0} \to Z\gamma$ search in the D2 category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-2 $gg \to X^{0} \to Z\gamma$ search in the ZMASS category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-2 $q\bar{q} \to X^{0} \to Z\gamma$ search in the BTAG category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-2 $q\bar{q} \to X^{0} \to Z\gamma$ search in the D2 category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The $m_{J\gamma}$ distributions of data events selected for the spin-2 $q\bar{q} \to X^{0} \to Z\gamma$ search in the ZMASS category. The background-only fit function shape is shown as the solid curve overlaid with a shaded band corresponding to statistical uncertainties in background parameters. Various signal shapes with cross-sections corresponding to expected limits obtained in this analysis are shown as dashed lines.

The 95% CL upper limits on $\sigma(pp\to X)\times B(X\to W/Z\gamma)$ as a function of $m_{X}$ for spin-0 $gg \to X^{0} \to Z\gamma$. The observed limits are shown as a solid black line and the expected ones are shown as a dashed line with the 1$\sigma$ (2$\sigma$) uncertainty band presented as the green (yellow) band.

The 95% CL upper limits on $\sigma(pp\to X)\times B(X\to W/Z\gamma)$ as a function of $m_{X}$ for spin-2 $gg \to X^{0} \to Z\gamma$. The observed limits are shown as a solid black line and the expected ones are shown as a dashed line with the 1$\sigma$ (2$\sigma$) uncertainty band presented as the green (yellow) band.

The 95% CL upper limits on $\sigma(pp\to X)\times B(X\to W/Z\gamma)$ as a function of $m_{X}$ for spin-2 $q\bar{q} \to X^{0} \to Z\gamma$. The observed limits are shown as a solid black line and the expected ones are shown as a dashed line with the 1$\sigma$ (2$\sigma$) uncertainty band presented as the green (yellow) band.

The 95% CL upper limits on $\sigma(pp\to X)\times B(X\to W/Z\gamma)$ as a function of $m_{X}$ for spin-1 $q\bar{q'}\to X^{\pm} \to W^{\pm}\gamma$. The observed limits are shown as a solid black line and the expected ones are shown as a dashed line with the 1$\sigma$ (2$\sigma$) uncertainty band presented as the green (yellow) band.

The jet mass distribution of large-$R$ jets originating from the hadronic decay of $W$ and $Z$ bosons produced from the decay of BSM bosons with mass $m_X = 2000$ GeV. The decays simulated are for the production models $q\bar{q'}\to X^{\pm} \to W^{\pm}\gamma$ with a spin-1 resonance $X^{\pm}$ and $gg\to X^0 \to Z\gamma$ with a spin-0 resonance $X^{0}$.

The photon $E_{T}$ selection criteria as a function of resonance mass for D2 categories. Selections for four different signal channels are shown together.

The photon $E_{T}$ selection criteria as a function of resonance mass for ZMASS or WMASS categories. Selections for four different signal channels are shown together.

Forward-rapidity $J/\psi$ —> $\mu^+\mu^-$ $d$+Au Nuclear Dependence at $\sqrt{s}$ = 200 GeV. (sys. A, B systematics are relative, i.e. they multiply the $R_{dAu}$ value. The sysA uncertainty includes statistical uncertainties as well as point-to-point uncorrelated systematics, sys. B represents uncertainties that are correlated from point to point, and sys. C represents uncertainties in the overall normalization.)

Mid rapidity $d$+Au —> $e^+e^-$ $J/\psi$ $R_{dAu}$ at $\sqrt{s}$=200 GeV. (Sys. A, B systematics are absolute, i.e. they add/subtract directly from $R_{dAu}$. The sys. A uncertainty includes statistical uncertainties as well as point-to-point uncorrelated systematics, sys. B represents uncertainties that are correlated from point to point, and sys. C represents uncertainties in the overall normalization.)

$d$+Au vs centrality vs $y$. (All uncertainties are absolute. The sys. A uncertainty includes both the statistical uncertainty and the point-to-point uncorrelated systematic, sys. B represents uncertainties that are correlated from point to point, and sys. C represents uncertainties in the overall normalization.)

Forward-rapidity $J/\psi$ —> $\mu^+\mu^-$ $d$+Au Nuclear Dependence at $\sqrt{s}$ = 200 GeV. (sys. A, B systematics are relative, i.e. they multiply the $R_{dAu}$ value. The sysA uncertainty includes statistical uncertainties as well as point-to-point uncorrelated systematics, sys. B represents uncertainties that are correlated from point to point, and sys. C represents uncertainties in the overall normalization.)

Mid rapidity $d$+Au —> $e^+e^-$ $J/\psi$ $R_{dAu}$ at $\sqrt{s}$=200 GeV. (Sys. A, B systematics are absolute, i.e. they add/subtract directly from $R_{dAu}$. The sys. A uncertainty includes statistical uncertainties as well as point-to-point uncorrelated systematics, sys. B represents uncertainties that are correlated from point to point, and sys. C represents uncertainties in the overall normalization.)

Forward-rapidity $J/\psi$ —> $\mu^+\mu^-$ $d$+Au Nuclear Dependence at $\sqrt{s}$ = 200 GeV. (sys. A, B systematics are relative, i.e. they multiply the $R_{CP}$ value. The sys. A uncertainty includes statistical uncertainties as well as point-to-point uncorrelated systematics, sys. B represents uncertainties that are correlated from point to point, and sys. C represents uncertainties in the overall normalization.)

Mid rapidity $d$+Au —> $e^+e^-$ $J/\psi$ $R_{CP}$ at $\sqrt{s}$=200 GeV. (Sys. A, B systematics are absolute, i.e. they add/subtract directly from $R_{CP}$. The sys. A uncertainty includes statistical uncertainties as well as point-to-point uncorrelated systematics, sys. B represents uncertainties that are correlated from point to point, and sys. C represents uncertainties in the overall normalization.)

Events in bins of the $E_{miss}^T$ in the SR.

Parametrization of the $A_{X}$ factor, defined as the fraction of diphoton resonances satisfying the fiducial acceptance at the particle level, as function of $m_{X}$ obtained from the narrow width simulated signal samples produced in gluon fusion.

The correction factor, $C_{X}$, defined as the ratio between the number of reconstructed signal events passing the analysis cuts and the number of signal events at the particle level generated within the fiducial volume, and acceptance correction factor, $A_{X}$, defined as the fraction of diphoton resonances satisfying the fiducial acceptance at the particle level. Both are computed for NWA spin-0 models as a function of $m_{X}$.

Effect of event selections on a scalar MC signal sample generated for $m_{X}$ = 15 GeV and on the data. For the MC sample, the efficiencies are shown after applying event weights and a truth level filter that requires two photons with $p^{\gamma\gamma}_{T}>40$ GeV; for the data, the absolute yields are shown. The initial yields for data include a trigger preselection that is the OR of a list of single photon and diphoton triggers. The "2 $loose$ photons" step includes the kinematic acceptance cuts.

Parameterization of the Double Sided Crystal Ball function parameters describing the scalar mass resolution model as a function of $m_{X}$ [GeV].

The reconstruction efficiency for caloDPJs produced by the decay of &gamma;_d into e⁺e^- or qq&#772;. The reconstruction efficiency is shown for &gamma;_d with 0&lt;|&eta;|&lt;1.1 as a function of their transverse momentum in events where the &gamma;_d L_xy is between 2&nbsp;m and 4&nbsp;m.

Observed 95% CL upper limits on the branching ratio (B) for the decay H&rarr;&nbsp;2&gamma;_d+X and ggF selection, for different &gamma;_d masses and a 125&nbsp;GeV Higgs boson. The limits are shown for the SR_2&mu;^ggF, assuming an FRVZ signal model. The hatched band denotes the region in which the branching ratio of H&rarr;&nbsp;2&gamma;_d+X is larger than unity.

Observed 95% CL upper limits on the branching ratio (B) for the decay H&rarr;&nbsp;2&gamma;_d+X and ggF selection, for different &gamma;_d masses and a 125&nbsp;GeV Higgs boson. The limits are shown for the SR_c+&mu;^ggF, assuming an FRVZ signal model. The hatched band denotes the region in which the branching ratio of H&rarr;&nbsp;2&gamma;_d+X is larger than unity.

Observed 95% CL upper limits on the branching ratio (B) for the decay H&rarr;&nbsp;2&gamma;_d+X and ggF selection, for different &gamma;_d masses and a 125&nbsp;GeV Higgs boson. The limits are shown for the SR_2c^ggF, assuming an FRVZ signal model. The hatched band denotes the region in which the branching ratio of H&rarr;&nbsp;2&gamma;_d+X is larger than unity.

Observed 95% CL upper limits on the branching ratio (B) for the decay H&rarr;&nbsp;2&gamma;_d+X and WH selection, for different &gamma;_d masses and a 125&nbsp;GeV Higgs boson. The limits are shown for the SR_c^WH, assuming an FRVZ signal model. The hatched band denotes the region in which the branching ratio of H&rarr;&nbsp;2&gamma;_d+X is larger than unity.

Observed 95% CL upper limits on the branching ratio (B) for the decay H&rarr;&nbsp;2&gamma;_d+X and WH selection, for different &gamma;_d masses and a 125&nbsp;GeV Higgs boson. The limits are shown for the SR_c+&mu;^WH, assuming an FRVZ signal model. The hatched band denotes the region in which the branching ratio of H&rarr;&nbsp;2&gamma;_d+X is larger than unity.

Observed 95% CL upper limits on the branching ratio (B) for the decay H&rarr;&nbsp;2&gamma;_d+X and WH selection, for different &gamma;_d masses and a 125&nbsp;GeV Higgs boson. The limits are shown for the SR_2c^WH, assuming an FRVZ signal model. The hatched band denotes the region in which the branching ratio of H&rarr;&nbsp;2&gamma;_d+X is larger than unity.

Observed 95% CL upper limits on the branching ratio (B) for the decay H&rarr;&nbsp;2&gamma;_d and ggF selection, for different &gamma;_d masses and a 125&nbsp;GeV Higgs boson. The limits are shown for the SR_2&mu;^ggF, assuming a HAHM signal model. The hatched band denotes the region in which the branching ratio of H&rarr;&nbsp;2&gamma;_d is larger than unity.

Observed 95% CL upper limits on the branching ratio (B) for the decay H&rarr;&nbsp;2&gamma;_d and ggF selection, for different &gamma;_d masses and a 125&nbsp;GeV Higgs boson. The limits are shown for the SR_c+&mu;^ggF, assuming a HAHM signal model. The hatched band denotes the region in which the branching ratio of H&rarr;&nbsp;2&gamma;_d is larger than unity.

Observed 95% CL upper limits on the branching ratio (B) for the decay H&rarr;&nbsp;2&gamma;_d and ggF selection, for different &gamma;_d masses and a 125&nbsp;GeV Higgs boson. The limits are shown for the SR_2c^ggF, assuming a HAHM signal model. The hatched band denotes the region in which the branching ratio of H&rarr;&nbsp;2&gamma;_d is larger than unity.

The CalRatio 2015--2016 trigger efficiency for the decay of &gamma;_d in e⁺e^- or qq&#772;. The trigger efficiency is shown for &gamma;_d with 0&lt;|&eta;|&lt;1.1 as function of the transverse decay length L_xy.

The CalRatio 2015--2016 trigger efficiency for the decay of &gamma;_d in e⁺e^- or qq&#772;. The trigger efficiency is shown for &gamma;_d with 0&lt;|&eta;|&lt;1.1 as function of the transverse momentum, in events where the &gamma;_d L_xy is between 2&nbsp;m and 4&nbsp;m.

The CalRatio 2017--2018 trigger efficiency for the decay of &gamma;_d in e⁺e^- or qq&#772;. The trigger efficiency is shown for &gamma;_d with 0&lt;|&eta;|&lt;1.1 as function of the transverse decay length L_xy.

The CalRatio 2017--2018 trigger efficiency for the decay of &gamma;_d in e⁺e^- or qq&#772;. The trigger efficiency is shown for &gamma;_d with 0&lt;|&eta;|&lt;1.1 as function of the transverse momentum, in events where the &gamma;_d L_xy is between 2&nbsp;m and 4&nbsp;m.

The narrow-scan 2015--2016 trigger efficiency for the decay of &gamma;_d in &mu;⁺&mu;^-. The trigger efficiency is shown for &gamma;_d with 0&lt;|&eta;|&lt;1.1 as function of the transverse decay length L_xy.

The narrow-scan 2015--2016 trigger efficiency for the decay of &gamma;_d in &mu;⁺&mu;^-. The trigger efficiency is shown for &gamma;_d with 0&lt;|&eta;|&lt;1.1 as function of the transverse momentum, in events where the &gamma;_d L_xy is below 6&nbsp;m.

The narrow-scan 2017--2018 trigger efficiency for the decay of &gamma;_d in &mu;⁺&mu;^-. The trigger efficiency is shown for &gamma;_d with 0&lt;|&eta;|&lt;1.1 as function of the transverse decay length L_xy.

The narrow-scan 2017--2018 trigger efficiency for the decay of &gamma;_d in &mu;⁺&mu;^-. The trigger efficiency is shown for &gamma;_d with 0&lt;|&eta;|&lt;1.1 as function of the transverse momentum, in events where the &gamma;_d L_xy is below 6&nbsp;m.

Efficiency of the cosmic-ray tagger as function of the &gamma;_d transverse decay length. The efficiency is calculated accepting the &mu;DPJs for which the cosmic-ray tagger score is &gt; 0.2 for each associated MS-only track.

Efficiency of the cosmic-ray tagger as function of the &gamma;_d transverse momentum. The efficiency is calculated accepting the &mu;DPJs for which the cosmic-ray tagger score is &gt; 0.2 for each associated MS-only track.

Efficiency of the QCD tagger as function of the &gamma;_d transverse decay length. The efficiency is calculated accepting the caloDPJs for which the QCD tagger score is &gt; 0.5.

Efficiency of the QCD tagger as function of the &gamma;_d transverse momentum. The efficiency is calculated accepting the caloDPJs for which the QCD tagger score is &gt; 0.5.

The extrapolated acceptance times efficiencies as a function of the &gamma;_d mean proper lifetime c&tau; for the FRVZ pp &rarr; H&rarr;&nbsp;2&gamma;_d+X process and a SM Higgs boson for the SR_2&mu;^ggF.

The extrapolated acceptance times efficiencies as a function of the &gamma;_d mean proper lifetime c&tau; for the FRVZ pp &rarr; H&rarr;&nbsp;2&gamma;_d+X process and a SM Higgs boson for the SR_c+&mu;^ggF.

The extrapolated acceptance times efficiencies as a function of the &gamma;_d mean proper lifetime c&tau; for the FRVZ pp &rarr; H&rarr;&nbsp;2&gamma;_d+X process and a SM Higgs boson for the SR_2c^ggF.

The extrapolated acceptance times efficiencies as a function of the &gamma;_d mean proper lifetime c&tau; for the FRVZ pp &rarr; H&rarr;&nbsp;2&gamma;_d+X process and a SM Higgs boson for the SR_c^WH.

The extrapolated acceptance times efficiencies as a function of the &gamma;_d mean proper lifetime c&tau; for the FRVZ pp &rarr; H&rarr;&nbsp;2&gamma;_d+X process and a SM Higgs boson for the SR_2c^WH.

The extrapolated acceptance times efficiencies as a function of the &gamma;_d mean proper lifetime c&tau; for the FRVZ pp &rarr; H&rarr;&nbsp;2&gamma;_d+X process and a SM Higgs boson for the SR_c+&mu;^WH.

The extrapolated acceptance times efficiencies for SR_2&mu;^ggF, as a function of the &gamma;_d mean proper lifetime c&tau; for the FRVZ pp &rarr; H&rarr;&nbsp;2&gamma;_d+X process and a Higgs-like scalar with a mass of 800&nbsp;GeV. The coloured bands represent the MC statistical uncertainties.

The extrapolated acceptance times efficiencies for SR_c+&mu;^ggF, as a function of the &gamma;_d mean proper lifetime c&tau; for the FRVZ pp &rarr; H&rarr;&nbsp;2&gamma;_d+X process and a Higgs-like scalar with a mass of 800&nbsp;GeV. The coloured bands represent the MC statistical uncertainties.

The extrapolated acceptance times efficiencies for SR_2c;^ggF, as a function of the &gamma;_d mean proper lifetime c&tau; for the FRVZ pp &rarr; H&rarr;&nbsp;2&gamma;_d+X process and a Higgs-like scalar with a mass of 800&nbsp;GeV. The coloured bands represent the MC statistical uncertainties.

The extrapolated acceptance times efficiencies for SR_2&mu;^ggF, as a function of the &gamma;_d mean proper lifetime c&tau; for the HAHM H&rarr;&nbsp;2&gamma;_d+X process and a SM Higgs boson for the ggF channel.

The extrapolated acceptance times efficiencies for SR_c+&mu;^ggF, as a function of the &gamma;_d mean proper lifetime c&tau; for the HAHM H&rarr;&nbsp;2&gamma;_d+X process and a SM Higgs boson for the ggF channel.

The extrapolated acceptance times efficiencies for SR_2c;^ggF, as a function of the &gamma;_d mean proper lifetime c&tau; for the HAHM H&rarr;&nbsp;2&gamma;_d+X process and a SM Higgs boson for the ggF channel.

The extrapolated acceptance times efficiencies for SR_2&mu;^ggF, as a function of the &gamma;_d mean proper lifetime c&tau; for the FRVZ pp &rarr; H&rarr;&nbsp;4&gamma;_d+X process and a Higgs-like scalar with a mass of 800&nbsp;GeV. The coloured bands represent the MC statistical uncertainties.

The extrapolated acceptance times efficiencies for SR_c+&mu;^ggF, as a function of the &gamma;_d mean proper lifetime c&tau; for the FRVZ pp &rarr; H&rarr;&nbsp;4&gamma;_d+X process and a Higgs-like scalar with a mass of 800&nbsp;GeV. The coloured bands represent the MC statistical uncertainties.

The extrapolated acceptance times efficiencies for SR_2c;^ggF, as a function of the &gamma;_d mean proper lifetime c&tau; for the FRVZ pp &rarr; H&rarr;&nbsp;4&gamma;_d+X process and a Higgs-like scalar with a mass of 800&nbsp;GeV. The coloured bands represent the MC statistical uncertainties.

The extrapolated acceptance times efficiencies for SR_2&mu;^ggF, as a function of the &gamma;_d mean proper lifetime c&tau; for the FRVZ pp &rarr; H&rarr;&nbsp;4&gamma;_d+X process and a SM Higgs boson. The coloured bands represent the MC statistical uncertainties.

The extrapolated acceptance times efficiencies for SR_c+&mu;^ggF, as a function of the &gamma;_d mean proper lifetime c&tau; for the FRVZ pp &rarr; H&rarr;&nbsp;4&gamma;_d+X process and a SM Higgs boson. The coloured bands represent the MC statistical uncertainties.

The extrapolated acceptance times efficiencies for SR_2c;^ggF, as a function of the &gamma;_d mean proper lifetime c&tau; for the FRVZ pp &rarr; H&rarr;&nbsp;4&gamma;_d+X process and a SM Higgs boson. The coloured bands represent the MC statistical uncertainties.

Observed limits at the 95% CL on the branching ratio (B), obtained from the SR_2&mu;^ggF, for the process pp &rarr; H&rarr;&nbsp;4&gamma;_d+X, for different &gamma;_d masses and a 125 GeV Higgs boson. The limits are shown assuming an FRVZ signal model.

Observed limits at the 95% CL on the branching ratio (B), obtained from the SR_c+&mu;^ggF, for the process pp &rarr; H&rarr;&nbsp;4&gamma;_d+X, for different &gamma;_d masses and a 125 GeV Higgs boson. The limits are shown assuming an FRVZ signal model.

Observed limits at the 95% CL on the branching ratio (B), obtained from the SR_2c^ggF, for the process pp &rarr; H&rarr;&nbsp;4&gamma;_d+X, for different &gamma;_d masses and a 125 GeV Higgs boson. The limits are shown assuming an FRVZ signal model.

Observed limits at the 95% CL on the cross section, obtained from the SR_2&mu;^ggF, for the process pp &rarr; H&rarr;&nbsp;2&gamma;_d+X, for different &gamma;_d masses and a Higgs-like scalar with a mass of 800&nbsp;GeV. The limits are shown assuming an FRVZ signal model.

Observed limits at the 95% CL on the cross section, obtained from the SR_c+&mu;^ggF, for the process pp &rarr; H&rarr;&nbsp;2&gamma;_d+X, for different &gamma;_d masses and a Higgs-like scalar with a mass of 800&nbsp;GeV. The limits are shown assuming an FRVZ signal model.

Observed limits at the 95% CL on the cross section, obtained from the SR_2c^ggF, for the process pp &rarr; H&rarr;&nbsp;2&gamma;_d+X and ggF selection, for different &gamma;_d masses and a Higgs-like scalar with a mass of 800&nbsp;GeV. The limits are shown assuming an FRVZ signal model.

Observed limits at the 95% CL on the cross section for the process pp &rarr; H&rarr;&nbsp;4&gamma;_d+X and ggF selection, for different &gamma;_d masses and a Higgs-like scalar with a mass of 800&nbsp;GeV. The limits are shown for the 2&mu; search channel, assuming an FRVZ signal model.

Observed limits at the 95% CL on the cross section for the process pp &rarr; H&rarr;&nbsp;4&gamma;_d+X and ggF selection, for different &gamma;_d masses and a Higgs-like scalar with a mass of 800&nbsp;GeV. The limits are shown for the c+&mu; search channel, assuming an FRVZ signal model.

Observed limits at the 95% CL on the cross section for the process pp &rarr; H&rarr;&nbsp;4&gamma;_d+X and ggF selection, for different &gamma;_d masses and a Higgs-like scalar with a mass of 800&nbsp;GeV. The limits are shown for the 2c search channel, assuming an FRVZ signal model.

Background-only fit to the observed numbers of events (black points) as a function of the diphoton invariant mass $m_{\gamma\gamma}$.

Expected and observed 95% CL upper limits on the signal fiducial cross section $\sigma_{\textrm{fid}}$, assuming 100% branching ratio for ALP decay into two photons, as a function of the hypothetical ALP mass $m_{\textrm{X}}$. The 1$\sigma$ and 2$\sigma$ confidence intervals are shown by the coloured bands.

Expected and observed 95% CL upper limits on the the ALP coupling constant, assuming 100% branching ratio for ALP decay into two photons, as a function of the hypothetical ALP mass $m_{\textrm{X}}$. The 1$\sigma$ and 2$\sigma$ confidence intervals are shown by the coloured bands. Contours of the ALP natural width $\Gamma$ are illustrated by the smooth blue solid lines.

Fiducial region acceptance of signal events for coupling constant f$^{-1}=0.05$ TeV$^{-1}$. EL, SD, and DD stand for the exclusive, single-dissociative, and double-dissociative processes. The acceptances are parameterised as $p_0 + p^{i}_{1} e^{p_{2}^{i} m_{\textrm{X}}} + p_{3}^{i} e^{p_{4}^{i} m_{\textrm{X}}} $ , where $i$ specifies the process type, exclusive, single-dissociative, or double-dissociative.

Correction factors $A_{\textrm{X}}/A_{0}$ and $C_{\textrm{X}}$ for the exclusive signal events with ALP coupling constant f$^{-1}$=0.05 TeV$^{-1}$, where acceptance $A_{\textrm{X}}$ and efficiency $C_{\textrm{X}}$ are defined in arXiv:2102.13405. $A_{\textrm{X}}$/$A_{0}$ and $C_{\textrm{X}}$ are parameterised as $p_0 + p^{i}_{1} e^{p_{2}^{i} m_{\textrm{X}}} + p_{3}^{i} e^{p_{4}^{i} m_{\textrm{X}}}$, where $i$ specifies the process type.

Simulated ALP signal cutflow for $m_{\rm{X}}$=300 GeV and $m_{\rm{X}}$=1200 GeV. The yields are normalized to a luminosity of 14.6 fb$^{-1}$. Selections are applied for each side after the selection on acoplanarity. The union of the sets of selected events is taken finally.

Comparison between data and background prediction before the fit (Pre-Fit) for the mass of the FCNC top-quark candidate in SR1. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The four FCNC LH signals are also shown separately, normalized to five times the cross-section corresponding to the most stringent observed branching ratio limits. The first (last) bin in all distributions includes the underflow (overflow). The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction before the fit (Pre-Fit) for the mass of the SM top-quark candidate in SR2. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The four FCNC LH signals are also shown separately, normalized to five times the cross-section corresponding to the most stringent observed branching ratio limits. The first (last) bin in all distributions includes the underflow (overflow). The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction before the fit (Pre-Fit) for the transverse momentum of the Z boson candidate in SR2. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The four FCNC LH signals are also shown separately, normalized to five times the cross-section corresponding to the most stringent observed branching ratio limits. The first (last) bin in all distributions includes the underflow (overflow). The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZu LH coupling extraction. The distribution is for the $D_{1}$ discriminant in the mass sideband CR1. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZu LH signals are also separately shown, normalized to 500 times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZu LH coupling extraction. The distribution is for the $D_{2}^{u}$ discriminant in the mass sideband CR2. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZu LH signals are also separately shown, normalized to 500 times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZu LH coupling extraction. The distribution is for the $D_{1}$ discriminant in SR1. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZu LH signals are also separately shown, normalized to 50 times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZu LH coupling extraction. The distribution is for the $D_{2}^{u}$ discriminant in SR2. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZu LH signals are also separately shown, normalized to 50 times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZc LH coupling extraction. The distribution is for the $D_{1}$ discriminant in the mass sideband CR1. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZc LH signals are also separately shown, normalized to 500 times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZc LH coupling extraction. The distribution is for the $D_{2}^{c}$ discriminant in the mass sideband CR2. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZc LH signals are also separately shown, normalized to 500 times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZc LH coupling extraction. The distribution is for the $D_{1}$ discriminant in SR1. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZc LH signals are also separately shown, normalized to 50 times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZc LH coupling extraction. The distribution is for the $D_{2}^{c}$ discriminant in SR2. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZc LH signals are also separately shown, normalized to 50 times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZu LH coupling extraction. The distribution is for the leading lepton $p_{T}$ in the $t\bar{t}$ CR. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZu LH signals are also shown separately, normalized to $10^{3}$ times the best fit of the signal yield. The last bin includes the overflow. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZu LH coupling extraction. The distribution is for the third lepton $p_{T}$ in the $t\bar{t}$ CR. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZu LH signals are also shown separately, normalized to $10^{3}$ times the best fit of the signal yield. The last bin includes the overflow. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZu LH coupling extraction. The distribution is for the leading lepton $p_{T}$ in the $t\bar{t}Z$ CR. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZu LH signals are also shown separately, normalized to $10^{3}$ times the best fit of the signal yield. The last bin includes the overflow. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Comparison between data and background prediction after the fit to data (Post-Fit) for the FCNC tZu LH coupling extraction. The distribution is for the $D_{1}$ discriminant in the $t\bar{t}Z$ CR. The uncertainty band includes both the statistical and systematic uncertainties in the background prediction. The FCNC tZu LH signals are also shown separately, normalized to $10^{3}$ times the best fit of the signal yield. The lower panels show the ratios of the data (Data) to the background prediction (Bkg.).

Measured differential cross section as a function of observable $ p_{T}^{ll} - p_{T}^{\gamma}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{ll} + p_{T}^{\gamma}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{ll} + p_{T}^{\gamma}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ \Delta R (l,l)$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ \Delta R (l,l)$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ N_{jets}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ N_{jets}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{Jet1}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{Jet1}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{Jet2}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{Jet2}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{Jet2}/p_{T}^{Jet1}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{Jet2}/p_{T}^{Jet1}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ m_{jj}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ m_{jj}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ m_{ll\gamma j}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ m_{ll\gamma j}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ H_{T}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ H_{T}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{\gamma} / \sqrt{H_{T}}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{\gamma} / \sqrt{H_{T}}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ \Delta \phi (Jet,\gamma)$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ \Delta \phi (Jet,\gamma)$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{ll\gamma j}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{ll\gamma j}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ \phi_{CS}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ \phi_{CS}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ \cos \theta_{CS}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ \cos \theta_{CS}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{ll\gamma} / m_{ll\gamma}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{ll\gamma} / m_{ll\gamma}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{ll\gamma} / m_{ll\gamma}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{ll\gamma} / m_{ll\gamma}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{ll\gamma} / m_{ll\gamma}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{ll\gamma} / m_{ll\gamma}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{ll} - p_{T}^{\gamma}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{ll} - p_{T}^{\gamma}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{ll} - p_{T}^{\gamma}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{ll} - p_{T}^{\gamma}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{ll} - p_{T}^{\gamma}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{ll} - p_{T}^{\gamma}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{ll\gamma j}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{ll\gamma j}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{ll\gamma j}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{ll\gamma j}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{ll\gamma j}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Measured differential cross section as a function of observable $ p_{T}^{ll\gamma j}$. Error on the measured cross-section include all the systematic uncertainties. SM predictions are produced with the event generators at particle level: Sherpa 2.2.4, Sherpa 2.2.11, MadGraph5_aMC@NLO, and MiNNLO$_{PS}$. Fixed order calculations results use MATRIX NNLO. Error represent statistical uncertainty and theoretical uncertainty (PDF and Scale variations).

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{\gamma} / \sqrt{H_{T}}$ (Fig. 8 (b))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{\gamma} / \sqrt{H_{T}}$ (Fig. 8 (b))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ H_{T}$ (Fig. 8 (a))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ H_{T}$ (Fig. 8 (a))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ \Delta \phi (Jet,\gamma)$ (Fig. 8 (c))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ \Delta \phi (Jet,\gamma)$ (Fig. 8 (c))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ \Delta R (l,l)$ (Fig. 5 (d))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ \Delta R (l,l)$ (Fig. 5 (d))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{ll} - p_{T}^{\gamma}$ (Fig. 5 (b))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{ll} - p_{T}^{\gamma}$ (Fig. 5 (b))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{ll} - p_{T}^{\gamma} \textrm{ in bin } p_{T}^{ll} + p_{T}^{\gamma} < 200 GeV$ (Fig. 11 (a))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{ll} - p_{T}^{\gamma} \textrm{ in bin } p_{T}^{ll} + p_{T}^{\gamma} < 200 GeV$ (Fig. 11 (a))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{ll} - p_{T}^{\gamma} \textrm{ in bin } 200 GeV < p_{T}^{ll} + p_{T}^{\gamma} < 300 GeV$ (Fig. 11 (b))