Two-dimensional exclusion limits for Dirac HNL pure electron coupling scenario

Two-dimensional exclusion limits for Majorana HNL pure muon coupling scenario

Two-dimensional exclusion limits for Dirac HNL pure muon coupling scenario

Two-dimensional exclusion limits for Majorana HNL pure tau lepton coupling scenario

Two-dimensional exclusion limits for Dirac HNL pure tau lepton coupling scenario

Two-dimensional exclusion limits for Majorana HNL mixed electron-muon coupling scenario

DTwo-dimensional exclusion limits for Dirac HNL mixed electron-muon coupling scenario

Two-dimensional exclusion limits for Majorana HNL mixed electron-tau lepton coupling scenario

Two-dimensional exclusion limits for Dirac HNL mixed electron-tau lepton coupling scenario

Two-dimensional exclusion limits for Majorana HNL mixed muon-tau lepton coupling scenario

Two-dimensional exclusion limits for Dirac HNL mixed muon-tau lepton coupling scenario

Two-dimensional exclusion limits for Majorana HNL democratic lepton coupling scenario

Two-dimensional exclusion limits for Dirac HNL democratic lepton coupling scenario

Lower limits on Majorana HNL mass for fixed proper decay length of $c\tau=0.1\,\text{mm}$

Lower limits on Dirac HNL mass for fixed proper decay length of $c\tau=0.1\,\text{mm}$

Lower limits on Majorana HNL mass for fixed proper decay length of $c\tau=1\,\text{mm}$

Lower limits on Dirac HNL mass for fixed proper decay length of $c\tau=1\,\text{mm}$

Lower limits on Majorana HNL proper decay length for fixed mass of 4.5 GeV

Lower limits on Dirac HNL proper decay length for fixed mass of 4.5 GeV

Lower limits on Majorana HNL proper decay length for fixed mass of 8 GeV

Lower limits on Dirac HNL proper decay length for fixed mass of 8 GeV

Isoscalar WIMP-nucleon elastic coupling limit for Operator 9

Isoscalar WIMP-nucleon elastic coupling limit for Operator 12

Isoscalar WIMP-nucleon elastic coupling limit for Operator 7

Isoscalar WIMP-nucleon elastic coupling limit for Operator 13

Isoscalar WIMP-nucleon elastic coupling limit for Operator 3

Isoscalar WIMP-nucleon elastic coupling limit for Operator 5

Isoscalar WIMP-nucleon elastic coupling limit for Operator 11

Isoscalar WIMP-nucleon elastic coupling limit for Operator 14

Isoscalar WIMP-nucleon elastic coupling limit for Operator 4

Isoscalar WIMP-nucleon elastic coupling limit for Operator 10

Isoscalar WIMP-nucleon elastic coupling limit for Operator 1

Isoscalar WIMP-nucleon elastic coupling limit for Operator 6

Isoscalar WIMP-nucleon elastic coupling limit for Operator 15

Isovector WIMP-nucleon elastic coupling limit for Operator 3

Isovector WIMP-nucleon elastic coupling limit for Operator 10

Isovector WIMP-nucleon elastic coupling limit for Operator 13

Isovector WIMP-nucleon elastic coupling limit for Operator 12

Isovector WIMP-nucleon elastic coupling limit for Operator 8

Isovector WIMP-nucleon elastic coupling limit for Operator 1

Isovector WIMP-nucleon elastic coupling limit for Operator 4

Isovector WIMP-nucleon elastic coupling limit for Operator 5

Isovector WIMP-nucleon elastic coupling limit for Operator 7

Isovector WIMP-nucleon elastic coupling limit for Operator 11

Isovector WIMP-nucleon elastic coupling limit for Operator 15

Isovector WIMP-nucleon elastic coupling limit for Operator 9

Isovector WIMP-nucleon elastic coupling limit for Operator 14

Isovector WIMP-nucleon elastic coupling limit for Operator 6

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 12

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 3

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 9

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 7

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 10

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 15

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 8

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 4

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 5

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 6

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 13

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 11

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 14

