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.

Bayesian upper exclusion limits on the ratio $g_{W'}/g_{W}$ for an SSM-like W$\prime$ boson are shown. The unity coupling ratio (blue dotted curve) corresponds to the SSM common benchmark. The 68 and 95% quantiles of the limits are represented by the green and yellow bands, respectively.

Bayesian lower exclusion limits on the NUGIM G(221) mixing angle $\cot(\theta_{E})$ are shown as a function of the W$\prime$ boson mass. The theoretically excluded region is shaded in grey. The 68 and 95% quantiles of the limits are represented by the green and yellow bands, respectively.

Bayesian upper exclusion limits at 95% CL on the product of the production cross section and branching fraction of a QBH in an associated $ au$ lepton and neutrino final state. Minimum threshold masses $m_{th}$ of up to 6.6 TeV are excluded at 95% CL. The observed limit (solid line) is obtained from the intersection with the LO QBH cross section (blue dotted curve). The 68 and 95% quantiles of the limits are represented by the green and yellow bands, respectively.

Bayesian upper limits at 95% CL on the cross section of the process $pp\rightarrow\tau\nu$ mediated via LQ exchange in the t-channel. The inner (green) band and the outer (yellow) band indicate the regions containing 68 and 95%, respectively, of the distribution of limits expected under the background-only hypothesis. The predicted LQ cross section at LO in the three coupling benchmark scenarios is depicted in different colors for $g_{U}=1$. The uncertainty bandy correspond to the sum in quadrature of PDF and scale variations. The first benchmark scenario considers only couplings to left-handed SM fermions (i.e. $\beta_{\text{R}}^{ij}=0$) and is referred to as "best fit LH". The second benchmark, referred to as "best fit LH+RH", considers $|\beta_{\text{R}}^{\text{b}\tau}|=1$ and all other $\beta_{\text{R}}^{ij}=0$. In the third "democratic" benchmark, equal couplings only to LH fermions are assumed, i.e. $\beta_{\text{L}}^{ij}=1$ and $\beta_{\text{R}}^{ij}=0$ for all $i$ and $j$.

Expected and observed lower limits of the LQ mass as a function of the coupling $g_{U}$ in the LH scenario. The blue band shows the 68% and 95% regions of $g_{U}$ preferred by the fit to the b anomalies data.

Expected and observed lower limits of the LQ mass as a function of the coupling $g_{U}$ in the LH+RH scenario. The blue band shows the 68% and 95% regions of $g_{U}$ preferred by the fit to the b anomalies data.

Expected and observed lower limits of the LQ mass as a function of the coupling $g_{U}$ in the democratic scenario. The blue band shows the 68% and 95% regions of $g_{U}$ preferred by the fit to the b anomalies data.

Bayesian upper exclusion limits at 95% CL on each of the Wilson coefficients described by the EFT model based on 2016-2018 data. The three different coupling types represent a left-handed vector coupling ($\epsilon^{cb}_{L}$), tensor-like coupling ($\epsilon^{cb}_{T}$), and scalar-tensor-like coupling ($\epsilon^{cb}_{S_{L}}$). The 68 and 95% quantiles of the limits are represented by the green and yellow bands, respectively.

Summary of exclusion limits (expected and observed) calculated at 95% CL for full Run-2 CMS data.

Background prediction and observed data yields in the signal region bins. The background yields are obtained from the background-only fit and serve as input to the simplified likelihood reinterpretation scheme. The naming of the bins is "year_binnumber", following the binning from Figure 4.

Matrix of covariance coefficients between signal region bins. The coefficients are obtained from the background-only fit and serve as input to the simplified likelihood reinterpretation scheme. The naming of the bins is "year_binnumber", following the binning used in Figure 4.

Predicted signal yields for the 2017 data-taking period, corresponding to 41.3 fb$^{-1}$, after the application of each search requirement (cumulative) for various signal hypothesis. The requirements listed are presented as total efficiencies w.r.t. the previous selection step.

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.

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

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>

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.

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))

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))

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} > 300 GeV$ (Fig. 11 (c))

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} > 300 GeV$ (Fig. 11 (c))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ m_{jj}$ (Fig. 7 (a))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ m_{jj}$ (Fig. 7 (a))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ m_{ll\gamma j}$ (Fig. 7 (b))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ m_{ll\gamma j}$ (Fig. 7 (b))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ N_{jets}$ (Fig. 6 (a))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ N_{jets}$ (Fig. 6 (a))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{Jet1}$ (Fig. 6 (b))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{Jet1}$ (Fig. 6 (b))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{Jet2}$ (Fig. 6 (c))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{Jet2}$ (Fig. 6 (c))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{Jet2}/p_{T}^{Jet1}$ (Fig. 6 (d))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{Jet2}/p_{T}^{Jet1}$ (Fig. 6 (d))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{ll}$ (Fig. 5 (a))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{ll}$ (Fig. 5 (a))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{ll\gamma j}$ (Fig. 8 (d))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{ll\gamma j}$ (Fig. 8 (d))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{ll\gamma j} \textrm{ in bin } p_{T}^{ll\gamma} < 50 GeV$ (Fig. 12 (a))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{ll\gamma j} \textrm{ in bin } p_{T}^{ll\gamma} < 50 GeV$ (Fig. 12 (a))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{ll\gamma j} \textrm{ in bin } 50 GeV < p_{T}^{ll\gamma} < 75 GeV$ (Fig. 12 (b))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{ll\gamma j} \textrm{ in bin } 50 GeV < p_{T}^{ll\gamma} < 75 GeV$ (Fig. 12 (b))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{ll\gamma j} \textrm{ in bin } p_{T}^{ll\gamma} > 75 GeV$ (Fig. 12 (c))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{ll\gamma j} \textrm{ in bin } p_{T}^{ll\gamma} > 75 GeV$ (Fig. 12 (c))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{ll\gamma} / m_{ll\gamma} \textrm{ in bin } 125 GeV < m_{ll\gamma} < 200 GeV$ (Fig. 10 (a))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{ll\gamma} / m_{ll\gamma} \textrm{ in bin } 125 GeV < m_{ll\gamma} < 200 GeV$ (Fig. 10 (a))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{ll\gamma} / m_{ll\gamma} \textrm{ in bin } 200 GeV < m_{ll\gamma} < 300 GeV$ (Fig. 10 (b))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{ll\gamma} / m_{ll\gamma} \textrm{ in bin } 200 GeV < m_{ll\gamma} < 300 GeV$ (Fig. 10 (b))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{ll\gamma} / m_{ll\gamma} \textrm{ in bin } m_{ll\gamma} > 300 GeV$ (Fig. 10 (c))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{ll\gamma} / m_{ll\gamma} \textrm{ in bin } m_{ll\gamma} > 300 GeV$ (Fig. 10 (c))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{ll} + p_{T}^{\gamma}$ (Fig. 5 (c))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ p_{T}^{ll} + p_{T}^{\gamma}$ (Fig. 5 (c))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ \cos \theta_{CS}$ (Fig. 9 (b))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ \cos \theta_{CS}$ (Fig. 9 (b))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ \phi_{CS}$ (Fig. 9 (a))

Correlation matrix of the statistical uncertainties in the measured fiducial cross section for the observable $ \phi_{CS}$ (Fig. 9 (a))

The observed and expected 95% CL upper limits on the minimal coupling vector LQ pair production cross-sections assuming B = 1 as a function of m$_\text{LQ}$.

The observed and expected 95% CL upper limits on the Yang-Mills coupling vector LQ pair production cross-sections assuming B = 1 as a function of m$_\text{LQ}$.

The observed and expected 95% CL upper limits on the branching ratio into charged leptons as a function of m$_\text{LQ}$ for the scalar LQ.

The observed and expected 95% CL upper limits on the branching ratio into charged leptons as a function of m$_\text{LQ}$ for the vector LQ minimal coupling scenario.

The observed and expected 95% CL upper limits on the branching ratio into charged leptons as a function of m$_\text{LQ}$ for the vector LQ Yang-Mills coupling scenario.

Breakdown of systematic uncertainties in the measured inclusive cross-section. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper.

Data bootstrap post unfolding for the inclusive cross-section. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb].

Differential absolute cross-section for $\textrm{p}_{\textrm{T}}^{l}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The last bin of the distribution contains the overflow.

Covariance Matrix (statistical uncertainties only) for the differential absolute cross-section for $\textrm{p}_{\textrm{T}}^{l}$. The last bin of the distribution contains the overflow.

Data bootstrap post unfolding for the differential absolute cross-section for $\textrm{p}_{\textrm{T}}^{l}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb/GeV]. The last bin of the distribution contains the overflow.

Differential absolute cross-section for $|\eta^{l}|$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper.

Covariance Matrix (statistical uncertainties only) for the differential absolute cross-section for $|\eta^{l}|$.

Data bootstrap post unfolding for the differential absolute cross-section for $|\eta^{l}|$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb].

