$R_{cp}$ as a function of $\eta$ for 1.5 < $p_T$ < 4.0 GeV/$c$ for different centrality classes. Systematic uncertainties which are point-to-point uncorrelated (sys-uncorr) and correlated (sys-corr) are shown.

Observed $m_{t\bar{t}}^{reco}$ distributions in data for SR$2b$, shown together with the result of the fit with the three-shape-parameter function. The error bars indicate the effect of the fit parameter uncertainty on the background prediction. The bin width of the distributions is chosen to be the same as that used in the background parameterization.

Observed and expected upper limits on the cross-section times branching fraction of the $Z'_{TC2} \rightarrow t\bar{t}$ as a function of the $Z'_{TC2}$ mass ($\Gamma/m = 1.2$%).

Fit parameter values obtained from fits to data in the two signal regions. The fit function is $f(x) = p_{0} (1-x)^{p_{1}} x^{p_{2} + p_{3} \log(x)}$, where $x = m^{reco}_{t\bar{t}}/\sqrt{s}$.

Distributions of $s_{\text{T}}$ in the single-tau one-bin SR. The stacked histograms show the various SM background contributions. The hatched band indicates the total statistical and systematic uncertainty of the SM background. The $t\bar{t}$ (2 real $\tau$) and $t\bar{t}$ (1 real $\tau$) as well as the single-top background contributions are scaled with the normalization factors obtained from the background-only fit. Minor backgrounds are grouped together and denoted as 'Other'. This includes $t\bar{t}$-fake, single top, and other top (di-tau channel) or $t\bar{t}$-fake, $t\bar{t}+H$, multiboson, and other top (single-tau channel). The overlaid dotted lines show the additional contributions for signal scenarios close to the expected exclusion contour with the particle type and the mass and $\beta$ parameters for the simplified models indicated in the legend. For the leptoquark signal model the shapes of the distributions for $\text{LQ}_{3}^{\text{d}}$ and $\text{LQ}_{3}^{\text{v}}$ (not shown) are similar to that of $\text{LQ}_{3}^{\text{u}}$. The rightmost bin includes the overflow.

Distributions of $m_{\text{T}}(\tau)$ in the single-tau one-bin SR. The stacked histograms show the various SM background contributions. The hatched band indicates the total statistical and systematic uncertainty of the SM background. The $t\bar{t}$ (2 real $\tau$) and $t\bar{t}$ (1 real $\tau$) as well as the single-top background contributions are scaled with the normalization factors obtained from the background-only fit. Minor backgrounds are grouped together and denoted as 'Other'. This includes $t\bar{t}$-fake, single top, and other top (di-tau channel) or $t\bar{t}$-fake, $t\bar{t}+H$, multiboson, and other top (single-tau channel). The overlaid dotted lines show the additional contributions for signal scenarios close to the expected exclusion contour with the particle type and the mass and $\beta$ parameters for the simplified models indicated in the legend. For the leptoquark signal model the shapes of the distributions for $\text{LQ}_{3}^{\text{d}}$ and $\text{LQ}_{3}^{\text{v}}$ (not shown) are similar to that of $\text{LQ}_{3}^{\text{u}}$. The rightmost bin includes the overflow.

Distributions of $\Sigma m_{\text{T}}(b_{1,2})$ in the single-tau $p_{\text{T}}(\tau)$-binned SR. The stacked histograms show the various SM background contributions. The hatched band indicates the total statistical and systematic uncertainty of the SM background. The $t\bar{t}$ (2 real $\tau$) and $t\bar{t}$ (1 real $\tau$) as well as the single-top background contributions are scaled with the normalization factors obtained from the background-only fit. Minor backgrounds are grouped together and denoted as 'Other'. This includes $t\bar{t}$-fake, single top, and other top (di-tau channel) or $t\bar{t}$-fake, $t\bar{t}+H$, multiboson, and other top (single-tau channel). The overlaid dotted lines show the additional contributions for signal scenarios close to the expected exclusion contour with the particle type and the mass and $\beta$ parameters for the simplified models indicated in the legend. For the leptoquark signal model the shapes of the distributions for $\text{LQ}_{3}^{\text{d}}$ and $\text{LQ}_{3}^{\text{v}}$ (not shown) are similar to that of $\text{LQ}_{3}^{\text{u}}$. The rightmost bin includes the overflow.