Isoscalar WIMP-nucleon inelastic coupling limit for Operator 1

Isovector WIMP-nucleon inelastic coupling limit for Operator 13

Isovector WIMP-nucleon inelastic coupling limit for Operator 1

Isovector WIMP-nucleon inelastic coupling limit for Operator 4

Isovector WIMP-nucleon inelastic coupling limit for Operator 7

Isovector WIMP-nucleon inelastic coupling limit for Operator 9

Isovector WIMP-nucleon inelastic coupling limit for Operator 3

Isovector WIMP-nucleon inelastic coupling limit for Operator 11

Isovector WIMP-nucleon inelastic coupling limit for Operator 5

Isovector WIMP-nucleon inelastic coupling limit for Operator 10

Isovector WIMP-nucleon inelastic coupling limit for Operator 8

Isovector WIMP-nucleon inelastic coupling limit for Operator 14

Isovector WIMP-nucleon inelastic coupling limit for Operator 15

Isovector WIMP-nucleon inelastic coupling limit for Operator 6

Isovector WIMP-nucleon inelastic coupling limit for Operator 12

The ratios between the hydrodynamic probes and their Taylor expansions. The ratios are plotted as functions of centrality in PbPb collisions at $\sqrt{s_{NN}}=5.02~\textrm{Te}\textrm{V}$. The measurement is performed with charged particles within the acceptance region.

The ratios between the hydrodynamic probes and their Taylor expansions. The ratios are plotted as functions of centrality in PbPb collisions at $\sqrt{s_{NN}}=5.02~\textrm{Te}\textrm{V}$. The measurement is performed with charged particles within the acceptance region.

The magnitudes of the measured skewness $\gamma_1^\textrm{exp}$ and "corrected" skewness $\gamma_{1,\textrm{corr}}^\textrm{exp}$ as functions of centrality in PbPb collisions at $\sqrt{s_{NN}}=5.02~\textrm{Te}\textrm{V}$. The measurement is performed with charged particles within the acceptance region.

The magnitudes of the measured kurtosis $\gamma_2^\textrm{exp}$ and "corrected" kurtosis $\gamma_{2,\textrm{corr}}^\textrm{exp}$ as functions of centrality in PbPb collisions at $\sqrt{s_{NN}}=5.02~\textrm{Te}\textrm{V}$. The measurement is performed with charged particles within the acceptance region.

The magnitudes of the measured superskewness $\gamma_3^\textrm{exp}$ and "corrected" superskewness $\gamma_{3,\textrm{corr}}^\textrm{exp}$ as functions of centrality in PbPb collisions at $\sqrt{s_{NN}}=5.02~\textrm{Te}\textrm{V}$. The measurement is performed with charged particles within the acceptance region.

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.

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>

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.

Transverse dilepton-$p_{T}^{miss}$ mass pre-fit distribution for selected events in SR2 of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Dilepton mass pre-fit distribution for selected events in SR3 of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Dilepton $p_T$ pre-fit distribution for selected events in SR3 of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Dijet $p_T$ pre-fit distribution for selected events in the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Lepton-dijet $p_T$ pre-fit distribution for selected events in the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Maximum $p_T$ of the visible particles divided by the missing $p_T$, for selected events in the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between jets, for selected events in the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between the lepton and missing $p_T$, for selected events in the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between the lepton and the dijet system, for selected events in the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta R$ distance between the leptons, for selected events in the non-prompt validation region of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between the di-lepton system and the missing $p_T$, for selected events in the non-prompt validation region of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between the lepton and the di-jet system, for selected events in the non-prompt validation region of the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between the lepton-dijet system and the missing $p_T$, for selected events in the non-prompt validation region of the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Leading lepton $p_T$ pre-fit distribution for selected events in the top control region of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Trailing lepton $p_T$ pre-fit distribution for selected events in the top control region of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Dilepton mass pre-fit distribution for selected events in the Drell-Yan control region of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Missing $p_T$ pre-fit distribution for selected events in the Drell-Yan control region of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Missing $p_T$ pre-fit distribution for selected events in the $WW$ control region of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Dilepton $p_T$ pre-fit distribution for selected events in the $WW$ control region of the di-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Minimum $p_T$ of the visible particles divided by the missing $p_T$, for selected events in the top control region of the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between jets, for selected events in the top control region of the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between the lepton and missing $p_T$, for selected events in the $W$+jets control region of the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Pre-fit distribution of the $\Delta \phi$ angle between the lepton and the dijet system, for selected events in the $W$+jets control region of the semi-leptonic channel for the full dataset. The error bars on the data points represent the statistical uncertainty of the data, and the error bars on the predicted yields represent the combined systematic and statistical uncertainty in each bin. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$. The last bin includes the overflow.