Differential absolute cross-section for $\textrm{E}^{e} + \textrm{E}^{\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The last bin of the distribution contains the overflow.

Covariance Matrix (statistical uncertainties only) for the differential absolute cross-section for $\textrm{E}^{e} + \textrm{E}^{\mu}$. The last bin of the distribution contains the overflow.

Data bootstrap post unfolding for the differential absolute cross-section for $\textrm{E}^{e} + \textrm{E}^{\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb/GeV]. The last bin of the distribution contains the overflow.

Differential absolute cross-section for $m^{e\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The last bin of the distribution contains the overflow.

Covariance Matrix (statistical uncertainties only) for the differential absolute cross-section for $m^{e\mu}$. The last bin of the distribution contains the overflow.

Data bootstrap post unfolding for the differential absolute cross-section for $m^{e\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb/GeV]. The last bin of the distribution contains the overflow.

Differential absolute cross-section for $\textrm{p}_{\textrm{T}}^{e} + \textrm{p}_{\textrm{T}}^{\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The last bin of the distribution contains the overflow.

Covariance Matrix (statistical uncertainties only) for the differential absolute cross-section for $\textrm{p}_{\textrm{T}}^{e} + \textrm{p}_{\textrm{T}}^{\mu}$. The last bin of the distribution contains the overflow.

Data bootstrap post unfolding for the differential absolute cross-section for $\textrm{p}_{\textrm{T}}^{e} + \textrm{p}_{\textrm{T}}^{\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb/GeV]. The last bin of the distribution contains the overflow.

Differential absolute cross-section for $\textrm{p}_{\textrm{T}}^{e\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The last bin of the distribution contains the overflow.

Covariance Matrix (statistical uncertainties only) for the differential absolute cross-section for $\textrm{p}_{\textrm{T}}^{e\mu}$. The last bin of the distribution contains the overflow.

Data bootstrap post unfolding for the differential absolute cross-section for $\textrm{p}_{\textrm{T}}^{e\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb/GeV]. The last bin of the distribution contains the overflow.

Differential absolute cross-section for $|\Delta\phi^{e\mu}|$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper.

Covariance Matrix (statistical uncertainties only) for the differential absolute cross-section for $|\Delta\phi^{e\mu}|$.

Data bootstrap post unfolding for the differential absolute cross-section for $|\Delta\phi^{e\mu}|$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb/rad].

Differential absolute cross-section for $|y^{e\mu}|$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper.

Covariance Matrix (statistical uncertainties only) for the differential absolute cross-section for $|y^{e\mu}|$.

Data bootstrap post unfolding for the differential absolute cross-section for $|y^{e\mu}|$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb].

Double differential absolute cross-section for $|y^{e\mu}|$ x $m^{e\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The highest invariant mass bin of the two distributions contains the overflows.

Covariance Matrix (statistical uncertainties only) for the double differential absolute cross-section for $|y^{e\mu}|$ x $m^{e\mu}$. The highest invariant mass bin of the two distributions contains the overflows.

Data bootstrap post unfolding for the double differential absolute cross-section for $|y^{e\mu}|$ x $m^{e\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb/GeV]. The highest invariant mass bin of the two distributions contains the overflows.

Double differential absolute cross-section for $|\Delta\phi^{e\mu}|$ x $m^{e\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The highest invariant mass bin of the two distributions contains the overflows.

Covariance Matrix (statistical uncertainties only) for the double differential absolute cross-section for $|\Delta\phi^{e\mu}|$ x $m^{e\mu}$. The highest invariant mass bin of the two distributions contains the overflows.

Data bootstrap post unfolding for the double differential absolute cross-section for $|\Delta\phi^{e\mu}|$ x $m^{e\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb/GeV/rad]. The highest invariant mass bin of the two distributions contains the overflows.

Double differential absolute cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{p}_{\textrm{T}}^{e\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The highest $\textrm{p}_{\textrm{T}}^{e\mu}$ bin contains the overflows.

Covariance Matrix (statistical uncertainties only) for the double differential absolute cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{p}_{\textrm{T}}^{e\mu}$. The highest $\textrm{p}_{\textrm{T}}^{e\mu}$ bin contains the overflows.

Data bootstrap post unfolding for the double differential absolute cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{p}_{\textrm{T}}^{e\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb/GeV/rad]. The highest $\textrm{p}_{\textrm{T}}^{e\mu}$ bin contains the overflows.

Double differential absolute cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{E}^{e} + \textrm{E}^{\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The highest $\textrm{E}^{e} + \textrm{E}^{\mu}$ bin contains the overflows.

Covariance Matrix (statistical uncertainties only) for the double differential absolute cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{E}^{e} + \textrm{E}^{\mu}$. The highest $\textrm{E}^{e} + \textrm{E}^{\mu}$ bin contains the overflows.

Data bootstrap post unfolding for the double differential absolute cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{E}^{e} + \textrm{E}^{\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb/GeV/rad]. The highest $\textrm{E}^{e} + \textrm{E}^{\mu}$ bin contains the overflows.

Differential normalised cross-section for $\textrm{p}_{\textrm{T}}^{l}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The last bin of the distribution contains the overflow.

Covariance Matrix (statistical uncertainties only) for the differential normalised cross-section for $\textrm{p}_{\textrm{T}}^{l}$. The last bin of the distribution contains the overflow.

Data bootstrap post unfolding for the differential normalised cross-section for $\textrm{p}_{\textrm{T}}^{l}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [1/GeV]. The last bin of the distribution contains the overflow.

Differential normalised cross-section for $|\eta^{l}|$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper.

Covariance Matrix (statistical uncertainties only) for the differential normalised cross-section for $|\eta^{l}|$.

Data bootstrap post unfolding for the differential normalised cross-section for $|\eta^{l}|$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry.

Differential normalised cross-section for $\textrm{E}^{e} + \textrm{E}^{\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The last bin of the distribution contains the overflow.

Covariance Matrix (statistical uncertainties only) for the differential normalised cross-section for $\textrm{E}^{e} + \textrm{E}^{\mu}$. The last bin of the distribution contains the overflow.

Data bootstrap post unfolding for the differential normalised cross-section for $\textrm{E}^{e} + \textrm{E}^{\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry. The last bin of the distribution contains the overflow.

Differential normalised cross-section for $m^{e\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The last bin of the distribution contains the overflow.

Covariance Matrix (statistical uncertainties only) for the differential normalised cross-section for $m^{e\mu}$. The last bin of the distribution contains the overflow.

Data bootstrap post unfolding for the differential normalised cross-section for $m^{e\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry. The last bin of the distribution contains the overflow.

Differential normalised cross-section for $\textrm{p}_{\textrm{T}}^{e} + \textrm{p}_{\textrm{T}}^{\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The last bin of the distribution contains the overflow.

Covariance Matrix (statistical uncertainties only) for the differential normalised cross-section for $\textrm{p}_{\textrm{T}}^{e} + \textrm{p}_{\textrm{T}}^{\mu}$. The last bin of the distribution contains the overflow.

Data bootstrap post unfolding for the differential normalised cross-section for $\textrm{p}_{\textrm{T}}^{e} + \textrm{p}_{\textrm{T}}^{\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [1/GeV]. The last bin of the distribution contains the overflow.

Differential normalised cross-section for $\textrm{p}_{\textrm{T}}^{e\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The last bin of the distribution contains the overflow.

Covariance Matrix (statistical uncertainties only) for the differential normalised cross-section for $\textrm{p}_{\textrm{T}}^{e\mu}$. The last bin of the distribution contains the overflow.

Data bootstrap post unfolding for the differential normalised cross-section for $\textrm{p}_{\textrm{T}}^{e\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [1/GeV]. The last bin of the distribution contains the overflow.

Differential normalised cross-section for $|\Delta\phi^{e\mu}|$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper.

Covariance Matrix (statistical uncertainties only) for the differential normalised cross-section for $|\Delta\phi^{e\mu}|$.

Data bootstrap post unfolding for the differential normalised cross-section for $|\Delta\phi^{e\mu}|$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [1/rad].

Differential normalised cross-section for $|y^{e\mu}|$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper.

Covariance Matrix (statistical uncertainties only) for the differential normalised cross-section for $|y^{e\mu}|$.

Data bootstrap post unfolding for the differential normalised cross-section for $|y^{e\mu}|$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry.

Double differential normalised cross-section for $|y^{e\mu}|$ x $m^{e\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The highest invariant mass bin of the two distributions contains the overflows.

Covariance Matrix (statistical uncertainties only) for the double differential normalised cross-section for $|y^{e\mu}|$ x $m^{e\mu}$. The highest invariant mass bin of the two distributions contains the overflows.

Data bootstrap post unfolding for the double differential normalised cross-section for $|y^{e\mu}|$ x $m^{e\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [1/GeV]. The highest invariant mass bin of the two distributions contains the overflows.

Double differential normalised cross-section for $|\Delta\phi^{e\mu}|$ x $m^{e\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The highest invariant mass bin of the two distributions contains the overflows.

Covariance Matrix (statistical uncertainties only) for the double differential normalised cross-section for $|\Delta\phi^{e\mu}|$ x $m^{e\mu}$. The highest invariant mass bin of the two distributions contains the overflows.

Data bootstrap post unfolding for the double differential normalised cross-section for $|\Delta\phi^{e\mu}|$ x $m^{e\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [1/GeV/rad]. The highest invariant mass bin of the two distributions contains the overflows.

Double differential normalised cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{p}_{\textrm{T}}^{e\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The highest $\textrm{p}_{\textrm{T}}^{e\mu}$ bin contains the overflows.

Covariance Matrix (statistical uncertainties only) for the double differential normalised cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{p}_{\textrm{T}}^{e\mu}$. The highest $\textrm{p}_{\textrm{T}}^{e\mu}$ bin contains the overflows.

Data bootstrap post unfolding for the double differential normalised cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{p}_{\textrm{T}}^{e\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [1/GeV/rad]. The highest $\textrm{p}_{\textrm{T}}^{e\mu}$ bin contains the overflows.

Double differential normalised cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{E}^{e} + \textrm{E}^{\mu}$. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper. The highest $\textrm{E}^{e} + \textrm{E}^{\mu}$ bin contains the overflows.

Covariance Matrix (statistical uncertainties only) for the double differential normalised cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{E}^{e} + \textrm{E}^{\mu}$. The highest $\textrm{E}^{e} + \textrm{E}^{\mu}$ bin contains the overflows.

Data bootstrap post unfolding for the double differential normalised cross-section for $|\Delta\phi^{e\mu}|$ x $\textrm{E}^{e} + \textrm{E}^{\mu}$. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [1/GeV/rad]. The highest $\textrm{E}^{e} + \textrm{E}^{\mu}$ bin contains the overflows.

The rejection factor (inverse of the efficiency) for jets that have the other origins (background jets) are shown. The background jet rejection factor is calculated using pre-selected large-$R$ jets in the sample of simulated $Z\rightarrow\nu\nu$ + jets events, dominated by initial state radiation (ISR) jets. The efficiency correction factors are applied on the background rejection for the $W_{qq}$/$Z_{qq}$-tagging. The systematic uncertainty is represented by the hashed bands.

The boson-tagging efficiency for jets arising from $Z/h$ bosons decaying into $b\bar{b}$ (signal jets). The signal jet efficiency of $Z_{bb}$/$h_{bb}$-tagging is evaluated using a sample of pre-selected large-$R$ jets ($p_{\textrm{T}}>200~\textrm{GeV}, |\eta|<2.0, m_{J} > 40~\textrm{GeV}$) in the simulated $(\tilde{W},\tilde{B})$ simplified model signal events with $\Delta m(\tilde{\chi}_{\textrm{heavy}},~\tilde{\chi}_{\textrm{light}}) \ge 400~\textrm{GeV}$. The jets are matched with generator-level $Z/h$-bosons by $\Delta R<1.0$ which decay into $b\bar{b}$. The systematic uncertainty is represented by the hashed bands.

The rejection factor (inverse of the efficiency) for jets that have the other origins (background jets) are shown. The background jet rejection factor is calculated using pre-selected large-$R$ jets in the sample of simulated $Z\rightarrow\nu\nu$ + jets events, dominated by initial state radiation (ISR) jets. As for the $Z_{bb}$/$h_{bb}$-tagging, the rejection is shown as the function of number of $b$- or $c$-quarks contained in the large-$R$ jet within $\Delta R<1.0$. The systematic uncertainty is represented by the hashed bands.

The signal jet efficiency and background jet rejection for the $W_{qq}$-/$Z_{qq}$-tagging working point used in the analysis in comparison with those for the nominal working points (''50% efficiency W/Z tagger'' in ATL-PHYS-PUB-2020-017). The definition of the efficiency and rejection follow the caption of Figure 4. Note that the efficiency for the nominal working point is not 50% since they are tuned using jet samples with different pre-selection.

The signal jet efficiency and background jet rejection for the $W_{qq}$-/$Z_{qq}$-tagging working point used in the analysis in comparison with those for the nominal working points (''50% efficiency W/Z tagger'' in ATL-PHYS-PUB-2020-017). The definition of the efficiency and rejection follow the caption of Figure 4. Note that the efficiency for the nominal working point is not 50% since they are tuned using jet samples with different pre-selection.

The signal jet efficiency and background jet rejection for the $W_{qq}$-/$Z_{qq}$-tagging working point used in the analysis in comparison with those for the nominal working points (''50% efficiency W/Z tagger'' in ATL-PHYS-PUB-2020-017). The definition of the efficiency and rejection follow the caption of Figure 4. Note that the efficiency for the nominal working point is not 50% since they are tuned using jet samples with different pre-selection.

The signal jet efficiency and background jet rejection for the $W_{qq}$-/$Z_{qq}$-tagging working point used in the analysis in comparison with those for the nominal working points (''50% efficiency W/Z tagger'' in ATL-PHYS-PUB-2020-017). The definition of the efficiency and rejection follow the caption of Figure 4. Note that the efficiency for the nominal working point is not 50% since they are tuned using jet samples with different pre-selection.

The total post-fit uncertainty in each of the SRs and VRs.

Summary of the observed data and predicted SM background in all SRs. The background prediction in SR-4Q (SR-2B2Q) is obtained by a background-only fit to CR0L-4Q (CR0L-2B2Q). The total systematic uncertainty on the background prediction is shown by the hatched area. Distributions of a few representative signals are overlaid. A few representative signals are overlaid. For the $(\tilde{W},\tilde{B})$ simplified model models, the label $(x,y)$ indicates $(m(\tilde{\chi}_{1}^{\pm}), m(\tilde{\chi}_{1}^{0}))=(x,y)~\textrm{GeV}$. The bottom panel shows the statistical significance of the discrepancy between the observed number of events and the SM expectation.

Number of observed data events and the SM backgrounds in the SRs and the CR0L bins. The SM backgrounds are predicted by the background-only fits. Note that the relative uncertainty on the expected yield in the CRs between the reducible backgrounds are identical, given that a common normalization factor is assigned for all of them in the fit. Yields for negligible backgrounds are indicated by "0.0000001±0.0000001" in the table.

Number of observed data events and the post-fit SM background prediction in the VR1L(1Y) bins and the corresponding CR1L(1Y) bins. "-" indicates negligibly small contribution. Note that the relative uncertainty on the expected yield in the CRs between the reducible backgrounds are identical, given that a common normalization factor is assigned for all of them in the fit. Yields for negligible backgrounds are indicated by "0.0000001±0.0000001" in the table.

$m_{\textrm{eff}}$ distribution in SR-4Q-VV. The post-fit SM background expectation using the background-only fit is shown in a histogram stack. Distributions of a few representative signals are overlaid. The bottom panels show the ratio of the observed data to the background prediction. The selection criteria on the variable shown by each plot is removed, while the arrow indicates the cut value to define the region. A few representative signals are overlaid. For the $(\tilde{W},\tilde{B})$ simplified model models, the label $(x,y)$ indicates $(m(\tilde{\chi}_{1}^{\pm}), m(\tilde{\chi}_{1}^{0}))=(x,y)~\text{GeV}$.

Distribution of $m(J_{1})$ in SR-4Q-VV. The post-fit SM background expectation using the background-only fit is shown in a histogram stack. Distributions of relevant signals are overlaid. The bottom panels show the ratio of the observed data to the background prediction. The selection criteria on the variable shown by each plot is removed, while the arrow indicates the cut value to define the region. A few representative signals are overlaid. For the $(\tilde{W},\tilde{B})$ simplified model models, the label $(x,y)$ indicates $(m(\tilde{\chi}_{1}^{\pm}), m(\tilde{\chi}_{1}^{0}))=(x,y)~\text{GeV}$.

Distribution of $D_2(J_{1})$ in SR-4Q-VV. For $D_2$, the cut value applied for $V_{qq}$-tagging ($D_{2,\text{cut}}$) is subtracted as the off-set so that $D_2-D_{2,\text{cut}}(J)<0$ represents $J$ passing the $D_2$ selection. The post-fit SM background expectation using the background-only fit is shown in a histogram stack. Distributions of relevant signals are overlaid. The bottom panels show the ratio of the observed data to the background prediction. The selection criteria on the variable shown by each plot is removed, while the arrow indicates the cut value to define the region. A few representative signals are overlaid. For the $(\tilde{W},\tilde{B})$ simplified model models, the label $(x,y)$ indicates $(m(\tilde{\chi}_{1}^{\pm}), m(\tilde{\chi}_{1}^{0}))=(x,y)~\text{GeV}$.

Distribution of $m(J_{2})$ in SR-4Q-VV. The post-fit SM background expectation using the background-only fit is shown in a histogram stack. Distributions of relevant signals are overlaid. The bottom panels show the ratio of the observed data to the background prediction. The selection criteria on the variable shown by each plot is removed, while the arrow indicates the cut value to define the region. A few representative signals are overlaid. For the $(\tilde{W},\tilde{B})$ simplified model models, the label $(x,y)$ indicates $(m(\tilde{\chi}_{1}^{\pm}), m(\tilde{\chi}_{1}^{0}))=(x,y)~\text{GeV}$.

Distribution of $D_2(J_{2})$ in SR-4Q-VV. For $D_2$, the cut value applied for $V_{qq}$-tagging ($D_{2,\text{cut}}$) is subtracted as the off-set so that $D_2-D_{2,\text{cut}}(J)<0$ represents $J$ passing the $D_2$ selection. The post-fit SM background expectation using the background-only fit is shown in a histogram stack. Distributions of relevant signals are overlaid. The bottom panels show the ratio of the observed data to the background prediction. The selection criteria on the variable shown by each plot is removed, while the arrow indicates the cut value to define the region. A few representative signals are overlaid. For the $(\tilde{W},\tilde{B})$ simplified model models, the label $(x,y)$ indicates $(m(\tilde{\chi}_{1}^{\pm}), m(\tilde{\chi}_{1}^{0}))=(x,y)~\text{GeV}$.

$m_{\textrm{T2}}$ distributions in SR-2B2Q-VZ. The post-fit SM background expectation using the background-only fit is shown in a histogram stack. Distributions of a few representative signals are overlaid. The bottom panels show the ratio of the observed data to the background prediction. The selection criteria on the variable shown by each plot is removed, while the arrow indicates the cut value to define the region. A few representative signals are overlaid. For the $(\tilde{W},\tilde{B})$ simplified model models, the label $(x,y)$ indicates $(m(\tilde{\chi}_{1}^{\pm}), m(\tilde{\chi}_{1}^{0}))=(x,y)~\text{GeV}$.

Distribution of $m(J_{bb})$ in SR-2B2Q-VZ. The post-fit SM background expectation using the background-only fit is shown in a histogram stack. Distributions of relevant signals are overlaid. The bottom panels show the ratio of the observed data to the background prediction. The selection criteria on the variable shown by each plot is removed, while the arrow indicates the cut value to define the region. A few representative signals are overlaid. For the $(\tilde{W},\tilde{B})$ simplified model models, the label $(x,y)$ indicates $(m(\tilde{\chi}_{1}^{\pm}), m(\tilde{\chi}_{1}^{0}))=(x,y)~\text{GeV}$.

Distribution of $m_{\textrm{eff}} $ in SR-2B2Q-VZ. The post-fit SM background expectation using the background-only fit is shown in a histogram stack. Distributions of relevant signals are overlaid. The bottom panels show the ratio of the observed data to the background prediction. The selection criteria on the variable shown by each plot is removed, while the arrow indicates the cut value to define the region. A few representative signals are overlaid. For the $(\tilde{W},\tilde{B})$ simplified model models, the label $(x,y)$ indicates $(m(\tilde{\chi}_{1}^{\pm}), m(\tilde{\chi}_{1}^{0}))=(x,y)~\text{GeV}$.

$m_{\textrm{T2}}$ distributions in SR-2B2Q-Vh. The post-fit SM background expectation using the background-only fit is shown in a histogram stack. Distributions of a few representative signals are overlaid. The bottom panels show the ratio of the observed data to the background prediction. The selection criteria on the variable shown by each plot is removed, while the arrow indicates the cut value to define the region. A few representative signals are overlaid. For the $(\tilde{W},\tilde{B})$ simplified model models, the label $(x,y)$ indicates $(m(\tilde{\chi}_{1}^{\pm}), m(\tilde{\chi}_{1}^{0}))=(x,y)~\text{GeV}$.

Distribution of $m(J_{bb})$ in SR-2B2Q-Vh. The post-fit SM background expectation using the background-only fit is shown in a histogram stack. Distributions of relevant signals are overlaid. The bottom panels show the ratio of the observed data to the background prediction. The selection criteria on the variable shown by each plot is removed, while the arrow indicates the cut value to define the region. A few representative signals are overlaid. For the $(\tilde{W},\tilde{B})$ simplified model models, the label $(x,y)$ indicates $(m(\tilde{\chi}_{1}^{\pm}), m(\tilde{\chi}_{1}^{0}))=(x,y)~\text{GeV}$.

Distribution of $m_{\textrm{eff}} $ in SR-2B2Q-Vh. The post-fit SM background expectation using the background-only fit is shown in a histogram stack. Distributions of relevant signals are overlaid. The bottom panels show the ratio of the observed data to the background prediction. The selection criteria on the variable shown by each plot is removed, while the arrow indicates the cut value to define the region. A few representative signals are overlaid. For the $(\tilde{W},\tilde{B})$ simplified model models, the label $(x,y)$ indicates $(m(\tilde{\chi}_{1}^{\pm}), m(\tilde{\chi}_{1}^{0}))=(x,y)~\text{GeV}$.

Exclusion limits for $(\tilde{W},\tilde{B})$ simplifiec model as a function of the produced wino mass ($m(\tilde{\chi}_{1}^{\pm}/\tilde{\chi}_{2}^{0})$) and the bino LSP mass ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed) and observed (solid) 95% CL exclusion limits on $(\tilde{W},\tilde{B})$ simplified models; C1C1-WW. The expected limit for the 1$\sigma$ down variation is not shown as no mass points could be excluded.

Exclusion limits for $(\tilde{W},\tilde{B})$ simplifiec model as a function of the produced wino mass ($m(\tilde{\chi}_{1}^{\pm}/\tilde{\chi}_{2}^{0})$) and the bino LSP mass ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed) and observed (solid) 95% CL exclusion limits on $(\tilde{W},\tilde{B})$ simplified models; C1C1-WW. The expected limit for the 1$\sigma$ down variation is not shown as no mass points could be excluded.

Exclusion limits for $(\tilde{W},\tilde{B})$ simplifiec model as a function of the produced wino mass ($m(\tilde{\chi}_{1}^{\pm}/\tilde{\chi}_{2}^{0})$) and the bino LSP mass ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed) and observed (solid) 95% CL exclusion limits on $(\tilde{W},\tilde{B})$ simplified models; C1C1-WW. The expected limit for the 1$\sigma$ down variation is not shown as no mass points could be excluded.

Exclusion limits for $(\tilde{W},\tilde{B})$ simplifiec model as a function of the produced wino mass ($m(\tilde{\chi}_{1}^{\pm}/\tilde{\chi}_{2}^{0})$) and the bino LSP mass ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed) and observed (solid) 95% CL exclusion limits on $(\tilde{W},\tilde{B})$ simplified models; C1C1-WW. The expected limit for the 1$\sigma$ down variation is not shown as no mass points could be excluded.

Exclusion limits for $(\tilde{W},\tilde{B})$ simplifiec model as a function of the produced wino mass ($m(\tilde{\chi}_{1}^{\pm}/\tilde{\chi}_{2}^{0})$) and the bino LSP mass ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed) and observed (solid) 95% CL exclusion limits on $(\tilde{W},\tilde{B})$ simplified models; C1C1-WW. The expected limit for the 1$\sigma$ down variation is not shown as no mass points could be excluded.

Exclusion limits for $(\tilde{W},\tilde{B})$ simplifiec model as a function of the produced wino mass ($m(\tilde{\chi}_{1}^{\pm}/\tilde{\chi}_{2}^{0})$) and the bino LSP mass ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed) and observed (solid) 95% CL exclusion limits on $(\tilde{W},\tilde{B})$ simplified models; C1N2-WZ.

Exclusion limits for $(\tilde{W},\tilde{B})$ simplifiec model as a function of the produced wino mass ($m(\tilde{\chi}_{1}^{\pm}/\tilde{\chi}_{2}^{0})$) and the bino LSP mass ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed) and observed (solid) 95% CL exclusion limits on $(\tilde{W},\tilde{B})$ simplified models; C1N2-WZ.

Exclusion limits for $(\tilde{W},\tilde{B})$ simplifiec model as a function of the produced wino mass ($m(\tilde{\chi}_{1}^{\pm}/\tilde{\chi}_{2}^{0})$) and the bino LSP mass ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed) and observed (solid) 95% CL exclusion limits on $(\tilde{W},\tilde{B})$ simplified models; C1N2-WZ.

Exclusion limits for $(\tilde{W},\tilde{B})$ simplifiec model as a function of the produced wino mass ($m(\tilde{\chi}_{1}^{\pm}/\tilde{\chi}_{2}^{0})$) and the bino LSP mass ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed) and observed (solid) 95% CL exclusion limits on $(\tilde{W},\tilde{B})$ simplified models; C1N2-WZ.

Exclusion limits for $(\tilde{W},\tilde{B})$ simplifiec model as a function of the produced wino mass ($m(\tilde{\chi}_{1}^{\pm}/\tilde{\chi}_{2}^{0})$) and the bino LSP mass ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed) and observed (solid) 95% CL exclusion limits on $(\tilde{W},\tilde{B})$ simplified models; C1N2-WZ.

Exclusion limits for $(\tilde{W},\tilde{B})$ simplifiec model as a function of the produced wino mass ($m(\tilde{\chi}_{1}^{\pm}/\tilde{\chi}_{2}^{0})$) and the bino LSP mass ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed) and observed (solid) 95% CL exclusion limits on $(\tilde{W},\tilde{B})$ simplified models; C1N2-WZ.

Exclusion limits for $(\tilde{W},\tilde{B})$ simplifiec model as a function of the produced wino mass ($m(\tilde{\chi}_{1}^{\pm}/\tilde{\chi}_{2}^{0})$) and the bino LSP mass ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed) and observed (solid) 95% CL exclusion limits on $(\tilde{W},\tilde{B})$ simplified models; C1N2-Wh.

Exclusion limits for $(\tilde{W},\tilde{B})$ simplifiec model as a function of the produced wino mass ($m(\tilde{\chi}_{1}^{\pm}/\tilde{\chi}_{2}^{0})$) and the bino LSP mass ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed) and observed (solid) 95% CL exclusion limits on $(\tilde{W},\tilde{B})$ simplified models; C1N2-Wh.

Exclusion limits for $(\tilde{W},\tilde{B})$ simplifiec model as a function of the produced wino mass ($m(\tilde{\chi}_{1}^{\pm}/\tilde{\chi}_{2}^{0})$) and the bino LSP mass ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed) and observed (solid) 95% CL exclusion limits on $(\tilde{W},\tilde{B})$ simplified models; C1N2-Wh.

Exclusion limits for $(\tilde{W},\tilde{B})$ simplifiec model as a function of the produced wino mass ($m(\tilde{\chi}_{1}^{\pm}/\tilde{\chi}_{2}^{0})$) and the bino LSP mass ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed) and observed (solid) 95% CL exclusion limits on $(\tilde{W},\tilde{B})$ simplified models; C1N2-Wh.

Exclusion limits for $(\tilde{W},\tilde{B})$ simplifiec model as a function of the produced wino mass ($m(\tilde{\chi}_{1}^{\pm}/\tilde{\chi}_{2}^{0})$) and the bino LSP mass ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed) and observed (solid) 95% CL exclusion limits on $(\tilde{W},\tilde{B})$ simplified models; C1N2-Wh.

Exclusion limits for $(\tilde{W},\tilde{B})$ simplifiec model as a function of the produced wino mass ($m(\tilde{\chi}_{1}^{\pm}/\tilde{\chi}_{2}^{0})$) and the bino LSP mass ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed) and observed (solid) 95% CL exclusion limits on $(\tilde{W},\tilde{B})$ simplified models; C1N2-Wh.

Exclusion limits for the $(\tilde{W},\tilde{B})$ and $(\tilde{H},\tilde{B})$ models shown as a function of the mass of wino/higgsino chargino ($m(\tilde{\chi}_{1}^{\pm})$) and the mass of bino LSP ($m(\tilde{\chi}_{1}^{0})$). Expected limits for various $\mathcal{B}(\tilde{\chi}_{2}^{0} \rightarrow Z \tilde{\chi}_{1}^{0})$ hypotheses overlaid. The outer and inner bundles correspond to the limits for the $(\tilde{W},\tilde{B})$ and $(\tilde{H},\tilde{B})$ models respectively.

Exclusion limits for the $(\tilde{W},\tilde{B})$ and $(\tilde{H},\tilde{B})$ models shown as a function of the mass of wino/higgsino chargino ($m(\tilde{\chi}_{1}^{\pm})$) and the mass of bino LSP ($m(\tilde{\chi}_{1}^{0})$). Observed limits for various $\mathcal{B}(\tilde{\chi}_{2}^{0} \rightarrow Z \tilde{\chi}_{1}^{0})$ hypotheses overlaid. The outer and inner bundles correspond to the limits for the $(\tilde{W},\tilde{B})$ and $(\tilde{H},\tilde{B})$ models respectively.

Exclusion limits for the $(\tilde{W},\tilde{B})$ and $(\tilde{H},\tilde{B})$ models shown as a function of the mass of wino/higgsino chargino ($m(\tilde{\chi}_{1}^{\pm})$) and the mass of bino LSP ($m(\tilde{\chi}_{1}^{0})$). Expected limits for various $\mathcal{B}(\tilde{\chi}_{2}^{0} \rightarrow Z \tilde{\chi}_{1}^{0})$ hypotheses overlaid. The outer and inner bundles correspond to the limits for the $(\tilde{W},\tilde{B})$ and $(\tilde{H},\tilde{B})$ models respectively.

Exclusion limits for the $(\tilde{W},\tilde{B})$ and $(\tilde{H},\tilde{B})$ models shown as a function of the mass of wino/higgsino chargino ($m(\tilde{\chi}_{1}^{\pm})$) and the mass of bino LSP ($m(\tilde{\chi}_{1}^{0})$). Observed limits for various $\mathcal{B}(\tilde{\chi}_{2}^{0} \rightarrow Z \tilde{\chi}_{1}^{0})$ hypotheses overlaid. The outer and inner bundles correspond to the limits for the $(\tilde{W},\tilde{B})$ and $(\tilde{H},\tilde{B})$ models respectively.

Exclusion limits for the $(\tilde{W},\tilde{B})$ and $(\tilde{H},\tilde{B})$ models shown as a function of the mass of wino/higgsino chargino ($m(\tilde{\chi}_{1}^{\pm})$) and the mass of bino LSP ($m(\tilde{\chi}_{1}^{0})$). Expected limits for various $\mathcal{B}(\tilde{\chi}_{2}^{0} \rightarrow Z \tilde{\chi}_{1}^{0})$ hypotheses overlaid. The outer and inner bundles correspond to the limits for the $(\tilde{W},\tilde{B})$ and $(\tilde{H},\tilde{B})$ models respectively.

Exclusion limits for the $(\tilde{W},\tilde{B})$ models shown as a function of the mass of wino/higgsino chargino ($m(\tilde{\chi}_{1}^{\pm})$) and the mass of bino LSP ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed) and observed (solid red) 95% CL exclusion limits for a representative branching ratio $\mathcal{B}(\tilde{\chi}_{2}^{0} \rightarrow Z\tilde{\chi}_{1}^{0})=50$%.

Exclusion limits for the $(\tilde{W},\tilde{B})$ models shown as a function of the mass of wino/higgsino chargino ($m(\tilde{\chi}_{1}^{\pm})$) and the mass of bino LSP ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed) and observed (solid red) 95% CL exclusion limits for a representative branching ratio $\mathcal{B}(\tilde{\chi}_{2}^{0} \rightarrow Z\tilde{\chi}_{1}^{0})=50$%.

Exclusion limits for the $(\tilde{W},\tilde{B})$ and $(\tilde{H},\tilde{B})$ models shown as a function of the mass of wino/higgsino chargino ($m(\tilde{\chi}_{1}^{\pm})$) and the mass of bino LSP ($m(\tilde{\chi}_{1}^{0})$). Observed limits for various $\mathcal{B}(\tilde{\chi}_{2}^{0} \rightarrow Z \tilde{\chi}_{1}^{0})$ hypotheses overlaid. The outer and inner bundles correspond to the limits for the $(\tilde{W},\tilde{B})$ and $(\tilde{H},\tilde{B})$ models respectively.

Exclusion limits for the $(\tilde{W},\tilde{B})$ models shown as a function of the mass of wino/higgsino chargino ($m(\tilde{\chi}_{1}^{\pm})$) and the mass of bino LSP ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed) and observed (solid red) 95% CL exclusion limits for a representative branching ratio $\mathcal{B}(\tilde{\chi}_{2}^{0} \rightarrow Z\tilde{\chi}_{1}^{0})=50$%.

