CL$_{s}$ value as a function of the branching fraction of $B^{+} \rightarrow K^{+} \nu \bar{\nu}$ for expected and observed signal yields and the corresponding upper limits at 90% confidence level (CL). The expected limit is derived for the background-only hypothesis. The observed limit is derived from a simultaneous fit to the on-resonance and off-resonance data, corresponding to an integrated luminosity of 63 fb$^{−1}$ and 9 fb$^{−1}$, respectively.

Signal efficiency as a function of the dineutrino invariant mass squared $q^{2}$ for events in the signal region (SR) ($\rm{BDT}_{1} > 0.9$ and $\rm{BDT}_{2} > 0.95$). The error bars indicate the statistical uncertainty.

Dijet signal efficiencies as a function of the signal mass and lifetime for events satisfying all event and vertex requirements, with corrections based on systematic differences in the vertex reconstruction efficiency between data and simulation.

The distribution of $d_{\mathrm{BV}}$ for $\geq$5-track one-vertex events in data and three simulated multijet signal samples each with a mass of 1600 GeV. The production cross section for each signal model is assumed to be the lower limit excluded by CMS-EXO-17-018, corresponding to values of 0.8, 0.25, and 0.15 fb for the samples with $c\tau =$ 0.3, 1.0, and 10 mm, respectively. The last bin includes the overflow events. This bin includes one event in data with a vertex with large $d_{\mathrm{BV}}$ that appears to arise from tracks originating from separate pp interaction vertices, consistent with background.

Distribution of the $x$-$y$ distances between vertices, $d_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution $d_{\mathrm{VV}}^{\kern 0.15em\mathrm{C}}$ (blue continuous line) is constructed from one-vertex events in data, and is normalized to the number of two-vertex events in data with two 3-track vertices. The two vertical red dashed lines separate the regions used in the fit.

Distribution of the $x$-$y$ distances between vertices, $d_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution $d_{\mathrm{VV}}^{\kern 0.15em\mathrm{C}}$ (blue continuous line) is constructed from one-vertex events in data, and is normalized to the number of two-vertex events in data which have exactly one 4-track vertex and one 3-track vertex. The two vertical red dashed lines separate the regions used in the fit.

Distribution of the $x$-$y$ distances between vertices, $d_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution $d_{\mathrm{VV}}^{\kern 0.15em\mathrm{C}}$ (blue continuous line) is constructed from one-vertex events in data, and is normalized to the number of two-vertex events in data with two 4-track vertices. The two vertical red dashed lines separate the regions used in the fit.

Distribution of the $x$-$y$ distances between vertices, $d_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution $d_{\mathrm{VV}}^{\kern 0.15em\mathrm{C}}$ (blue continuous line) is constructed from one-vertex events in data, and is normalized using $\geq$5-track one-vertex event information. The two vertical red dashed lines separate the regions used in the fit.

Predicted yields for the background-only normalized template, predicted yields for three simulated multijet signals each with a mass of 1600 GeV, and the observed yield in each $d_{\mathrm{VV}}$ bin. The production cross section for each signal model is assumed to be the lower limit excluded by CMS-EXO-17-018, corresponding to values of 0.8, 0.25, and 0.15 fb for samples with $c\tau =$ 0.3, 1.0, and 10 mm, respectively. The uncertainties in the signal yields and the systematic uncertainties in the background prediction reflect the systematic uncertainties given in the text.

Observed 95% CL upper limits on the product of cross section and branching fraction squared for the multijet signals, as a function of mass and $c\tau$. The overlaid mass-lifetime exclusion curves assume pair-production cross sections for the neutralino (red) and gluino (purple) with 100% branching fraction to each model's respective decay mode specified. The solid black (dashed colored) lines represent the observed (median expected) limits at 95% CL. The thin black lines represent the variation of the observed limit within theoretical uncertainties of the signal cross section. The thin dashed colored lines represent the region containing 68% of the expected limit distribution under the background-only hypothesis. The observed limits from the CMS displaced jets search (CMS-EXO-19-021) are also shown in teal for comparison.

