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

Observed (solid line) value of $-2\ln\Lambda$ as a function of $\kappa_{2V}$ for the combined fit, when all other coupling modifiers are fixed to their SM predictions.

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

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

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

Likelihood contours at 68% (solid line) and 95% (dashed line) CL in the ($c_{hhh}$, $c_{gghh}$) parameter space, when all other Wilson coefficients are fixed to their SM values. The corresponding expected contours are shown by the inner and outer shaded regions. The SM prediction is indicated by the star, while the observed best-fit value is denoted by the black cross. The best-fit value corresponds to the very first entry in the table.

Likelihood contours at 68% (solid line) and 95% (dashed line) CL in the ($c_{hhh}$, $c_{gghh}$) parameter space, when all other Wilson coefficients are fixed to their SM values. The corresponding expected contours are shown by the inner and outer shaded regions. The SM prediction is indicated by the star, while the observed best-fit value is denoted by the black cross.

Likelihood contours at 68% (solid line) and 95% (dashed line) CL in the ($c_{hhh}$, $c_{tthh}$) parameter space, when all other Wilson coefficients are fixed to their SM values. The corresponding expected contours are shown by the inner and outer shaded regions. The SM prediction is indicated by the star, while the observed best-fit value is denoted by the black cross. The best-fit value corresponds to the very first entry in the table.

Likelihood contours at 68% (solid line) and 95% (dashed line) CL in the ($c_{hhh}$, $c_{tthh}$) parameter space, when all other Wilson coefficients are fixed to their SM values. The corresponding expected contours are shown by the inner and outer shaded regions. The SM prediction is indicated by the star, while the observed best-fit value is denoted by the black cross.

Likelihood contours at 68% (solid line) and 95% (dashed line) CL in the ($c_{H}$, $c_{H◻}$) parameter space, when all other Wilson coefficients are fixed to their SM values. The corresponding expected contours are shown by the inner and outer shaded regions. The SM prediction is indicated by the star, while the observed best-fit value is denoted by the black cross. The best-fit value corresponds to the very first entry in the table.

Likelihood contours at 68% (solid line) and 95% (dashed line) CL in the ($c_{H}$, $c_{H◻}$) parameter space, when all other Wilson coefficients are fixed to their SM values. The corresponding expected contours are shown by the inner and outer shaded regions. The SM prediction is indicated by the star, while the observed best-fit value is denoted by the black cross.

The acceptance times efficiency [in %] for the signal ggF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category for the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ channel. The other coupling modifiers affecting ggF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal ggF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ SLT channel. The other coupling modifiers affecting ggF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal ggF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ LTT channel. The other coupling modifiers affecting ggF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category for the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ channel. The other coupling modifiers affecting VBF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ SLT channel. The other coupling modifiers affecting VBF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ LTT channel. The other coupling modifiers affecting VBF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{2V}$ in each analysis category for the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ channel. The other coupling modifiers affecting VBF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{2V}$ in each analysis category for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ SLT channel. The other coupling modifiers affecting VBF $HH$ production are fixed to their SM predictions.

The acceptance times efficiency [in %] for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{2V}$ in each analysis category for the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ LTT channel. The other coupling modifiers affecting VBF $HH$ production are fixed to their SM predictions.

Expected (dashed line) value of $-2\ln\Lambda$ as a function of $\kappa_{\lambda}$ for the fit to the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

Expected (dashed line) value of $-2\ln\Lambda$ as a function of $\kappa_{\lambda}$ for the fit to the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

Observed (solid line) value of $-2\ln\Lambda$ as a function of $\kappa_{\lambda}$ for the fit to the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

Observed (solid line) value of $-2\ln\Lambda$ as a function of $\kappa_{\lambda}$ for the fit to the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

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

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