Exclusion limits for the $(\tilde{W},\tilde{B})$ models shown as a function of the mass of wino/higgsino chargino ($m(\tilde{\chi}_{1}^{\pm})$) and the mass of bino LSP ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed) and observed (solid red) 95% CL exclusion limits for a representative branching ratio $\mathcal{B}(\tilde{\chi}_{2}^{0} \rightarrow Z\tilde{\chi}_{1}^{0})=50$%.

Exclusion limits for the $(\tilde{W},\tilde{B})$ and $(\tilde{H},\tilde{B})$ models shown as a function of the mass of wino/higgsino chargino ($m(\tilde{\chi}_{1}^{\pm})$) and the mass of bino LSP ($m(\tilde{\chi}_{1}^{0})$). Expected limits for various $\mathcal{B}(\tilde{\chi}_{2}^{0} \rightarrow Z \tilde{\chi}_{1}^{0})$ hypotheses overlaid. The outer and inner bundles correspond to the limits for the $(\tilde{W},\tilde{B})$ and $(\tilde{H},\tilde{B})$ models respectively.

Exclusion limits for the $(\tilde{W},\tilde{B})$ and $(\tilde{H},\tilde{B})$ models shown as a function of the mass of wino/higgsino chargino ($m(\tilde{\chi}_{1}^{\pm})$) and the mass of bino LSP ($m(\tilde{\chi}_{1}^{0})$). Observed limits for various $\mathcal{B}(\tilde{\chi}_{2}^{0} \rightarrow Z \tilde{\chi}_{1}^{0})$ hypotheses overlaid. The outer and inner bundles correspond to the limits for the $(\tilde{W},\tilde{B})$ and $(\tilde{H},\tilde{B})$ models respectively.

