The distribution of the fit discriminant mTVH in the 3+ b-tag signal region of the vvbb channel for the mH=300 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The m(tt) distribution in the Lhi3_Zin region of the lltt channel. <br><br><a href="?table=overview">return to overview</a>

The m(bb) distribution in the 2 b-tag 0L region of the vvbb channel. <br><br><a href="?table=overview">return to overview</a>

The m(bb) distribution in the 3+ b-tag 0L region of the vvbb channel. <br><br><a href="?table=overview">return to overview</a>

Expected and observed upper limits at 95% CL on sigma(gg->A)*B(A->ZH)*B(H->tt) in the (mA,mH) plane. The limits are shown for tanbeta=1. The tanbeta value is relevant only for the choice of the A boson width. <br><a href="?table=overview">return to overview</a>

Expected and observed upper limits at 95% CL on sigma(bbA)*B(A->ZH)*B(H->tt)$ in the (mA,mH) plane. The limits are shown for tanbeta=10. The tanbeta value is relevant only for the choice of the A boson width. <br><br><a href="?table=overview">return to overview</a>

Expected and observed upper limits at 95% CL on sigma(gg->A)*B(A->ZH)*B(H->bb) in the (mA,mH) plane. The limits are shown for tanbeta=1. The tanbeta value is relevant only for the choice of the A boson width. <br><br><a href="?table=overview">return to overview</a>

Expected and observed upper limits at 95% CL on sigma(bbA)*B(A->ZH)*B(H->bb) in the (mA,mH) plane. The limits are shown for tanbeta=10. The tanbeta value is relevant only for the choice of the A boson width. <br><br><a href="?table=overview">return to overview</a>

Expected and observed upper limits at 95% CL on sigma(gg->A)*B(A->ZH)*B(H->tt) in he (mA,mH) plane. The limits are shown for tanbeta=0.5. The tanbeta value is relevant only for the choice of the A boson width. <br><a href="?table=overview">return to overview</a>

Expected and observed upper limits at 95% CL on sigma(gg->A)*B(A->ZH)*B(H->tt) in the (mA,mH) plane. The limits are shown for tanbeta=5. The tanbeta value is relevant only for the choice of the A boson width. <br><a href="?table=overview">return to overview</a>

Expected and observed upper limits at 95% CL on sigma(bbA)*B(A->ZH)*B(H->tt) in the (mA,mH) plane. The limits are shown for tanbeta=1. The tanbeta value is relevant only for the choice of the A boson width. <br><br><a href="?table=overview">return to overview</a>

Expected and observed upper limits at 95% CL on sigma(bbA)*B(A->ZH)*B(H->tt)$ in the (mA,mH) plane. The limits are shown for tanbeta=5. The tanbeta value is relevant only for the choice of the A boson width. <br><br><a href="?table=overview">return to overview</a>

Expected and observed upper limits at 95% CL on sigma(gg->A)*B(A->ZH)*B(H->bb) in the (mA,mH) plane. The limits are shown for tanbeta=0.5. The tanbeta value is relevant only for the choice of the A boson width. <br><br><a href="?table=overview">return to overview</a>

Expected and observed upper limits at 95% CL on sigma(gg->A)*B(A->ZH)*B(H->bb) in the (mA,mH) plane. The limits are shown for tanbeta=5. The tanbeta value is relevant only for the choice of the A boson width. <br><br><a href="?table=overview">return to overview</a>

Expected and observed upper limits at 95% CL on sigma(bbA)*B(A->ZH)*B(H->bb) in the (mA,mH) plane. The limits are shown for tanbeta=1. The tanbeta value is relevant only for the choice of the A boson width. <br><br><a href="?table=overview">return to overview</a>

Expected and observed upper limits at 95% CL on sigma(bbA)*B(A->ZH)*B(H->bb) in the (mA,mH) plane. The limits are shown for tanbeta=5. The tanbeta value is relevant only for the choice of the A boson width. <br><br><a href="?table=overview">return to overview</a>

Expected and observed upper limits at 95% CL on sigma(bbA)*B(A->ZH)*B(H->bb) in the (mA,mH) plane. The limits are shown for tanbeta=20. The tanbeta value is relevant only for the choice of the A boson width. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant m(lltt)-m(tt) in the signal region of the lltt channel for the mH=450 GeV hypothesis with the bbA signal shown for comparison. <br><br><a href="?table=overview">return to overview</a>

The m(tt) distribution in the L3hi_Zin region of the lltt channel for the mH=450 GeV hypothesis with the bbA signal shown for comparison. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant m(lltt)-m(tt) in the signal region of the lltt channel for the mH=350 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant m(lltt)-m(tt) in the signal region of the lltt channel for the mH=400 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant m(lltt)-m(tt) in the signal region of the lltt channel for the mH=500 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant m(lltt)-m(tt) in the signal region of the lltt channel for the mH=550 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant m(lltt)-m(tt) in the signal region of the lltt channel for the mH=600 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant m(lltt)-m(tt) in the signal region of the lltt channel for the mH=700 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant m(lltt)-m(tt) in the signal region of the lltt channel for the mH=800 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 2 b-tag signal region of the vvbb channel for the mH=130 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 2 b-tag signal region of the vvbb channel for the mH=150 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 2 b-tag signal region of the vvbb channel for the mH=200 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 2 b-tag signal region of the vvbb channel for the mH=250 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 2 b-tag signal region of the vvbb channel for the mH=350 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 2 b-tag signal region of the vvbb channel for the mH=400 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 2 b-tag signal region of the vvbb channel for the mH=450 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 2 b-tag signal region of the vvbb channel for the mH=500 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 2 b-tag signal region of the vvbb channel for the mH=600 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 2 b-tag signal region of the vvbb channel for the mH=700 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 2 b-tag signal region of the vvbb channel for the mH=800 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH i+ the 3p b-tag signal region of the vvbb channel for the mH=130 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 3+ b-tag signal region of the vvbb channel for the mH=150 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 3+ b-tag signal region of the vvbb channel for the mH=200 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 3+ b-tag signal region of the vvbb channel for the mH=250 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 3+ b-tag signal region of the vvbb channel for the mH=350 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 3+ b-tag signal region of the vvbb channel for the mH=400 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 3+ b-tag signal region of the vvbb channel for the mH=450 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 3+ b-tag signal region of the vvbb channel for the mH=500 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 3+ b-tag signal region of the vvbb channel for the mH=600 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 3+ b-tag signal region of the vvbb channel for the mH=700 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 3+ b-tag signal region of the vvbb channel for the mH=800 GeV hypothesis. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 2 b-tag 2L region of the vvbb channel. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 2 b-tag em region of the vvbb channel. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 3+ b-tag 2L region of the vvbb channel. <br><br><a href="?table=overview">return to overview</a>

The distribution of the fit discriminant mTVH in the 3+ b-tag em region of the vvbb channel. <br><br><a href="?table=overview">return to overview</a>

Local significance for lltt ggF signals for tanbeta=1. <br><br><a href="?table=overview">return to overview</a>

Local significance for lltt bbA signals for tanbeta=20. <br><br><a href="?table=overview">return to overview</a>

Local significance for vvbb ggF signals for tanbeta=1. <br><br><a href="?table=overview">return to overview</a>

Local significance for vvbb bbA signals for tanbeta=20. <br><br><a href="?table=overview">return to overview</a>

Signal acceptance times efficiency for lltt ggF signals. <br><br><a href="?table=overview">return to overview</a>

Signal acceptance times efficiency for lltt bbA signals. <br><br><a href="?table=overview">return to overview</a>

Signal acceptance times efficiency for vvbb ggF signals. <br><br><a href="?table=overview">return to overview</a>

Signal acceptance times efficiency for vvbb bbA signals. <br><br><a href="?table=overview">return to overview</a>

The pT distribution of the third lepton in the L3hi_Zin region of the lltt channel. <br><br><a href="?table=overview">return to overview</a>

The eta distribution of the H boson candidate in the rest frame of the A boson candidate in the L3hi_Zin region of the lltt channel. <br><br><a href="?table=overview">return to overview</a>

The ETmiss distribution in the 2 b-tag 0L region of the vvbb channel. <br><br><a href="?table=overview">return to overview</a>

The m_top_near distribution in the 2 b-tag 0L region of the vvbb channel. <br><br><a href="?table=overview">return to overview</a>

The ETmiss distribution in the 3+ b-tag 0L region of the vvbb channel. <br><br><a href="?table=overview">return to overview</a>

The m_top_near distribution in the 3+ b-tag 0L region of the vvbb channel. <br><br><a href="?table=overview">return to overview</a>

Distribution of $m_{\mathrm{JJ}}$ in the SR (black points) compared to the background-only expectation, where the shape is derived from the CR and the normalisation is fitted in the SR.

Observed and expected limits at 95% CL on the cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model A (see the text for details). The darker and lighter shaded bands around the expected limits represent the $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty range, respectively.

The predicted theory cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model A (see the text for details). In the theory model, a universal coupling of the $Z^\prime$ to each SM quark $g_q$ is assumed and the dark quarks are all represented by one species with an effective coupling $g_{q_d}$ to the $Z^\prime$. For model A, the predicted cross-section is shown for $g_q=0.05$ and $g_{q_d} = 0.2$.

Observed and expected limits at 95% CL on the cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model B (see the text for details). The darker and lighter shaded bands around the expected limits represent the $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty range, respectively.

The predicted theory cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model B (see the text for details). In the theory model, a universal coupling of the $Z^\prime$ to each SM quark $g_q$ is assumed and the dark quarks are all represented by one species with an effective coupling $g_{q_d}$ to the $Z^\prime$. For model B, the predicted cross-section is shown for $g_q=0.15$ and $g_{q_d} = 0.5$.

Observed and expected limits at 95% CL on the cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model C (see the text for details). The darker and lighter shaded bands around the expected limits represent the $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty range, respectively.

