Upper limit on the WIMP-nucleus scattering cross section from the single-scatter analysis.

Total acceptance evaluated on a simulated dataset of MIMP events. The mean and upper and lower bounds of the 1 sigma uncertainty band are provided.

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.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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.

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.

Isoscalar WIMP-nucleon elastic coupling limit for Operator 9

Isoscalar WIMP-nucleon elastic coupling limit for Operator 12

Isoscalar WIMP-nucleon elastic coupling limit for Operator 7

Isoscalar WIMP-nucleon elastic coupling limit for Operator 13

Isoscalar WIMP-nucleon elastic coupling limit for Operator 3

Isoscalar WIMP-nucleon elastic coupling limit for Operator 5

Isoscalar WIMP-nucleon elastic coupling limit for Operator 11

Isoscalar WIMP-nucleon elastic coupling limit for Operator 14

Isoscalar WIMP-nucleon elastic coupling limit for Operator 4

Isoscalar WIMP-nucleon elastic coupling limit for Operator 10

Isoscalar WIMP-nucleon elastic coupling limit for Operator 1

Isoscalar WIMP-nucleon elastic coupling limit for Operator 6

Isoscalar WIMP-nucleon elastic coupling limit for Operator 15

Isovector WIMP-nucleon elastic coupling limit for Operator 3

Isovector WIMP-nucleon elastic coupling limit for Operator 10

Isovector WIMP-nucleon elastic coupling limit for Operator 13

Isovector WIMP-nucleon elastic coupling limit for Operator 12

Isovector WIMP-nucleon elastic coupling limit for Operator 8

Isovector WIMP-nucleon elastic coupling limit for Operator 1

Isovector WIMP-nucleon elastic coupling limit for Operator 4

Isovector WIMP-nucleon elastic coupling limit for Operator 5

Isovector WIMP-nucleon elastic coupling limit for Operator 7

Isovector WIMP-nucleon elastic coupling limit for Operator 11

Isovector WIMP-nucleon elastic coupling limit for Operator 15

Isovector WIMP-nucleon elastic coupling limit for Operator 9

Isovector WIMP-nucleon elastic coupling limit for Operator 14

Isovector WIMP-nucleon elastic coupling limit for Operator 6

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 12

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 3

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 9

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 7

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 10

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 15

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 8

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 4

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 5

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 6

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 13

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 11

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 14

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 1

Isovector WIMP-nucleon inelastic coupling limit for Operator 13

Isovector WIMP-nucleon inelastic coupling limit for Operator 1

Isovector WIMP-nucleon inelastic coupling limit for Operator 4

Isovector WIMP-nucleon inelastic coupling limit for Operator 7

Isovector WIMP-nucleon inelastic coupling limit for Operator 9

Isovector WIMP-nucleon inelastic coupling limit for Operator 3

Isovector WIMP-nucleon inelastic coupling limit for Operator 11

Isovector WIMP-nucleon inelastic coupling limit for Operator 5

Isovector WIMP-nucleon inelastic coupling limit for Operator 10

Isovector WIMP-nucleon inelastic coupling limit for Operator 8

Isovector WIMP-nucleon inelastic coupling limit for Operator 14

Isovector WIMP-nucleon inelastic coupling limit for Operator 15

Isovector WIMP-nucleon inelastic coupling limit for Operator 6

Isovector WIMP-nucleon inelastic coupling limit for Operator 12

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

Values of $\Delta Z$ for the discovery sensitivity, as defined in the text, as a function of the invariant mass $\textit{m}$. The j+b-jet invariant mass distribution is calculated in the 10 pb AR. Positive percentages indicate improvements in sensitivity. Horizontal dashed lines are drawn at 100% and 200% to guide the eye. The five benchmark BSM models are (1) charged Higgs boson production in association with a top quark, $tbH^{+}$ with $H^{+} \rightarrow t\bar{b}$; (2) a Kaluza-Klein gauge boson, $W_{KK}$, with the SM $W$ boson and a radion $\phi$; (3) a $Z'$ boson decaying to a composite lepton $E$ and $\ell$, with $E \rightarrow Z\ell$; (4) the sequential standard model $W' \rightarrow W Z' \rightarrow \ell\nu q\bar{q}$; (5) a simplified dark-matter model with an axial-vector mediator $Z' \rightarrow q\bar{q}$, where one of the quarks radiates a $W$ boson decaying to $\ell\nu$. The multiple markers shown for the composite-lepton model at the same invariant mass values correspond to different composite lepton ($E$) masses between 0.25 and 3.5 TeV. The center positions of the markers are set to the masses of the corresponding heavy particles.

Values of $\Delta Z$ for the discovery sensitivity, as defined in the text, as a function of the invariant mass $\textit{m}$. The 2b-jet invariant mass distribution is calculated in the 10 pb AR. Positive percentages indicate improvements in sensitivity. Horizontal dashed lines are drawn at 100% and 200% to guide the eye. The five benchmark BSM models are (1) charged Higgs boson production in association with a top quark, $tbH^{+}$ with $H^{+} \rightarrow t\bar{b}$; (2) a Kaluza-Klein gauge boson, $W_{KK}$, with the SM $W$ boson and a radion $\phi$; (3) a $Z'$ boson decaying to a composite lepton $E$ and $\ell$, with $E \rightarrow Z\ell$; (4) the sequential standard model $W' \rightarrow W Z' \rightarrow \ell\nu q\bar{q}$; (5) a simplified dark-matter model with an axial-vector mediator $Z' \rightarrow q\bar{q}$, where one of the quarks radiates a $W$ boson decaying to $\ell\nu$. The multiple markers shown for the composite-lepton model at the same invariant mass values correspond to different composite lepton ($E$) masses between 0.25 and 3.5 TeV. The center positions of the markers are set to the masses of the corresponding heavy particles.