Observed 95% CL upper limits on the product of cross section and branching fraction squared for the multijet signals, as a function of mass and $c\tau$. The overlaid mass-lifetime exclusion curves assume pair-production cross sections for the neutralino (red) and gluino (purple) with 100% branching fraction to each model's respective decay mode specified. The solid black (dashed colored) lines represent the observed (median expected) limits at 95% CL. The thin black lines represent the variation of the observed limit within theoretical uncertainties of the signal cross section. The thin dashed colored lines represent the region containing 68% of the expected limit distribution under the background-only hypothesis. The observed limits from the CMS displaced jets search (CMS-EXO-19-021) are also shown in teal for comparison.

Observed 95% CL upper limits on the product of cross section and branching fraction squared for the multijet signals, as a function of mass and $c\tau$. The overlaid mass-lifetime exclusion curves assume pair-production cross sections for the neutralino (red) and gluino (purple) with 100% branching fraction to each model's respective decay mode specified. The solid black (dashed colored) lines represent the observed (median expected) limits at 95% CL. The thin black lines represent the variation of the observed limit within theoretical uncertainties of the signal cross section. The thin dashed colored lines represent the region containing 68% of the expected limit distribution under the background-only hypothesis. The observed limits from the CMS displaced jets search (CMS-EXO-19-021) are also shown in teal for comparison.

Observed 95% CL upper limits on the product of cross section and branching fraction squared for the dijet signals, as a function of mass and $c\tau$. The overlaid mass-lifetime exclusion curves assume pair-production cross sections for the top squark with 100% branching fraction to each model's respective decay mode specified. The solid black (dashed colored) lines represent the observed (median expected) limits at 95% CL. The thin black lines represent the variation of the observed limit within theoretical uncertainties of the signal cross section. The thin dashed colored lines represent the region containing 68% of the expected limit distribution under the background-only hypothesis. The observed limits from the CMS displaced jets search (CMS-EXO-19-021) are also shown in teal for comparison.

Observed 95% CL upper limits on the product of cross section and branching fraction squared for the dijet signals, as a function of mass and $c\tau$. The overlaid mass-lifetime exclusion curves assume pair-production cross sections for the top squark with 100% branching fraction to each model's respective decay mode specified. The solid black (dashed colored) lines represent the observed (median expected) limits at 95% CL. The thin black lines represent the variation of the observed limit within theoretical uncertainties of the signal cross section. The thin dashed colored lines represent the region containing 68% of the expected limit distribution under the background-only hypothesis. The observed limits from the CMS displaced jets search (CMS-EXO-19-021) are also shown in teal for comparison.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of mass for multijet signals, for a fixed $c\tau$ of 300um in the full Run-2 data set. The neutralino and gluino pair production cross sections are shown for the multijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of mass for dijet signals, for a fixed $c\tau$ of 300um in the full Run-2 data set. The top squark pair-production cross section is shown for the dijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of mass for multijet signals, for a fixed $c\tau$ of 1 mm in the full Run-2 data set. The neutralino and gluino pair production cross sections are shown for the multijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of mass for dijet signals, for a fixed $c\tau$ of 1 mm in the full Run-2 data set. The top squark pair-production cross section is shown for the dijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of mass for multijet signals, for a fixed $c\tau$ of 10 mm in the full Run-2 data set. The neutralino and gluino pair production cross sections are shown for the multijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of mass for dijet signals, for a fixed $c\tau$ of 10 mm in the full Run-2 data set. The top squark pair-production cross section is shown for the dijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of $c\tau$ for multijet signals, for a fixed mass of 800 GeV in the full Run-2 data set. The neutralino and gluino pair production cross sections are shown for the multijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of $c\tau$ for dijet signals, for a fixed mass of 800 GeV in the full Run-2 data set. The top squark pair-production cross section is shown for the dijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of $c\tau$ for multijet signals, for a fixed mass of 1600 GeV in the full Run-2 data set. The neutralino and gluino pair production cross sections are shown for the multijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of $c\tau$ for dijet signals, for a fixed mass of 1600 GeV in the full Run-2 data set. The top squark pair-production cross section is shown for the dijet signals.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of $c\tau$ for multijet signals, for a fixed mass of 2400 GeV in the full Run-2 data set. The neutralino and gluino pair production cross sections are shown for the multijet signals. For $m$ = 2400 GeV, the expected neutralino cross section is $\approx 8\times 10^{-5}$ fb and is not shown.