Distributions of $p_{\text{T}}(\tau)$ in the single-tau $p_{\text{T}}(\tau)$-binned SR. The stacked histograms show the various SM background contributions. The hatched band indicates the total statistical and systematic uncertainty of the SM background. The $t\bar{t}$ (2 real $\tau$) and $t\bar{t}$ (1 real $\tau$) as well as the single-top background contributions are scaled with the normalization factors obtained from the background-only fit. Minor backgrounds are grouped together and denoted as 'Other'. This includes $t\bar{t}$-fake, single top, and other top (di-tau channel) or $t\bar{t}$-fake, $t\bar{t}+H$, multiboson, and other top (single-tau channel). The overlaid dotted lines show the additional contributions for signal scenarios close to the expected exclusion contour with the particle type and the mass and $\beta$ parameters for the simplified models indicated in the legend. For the leptoquark signal model the shapes of the distributions for $\text{LQ}_{3}^{\text{d}}$ and $\text{LQ}_{3}^{\text{v}}$ (not shown) are similar to that of $\text{LQ}_{3}^{\text{u}}$. The rightmost bin includes the overflow.

Observed event yields in data ('Observed') and expected event yields for SM background processes obtained from the background-only fit ('Total bkg.' and rows below) in the signal regions of the di-tau and single-tau channels. The quoted uncertainties include both the statistical and systematic uncertainties and are truncated at zero yield. By construction, no $t\bar{t}$ (2 real $\tau$) events can pass the selections in the single-tau channel. As the individual uncertainties are correlated, they do not add in quadrature to equal the total background uncertainty.

From left to right: upper limits at the 95% confidence level (CL) on the visible cross section ($\sigma_\text{vis}$) and on the number of signal events ($S_{\text{obs}}^{95}$). The third column ($S_{\text{exp}}^{95}$) shows the upper limit at the 95% CL on the number of signal events, given the expected number (and $\pm 1\,\sigma$ excursions on the expectation) of background events. The last two columns indicate the confidence level observed for the background-only hypothesis ($\text{CL}_{b}$), the discovery $p$-value ($p(s=0)$) and the significance ($Z$). In the di-tau SR, where fewer events are observed than predicted by the fitted background estimate, the $p$-value is capped at 0.5.

Expected and observed exclusion contours at the 95% confidence level for the vector third-generation leptoquark signal model, as a function of the mass $m(\text{LQ}_{3}^{\text{v}})$ and the branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ into a quark and a charged lepton. The plot shows the exclusion contour for the minimal-coupling scenario. The limits are derived from the binned single-tau signal region.

Expected and observed exclusion contours at the 95% confidence level for the vector third-generation leptoquark signal model, as a function of the mass $m(\text{LQ}_{3}^{\text{v}})$ and the branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ into a quark and a charged lepton. The plot shows the exclusion contour for the minimal-coupling scenario. The limits are derived from the binned single-tau signal region.

Expected and observed exclusion contours at the 95% confidence level for the vector third-generation leptoquark signal model, as a function of the mass $m(\text{LQ}_{3}^{\text{v}})$ and the branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ into a quark and a charged lepton. The plot shows the exclusion contour for vector leptoquarks with additional gauge couplings. The limits are derived from the binned single-tau signal region.

Expected and observed exclusion contours at the 95% confidence level for the vector third-generation leptoquark signal model, as a function of the mass $m(\text{LQ}_{3}^{\text{v}})$ and the branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ into a quark and a charged lepton. The plot shows the exclusion contour for vector leptoquarks with additional gauge couplings. The limits are derived from the binned single-tau signal region.

Exclusion contours at the 95% confidence level for the stop-stau signal model as a function of the masses of the top squark $m(\tilde{t}_{1})$ and of the tau slepton $m(\tilde{\tau}_{1})$. Expected and observed limits are shown for the present search in comparison to observed limits from previous ATLAS analyses based on data from Run-1 of the LHC at $\sqrt{s} = 8$ TeV [Eur. Phys. J. C 76 (2016)] and on a partial dataset from Run 2 at $\sqrt{s} = 13$ TeV [Phys. Rev. D 98 (2018) 032008]. The green band indicates the limit on the mass of the tau slepton (for a massless LSP) from the LEP experiments.