Values of $\Delta Z$ for the discovery sensitivity, as defined in the text, as a function of the invariant mass $\textit{m}$. The j+e invariant mass distribution is calculated in the 10 pb AR. Positive percentages indicate improvements in sensitivity. Horizontal dashed lines are drawn at 100% and 200% to guide the eye. The five benchmark BSM models are (1) charged Higgs boson production in association with a top quark, $tbH^{+}$ with $H^{+} \rightarrow t\bar{b}$; (2) a Kaluza-Klein gauge boson, $W_{KK}$, with the SM $W$ boson and a radion $\phi$; (3) a $Z'$ boson decaying to a composite lepton $E$ and $\ell$, with $E \rightarrow Z\ell$; (4) the sequential standard model $W' \rightarrow W Z' \rightarrow \ell\nu q\bar{q}$; (5) a simplified dark-matter model with an axial-vector mediator $Z' \rightarrow q\bar{q}$, where one of the quarks radiates a $W$ boson decaying to $\ell\nu$. The multiple markers shown for the composite-lepton model at the same invariant mass values correspond to different composite lepton ($E$) masses between 0.25 and 3.5 TeV. The center positions of the markers are set to the masses of the corresponding heavy particles.

Values of $\Delta Z$ for the discovery sensitivity, as defined in the text, as a function of the invariant mass $\textit{m}$. The b-jet+e invariant mass distribution is calculated in the 10 pb AR. Positive percentages indicate improvements in sensitivity. Horizontal dashed lines are drawn at 100% and 200% to guide the eye. The five benchmark BSM models are (1) charged Higgs boson production in association with a top quark, $tbH^{+}$ with $H^{+} \rightarrow t\bar{b}$; (2) a Kaluza-Klein gauge boson, $W_{KK}$, with the SM $W$ boson and a radion $\phi$; (3) a $Z'$ boson decaying to a composite lepton $E$ and $\ell$, with $E \rightarrow Z\ell$; (4) the sequential standard model $W' \rightarrow W Z' \rightarrow \ell\nu q\bar{q}$; (5) a simplified dark-matter model with an axial-vector mediator $Z' \rightarrow q\bar{q}$, where one of the quarks radiates a $W$ boson decaying to $\ell\nu$. The multiple markers shown for the composite-lepton model at the same invariant mass values correspond to different composite lepton ($E$) masses between 0.25 and 3.5 TeV. The center positions of the markers are set to the masses of the corresponding heavy particles.

Values of $\Delta Z$ for the discovery sensitivity, as defined in the text, as a function of the invariant mass $\textit{m}$. The j+$\gamma$ invariant mass distribution is calculated in the 10 pb AR. Positive percentages indicate improvements in sensitivity. Horizontal dashed lines are drawn at 100% and 200% to guide the eye. The five benchmark BSM models are (1) charged Higgs boson production in association with a top quark, $tbH^{+}$ with $H^{+} \rightarrow t\bar{b}$; (2) a Kaluza-Klein gauge boson, $W_{KK}$, with the SM $W$ boson and a radion $\phi$; (3) a $Z'$ boson decaying to a composite lepton $E$ and $\ell$, with $E \rightarrow Z\ell$; (4) the sequential standard model $W' \rightarrow W Z' \rightarrow \ell\nu q\bar{q}$; (5) a simplified dark-matter model with an axial-vector mediator $Z' \rightarrow q\bar{q}$, where one of the quarks radiates a $W$ boson decaying to $\ell\nu$. The multiple markers shown for the composite-lepton model at the same invariant mass values correspond to different composite lepton ($E$) masses between 0.25 and 3.5 TeV. The center positions of the markers are set to the masses of the corresponding heavy particles.

Values of $\Delta Z$ for the discovery sensitivity, as defined in the text, as a function of the invariant mass $\textit{m}$. The j+$\mu$ invariant mass distribution is calculated in the 10 pb AR. Positive percentages indicate improvements in sensitivity. Horizontal dashed lines are drawn at 100% and 200% to guide the eye. The five benchmark BSM models are (1) charged Higgs boson production in association with a top quark, $tbH^{+}$ with $H^{+} \rightarrow t\bar{b}$; (2) a Kaluza-Klein gauge boson, $W_{KK}$, with the SM $W$ boson and a radion $\phi$; (3) a $Z'$ boson decaying to a composite lepton $E$ and $\ell$, with $E \rightarrow Z\ell$; (4) the sequential standard model $W' \rightarrow W Z' \rightarrow \ell\nu q\bar{q}$; (5) a simplified dark-matter model with an axial-vector mediator $Z' \rightarrow q\bar{q}$, where one of the quarks radiates a $W$ boson decaying to $\ell\nu$. The multiple markers shown for the composite-lepton model at the same invariant mass values correspond to different composite lepton ($E$) masses between 0.25 and 3.5 TeV. The center positions of the markers are set to the masses of the corresponding heavy particles.

Values of $\Delta Z$ for the discovery sensitivity, as defined in the text, as a function of the invariant mass $\textit{m}$. The b-jet+$\mu$ invariant mass distribution is calculated in the 10 pb AR. Positive percentages indicate improvements in sensitivity. Horizontal dashed lines are drawn at 100% and 200% to guide the eye. The five benchmark BSM models are (1) charged Higgs boson production in association with a top quark, $tbH^{+}$ with $H^{+} \rightarrow t\bar{b}$; (2) a Kaluza-Klein gauge boson, $W_{KK}$, with the SM $W$ boson and a radion $\phi$; (3) a $Z'$ boson decaying to a composite lepton $E$ and $\ell$, with $E \rightarrow Z\ell$; (4) the sequential standard model $W' \rightarrow W Z' \rightarrow \ell\nu q\bar{q}$; (5) a simplified dark-matter model with an axial-vector mediator $Z' \rightarrow q\bar{q}$, where one of the quarks radiates a $W$ boson decaying to $\ell\nu$. The multiple markers shown for the composite-lepton model at the same invariant mass values correspond to different composite lepton ($E$) masses between 0.25 and 3.5 TeV. The center positions of the markers are set to the masses of the corresponding heavy particles.