Observed and expected 95% CL upper limits on the product of cross section and branching fraction squared, as a function of $c\tau$ for dijet signals, for a fixed mass of 2400 GeV in the full Run-2 data set. The top squark pair-production cross section is shown for the dijet signals.

Data-to-simulation efficiency correction factors, shown for multijet and dijet signal topologies in several ranges of $c\tau$. Note that these correction factors account for the two long-lived particles in the simulated events, and are therefore the total correction factors used to scale event yields rather than the correction factors one would apply to individual vertices.

Distribution of the azimuthal angle between vertices, $\Delta\phi_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution (blue continuous line) is constructed from 3-track one-vertex events in data, and is normalized to the number of 3-track two-vertex events in data.

Distribution of the azimuthal angle between vertices, $\Delta\phi_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution (blue continuous line) is constructed from 4-track and 3-track one-vertex events in data, and is normalized to the number of two-vertex events in data which have exactly one 4-track vertex and one 3-track vertex.

Distribution of the azimuthal angle between vertices, $\Delta\phi_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution (blue continuous line) is constructed from 4-track one-vertex events in data, and is normalized to the number of 4-track two-vertex events in data.

Distribution of the azimuthal angle between vertices, $\Delta\phi_{\mathrm{VV}}$, for 2017 and 2018 data. The background distribution (blue continuous line) is constructed from $\geq$5-track one-vertex events in data, and is normalized using one-vertex event information. No $\geq$5-track two-vertex data events pass the selection.

95% CL observed upper limits on the Yukawa couplings

The number of background events obtained from the 1- and 2-leg predictions derived from data, together with the direct observation in data, in bins in ChF, where either the leading or subleading jet has a ChF within the bin edges, and both have a ChF below the upper bin threshold. The bottom panel shows the ratios of the observation in data to the 1-leg and the 2-leg background predictions.

The expected and observed 95% CL upper limits on the production cross section for SIMPs with masses between 1 and 1000 GeV, with the assumption that the SIMP interaction in the detector can be approximated as neutron-like. The theoretical prediction of a simplified model incorporating this approximation and including a scalar mediator with couplings g_chi = -1 and g_q = 1 is also shown (red line). For masses above 100 GeV, where the modelling of the SIMP-nucleon interaction becomes more speculative, the obtained cross section upper limits are increasingly uncertain (shaded area).

The best-fit predicted and observed distributions of the combined NN output in the high-$p_\text{T}$-SR for the $\mu\tau_e$ channel. The first and last bin include underflow and overflow events, respectively.

The best-fit predicted and observed distributions of the collinear mass $m_\text{coll}$ in the high-$p_\text{T}$-SR for the $e\tau_\mu$ channel. The first and last bin include underflow and overflow events, respectively.

The best-fit predicted and observed distributions of the collinear mass $m_\text{coll}$ in the high-$p_\text{T}$-SR for the $\mu\tau_e$ channel. The first and last bin include underflow and overflow events, respectively.

Observed and expected upper limits on $\mathcal{B}(Z\rightarrow\ell\tau)$ with leptonically decaying $\tau$ at 95% confidence level.

Observed and expected upper limits on $\mathcal{B}(Z\rightarrow\ell\tau)$ at 95% confidence level, combination of hadronically and leptonically decaying $\tau$ channels.

The best-fit predicted and observed distributions of the collinear mass $m_\text{coll}$ in the CRZ$\tau\tau$ for the $e\tau_\mu$ channel. The first and last bin include underflow and overflow events, respectively.

The best-fit predicted and observed distributions of the collinear mass $m_\text{coll}$ in the CRZ$\tau\tau$ for the $\mu\tau_e$ channel. The first and last bin include underflow and overflow events, respectively.

