Comparison of the shapes extrapolated with a linear or quadratic extrapolation of the QCD contribution in the W+ signal region.

Comparison of the shapes extrapolated with a linear or quadratic extrapolation of the QCD contribution in the W- signal region.

Post-fit distributions of the transverse mass in the W+ signal region.

Post-fit distributions of the transverse mass in the W- signal region.

Post-fit distributions of the dimuon mass in the Z boson signal region.

Example plot for the extrapolation of the QCD contribution from a region with inverted relative isolation criterion to the W+ signal region in the transverse mass bin 12 < mT < 18 GeV. The values corresponding to the smallest relative muon isolation are extrapolated with a linear (smaller value) and quadratic (larger value) function.

Example plot for the extrapolation of the QCD contribution from a region with inverted relative isolation criterion to the W- signal region in the transverse mass bin 12 < mT < 18 GeV.The values corresponding to the smallest relative muon isolation are extrapolated with a linear (smaller value) and quadratic (larger value) function.

Example plot for the extrapolation of the QCD contribution from a region with inverted relative isolation criterion to the W+ signal region in the transverse mass bin 102 < mT < 108 GeV.The values corresponding to the smallest relative muon isolation are extrapolated with a linear (larger value) and quadratic (smaller value) function.

Example plot for the extrapolation of the QCD contribution from a region with inverted relative isolation criterion to the W- signal region in the transverse mass bin 102 < mT < 108 GeV.The values corresponding to the smallest relative muon isolation are extrapolated with a linear (smaller value) and quadratic (larger value) function.

The table compares theoretical and measured values of the product of the fiducial cross sections and branching fractions. The theoretical predictions are generated with DYTurbo in combination with three different PDF sets. All values are normalized to the corresponding measured value. The relative uncertainties on the measured values are given separately with and without the contribution from the luminosity uncertainty.

The table compares theoretical and measured values of the product of the total cross sections and branching fractions. The theoretical predictions are generated with DYTurbo in combination with three different PDF sets. All values are normalized to the corresponding measured value. The relative uncertainties on the measured values are given separately with and without the contribution from the luminosity uncertainty.

The table compares theoretical and measured ratios of the product of the fiducial cross sections and branching fractions. The theoretical predictions are generated with DYTurbo in combination with three different PDF sets. All values are normalized to the corresponding measured value.

The table compares theoretical and measured ratios of the product of the total cross sections and branching fractions. The theoretical predictions are generated with DYTurbo in combination with three different PDF sets. All values are normalized to the corresponding measured value.

Comparison of the measured and theoretical cross sections for Z, W+, W-, and W production at different center-of-mass energies. The theoretical predictions are generated with DYTurbo in combination with the NNPDF 3.1 PDF set.

Upper limits at $95\%$ CL on the LQ-fermion couplings, $|y_{e u}|$, versus $S_{e u}$ mass for scalar LQs coupled to electrons.

Observed and Expected UL exclusions on the $BR(H\to S)$ of generic signals with $m_{A'} = 1.0\;GeV$ and $BR(A' \rightarrow \pi\pi) = 1.0$.

The observed data in the dielectron channel and the fitted signal-plus-background templates, shown for the $S_{e d}$ scenario with a candidate LQ mass of 2.5 TeV. Distributions of events are binned in the reconstructed dilepton mass, rapidity, and cosine theta.

Observed and Expected UL exclusions on the number of 'signal' yields in sideband A from the model agnostic fit, assuming no 'SUEP-like' signal contamination in the CRs.

Upper limits at $95\%$ CL on the LQ-fermion couplings, $|y_{e d}|$, versus $S_{e d}$ mass for scalar LQs coupled to electrons.

Postfit values for SR in the B-Only CR+SR and CR-Only fits with the prefit signal overlayed for reference

The observed data in the dimuon channel and the fitted signal-plus-background templates, shown for the $S_{\mu u}$ scenario with a candidate LQ mass of 2.5 TeV. Distributions of events are binned in the reconstructed dilepton mass, rapidity, and cosine theta.

Postfit values for PR in the B-Only CR+SR and CR-Only fits with the prefit signal overlayed for reference

Upper limits at $95\%$ CL on the LQ-fermion couplings, $|y_{\mu u}|$, versus $S_{\mu u}$ mass for scalar LQs coupled to muons.

Postfit values for CRDY in the B-Only CR+SR and CR-Only fits with the prefit signal overlayed for reference

The observed data in the dimuon channel and the fitted signal-plus-background templates, shown for the $S_{\mu d}$ scenario with a candidate LQ mass of 2.5 TeV. Distributions of events are binned in the reconstructed dilepton mass, rapidity, and cosine theta.