Exclusion limits for the $(\tilde{W},\tilde{B})$ and $(\tilde{H},\tilde{B})$ models shown as a function of the mass of wino/higgsino chargino ($m(\tilde{\chi}_{1}^{\pm})$) and the mass of bino LSP ($m(\tilde{\chi}_{1}^{0})$). Expected limits for various $\mathcal{B}(\tilde{\chi}_{2}^{0} \rightarrow Z \tilde{\chi}_{1}^{0})$ hypotheses overlaid. The outer and inner bundles correspond to the limits for the $(\tilde{W},\tilde{B})$ and $(\tilde{H},\tilde{B})$ models respectively.

Exclusion limits for the $(\tilde{W},\tilde{B})$ and $(\tilde{H},\tilde{B})$ models shown as a function of the mass of wino/higgsino chargino ($m(\tilde{\chi}_{1}^{\pm})$) and the mass of bino LSP ($m(\tilde{\chi}_{1}^{0})$). Observed limits for various $\mathcal{B}(\tilde{\chi}_{2}^{0} \rightarrow Z \tilde{\chi}_{1}^{0})$ hypotheses overlaid. The outer and inner bundles correspond to the limits for the $(\tilde{W},\tilde{B})$ and $(\tilde{H},\tilde{B})$ models respectively.