Exclusion contours at the 95% confidence level for the stop-stau signal model as a function of the masses of the top squark $m(\tilde{t}_{1})$ and of the tau slepton $m(\tilde{\tau}_{1})$. Expected and observed limits are shown for the present search in comparison to observed limits from previous ATLAS analyses based on data from Run-1 of the LHC at $\sqrt{s} = 8$ TeV [Eur. Phys. J. C 76 (2016)] and on a partial dataset from Run 2 at $\sqrt{s} = 13$ TeV [Phys. Rev. D 98 (2018) 032008]. The green band indicates the limit on the mass of the tau slepton (for a massless LSP) from the LEP experiments.

Expected and observed exclusion contours at the 95% confidence level for the scalar third-generation leptoquark signal model, as a function of the mass $m(\text{LQ}_{3}^{\text{u}})$ and the branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow q\ell)$ into a quark and a charged lepton. The plot shows the exclusion contour for up-type leptoquarks $\text{LQ}_{3}^{\text{u}})$ with charge $+2/3e$. The limits are derived from the binned single-tau signal region. Shown in gray for comparison are the observed exclusion-limit contours from the previous ATLAS publication that targets the same leptoquark models but is based on a subset of the Run-2 data [JHEP 06 (2019) 144]. In this previous publication five different analyses are considered that target not only the final state studied here but also the final states that correspond to a branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow q\ell)$ of 0 or 1, leading to the concave shapes of the gray exclusion contours.

Expected and observed exclusion contours at the 95% confidence level for the scalar third-generation leptoquark signal model, as a function of the mass $m(\text{LQ}_{3}^{\text{u}})$ and the branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow q\ell)$ into a quark and a charged lepton. The plot shows the exclusion contour for up-type leptoquarks $\text{LQ}_{3}^{\text{u}})$ with charge $+2/3e$. The limits are derived from the binned single-tau signal region. Shown in gray for comparison are the observed exclusion-limit contours from the previous ATLAS publication that targets the same leptoquark models but is based on a subset of the Run-2 data [JHEP 06 (2019) 144]. In this previous publication five different analyses are considered that target not only the final state studied here but also the final states that correspond to a branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow q\ell)$ of 0 or 1, leading to the concave shapes of the gray exclusion contours.

Expected and observed exclusion contours at the 95% confidence level for the scalar third-generation leptoquark signal model, as a function of the mass $m(\text{LQ}_{3}^{\text{d}})$ and the branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow q\ell)$ into a quark and a charged lepton. The plot shows the exclusion contour for down-type leptoquarks $\text{LQ}_{3}^{\text{d}})$ with charge $-1/3e$. The limits are derived from the binned single-tau signal region. Shown in gray for comparison are the observed exclusion-limit contours from the previous ATLAS publication that targets the same leptoquark models but is based on a subset of the Run-2 data [JHEP 06 (2019) 144]. In this previous publication five different analyses are considered that target not only the final state studied here but also the final states that correspond to a branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow q\ell)$ of 0 or 1, leading to the concave shapes of the gray exclusion contours.

Expected and observed exclusion contours at the 95% confidence level for the scalar third-generation leptoquark signal model, as a function of the mass $m(\text{LQ}_{3}^{\text{d}})$ and the branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow q\ell)$ into a quark and a charged lepton. The plot shows the exclusion contour for down-type leptoquarks $\text{LQ}_{3}^{\text{d}})$ with charge $-1/3e$. The limits are derived from the binned single-tau signal region. Shown in gray for comparison are the observed exclusion-limit contours from the previous ATLAS publication that targets the same leptoquark models but is based on a subset of the Run-2 data [JHEP 06 (2019) 144]. In this previous publication five different analyses are considered that target not only the final state studied here but also the final states that correspond to a branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow q\ell)$ of 0 or 1, leading to the concave shapes of the gray exclusion contours.

Upper limits on the signal cross section at the 95 % confidence level for the stop-stau signal model.

Upper limits on the signal cross section at the 95 % confidence level for the scalar third-generation leptoquark signal model with up-type leptoquarks.

Upper limits on the signal cross section at the 95 % confidence level for the scalar third-generation leptoquark signal model with down-type leptoquarks.

Upper limits on the signal cross section at the 95 % confidence level for the vector third-generation leptoquark signal model with minimal coupling (MC).

Upper limits on the signal cross section at the 95 % confidence level for the vector third-generation leptoquark signal model with additional gauge couplings (YM).

Acceptance of the one-bin signal region of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$.