Values of $\Delta Z$ for the discovery sensitivity, as defined in the text, as a function of the invariant mass $\textit{m}$. The b-jet+$\gamma$ invariant mass distribution is calculated in the 10 pb AR. Positive percentages indicate improvements in sensitivity. Horizontal dashed lines are drawn at 100% and 200% to guide the eye. The five benchmark BSM models are (1) charged Higgs boson production in association with a top quark, $tbH^{+}$ with $H^{+} \rightarrow t\bar{b}$; (2) a Kaluza-Klein gauge boson, $W_{KK}$, with the SM $W$ boson and a radion $\phi$; (3) a $Z'$ boson decaying to a composite lepton $E$ and $\ell$, with $E \rightarrow Z\ell$; (4) the sequential standard model $W' \rightarrow W Z' \rightarrow \ell\nu q\bar{q}$; (5) a simplified dark-matter model with an axial-vector mediator $Z' \rightarrow q\bar{q}$, where one of the quarks radiates a $W$ boson decaying to $\ell\nu$. The multiple markers shown for the composite-lepton model at the same invariant mass values correspond to different composite lepton ($E$) masses between 0.25 and 3.5 TeV. The center positions of the markers are set to the masses of the corresponding heavy particles.

The 95% CL upper limits on the cross section times acceptance ($A$), efficiency ($\epsilon$), and branching ratio ($B$) for Gaussian-shaped signals with different signal widths. The limits are calculated for events with at least one lepton with pT > 60 GeV in the 10 pb AR. Two width hypotheses are shown; $\sigma X / mX = 0$ and $\sigma X / mX = 0.15$. In both cases, the detector resolution for jets is included in the simulation of signal samples. The $\pm1 \sigma$ and $\pm2 \sigma$ bands around the expected limit are shown for $\sigma X / mX = 0$ signals. Mass points are spaced 5% apart, relative to the preceding point, starting at 0.3 TeV.

Distributions of the anomaly score from the AE for data and five benchmark BSM models. Their legends, from top to bottom, are; (1) charged Higgs boson production in association with a top quark, $tbH^{+}$ with $H^{+} \rightarrow t\bar{b}$, with the $H^{+}$ mass of 2 TeV; (2) a 6 TeV Kaluza-Klein gauge boson, $W_{KK}$, with the SM $W$ boson and a radion $\phi$; (3) a 4 TeV $Z'$ boson decaying to a composite lepton $E$ with a mass of 0.5 TeV; (4) the sequential standard model $W' \rightarrow WZ' \rightarrow \ell\nu q\bar{q}$, with the $W'$ mass of 6.2 TeV and the $Z'$ mass of 6 TeV; (5) a simplified dark-matter model with an axial-vector mediator $Z' \rightarrow q\bar{q}$, where one of the quarks radiates a $W$ boson decaying to $\ell\nu$ with a mass of 6 TeV. The distributions for the BSM models are smoothed to remove fluctuations due to low MC event counts. The vertical lines indicate the start of the three anomaly regions (ARs). The labels of the three ARs indicate the visible cross section for hypothetical processes yielding the same number of events as observed in the 140 $fb^{-1}$ dataset. The AE is applied to preselected events without any requirements on invariant mass distributions.

Invariant mass distributions of $j+Y$ for $M_{jY}$ > 0.3 TeV after preselection along with the fit from Eq.(1). The fit is represented by red lines, and the associated uncertainties are indicated by yellow bands. The lower panels show the bin-by-bin significances of deviations from the fit, calculated as $(d_{\textit{i}} - f_{\textit{i}}) / \delta_{\textit{i}}$, where $d_{\textit{i}}$ is the data yield, $f_{\textit{i}}$ is the fit value, and $\delta_{i}$ is the uncertainty of data in the $\textit{i}$-th bin

Distributions of the anomaly score for data and several anomaly scenarios. The example BSM model (shown with the dashed blue lines) is the sequential standard model $W' \rightarrow WZ' \rightarrow \ell\nu q\bar{q}$. The mass of $W'$ is set to 2.2 TeV and the mass of the $Z'$ is set to 2 TeV. This model leads to the final state of one lepton, two jets, and small missing transverse energy that is similar to the SM backgrounds. The other histograms represent events from the same model with artificial modifications to represent anomalous events; (1) "Anomaly 1" is the case where all jets beyond the second jet ($N_{j}$ > 2) are replaced with photons; (2) "Anomaly 2" is the case where all jets beyond the second jet are replaced with $b$-jets; (3) "Anomaly 3" is the case of low-multiplicity events where all jets beyond the second jet are removed. The histograms are normalized to unit. The left peak near log(Loss) = -9 visible in Anomaly 1 and 2 is from the events without a third jet.

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

The observed 90% C.L upper limit on neutrino effective millicharge (\delta_{Q}[e_{0}]) in SR1. The data contains observed upper limit, median sensitivity and 1\sigma and 2\sigma sensitivity range.

The observerd 90% C.L upper limit on solar axion g_{ae} coupling constant in valid across the mass range for SR1. The data contains observed upper limit, median sensitivity and 1\sigma and 2\sigma sensitivity range as a function of Solar Axion Mass .

The observed 90% C.L. upper limit on axion like particle (ALP) coupling constant (g_{ae}) in SR1. ALP mass spanning from 1 keV to 17 keV is considered. The data contains observed upper limit, median sensitivity and 1\sigma and 2\sigma sensitivity range as a function of ALP Mass.

The observed 90% C.L. upper limit on hidden photon (HP) coupling constant squared (\kappa^2) in SR1. HP mass spanning from 1 keV to 17 keV is considered. The data contains observed upper limit, median sensitivity and 1\sigma and 2\sigma sensitivity range as a function of HP Mass.