Exclusion limits for the $(\tilde{W},\tilde{B})$ and $(\tilde{H},\tilde{B})$ models shown as a function of the mass of wino/higgsino chargino ($m(\tilde{\chi}_{1}^{\pm})$) and the mass of bino LSP ($m(\tilde{\chi}_{1}^{0})$). Expected limits for various $\mathcal{B}(\tilde{\chi}_{2}^{0} \rightarrow Z \tilde{\chi}_{1}^{0})$ hypotheses overlaid. The outer and inner bundles correspond to the limits for the $(\tilde{W},\tilde{B})$ and $(\tilde{H},\tilde{B})$ models respectively. The expected limit for the 1$\sigma$ down variation is not shown as no mass points could be excluded.

Exclusion limits for the $(\tilde{H},\tilde{B})$ models shown as a function of the mass of wino/higgsino chargino ($m(\tilde{\chi}_{1}^{\pm})$) and the mass of bino LSP ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed) and observed (solid red) 95% CL exclusion limits for a representative branching ratio $\mathcal{B}(\tilde{\chi}_{2}^{0} \rightarrow Z\tilde{\chi}_{1}^{0})=50$%. The expected limit for the 1$\sigma$ down variation is not shown as no mass points could be excluded.

Exclusion limits for the $(\tilde{W},\tilde{B})$ and $(\tilde{H},\tilde{B})$ models shown as a function of the mass of wino/higgsino chargino ($m(\tilde{\chi}_{1}^{\pm})$) and the mass of bino LSP ($m(\tilde{\chi}_{1}^{0})$). Observed limits for various $\mathcal{B}(\tilde{\chi}_{2}^{0} \rightarrow Z \tilde{\chi}_{1}^{0})$ hypotheses overlaid. The outer and inner bundles correspond to the limits for the $(\tilde{W},\tilde{B})$ and $(\tilde{H},\tilde{B})$ models respectively. The expected limit for the 1$\sigma$ down variation is not shown as no mass points could be excluded.

Exclusion limits for the $(\tilde{H},\tilde{B})$ models shown as a function of the mass of wino/higgsino chargino ($m(\tilde{\chi}_{1}^{\pm})$) and the mass of bino LSP ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed) and observed (solid red) 95% CL exclusion limits for a representative branching ratio $\mathcal{B}(\tilde{\chi}_{2}^{0} \rightarrow Z\tilde{\chi}_{1}^{0})=50$%. The expected limit for the 1$\sigma$ down variation is not shown as no mass points could be excluded.

Exclusion limits for the $(\tilde{H},\tilde{B})$ models shown as a function of the mass of wino/higgsino chargino ($m(\tilde{\chi}_{1}^{\pm})$) and the mass of bino LSP ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed) and observed (solid red) 95% CL exclusion limits for a representative branching ratio $\mathcal{B}(\tilde{\chi}_{2}^{0} \rightarrow Z\tilde{\chi}_{1}^{0})=50$%. The expected limit for the 1$\sigma$ down variation is not shown as no mass points could be excluded.