Efficiency of the one-bin signal region of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$.

Efficiency of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$.

Efficiency of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$.

Efficiency of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the signal region of the di-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$.

Efficiency of the signal region of the di-tau channel for pair production of up-type leptoquarks $\text{LQ}_{3}^{\text{u}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{u}} \rightarrow b\tau)$ of 0 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the one-bin signal region of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$.

Efficiency of the one-bin signal region of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow t\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$.

Efficiency of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow t\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$.

Efficiency of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow t\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$.

Efficiency of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow t\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the signal region of the di-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$.

Efficiency of the signal region of the di-tau channel for pair production of down-type leptoquarks $\text{LQ}_{3}^{\text{d}}$. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{d}} \rightarrow t\tau)$ of 0 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the one-bin signal region of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario.

Efficiency of the one-bin signal region of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario.

Efficiency of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario.

Efficiency of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario.

Efficiency of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the signal region of the di-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario.

Efficiency of the signal region of the di-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ in the minimal-coupling scenario. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the one-bin signal region of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings.

Efficiency of the one-bin signal region of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings.

Efficiency of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings.

Efficiency of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings.

Efficiency of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 or 1 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the signal region of the di-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings.

Efficiency of the signal region of the di-tau channel for pair production of vector leptoquarks $\text{LQ}_{3}^{\text{v}}$ with additional gauge couplings. The plot does not show efficiencies for a branching fraction $B(\text{LQ}_{3}^{\text{v}} \rightarrow b\tau)$ of 0 because here the acceptance at generator level becomes zero and the efficiency is thus undefined.

Acceptance of the one-bin signal region of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Efficiency of the one-bin signal region of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Acceptance of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Efficiency of the first bin of the multi-bin signal region (50 GeV $< p_{\text{T}}(\tau) <$ 100 GeV) of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Acceptance of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Efficiency of the middle bin of the multi-bin signal region (100 GeV $< p_{\text{T}}(\tau) <$ 200 GeV) of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Acceptance of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Efficiency of the last bin of the multi-bin signal region (200 GeV $< p_{\text{T}}(\tau)$) of the single-tau channel for pair production of top squarks with decays via tau sleptons.

Acceptance of the signal region of the di-tau channel for pair production of top squarks with decays via tau sleptons.

Efficiency of the signal region of the di-tau channel for pair production of top squarks with decays via tau sleptons.

Cutflow for the benchmark signal model $m(\tilde{t}_{1}) = 1350$ GeV, $m(\tilde{\tau}_{1}) = 1090$ GeV for the di-tau SR. The simulated sample contains 30,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the di-tau channel.

Cutflow for the benchmark signal model $m(\tilde{t}_{1}) = 1350$ GeV, $m(\tilde{\tau}_{1}) = 1090$ GeV for the single-tau one-bin SR. The simulated sample contains 30,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\tilde{t}_{1}) = 1350$ GeV, $m(\tilde{\tau}_{1}) = 1090$ GeV for the single-tau multi-bin SR. The simulated sample contains 30,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{u}}) = 1.2$ TeV, $\beta = 0.5$ for the di-tau SR. The simulated sample contains 210,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the di-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{u}}) = 1.2$ TeV, $\beta = 0.5$ for the single-tau one-bin SR. The simulated sample contains 210,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{u}}) = 1.2$ TeV, $\beta = 0.5$ for the single-tau multi-bin SR. The simulated sample contains 210,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{d}}) = 1.2$ TeV, $\beta = 0.5$ for the di-tau SR. The simulated sample contains 210,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the di-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{d}}) = 1.2$ TeV, $\beta = 0.5$ for the single-tau one-bin SR. The simulated sample contains 210,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{d}}) = 1.2$ TeV, $\beta = 0.5$ for the single-tau multi-bin SR. The simulated sample contains 210,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{v}}) = 1.4$ TeV, $\beta = 0.5$ in the minimal-coupling scenario for the di-tau SR. The simulated sample contains 50,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the di-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{v}}) = 1.4$ TeV, $\beta = 0.5$ in the minimal-coupling scenario for the single-tau one-bin SR. The simulated sample contains 50,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{v}}) = 1.4$ TeV, $\beta = 0.5$ in the minimal-coupling scenario for the single-tau multi-bin SR. The simulated sample contains 50,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{v}}) = 1.4$ TeV, $\beta = 0.5$ in the Yang--Mills scenario for the di-tau SR. The simulated sample contains 50,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the di-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{v}}) = 1.4$ TeV, $\beta = 0.5$ in the Yang--Mills scenario for the single-tau one-bin SR. The simulated sample contains 50,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Cutflow for the benchmark signal model $m(\text{LQ}_{3}^{\text{v}}) = 1.4$ TeV, $\beta = 0.5$ in the Yang--Mills scenario for the single-tau multi-bin SR. The simulated sample contains 50,000 raw MC events. Weighted event yields are reported, normalized to an integrated luminosity of 139 fb$^{-1}$. 'Preselection' refers to the preselection for the single-tau channel.