Postfit values for CRTT in the B-Only CR+SR and CR-Only fits with the prefit signal overlayed for reference

Upper limits at $95\%$ CL on the LQ-fermion couplings, $|y_{\mu d}|$, versus $S_{\mu d}$ mass for scalar LQs coupled to muons.

Per-selection efficiencies on Generic - $m_{A'} = 1.0\;GeV$ $BR(A' \rightarrow \pi\pi) = 1.0$ Hadronic - $m_{A'} = 0.7\;GeV$ and $BR(A' \rightarrow ee) = BR(A' \rightarrow \mu\mu) = 0.15$ $BR(A' \rightarrow \pi\pi) = 0.7$ Leptonic - $m_{A'} = 0.5\;GeV$ and $BR(A' \rightarrow ee) = BR(A' \rightarrow \mu\mu) = 0.2$ and $BR(A' \rightarrow \pi\pi) = 0.6$ samples. Selections: 2 OSSF leptons, one SUEP cluster, $60\leq Z_{m}\leq 120$ GeV, $p_{T}^Z\geq 25$ GeV, veto events with b-tagged jets, $p_{T}^{SUEP} > 60\;GeV$, and $dR^{Ak4}_{Ak15} < 1.5$. Efficiencies are the ratio of the histogram size at the selection step to the previous step (except twoLepton vs. initial). The final column is the total cumulative efficiency (dR vs. initial). Note that these selection numbers were all calculated using samples generated with 600,000 events for each signal point considered.

The observed data in the dielectron channel and the fitted signal-plus-background templates, shown for the $V_{e u}$ scenario with a candidate LQ mass of 2.5 TeV. Distributions of events are binned in the reconstructed dilepton mass, rapidity, and cosine theta.

Interpolated observed and expected UL exclusions on $BR(H\to S) = 1$ for generic signals with $m_{A'} = 1.0\;GeV$ and $BR(A' \rightarrow \pi\pi) = 1.0$. The discreate heatmap points can be found in the UL Generic (Observed & Expected) figure in this sandbox

Upper limits at $95\%$ CL on the LQ-fermion couplings, $|g_{e u}|$, versus $V_{e u}$ mass for vector LQs coupled to electrons.

The observed data in the dielectron channel and the fitted signal-plus-background templates, shown for the $V_{e d}$ scenario with a candidate LQ mass of 2.5 TeV. Distributions of events are binned in the reconstructed dilepton mass, rapidity, and cosine theta.

Upper limits at $95\%$ CL on the LQ-fermion couplings, $|g_{e d}|$, versus $V_{e d}$ mass for vector LQs coupled to electrons.

The observed data in the dimuon channel and the fitted signal-plus-background templates, shown for the $V_{\mu u}$ scenario with a candidate LQ mass of 2.5 TeV. Distributions of events are binned in the reconstructed dilepton mass, rapidity, and cosine theta.

Upper limits at $95\%$ CL on the LQ-fermion couplings, $|g_{\mu u}|$, versus $V_{\mu u}$ mass for vector LQs coupled to muons.

The observed data in the dimuon channel and the fitted signal-plus-background templates, shown for the $V_{\mu d}$ scenario with a candidate LQ mass of 2.5 TeV. Distributions of events are binned in the reconstructed dilepton mass, rapidity, and cosine theta.

Upper limits at $95\%$ CL on the LQ-fermion couplings, $|g_{\mu d}|$, versus $V_{\mu d}$ mass for vector LQs coupled to muons.

Observed and expected 95% CL upper limits on the cross section for the resonant production of a new spin-2 particle $\mathrm{X}^{(2)}$ which decays to Higgs boson pairs, $\sigma(\mathrm{pp} \to \mathrm{X}^{(2)} \to \mathrm{HH})$, given for different values of $m_\mathrm{X}$ in the range 260-1000 GeV.

Observed and expected 95% CL upper limits on the cross section for the resonant production of a new scalar particle $\mathrm{X}$ which decays to a SM Higgs boson and a new scalar particle $\mathrm{Y}$ which subsequently decay to a pair of photons and a pair of tau leptons, $\sigma(\mathrm{pp} \to \mathrm{X} \to \mathrm{YH} \to \gamma\gamma\tau\tau)$. The limits are shown for different values of $m_\mathrm{Y}$ in the range 50-800 GeV and at particular values of $m_\mathrm{X}$ in the range 300-1000 GeV.