The observed 90% C.L. upper limit for the spin-independent (cm^2) WIMP cross-section vs. WIMP mass (GeV/c^2), using the Migdal signal pathway. The result is obtained for WIMPs undergoing the Migdal effect from 0.5 to 9 GeV/c2. The data contains observed upper limit, median sensitivity and 1\sigma and 2\sigma sensitivity range as a function of WIMP Mass.

The observed 90% C.L. upper limit for the spin-dependent WIMP cross-section (cm^2) vs. WIMP mass (GeV/c^2) for coupling on neutrons. The result is obtained for WIMPs undergoing the Migdal effect from 0.5 to 9 GeV/c2. The data contains observed upper limit, median sensitivity and 1\sigma and 2\sigma sensitivity range as a function of WIMP Mass.

The observed 90% C.L. upper limit for the spin-dependent WIMP cross-section (cm^2) vs. WIMP mass (GeV/c^2) for coupling on protons. The result is obtained for WIMPs undergoing the Migdal effect from 0.5 to 9 GeV/c2. The data contains observed upper limit, median sensitivity and 1\sigma and 2\sigma sensitivity range as a function of WIMP Mass.

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.

Full experimental (statistical and systematic) correlations (in \%) of the normalized partial decay rates over the $\overline{B}^0\to D^{*+} e^- \bar\nu_e$ and $\overline{B}^0\to D^{*+} \mu^- \bar\nu_\mu$ decays.

Measured partial decay rates $\Delta\Gamma$ (in units of $10^{-15}$ GeV) using the matrix inversion unfolding method

Average of normalized decay rates that are determined using the matrix inversion unfolding method

Full experimental (statistical and systematic) correlations (in \%) of the partial decay rates determined with the matrix inversion unfolding method for the $\overline{B}^0\to D^{*+} e^- \bar\nu_e$ and $\overline{B}^0\to D^{*+} \mu^- \bar\nu_\mu$ decays.

Full experimental (statistical and systematic) correlations (in \%) of the normalized partial decay rates determined using the matrix inversion unfolding method

Summary of results for cross-section of non-prompt $\psi(2S)$ decaying to a muon pair for 13 TeV data in fb/GeV. Uncertainties are statistical and systematic, respectively.

Summary of results for non-prompt fraction of $J/\psi$ decaying to a muon pair as a percentage for 13 TeV data. Uncertainties are statistical and systematic, respectively.

Summary of results for non-prompt fraction of $\psi(2S)$ decaying to a muon pair as a percentage for 13 TeV data. Uncertainties are statistical and systematic, respectively.

Summary of results for prompt $\psi(2S)$ over $J/\psi$ ratio as a percentage for 13 TeV data. Uncertainties are statistical and systematic, respectively.

Summary of results for non-prompt $\psi(2S)$ over $J/\psi$ ratio as a percentage for 13 TeV data. Uncertainties are statistical and systematic, respectively.

Correction factors for an assumed angular dependence $\propto 1+\lambda_{\theta} \cos^2\theta $ in the helicity frame, for several values of the parameter $\lambda_{\theta}$. The correction factors were found to be the same (within 1% - 2%) for $J\psi$ and $\psi(2S)$ mesons, for prompt and non-prompt production mechanisms, and for the three rapidity intervals considered in this paper.

Data selection efficiency as a function of nuclear recoil energy

Data points used in analysis in log_10(S2)-S1 space

The average tetrajet invariant mass distributions in data, along with the fitted background estimates for 0.16 < $\alpha$ < 0.18.

The average tetrajet invariant mass distributions in data, along with the fitted background estimates for 0.18 < $\alpha$ < 0.20.

The average tetrajet invariant mass distributions in data, along with the fitted background estimates for 0.20 < $\alpha$ < 0.22.

The average tetrajet invariant mass distributions in data, along with the fitted background estimates for 0.22 < $\alpha$ < 0.24.

The average tetrajet invariant mass distributions in data, along with the fitted background estimates for 0.24 < $\alpha$ < 0.26.

The average tetrajet invariant mass distributions in data, along with the fitted background estimates for 0.26 < $\alpha$ < 0.28.

The average tetrajet invariant mass distributions in data, along with the fitted background estimates for 0.28 < $\alpha$ < 0.30.

The average tetrajet invariant mass distributions in data, along with the fitted background estimates for 0.30 < $\alpha$ < 0.32.

The average tetrajet invariant mass distributions in data, along with the fitted background estimates for 0.32 < $\alpha$ < 0.34.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.10 < $\alpha$ < 0.12.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.12 < $\alpha$ < 0.14.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.14 < $\alpha$ < 0.16.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.16 < $\alpha$ < 0.18.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.18 < $\alpha$ < 0.20.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.20 < $\alpha$ < 0.22.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.22 < $\alpha$ < 0.24.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.24 < $\alpha$ < 0.26.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.26 < $\alpha$ < 0.28.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.28 < $\alpha$ < 0.30.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.30 < $\alpha$ < 0.32.

The average dijet invariant mass distributions in data, along with the fitted background estimates for 0.32 < $\alpha$ < 0.34.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.10 < $\alpha$ < 0.12.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.12 < $\alpha$ < 0.14.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.14 < $\alpha$ < 0.16.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.16 < $\alpha$ < 0.18.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.18 < $\alpha$ < 0.20.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.20 < $\alpha$ < 0.22.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.22 < $\alpha$ < 0.24.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.24 < $\alpha$ < 0.26.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.26 < $\alpha$ < 0.28.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.28 < $\alpha$ < 0.30.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.30 < $\alpha$ < 0.32.

The expected and observed limits on $m_{4j}$ for signal templates using a 4-parameter fit function for 0.32 < $\alpha$ < 0.34.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.10 < $\alpha$ < 0.12.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.12 < $\alpha$ < 0.14.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.14 < $\alpha$ < 0.16.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.16 < $\alpha$ < 0.18.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.18 < $\alpha$ < 0.20.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.20 < $\alpha$ < 0.22.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.22 < $\alpha$ < 0.24.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.24 < $\alpha$ < 0.26.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.26 < $\alpha$ < 0.28.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.28 < $\alpha$ < 0.30.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.30 < $\alpha$ < 0.32.