The predicted theory cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model C (see the text for details). In the theory model, a universal coupling of the $Z^\prime$ to each SM quark $g_q$ is assumed and the dark quarks are all represented by one species with an effective coupling $g_{q_d}$ to the $Z^\prime$. For model C, the predicted cross-section is shown for $g_q=0.15$ and $g_{q_d}= 0.5$.

Observed and expected limits at 95% CL on the cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model D (see the text for details). The darker and lighter shaded bands around the expected limits represent the $\pm 1 \sigma$ and $\pm 2 \sigma$ uncertainty range, respectively.

The predicted theory cross-section times branching ratio for the production of dark quarks through a $Z^\prime$ as a function of the $Z^\prime$ mass for model D (see the text for details). In the theory model, a universal coupling of the $Z^\prime$ to each SM quark $g_q$ is assumed and the dark quarks are all represented by one species with an effective coupling $g_{q_d}$ to the $Z^\prime$. For model D, the predicted cross-section is shown for $g_q=0.05$ and $g_{q_d} = 0.2$.

Relative efficiency (in %) of each stage of the analysis selection with respect to the previous one for models A through D and $m_{Z^\prime}=2.5$ TeV. Approximately 90000 events were generated per signal mass hypothesis.

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

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

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

Likelihood observed contours at 68% CL (solid line) and 95% CL (dashed line) in the $c_{gghh}$ versus $c_{hhh}$ HEFT parameter space, with the remaining coefficient fixed to its SM value. The best-fit value, denoted by a cross in the plot, corresponds to the very first entry in the table.

Likelihood expected contours at 68% CL (teal shaded region) 95% CL (yellow shaded region) in the $c_{gghh}$ versus $c_{hhh}$ HEFT parameter space, with the remaining coefficient fixed to its SM value. In the plot, the SM prediction is indicated by the star.

Likelihood observed contours at 68% CL (solid line) and 95% CL (dashed line) in the $c_{tthh}$ versus $c_{hhh}$ HEFT parameter space, with the remaining coefficient fixed to its SM value. The best-fit value, denoted by a cross in the plot, corresponds to the very first entry in the table.

Likelihood exptected contours at 68% CL (teal shaded region) and 95% CL (yellow shaded region) in the $c_{tthh}$ versus $c_{hhh}$ HEFT parameter space, with the remaining coefficient fixed to its SM value. In the plot, the SM prediction is indicated by the star.

The observed (filled circles) and expected (hollow circles) 95% CL upper limits on the $HH$ ggF production cross-section in the Standard Model and for seven HEFT benchmark points. The teal and yellow colored bands indicate the $\pm 1\sigma$ and $\pm 2\sigma$ variations on the expected limit due to statistical and systematic uncertainties. The predicted cross-sections of each of the models under consideration are shown by the red crosses. The contribution from VBF production to the total $HH$ production cross-section is neglected.

Likelihood observed contours at 68% CL (solid line) and 95% CL (dashed line) in the $c_{{H}\boxed{}}$ versus $c_{H}$ SMEFT parameter space, with the remaining coefficient fixed to its SM value. The best-fit value, denoted by a cross in the plot, corresponds to the very first entry in the table.

Likelihood expected contours at 68% CL (teal shaded region) 95% CL (yellow shaded region) in the $c_{{H}\boxed{}}$ versus $c_{H}$ SMEFT parameter space, with the remaining coefficient fixed to its SM value. In the plot, the SM prediction is indicated by the star.

The acceptance times efficiency for the signal ggF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category. The other coupling modifiers affecting ggF $HH$ production are fixed to their SM predictions. The bands indicate the simulation statistical uncertainty.

The acceptance times efficiency for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ in each analysis category. When varying $\kappa_{\lambda}$, the other coupling modifiers affecting the VBF $HH$ production mode are set to their SM predictions. The bands indicate the simulation statistical uncertainty.

The acceptance times efficiency for the signal VBF $HH$ process as a function of the coupling modifier $\kappa_{2V}$ in each analysis category. When varying $\kappa_{2V}$, the other coupling modifiers affecting the VBF $HH$ production mode are set to their SM predictions. The bands indicate the simulation statistical uncertainty.

Expected yields of the signal $HH$ process as a function of the coupling modifier $\kappa_{\lambda}$ after applying the analysis preselection and in each analysis category after the final selection. The $HH$ yield is obtained considering both the $ggF$ and $VBF$ contributions. The other coupling modifiers affecting $HH$ production are fixed to their SM predictions. The bands indicate the simulation statistical uncertainty.

The acceptance times efficiency for the signal ggF $HH$ process in each analysis category as a function of the $c_{hhh}$ HEFT coefficients. The dashed lines denote values that are excluded at 95% CL. The bottom panels show the efficiency of the sum of all analysis categories.

The acceptance times efficiency for the signal ggF $HH$ process in each analysis category as a function of the $c_{tthh}$ HEFT coefficients. The dashed lines denote values that are excluded at 95% CL. The bottom panels show the efficiency of the sum of all analysis categories.

The acceptance times efficiency for the signal ggF $HH$ process in each analysis category as a function of the $c_{gghh}$ HEFT coefficients. The dashed lines denote values that are excluded at 95% CL. The bottom panels show the efficiency of the sum of all analysis categories.

The acceptance times efficiency for the signal ggF $HH$ process in each analysis category as a function of the $c_{H}$ SMEFT coefficients. The dashed lines denote values that are excluded at 95% CL. The bottom panels show the efficiency of the sum of all analysis categories.

The acceptance times efficiency for the signal ggF $HH$ process in each analysis category as a function of the $c_{{H}\boxed{}}$ SMEFT coefficients. The dashed lines denote values that are excluded at 95% CL. The bottom panels show the efficiency of the sum of all analysis categories.

The yield of the signal ggF $HH$ process in each analysis category as a function of the $c_{hhh}$ HEFT coefficients. The dashed lines denote values that are excluded at 95% CL. The bottom panels show the efficiency of the sum of all analysis categories.

The yield of the signal ggF $HH$ process in each analysis category as a function of the $c_{tthh}$ HEFT coefficients. The dashed lines denote values that are excluded at 95% CL. The bottom panels show the efficiency of the sum of all analysis categories.

The yield of the signal ggF $HH$ process in each analysis category as a function of the $c_{gghh}$ HEFT coefficients. The dashed lines denote values that are excluded at 95% CL. The bottom panels show the efficiency of the sum of all analysis categories.

The yield of the signal ggF $HH$ process in each analysis category as a function of the $c_{H}$ SMEFT coefficients. The dashed lines denote values that are excluded at 95% CL. The bottom panels show the efficiency of the sum of all analysis categories.

The yield of the signal ggF $HH$ process in each analysis category as a function of the $c_{{H}\boxed{}}$ SMEFT coefficients. The dashed lines denote values that are excluded at 95% CL. The bottom panels show the efficiency of the sum of all analysis categories.

Distributions of the diphoton invariant mass for events in data (black dots with error bars) compared to the sum of the expected signal and backgrounds (histograms) in the low mass 1 region.

Distributions of the diphoton invariant mass for events in data (black dots with error bars) compared to the sum of the expected signal and backgrounds (histograms) in the low mass 2 region.

Distributions of the diphoton invariant mass for events in data (black dots with error bars) compared to the sum of the expected signal and backgrounds (histograms) in the low mass 3 region.

Distributions of the diphoton invariant mass for events in data (black dots with error bars) compared to the sum of the expected signal and backgrounds (histograms) in the low mass 4 region.

Distributions of the diphoton invariant mass for events in data (black dots with error bars) compared to the sum of the expected signal and backgrounds (histograms) in the high mass 1 region.

Distributions of the diphoton invariant mass for events in data (black dots with error bars) compared to the sum of the expected signal and backgrounds (histograms) in the high mass 2 region.

Distributions of the diphoton invariant mass for events in data (black dots with error bars) compared to the sum of the expected signal and backgrounds (histograms) in the high mass 3 region.

Observed (solid line) value of $-2\ln\Lambda$ as a function of the $c_{hhh}$ HEFT coefficients when all other coupling modifiers are fixed to the SM value

Expected (dashed line) value of $-2\ln\Lambda$ as a function of the $c_{hhh}$ HEFT coefficients when all other coupling modifiers are fixed to the SM value

Observed (solid line) value of $-2\ln\Lambda$ as a function of the $c_{gghh}$ HEFT coefficients when all other coupling modifiers are fixed to the SM value

Expected (dashed line) value of $-2\ln\Lambda$ as a function of the $c_{gghh}$ HEFT coefficients when all other coupling modifiers are fixed to the SM value

Observed (solid line) value of $-2\ln\Lambda$ as a function of the $c_{tthh}$ HEFT coefficients when all other coupling modifiers are fixed to the SM value

Expected (dashed line) value of $-2\ln\Lambda$ as a function of the $c_{tthh}$ HEFT coefficients when all other coupling modifiers are fixed to the SM value

Observed (solid line) value of $-2\ln\Lambda$ as a function of the $c_{H}$ SMEFT coefficients when all other coupling modifiers are fixed to the SM value.

Expected (dashed line) value of $-2\ln\Lambda$ as a function of the $c_{H}$ SMEFT coefficients when all other coupling modifiers are fixed to the SM value.

Observed (solid line) value of $-2\ln\Lambda$ as a function of the $c_{{H}\boxed{}}$ SMEFT coefficients when all other coupling modifiers are fixed to the SM value.

Expected (dashed line) value of $-2\ln\Lambda$ as a function of the $c_{{H}\boxed{}}$ SMEFT coefficients when all other coupling modifiers are fixed to the SM value.

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

Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64822.9 $\pm$ 254.6, NPowHeg truth =634276.

Unfolded dR distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64273.2 $\pm$ 253.5, NPowHeg truth =338671.

Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 65362.4 $\pm$ 255.7 , NPowHeg truth =634214.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63486.8 $\pm$ 252.0, NPowHeg truth =333348.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64432.6 $\pm$ 253.8, NPowHeg truth =624059.