Measured pseudorapidity dependence of DN/DETARAP for NSD collisions at a centre-of-mass energy of 2360 GeV.

Charged particle pseudorapidiy densities in the pseudorapidity region -0.5 to 0.5 for INEL collisions.

Charged particle pseudorapidiy densities in the pseudorapidity region -0.5 to 0.5 for NSD collisions.

Charged particle pseudorapidiy densities in the pseudorapidity region -0.5 to 0.5 for INEL>0 collisions.

Multiplicity distribution normalized to the bin width in the pseudorapidity region -0.5 to 0.5 for NSD collisions at a centre-of-mass energy of 900 GeV. Note that the statistical as well as the systematic uncertainties are strongly correlated between neighbouring points. See text of paper for details.

Multiplicity distribution normalized to the bin width in the pseudorapidity region -1.0 to 1.0 for NSD collisions at a centre-of-mass energy of 900 GeV. Note that the statistical as well as the systematic uncertainties are strongly correlated between neighbouring points. See text of paper for details.

Multiplicity distribution normalized to the bin width in the pseudorapidity region -1.3 to 1.3 for NSD collisions at a centre-of-mass energy of 900 GeV. Note that the statistical as well as the systematic uncertainties are strongly correlated between neighbouring points. See text of paper for details.

Multiplicity distribution normalized to the bin width in the pseudorapidity region -0.5 to 0.5 for INEL collisions at a centre-of-mass energy of 900 GeV. Note that the statistical as well as the systematic uncertainties are strongly correlated between neighbouring points. See text of paper for details.

Multiplicity distribution normalized to the bin width in the pseudorapidity region -1.0 to 1.0 for INEL collisions at a centre-of-mass energy of 900 GeV. Note that the statistical as well as the systematic uncertainties are strongly correlated between neighbouring points. See text of paper for details.

Multiplicity distribution normalized to the bin width in the pseudorapidity region -1.3 to 1.3 for INEL collisions at a centre-of-mass energy of 900 GeV. Note that the statistical as well as the systematic uncertainties are strongly correlated between neighbouring points. See text of paper for details.

Multiplicity distribution normalized to the bin width in the pseudorapidity region -0.5 to 0.5 for NSD collisions at a centre-of-mass energy of 2360 GeV. Note that the statistical as well as the systematic uncertainties are strongly correlated between neighbouring points. See text of paper for details.

Multiplicity distribution normalized to the bin width in the pseudorapidity region -1.0 to 1.0 for NSD collisions at a centre-of-mass energy of 2360 GeV. Note that the statistical as well as the systematic uncertainties are strongly correlated between neighbouring points. See text of paper for details.

Multiplicity distribution normalized to the bin width in the pseudorapidity region -1.3 to 1.3 for NSD collisions at a centre-of-mass energy of 2360 GeV. Note that the statistical as well as the systematic uncertainties are strongly correlated between neighbouring points. See text of paper for details.

Multiplicity distribution normalized to the bin width in the pseudorapidity region -0.5 to 0.5 for INEL collisions at a centre-of-mass energy of 2360 GeV. Note that the statistical as well as the systematic uncertainties are strongly correlated between neighbouring points. See text of paper for details.

Multiplicity distribution normalized to the bin width in the pseudorapidity region -1.0 to 1.0 for INEL collisions at a centre-of-mass energy of 2360 GeV. Note that the statistical as well as the systematic uncertainties are strongly correlated between neighbouring points. See text of paper for details.

Multiplicity distribution normalized to the bin width in the pseudorapidity region -1.3 to 1.3 for INEL collisions at a centre-of-mass energy of 2360 GeV. Note that the statistical as well as the systematic uncertainties are strongly correlated between neighbouring points. See text of paper for details.