95% CL exclusion limits for the $(\tilde{W},\tilde{H})$ models. The limits are projected on to a two-dimensional plane of the physical electroweakino masses $(m(\tilde{\chi}_{2}^{\pm}),m(\tilde{\chi}_{1}^{0}))$ representing $(m(\tilde{\chi}_{\textrm{heavy}}),m(\tilde{\chi}_{\textrm{light}}))$.

95% CL exclusion limits for the $(\tilde{W},\tilde{H})$ models. The limits are projected on to a two-dimensional plane of the physical electroweakino masses $(m(\tilde{\chi}_{2}^{\pm}),m(\tilde{\chi}_{1}^{0}))$ representing $(m(\tilde{\chi}_{\textrm{heavy}}),m(\tilde{\chi}_{\textrm{light}}))$.

95% CL exclusion limits for the $(\tilde{W},\tilde{H})$ models. The limits are projected on to a two-dimensional plane of the physical electroweakino masses $(m(\tilde{\chi}_{2}^{\pm}),m(\tilde{\chi}_{1}^{0}))$ representing $(m(\tilde{\chi}_{\textrm{heavy}}),m(\tilde{\chi}_{\textrm{light}}))$.

95% CL exclusion limits for the $(\tilde{W},\tilde{H})$ models. The limits are projected on to a two-dimensional plane of the physical electroweakino masses $(m(\tilde{\chi}_{2}^{\pm}),m(\tilde{\chi}_{1}^{0}))$ representing $(m(\tilde{\chi}_{\textrm{heavy}}),m(\tilde{\chi}_{\textrm{light}}))$.

95% CL exclusion limits for the $(\tilde{W},\tilde{H})$ models. The limits are projected on to a two-dimensional plane of the physical electroweakino masses $(m(\tilde{\chi}_{2}^{\pm}),m(\tilde{\chi}_{1}^{0}))$ representing $(m(\tilde{\chi}_{\textrm{heavy}}),m(\tilde{\chi}_{\textrm{light}}))$.

95% CL exclusion limits for the $(\tilde{W},\tilde{H})$ models. The limits are projected on to a two-dimensional plane of the physical electroweakino masses $(m(\tilde{\chi}_{2}^{\pm}),m(\tilde{\chi}_{1}^{0}))$ representing $(m(\tilde{\chi}_{\textrm{heavy}}),m(\tilde{\chi}_{\textrm{light}}))$.

95% CL exclusion limits for the $(\tilde{H},\tilde{W})$ models. The limits are projected on to a two-dimensional plane of the physical electroweakino masses $(m(\tilde{\chi}_{2}^{\pm}),m(\tilde{\chi}_{1}^{0}))$ representing $(m(\tilde{\chi}_{\textrm{heavy}}),m(\tilde{\chi}_{\textrm{light}}))$. The expected limit for the 1$\sigma$ down variation is not shown as no mass points could be excluded.

95% CL exclusion limits for the $(\tilde{H},\tilde{W})$ models. The limits are projected on to a two-dimensional plane of the physical electroweakino masses $(m(\tilde{\chi}_{2}^{\pm}),m(\tilde{\chi}_{1}^{0}))$ representing $(m(\tilde{\chi}_{\textrm{heavy}}),m(\tilde{\chi}_{\textrm{light}}))$. The expected limit for the 1$\sigma$ down variation is not shown as no mass points could be excluded.

95% CL exclusion limits for the $(\tilde{H},\tilde{W})$ models. The limits are projected on to a two-dimensional plane of the physical electroweakino masses $(m(\tilde{\chi}_{2}^{\pm}),m(\tilde{\chi}_{1}^{0}))$ representing $(m(\tilde{\chi}_{\textrm{heavy}}),m(\tilde{\chi}_{\textrm{light}}))$. The expected limit for the 1$\sigma$ down variation is not shown as no mass points could be excluded.

95% CL exclusion limits for the $(\tilde{H},\tilde{W})$ models. The limits are projected on to a two-dimensional plane of the physical electroweakino masses $(m(\tilde{\chi}_{2}^{\pm}),m(\tilde{\chi}_{1}^{0}))$ representing $(m(\tilde{\chi}_{\textrm{heavy}}),m(\tilde{\chi}_{\textrm{light}}))$. The expected limit for the 1$\sigma$ down variation is not shown as no mass points could be excluded.

95% CL exclusion limits for the $(\tilde{H},\tilde{W})$ models. The limits are projected on to a two-dimensional plane of the physical electroweakino masses $(m(\tilde{\chi}_{2}^{\pm}),m(\tilde{\chi}_{1}^{0}))$ representing $(m(\tilde{\chi}_{\textrm{heavy}}),m(\tilde{\chi}_{\textrm{light}}))$. The expected limit for the 1$\sigma$ down variation is not shown as no mass points could be excluded.

95% CL exclusion limits for the $(\tilde{W},\tilde{H})$ models. The limits are projected on to a two-dimensional plane as a function of the wino/higgsino mass parameters: $(M_2,\mu)$. For the limits shown on the $(M_2,\mu)$ plane, the excluded regions are indicated by the area inside the contours. The round excluded area in the top part corresponds to the excluded parameter space in the $(\tilde{W},~\tilde{H})$ model $(M_{2} > |\mu|)$, while the two small areas in at the bottom are that in the $(\tilde{H},~\tilde{W})$ model $(M_{2} < |\mu|)$.

95% CL exclusion limits for the $(\tilde{W},\tilde{H})$ models. The limits are projected on to a two-dimensional plane as a function of the wino/higgsino mass parameters: $(M_2,\mu)$. For the limits shown on the $(M_2,\mu)$ plane, the excluded regions are indicated by the area inside the contours. The expected limits for the 1$\sigma$ down variation are not shown as no mass points could be excluded. The round excluded area in the top part corresponds to the excluded parameter space in the $(\tilde{W},~\tilde{H})$ model $(M_{2} > |\mu|)$, while the two small areas in at the bottom are that in the $(\tilde{H},~\tilde{W})$ model $(M_{2} < |\mu|)$.

95% CL exclusion limits for the $(\tilde{W},\tilde{H})$ models. The limits are projected on to a two-dimensional plane as a function of the wino/higgsino mass parameters: $(M_2,\mu)$. For the limits shown on the $(M_2,\mu)$ plane, the excluded regions are indicated by the area inside the contours. The expected limits for the 1$\sigma$ down variation are not shown as no mass points could be excluded. The round excluded area in the top part corresponds to the excluded parameter space in the $(\tilde{W},~\tilde{H})$ model $(M_{2} > |\mu|)$, while the two small areas in at the bottom are that in the $(\tilde{H},~\tilde{W})$ model $(M_{2} < |\mu|)$.

95% CL exclusion limits for the $(\tilde{W},\tilde{H})$ models. The limits are projected on to a two-dimensional plane as a function of the wino/higgsino mass parameters: $(M_2,\mu)$. For the limits shown on the $(M_2,\mu)$ plane, the excluded regions are indicated by the area inside the contours. The round excluded area in the top part corresponds to the excluded parameter space in the $(\tilde{W},~\tilde{H})$ model $(M_{2} > |\mu|)$, while the two small areas in at the bottom are that in the $(\tilde{H},~\tilde{W})$ model $(M_{2} < |\mu|)$.

95% CL exclusion limits for the $(\tilde{W},\tilde{H})$ models. The limits are projected on to a two-dimensional plane as a function of the wino/higgsino mass parameters: $(M_2,\mu)$. For the limits shown on the $(M_2,\mu)$ plane, the excluded regions are indicated by the area inside the contours. The expected limits for the 1$\sigma$ down variation are not shown as no mass points could be excluded. The round excluded area in the top part corresponds to the excluded parameter space in the $(\tilde{W},~\tilde{H})$ model $(M_{2} > |\mu|)$, while the two small areas in at the bottom are that in the $(\tilde{H},~\tilde{W})$ model $(M_{2} < |\mu|)$.

95% CL exclusion limits for the $(\tilde{W},\tilde{H})$ models. The limits are projected on to a two-dimensional plane as a function of the wino/higgsino mass parameters: $(M_2,\mu)$. For the limits shown on the $(M_2,\mu)$ plane, the excluded regions are indicated by the area inside the contours. The expected limits for the 1$\sigma$ down variation are not shown as no mass points could be excluded. The round excluded area in the top part corresponds to the excluded parameter space in the $(\tilde{W},~\tilde{H})$ model $(M_{2} > |\mu|)$, while the two small areas in at the bottom are that in the $(\tilde{H},~\tilde{W})$ model $(M_{2} < |\mu|)$.

95% CL exclusion limits for the $(\tilde{H},\tilde{W})$ models. The limits are projected on to a two-dimensional plane as a function of the wino/higgsino mass parameters: $(M_2,\mu)$. For the limits shown on the $(M_2,\mu)$ plane, the excluded regions are indicated by the area inside the contours. The round excluded area in the top part corresponds to the excluded parameter space in the $(\tilde{W},~\tilde{H})$ model $(M_{2} > |\mu|)$, while the two small areas in at the bottom are that in the $(\tilde{H},~\tilde{W})$ model $(M_{2} < |\mu|)$.

95% CL exclusion limits for the $(\tilde{H},\tilde{W})$ models. The limits are projected on to a two-dimensional plane as a function of the wino/higgsino mass parameters: $(M_2,\mu)$. For the limits shown on the $(M_2,\mu)$ plane, the excluded regions are indicated by the area inside the contours. The expected limits for the 1$\sigma$ down variation are not shown as no mass points could be excluded. The round excluded area in the top part corresponds to the excluded parameter space in the $(\tilde{W},~\tilde{H})$ model $(M_{2} > |\mu|)$, while the two small areas in at the bottom are that in the $(\tilde{H},~\tilde{W})$ model $(M_{2} < |\mu|)$.