The expected and observed limits on $m_{\langle 2j \rangle}$ for signal templates using a 4-parameter fit function for 0.32 < $\alpha$ < 0.34.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.10 < $\alpha$ < 0.12.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.12 < $\alpha$ < 0.14.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.14 < $\alpha$ < 0.16.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.16 < $\alpha$ < 0.18.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.18 < $\alpha$ < 0.20.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.20 < $\alpha$ < 0.22.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.22 < $\alpha$ < 0.24.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.24 < $\alpha$ < 0.26.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.26 < $\alpha$ < 0.28.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.28 < $\alpha$ < 0.30.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.30 < $\alpha$ < 0.32.

The expected and observed limits on $m_{4j}$ for Gaussian signal templates using a 4-parameter fit function for 0.32 < $\alpha$ < 0.34.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.10 < $\alpha$ < 0.12.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.12 < $\alpha$ < 0.14.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.14 < $\alpha$ < 0.16.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.16 < $\alpha$ < 0.18.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.18 < $\alpha$ < 0.20.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.20 < $\alpha$ < 0.22.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.22 < $\alpha$ < 0.24.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.24 < $\alpha$ < 0.26.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.26 < $\alpha$ < 0.28.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.28 < $\alpha$ < 0.30.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.30 < $\alpha$ < 0.32.

The expected and observed limits on $m_{\langle 2j \rangle}$ for Gaussian signal templates using a 4-parameter fit function for 0.32 < $\alpha$ < 0.34.

The product of the analysis acceptance and selection efficiency for all analysis selection criteria is shown as a function of $m_Y$ and $m_{X}/m_{Y}$.

The cross-section of the signal samples is shown as a function of $m_Y$ and $m_{X}/m_{Y}$.

Signal region detector-level distribution for the observable $p_{\mathrm{T}}^{\mathrm{lead}\, \ell}$.

Signal region detector-level distribution for the observable $m_{e\mu}$.

Signal region detector-level distribution for the observable $p_{\mathrm{T}}^{e\mu}$.

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|y_{e\mu}|$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|y_{e\mu}|$

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|\Delta \phi(e \mu)|$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|\Delta \phi(e \mu)|$

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $ \cos\theta^{\ast}$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $ \cos\theta^{\ast}$

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{\mathrm{lead}\, \ell}$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{\mathrm{lead}\, \ell}$

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $m_{e\mu}$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $m_{e\mu}$

Measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{e\mu}$

Relative systematic uncertainties for the fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{e\mu}$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|y_{e\mu}|$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|y_{e\mu}|$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|\Delta \phi(e \mu)|$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $|\Delta \phi(e \mu)|$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $ \cos\theta^{\ast}$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $ \cos\theta^{\ast}$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{\mathrm{lead}\, \ell}$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{\mathrm{lead}\, \ell}$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $m_{e\mu}$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $m_{e\mu}$

The statistical correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{e\mu}$

The total correlation coefficients (in percentage) between bins for the measured fiducial differential cross-section of $WW \rightarrow e^{\pm}\nu\mu^{\mp}\nu$ production for the observable $p_{\mathrm{T}}^{e\mu}$

High $p_T$ ($p_T$ > 6 GeV/$c$) integrated $v_2$ vs $N_{part}$ for (a) $\pi^0$ (b) inclusive photon, and (c) direct photon. Results are shown with both reaction plane detectors$:$ (solid black circles) BBC and (solid red squares) RXN. Each point represents a 10% wide centrality bin from 60%-0%.

Inclusive unlike-sign $\mu \mu$ pair yield from $p$+$p$ collisions at $\sqrt{s}$ = 200 GeV as a function of mass in different $p_T$ slices for the south and north muon arms and the ratio of data to expected sources.

Inclusive unlike-sign $\mu \mu$ pair yield from $p$+$p$ collisions at $\sqrt{s}$ = 200 GeV as a function of mass in different $p_T$ slices for the south and north muon arms and the ratio of data to expected sources.

Inclusive like-sign $\mu \mu$ pair yield from $p$+$p$ collisions at $\sqrt{s}$ = 200 GeV as a function of mass in different $p_T$ slices for the south and north muon arms and the ratio of data to expected sources.

Inclusive like-sign $\mu \mu$ pair yield from $p$+$p$ collisions at $\sqrt{s}$ = 200 GeV as a function of mass in different $p_T$ slices for the south and north muon arms and the ratio of data to expected sources.

Inclusive like-sign $\mu \mu$ pair yield from $p$+$p$ collisions at $\sqrt{s}$ = 200 GeV as a function of mass in different $p_T$ slices for the south and north muon arms and the ratio of data to expected sources.

The $\mu \mu$ pair data in the charm mass region as a function of $\Delta \phi$.

The $\mu \mu$ pair data in the charm mass region as a function of $\Delta \phi.$

The $\mu \mu$ pair data in the charm mass region as a function of pair $p_T$.

The $\mu \mu$ pair data in the charm mass region as a function of pair $p_T$

The like-sign $\mu \mu$ pair data in the bottom mass region as a function of $\Delta \phi$.

The like-sign $\mu \mu$ pair data in the bottom mass region as a function of $\Delta \phi$.

The like-sign $\mu \mu$ pair data in the bottom mass region as a function of pair $p_T$.

The like-sign $\mu \mu$ pair data in the bottom mass region as a function of pair $p_T$.

The unlike-sign $\mu \mu$ pair data used to determine the Drell-Yan contribution as a function of mass.

The unlike-sign $\mu \mu$ pair data used to determine the Drell-Yan contribution as a function of mass.

The unlike-sign $\mu \mu$ pair data used to determine the Drell-Yan contribution as a function of pair $p_T$.

The unlike-sign $\mu \mu$ pair data used to determine the Drell-Yan contribution as a function of pair $p_T$.