The best-fit predicted and observed distributions of the combined NN output in the CRZ$\tau\tau$ for the $e\tau_\mu$ channel. The first and last bin include underflow and overflow events, respectively.

The best-fit predicted and observed distributions of the combined NN output in the CRZ$\tau\tau$ for the $\mu\tau_e$ channel. The first and last bin include underflow and overflow events, respectively.

The best-fit predicted and observed distributions of the leading lepton transverse momentum $p_\text{T}(e)$ in the high-$p_\text{T}$-SR for the $e\tau_\mu$ channel. The first and last bin include underflow and overflow events, respectively.

The best-fit predicted and observed distributions of the leading lepton transverse momentum $p_\text{T}(\mu)$ in the high-$p_\text{T}$-SR for the $\mu\tau_e$ channel. The first and last bin include underflow and overflow events, respectively.

The best-fit predicted and observed distributions of the missing transverse momentum $E^\text{miss}_\text{T}$ in the low-$p_\text{T}$-SR for the $e\tau_\mu$ channel. The first and last bin include underflow and overflow events, respectively.

The best-fit predicted and observed distributions of the missing transverse momentum $E^\text{miss}_\text{T}$ in the low-$p_\text{T}$-SR for the $\mu\tau_e$ channel. The first and last bin include underflow and overflow events, respectively.

Distribution of the transverse momentum of the diphoton system for the $\mathrm{Z}\gamma\gamma$ muon channel. The predicted yields are shown with their pre-fit normalisations. The observed data, the expected signal contribution and the background estimates are presented with error bars showing the corresponding statistical uncertainties.

Distribution of the transverse momentum of the diphoton system, obtained in the control region enriched in misidentified photons for the $\mathrm{W}\gamma\gamma$ electron channel. The predicted yields are shown with their pre-fit normalisations. The observed data, the expected signal contribution and the background estimates are presented with error bars showing the corresponding statistical uncertainties.

Distribution of the transverse momentum of the diphoton system, obtained in the control region enriched in misidentified photons for the $\mathrm{Z}\gamma\gamma$ electron channel. The predicted yields are shown with their pre-fit normalisations. The observed data, the expected signal contribution and the background estimates are presented with error bars showing the corresponding statistical uncertainties.

The measured $\mathrm{W}\gamma\gamma$ and $\mathrm{Z}\gamma\gamma$ fiducial cross sections and corresponding theoretical predictions from MadGraph5_aMC@NLO. The measured yields in the electron and muon channels are extrapolated to a common fiducial phase space determined from simulated signal events at the generated particle level. Generated particles are considered stable if their mean decay length is larger than $1\,\mathrm{cm}$. Generated leptons are required to have a $\mathrm{p}_\mathrm{t}>15\,\mathrm{GeV}$ and $|\eta|<2.5$. The momenta of photons in a cone of $\Delta R=0.1$, the same cone size as the one applied to reconstructed data, are added to the charged lepton momentum to correct for final-state radiation. Generated photons are required to have $\mathrm{p}_\mathrm{t}>15\,\mathrm{GeV}$ and $|\eta|<2.5$. Additionally, the candidate photons are required to have no selected leptons or photons in a cone of radius $\Delta R=0.4$ and no other stable particles, apart from photons and neutrinos, in a cone of radius $\Delta R=0.1$. Events are then selected in the $\mathrm{W}\gamma\gamma$ channel by requiring exactly one electron (muon) with $\mathrm{p}_\mathrm{t}>30\,\mathrm{GeV}$ and at least two photons with $\mathrm{p}_\mathrm{t}>20\,\mathrm{GeV}$. Events are selected in the $\mathrm{Z}\gamma\gamma$ channel by requiring two electrons (muons), at least one of them with $\mathrm{p}_\mathrm{t}>30\,\mathrm{GeV}$, and not less than two photons, each of them with $\mathrm{p}_\mathrm{t}>20\,\mathrm{GeV}$. Additionally, the invariant mass of the dilepton system is required to be $m_{\ell\ell}>55\,\mathrm{GeV}$.