95% CL exclusion limits for the $(\tilde{H},\tilde{W})$ models. The limits are projected on to a two-dimensional plane as a function of the wino/higgsino mass parameters: $(M_2,\mu)$. For the limits shown on the $(M_2,\mu)$ plane, the excluded regions are indicated by the area inside the contours. The round excluded area in the top part corresponds to the excluded parameter space in the $(\tilde{W},~\tilde{H})$ model $(M_{2} > |\mu|)$, while the two small areas in at the bottom are that in the $(\tilde{H},~\tilde{W})$ model $(M_{2} < |\mu|)$.

95% CL exclusion limits for the $(\tilde{H},\tilde{W})$ models. The limits are projected on to a two-dimensional plane as a function of the wino/higgsino mass parameters: $(M_2,\mu)$. For the limits shown on the $(M_2,\mu)$ plane, the excluded regions are indicated by the area inside the contours. The expected limits for the 1$\sigma$ down variation are not shown as no mass points could be excluded. The round excluded area in the top part corresponds to the excluded parameter space in the $(\tilde{W},~\tilde{H})$ model $(M_{2} > |\mu|)$, while the two small areas in at the bottom are that in the $(\tilde{H},~\tilde{W})$ model $(M_{2} < |\mu|)$.

95% CL exclusion limits for the $(\tilde{H},\tilde{W})$ models. The limits are projected on to a two-dimensional plane as a function of the wino/higgsino mass parameters: $(M_2,\mu)$. For the limits shown on the $(M_2,\mu)$ plane, the excluded regions are indicated by the area inside the contours. The expected limits for the 1$\sigma$ down variation are not shown as no mass points could be excluded. The round excluded area in the top part corresponds to the excluded parameter space in the $(\tilde{W},~\tilde{H})$ model $(M_{2} > |\mu|)$, while the two small areas in at the bottom are that in the $(\tilde{H},~\tilde{W})$ model $(M_{2} < |\mu|)$.

Expected (dashed) and observed (solid red) 95% CL exclusion limit derived for the $(\tilde{H},\tilde{G})$ model, as a function of the lightest higgsino mass ($m(\tilde{\chi}_{1}^{0})$) and the branching ratio $\mathcal{B}(\tilde{\chi}_{1}^{0} \to Z\tilde{G}) (=1-\mathcal{B}(\tilde{\chi}_{1}^{0} \to h\tilde{G}))$. The excluded region is indicated by the area inside the contour.

Expected (dashed) and observed (solid red) 95% CL exclusion limit derived for the $(\tilde{H},\tilde{G})$ model, as a function of the lightest higgsino mass ($m(\tilde{\chi}_{1}^{0})$) and the branching ratio $\mathcal{B}(\tilde{\chi}_{1}^{0} \to Z\tilde{G}) (=1-\mathcal{B}(\tilde{\chi}_{1}^{0} \to h\tilde{G}))$. The excluded region is indicated by the area inside the contour.

Expected (dashed) and observed (solid red) 95% CL exclusion limit derived for the $(\tilde{H},\tilde{G})$ model, as a function of the lightest higgsino mass ($m(\tilde{\chi}_{1}^{0})$) and the branching ratio $\mathcal{B}(\tilde{\chi}_{1}^{0} \to Z\tilde{G}) (=1-\mathcal{B}(\tilde{\chi}_{1}^{0} \to h\tilde{G}))$. The excluded region is indicated by the area inside the contour.

Expected (dashed) and observed (solid red) 95% CL exclusion limit derived for the $(\tilde{H},\tilde{G})$ model, as a function of the lightest higgsino mass ($m(\tilde{\chi}_{1}^{0})$) and the branching ratio $\mathcal{B}(\tilde{\chi}_{1}^{0} \to Z\tilde{G}) (=1-\mathcal{B}(\tilde{\chi}_{1}^{0} \to h\tilde{G}))$. The excluded region is indicated by the area inside the contour.

Expected (dashed) and observed (solid red) 95% CL exclusion limit derived for the $(\tilde{H},\tilde{G})$ model, as a function of the lightest higgsino mass ($m(\tilde{\chi}_{1}^{0})$) and the branching ratio $\mathcal{B}(\tilde{\chi}_{1}^{0} \to Z\tilde{G}) (=1-\mathcal{B}(\tilde{\chi}_{1}^{0} \to h\tilde{G}))$. The excluded region is indicated by the area inside the contour.

Expected (dashed) and observed (solid red) 95% CL exclusion limit derived for the $(\tilde{H},\tilde{G})$ model, as a function of the lightest higgsino mass ($m(\tilde{\chi}_{1}^{0})$) and the branching ratio $\mathcal{B}(\tilde{\chi}_{1}^{0} \to Z\tilde{G}) (=1-\mathcal{B}(\tilde{\chi}_{1}^{0} \to h\tilde{G}))$. The excluded region is indicated by the area inside the contour.

95% CL exclusion limits for the $(\tilde{H},\tilde{a})$ model as the function of axino mass ($m(\tilde{a})$) and the lightest higgsino ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed lines) and observed (solid lines) limits with various $\mathcal{B}(\tilde{\chi}_{1}^{0} \to Z\tilde{a}) (=1-\mathcal{B}(\tilde{\chi}_{1}^{0} \to h\tilde{a}))$ hypotheses. No expected limit is derived for the case with $\mathcal{B}(\tilde{\chi}_{1}^{0} \to Z\tilde{a})=100$% as no mass point on the plane can be excluded.

95% CL exclusion limits for the $(\tilde{H},\tilde{a})$ model as the function of axino mass ($m(\tilde{a})$) and the lightest higgsino ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed line) and observed (solid red line) limits calculated for $\mathcal{B}(\tilde{\chi}_{1}^{0} \to Z\tilde{a})=100$%.

95% CL exclusion limits for the $(\tilde{H},\tilde{a})$ model as the function of axino mass ($m(\tilde{a})$) and the lightest higgsino ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed line) and observed (solid red line) limits calculated for $\mathcal{B}(\tilde{\chi}_{1}^{0} \to Z\tilde{a})=100$%.

95% CL exclusion limits for the $(\tilde{H},\tilde{a})$ model as the function of axino mass ($m(\tilde{a})$) and the lightest higgsino ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed lines) and observed (solid lines) limits with various $\mathcal{B}(\tilde{\chi}_{1}^{0} \to Z\tilde{a}) (=1-\mathcal{B}(\tilde{\chi}_{1}^{0} \to h\tilde{a}))$ hypotheses. No expected limit is derived for the case with $\mathcal{B}(\tilde{\chi}_{1}^{0} \to Z\tilde{a})=100$% as no mass point on the plane can be excluded.

95% CL exclusion limits for the $(\tilde{H},\tilde{a})$ model as the function of axino mass ($m(\tilde{a})$) and the lightest higgsino ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed line) and observed (solid red line) limits calculated for $\mathcal{B}(\tilde{\chi}_{1}^{0} \to Z\tilde{a})=100$%.

95% CL exclusion limits for the $(\tilde{H},\tilde{a})$ model as the function of axino mass ($m(\tilde{a})$) and the lightest higgsino ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed line) and observed (solid red line) limits calculated for $\mathcal{B}(\tilde{\chi}_{1}^{0} \to Z\tilde{a})=100$%.

95% CL exclusion limits for the $(\tilde{H},\tilde{a})$ model as the function of axino mass ($m(\tilde{a})$) and the lightest higgsino ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed lines) and observed (solid lines) limits with various $\mathcal{B}(\tilde{\chi}_{1}^{0} \to Z\tilde{a}) (=1-\mathcal{B}(\tilde{\chi}_{1}^{0} \to h\tilde{a}))$ hypotheses. No expected limit is derived for the case with $\mathcal{B}(\tilde{\chi}_{1}^{0} \to Z\tilde{a})=75$% as no mass point on the plane can be excluded.

95% CL exclusion limits for the $(\tilde{H},\tilde{a})$ model as the function of axino mass ($m(\tilde{a})$) and the lightest higgsino ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed lines) and observed (solid lines) limits with various $\mathcal{B}(\tilde{\chi}_{1}^{0} \to Z\tilde{a}) (=1-\mathcal{B}(\tilde{\chi}_{1}^{0} \to h\tilde{a}))$ hypotheses. No expected limit is derived for the case with $\mathcal{B}(\tilde{\chi}_{1}^{0} \to Z\tilde{a})=75$% as no mass point on the plane can be excluded.

95% CL exclusion limits for the $(\tilde{H},\tilde{a})$ model as the function of axino mass ($m(\tilde{a})$) and the lightest higgsino ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed lines) and observed (solid lines) limits with various $\mathcal{B}(\tilde{\chi}_{1}^{0} \to Z\tilde{a}) (=1-\mathcal{B}(\tilde{\chi}_{1}^{0} \to h\tilde{a}))$ hypotheses. No expected limit is derived for the case with $\mathcal{B}(\tilde{\chi}_{1}^{0} \to Z\tilde{a})=50$% as no mass point on the plane can be excluded.

95% CL exclusion limits for the $(\tilde{H},\tilde{a})$ model as the function of axino mass ($m(\tilde{a})$) and the lightest higgsino ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed lines) and observed (solid lines) limits with various $\mathcal{B}(\tilde{\chi}_{1}^{0} \to Z\tilde{a}) (=1-\mathcal{B}(\tilde{\chi}_{1}^{0} \to h\tilde{a}))$ hypotheses. No expected limit is derived for the case with $\mathcal{B}(\tilde{\chi}_{1}^{0} \to Z\tilde{a})=50$% as no mass point on the plane can be excluded.