Summary of arm-averaged relative systematic uncertainties for the total bottom cross section $\sigma _{b\bar{b}}$, the differential Drell-Yan cross section $d^{2}\sigma _{DY\implies \mu\mu}/dmdy$, and the $b\bar{b}$ ($c\bar{c}$) differential yields $dN_{b\bar{b}(c\bar{c})\implies \mu\mu}/d\Delta\phi$. The systematic uncertainty type is indicated in the second column and is applicable only to the differential measurements. The uncertainties for the differential measurements vary with azimuthal opening angle, pair $p_{T}$ , or mass. Asymmetric uncertainties are quoted in bracketed values. For the $c\bar{c}$ measurement, the regions $\Delta\phi$ < $\pi$/2, $p_{T}$ < 0.5 GeV/$c$ and $p_{T}$ > 2.0 GeV/$c$ are excluded because the yield approaches zero and relative systematic uncertainties diverge. With these regions excluded, the difference between the systematic uncertainties of all measurements for the south and north muon arms differs by no greater than 2% for all systematic uncertainties sources.

The differential yield $dN/d\phi$ of unlike-sign muon pairs from charm with mass 1.5-2.5 GeV/$c^2$ in the ideal muon arm acceptance, as a function of the pair azimuthal opening angle. Statistical (stat.) and types B and C systematic errors (sys. B & sys. C) are all $10^{-9}\times$[rad$^{-1}$].

The differential yield $dN/d\phi$ of like-sign muon pairs from bottom with mass 3.5-10.0 GeV/$c^2$ in the ideal muon arm acceptance, as a function of the pair azimuthal opening angle. Statistical (stat.) and types B and C systematic errors (sys. B & sys. C) are all $10^{-9}\times$[rad$^{-1}$].

The differential yield $dN/dp_{T}$ of unlike-sign muon pairs from bottom with mass 1.5-2.5 GeV/$c^2$ in the ideal muon arm acceptance, as a function of the pair transverse momentum. Statistical (stat.) and types B and C systematic errors (sys. B & sys. C) are all $10^{-9}\times$[$c$/GeV].

The differential yield $dN/dp_{T}$ of like-sign muon pairs from bottom with mass 3.5-10.0 GeV/$c^2$ in the ideal muon arm acceptance, as a function of the pair transverse momentum. Statistical (stat.) and types B and C systematic errors (sys. B & sys. C) are all $10^{-9}\times$[$c$/GeV].

$\sigma_{b\bar{b}}$ from fit using different models. Only statistical uncertainties are shown.

The differential Drell-Yan cross section $\frac{d^{2}\sigma}{dmdy}$ as a function of the muon pair mass, where the muon pair rapidity |$y\mu\mu$| is between 1.2 and 2.2. Statistical (stat.) and types B and C systematic errors (sys. B & sys. C) are all $10^{-9}\times$[pb $c^2$/GeV]. In the data table below, the value of mass correspond to the center of the stated range.

The differential Drell-Yan cross section $\frac{1}{2\pi p_{T}} \frac{d^{2}\sigma}{dydp_{T}}$ as a function of the muon pair transverse momentum, where the muon pair mass $m_{\mu\mu}$ is between 4.8 and 8.2 GeV/$c^2$ and the muon pair rapidity |$y\mu\mu$| is between 1.2 and 2.2. Statistical (stat.) and types B and C systematic errors (sys. B & sys. C) are all $10^{-9}\times$[pb $c^2$/GeV]. In the data table below, the value of p_T correspond to the center of the stated range.

The relative nuclear-modification factor for U+U and Au+Au as a function of collision centrality.

The ratio of the invariant yield for U+U to that for Au+Au, as a function of collision centrality.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Expected 95% CL exclusion contours in the $(m(W_{R}), m(N_{R}))$ plane in the muon channel for boosted.

Observed 95% CL exclusion contours for the assumption of Dirac type in the $(m(W_{R}), m(N_{R}))$ plane in the electron channel for boosted.

Expected 95% CL exclusion contours for the assumption of Dirac type in the $(m(W_{R}), m(N_{R}))$ plane in the electron channel for boosted.

Observed 95% CL exclusion contours for the assumpition of Dirac type in the $(m(W_{R}), m(N_{R}))$ plane in the muon channel for boosted.

Expected 95% CL exclusion contours for the assumption of Dirac type in the $(m(W_{R}), m(N_{R}))$ plane in the muon channel for boosted.

Observed 95% CL exclusion contours for the assumption of Majorana type in the $(m(W_{R}), m(N_{R}))$ plane in the electron channel for resolved.

Expected 95% CL exclusion contours for the assumption of Majorana type in the $(m(W_{R}), m(N_{R}))$ plane in the electron channel for resolved.

Observed 95% CL exclusion contours for the assumption of Majorana type in the $(m(W_{R}), m(N_{R}))$ plane in the muon channel for resolved.

Expected 95% CL exclusion contours for the assumption of Majorana type in the $(m(W_{R}), m(N_{R}))$ plane in the muon channel for resolved.

Observed 95% CL exclusion contours for the assumption of Dirac type in the $(m(W_{R}), m(N_{R}))$ plane in the electron channel for resolved.

Expected 95% CL exclusion contours for the assumption of Dirac type in the $(m(W_{R}), m(N_{R}))$ plane in the electron channel for resolved.

Observed 95% CL exclusion contours for the assumption of Dirac type in the $(m(W_{R}), m(N_{R}))$ plane in the muon channel for resolved.

Expected 95% CL exclusion contours for the assumption of Dirac type in the $(m(W_{R}), m(N_{R}))$ plane in the muon channel for resolved.

The $m_{eJ}$ distribution in the boosted one electron channel.

The $m_{eeJ}$ distribution in the boosted two electron channel.

The $m_{\mu\mu J}$ distribution in the boosted two muon channel.

The $m_{jj}$ distribution in the resolved electron channel.

The $m_{jj}$ distribution in the resolved muon channel.

The $m_{eejj}$ distribution in the resolved electron channel.