Mean charged particle multiplicities in 3 pseudorapidity intervals in P P NSD collisions at 900 and 2360 GeV.

Mean CQ moments of the multiplicity distributions for the pseudorapidity range -0.5 to 0.5 in P P NSD collisions at centre-of-mass energies 900 and 2360 GeV.

Mean CQ moments of the multiplicity distributions for the pseudorapidity range -1.0 to 1.0 in P P NSD collisions at centre-of-mass energies 900 and 2360 GeV.

Mean CQ moments of the multiplicity distributions for the pseudorapidity range -1.3 to 1.3 in P P NSD collisions at centre-of-mass energies 900 and 2360 GeV.

The $\Lambda_{\rm {c}}^{+}$/${\rm D}^{0}$ ratio measured in p-Pb collisions at $\sqrt{s_{\rm {NN}}} = 5.02$ TeV in the rapidity interval $-0.96 \lt y \lt 0.04$ as a function of $p_{\rm {T}}$.

The $\Lambda_{\rm {c}}^{+}$/${\rm D}^{0}$ ratio measured in pp collisions at $\sqrt{s} = 7$ TeV in the transverse momentum interval $ 1 \lt p_{\rm T} \lt 8$ GeV/c and in the rapidity interval $|y|<0.5$.

The $\Lambda_{\rm {c}}^{+}$/${\rm D}^{0}$ ratio measured in p-Pb collisions at $\sqrt{s_{\rm {NN}}} = 5.02$ TeV in the transverse momentum interval $2 \lt p_{\rm T} \lt 12$ GeV/c and in the rapidity interval $-0.96 \lt y \lt 0.04$.

The nuclear modification factor $R_\mathrm{pPb}$ of prompt $\Lambda_{\rm {c}}^{+}$ baryons in p-Pb collisions at $\sqrt{s_{\rm {NN}}} = 5.02$ TeV in the rapidity interval $ -0.96\lt y \lt 0.04$.

The expected (pre-fit) and observed numbers of events for an integrated luminosity of 139 fb$^{-1}$ at $\sqrt{s}$=13 TeV in the signal region 115 < $m_{4l}$< 130 GeV and sideband region $m_{4l}$ in 105-115 GeV or 130-160 GeV (350 GeV for $tXX$-enriched) in each reconstructed event category assuming the SM Higgs boson signal with a mass $m_{H}$= 125 GeV. Combined statistical and systematic uncertainties are included for the predictions. Expected contributions that are below 0.2% of the total yield in each reconstructed event category are not shown and replaced by -.

The expected and the observed (post-fit) the jet multiplicity distribution after the inclusive event selection for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed (post-fit) the four-lepton transverse momentum distribution for events with zero jets, after the inclusive event selection for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed (post-fit) the four-lepton transverse momentum distribution for events with one jet, after the inclusive event selection for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed (post-fit) the four-lepton transverse momentum distribution for events with at least two jets, after the inclusive event selection for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed (post-fit) the dijet invariant mass distribution for events with at least two jets, after the inclusive event selection for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed $NN_{ggF}$ output (post-fit) distribution in 0$j$-$p_{T}^{4l}$-Low category for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed $NN_{ggF}$ output (post-fit) distribution in 0$j$-$p_{T}^{4l}$-Med category for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed $NN_{VBF}$ output (post-fit) distribution in 1$j$-$p_{T}^{4l}$-Low category with $NN_{ZZ}$<0.25 for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed $NN_{ZZ}$ output (post-fit) distribution in 1$j$-$p_{T}^{4l}$-Low category with $NN_{ZZ}$>0.25 for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed $NN_{VBF}$ output (post-fit) distribution in 1$j$-$p_{T}^{4l}$-Med category with $NN_{ZZ}$<0.25 for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed $NN_{ZZ}$ output (post-fit) distribution in 1$j$-$p_{T}^{4l}$-Med category with $NN_{ZZ}$>0.25 for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed $NN_{VBF}$ output (post-fit) distribution in 1$j$-$p_{T}^{4l}$-High category for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed $NN_{VBF}$ output (post-fit) distribution in 2$j$ category with $NN_{VH}$<0.2 for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed $NN_{VH}$ output (post-fit) distribution in 2$j$ category with $NN_{VH}$>0.2 for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed $NN_{VBF}$ output (post-fit) distribution in 2$j$-BSM-like category for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed $NN_{ttH}$ output (post-fit) distribution in $ttH$-Had-enriched category with $NN_{tXX}$<0.4 for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed $NN_{tXX}$ output (post-fit) distribution in $ttH$-Had-enriched category with $NN_{tXX}$>0.4 for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed $NN_{ttH}$ output (post-fit) distribution in $VH$-Lep-enriched category for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed events in the categories where no NN discriminant is used for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected and the observed events in the side bands used to constraint the $ZZ^{*}$ and $tXX$ backgrounds for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13TeV. The SM Higgs boson signal is assumed to have a mass $m_{H}$=125 GeV.