Unfolded $M(l^{+}\gamma)$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57095.9 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =76327.

Unfolded $M(l^{-}\gamma)$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57130 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =763273.

Unfolded $M(l^{+}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58310.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58320.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Unfolded dR distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57705.4 $\pm$ 254 , NPowHeg truth =301655, N Sherpa truth =763259.

Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58869.2 $\pm$ 256, NPowHeg truth =564660, NS herpa truth =813995.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57358.6 $\pm$ 254 , NPowHeg truth =298662, N Sherpa truth =756609.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58226.4 $\pm$ 254, NPowHeg truth =558949, N Sherpa truth =806988.

Unfolded $M(l^{+}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63717.4 $\pm$ 252.4, NPowHeg truth =338714.

Unfolded $M(l^{-}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63855.8 $\pm$ 252.7 , NPowHeg truth =338708.

Unfolded $M(l^{+}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64809.8 $\pm$ 254.6, NPowHeg truth =634285.

Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64822.9 $\pm$ 254.6, NPowHeg truth =634276.

Unfolded dR distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64273.2 $\pm$ 253.5, NPowHeg truth =338671.

Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 65362.4 $\pm$ 255.7 , NPowHeg truth =634214.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63486.8 $\pm$ 252.0, NPowHeg truth =333348.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64432.6 $\pm$ 253.8, NPowHeg truth =624059.

Unfolded $M(l^{+}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57095.9 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =76327.

Unfolded $M(l^{-}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57130 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =763273.

Unfolded $M(l^{+}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58310.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58320.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Unfolded dR distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57705.4 $\pm$ 254 , NPowHeg truth =301655, N Sherpa truth =763259.

Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58869.2 $\pm$ 256, NPowHeg truth =564660, NS herpa truth =813995.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57358.6 $\pm$ 254 , NPowHeg truth =298662, N Sherpa truth =756609.

Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58226.4 $\pm$ 254, NPowHeg truth =558949, N Sherpa truth =806988.

Truth $P_{T}^{Z}$ distribution after reweghting to data.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63717.4 $\pm$ 252.4, NPowHeg truth =338714.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63855.8 $\pm$ 252.7 , NPowHeg truth =338708.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64809.8 $\pm$ 254.6, NPowHeg truth =634285.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64822.9 $\pm$ 254.6, NPowHeg truth =634276.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64273.2 $\pm$ 253.5, NPowHeg truth =338671.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 65362.4 $\pm$ 255.7 , NPowHeg truth =634214.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63486.8 $\pm$ 252.0, NPowHeg truth =333348.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64432.6 $\pm$ 253.8, NPowHeg truth =624059.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57095.9 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =76327.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57130 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =763273.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58310.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58320.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57705.4 $\pm$ 254 , NPowHeg truth =301655, N Sherpa truth =763259.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58869.2 $\pm$ 256, NPowHeg truth =564660, NS herpa truth =813995.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57358.6 $\pm$ 254 , NPowHeg truth =298662, N Sherpa truth =756609.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with dressed leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58226.4 $\pm$ 254, NPowHeg truth =558949, N Sherpa truth =806988.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63717.4 $\pm$ 252.4, NPowHeg truth =338714.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63855.8 $\pm$ 252.7 , NPowHeg truth =338708.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64809.8 $\pm$ 254.6, NPowHeg truth =634285.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64822.9 $\pm$ 254.6, NPowHeg truth =634276.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64273.2 $\pm$ 253.5, NPowHeg truth =338671.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 65362.4 $\pm$ 255.7 , NPowHeg truth =634214.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 63486.8 $\pm$ 252.0, NPowHeg truth =333348.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>20$ GeV. Nexp.un f. = 64432.6 $\pm$ 253.8, NPowHeg truth =624059.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57095.9 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =76327.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57130 $\pm$ 254 , NPowHeg truth =301666, N Sherpa truth =763273.

Combined Covariance Matrix for Unfolded $M(l^{+}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58310.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Combined Covariance Matrix for Unfolded $M(l^{-}\gamma)$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58320.9 $\pm$ 245, NPowHeg truth =564672, N Sherpa truth =814016.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57705.4 $\pm$ 254 , NPowHeg truth =301655, N Sherpa truth =763259.

Combined Covariance Matrix for Unfolded dR distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58869.2 $\pm$ 256, NPowHeg truth =564660, NS herpa truth =813995.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to ee\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 57358.6 $\pm$ 254 , NPowHeg truth =298662, N Sherpa truth =756609.

Combined Covariance Matrix for Unfolded $P_{T}^{\gamma}$ distribution for $Z \to \mu\mu\gamma$ process with bare leptons and bkg subtraction. $M_{ll}>45$ GeV. Nexp.un f. = 58226.4 $\pm$ 254, NPowHeg truth =558949, N Sherpa truth =806988.

Unfolded $M(l^{+}\gamma\gamma)$ distribution for $Z \to ll\gamma\gamma$ process with dressed leptons.

Unfolded $M(l^{+}\gamma\gamma)$ distribution for $Z \to ll\gamma\gamma$ process with bare leptons.

Unfolded dR distribution for $Z \to ll\gamma\gamma$ process with dressed leptons for leading photon.

Unfolded dR distribution for $Z \to ll\gamma\gamma$ process with bare leptons for leading photon.

Unfolded dR distribution for $Z \to ll\gamma\gamma$ process with dressed leptons for second photon.

Unfolded dR distribution for $Z \to ll\gamma\gamma$ process with bare leptons for second photon.

Unfolded dRgg distribution for $Z \to ll\gamma\gamma$ process with dressed leptons.

Unfolded dRgg distribution for $Z \to ll\gamma\gamma$ process with bare leptons.

Unfolded $P_{T}^{\gamma 1}$ distribution for $Z \to ll\gamma\gamma$ process with dressed leptons for leading photon.

Unfolded $P_{T}^{\gamma 1}$ distribution for $Z \to ll\gamma\gamma$ process with bare leptons for leading photon.

Unfolded $P_{T}^{\gamma 2}$ distribution for $Z \to ll\gamma\gamma$ process with dressed leptons for second photon.

Unfolded $P_{T}^{\gamma 2}$ distribution for $Z \to ll\gamma\gamma$ process with bare leptons for second photon.

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

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

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

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

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

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

Measured $p_T$ cross sections in full-lepton phase space for 1.2 < |y| < 1.6.

Measured $p_T$ cross sections in full-lepton phase space for 1.6 < |y| < 2.0.

Measured $p_T$ cross sections in full-lepton phase space for 2.0 < |y| < 2.4.

Measured $p_T$ cross sections in full-lepton phase space for 2.4 < |y| < 2.8.

Measured $p_T$ cross sections in full-lepton phase space for 2.8 < |y| < 3.2.

Measured cross sections in full-lepton phase space as a function of |y|.

Normalised measured $p_T$ cross sections in full-lepton phase space for |y| < 1.6.

Predicted and observed yields as a function of $m_{4\ell}$ in the VBS-Suppressed region. Overflow events are included in the last bin of the distribution.

Relative systematic uncertainties in the unfolded differential cross section for 4ljj production in the VBS-Enhanced region as a function of m_{jj}

Relative systematic uncertainties in the unfolded differential cross section for 4ljj production in the VBS-Enhanced region as a function of m_{4\ell}

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $m_{jj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $m_{jj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $m_{4\ell}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $m_{4\ell}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $ p_{T,4l}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $ p_{T,4l}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $\Delta\phi_{jj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $\Delta\phi_{jj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $|\Delta y_{jj}|$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $|\Delta y_{jj}|$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $cos(\theta^{*}_{12})$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $cos(\theta^{*}_{12})$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $cos(\theta^{*}_{34})$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $cos(\theta^{*}_{34})$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $p_{T,jj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $p_{T,jj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $p_{T,4ljj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $p_{T,4ljj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Enhanced region as a function of $S_{T,4ljj}$. Overflow events are included in the last bin of the distribution.

Differential cross-sections for inclusive 4ljj production in the VBS-Suppressed region as a function of $S_{T,4ljj}$. Overflow events are included in the last bin of the distribution.

Statistical correlation between different bins of the unfolded cross-section for all the measured observables in the VBS-Enhanced region.

Statistical correlation between different bins of the unfolded cross-section for all the measured observables in the VBS-Suppressed region.

Distribution (events/100 GeV) of the reconstructed $m_{tb}$ for data and backgrounds in the 0-lepton channel's the signal region 3 after the background-only fit to data. The systematics uncertainty is shown for the post-fit background sum, including the background statistical uncertainty. The individual background components are obtained after the fit, too. There are also the pre-fit background sum and the expected signal distribution. The distribution of the $W^{\prime}$ boson signal for a mass of 3 TeV is normalised to the predicted cross-section. The last bin in each distribution includes overflow.

Distribution (events/100 GeV) of the reconstructed $m_{tb}$ for data and backgrounds in the 0-lepton channel's validation region before the fit to data. The systematics uncertainty is shown for the pre-fit background sum, including the background statistical uncertainty. The individual background components are obtained before the fit, too. There is also the expected signal distribution. The distribution of the $W^{\prime}$ boson signal for a mass of 3 TeV is normalised to the predicted cross-section. The last bin in each distribution includes overflow.

Reconstructed $m_{tb}$ distributions (events/50 GeV) for data and backgrounds in the 1-lepton channel's 2j1b signal region after the background-only fit to data. The systematics uncertainty is shown for the post-fit background sum, including the background statistical uncertainty. The individual background components are obtained after the fit, too. There are also the pre-fit background sum and the expected signal distribution. The distribution of the $W^{\prime}$ boson signal for a mass of 3 TeV is normalised to the predicted cross-section. The last bin in each distribution includes overflow.

Reconstructed $m_{tb}$ distributions (events/50 GeV) for data and backgrounds in the 1-lepton channel's 3j1b signal region after the background-only fit to data. The systematics uncertainty is shown for the post-fit background sum, including the background statistical uncertainty. The individual background components are obtained after the fit, too. There are also the pre-fit background sum and the expected signal distribution. The distribution of the $W^{\prime}$ boson signal for a mass of 3 TeV is normalised to the predicted cross-section. The last bin in each distribution includes overflow.

Reconstructed $m_{tb}$ distributions (events/50 GeV) for data and backgrounds in the 1-lepton channel's 2j2b signal region after the background-only fit to data. The systematics uncertainty is shown for the post-fit background sum, including the background statistical uncertainty. The individual background components are obtained after the fit, too. There are also the pre-fit background sum and the expected signal distribution. The distribution of the $W^{\prime}$ boson signal for a mass of 3 TeV is normalised to the predicted cross-section. The last bin in each distribution includes overflow.

Reconstructed $m_{tb}$ distributions (events/50 GeV) for data and backgrounds in the 1-lepton channel's 3j2b signal region after the background-only fit to data. The systematics uncertainty is shown for the post-fit background sum, including the background statistical uncertainty. The individual background components are obtained after the fit, too. There are also the pre-fit background sum and the expected signal distribution. The distribution of the $W^{\prime}$ boson signal for a mass of 3 TeV is normalised to the predicted cross-section. The last bin in each distribution includes overflow.

Reconstructed $m_{tb}$ distributions (events/50 GeV) for data and backgrounds in the 1-lepton channel's 2j1b control region after the background-only fit to data. The systematics uncertainty is shown for the post-fit background sum, including the background statistical uncertainty. The individual background components are obtained after the fit, too. There are also the pre-fit background sum and the expected signal distribution. The distribution of the $W^{\prime}$ boson signal for a mass of 3 TeV is normalised to the predicted cross-section. The last bin in each distribution includes overflow.

Reconstructed $m_{tb}$ distributions (events/50 GeV) for data and backgrounds in the 1-lepton channel's 3j1b control region after the background-only fit to data. The systematics uncertainty is shown for the post-fit background sum, including the background statistical uncertainty. The individual background components are obtained after the fit, too. There are also the pre-fit background sum and the expected signal distribution. The distribution of the $W^{\prime}$ boson signal for a mass of 3 TeV is normalised to the predicted cross-section. The last bin in each distribution includes overflow.

Reconstructed $m_{tb}$ distributions (events/50 GeV) for data and backgrounds in the 1-lepton channel's 2j1b validation region after the background-only fit to data. The systematics uncertainty is shown for the post-fit background sum, including the background statistical uncertainty. The individual background components are obtained after the fit, too. There are also the pre-fit background sum and the expected signal distribution. The distribution of the $W^{\prime}$ boson signal for a mass of 3 TeV is normalised to the predicted cross-section. The last bin in each distribution includes overflow.

Reconstructed $m_{tb}$ distributions (events/50 GeV) for data and backgrounds in the 1-lepton channel's 3j1b validation region after the background-only fit to data. The systematics uncertainty is shown for the post-fit background sum, including the background statistical uncertainty. The individual background components are obtained after the fit, too. There are also the pre-fit background sum and the expected signal distribution. The distribution of the $W^{\prime}$ boson signal for a mass of 3 TeV is normalised to the predicted cross-section. The last bin in each distribution includes overflow.

Observed and expected 95% CL limits on the cross-section times branching ratio for the production of a $W^{\prime}$ boson with decay to tb and left-handed couplings as a function of the mass of the $W^{\prime}$ boson and a coupling value of $g^{\prime}/g$=0.5. The observed and expected limits are derived using a linear interpolation between simulated signal mass hypotheses. The uncertainty on the theory prediction includes components from factorisation and renormalisation scales, PDFs, the strong coupling constant, and the top-quark mass.

Observed and expected 95% CL limits on the cross-section times branching ratio for the production of a $W^{\prime}$ boson with decay to tb and left-handed couplings as a function of the mass of the $W^{\prime}$ boson and a coupling value of $g^{\prime}/g$=1.0. The observed and expected limits are derived using a linear interpolation between simulated signal mass hypotheses. The uncertainty on the theory prediction includes components from factorisation and renormalisation scales, PDFs, the strong coupling constant, and the top-quark mass.

Observed and expected 95% CL limits on the cross-section times branching ratio for the production of a $W^{\prime}$ boson with decay to tb and left-handed couplings as a function of the mass of the $W^{\prime}$ boson and a coupling value of $g^{\prime}/g$=2.0. The observed and expected limits are derived using a linear interpolation between simulated signal mass hypotheses. The uncertainty on the theory prediction includes components from factorisation and renormalisation scales, PDFs, the strong coupling constant, and the top-quark mass.

Observed and expected 95% CL limits on the cross-section times branching ratio for the production of a $W^{\prime}$ boson with decay to tb and right-handed couplings as a function of the mass of the $W^{\prime}$ boson and a coupling value of $g^{\prime}/g$=0.5. The observed and expected limits are derived using a linear interpolation between simulated signal mass hypotheses. The uncertainty on the theory prediction includes components from factorisation and renormalisation scales, PDFs, the strong coupling constant, and the top-quark mass.

Observed and expected 95% CL limits on the cross-section times branching ratio for the production of a $W^{\prime}$ boson with decay to tb and right-handed couplings as a function of the mass of the $W^{\prime}$ boson and a coupling value of $g^{\prime}/g$=1.0. The observed and expected limits are derived using a linear interpolation between simulated signal mass hypotheses. The uncertainty on the theory prediction includes components from factorisation and renormalisation scales, PDFs, the strong coupling constant, and the top-quark mass.

Observed and expected 95% CL limits on the cross-section times branching ratio for the production of a $W^{\prime}$ boson with decay to tb and right-handed couplings as a function of the mass of the $W^{\prime}$ boson and a coupling value of $g^{\prime}/g$=2.0. The observed and expected limits are derived using a linear interpolation between simulated signal mass hypotheses. The uncertainty on the theory prediction includes components from factorisation and renormalisation scales, PDFs, the strong coupling constant, and the top-quark mass.

95% CL limits on the coupling value ($g^{\prime}/g$) and the $W^{\prime}$-boson mass for left-handed $W^{\prime}$-boson couplings. The area above the line is excluded.

95% CL limits on the coupling value ($g^{\prime}/g$) and the $W^{\prime}$-boson mass for right-handed $W^{\prime}$-boson couplings. The area above the line is excluded.

The product of fiducial acceptance and selection efficiency for the three signal regions of the 0-lepton channel as a function of the mass of the $W^{\prime}$ boson with left-handed chirality. The $W^{\prime}$ boson coupling strength is set to $g^{\prime}/g=1$.

The product of fiducial acceptance and selection efficiency for the three signal regions of the 0-lepton channel as a function of the mass of the $W^{\prime}$ boson with left-handed chirality. The $W^{\prime}$ boson coupling strength is set to $g^{\prime}/g=0.5$.

The product of fiducial acceptance and selection efficiency for the three signal regions of the 0-lepton channel as a function of the mass of the $W^{\prime}$ boson with left-handed chirality. The $W^{\prime}$ boson coupling strength is set to $g^{\prime}/g=2$.

The product of fiducial acceptance and selection efficiency for the four signal regions of the 1-lepton channel as a function of the mass of the $W^{\prime}$ boson with left-handed chirality. The $W^{\prime}$ boson coupling strength is set to $g^{\prime}/g=1$.

The product of fiducial acceptance and selection efficiency for the four signal regions of the 1-lepton channel as a function of the mass of the $W^{\prime}$ boson with left-handed chirality. The $W^{\prime}$ boson coupling strength is set to $g^{\prime}/g=0.5$.

The product of fiducial acceptance and selection efficiency for the four signal regions of the 1-lepton channel as a function of the mass of the $W^{\prime}$ boson with left-handed chirality. The $W^{\prime}$ boson coupling strength is set to $g^{\prime}/g=2$.

The product of fiducial acceptance and selection efficiency for the three signal regions of the 0-lepton channel as a function of the mass of the $W^{\prime}$ boson with right-handed chirality. The $W^{\prime}$ boson coupling strength is set to $g^{\prime}/g=1$.

The product of fiducial acceptance and selection efficiency for the three signal regions of the 0-lepton channel as a function of the mass of the $W^{\prime}$ boson with right-handed chirality. The $W^{\prime}$ boson coupling strength is set to $g^{\prime}/g=0.5$.

The product of fiducial acceptance and selection efficiency for the three signal regions of the 0-lepton channel as a function of the mass of the $W^{\prime}$ boson with right-handed chirality. The $W^{\prime}$ boson coupling strength is set to $g^{\prime}/g=2$.

The product of fiducial acceptance and selection efficiency for the four signal regions of the 1-lepton channel as a function of the mass of the $W^{\prime}$ boson with right-handed chirality. The $W^{\prime}$ boson coupling strength is set to $g^{\prime}/g=1$.

The product of fiducial acceptance and selection efficiency for the four signal regions of the 1-lepton channel as a function of the mass of the $W^{\prime}$ boson with right-handed chirality. The $W^{\prime}$ boson coupling strength is set to $g^{\prime}/g=0.5$.

The product of fiducial acceptance and selection efficiency for the four signal regions of the 1-lepton channel as a function of the mass of the $W^{\prime}$ boson with right-handed chirality. The $W^{\prime}$ boson coupling strength is set to $g^{\prime}/g=2$.

Observed 95% CL upper limits on the cross section for all masses and charges of photon-fusion pair-produced spin-1/2 monopoles as a function of mass for magnetic charges $|g|=1g_D$ and $|g|=2g_D$.

Observed 95% CL upper limits on the cross section for all masses and charges of Drell-Yan spin-0 HECOs production as a function of mass for various values of electric charge in the range $20\le|z|\le100$.

Observed 95% CL upper limits on the cross section for all masses and charges of Drell-Yan spin-1/2 HECOs production as a function of mass for various values of electric charge in the range $20\le|z|\le100$.

Observed 95% CL upper limits on the cross section for all masses and charges of photon-fusion pair-produced spin-0 HECOs as a function of mass for various values of electric charge in the range $20\le|z|\le100$.

Observed 95% CL upper limits on the cross section for all masses and charges of photon-fusion pair-produced spin-1/2 HECOs as a function of mass for various values of electric charge in the range $20\le|z|\le100$.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=1g_\textrm{D}$ monopoles of mass 200 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=1g_\textrm{D}$ monopoles of mass 500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=1g_\textrm{D}$ monopoles of mass 1000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=1g_\textrm{D}$ monopoles of mass 1500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=1g_\textrm{D}$ monopoles of mass 2000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=1g_\textrm{D}$ monopoles of mass 2500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=1g_\textrm{D}$ monopoles of mass 3000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=1g_\textrm{D}$ monopoles of mass 4000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=2g_\textrm{D}$ monopoles of mass 200 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=2g_\textrm{D}$ monopoles of mass 500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=2g_\textrm{D}$ monopoles of mass 1000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=2g_\textrm{D}$ monopoles of mass 1500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=2g_\textrm{D}$ monopoles of mass 2000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=2g_\textrm{D}$ monopoles of mass 2500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=2g_\textrm{D}$ monopoles of mass 3000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for $g=2g_\textrm{D}$ monopoles of mass 4000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=20$ of mass 200 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=20$ of mass 500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=20$ of mass 1000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=20$ of mass 1500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=20$ of mass 2000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=20$ of mass 2500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=20$ of mass 3000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=20$ of mass 4000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=40$ of mass 200 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=40$ of mass 500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=40$ of mass 1000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=40$ of mass 1500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=40$ of mass 2000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=40$ of mass 2500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=40$ of mass 3000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=40$ of mass 4000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=60$ of mass 200 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=60$ of mass 500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=60$ of mass 1000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=60$ of mass 1500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=60$ of mass 2000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=60$ of mass 2500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=60$ of mass 3000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=60$ of mass 4000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=80$ of mass 200 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=80$ of mass 500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=80$ of mass 1000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=80$ of mass 1500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=80$ of mass 2000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=80$ of mass 2500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=80$ of mass 3000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=80$ of mass 4000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=100$ of mass 200 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=100$ of mass 500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=100$ of mass 1000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=100$ of mass 1500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=100$ of mass 2000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=100$ of mass 2500 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=100$ of mass 3000 GeV.

Selection efficiency as a function of transverse kinetic energy $E^\text{kin}_\text{T}=E_\text{kin}\sin\theta$ and pseudorapidity $|\eta|$ for HECOs of charge $|z|=100$ of mass 4000 GeV.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Event yields after a background-only fit and data in the three signal regions. All the backgrounds except the fake $\tau_{\mathrm{had}}$ are estimated using the real $\tau_{\mathrm{had}}$ contributions from MC. The fake $\tau_{\mathrm{had}}$ background contribution is data-driven. The top-quark, $V$+jets, and small backgrounds are included in "Others". The total uncertainties include both systematic and statistical. In the $WW1\ell2\tau_{\mathrm{had}}$ channel, the smaller uncertainty in the total background than the fake $\tau_{\mathrm{had}}$ is due to the anticorrelations between the nuisance parameters of the fake $\tau_{\mathrm{had}}$ and diboson backgrounds obtained in the fit.

Uncertainties on $\sigma(pp\to X\to SH)$ obtained from the combined fit to the data for the least and most sensitive mass points. The uncertainties are symmetrised and grouped into the categories described in the systematic section.

Summary of observed and expected 95\% CL upper limits on the production cross-section of $\sigma(pp\to X\to SH)$ in fb in the combination of $WW$+$ZZ$, and in the $WW$ and $ZZ$ individual channels. The best expected limit is set to 85 fb for $m_{X}$ = 1250 GeV and $_{S}$ = 300 GeV. The combination limits gain about 26--53% over the best individual limit from the $WW1\ell2\tau_{\mathrm{had}}$ channel.

Observed and expected upper limits at 95% CL on the production cross-section times branching fraction as a function of $m_{A'}$ at dark Higgs boson mass of 50 GeV

Observed and expected upper limits at 95% CL on the production cross-section times branching fraction as a function of $m_{A'}$ at dark Higgs boson mass of 60 GeV

Observed and expected upper limits at 95% CL on the production cross-section times branching fraction as a function of $m_{A'}$ at dark Higgs boson mass of 70 GeV

Observed 90% CL upper limits on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 20 GeV

Observed 90% CL upper limits on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 30 GeV

Observed 90% CL upper limits on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 40 GeV

Observed 90% CL upper limits on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 50 GeV

Observed 90% CL upper limits on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 60 GeV

Observed 90% CL upper limits on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 70 GeV

Observed 90% CL upper limits on $\epsilon^{2}$, assuming of $\alpha_{D}=0.1$, as a function of $m_{A'}$ at dark Higgs boson mass of 20 GeV

Observed 90% CL upper limits on $\epsilon^{2}$, assuming of $\alpha_{D}=0.1$, as a function of $m_{A'}$ at dark Higgs boson mass of 30 GeV

Observed 90% CL upper limits on $\epsilon^{2}$, assuming of $\alpha_{D}=0.1$, as a function of $m_{A'}$ at dark Higgs boson mass of 40 GeV

Observed 90% CL upper limits on $\epsilon^{2}$, assuming of $\alpha_{D}=0.1$, as a function of $m_{A'}$ at dark Higgs boson mass of 50 GeV

Observed 90% CL upper limits on $\epsilon^{2}$, assuming of $\alpha_{D}=0.1$, as a function of $m_{A'}$ at dark Higgs boson mass of 60 GeV

Observed 90% CL upper limits on $\epsilon^{2}$, assuming of $\alpha_{D}=0.1$, as a function of $m_{A'}$ at dark Higgs boson mass of 70 GeV

Observed and expected upper limits at 95% CL on the branching fraction as a function of $m_{A'}$ at dark Higgs boson mass of 20 GeV

Observed and expected upper limits at 95% CL on the branching fraction as a function of $m_{A'}$ at dark Higgs boson mass of 30 GeV

Observed and expected upper limits at 95% CL on the branching fraction as a function of $m_{A'}$ at dark Higgs boson mass of 40 GeV

Observed and expected upper limits at 95% CL on the branching fraction as a function of $m_{A'}$ at dark Higgs boson mass of 50 GeV

Observed and expected upper limits at 95% CL on the branching fraction as a function of $m_{A'}$ at dark Higgs boson mass of 60 GeV

Observed and expected upper limits at 95% CL on the branching fraction as a function of $m_{A'}$ at dark Higgs boson mass of 70 GeV

Observed and expected upper limits at 95% CL on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 20 GeV

Observed and expected upper limits at 95% CL on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 30 GeV

Observed and expected upper limits at 95% CL on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 40 GeV

Observed and expected upper limits at 95% CL on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 50 GeV

Observed and expected upper limits at 95% CL on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 60 GeV

Observed and expected upper limits at 95% CL on $\alpha_{D}\epsilon^{2}$ as a function of $m_{A'}$ at dark Higgs boson mass of 70 GeV

Events in bins of the $E_{miss}^T$ in the SR.

Di-muon mass distribution for the $\mu\mu\mu$ signal region for data and expected background. The expected signal distribution for $m_a = 35$ GeV is shown assuming $\sigma(t\bar{t}a)\times \text{Br}(a\rightarrow\mu\mu) = 4$ fb. Rare background processes include $ZZ+$jets, $WWZ$, $WZZ$, $ZZZ$, and $t\bar{t}t\bar{t}$.

Fit of the Double Crystal Ball width as a function of the $a$-boson mass in the $e\mu\mu$ channel.

Expected and observed 95% CL upper limit on signal cross section as a function of $m_a$ for $t\bar{t}a$.

Expected and observed 95% CL upper limit on signal cross section as a function of $m_a$ for $H^{+}\rightarrow aW$, considering a top pair cross section of $\sigma(pp\rightarrow t\bar{t})= 833$ pb for the charged Higgs boson mass, $m_{H^{+}} = 120$ GeV.

Expected and observed 95% CL upper limit on signal cross section as a function of $m_a$ for $H^{+}\rightarrow aW$, considering a top pair cross section of $\sigma(pp\rightarrow t\bar{t})= 833$ pb for the charged Higgs boson mass, $m_{H^{+}} = 140$ GeV.

Expected and observed 95% CL upper limit on signal cross section as a function of $m_a$ for $H^{+}\rightarrow aW$, considering a top pair cross section of $\sigma(pp\rightarrow t\bar{t})= 833$ pb for the charged Higgs boson mass, $m_{H^{+}} = 160$ GeV.

Comparison of the expected background with data in a region dominated by non-prompt background. The region is defined exactly like the $e\mu\mu$ signal region but requiring the two selected muons to have the same electrical charge.

Scan of the observed p-value as a function of $m_a$ for the background-only hypothesis.

Scan of the observed p-value as a function of $m_a$ for the background-only hypothesis for $e\mu\mu$ channel.

Scan of the observed p-value as a function of $m_a$ for the background-only hypothesis for $\mu\mu\mu$ channel.

Acceptance times efficiency for the $t\bar{t}a$, $a\rightarrow \mu\mu$ signal as a function of the $m_a$ hypothesis.

Acceptance times efficiency for the $t \rightarrow bH^+$, $H^+ \rightarrow W^+ a$, $a\rightarrow \mu\mu$ signal as a function of the $m_a$ and $m_{H^+}$ hypotheses, for the $e\mu\mu$ signal region. The upper left region of the plot corresponds to signals for which on-shell decays of $H^+ \rightarrow Wa$ are not possible. The CP conjugate process is also included.

Acceptance times efficiency for the $t \rightarrow bH^+ $, $H^+ \rightarrow W^+ a$, $a\rightarrow \mu\mu$ signal as a function of the $m_a$ and $m_{H^+}$ hypotheses, for the $\mu\mu\mu$ signal region. The upper left region of the plot corresponds to signals for which on-shell decays of $H^+ \rightarrow Wa$ are not possible. The CP conjugate process is also included.

Expected lower limits on $\tan\beta$ for the $t\rightarrow H^{+}b, H^{+}\rightarrow aW, a\rightarrow\mu\mu$ signal in the context of the 2HDM as a function of $m_{H^{+}}$ and $m_a$.

Observed lower limits on $\tan\beta$ for the $t\rightarrow H^{+}b, H^{+}\rightarrow aW, a\rightarrow\mu\mu$ signal in the context of the 2HDM as a function of $m_{H^{+}}$ and $m_a$.

Cutflow for two signal mass points for $t\bar{t}a$, $m_{a} = 20$ GeV and $m_{a} = 60$ GeV, and one signal mass point for $H^+ \rightarrow Wa$ for the electron channel $e\mu\mu$. For the signal yields, a cross section times branching ratio of 1 fb is assumed. The yields are presented before the profile likelihood fit.

Cutflow for the dominant backgrounds estimated from simulation ($t\bar{t}Z$, $WZ$, $t\bar{t}H$) and data, for the electron channel $e\mu\mu$. The yields are presented before the profile likelihood fit.

Cutflow for the signal mass point $t\bar{t}a$, $m_{a} = 20$ GeV, as well as the dominant backgrounds estimated from simulation ($t\bar{t}Z$, $WZ$, $t\bar{t}H$) and data, for the corresponding signal mass hypothesis, for the muon channel $\mu\mu\mu$. For the signal yields, a cross section times branching ratio of 1 fb is assumed. The yields are presented before the profile likelihood fit.

Cutflow for the signal mass point $t\bar{t}a$, $m_{a} = 60$ GeV, as well as the dominant backgrounds estimated from simulation ($t\bar{t}Z$, $WZ$, $t\bar{t}H$) and data, for the corresponding signal mass hypothesis, for the muon channel $\mu\mu\mu$. For the signal yields, a cross section times branching ratio of 1 fb is assumed. The yields are presented before the profile likelihood fit.

Cutflow for $H^+ \rightarrow W a$ with $m_{H^{+}} = 120$ GeV, $m_{a} = 30$ GeV for the muon channel $\mu\mu\mu$. A cross section times branching ratio of 1 fb is assumed. The yields are presented before the profile likelihood fit.

Total yields in the $e\mu\mu$ and $\mu\mu\mu$ signal regions for the expected background and data after the profile likelihood fit for a $m_a = 30$ GeV signal mass hypothesis.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Acceptance $\times$ efficiency for bSR1e.

Acceptance $\times$ efficiency for bSR2e.

Acceptance $\times$ efficiency for bSR2mu.

Acceptance $\times$ efficiency for rSROS2e.

Acceptance $\times$ efficiency for rSROS2mu.

Acceptance $\times$ efficiency for rSRSS2e.

Acceptance $\times$ efficiency for rSRSS2mu.

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

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

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

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

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

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

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

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

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

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

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

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

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

Signal selection efficiency ($\epsilon$) times acceptance ($A$) as a function of $H^{\pm}$. The estimated $\epsilon\times A$ arises from the lepton selection and triggering (∼30%) as well as jet selection and flavour tagging (∼10% or lower). The decrease of $\epsilon\times A$ for $m_{H^{\pm}}$ = 120 GeV and higher is expected from the kinematic constraint on the $H^{\pm}$ decay products due to the top-quark mass.

Fit results of the simultaneous measurements of the $H\to e\tau$ and $H\to \mu\tau$ signals (2POI) showing best-fit values of the LFV branching ratios of the Higgs boson $\hat{B}$($H\to e\tau$). The results from standalone channel/categories fits are compared with the results of the combined fit.

Fit results of the simultaneous measurements of the $H\to e\tau$ and $H\to \mu\tau$ signals (2POI) showing upper limits at 95% C.L. on the LFV branching ratios of the Higgs boson $H\to \mu\tau$. The results from standalone channel/categories fits are compared with the results of the combined fit.

Fit results of the simultaneous measurements of the $H\to e\tau$ and $H\to \mu\tau$ signals (2POI) showing upper limits at 95% C.L. on the LFV branching ratios of the Higgs boson $H\to \mu\tau$. The results from standalone channel/categories fits are compared with the results of the combined fit.

Fit results of the simultaneous measurements of the $H\to e\tau$ and $H\to \mu\tau$ signals (2POI) showing best-fit values of the LFV branching ratios of the Higgs boson $\hat{B}$($H\to \mu\tau$). The results from standalone channel/categories fits are compared with the results of the combined fit.

Fit results of the simultaneous measurements of the $H\to e\tau$ and $H\to \mu\tau$ signals (2POI) showing best-fit values of the LFV branching ratios of the Higgs boson $\hat{B}$($H\to \mu\tau$). The results from standalone channel/categories fits are compared with the results of the combined fit.

Fit results of the independent searches (1 POI) showing upper limits at 95% C.L. on the LFV branching ratios of the Higgs boson $H\to e\tau$. The results from standalone channel/categories fits are compared with the results of the combined fit.

Fit results of the independent searches (1 POI) showing upper limits at 95% C.L. on the LFV branching ratios of the Higgs boson $H\to e\tau$. The results from standalone channel/categories fits are compared with the results of the combined fit.

Fit results of the independent searches (1 POI) showing best-fit values of the LFV branching ratios of the Higgs boson $\hat{B}$($H\to e\tau$). The results from standalone channel/categories fits are compared with the results of the combined fit.

Fit results of the independent searches (1 POI) showing best-fit values of the LFV branching ratios of the Higgs boson $\hat{B}$($H\to e\tau$). The results from standalone channel/categories fits are compared with the results of the combined fit.

Fit results of the independent searches (1 POI) showing upper limits at 95% C.L. on the LFV branching ratios of the Higgs boson $H\to \mu\tau$. The results from standalone channel/categories fits are compared with the results of the combined fit.

Fit results of the independent searches (1 POI) showing upper limits at 95% C.L. on the LFV branching ratios of the Higgs boson $H\to \mu\tau$. The results from standalone channel/categories fits are compared with the results of the combined fit.

Fit results of the independent searches (1 POI) showing best-fit values of the LFV branching ratios of the Higgs boson $\hat{B}$($H\to \mu\tau$). The results from standalone channel/categories fits are compared with the results of the combined fit.

Fit results of the independent searches (1 POI) showing best-fit values of the LFV branching ratios of the Higgs boson $\hat{B}$($H\to \mu\tau$). The results from standalone channel/categories fits are compared with the results of the combined fit.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{\ell'}$ channel for $e\tau_{\mu}$ final state in the non-VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{\ell'}$ channel for $e\tau_{\mu}$ final state in the non-VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{\ell'}$ channel for $\mu\tau_{e}$ final state in the non-VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{\ell'}$ channel for $\mu\tau_{e}$ final state in the non-VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{\ell'}$ channel for $e\tau_{\mu}$ final state in the VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{\ell'}$ channel for $e\tau_{\mu}$ final state in the VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{\ell'}$ channel for $\mu\tau_{e}$ final state in the VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{\ell'}$ channel for $\mu\tau_{e}$ final state in the VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{had}$ channel for $e\tau_{had}$ final state in the non-VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{had}$ channel for $e\tau_{had}$ final state in the non-VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{had}$ channel for $\mu\tau_{had}$ final state in the non-VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{had}$ channel for $\mu\tau_{had}$ final state in the non-VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{had}$ channel for $e\tau_{had}$ final state in the VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{had}$ channel for $e\tau_{had}$ final state in the VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{had}$ channel for $\mu\tau_{had}$ final state in the VBF category.

BDT score distribution, after the simultaneous fit of the $H\to e\tau$ and $H\to \mu\tau$ signals, obtained by fitting the data of the MC-template $\ell\tau_{had}$ channel for $\mu\tau_{had}$ final state in the VBF category.

NN score distribution, after an independent fit of the $H\to \mu\tau$ signal, obtained by fitting the data of the Symmetry $\ell\tau_{\ell'}$ channel, for the $\mu\tau_{e}$ final state in the non-VBF category.

NN score distribution, after an independent fit of the $H\to \mu\tau$ signal, obtained by fitting the data of the Symmetry $\ell\tau_{\ell'}$ channel, for the $\mu\tau_{e}$ final state in the non-VBF category.

NN score distribution, after an independent fit of the $H\to \mu\tau$ signal, obtained by fitting the data of the Symmetry $\ell\tau_{\ell'}$ channel, for the $\mu\tau_{e}$ final state in the VBF category.

NN score distribution, after an independent fit of the $H\to \mu\tau$ signal, obtained by fitting the data of the Symmetry $\ell\tau_{\ell'}$ channel, for the $\mu\tau_{e}$ final state in the VBF category.

Best-fit $\mathcal{B}(H\to \mu\tau)- \mathcal{B}(H\to e\tau)$ values, given in %, obtained from the standalone fits of the $\ell\tau_\ell'$ final state with either background estimation method. For the MC-template method, the correlated uncertainties are fixed to their best-fit values and the uncertainties are re-derived considering only the uncorrelated sources. For the Symmetry method, the full uncertainty is considered.

Best-fit $\mathcal{B}(H\to \mu\tau)- \mathcal{B}(H\to e\tau)$ values, given in %, obtained from the standalone fits of the $\ell\tau_\ell'$ final state with either background estimation method. For the MC-template method, the correlated uncertainties are fixed to their best-fit values and the uncertainties are re-derived considering only the uncorrelated sources. For the Symmetry method, the full uncertainty is considered.

Best-fit value of the branching ratios $\mathcal{B}(H\to e\tau)$ and $\mathcal{B}(H\to \mu\tau)$, given in %, and likelihood contours at 68% and 95% C.L. Obtained from the simultaneous fit of $H\to e\tau$ and $H\to \mu\tau$ signals based on the MC-template method, compared to the SM expectation.

Best-fit value of the branching ratios $\mathcal{B}(H\to e\tau)$ and $\mathcal{B}(H\to \mu\tau)$, given in %, and likelihood contours at 68% and 95% C.L. Obtained from the simultaneous fit of $H\to e\tau$ and $H\to \mu\tau$ signals based on the MC-template method, compared to the SM expectation.

The performance of jet energy resolution (JER) for jets with |y| < 2.1 evaluated as a function of pT_truth in different centrality bins. Simulated hard scatter events were overlaid onto events from a dedicated sample of minimum-bias Xe+Xe data. The fit parameters are listed in a sperate table (Extras 1)

The relative magnitude of systematic uncertainties for per-pair normalized xJ distributions in 0-10% Xe+Xe centrality

The relative magnitude of systematic uncertainties for absolutely normalized xJ distributions in 0-10% Xe+Xe centrality

The relative magnitude of systematic uncertainties for rho distributions for leading jets in 0-10% Xe+Xe centrality

Per-pair normalized xJ distribution evaluated in four centrality intervals and given pT1 interval.

Per-pair normalized xJ distribution evaluated in four centrality intervals and given pT1 interval.

Per-pair normalized xJ distribution evaluated in four centrality intervals and given pT1 interval.

Absolutely normalized xJ distribution evaluated in four centrality intervals and given pT1 interval.

Absolutely normalized xJ distribution evaluated in four centrality intervals and given pT1 interval.

Absolutely normalized xJ distribution evaluated in four centrality intervals and given pT1 interval.

Per-pair normalized xJ distribution in Xe+Xe collisions.

Per-pair normalized xJ distribution in Pb+Pb collisions.

Per-pair normalized xJ distribution in Xe+Xe collisions.

Per-pair normalized xJ distribution in Pb+Pb collisions.

Per-pair normalized xJ distribution in Xe+Xe collisions.

Per-pair normalized xJ distribution in Pb+Pb collisions.

Absolutely normalized xJ distribution in Xe+Xe collisions.

Absolutely normalized xJ distribution in Xe+Xe collisions.

Absolutely normalized xJ distribution in Xe+Xe collisions.

Absolutely normalized xJ distribution in Xe+Xe collisions.

Absolutely normalized xJ distribution in Xe+Xe collisions.

Absolutely normalized xJ distribution in Xe+Xe collisions.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of leading jet pT in the same centrality intervals.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of leading jet pT in the same centrality intervals.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of leading jet pT in the same centrality intervals.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of subleading jet pT in the same centrality intervals.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of subleading jet pT in the same centrality intervals.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of subleading jet pT in the same centrality intervals.

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of leading jet pT in the same SUM ETFCal intervals (selecting equivalent event activity)

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of leading jet pT in the same SUM ETFCal intervals (selecting equivalent event activity)

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of leading jet pT in the same SUM ETFCal intervals (selecting equivalent event activity)

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of subleading jet pT in the same SUM ETFCal intervals (selecting equivalent event activity)

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of subleading jet pT in the same SUM ETFCal intervals (selecting equivalent event activity)

The ratios of Xe+Xe and Pb+Pb pair nuclear-modification factors, rho, evaluated as a function of subleading jet pT in the same SUM ETFCal intervals (selecting equivalent event activity)

Parameter a,b, and c from JER fits in Figure 1b.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Per-pair normalized xJ distribution in Xe+Xe collisions for selected Pb+Pb centrality and pT1 bin.

Distribution of the NN output in the $\geq$1fj$\,$SR in data and the expected contribution of the signal and background processes after the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions considering the correlations of the uncertainties as obtained by the fit.

Distribution of the NN output in the $t\bar{t}\gamma\,$CR in data and the expected contribution of the signal and background processes after the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions considering the correlations of the uncertainties as obtained by the fit.

Total event yield in the $W\gamma\,$CR in data and the expected contribution of the signal and background processes after the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions considering the correlations of the uncertainties as obtained by the fit.

Distribution of the scalar sum of the jet transverse momenta in the 0fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions. The first and last bins include the underflow and overflow, respectively.

Distribution of the $\eta$ of the $b$-tagged jet in the 0fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions.

Distribution of the reconstructed top-quark mass in the 0fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions. The first and last bins include the underflow and overflow, respectively.

Distribution of the $p_{\mathrm{T}}$ of the top-quark+photon system in the 0fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions and the last bin includes the overflow.

Distribution of the photon $p_{\mathrm{T}}$ in the 0fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions and the last bin includes the overflow.

Distribution of the photon $\eta$ in the 0fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions.

Distribution of the scalar sum of the jet transverse momenta in the $\geq$1fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions and the last bin includes the overflow.

Distribution of the invariant mass of the $b$-tagged jet and the highest-$p_{\mathrm{T}}$ forward jet in the $\geq$1fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions and the last bin includes the overflow.

Distribution of $p_{\mathrm{T}}$ of the highest-$p_{\mathrm{T}}$ forward jet in the $\geq$1fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions and the last bin includes the overflow.

Distribution of the difference in $\eta$ between the highest-$p_{\mathrm{T}}$ forward jet and the photon in the $\geq$1fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions and the last bin includes the overflow.

Distribution of the energy of the system formed by the highest-$p_{\mathrm{T}}$ forward jet and the photon in the $\geq$1fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions. The first and last bins include the underflow and overflow, respectively.

Distribution of the $\eta$ of the $b$-tagged jet in the $\geq$1fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions.

Distribution of the reconstructed top-quark mass in the $\geq$1fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions. The first and last bins include the underflow and overflow, respectively.

Distribution of the $p_{\mathrm{T}}$ of the top-quark and photon system in the $\geq$1fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions and the last bin includes the overflow.

Distribution of the photon $p_{\mathrm{T}}$ in the $\geq$1fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions and the last bin includes the overflow.

Distribution of the photon $\eta$ in the $\geq$1fj$\,$SR in data and for the sum of all processes expectations before the profile-likelihood fit. The "Total" column corresponds to the sum of the expected contributions from the signal and background processes. The uncertainty represents the sum of statistical and systematic uncertainties in the signal and background predictions.

Ordered list of the 30 systematic uncertainties with the largest impact on the measured signal normalisation in the fit to data in the parton-level measurement considered as nuisance parameters (NPs) in the profile-likelihood fit. The column "NP value, error" corresponds to the nominal best-fit values and the corresponding uncertainties. The impact of each NP, $\Delta\sigma$/$\sigma_{\mathrm{pred}}$, is computed by comparing the nominal best-fit value of the POI ($\sigma$/$\sigma_{\mathrm{pred}}$) with the result of the fits when fixing the considered NP to its best-fit value shifted by its pre-fit and post-fit uncertainties. The corresponding impacts are listed in the "POI impact prefit high/low" and "POI impact high/low" columns, respectively. The "MC stat." NPs represent the MC statistical uncertainty and they enter the likelihood with a Poisson term, while all the other NPs enter the likelihood via a Gaussian term.

Ordered list of the 30 systematic uncertainties with the largest impact on the measured signal normalisation in the fit to data in the particle-level measurement considered as nuisance parameters (NPs) in the profile-likelihood fit. The column "NP value, error" corresponds to the nominal best-fit values and the corresponding uncertainties. The impact of each NP, $\Delta\sigma$/$\sigma_{\mathrm{pred}}$, is computed by comparing the nominal best-fit value of the POI ($\sigma$/$\sigma_{\mathrm{pred}}$) with the result of the fits when fixing the considered NP to its best-fit value shifted by its pre-fit and post-fit uncertainties. The corresponding impacts are listed in the "POI impact prefit high/low" and "POI impact high/low" columns, respectively. The "MC stat." NPs represent the MC statistical uncertainty and they enter the likelihood with a Poisson term, while all the other NPs enter the likelihood via a Gaussian term.

This table lists the kinematic requirements on parton-level objects used to define of the fiducial phase space for the parton-level measurement. Frixione isolation ($\href{https://arxiv.org/abs/hep-ph/9801442}{\text{hep-ph/9801442}}$) with a chosen radius of $\Delta R = 0.2$ is applied to photons ($\gamma$). The measured fiducial parton-level cross section is $\sigma_{tq\gamma}\times\mathcal{B}(t\rightarrow l\nu b) = 688\pm 23(\text{stat.})^{+75}_{-71}(\text{syst.})\,$fb.

This table lists the kinematic requirements on particle-level objects used to define of the fiducial phase space for the particle-level measurement. The particle level objects are photons ($\gamma$) not from a hadron decay, neutrinos not from a hadron decay ($\nu$), prompt electrons and muons ($\ell$) "dressed" by adding close-by ($\Delta R < 0.1$) photons, and anti-$k_t$ $R = 0.4$ jets built from stable particles ($\tau > 30\,$ps) and tau leptons excluding neutrinos and prompt dressed muons. Jets are $b$-tagged ($b$-jet) using ghost-matched $b$-hadrons with $p_{\text{T}} > 5\,$GeV. Apart from the kinematic requirements, isolation and overlap removal criteria are applied. Jets within $\Delta R = 0.4$ of a photon are removed if the $p_{\text{T}}$ of charged particles within $\Delta R = 0.3$ of the photon is smaller than $10\,\%$ of its $p_{\text{T}}$. Jets within $\Delta R = 0.4$ of a lepton are removed. Events are removed where a photon is close ($\Delta R < 0.4$) to a lepton or a surviving jet. The measured fiducial particle-level cross section is $\sigma_{tq\gamma}\times\mathcal{B}(t\rightarrow l\nu b)+\sigma_{t(\rightarrow l\nu b\gamma)q} = 303\pm 9(\text{stat.})^{+33}_{-32}(\text{syst.})\,$fb.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.56<|\eta^{\gamma}|<1.81$ and photon isolation cone radius $R=0.4$.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.81<|\eta^{\gamma}|<2.01$ and photon isolation cone radius $R=0.4$.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $2.01<|\eta^{\gamma}|<2.37$ and photon isolation cone radius $R=0.4$.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $|\eta^{\gamma}|<0.6$ and photon isolation cone radius $R=0.2$.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $0.6<|\eta^{\gamma}|<0.8$ and photon isolation cone radius $R=0.2$.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $0.8<|\eta^{\gamma}|<1.37$ and photon isolation cone radius $R=0.2$.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.56<|\eta^{\gamma}|<1.81$ and photon isolation cone radius $R=0.2$.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.81<|\eta^{\gamma}|<2.01$ and photon isolation cone radius $R=0.2$.

Measured cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $2.01<|\eta^{\gamma}|<2.37$ and photon isolation cone radius $R=0.2$.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $|\eta^{\gamma}|<0.6$ and isolation cone radius $0.4$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $0.6<|\eta^{\gamma}|<0.8$ and isolation cone radius $0.4$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $0.8<|\eta^{\gamma}|<1.37$ and isolation cone radius $0.4$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.56<|\eta^{\gamma}|<1.81$ and isolation cone radius $0.4$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.81<|\eta^{\gamma}|<2.01$ and isolation cone radius $0.4$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $2.01<|\eta^{\gamma}|<2.37$ and isolation cone radius $0.4$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $|\eta^{\gamma}|<0.6$ and isolation cone radius $0.2$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $0.6<|\eta^{\gamma}|<0.8$ and isolation cone radius $0.2$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $0.8<|\eta^{\gamma}|<1.37$ and isolation cone radius $0.2$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.56<|\eta^{\gamma}|<1.81$ and isolation cone radius $0.2$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.81<|\eta^{\gamma}|<2.01$ and isolation cone radius $0.2$ at NNLO QCD.

Predicted cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $2.01<|\eta^{\gamma}|<2.37$ and isolation cone radius $0.2$ at NNLO QCD.

Measured ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $|\eta^{\gamma}|<0.6$.

Measured ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $0.6<|\eta^{\gamma}|<0.8$.

Measured ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $0.8<|\eta^{\gamma}|<1.37$.

Measured ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $1.56<|\eta^{\gamma}|<1.81$.

Measured ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $1.81<|\eta^{\gamma}|<2.01$.

Measured ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $2.01<|\eta^{\gamma}|<2.37$.

Predicted ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $|\eta^{\gamma}|<0.6$ at NNLO QCD.

Predicted ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $0.6<|\eta^{\gamma}|<0.8$ at NNLO QCD.

Predicted ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $0.8<|\eta^{\gamma}|<1.37$ at NNLO QCD.

Predicted ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $1.56<|\eta^{\gamma}|<1.81$ at NNLO QCD.

Predicted ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $1.81<|\eta^{\gamma}|<2.01$ at NNLO QCD.

Predicted ratio of the differential cross sections for inclusive isolated-photon production for $R=0.2$ and $R=0.4$ as a function of $E_{\rm T}^{\gamma}$ for $2.01<|\eta^{\gamma}|<2.37$ at NNLO QCD.

Measured fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $|\eta^{\gamma}|<0.6$.

Measured fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $0.6<|\eta^{\gamma}|<0.8$.

Measured fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $0.8<|\eta^{\gamma}|<1.37$.

Measured fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $1.56<|\eta^{\gamma}|<1.81$.

Measured fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $1.81<|\eta^{\gamma}|<2.01$.

Measured fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $2.01<|\eta^{\gamma}|<2.37$.

Predicted fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $|\eta^{\gamma}|<0.6$ at NNLO QCD.

Predicted fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $0.6<|\eta^{\gamma}|<0.8$ at NNLO QCD.

Predicted fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $0.8<|\eta^{\gamma}|<1.37$ at NNLO QCD.

Predicted fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $1.56<|\eta^{\gamma}|<1.81$ at NNLO QCD.

Predicted fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $1.81<|\eta^{\gamma}|<2.01$ at NNLO QCD.

Predicted fiducial integrated cross section for inclusive isolated-photon production as a function of $R$ for $2.01<|\eta^{\gamma}|<2.37$ at NNLO QCD.

Two-dimensional distribution of the invariant mass $m_{DV}$ and the track multiplicity in the High-pT jet SR for observed data events

Two-dimensional distribution of the invariant mass $m_{DV}$ and the track multiplicity in the High-pT jet SR for expected signal events in the strong gluino pair pair production model with m(gluino)=1.8 TeV, m(chi0)=0.2 TeV, tau(chi0)=0.1 ns

Two-dimensional distribution of the invariant mass $m_{DV}$ and the track multiplicity in the Trackless jet SR for observed data events

Two-dimensional distribution of the invariant mass $m_{DV}$ and the track multiplicity in the Trackless jet SR for expected signal events in the electroweak pair production model

Expected exclusion limits at 95% CL on the lifetime and mass of the neutralino in electroweakino pair production models

Expected (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in electroweakino pair production models

Expected (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in electroweakino pair production models

Observed exclusion limits at 95% CL on the lifetime and mass of the neutralino in electroweakino pair production models

Observed (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in electroweakino pair production models

Observed (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in electroweakino pair production models

Expected exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.4 TeV

Expected (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.4 TeV

Expected (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.4 TeV

Observed exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.4 TeV

Observed (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.4 TeV

Observed (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.4 TeV

Exclusion limits at 95% CL on the production cross section in the electroweak pair production model.

Exclusion limits at 95% CL on the production cross section in the strong gluino pair production models and m(gluino)=2.4 TeV

Expected exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.0 TeV

Expected (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.0 TeV

Expected (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.0 TeV

Observed exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.0 TeV

Observed (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.0 TeV

Observed (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.0 TeV

Expected exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.2 TeV

Expected (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.2 TeV

Expected (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.2 TeV

Observed exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.2 TeV

Observed (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.2 TeV

Observed (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the neutralino in strong gluino pair production models and m(gluino)=2.2 TeV

Expected exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=50 GeV

Expected (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=50 GeV

Expected (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=50 GeV

Observed exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=50 GeV

Observed (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=50 GeV

Observed (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=50 GeV

Expected exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=450 GeV

Expected (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=450 GeV

Expected (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=450 GeV

Observed exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=450 GeV

Observed (+1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=450 GeV

Observed (-1 sigma) exclusion limits at 95% CL on the lifetime and mass of the gluino in strong gluino pair production models and m(chi0)=450 GeV

Expected exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.01 ns

Expected (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.01 ns

Expected (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.01 ns

Observed exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.01 ns

Observed (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.01 ns

Observed (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.01 ns

Expected exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.1 ns

Expected (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.1 ns

Expected (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.1 ns

Observed exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.1 ns

Observed (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.1 ns

Observed (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=0.1 ns

Expected exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=1 ns

Expected (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=1 ns

Expected (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=1 ns

Observed exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=1 ns

Observed (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=1 ns

Observed (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=1 ns

Expected exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=10 ns

Expected (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=10 ns

Expected (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=10 ns

Observed exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=10 ns

Observed (+1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=10 ns

Observed (-1 sigma) exclusion limits at 95% CL on the mass of the gluino and neutralino in strong gluino pair production models and tau(chi0)=10 ns

Exclusion limits at 95% CL on the production cross section in the strong gluino pair production models and m($ ilde{\chi}^0_1$)=1.25 TeV

Acceptance cutflow for the High-pT SR for representative points in the strong gluino pair production model. See additional resources for more information.

Acceptance cutflow for the Trackless SR for representative points in the electroweak pair production model. See additional resources for more information.

Acceptance cutflow for the Trackless SR for representative points in the electroweak pair production model with heavy-flavor quarks final state. See additional resources for more information.

Acceptance cutflow for the High-pT SR for representative points in the electroweak pair production model with heavy-flavor quarks final state. See additional resources for more information.

Reinterpretation Material: Event-level Efficiency for HighPt SR selections, R < 1150 mm

Reinterpretation Material: Event-level Efficiency for HighPt SR selections, R [1150, 3870] mm

Reinterpretation Material: Event-level Efficiency for HighPt SR selections, R > 3870 mm

Reinterpretation Material: Event-level Efficiency for Trackless SR selections, R < 1150 mm

Reinterpretation Material: Event-level Efficiency for Trackless SR selections, R [1150, 3870] mm

Reinterpretation Material: Event-level Efficiency for Trackless SR selections, R > 3870 mm

Reinterpretation Material: Vertex-level Efficiency for R < 22 mm

Reinterpretation Material: Vertex-level Efficiency for R [22, 25] mm

Reinterpretation Material: Vertex-level Efficiency for R [25, 29] mm

Reinterpretation Material: Vertex-level Efficiency for R [29, 38] mm

Reinterpretation Material: Vertex-level Efficiency for R [38, 46] mm

Reinterpretation Material: Vertex-level Efficiency for R [46, 73] mm

Reinterpretation Material: Vertex-level Efficiency for R [73, 84] mm

Reinterpretation Material: Vertex-level Efficiency for R [84, 111] mm

Reinterpretation Material: Vertex-level Efficiency for R [111, 120] mm

Reinterpretation Material: Vertex-level Efficiency for R [120, 145] mm

Reinterpretation Material: Vertex-level Efficiency for R [145, 180] mm

Reinterpretation Material: Vertex-level Efficiency for R [180, 300] mm

Cutflow (acceptance x efficiency) for the High-pT SR for representative points in the strong gluino pair production model. See additional resources for more information.

Cutflow (acceptance x efficiency) for the Trackless SR for representative points in the electroweak pair production model. See additional resources for more information.

Cutflow (acceptance x efficiency) for the Trackless SR for representative points in the electroweak pair production model with heavy-flavor quarks. See additional resources for more information.

Cutflow (acceptance x efficiency) for the High-pT SR for representative points in the electroweak pair production model with heavy-flavor quarks. See additional resources for more information.