The $m_{\mu\mu jj}$ distribution in the resolved muon channel.

The $h_{T}$ distribution in the resolved electron channel.

The $h_{T}$ distribution in the resolved muon channel.

Acceptance $\times$ efficiency for bSR1e.

Acceptance $\times$ efficiency for bSR2e.

Acceptance $\times$ efficiency for bSR2mu.

Acceptance $\times$ efficiency for rSROS2e.

Acceptance $\times$ efficiency for rSROS2mu.

Acceptance $\times$ efficiency for rSRSS2e.

Acceptance $\times$ efficiency for rSRSS2mu.

Expected and observed 95% CL upper limits for $\sigma(pp\to W_{R})\times B(W_{R}\to eeqq)$ as a function of $m(W_{R})$ for the special case of $m(W_{R}) = 2\times m(N_{R})$. The theoretical cross section corresponds to either the Dirac or Majorana neutrino interpretation by ignoring the charge of the leptons.

Expected and observed 95% CL upper limits for $\sigma(pp\to W_{R})\times B(W_{R} \to \mu\mu qq)$ as a function of $m(W_{R})$ for the special case of $m(W_{R}) = 2\times m(N_{R})$. The theoretical cross section corresponds to either the Dirac or Majorana neutrino interpretation by ignoring the charge of the leptons.

Expected and observed 95% CL upper limits for $\sigma(pp\to W_{R})\times B(W_{R} \to eeqq)$ as a function of $m(W_{R})$ for the special case of $m(N_{R})=200$ GeV. The theoretical cross section corresponds to either the Dirac or Majorana neutrino interpretation by ignoring the charge of the leptons.

Expected and observed 95% CL upper limits for $\sigma(pp\to W_{R})\times B(W_{R} \to \mu \mu qq)$ as a function of $m(W_{R})$ for the special case of $m(N_{R})=200$ GeV. The theoretical cross section corresponds to either the Dirac or Majorana neutrino interpretation by ignoring the charge of the leptons.

Expected and observed 95% CL upper limits for $\sigma(pp \to W_{R}) \times B(W_{R} \to eeqq)$ as a function of $m(W_{R})$ for the Majorana neutrino interpretations for the special case of $m(W_{R}) = 2\times m(N_{R})$.

Expected and observed 95% CL upper limits for $\sigma(pp \to W_{R}) \times B(W_{R} \to \mu\mu qq)$ as a function of $m(W_{R})$ for the Majorana neutrino interpretations for the special case of $m(W_{R}) = 2 \times m(N_{R})$.

Expected and observed 95% CL upper limits for $\sigma(pp \to W_{R}) \times B(W_{R} \to eeqq)$ as a function of $m(W_{R})$ for the Dirac neutrino interpretations for the special case of $m(W_{R}) = 2 \times m(N_{R})$.

Expected and observed 95% CL upper limits for $\sigma(pp \to W_{R}) \times B(W_{R} \to \mu \mu qq)$ as a function of $m(W_{R})$ for the Dirac neutrino interpretations for the special case of $m(W_{R}) = 2 \times m(N_{R})$.

The values of $\langle \mathrm{d}N_\mathrm{ch} / \mathrm{d}\eta \rangle$ averaged over $|\eta|<0.5$ for the INEL, NSD, and INEL>0 event classes at $\sqrt{s} = 5.02$ TeV

Pseudorapidity density distributions of charged particles, $\mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta$, in pp collisions at $\sqrt{s} = $ 5.02 TeV for the INEL>0 with momentum threshold $p_{T}>0.15 \mathrm{GeV}/c$ ($\mathrm{INEL>0}_{p_{\mathrm{T}}>0.15}^{|\eta|<0.8}$).

Pseudorapidity density distributions of charged particles, $\mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta$, in pp collisions at $\sqrt{s} = $ 5.02 TeV for the INEL>0 with momentum threshold $p_{T}>0.5 \mathrm{GeV}/c$ ($\mathrm{INEL>0}_{p_{\mathrm{T}}>0.5}^{|\eta|<0.8}$).

Pseudorapidity density distributions of charged particles, $\mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta$, in pp collisions at $\sqrt{s} = $ 5.02 TeV for the INEL>0 with momentum threshold $p_{T}>1 \mathrm{GeV}/c$ ($\mathrm{INEL>0}_{p_{\mathrm{T}}>1}^{|\eta|<0.8}$).

Pseudorapidity density distributions of charged particles, $\mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta$, in pp collisions at $\sqrt{s} = $ 5.02 TeV for the INEL>0 with momentum threshold $p_{T}>2 \mathrm{GeV}/c$ ($\mathrm{INEL>0}_{p_{\mathrm{T}}>2}^{|\eta|<0.8}$).

Pseudorapidity density distributions of charged particles, $\mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta$, in pp collisions at $\sqrt{s} = $ 13 TeV for the INEL>0 with momentum threshold $p_{T}>0.15 \mathrm{GeV}/c$ ($\mathrm{INEL>0}_{p_{\mathrm{T}}>0.15}^{|\eta|<0.8}$).

Pseudorapidity density distributions of charged particles, $\mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta$, in pp collisions at $\sqrt{s} = $ 13 TeV for the INEL>0 with momentum threshold $p_{T}>0.5 \mathrm{GeV}/c$ ($\mathrm{INEL>0}_{p_{\mathrm{T}}>0.5}^{|\eta|<0.8}$).

Pseudorapidity density distributions of charged particles, $\mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta$, in pp collisions at $\sqrt{s} = $ 13 TeV for the INEL>0 with momentum threshold $p_{T}>1 \mathrm{GeV}/c$ ($\mathrm{INEL>0}_{p_{\mathrm{T}}>1}^{|\eta|<0.8}$).