Observed and expected 95% CL upper limits on the cross section for the resonant production of a new scalar particle $\mathrm{X}$ which decays to a SM Higgs boson and a new scalar particle $\mathrm{Y}$ which subsequently decay to a pair of photons and a pair of tau leptons, $\sigma(\mathrm{pp} \to \mathrm{X} \to \mathrm{YH} \to \gamma\gamma\tau\tau)$. The limits are shown for different values of $m_\mathrm{X}$ in the range 300-1000 GeV and at particular values of $m_\mathrm{Y}$ in the range 50-800 GeV.

Observed and expected 95% CL upper limits on the cross section for the resonant production of a new scalar particle $\mathrm{X}$ which decays to a SM Higgs boson and a new scalar particle $\mathrm{Y}$ which subsequently decay to a pair of photons and a pair of tau leptons, $\sigma(\mathrm{pp} \to \mathrm{X} \to \mathrm{YH} \to \gamma\gamma\tau\tau)$. The limits are shown for different values of $m_\mathrm{X}$ in the range 300-1000 GeV and $m_\mathrm{Y}$ in the range 50-800 GeV.

Observed and expected 95% CL upper limits on the cross section for the resonant production of a new scalar particle $\mathrm{X}$ which decays to a SM Higgs boson and a new scalar particle $\mathrm{Y}$, multiplied by the $\mathrm{Y} \to \gamma\gamma$ branching fraction, $\sigma(\mathrm{pp} \to \mathrm{X} \to \mathrm{YH})B(\mathrm{Y} \to \gamma\gamma)$. The limits are shown for different values of $m_\mathrm{Y}$ in the range 70-125 GeV and at particular values of $m_\mathrm{X}$ in the range 300-1000 GeV.

Observed and expected 95% CL upper limits on the cross section for the resonant production of a new scalar particle $\mathrm{X}$ which decays to a SM Higgs boson and a new scalar particle $\mathrm{Y}$, multiplied by the $\mathrm{Y} \to \gamma\gamma$ branching fraction, $\sigma(\mathrm{pp} \to \mathrm{X} \to \mathrm{YH})B(\mathrm{Y} \to \gamma\gamma)$. The limits are shown for different values of $m_\mathrm{X}$ in the range 300-800 GeV and at particular values of $m_\mathrm{Y}$ in the range 70-125 GeV.

Observed and expected 95% CL upper limits on the cross section for the resonant production of a new scalar particle $\mathrm{X}$ which decays to a SM Higgs boson and a new scalar particle $\mathrm{Y}$, multiplied by the $\mathrm{Y} \to \gamma\gamma$ branching fraction, $\sigma(\mathrm{pp} \to \mathrm{X} \to \mathrm{YH})B(\mathrm{Y} \to \gamma\gamma)$. The limits are shown for different values of $m_\mathrm{X}$ in the range 300-1000 GeV and $m_\mathrm{Y}$ in the range 70-125 GeV.

Observed and expected 95% CL upper limits on the cross section for the resonant production of a new scalar particle $\mathrm{X}$ which decays to a SM Higgs boson and a new scalar particle $\mathrm{Y}$, multiplied by the $\mathrm{Y} \to \gamma\gamma$ branching fraction, $\sigma(\mathrm{pp} \to \mathrm{X} \to \mathrm{YH})B(\mathrm{Y} \to \gamma\gamma)$. The limits are shown for different values of $m_\mathrm{Y}$ in the range 125-800 GeV and at particular values of $m_\mathrm{X}$ in the range 300-1000 GeV.

Observed and expected 95% CL upper limits on the cross section for the resonant production of a new scalar particle $\mathrm{X}$ which decays to a SM Higgs boson and a new scalar particle $\mathrm{Y}$, multiplied by the $\mathrm{Y} \to \gamma\gamma$ branching fraction, $\sigma(\mathrm{pp} \to \mathrm{X} \to \mathrm{YH})B(\mathrm{Y} \to \gamma\gamma)$. The limits are shown for different values of $m_\mathrm{X}$ in the range 300-800 GeV and at particular values of $m_\mathrm{Y}$ in the range 125-800 GeV.

Observed and expected 95% CL upper limits on the cross section for the resonant production of a new scalar particle $\mathrm{X}$ which decays to a SM Higgs boson and a new scalar particle $\mathrm{Y}$, multiplied by the $\mathrm{Y} \to \gamma\gamma$ branching fraction, $\sigma(\mathrm{pp} \to \mathrm{X} \to \mathrm{YH})B(\mathrm{Y} \to \gamma\gamma)$. The limits are shown for different values of $m_\mathrm{X}$ in the range 300-1000 GeV and $m_\mathrm{Y}$ in the range 125-800 GeV.