Expected and observed 95% confidence level intervals for the different anomalous couplings parameters in both the $\ell\nu\gamma\gamma$ and $\ell\ell\gamma\gamma$ channels. All parameters are fixed to their SM values except the one that is fitted. No unitarity regularisation scheme is applied. Intervals with a NAN value are not measured in the corresponding channels.

Differential cross sections as function of transverse momentum imbalance between Z boson and dijet for Z+ ≥ 2 jets events.

Differential cross sections as function of transverse momentum imbalance between leading and subleading jet for Z+ ≥ 2 jets events.

Area normalized distribution as function of Delta Phi between Z boson and the leading jet for Z+ ≥ 1 jet events.

Area normalized distribution as function of transverse momentum imbalance between Z boson and the leading jet for Z+ ≥ 1 jet events.

Area normalized distribution as function of Delta Phi between Z boson and dijet for Z+ ≥ 2 jets events.

Area normalized distribution as function of transverse momentum imbalance between Z boson and dijet for Z+ ≥ 2 jets events.

Area normalized distribution as function of transverse momentum imbalance between leading and subleading jet for Z+ ≥ 2 jets events

Response matrix for Delta Phi between Z boson and the leading jet for Z+ ≥ 1 jet events.

Response matrix for transverse momentum imbalance between Z boson and the leading jet for Z+ ≥ 1 jet events.

Response matrix for Delta Phi between Z boson and dijet for Z+ ≥ 2 jets events.

Response matrix for transverse momentum imbalance between Z boson and dijet for Z+ ≥ 2 jets events.

Response matrix for transverse momentum imbalance between the leading jet and subleading jet for Z+ ≥ 2 jets events.

Correlation matrix for Delta Phi between Z boson and the leading jet for Z+ ≥ 1 jet events (for absolute differential cross section measurements).

Correlation matrix for transverse momentum imbalance between Z boson and the leading jet for Z+ ≥ 1 jet events (for absolute differential cross section measurements).

Correlation matrix for Delta Phi between Z boson and dijet for Z+ ≥ 2 jets events (for absolute differential cross section measurements).

Correlation matrix for transverse momentum imbalance between Z boson and dijet for Z+ ≥ 2 jets events (for absolute differential cross section measurements).

Correlation matrix for transverse momentum imbalance between the leading jet and subleading jet for Z+ ≥ 2 jets events (for absolute differential cross section measurements).

Correlation matrix for Delta Phi between Z boson and the leading jet for Z+ ≥ 1 jet events (for normalized differential cross section measurements).

Correlation matrix for transverse momentum imbalance between Z boson and the leading jet for Z+ ≥ 1 jet events (for normalized differential cross section measurements).

Correlation matrix for Delta Phi between Z boson and dijet for Z+ ≥ 2 jets events (for normalized differential cross section measurements).

Correlation matrix for transverse momentum imbalance between Z boson and dijet for Z+ ≥ 2 jets events (for normalized differential cross section measurements).

Correlation matrix for transverse momentum imbalance between the leading jet and subleading jet for Z+ ≥ 2 jets events (for normalized differential cross section measurements).

Expected and observed exclusion limits for heavy higgs boson mass 360 GeV

Expected and observed exclusion limits for heavy higgs boson mass 400 GeV

Expected and observed exclusion limits for heavy higgs boson mass 450 GeV

Expected and observed exclusion limits for heavy higgs boson mass 500 GeV

Expected and observed exclusion limits for heavy higgs boson mass 550 GeV

Expected and observed exclusion limits for heavy higgs boson mass 600 GeV

Expected and observed exclusion limits for heavy higgs boson mass 700 GeV

Expected and observed exclusion limits for heavy higgs boson mass 800 GeV

Expected and observed exclusion limits for heavy higgs boson mass 900 GeV

Expected and observed exclusion limits for heavy higgs boson mass 1000 GeV

Expected and observed exclusion limits for heavy higgs boson mass 1200 GeV

Expected and observed exclusion limits for heavy higgs boson mass 1400 GeV

Expected and observed exclusion limits for heavy higgs boson mass 1600 GeV

Expected and observed exclusion limits for heavy higgs boson mass 1800 GeV

Expected and observed exclusion limits for heavy higgs boson mass 2000 GeV