Pseudorapidity density distributions of charged particles, $\mathrm{d}N_{\mathrm{ch}}/\mathrm{d}\eta$, in pp collisions at $\sqrt{s} = $ 13 TeV for the INEL>0 with momentum threshold $p_{T}>2 \mathrm{GeV}/c$ ($\mathrm{INEL>0}_{p_{\mathrm{T}}>2}^{|\eta|<0.8}$).

Lévy scale parameter $R$ versus average $m_T$ of the pair. The graphical representation of statistical and systematic uncertainties is the same as in Fig. 4.

Lévy index parameter $\alpha$ versus average $m_T$ of the pair. Statistical and systematic uncertainties are indicated similarly to Fig. 4. The horizontal line, %\aplha$ = 1.207, represents the 0%-30% centrality average value of %\aplha$.

Inverse square of the Lévy scale parameter 1/$R^2$ versus average $m_T$ of the pair. Statistical and systematic uncertainties shown as bars and boxes, respectively.

New scale parameter $\hat{R}$ versus average $m_T$ of the pair, with a linear fit. Statistical and systematic uncertainties shown as bars and boxes, respectively.

Normalized correlation strength parameter $\lambda/\lambda_{max}$ versus average $m_T$ of the pair. The data are compared with parameter scans from Refs. [58, 59] using different values of in-medium $\eta^{\prime}$ mass $m^{*}_{\eta^{\prime}}$ and slope parameter $B^{-1}_{\eta^{\prime}}$ . A best fit with Eq. (63) and the resulting $H$ and $\sigma$ parameters are also shown.

Normalized correlation strength parameter $\lambda/\lambda_{max}$ versus average $m_T$ of the pair. The data are compared with parameter scans from Refs. [58, 59] using different values of in-medium $\eta^{\prime}$ mass $m^{*}_{\eta^{\prime}}$ and slope parameter $B^{-1}_{\eta^{\prime}}$ . A best fit with Eq. (63) and the resulting $H$ and $\sigma$ parameters are also shown.

Normalized correlation strength parameter $\lambda/\lambda_{max}$ versus average $m_T$ of the pair. The data are compared with parameter scans from Refs. [58, 59] using different values of in-medium $\eta^{\prime}$ mass $m^{*}_{\eta^{\prime}}$ and slope parameter $B^{-1}_{\eta^{\prime}}$ . A best fit with Eq. (63) and the resulting $H$ and $\sigma$ parameters are also shown.

Normalized correlation strength parameter $\lambda/\lambda_{max}$ versus average $m_T$ of the pair. The data are compared with parameter scans from Refs. [58, 59] using different values of in-medium $\eta^{\prime}$ mass $m^{*}_{\eta^{\prime}}$ and slope parameter $B^{-1}_{\eta^{\prime}}$ . A best fit with Eq. (63) and the resulting $H$ and $\sigma$ parameters are also shown.

Normalized correlation strength parameter $\lambda/\lambda_{max}$ versus average $m_T$ of the pair. The data are compared with parameter scans from Refs. [58, 59] using different values of in-medium $\eta^{\prime}$ mass $m^{*}_{\eta^{\prime}}$ and slope parameter $B^{-1}_{\eta^{\prime}}$ . A best fit with Eq. (63) and the resulting $H$ and $\sigma$ parameters are also shown.

Dependence of $\overline{p}$-$\overline{d}$ correlation on pseudorapidity acceptance of $\overline{p}$ and $\overline{d}$ selection in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV. Results are for 20.0--30.0$\%$ collision centrality.

Dependence of $\overline{p}$-$\overline{d}$ correlation on pseudorapidity acceptance of $\overline{p}$ and $\overline{d}$ selection in Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}$ = 5.02 TeV. Results are for 50.0--60.0$\%$ collision centrality.

HBT parameters of charge-combined kaon pairs, shown as value $\pm$ statistical uncertainty [absolute value] $\pm$ systematic uncertainty [%] for the centrality bins shown in Fig. 3.

$R_o/R_s$ of positive pion pairs shown as value $\pm$ statistical uncertainty [absolute value] $\pm$ systematic uncertainty [%] for the centrality bins shown in Fig. 4.

$R_o/R_s$ of negative pion pairs shown as value $\pm$ statistical uncertainty [absolute value] $\pm$ systematic uncertainty [%] for the centrality bins shown in Fig. 4.

$R_o/R_s$ of charge-combined kaon pairs shown as value $\pm$ statistical uncertainty [absolute value] $\pm$ systematic uncertainty [%] for the centrality bins shown in Fig. 4.

Azimuthal angle dependence of HBT radii of charged kaons, shown as value $\pm$ statistical uncertainty $\pm$ systematic uncertainty for the centrality bins shown in Fig. 8.

Azimuthal angle dependence of HBT radii of charged pions, shown as value $\pm$ statistical uncertainty $\pm$ systematic uncertainty for the 0%–20% and 20%–60% centrality bins.

Oscillation amplitudes relative to the event plane for charged pions, shown as value $\pm$ statistical uncertainty $\pm$ systematic uncertainty for the 0%–20% and 20%–60% centrality bins.

Oscillation amplitudes relative to the event plane for charged kaons, shown as value $\pm$ statistical uncertainty $\pm$ systematic uncertainty for the 0%–20% and 20%–60% centrality bins.

Summary of extracted parameters in the blast-wave model fit for two fitting conditions. The $T_f$ and $\rho_0$ parameters were obtained by fits to $p_T$ spectra, and the other parameters were obtained by a simultaneous fit to $v_2$ and the HBT radii.

Distribution of ${E}_{T}^{miss}$ after the fit of the multijet backgrounds, in the electron channel, in the $t\bar{t}$ 3-jets VR, without applying the cut on ${E}_{T}^{miss}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Distribution of ${E}_{T}^{miss}$ after the fit of the multijet backgrounds, in the electron channel, in the $t\bar{t}$ 4-jets VR, without applying the cut on ${E}_{T}^{miss}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Distribution of $m_{T}^{W}$ after the fit of the multijet backgrounds, in the muon channel, in the signal region, without applying the cut on $m_{T}^{W}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.