The expected SM cross section $(\sigma \cdot BR)_{SM}$, the observed value of $(\sigma \cdot BR)$, and their ratio $(\sigma \cdot BR)/(\sigma \cdot BR)_{SM}$ for the inclusive production and for each Production Mode and Reduced Stage-1.1 production bin for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13 TeV.

The correlation matrices between the measured cross-sections and the $ZZ$ and $tXX$ normalisation factors for Production Mode Stage.

The correlation matrices between the measured cross-sections and the $ZZ$ and $tXX$ normalisation factors for Reduced Stage 1.1.

Likelihood contours at 68% CLs in the (ggF,VBF) plane. The observed best fit value is the first value of the table.

Likelihood contours at 95% CLs in the (ggF,VBF) plane. The observed best fit value is the first value of the table.

Likelihood contours at 68% CLs in the (ggF,VH) plane. The observed best fit value is the first value of the table.

Likelihood contours at 95% CLs in the (ggF,VH) plane. The observed best fit value is the first value of the table.

Likelihood contours at 68% CLs in the (VBF,VH) plane. The observed best fit value is the first value of the table.

Likelihood contours at 95% CLs in the (VBF,VH) plane. The observed best fit value is the first value of the table.

Likelihood contours at 68% CLs in the ($gg2H-0j-p_T^H-Low$,$gg2H-0j-p_T^H-High$) plane. The observed best fit value is the first value of the table.

Likelihood contours at 95% CLs in the ($gg2H-0j-p_T^H-Low$,$gg2H-0j-p_T^H-High$) plane. The observed best fit value is the first value of the table.

Likelihood contours at 68% CLs in the ($\kappa_{V}$vs.$\kappa_{F}$) plane. The observed best fit value is the first value of the table.

Likelihood contours at 95% CLs in the ($\kappa_{V}$vs.$\kappa_{F}$) plane. The observed best fit value is the first value of the table.

The expected signal yield ratio for chosen CP-even and CP-odd EFT parameter values together with the corresponding cross-section measurement in each production bin of Reduced Stage 1.1. The parameter values correspond approximately to the expected confidence intervals at the 68% CLs obtained from the statistical interpretation of data.

The expected and observed confidence intervals at 68% and 95% CLs of the SMEFT Wilson coefficients for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13 TeV. For the coupling with two best fit values, the two values have been splitted in the two different columns 'Two Best-fit value - first' and 'Two Best-fit value - second', and the symbol '-' has been inserted in the missing fields.

Expected 2D-fit likelihood curves at the 95% CLs for the $c_{HW}$ versus $c_{HB}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Observed 2D-fit likelihood curves at the 95% CLs for the $c_{HW}$ versus $c_{HB}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Expected 2D-fit likelihood curves at the 95% CLs for the $c_{HB}$ versus $c_{HG}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Observed 2D-fit likelihood curves at the 95% CLs for the $c_{HB}$ versus $c_{HG}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Expected 2D-fit likelihood curves at the 95% CLs for the $c_{uH}$ versus $c_{HG}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Observed 2D-fit likelihood curves at the 95% CLs for the $c_{uH}$ versus $c_{HG}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

The best-fit values and the corresponding deviation from the SM prediction obtained from two-dimensional likelihood scans of the CP-odd BSM coupling parameters with 139fb$^{-1}$ data at $\sqrt{s}$=13 TeV. For the coupling with two best fit values, the two values have been splitted in the two different columns 'Two Observed best-fit first/second_POI - first' and 'Two Observed best-fit first/second_POI - second', and the symbol '-' has been inserted in the missing fields.