Expected and observed exclusion limits for heavy higgs boson mass 2500 GeV

Expected and observed exclusion limits for heavy higgs boson mass 3000 GeV

Observed and expected (background-only fitted) invariant mass spectra of Wgamma events. The fitted function is ${ d N}/{ d m} = p_{0} * (m/\sqrt{s})^{p_{1} + p_{2} * \log(m/\sqrt{s}) + p_{3} * \log^{2}(m/\sqrt{s})}$

Expected and observed 95% CL upper limits on the product of the cross section and branching fraction for narrow scalar Wgamma resonances. Limits are compared to predicted cross sections for the heavy scalar triplet model described in arXiv:1912.08234

Expected and observed 95% CL upper limits on the product of the cross section and branching fraction for broad scalar Wgamma resonances.

Expected and observed 95% CL upper limits on the product of the cross section and branching fraction for narrow vector Wgamma resonances. Limits are compared to predicted cross sections for the heavy vector triplet model described in arXiv:1912.08234

Expected and observed 95% CL upper limits on the product of the cross section and branching fraction for broad vector Wgamma resonances.

Expected and observed model-independent 95% CL upper limits on the product of the cross section, branching fraction and signal acceptance for general Wgamma resonances.

Expected and observed model-independent 95% CL upper limits on the product of the cross section, branching fraction, signal acceptance and W tagging efficiency for general Jgamma resonances.

The measured single-differential cross sections in leading jet transverse momenta for the pure electroweak Z$\gamma$jj production. The total uncertainty includes the stastical uncertianty and the systematic uncertainty, while the uncertainty of the predicted results is the theoretical uncertainty from the MadGraph5_aMC@NLO. The last bin includes overflow events.

The measured single-differential cross sections in leading lepton transverse momenta for the pure electroweak Z$\gamma$jj production. The total uncertainty includes the stastical uncertianty and the systematic uncertainty, while the uncertainty of the predicted results is the theoretical uncertainty from the MadGraph5_aMC@NLO. The last bin includes overflow events.

The measured double-differential cross sections as functions of $m_{\mathrm{jj}}$ and $|\Delta\eta_{\mathrm{jj}}|$ for the pure electroweak Z$\gamma$jj production. The total uncertainty includes the stastical uncertianty and the systematic uncertainty, while the uncertainty of the predicted results is the theoretical uncertainty from the MadGraph5_aMC@NLO. The last bin includes overflow events.

The measured single-differential cross sections in photon transverse momenta for the combined QCD-induced and electroweak Z$\gamma$jj production. The total uncertainty includes the stastical uncertianty and the systematic uncertainty, while the uncertainty of the predicted results is the theoretical uncertainty from the MadGraph5_aMC@NLO. The last bin includes overflow events.

The measured single-differential cross sections in leading jet transverse momenta for the combined QCD-induced and electroweak Z$\gamma$jj production. The total uncertainty includes the stastical uncertianty and the systematic uncertainty, while the uncertainty of the predicted results is the theoretical uncertainty from the MadGraph5_aMC@NLO. The last bin includes overflow events.

The measured single-differential cross sections in leading lepton transverse momenta for the combined QCD-induced and electroweak Z$\gamma$jj production. The total uncertainty includes the stastical uncertianty and the systematic uncertainty, while the uncertainty of the predicted results is the theoretical uncertainty from the MadGraph5_aMC@NLO. The last bin includes overflow events.

The measured double-differential cross sections as functions of $m_{\mathrm{jj}}$ and $|\Delta\eta_{\mathrm{jj}}|$ for the combined QCD-induced and electroweak Z$\gamma$jj production. The total uncertainty includes the stastical uncertianty and the systematic uncertainty, while the uncertainty of the predicted results is the theoretical uncertainty from the MadGraph5_aMC@NLO. The last bin includes overflow events.

The expected and observed limits on the aQGC parameters at 95% confidence level. The last column presents the scattering energy values for which the amplitude would violate unitarity for the observed value of the aQGC parameter. All coupling parameter limits are set in TeV$^{-4}$, whereas the unitarity bounds are in TeV.