Unrolled $m_{ll}-m_{T}^{l min, p_{T}^{miss}}$ post-fit distributions in the di-leptonic channel in SR1 for the full data set. The histogram bins are spaced uniformly. Each group of five bins (from left to right) corresponds to $m_{T}^{l min, p_{T}^{miss}}$ distribution in a $m_{ll}$ region, placed in ascending order. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$.

Unrolled $m_{ll}-m_{T}^{l min, p_{T}^{miss}}$ post-fit distributions in the di-leptonic channel in SR2 for the full data set. The histogram bins are spaced uniformly. Each group of five bins (from left to right) corresponds to $m_{T}^{l min, p_{T}^{miss}}$ distribution in a $m_{ll}$ region, placed in ascending order. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$.

Unrolled $m_{ll}-m_{T}^{l min, p_{T}^{miss}}$ post-fit distributions in the di-leptonic channel in SR3 for the full data set. The histogram bins are spaced uniformly. Each group of five bins (from left to right) corresponds to $m_{T}^{l min, p_{T}^{miss}}$ distribution in a $m_{ll}$ region, placed in ascending order. The signal prediction represents the mass point: $m_s = 160 GeV, m_{\chi} = 100 GeV, m_{Z'} = 500 GeV$.

Post-fit BDT distributions in the semi-leptonic channel for the full data set in the top control region. The signal prediction represents the mass point: $m_{s} = 160$ GeV, $m_{\chi} = 100$ GeV, $m_{Z'} = 500$ GeV

Post-fit BDT distributions in the semi-leptonic channel for the full data set in the $W$+jets control region. The signal prediction represents the mass point: $m_{s} = 160$ GeV, $m_{\chi} = 100$ GeV, $m_{Z'} = 500$ GeV

Post-fit BDT distributions in the semi-leptonic channel in the 2016 signal region. The signal prediction represents the mass point: $m_{s} = 160$ GeV, $m_{\chi} = 100$ GeV, $m_{Z'} = 500$ GeV

Post-fit BDT distributions in the semi-leptonic channel in the 2017+2018 signal region. The signal prediction represents the mass point: $m_{s} = 160$ GeV, $m_{\chi} = 100$ GeV, $m_{Z'} = 500$ GeV

Contour of the 95% CL observed upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 100 GeV$.

Contour of the 95% CL expected + 2 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 100 GeV$.

Contour of the 95% CL expected + 1 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 100 GeV$.

Contour of the 95% CL expected upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 100 GeV$.

Contour of the 95% CL expected - 1 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 100 GeV$.

Contour of the 95% CL expected - 2 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 100 GeV$.

Contour of the 95% CL observed upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 150 GeV$.

Contour of the 95% CL expected + 2 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 150 GeV$.

Contour of the 95% CL expected + 1 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 150 GeV$.

Contour of the 95% CL expected upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 150 GeV$.

Contour of the 95% CL expected - 1 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 150 GeV$.

Contour of the 95% CL expected - 2 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 150 GeV$.

Contour of the 95% CL observed upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 200 GeV$.

Contour of the 95% CL expected + 2 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 200 GeV$.

Contour of the 95% CL expected + 1 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 200 GeV$.

Contour of the 95% CL expected upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 200 GeV$.

Contour of the 95% CL expected - 1 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 200 GeV$.

Contour of the 95% CL expected - 2 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 200 GeV$.

Contour of the 95% CL observed upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 300 GeV$.

Contour of the 95% CL expected + 2 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 300 GeV$.

Contour of the 95% CL expected + 1 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 300 GeV$.

Contour of the 95% CL expected upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 300 GeV$.

Contour of the 95% CL expected - 1 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 300 GeV$.

Contour of the 95% CL expected - 2 std. upper limits representing $\sigma / \sigma_{theory} = 1$ for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 300 GeV$.

The $\sigma / \sigma_{theory}$ 95% CL upper limit values for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 100 GeV$.

The $\sigma / \sigma_{theory}$ 95% CL upper limit values for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 150 GeV$.

The $\sigma / \sigma_{theory}$ 95% CL upper limit values for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 200 GeV$.

The $\sigma / \sigma_{theory}$ 95% CL upper limit values for the dark Higgs model in the ($m_{s}$, $m_{Z'}$) plane for $m_{\chi} = 300 GeV$.

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.

CF vs qinv

CF vs qinv

CF vs qinv

CF vs qinv

CF vs qinv

CF vs qinv

CF vs qinv

CF vs qinv

CF vs qinv

CF vs qinv

CF vs qinv

CF vs qinv

CF vs qinv

CF vs qinv

CF vs qinv

CF vs qinv

CF vs qinv

CF vs qinv

CF vs qinv

CF vs qinv

CF vs qinv

CF vs qinv

CF vs qinv

CF vs qinv

CF vs qinv

CF vs qinv

CF vs qinv

CF vs qinv

CF vs qinv

CF vs qinv

CF vs qinv

CF vs qinv

CF vs qinv

CF vs qinv

CF vs qinv

CF vs qinv

CF vs qinv

CF vs qinv

CF vs qinv

Rinv vs mT

Rinv vs mT

Rinv vs mT

Rinv vs mT

Rinv vs mT

Rinv vs mT

Rinv vs mT

Rinv vs mT

Rinv vs mT

Rinv vs mT

Rinv vs mT

Rinv vs mT

Rinv vs mT

Rinv vs mT

Rinv vs mT

Rinv vs mT

Rinv vs mT

Rinv vs mT

Rinv vs mT

Rinv vs mT

lambda vs mT

lambda vs mT

lambda vs mT

lambda vs mT

lambda vs mT

lambda vs mT

lambda vs mT

lambda vs mT

lambda vs mT

lambda vs mT

lambda vs mT

lambda vs mT

$R_\mathrm{T}$ distribution using events with trigger particles $5<p_\mathrm{T}^\mathrm{trig}<40~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ in pp collisions at $\sqrt{s}=13~\mathrm{TeV}$

$p_\mathrm{T}$ as a function of $R_\mathrm{T}$ in the toward region for different $R_\mathrm{T}$ intervals using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for pp collisions at $\sqrt{s}=5.02~\mathrm{TeV}$

$p_\mathrm{T}$ as a function of $R_\mathrm{T}$ in the away region for different $R_\mathrm{T}$ intervals using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for pp collisions at $\sqrt{s}=5.02~\mathrm{TeV}$

$p_\mathrm{T}$ as a function of $R_\mathrm{T}$ in the transverse region for different $R_\mathrm{T}$ intervals using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for pp collisions at $\sqrt{s}=5.02~\mathrm{TeV}$

$p_\mathrm{T}$ as a function of $R_\mathrm{T}$ in the toward region for different $R_\mathrm{T}$ intervals using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for p-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$

$p_\mathrm{T}$ as a function of $R_\mathrm{T}$ in the away region for different $R_\mathrm{T}$ intervals using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for p-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$

$p_\mathrm{T}$ as a function of $R_\mathrm{T}$ in the transverse region for different $R_\mathrm{T}$ intervals using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for p-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$

$p_\mathrm{T}$ as a function of $R_\mathrm{T}$ in the toward region for different $R_\mathrm{T}$ intervals using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$

$p_\mathrm{T}$ as a function of $R_\mathrm{T}$ in the away region for different $R_\mathrm{T}$ intervals using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$

$p_\mathrm{T}$ as a function of $R_\mathrm{T}$ in the transverse region for different $R_\mathrm{T}$ intervals using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$

Average $p_\mathrm{T}$ as a function of $R_\mathrm{T}$ in the toward region using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for pp collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$

Average $p_\mathrm{T}$ as a function of $R_\mathrm{T}$ in the toward region using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for p-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$

Average $p_\mathrm{T}$ as a function of $R_\mathrm{T}$ in the toward region using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$

Average $p_\mathrm{T}$ as a function of $R_\mathrm{T}$ in the away region using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for pp collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$

Average $p_\mathrm{T}$ as a function of $R_\mathrm{T}$ in the away region using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for p-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$

Average $p_\mathrm{T}$ as a function of $R_\mathrm{T}$ in the away region using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$

Average $p_\mathrm{T}$ as a function of $R_\mathrm{T}$ in the transverse region using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for pp collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$

Average $p_\mathrm{T}$ as a function of $R_\mathrm{T}$ in the transverse region using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for p-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$

Average $p_\mathrm{T}$ as a function of $R_\mathrm{T}$ in the transverse region using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$

$p_\mathrm{T}$-integrated yield normalised to $<N_\mathrm{ch}^\mathrm{T}>$ as a function of $R_\mathrm{T}$ in the toward region using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for pp collisions at $\sqrt{s}=5.02~\mathrm{TeV}$

$p_\mathrm{T}$-integrated yield normalised to $<N_\mathrm{ch}^\mathrm{T}>$ as a function of $R_\mathrm{T}$ in the away region using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for pp collisions at $\sqrt{s}=5.02~\mathrm{TeV}$

$p_\mathrm{T}$-integrated yield normalised to $<N_\mathrm{ch}^\mathrm{T}>$ as a function of $R_\mathrm{T}$ in the toward region using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for p-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$

$p_\mathrm{T}$-integrated yield normalised to $<N_\mathrm{ch}^\mathrm{T}>$ as a function of $R_\mathrm{T}$ in the away region using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for p-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$

$p_\mathrm{T}$-integrated yield normalised to $<N_\mathrm{ch}^\mathrm{T}>$ as a function of $R_\mathrm{T}$ in the toward region using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$

$p_\mathrm{T}$-integrated yield normalised to $<N_\mathrm{ch}^\mathrm{T}>$ as a function of $R_\mathrm{T}$ in the away region using events with trigger particles $8<p_\mathrm{T}^\mathrm{trig}<15~\mathrm{GeV}/c$ in the pseudorapidity range of $|\eta|<0.8$ and with $p_\mathrm{T}>0.5~\mathrm{GeV}/c$ for Pb-Pb collisions at $\sqrt{s_\mathrm{NN}}=5.02~\mathrm{TeV}$

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

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

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

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

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

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

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

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

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

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

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

Observed upper limits at 90% CL on the square of the kinetic mixing coefficient $\epsilon$ in the minimal model of a dark photon

Observed upper limits at 90% CL on the mixing angle $\theta_H$ for the type IV 2HDM+S scenario with fixed $tan\beta=0.5$

The combined efficiency of the dimuon scouting trigger and the MVA muon selection, averaged between 2017 and 2018, weighted by the integrated luminosity of each year. The trigger and selection efficiencies are also included in the table.

Theory cross-section times branching fraction to muons times acceptance for the dark photon and 2HDM+S models.

'$\overline{\Lambda}$ $D_{LL}$ as a function of $p_{T}$ at $0 < \eta_{\Lambda(\overline{\Lambda})} < 1.2$'

'$\Lambda$ $D_{LL}$ as a function of $p_{T}$ at $-1.2 < \eta_{\Lambda(\overline{\Lambda})} < 0$'

'$\overline{\Lambda}$ $D_{LL}$ as a function of $p_{T}$ at $-1.2 < \eta_{\Lambda(\overline{\Lambda})} < 0$'

'$\Lambda$ $D_{LL}$ as a function of $p_{T}$ at $0 < \eta_{\Lambda(\overline{\Lambda})} < 1.2$'

'$\overline{\Lambda}$ $D_{LL}$ as a function of $p_{T}$ at $0 < \eta_{\Lambda(\overline{\Lambda})} < 1.2$'

'Two-year-combined $\Lambda$ $D_{LL}$ as a function of $p_{T}$ at $0 < \eta_{\Lambda(\overline{\Lambda})} < 1.2$'

'Two-year-combined $\overline{\Lambda}$ $D_{LL}$ as a function of $p_{T}$ at $0 < \eta_{\Lambda(\overline{\Lambda})} < 1.2$'

'Two-year-combined $\Lambda+\overline{\Lambda}$ $D_{LL}$ as a function of $p_{T}$ at $0 < \eta_{\Lambda(\overline{\Lambda})} < 1.2$'

'Two-year-combined $\Lambda$ $D_{LL}$ as a function of $p_{T}$ at $-1.2 < \eta_{\Lambda(\overline{\Lambda})} < 0$'

'Two-year-combined $\overline{\Lambda}$ $D_{LL}$ as a function of $p_{T}$ at $-1.2 < \eta_{\Lambda(\overline{\Lambda})} < 0$'

'$\Lambda$ $D_{LL}$ as a function of z at $0 < \eta_{jet} < 1.0$'

'$\overline{\Lambda}$ $D_{LL}$ as a function of z at $0 < \eta_{jet} < 1.0$'

'$\Lambda$ $D_{LL}$ as a function of z at $-1.0 < \eta_{jet} < 0$'

'$\overline{\Lambda}$ $D_{LL}$ as a function of z at $-1.0 < \eta_{jet} < 0$'

'particle jet pt as a function of $\Lambda$ z'

'particle jet pt as a function of $\overline{\Lambda}$ z'

'$\Lambda$ $D_{TT}$ as a function of $p_{T}$ at $0 < \eta_{\Lambda(\overline{\Lambda})} < 1.2$'

'$\overline{\Lambda}$ $D_{TT}$ as a function of $p_{T}$ at $0 < \eta_{\Lambda(\overline{\Lambda})} < 1.2$'

'$\Lambda$ $D_{TT}$ as a function of $p_{T}$ at $-1.2 < \eta_{\Lambda(\overline{\Lambda})} < 0$'

'$\overline{\Lambda}$ $D_{TT}$ as a function of $p_{T}$ at $-1.2 < \eta_{\Lambda(\overline{\Lambda})} < 0$'

'$\Lambda$ $D_{TT}$ as a function of $p_{T}$ at $0 < \eta_{\Lambda(\overline{\Lambda})} < 1.2$'

'$\overline{\Lambda}$ $D_{TT}$ as a function of $p_{T}$ at $0 < \eta_{\Lambda(\overline{\Lambda})} < 1.2$'

'Two-year-combined $\Lambda$ $D_{TT}$ as a function of $p_{T}$ at $0 < \eta_{\Lambda(\overline{\Lambda})} < 1.2$'

'Two-year-combined $\overline{\Lambda}$ $D_{TT}$ as a function of $p_{T}$ at $0 < \eta_{\Lambda(\overline{\Lambda})} < 1.2$'

'Two-year-combined $\Lambda$ $D_{TT}$ as a function of $p_{T}$ at $-1.2 < \eta_{\Lambda(\overline{\Lambda})} < 0$'

'Two-year-combined $\overline{\Lambda}$ $D_{TT}$ as a function of $p_{T}$ at $-1.2 < \eta_{\Lambda(\overline{\Lambda})} < 0$'

'$\Lambda$ $D_{TT}$ as a function of z at $0 < \eta_{jet} < 1.0$'

'$\overline{\Lambda}$ $D_{TT}$ as a function of z at $0 < \eta_{jet} < 1.0$'

'$\Lambda$ $D_{TT}$ as a function of z at $-1.0 < \eta_{jet} < 0$'

'$\overline{\Lambda}$ $D_{TT}$ as a function of z at $-1.0 < \eta_{jet} < 0$'

'particle jet pt as a function of $\Lambda$ z'

'particle jet pt as a function of $\overline{\Lambda}$ z'

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.

Correlation matrix of the unfolding uncertainty between the inclusive-jet measurement and the previous dijet measurement. Correlations are given in percent.

Summary of different determinations of $\alpha_s(M_Z^2)$ at NNLO or higher order, adapted from PDG, see references therein. The red points are included in the PDG world average. The averages from each sub-field are shown as yellow bands and the world average as a blue band. A recent measurement from CMS using jet cross sections and the latest determination from HERAPDF, which are not yet included in the world average, are shown in green. The current determination, assuming half-correlated and half-uncorrelated scale uncertainties, is shown in black.

Value of the strong coupling $\alpha_s(\mu^2)$ as a function of the scale $\mu$. The data points indicate determinations from measurements that were performed close to the indicated scale. The uncertainties represent the full uncertainty of each determination. All depicted results were obtained at least at NNLO. They are based on data from $e^+e^-$, $ep$ and $pp$ collisions, as well as from $\tau$ lepton decays and quarkonium states. The solid blue line shows the PDG world average. Also shown are the $\alpha_s(M_Z^2)$ values corresponding to each data point.

Normalized $\Delta\it{v}_{2}$ slope of protons as a function of centrality in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

Normalized $\Delta\it{v}_{3}$ slope of charged hadrons as a function of centrality in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

Normalized $\Delta\it{v}_{3}$ slope of kaons as a function of centrality in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

Normalized $\Delta\it{v}_{3}$ slope of pions as a function of centrality in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

Normalized $\Delta\it{v}_{3}$ slope of protons as a function of centrality in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

Delta Integral Covariance as function of $v_2$ in 0-10% Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

Delta Integral Covariance as function of $v_2$ in 10-20% Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

Delta Integral Covariance as function of $v_2$ in 20-30% Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

Delta Integral Covariance as function of $v_2$ in 30-40% Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

Delta Integral Covariance as function of $v_2$ in 40-50% Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

Delta Integral Covariance as function of $v_2$ in 50-60% Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

$f_{CMW}$ as function of CENTRALITY in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

Non-prompt J/$\psi$ fraction as a function of average number of participants in Pb-Pb at 5.02 TeV

Non-prompt J/$\psi$ yields as a function of transverse momentum in Pb-Pb at 5.02 TeV, centrality 0-10%

Prompt J/$\psi$ yields as a function of transverse momentum in Pb-Pb at 5.02 TeV, centrality 0-10%

Prompt J/$\psi$ nuclear modification factor as a function of transverse momentum in Pb-Pb at 5.02 TeV, centrality 0-10%. Global uncertainty 2.3%

Prompt J/$\psi$ nuclear modification factor as a function of transverse momentum in Pb-Pb at 5.02 TeV, centrality 10-30%. Global uncertainty 2.3%

Prompt J/$\psi$ nuclear modification factor as a function of transverse momentum in Pb-Pb at 5.02 TeV, centrality 30-50%. Global uncertainty 2.7%

Prompt J/$\psi$ nuclear modification factor as a function of average number of participants in Pb-Pb at 5.02 TeV. Global uncertainty 2.2%

Non-prompt J/$\psi$ nuclear modification factor as a function of transverse momentum in Pb-Pb at 5.02 TeV, centrality 0-10%. Global uncertainty 2.3%

Non-prompt J/$\psi$ nuclear modification factor as a function of transverse momentum in Pb-Pb at 5.02 TeV,, centrality 10-30%. Global uncertainty 2.3%

Non-prompt J/$\psi$ nuclear modification factor as a function of transverse momentum in Pb-Pb at 5.02 TeV, centrality 30-50%. Global uncertainty 2.7%

Non-prompt J/$\psi$ nuclear modification factor as a function of average number of participants in Pb-Pb at 5.02 TeV. Global uncertainty 2.7%

Prompt J/$\psi$ yields as a function of transverse momentum in Pb-Pb at 5.02 TeV, centrality 30-50%

Non-prompt J/$\psi$ yields as a function of transverse momentum in Pb-Pb at 5.02 TeV, centrality 30-50%

$p_{\rm T}$-distributions of $\rm{K}^{*}$ (average of particle and anti-particle) meson measured in Pb-Pb collisions at \sqrt{s_{NN}}$ = 5.02 TeV for 40-60\% centrality.

$p_{\rm T}$-distributions of $\rm{K}^{*}$ (average of particle and anti-particle) meson measured in Pb-Pb collisions at \sqrt{s_{NN}}$ = 5.02 TeV for 60-80\% centrality.

$p_{\rm T}$-integrated yield d$N$/d$y$ of K$^{*}$ (average of particle and anti-particle) meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

Mean $p_{\rm T}$ of K$^{*}$ (average of particle and anti-particle) meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

$p_{\rm T}$-integrated yield ratio of $\rm{K}^{*}$ (average of particleand anti-particle) to K (average of particle and anti-particle) measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV.

$p_{\rm T}$ differential ratio $(\rm{K}^{*\pm} + \rm{K}^{*\mp}) / (\rm{K}^{+}+\rm{K}^{-})$ measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV for 0-10\% centrality interval.

$p_{\rm T}$ differential ratio $(\rm{K}^{*\pm} + \rm{K}^{*\mp}) / (\rm{K}^{+}+\rm{K}^{-})$ measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV for 60-80\% centrality interval.

$p_{\rm T}$ differential ratio $(\rm{K}^{*\pm} + \rm{K}^{*\mp}) / (\rm{pi}^{+}+\rm{pi}^{-})$ measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV for 0-10\% centrality interval.

$p_{\rm T}$ differential ratio $(\rm{K}^{*\pm} + \rm{K}^{*\mp}) / (\rm{pi}^{+}+\rm{pi}^{-})$ measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV for 60-80\% centrality interval.

$p_{T}$-dependent nuclear modification factor of $\rm{K}^{*}$ (average of particle and anti-particle) meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV for 0-10\% centrality.

$p_{T}$-dependent nuclear modification factor of $\rm{K}^{*}$ (average of particle and anti-particle) meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV for 10-20\% centrality.

$p_{T}$-dependent nuclear modification factor of $\rm{K}^{*}$ (average of particle and anti-particle) meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV for 20-40\% centrality.

$p_{T}$-dependent nuclear modification factor of $\rm{K}^{*}$ (average of particle and anti-particle) meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV for 40-60\% centrality.

$p_{T}$-dependent nuclear modification factor of $\rm{K}^{*}$ (average of particle and anti-particle) meson measured in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV for 60-80\% centrality.

Statistical correlation matrix for the one-dimensional $y^{Z}$ measurement.

Statistical correlation matrix for the one-dimensional $p_{T}^{Z}$ measurement.

Statistical correlation matrix for the one-dimensional $\phi_{\eta}^{*}$ measurement.

Correlation matrix for the efficiency uncertainty of the one-dimensional $y^{Z}$ measurement.

Correlation matrix for the efficiency uncertainty of the one-dimensional $p_{T}^{Z}$ measurement.

Correlation matrix for the efficiency uncertainty of the one-dimensional $\phi_{\eta}^{*}$ measurement.

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