Observed (solid line) value of $-2\ln\Lambda$ as a function of $\kappa_{2V}$ for the fit to the $\tau_{\mathrm{had}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

Observed (solid line) value of $-2\ln\Lambda$ as a function of $\kappa_{2V}$ for the fit to the $\tau_{\mathrm{lep}}\tau_{\mathrm{had}}$ channel, when all other coupling modifiers are fixed to their SM predictions.

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

Expected CLs values in the High mA SR, mA vs ma signal grid.

Observed CLs values in the High mA SR, mA vs ma signal grid.

Expected CLs values in the Low mA SR, mA vs ma signal grid.

Observed CLs values in the High mA SR, mA vs ma signal grid.

CLs+b values in the Low mA SR, mA vs tanB signal grid.

CLs+b values in the High mA SR, mA vs ma signal grid.

CLs+b values in the Low mA SR, mA vs ma signal grid.

Cut flow of the 2HDM+a signal points, gluon–gluon fusion production, Low mA SR. tanB = 1, $sin\theta$ = 0.35. The first two entries in the tables are number of raw MC events, third entry is theoretical prediction, and all other lines include the correct weights. Note that during the generation Higgs boson’s branching ratio to taus has been set to 1. An additional factor of 0.0627 is used to account for that, starting from the ‘Initial’ entry.

Cut flow of the 2HDM+a signal points, gluon–gluon fusion production, High mA SR. tanB = 1, $sin\theta$ = 0.35. The first two entries in the tables are number of raw MC events, third entry is theoretical prediction, and all other lines include the correct weights. Note that during the generation Higgs boson’s branching ratio to taus has been set to 1. An additional factor of 0.0627 is used to account for that, starting from the ‘Initial’ entry.

Cut flow of the 2HDM+a signal points, bb annihilation production, Low mA SR. tanB = 1, $sin\theta$ = 0.35. The first two entries in the tables are number of raw MC events, third entry is theoretical prediction, and all other lines include the correct weights. Note that during the generation Higgs boson’s branching ratio to taus has been set to 1. An additional factor of 0.0627 is used to account for that, starting from the ‘Initial’ entry.

Cut flow of the 2HDM+a signal points, bb annihilation production, High mA SR. tanB = 1, $sin\theta$ = 0.35. The first two entries in the tables are number of raw MC events, third entry is theoretical prediction, and all other lines include the correct weights. Note that during the generation Higgs boson’s branching ratio to taus has been set to 1. An additional factor of 0.0627 is used to account for that, starting from the ‘Initial’ entry.

A comparison of the observed and expected yields in the four bins of the Low mA SR.

A comparison of the observed and expected yields in the two bins of the High mA SR.

Expected exclusion contours at 95% CL as a function of mA and ma.

Observed exclusion contours at 95% CL as a function of mA and ma.

Expected +- 1sigma exclusion contours at 95% CL as a function of mA and ma.

Expected exclusion contours at 95% CL as a function of mA and ma.

Observed exclusion contours at 95% CL as a function of mA and ma.

Expected +- 1sigma exclusion contours at 95% CL as a function of mA and ma.

Acceptance, High mA SR, mA vs tanB grid, 400-750 GeV, bb prod

Acceptance, High mA SR, mA vs tanB grid, >750 GeV, bb prod

Acceptance, Low mA SR, mA vs tanB grid, 100-250 GeV, bb prod

Acceptance, Low mA SR, mA vs tanB grid, 250-400 GeV, bb prod

Acceptance, Low mA SR, mA vs tanB grid, 400-550 GeV, bb prod

Acceptance, Low mA SR, mA vs tanB grid, >550 GeV, bb prod

Acceptance, High mA SR, mA vs ma grid, 400-750 GeV, bb prod

Acceptance, High mA SR, mA vs ma grid, >750 GeV, bb prod

Acceptance, Low mA SR, mA vs ma grid, 100-250 GeV, bb prod

Acceptance, Low mA SR, mA vs ma grid, 250-400 GeV, bb prod

Acceptance, Low mA SR, mA vs ma grid, 400-550 GeV, bb prod

Acceptance, Low mA SR, mA vs ma grid, >550 GeV, bb prod

Acceptance, High mA SR, mA vs tanB grid, 400-750 GeV, ggF prod

Acceptance, High mA SR, mA vs tanB grid, >750 GeV, ggF prod

Acceptance, Low mA SR, mA vs tanB grid, 100-250 GeV, ggF prod

Acceptance, Low mA SR, mA vs tanB grid, 250-400 GeV, ggF prod

Acceptance, Low mA SR, mA vs tanB grid, 400-550 GeV, ggF prod

Acceptance, Low mA SR, mA vs tanB grid, >550 GeV, ggF prod

Acceptance, High mA SR, mA vs ma grid, 400-750 GeV, ggF prod

Acceptance, High mA SR, mA vs ma grid, >750 GeV, ggF prod

Acceptance, Low mA SR, mA vs ma grid, 100-250 GeV, ggF prod

Acceptance, Low mA SR, mA vs ma grid, 250-400 GeV, ggF prod

Acceptance, Low mA SR, mA vs ma grid, 400-550 GeV, ggF prod

Acceptance, Low mA SR, mA vs ma grid, >550 GeV, ggF prod

Efficiency, High mA SR, mA vs tanB grid, 400-750 GeV, bb prod

Efficiency, High mA SR, mA vs tanB grid, >750 GeV, bb prod

Efficiency, Low mA SR, mA vs tanB grid, 100-250 GeV, bb prod

Efficiency, Low mA SR, mA vs tanB grid, 250-400 GeV, bb prod

Efficiency, Low mA SR, mA vs tanB grid, 400-550 GeV, bb prod

Efficiency, Low mA SR, mA vs tanB grid, >550 GeV, bb prod

Efficiency, High mA SR, mA vs ma grid, 400-750 GeV, bb prod

Efficiency, High mA SR, mA vs ma grid, >750 GeV, bb prod

Efficiency, Low mA SR, mA vs ma grid, 100-250 GeV, bb prod

Efficiency, Low mA SR, mA vs ma grid, 250-400 GeV, bb prod

Efficiency, Low mA SR, mA vs ma grid, 400-550 GeV, bb prod

Efficiency, Low mA SR, mA vs ma grid, >550 GeV, bb prod

Efficiency, High mA SR, mA vs tanB grid, 400-750 GeV, ggF prod

Efficiency, High mA SR, mA vs tanB grid, >750 GeV, ggF prod

Efficiency, Low mA SR, mA vs tanB grid, 100-250 GeV, ggF prod

Efficiency, Low mA SR, mA vs tanB grid, 250-400 GeV, ggF prod

Efficiency, Low mA SR, mA vs tanB grid, 400-550 GeV, ggF prod

Efficiency, Low mA SR, mA vs tanB grid, >550 GeV, ggF prod

Efficiency, High mA SR, mA vs ma grid, 400-750 GeV, ggF prod

Efficiency, High mA SR, mA vs ma grid, >750 GeV, ggF prod

Efficiency, Low mA SR, mA vs ma grid, 100-250 GeV, ggF prod

Efficiency, Low mA SR, mA vs ma grid, 250-400 GeV, ggF prod

Efficiency, Low mA SR, mA vs ma grid, 400-550 GeV, ggF prod

Efficiency, Low mA SR, mA vs ma grid, >550 GeV, ggF prod

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Cutflow table for a pair produced coloron of 1500 GeV decaying into two quarks.

The observed number of data, background and top squark signal events in each of the signal regions of the inclusive selection

The observed number of data, background and coloron signal events in each of the signal regions of the inclusive selection

The observed number of data, background and top squark signal events in each of the signal regions of the b-tagged selection

Signal acceptance and efficiency (in %) as a function of M(STOP), before mass windows

Signal acceptance (in %) and efficiency as a function of M(STOP), after mass windows

Signal acceptance and efficiency (in %) as a function of M(RHO), before mass windows

Signal acceptance and efficiency (in %) as a function of M(RHO), after mass windows

Signal acceptance (in %) and efficiency as a function of M(STOP), before mass windows

Signal acceptance (in %) and efficiency as a function of M(STOP), after mass windows

Cross section excluded at 95% CL as a function of the top squark mass, for a pair produced top squark with decays into a pair of light-quarks.

Cross section excluded at 95% CL as a function of the cooron mass, for a pair produced coloron with decays into a pair of light-quarks.

Cross section excluded at 95% CL as a function of the top squark mass, for a pair produced top squark with decays into a b- and an s-quark.

Observed exclusion limit contour in the (BRe,BRtau) plane for a 700 GeV stop. All limits are computed at 95% CL.

Expected exclusion limit contour in the (BRe,BRtau) plane for a 800 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit contour in the (BRe,BRtau) plane for a 800 GeV stop. All limits are computed at 95% CL.

Expected exclusion limit contour in the (BRe,BRtau) plane for a 900 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit contour in the (BRe,BRtau) plane for a 900 GeV stop. All limits are computed at 95% CL.

Expected exclusion limit contour in the (BRe,BRtau) plane for a 1000 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit contour in the (BRe,BRtau) plane for a 1000 GeV stop. All limits are computed at 95% CL.

Expected exclusion limit contour in the (BRe,BRtau) plane for a 1050 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit contour in the (BRe,BRtau) plane for a 1050 GeV stop. All limits are computed at 95% CL.

Expected exclusion limit contour in the (BRe,BRtau) plane for a 1100 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit contour in the (BRe,BRtau) plane for a 1100 GeV stop. All limits are computed at 95% CL.

Expected exclusion limit contour in the (BRe,BRtau) plane for a 1150 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit contour in the (BRe,BRtau) plane for a 1150 GeV stop. All limits are computed at 95% CL.

Expected exclusion limit contour in the (BRe,BRtau) plane for a 1200 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit contour in the (BRe,BRtau) plane for a 1200 GeV stop. All limits are computed at 95% CL.

Expected exclusion limit contour in the (BRe,BRtau) plane for a 1250 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit contour in the (BRe,BRtau) plane for a 1250 GeV stop. All limits are computed at 95% CL.

Expected exclusion limit contour in the (BRe,BRtau) plane for a 1300 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit contour in the (BRe,BRtau) plane for a 1300 GeV stop. All limits are computed at 95% CL.

Expected exclusion limit contour in the (BRe,BRtau) plane for a 1350 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit contour in the (BRe,BRtau) plane for a 1350 GeV stop. All limits are computed at 95% CL.

Expected exclusion limit contour in the (BRe,BRtau) plane for a 1400 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit contour in the (BRe,BRtau) plane for a 1400 GeV stop. All limits are computed at 95% CL.

Expected exclusion limit contour in the (BRe,BRtau) plane for a 1450 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit contour in the (BRe,BRtau) plane for a 1450 GeV stop. All limits are computed at 95% CL.

Expected exclusion limit contour in the (BRe,BRtau) plane for a 1500 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit contour in the (BRe,BRtau) plane for a 1500 GeV stop. All limits are computed at 95% CL.

Signal acceptance (in %) in the (BRe,BRtau) plane for a 800 GeV stop, for the SR800 signal region.

Signal acceptance (in %) in the (BRe,BRtau) plane for a 800 GeV stop, for the SR1100 signal region.

Signal efficiency (in %) in the (BRe,BRtau) plane for a 800 GeV stop, for the SR800 signal region.

Signal efficiency (in %) in the (BRe,BRtau) plane for a 800 GeV stop, for the SR1100 signal region.

Signal acceptance (in %) in the (BRe,BRtau) plane for a 1200 GeV stop, for the SR800 signal region.

Signal acceptance (in %) in the (BRe,BRtau) plane for a 1200 GeV stop, for the SR1100 signal region.

Signal efficiency (in %) in the (BRe,BRtau) plane for a 1200 GeV stop, for the SR800 signal region.

Signal efficiency (in %) in the (BRe,BRtau) plane for a 1200 GeV stop, for the SR1100 signal region.

Signal acceptance (in %) in the (BRe,BRtau) plane for a 1500 GeV stop, for the SR800 signal region.

Signal acceptance (in %) in the (BRe,BRtau) plane for a 1500 GeV stop, for the SR1100 signal region.

Signal efficiency (in %) in the (BRe,BRtau) plane for a 1500 GeV stop, for the SR800 signal region.

Signal efficiency (in %) in the (BRe,BRtau) plane for a 1500 GeV stop, for the SR1100 signal region.

$m_{bl}^{0}$ distribution in SR800. All selection criteria are applied, except the selection on the variable that is displayed in the plot. The SM backgrounds are normalized to the values determined in the background-only fit. The last bin includes overflows.

$m_{bl}^\mathrm{asym}$ distribution in SR800. All selection criteria are applied, except the selection on the variable that is displayed in the plot. The SM backgrounds are normalized to the values determined in the background-only fit. The last bin includes overflows.

$H_\mathrm{T}$ distribution in SR800. All selection criteria are applied, except the selection on the variable that is displayed in the plot. The SM backgrounds are normalized to the values determined in the background-only fit. The last bin includes overflows.

$m_{ll}$ distribution in SR800. All selection criteria are applied, except the selection on the variable that is displayed in the plot. The SM backgrounds are normalized to the values determined in the background-only fit. The last bin includes overflows.

$m_{bl}^{1}$(rej) distribution in SR800. All selection criteria are applied, except the selection on the variable that is displayed in the plot. The SM backgrounds are normalized to the values determined in the background-only fit. The last bin includes overflows.

Full list of event selections and MC generator-weighted yields and efficiencies in the inclusive SR800 and SR1100 signal regions for several signal samples of varying stop mass with decay into b-electron, b-muon or b-tau at 1/3 branching ratio.

Observed exclusion limit in the (BRe,BRtau) plane on the cross section for a 600 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit in the (BRe,BRtau) plane on the cross section for a 700 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit in the (BRe,BRtau) plane on the cross section for a 800 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit in the (BRe,BRtau) plane on the cross section for a 900 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit in the (BRe,BRtau) plane on the cross section for a 1000 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit in the (BRe,BRtau) plane on the cross section for a 1050 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit in the (BRe,BRtau) plane on the cross section for a 1100 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit in the (BRe,BRtau) plane on the cross section for a 1150 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit in the (BRe,BRtau) plane on the cross section for a 1200 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit in the (BRe,BRtau) plane on the cross section for a 1250 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit in the (BRe,BRtau) plane on the cross section for a 1300 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit in the (BRe,BRtau) plane on the cross section for a 1350 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit in the (BRe,BRtau) plane on the cross section for a 1400 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit in the (BRe,BRtau) plane on the cross section for a 1450 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit in the (BRe,BRtau) plane on the cross section for a 1500 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit in the (BRe,BRtau) plane on the cross section for a 1550 GeV stop. All limits are computed at 95% CL.

Observed exclusion limit in the (BRe,BRtau) plane on the cross section for a 1600 GeV stop. All limits are computed at 95% CL.

The chosen signal region in the (BRe,BRtau) plane with the best expected exclusion on the cross section for a 1350 GeV stop. All limits are computed at 95% CL.

The chosen signal region in the (BRe,BRtau) plane with the best expected exclusion on the cross section for a 1400 GeV stop. All limits are computed at 95% CL.

The chosen signal region in the (BRe,BRtau) plane with the best expected exclusion on the cross section for a 1450 GeV stop. All limits are computed at 95% CL.

The chosen signal region in the (BRe,BRtau) plane with the best expected exclusion on the cross section for a 1500 GeV stop. All limits are computed at 95% CL.

The chosen signal region in the (BRe,BRtau) plane with the best expected exclusion on the cross section for a 1550 GeV stop. All limits are computed at 95% CL.

The chosen signal region in the (BRe,BRtau) plane with the best expected exclusion on the cross section for a 1600 GeV stop. All limits are computed at 95% CL.

The chosen signal region in the (BRe,BRtau) plane with the best expected exclusion on the cross section for a 600 GeV stop. All limits are computed at 95% CL.

The chosen signal region in the (BRe,BRtau) plane with the best expected exclusion on the cross section for a 700 GeV stop. All limits are computed at 95% CL.

The chosen signal region in the (BRe,BRtau) plane with the best expected exclusion on the cross section for a 800 GeV stop. All limits are computed at 95% CL.

The chosen signal region in the (BRe,BRtau) plane with the best expected exclusion on the cross section for a 900 GeV stop. All limits are computed at 95% CL.

The chosen signal region in the (BRe,BRtau) plane with the best expected exclusion on the cross section for a 1000 GeV stop. All limits are computed at 95% CL.

The chosen signal region in the (BRe,BRtau) plane with the best expected exclusion on the cross section for a 1050 GeV stop. All limits are computed at 95% CL.

The chosen signal region in the (BRe,BRtau) plane with the best expected exclusion on the cross section for a 1100 GeV stop. All limits are computed at 95% CL.

The chosen signal region in the (BRe,BRtau) plane with the best expected exclusion on the cross section for a 1150 GeV stop. All limits are computed at 95% CL.

The chosen signal region in the (BRe,BRtau) plane with the best expected exclusion on the cross section for a 1200 GeV stop. All limits are computed at 95% CL.

The chosen signal region in the (BRe,BRtau) plane with the best expected exclusion on the cross section for a 1250 GeV stop. All limits are computed at 95% CL.

The chosen signal region in the (BRe,BRtau) plane with the best expected exclusion on the cross section for a 1300 GeV stop. All limits are computed at 95% CL.

Dijet mass distributions for data in the (a) WW signal region.

Dijet mass distributions for data in the (b) WZ signal region.

Dijet mass distributions for data in the (c) ZZ signal region.

Dijet mass distributions for data in the (d) WZ+WW signal region.

Dijet mass distributions for data in the (e) WW+ZZ signal region.

Upper limits at the 95% CL on the cross section times branching ratio for WW+WZ production as a function of V' mass

Upper limits at the 95% CL on the cross section times branching ratio for WW+ZZ production as a function of GKK mass for the bulk RS model with k/M̄Pl=1.

Upper limits at the 95% CL on the cross section times branching ratio for (c) WW+ZZ production as a function of scalar mass.

Upper limits at the 95% CL on the cross section times branching ratio for WW production as a function of V' mass

Upper limits at the 95% CL on the cross section times branching ratio for WZ production as a function of V' mass

Upper limits at the 95% CL on the cross section times branching ratio for WW production as a function of GKK mass for the bulk RS model with k/M̄Pl=0.5.

Upper limits at the 95% CL on the cross section times branching ratio for ZZ production as a function of GKK mass for the bulk RS model with k/M̄Pl=0.5.

Upper limits at the 95% CL on the cross section times branching ratio for WW production as a function of GKK mass for the bulk RS model with k/M̄Pl=1.

Upper limits at the 95% CL on the cross section times branching ratio for ZZ production as a function of GKK mass for the bulk RS model with k/M̄Pl=1.

Upper limits at the 95% CL on the cross section times branching ratio for WW+ZZ production as a function of GKK mass for the bulk RS model with k/M̄Pl=0.5.

Upper limits at the 95% CL on the cross section times branching ratio for WW production as a function of scalar mass.

Upper limits at the 95% CL on the cross section times branching ratio for ZZ production as a function of scalar mass.

Observed 95% CL limit for the pMSSM grid when the signal cross section is decreased by one standard deviation.

Expected 95% CL limit for the pMSSM grid.

Expected 95% CL limit for the pMSSM grid with an up variation of the uncertainties.

Expected 95% CL limit for the pMSSM grid with a down variation of the uncertainties.

Observed 95% CL limit for the 2Step grid.

Observed 95% CL limit for the 2Step grid when the signal cross section is increased by one standard deviation.

Observed 95% CL limit for the 2Step grid when the signal cross section is decreased by one standard deviation.

Expected 95% CL limit for the 2Step grid.

Expected 95% CL limit for the 2Step grid with an up variation of the uncertainties.

Expected 95% CL limit for the 2Step grid with a down variation of the uncertainties.

Observed 95% CL limit for the gtt off-shell grid.

Observed 95% CL limit for the gtt off-shell grid when the signal cross section is increased by one standard deviation.

Observed 95% CL limit for the gtt off-shell grid when the signal cross section is decreased by one standard deviation.

Expected 95% CL limit for the gtt off-shell grid.

Expected 95% CL limit for the gtt off-shell grid with an up variation of the uncertainties.

Expected 95% CL limit for the gtt off-shell grid with a down variation of the uncertainties.

Observed 95% CL limit for the RPV grid.

Observed 95% CL limit for the RPV grid when the signal cross section is increased by one standard deviation.

Observed 95% CL limit for the RPV grid when the signal cross section is decreased by one standard deviation.

Expected 95% CL limit for the RPV grid.

Expected 95% CL limit for the RPV grid with an up variation of the uncertainties.

Expected 95% CL limit for the RPV grid with a down variation of the uncertainties.

Number of signal events expected for 36.1 fb$^{-1}$ at different stages of the event selection for the signal region SR-7j80-0b in a pMSSM inspired model where m($\tilde{g}$) = 1400 GeV and m($\tilde{\chi}_{0}^{1}$) = 200 GeV.

Number of signal events expected for 36.1 fb$^{-1}$ at different stages of the event selection for the signal region SR-7j80-1b in a pMSSM inspired model where m($\tilde{g}$) = 1400 GeV and m($\tilde{\chi}_{0}^{1}$) = 200 GeV.

Number of signal events expected for 36.1 fb$^{-1}$ at different stages of the event selection for the signal region SR-7j80-2b in a pMSSM inspired model where m($\tilde{g}$) = 1400 GeV and m($\tilde{\chi}_{0}^{1}$) = 200 GeV.

Number of signal events expected for 36.1 fb$^{-1}$ at different stages of the event selection for the signal region SR-8j80-0b in a pMSSM inspired model where m($\tilde{g}$) = 1400 GeV and m($\tilde{\chi}_{0}^{1}$) = 200 GeV.

Number of signal events expected for 36.1 fb$^{-1}$ at different stages of the event selection for the signal region SR-8j80-1b in a pMSSM inspired model where m($\tilde{g}$) = 1400 GeV and m($\tilde{\chi}_{0}^{1}$) = 200 GeV.

Number of signal events expected for 36.1 fb$^{-1}$ at different stages of the event selection for the signal region SR-8j80-2b in a pMSSM inspired model where m($\tilde{g}$) = 1400 GeV and m($\tilde{\chi}_{0}^{1}$) = 200 GeV.

Number of signal events expected for 36.1 fb$^{-1}$ at different stages of the event selection for the signal region SR-9j80-0b in a pMSSM inspired model where m($\tilde{g}$) = 1400 GeV and m($\tilde{\chi}_{0}^{1}$) = 200 GeV.

Number of signal events expected for 36.1 fb$^{-1}$ at different stages of the event selection for the signal region SR-9j80-1b in a pMSSM inspired model where m($\tilde{g}$) = 1400 GeV and m($\tilde{\chi}_{0}^{1}$) = 200 GeV.

Number of signal events expected for 36.1 fb$^{-1}$ at different stages of the event selection for the signal region SR-9j80-2b in a pMSSM inspired model where m($\tilde{g}$) = 1400 GeV and m($\tilde{\chi}_{0}^{1}$) = 200 GeV.

Number of signal events expected for 36.1 fb$^{-1}$ at different stages of the event selection for the signal region SR-8j50-0b in a pMSSM inspired model where m($\tilde{g}$) = 1400 GeV and m($\tilde{\chi}_{0}^{1}$) = 200 GeV.

Number of signal events expected for 36.1 fb$^{-1}$ at different stages of the event selection for the signal region SR-8j50-0b in a pMSSM inspired model where m($\tilde{g}$) = 1400 GeV and m($\tilde{\chi}_{0}^{1}$) = 200 GeV.

Number of signal events expected for 36.1 fb$^{-1}$ at different stages of the event selection for the signal region SR-8j50-1b in a pMSSM inspired model where m($\tilde{g}$) = 1400 GeV and m($\tilde{\chi}_{0}^{1}$) = 200 GeV.

Number of signal events expected for 36.1 fb$^{-1}$ at different stages of the event selection for the signal region SR-8j50-2b in a pMSSM inspired model where m($\tilde{g}$) = 1400 GeV and m($\tilde{\chi}_{0}^{1}$) = 200 GeV.

Number of signal events expected for 36.1 fb$^{-1}$ at different stages of the event selection for the signal region SR-9j50-0b in a pMSSM inspired model where m($\tilde{g}$) = 1400 GeV and m($\tilde{\chi}_{0}^{1}$) = 200 GeV.

Number of signal events expected for 36.1 fb$^{-1}$ at different stages of the event selection for the signal region SR-9j50-0b in a pMSSM inspired model where m($\tilde{g}$) = 1400 GeV and m($\tilde{\chi}_{0}^{1}$) = 200 GeV.

Number of signal events expected for 36.1 fb$^{-1}$ at different stages of the event selection for the signal region SR-9j50-1b in a pMSSM inspired model where m($\tilde{g}$) = 1400 GeV and m($\tilde{\chi}_{0}^{1}$) = 200 GeV.

Number of signal events expected for 36.1 fb$^{-1}$ at different stages of the event selection for the signal region SR-9j50-2b in a pMSSM inspired model where m($\tilde{g}$) = 1400 GeV and m($\tilde{\chi}_{0}^{1}$) = 200 GeV.

Number of signal events expected for 36.1 fb$^{-1}$ at different stages of the event selection for the signal region SR-10j50-0b in a pMSSM inspired model where m($\tilde{g}$) = 1400 GeV and m($\tilde{\chi}_{0}^{1}$) = 200 GeV.

Number of signal events expected for 36.1 fb$^{-1}$ at different stages of the event selection for the signal region SR-10j50-0b in a pMSSM inspired model where m($\tilde{g}$) = 1400 GeV and m($\tilde{\chi}_{0}^{1}$) = 200 GeV.

Number of signal events expected for 36.1 fb$^{-1}$ at different stages of the event selection for the signal region SR-10j50-1b in a pMSSM inspired model where m($\tilde{g}$) = 1400 GeV and m($\tilde{\chi}_{0}^{1}$) = 200 GeV.

Number of signal events expected for 36.1 fb$^{-1}$ at different stages of the event selection for the signal region SR-10j50-2b in a pMSSM inspired model where m($\tilde{g}$) = 1400 GeV and m($\tilde{\chi}_{0}^{1}$) = 200 GeV.

Number of signal events expected for 36.1 fb$^{-1}$ at different stages of the event selection for the signal region SR-11j50-0b in a pMSSM inspired model where m($\tilde{g}$) = 1400 GeV and m($\tilde{\chi}_{0}^{1}$) = 200 GeV.

Number of signal events expected for 36.1 fb$^{-1}$ at different stages of the event selection for the signal region SR-11j50-1b in a pMSSM inspired model where m($\tilde{g}$) = 1400 GeV and m($\tilde{\chi}_{0}^{1}$) = 200 GeV.

Number of signal events expected for 36.1 fb$^{-1}$ at different stages of the event selection for the signal region SR-11j50-2b in a pMSSM inspired model where m($\tilde{g}$) = 1400 GeV and m($\tilde{\chi}_{0}^{1}$) = 200 GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region SR-8j50-0b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region SR-8j50-1b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region SR-8j50-2b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region SR-9j50-0b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region SR-9j50-1b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region SR-9j50-2b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region SR-10j50-0b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region SR-10j50-1b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region SR-10j50-2b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region SR-11j50-0b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region SR-11j50-1b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region SR-11j50-2b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region SR-7j80-0b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region SR-7j80-1b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region SR-7j80-2b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region SR-8j80-0b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region SR-8j80-1b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region SR-8j80-2b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region SR-9j80-0b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region SR-9j80-1b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region SR-9j80-2b. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region SR-8j50-0b-MJ340. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region SR-8j50-0b-MJ500. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region SR-9j50-0b-MJ340. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region SR-9j50-0b-MJ500. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region SR-10j50-0b-MJ340. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

$E_{\mathrm{T}}^{\mathrm{miss}} / \sqrt{H_{\mathrm{T}}}$ distribution in signal region SR-10j50-0b-MJ500. Two benchmark signal models are overlaid on the plot for comparison. Labelled `pMSSM' and `2-step', they show signal distributions from the example SUSY models (as described in the paper): a pMSSM slice model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{\pm}}$) = (1300, 200) GeV and a cascade decay model with ($m \tilde{g}$, $m \tilde{\chi_{1}^{0}}$) = (1300, 200) GeV.

Degree of multijet closure for signal and vaidation regions (prior to the leptonic background fit) for the flavour stream with no b-jet requirement and a minimum transverse momentum of 50 GeV. The solid lines are the pre-fit predicted numbers of events and the points are the observed numbers. The blue hatched band shows only the statistical (MC and data) uncertainty on the background estimate. The template closure uncertainty for each SR bin is given by the maximal deviation of data from prediction in any non-SR bin to its left on this plot (although those for 80 GeV regions are independent of deviations in 50 GeV regions).

Degree of multijet closure for signal and vaidation regions (prior to the leptonic background fit) for the flavour stream with one inclusive b-jet required and a minimum transverse momentum of 50 GeV. The solid lines are the pre-fit predicted numbers of events and the points are the observed numbers. The blue hatched band shows only the statistical (MC and data) uncertainty on the background estimate. The template closure uncertainty for each SR bin is given by the maximal deviation of data from prediction in any non-SR bin to its left on this plot (although those for 80 GeV regions are independent of deviations in 50 GeV regions).

Degree of multijet closure for signal and vaidation regions (prior to the leptonic background fit) for the flavour stream with two inclusive b-jets required and a minimum transverse momentum of 50 GeV. The solid lines are the pre-fit predicted numbers of events and the points are the observed numbers. The blue hatched band shows only the statistical (MC and data) uncertainty on the background estimate. The template closure uncertainty for each SR bin is given by the maximal deviation of data from prediction in any non-SR bin to its left on this plot (although those for 80 GeV regions are independent of deviations in 50 GeV regions).

Degree of multijet closure for signal and vaidation regions (prior to the leptonic background fit) for the flavour stream with no b-jet requirement and a minimum transverse momentum of 80 GeV. The solid lines are the pre-fit predicted numbers of events and the points are the observed numbers. The blue hatched band shows only the statistical (MC and data) uncertainty on the background estimate. The template closure uncertainty for each SR bin is given by the maximal deviation of data from prediction in any non-SR bin to its left on this plot (although those for 80 GeV regions are independent of deviations in 50 GeV regions).

Degree of multijet closure for signal and vaidation regions (prior to the leptonic background fit) for the flavour stream with one inclusive b-jet required and a minimum transverse momentum of 80 GeV. The solid lines are the pre-fit predicted numbers of events and the points are the observed numbers. The blue hatched band shows only the statistical (MC and data) uncertainty on the background estimate. The template closure uncertainty for each SR bin is given by the maximal deviation of data from prediction in any non-SR bin to its left on this plot (although those for 80 GeV regions are independent of deviations in 50 GeV regions).

Degree of multijet closure for signal and vaidation regions (prior to the leptonic background fit) for the flavour stream with two inclusive b-jets required and a minimum transverse momentum of 80 GeV. The solid lines are the pre-fit predicted numbers of events and the points are the observed numbers. The blue hatched band shows only the statistical (MC and data) uncertainty on the background estimate. The template closure uncertainty for each SR bin is given by the maximal deviation of data from prediction in any non-SR bin to its left on this plot (although those for 80 GeV regions are independent of deviations in 50 GeV regions).

Degree of multijet closure for signal and vaidation regions (prior to the leptonic background fit) for the fat-jet stream with MJSigma above 340 GeV. The solid lines are the pre-fit predicted numbers of events and the points are the observed numbers. The blue hatched band shows only the statistical (MC and data) uncertainty on the background estimate. The template closure uncertainty for each SR bin is given by the maximal deviation of data from prediction in any non-SR bin to its left on this plot (although those for 80 GeV regions are independent of deviations in 50 GeV regions).

Degree of multijet closure for signal and vaidation regions (prior to the leptonic background fit) for the fat-jet stream with MJSigma above 500 GeV. The solid lines are the pre-fit predicted numbers of events and the points are the observed numbers. The blue hatched band shows only the statistical (MC and data) uncertainty on the background estimate. The template closure uncertainty for each SR bin is given by the maximal deviation of data from prediction in any non-SR bin to its left on this plot (although those for 80 GeV regions are independent of deviations in 50 GeV regions).

The best-expected signal region and the corresponding best-observed and best-expected CLs values for the 2Step grid.

The best-expected signal region and the corresponding best-observed and best-expected CLs values for the 2Step grid.

The best-expected signal region and the corresponding best-observed and best-expected CLs values for the 2Step grid.

The best-expected signal region and the corresponding best-observed and best-expected CLs values for the pMSSM grid.

The best-expected signal region and the corresponding best-observed and best-expected CLs values for the pMSSM grid.

The best-expected signal region and the corresponding best-observed and best-expected CLs values for the pMSSM grid.

The best-expected signal region and the corresponding best-observed and best-expected CLs values for the RPV grid.

The best-expected signal region and the corresponding best-observed and best-expected CLs values for the RPV grid.

The best-expected signal region and the corresponding best-observed and best-expected CLs values for the RPV grid.

The best-expected signal region and the corresponding best-observed and best-expected CLs values for the gtt off-shell grid.

The best-expected signal region and the corresponding best-observed and best-expected CLs values for the gtt off-shell grid.

The best-expected signal region and the corresponding best-observed and best-expected CLs values for the gtt off-shell grid.

95% CLs observed upper limit on model cross-section (in fb) for 2Step signal points for the best-expected signal region.

95% CLs observed upper limit on model cross-section (in fb) for RPV signal points for the best-expected signal region.

95% CLs observed upper limit on model cross-section (in fb) for gtt off-shell signal points for the best-expected signal region.

Performance of the SR-8j50-0b for the 2Step grid: fractional acceptance; fractional efficiency.

Performance of the SR-8j50-0b-MJ340 for the 2Step grid: fractional acceptance; fractional efficiency.

Performance of the SR-8j50-0b-MJ500 for the 2Step grid: fractional acceptance; fractional efficiency.

Performance of the SR-8j50-1b for the 2Step grid: fractional acceptance; fractional efficiency.

Performance of the SR-8j50-2b for the 2Step grid: fractional acceptance; fractional efficiency.

Performance of the SR-9j50-0b for the 2Step grid: fractional acceptance; fractional efficiency.

Performance of the SR-9j50-0b-MJ340 for the 2Step grid: fractional acceptance; fractional efficiency.

Performance of the SR-9j50-0b-MJ500 for the 2Step grid: fractional acceptance; fractional efficiency.

Performance of the SR-9j50-1b for the 2Step grid: fractional acceptance; fractional efficiency.

Performance of the SR-9j50-2b for the 2Step grid: fractional acceptance; fractional efficiency.

Performance of the SR-10j50-0b for the 2Step grid: fractional acceptance; fractional efficiency.

Performance of the SR-10j50-0b-MJ340 for the 2Step grid: fractional acceptance; fractional efficiency.

Performance of the SR-10j50-0b-MJ500 for the 2Step grid: fractional acceptance; fractional efficiency.

Performance of the SR-10j50-1b for the 2Step grid: fractional acceptance; fractional efficiency.

Performance of the SR-10j50-2b for the 2Step grid: fractional acceptance; fractional efficiency.

Performance of the SR-11j50-0b for the 2Step grid: fractional acceptance; fractional efficiency.

Performance of the SR-11j50-1b for the 2Step grid: fractional acceptance; fractional efficiency.

Performance of the SR-11j50-2b for the 2Step grid: fractional acceptance; fractional efficiency.

Performance of the SR-7j80-0b for the 2Step grid: fractional acceptance; fractional efficiency.

Performance of the SR-7j80-1b for the 2Step grid: fractional acceptance; fractional efficiency.

Performance of the SR-7j80-2b for the 2Step grid: fractional acceptance; fractional efficiency.

Performance of the SR-8j80-0b for the 2Step grid: fractional acceptance; fractional efficiency.

Performance of the SR-8j80-1b for the 2Step grid: fractional acceptance; fractional efficiency.

Performance of the SR-8j80-2b for the 2Step grid: fractional acceptance; fractional efficiency.

Performance of the SR-9j80-0b for the 2Step grid: fractional acceptance; fractional efficiency.

Performance of the SR-9j80-1b for the 2Step grid: fractional acceptance; fractional efficiency.

Performance of the SR-9j80-2b for the 2Step grid: fractional acceptance; fractional efficiency.

Upper limits on the fiducial cross section times branching ratio to two photons at centre-of-mass energy of 13 TeV of a narrow-width (Γ_X = 4 MeV) spin-0 resonance as a function of its mass m_X.

Upper limits on the production cross section times branching ratio to two photons at centre-of-mass energy of 13 TeV of the lightest KK graviton as a function of its mass m_G* for k/MPl=0.01.

Upper limits on the production cross section times branching ratio to two photons at centre-of-mass energy of 13 TeV of the lightest KK graviton as a function of its mass m_G* for k/MPl=0.10. Only the limits derived with the asymptotic formulae are quoted here.

Upper limits on the production cross section times branching ratio to two photons at centre-of-mass energy of 13 TeV of the lightest KK graviton as a function of its mass m_G* for k/MPl=0.20.

Upper limits on the production cross section times branching ratio to two photons at centre-of-mass energy of 13 TeV of the lightest KK graviton as a function of its mass m_G* for k/MPl=0.30.

Distribution of the diphoton invariant mass for data events passing the spin-0 selection. The data points are plotted at the centre of each bin. The error bars indicate statistical uncertainties. The distribution consists of not only genuine diphoton events, but also events where at least one photon is faked by a jet. There is no data event above 2700 GeV.

Distribution of the diphoton invariant mass for data events passing the spin-2 selection. The data points are plotted at the centre of each bin. The error bars indicate statistical uncertainties. The distribution consists of not only genuine diphoton events, but also events where at least one photon is faked by a jet. There is no data event above 2700 GeV.

Fiducial acceptance and product of the fiducial acceptance times the combined reconstruction and identification efficiency as a function of the mass for a spin-0 resonant signal with narrow width under the spin-0 selection. Only MC statistical uncertainties are included.

Fiducial acceptance and product of the fiducial acceptance times the combined reconstruction and identification efficiency as a function of the mass for a spin-2 resonant signal with k/MPl=0.1 under the spin-2 selection. The definition of fiducial region is provided in JHEP 09 (2016) 001. Only MC statistical uncertainties are included.

Fiducial acceptance and product of the fiducial acceptance times the combined reconstruction and identification efficiency as a function of the mass for an ADD excess (GRW formalism), calculated as the difference between the sum of ADD and SM contributions (including their interference) and the SM contribution, under the spin-2 selection in the signal region M(gamma gamma)>2240 GeV. The definition of fiducial region is provided in JHEP 09 (2016) 001, and M(gamma gamma)>2240 GeV is required in addition. Only MC statistical uncertainties are included.

Selection acceptance times efficiency as a function of WPRIME mass at truth level for left- and right-handed WPRIME MC. The total efficiency curves correspond to the sum of the efficiencies of the one b-tag and two b-tag categories.

Cutflow (efficiency with respect to total number of events in %) for several WPRIME masses in the left-handed and in the right-handed model for hadronic top-quark decays.

Overview of the signal acceptance times efficiency (A x Eff) for the hadronic top-quark decay for both categories, the predicted cross section times branching ratio to TOP BOTTOM, as well as the observed 95% CL limit on the cross section times branching ratio for several WPRIME masses in the left-handed and in the right-handed model. The expected signal event yields in the two categories for 20.3 fb-1 of 8 TeV data are also shown.