Expected 2D-fit likelihood curves at the 95% CLs for the $c_{H\tilde{W}}$ versus $c_{H\tilde{B}}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Observed 2D-fit likelihood curves at the 95% CLs for the $c_{H\tilde{W}}$ versus $c_{H\tilde{B}}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Expected 2D-fit likelihood curves at the 95% CLs for the $c_{H\tilde{B}}$ versus $c_{H\tilde{G}}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Observed 2D-fit likelihood curves at the 95% CLs for the $c_{H\tilde{B}}$ versus $c_{H\tilde{G}}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Expected 2D-fit likelihood curves at the 95% CLs for the $c_{\tilde{u}H}$ versus $c_{H\tilde{G}}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Observed 2D-fit likelihood curves at the 95% CLs for the $c_{\tilde{u}H}$ versus $c_{H\tilde{G}}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Signal acceptance times efficiency, expressed as a percentage, obtained from the ratio of the number of simulated signal events satisfying the event selection criteria in each reconstructed category to the total number of events generated in the phase space specified by a given Production model bin.

Signal acceptance times efficiency, expressed as a percentage, obtained from the ratio of the number of simulated signal events satisfying the event selection criteria in each reconstructed category to the total number of events generated in the phase space specified by a given Reduced Stage-1.1 production bin.

Observed covariance matrix for the different parameters of interest when fitting the for Production Mode Stage.

Observed covariance matrix for the different parameters of interest when fitting the for Reduced Stage 1.1.

The expected SM cross section $(\sigma\cdot BR)_{SM}$, the observed value of $(\sigma\cdot BR)$, and their ratio $(\sigma\cdot BR)/(\sigma\cdot BR)_{SM}$ for the alternate Reduced Stage-1.1 scheme for an integrated luminosity of 139fb$^{-1}$ at $\sqrt{s}$=13 TeV.

The correlation matrices between the measured cross-sections and the ZZ and $tXX$ normalisation factors for alternative Reduced Stage 1.1.

Expected 2D-fit likelihood curves at the 95% CLs for the $c_{HW}$ versus $c_{HWB}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Observed 2D-fit likelihood curves at the 95% CLs for the $c_{HW}$ versus $c_{HWB}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Expected 2D-fit likelihood curves at the 95% CLs for the $c_{HB}$ versus $c_{HWB}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Observed 2D-fit likelihood curves at the 95% CLs for the $c_{HB}$ versus $c_{HWB}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Expected 2D-fit likelihood curves at the 95% CLs for the $c_{HW}$ versus $c_{HG}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Observed 2D-fit likelihood curves at the 95% CLs for the $c_{HW}$ versus $c_{HG}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Expected 2D-fit likelihood curves at the 95% CLs for the $c_{HWB}$ versus $c_{HG}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Observed 2D-fit likelihood curves at the 95% CLs for the $c_{HWB}$ versus $c_{HG}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Expected 2D-fit likelihood curves at the 95% CLs for the $c_{H\tilde{W}}$ versus $c_{H\tilde{W}B}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Observed 2D-fit likelihood curves at the 95% CLs for the $c_{H\tilde{W}}$ versus $c_{H\tilde{W}B}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Expected 2D-fit likelihood curves at the 95% CLs for the $c_{H\tilde{B}}$ versus $c_{H\tilde{W}B}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Observed 2D-fit likelihood curves at the 95% CLs for the $c_{H\tilde{B}}$ versus $c_{H\tilde{W}B}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Expected 2D-fit likelihood curves at the 95% CLs for the $c_{H\tilde{W}}$ versus $c_{H\tilde{G}}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Observed 2D-fit likelihood curves at the 95% CLs for the $c_{H\tilde{W}}$ versus $c_{H\tilde{G}}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Expected 2D-fit likelihood curves at the 95% CLs for the $c_{H\tilde{W}B}$ versus $c_{H\tilde{G}}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Observed 2D-fit likelihood curves at the 95% CLs for the $c_{H\tilde{W}B}$ versus $c_{H\tilde{G}}$ coupling parameters at an integrated luminosity of 139 fb$^{-1}$ and $\sqrt{s}$=13 TeV. Except for the two fitted Wilson coefficients, all others are set to zero.

Summary of results for bottom muon v2 as a function of pT for different centrality. Uncertainties are statistical and systematic, respectively.

Summary of results for charm muon v3 as a function of pT for different centrality. Uncertainties are statistical and systematic, respectively.

Summary of results for bottom muon v3 as a function of pT for different centrality. Uncertainties are statistical and systematic, respectively.