Data vs. the ABCD method background prediction for 2017 in $\mathrm{SR_{b}^{mm}}$. Events are divided by the bin width, hence fractional data counts. Error bars show statistical uncertainties of data. The blue band shows the propagated uncertainty of all individual fit variations in a given bin, which we consider to be uncorrelated. The lower panels show the ratio of the observed data to the background estimation.

Data vs. the ABCD method background prediction for 2018 in $\mathrm{SR_{b}^{mm}}$. Events are divided by the bin width, hence fractional data counts. Error bars show statistical uncertainties of data. The blue band shows the propagated uncertainty of all individual fit variations in a given bin, which we consider to be uncorrelated. The lower panels show the ratio of the observed data to the background estimation.

Data vs. the ABCD method background prediction for 2016 in $\mathrm{SR_{b+\textrm{j}/b}^{mm}}$. Events are divided by the bin width, hence fractional data counts. Error bars show statistical uncertainties of data. The blue band shows the propagated uncertainty of all individual fit variations in a given bin, which we consider to be uncorrelated. The lower panels show the ratio of the observed data to the background estimation.

Data vs. the ABCD method background prediction for 2017 in $\mathrm{SR_{b+\textrm{j}/b}^{mm}}$. Events are divided by the bin width, hence fractional data counts. Error bars show statistical uncertainties of data. The blue band shows the propagated uncertainty of all individual fit variations in a given bin, which we consider to be uncorrelated. The lower panels show the ratio of the observed data to the background estimation.

Data vs. the ABCD method background prediction for 2018 in $\mathrm{SR_{b+\textrm{j}/b}^{mm}}$. Events are divided by the bin width, hence fractional data counts. Error bars show statistical uncertainties of data. The blue band shows the propagated uncertainty of all individual fit variations in a given bin, which we consider to be uncorrelated. The lower panels show the ratio of the observed data to the background estimation.

The asymptotic 95% CL limits on the product of cross section ($\sigma$), branching fraction into a pair of muons ($\mathcal{B}$) to dimuons, acceptance (A), and efficiency ($\epsilon$) in $\mathrm{SR_{b}^{mm}}$. A combination of both limits depends on the relative acceptance contributions in each region and is omitted here.

The asymptotic 95% CL limits on the product of cross section ($\sigma$), branching fraction into a pair of muons ($\mathcal{B}$) to dimuons, acceptance (A), and efficiency ($\epsilon$) in $\mathrm{SR_{b+\textrm{j}/b}^{mm}}$. A combination of both limits depends on the relative acceptance contributions in each region and is omitted here.

Cutflow for signal samples in 2016. Samples are divided by generated mass, signal region, and the b to s quark flavor changing coupling of the $Z^\prime$, dbs. For some mass values, multiple values of this coupling strength have been generated. It is described in this minimal lagrangian: $\mathcal{L} \supset Z^{\prime}_{\eta} \Big[ g_{\mu}\, \bar{\mu}\, \gamma^{\eta}\, \mu + \big( g_b\, \delta_{bs}\, \bar{s}\, \gamma^{\eta}\, P_{L}\, b + \text{h.c.} \big) \Big]$. The cutflow shows the selection into signal regions, and then the effects of specialized HTLT and RelMET selections, which grow with greater mass.

Cutflow for signal samples in 2017. Samples are divided by generated mass, signal region, and the b to s quark flavor changing coupling of the $Z^\prime$, dbs. For some mass values, multiple values of this coupling strength have been generated. It is described in this minimal lagrangian: $\mathcal{L} \supset Z^{\prime}_{\eta} \Big[ g_{\mu}\, \bar{\mu}\, \gamma^{\eta}\, \mu + \big( g_b\, \delta_{bs}\, \bar{s}\, \gamma^{\eta}\, P_{L}\, b + \text{h.c.} \big) \Big]$. The cutflow shows the selection into signal regions, and then the effects of specialized HTLT and RelMET selections, which grow with greater mass.

Cutflow for signal samples in 2018. Samples are divided by generated mass, signal region, and the b to s quark flavor changing coupling of the $Z^\prime$, dbs. For some mass values, multiple values of this coupling strength have been generated. It is described in this minimal lagrangian: $\mathcal{L} \supset Z^{\prime}_{\eta} \Big[ g_{\mu}\, \bar{\mu}\, \gamma^{\eta}\, \mu + \big( g_b\, \delta_{bs}\, \bar{s}\, \gamma^{\eta}\, P_{L}\, b + \text{h.c.} \big) \Big]$. The cutflow shows the selection into signal regions, and then the effects of specialized HTLT and RelMET selections, which grow with greater mass.