95% CL exclusion limits for the $(\tilde{H},\tilde{a})$ model as the function of axino mass ($m(\tilde{a})$) and the lightest higgsino ($m(\tilde{\chi}_{1}^{0})$). Expected (dashed lines) and observed (solid lines) limits with various $\mathcal{B}(\tilde{\chi}_{1}^{0} \to Z\tilde{a}) (=1-\mathcal{B}(\tilde{\chi}_{1}^{0} \to h\tilde{a}))$ hypotheses. No expected limit is derived for the case with $\mathcal{B}(\tilde{\chi}_{1}^{0} \to Z\tilde{a})=25$% as no mass point on the plane can be excluded.

The $\tilde{\chi}_{\textrm{heavy}}$ branching ratios in the $(\tilde{W},\tilde{H})$ model ($\tan\beta=10, \mu>0$). The branching ratios of wino-like chargino ($\tilde{\chi}_{2}^{\pm}$) and wino-like neutralino ($\tilde{\chi}_{3}^{0}$) are shown as the function of the mass of $\tilde{\chi}_{2}^{\pm}$ and that of the lightest higgsino-like neutralino ($\tilde{\chi}_{1}^{0}$), representing $m(\tilde{\chi}_{\textrm{light}})$ and $m(\tilde{\chi}_{\textrm{heavy}})$ respectively. The branchings into the two neutral higgsinos are summed up in the presentation i.e. $\mathcal{B}(\dots \rightarrow \tilde{\chi}_{1,2}^{0}) := \mathcal{B}(\dots \rightarrow \tilde{\chi}_{1}^{0})+\mathcal{B}(\dots \rightarrow \tilde{\chi}_{2}^{0})$.

The $\tilde{\chi}_{\textrm{heavy}}$ branching ratios in the $(\tilde{W},\tilde{H})$ model ($\tan\beta=10, \mu>0$). The branching ratios of wino-like chargino ($\tilde{\chi}_{2}^{\pm}$) and wino-like neutralino ($\tilde{\chi}_{3}^{0}$) are shown as the function of the mass of $\tilde{\chi}_{2}^{\pm}$ and that of the lightest higgsino-like neutralino ($\tilde{\chi}_{1}^{0}$), representing $m(\tilde{\chi}_{\textrm{light}})$ and $m(\tilde{\chi}_{\textrm{heavy}})$ respectively. The branchings into the two neutral higgsinos are summed up in the presentation i.e. $\mathcal{B}(\dots \rightarrow \tilde{\chi}_{1,2}^{0}) := \mathcal{B}(\dots \rightarrow \tilde{\chi}_{1}^{0})+\mathcal{B}(\dots \rightarrow \tilde{\chi}_{2}^{0})$.

The $\tilde{\chi}_{\textrm{heavy}}$ branching ratios in the $(\tilde{W},\tilde{H})$ model ($\tan\beta=10, \mu>0$). The branching ratios of wino-like chargino ($\tilde{\chi}_{2}^{\pm}$) and wino-like neutralino ($\tilde{\chi}_{3}^{0}$) are shown as the function of the mass of $\tilde{\chi}_{2}^{\pm}$ and that of the lightest higgsino-like neutralino ($\tilde{\chi}_{1}^{0}$), representing $m(\tilde{\chi}_{\textrm{light}})$ and $m(\tilde{\chi}_{\textrm{heavy}})$ respectively. The branchings into the two neutral higgsinos are summed up in the presentation i.e. $\mathcal{B}(\dots \rightarrow \tilde{\chi}_{1,2}^{0}) := \mathcal{B}(\dots \rightarrow \tilde{\chi}_{1}^{0})+\mathcal{B}(\dots \rightarrow \tilde{\chi}_{2}^{0})$.

The $\tilde{\chi}_{\textrm{heavy}}$ branching ratios in the $(\tilde{W},\tilde{H})$ model ($\tan\beta=10, \mu>0$). The branching ratios of wino-like chargino ($\tilde{\chi}_{2}^{\pm}$) and wino-like neutralino ($\tilde{\chi}_{3}^{0}$) are shown as the function of the mass of $\tilde{\chi}_{2}^{\pm}$ and that of the lightest higgsino-like neutralino ($\tilde{\chi}_{1}^{0}$), representing $m(\tilde{\chi}_{\textrm{light}})$ and $m(\tilde{\chi}_{\textrm{heavy}})$ respectively. The branchings into the two neutral higgsinos are summed up in the presentation i.e. $\mathcal{B}(\dots \rightarrow \tilde{\chi}_{1,2}^{0}) := \mathcal{B}(\dots \rightarrow \tilde{\chi}_{1}^{0})+\mathcal{B}(\dots \rightarrow \tilde{\chi}_{2}^{0})$.

The $\tilde{\chi}_{\textrm{heavy}}$ branching ratios in the $(\tilde{W},\tilde{H})$ model ($\tan\beta=10, \mu>0$). The branching ratios of wino-like chargino ($\tilde{\chi}_{2}^{\pm}$) and wino-like neutralino ($\tilde{\chi}_{3}^{0}$) are shown as the function of the mass of $\tilde{\chi}_{2}^{\pm}$ and that of the lightest higgsino-like neutralino ($\tilde{\chi}_{1}^{0}$), representing $m(\tilde{\chi}_{\textrm{light}})$ and $m(\tilde{\chi}_{\textrm{heavy}})$ respectively. The branchings into the two neutral higgsinos are summed up in the presentation i.e. $\mathcal{B}(\dots \rightarrow \tilde{\chi}_{1,2}^{0}) := \mathcal{B}(\dots \rightarrow \tilde{\chi}_{1}^{0})+\mathcal{B}(\dots \rightarrow \tilde{\chi}_{2}^{0})$.

The $\tilde{\chi}_{\textrm{heavy}}$ branching ratios in the $(\tilde{W},\tilde{H})$ model ($\tan\beta=10, \mu>0$). The branching ratios of wino-like chargino ($\tilde{\chi}_{2}^{\pm}$) and wino-like neutralino ($\tilde{\chi}_{3}^{0}$) are shown as the function of the mass of $\tilde{\chi}_{2}^{\pm}$ and that of the lightest higgsino-like neutralino ($\tilde{\chi}_{1}^{0}$), representing $m(\tilde{\chi}_{\textrm{light}})$ and $m(\tilde{\chi}_{\textrm{heavy}})$ respectively. The branchings into the two neutral higgsinos are summed up in the presentation i.e. $\mathcal{B}(\dots \rightarrow \tilde{\chi}_{1,2}^{0}) := \mathcal{B}(\dots \rightarrow \tilde{\chi}_{1}^{0})+\mathcal{B}(\dots \rightarrow \tilde{\chi}_{2}^{0})$.

The $\tilde{\chi}_{\textrm{heavy}}$ branching ratios in the $(\tilde{H},\tilde{W})$ model ($\tan\beta=10, \mu>0$). The branching ratios of higgsino-like chargino ($\tilde{\chi}_{2}^{\pm}$) and higgsino-like neutralinos ($\tilde{\chi}_{1}^{0}$, $\tilde{\chi}_{3}^{0}$) as the function of the mass of $\tilde{\chi}_{2}^{\pm}$ and that of wino-like neutralino ($\tilde{\chi}_{1}^{0}$), representing $m(\tilde{\chi}_{\textrm{light}})$ and $m(\tilde{\chi}_{\textrm{heavy}})$ respectively.

The $\tilde{\chi}_{\textrm{heavy}}$ branching ratios in the $(\tilde{H},\tilde{W})$ model ($\tan\beta=10, \mu>0$). The branching ratios of higgsino-like chargino ($\tilde{\chi}_{2}^{\pm}$) and higgsino-like neutralinos ($\tilde{\chi}_{1}^{0}$, $\tilde{\chi}_{3}^{0}$) as the function of the mass of $\tilde{\chi}_{2}^{\pm}$ and that of wino-like neutralino ($\tilde{\chi}_{1}^{0}$), representing $m(\tilde{\chi}_{\textrm{light}})$ and $m(\tilde{\chi}_{\textrm{heavy}})$ respectively.

The $\tilde{\chi}_{\textrm{heavy}}$ branching ratios in the $(\tilde{H},\tilde{W})$ model ($\tan\beta=10, \mu>0$). The branching ratios of higgsino-like chargino ($\tilde{\chi}_{2}^{\pm}$) and higgsino-like neutralinos ($\tilde{\chi}_{1}^{0}$, $\tilde{\chi}_{3}^{0}$) as the function of the mass of $\tilde{\chi}_{2}^{\pm}$ and that of wino-like neutralino ($\tilde{\chi}_{1}^{0}$), representing $m(\tilde{\chi}_{\textrm{light}})$ and $m(\tilde{\chi}_{\textrm{heavy}})$ respectively.

The $\tilde{\chi}_{\textrm{heavy}}$ branching ratios in the $(\tilde{H},\tilde{W})$ model ($\tan\beta=10, \mu>0$). The branching ratios of higgsino-like chargino ($\tilde{\chi}_{2}^{\pm}$) and higgsino-like neutralinos ($\tilde{\chi}_{1}^{0}$, $\tilde{\chi}_{3}^{0}$) as the function of the mass of $\tilde{\chi}_{2}^{\pm}$ and that of wino-like neutralino ($\tilde{\chi}_{1}^{0}$), representing $m(\tilde{\chi}_{\textrm{light}})$ and $m(\tilde{\chi}_{\textrm{heavy}})$ respectively.

The $\tilde{\chi}_{\textrm{heavy}}$ branching ratios in the $(\tilde{H},\tilde{W})$ model ($\tan\beta=10, \mu>0$). The branching ratios of higgsino-like chargino ($\tilde{\chi}_{2}^{\pm}$) and higgsino-like neutralinos ($\tilde{\chi}_{1}^{0}$, $\tilde{\chi}_{3}^{0}$) as the function of the mass of $\tilde{\chi}_{2}^{\pm}$ and that of wino-like neutralino ($\tilde{\chi}_{1}^{0}$), representing $m(\tilde{\chi}_{\textrm{light}})$ and $m(\tilde{\chi}_{\textrm{heavy}})$ respectively.