Measured values of $R(\tau/e)$ in the different $p_T$ bins. The first row gives the main inclusive measurement of $R(\tau/e)$.

Breakdown of systematic uncertainties in the fit to data. Different individual components used in the fit are combined into categories. The size of the impact each group of uncertainties has on $R(\tau/e)$, which is shown in the second column, is assessed by subtracting in quadrature the uncertainty from the nominal fit and the fit with the relevant nuisance parameters fixed to their post-fit values. The group 'Other sources' includes the uncertainties due to trigger, JVT, luminosity, and pile-up.

Distributions of $|d_0|$ for the $\mu e$ channel of the nominal sample in $7 < p_T < 10$ GeV bin. The data are represented by black markers. The different contributions to the nominal sample, taken from simulation after the fit (`Post-Fit'), are shown by stacked histograms. The bottom panels show the ratio of the data to the simulation prediction after the fit. The systematic uncertainties in the prediction with the constraints from the fit and the correlation between the nuisance parameters applied are represented by the blue hatched bands.

Distributions of $|d_0|$ for the $\mu e$ channel of the nominal sample in $10 < p_T < 20$ GeV bin. The data are represented by black markers. The different contributions to the nominal sample, taken from simulation after the fit (`Post-Fit'), are shown by stacked histograms. The bottom panels show the ratio of the data to the simulation prediction after the fit. The systematic uncertainties in the prediction with the constraints from the fit and the correlation between the nuisance parameters applied are represented by the blue hatched bands.

Distributions of $|d_0|$ for the $\mu e$ channel of the nominal sample in $20 < p_T < 250$ GeV bin. The data are represented by black markers. The different contributions to the nominal sample, taken from simulation after the fit (`Post-Fit'), are shown by stacked histograms. The bottom panels show the ratio of the data to the simulation prediction after the fit. The systematic uncertainties in the prediction with the constraints from the fit and the correlation between the nuisance parameters applied are represented by the blue hatched bands.

Distributions of $|d_0|$ for the $e e$ channel of the nominal sample in $5 < p_T < 10$ GeV bin. The data are represented by black markers. The different contributions to the nominal sample, taken from simulation after the fit (`Post-Fit'), are shown by stacked histograms. The bottom panels show the ratio of the data to the simulation prediction after the fit. The systematic uncertainties in the prediction with the constraints from the fit and the correlation between the nuisance parameters applied are represented by the blue hatched bands.

Distributions of $|d_0|$ for the $e e$ channel of the nominal sample in $10 < p_T < 20$ GeV bin. The data are represented by black markers. The different contributions to the nominal sample, taken from simulation after the fit (`Post-Fit'), are shown by stacked histograms. The bottom panels show the ratio of the data to the simulation prediction after the fit. The systematic uncertainties in the prediction with the constraints from the fit and the correlation between the nuisance parameters applied are represented by the blue hatched bands.

Distributions of $|d_0|$ for the $e e$ channel of the nominal sample in $20 < p_T < 250$ GeV bin. The data are represented by black markers. The different contributions to the nominal sample, taken from simulation after the fit (`Post-Fit'), are shown by stacked histograms. The bottom panels show the ratio of the data to the simulation prediction after the fit. The systematic uncertainties in the prediction with the constraints from the fit and the correlation between the nuisance parameters applied are represented by the blue hatched bands.

Statistical covariance between bins of Table 2

Mean multiplicity of $K^{0}_{S}$ per unit $(\eta, \phi)$ in the transverse region vs. leading-jet $p_{T}$

Statistical covariance between bins of Table 3

Mean scalar sum-$p_{T}$ of $K^{0}_{S}$ per unit $(\eta, \phi)$ in the away region vs. leading-jet $p_{T}$

Statistical covariance between bins of Table 4

Mean scalar sum-$p_{T}$ of $K^{0}_{S}$ per unit $(\eta, \phi)$ in the towards region vs. leading-jet $p_{T}$

Statistical covariance between bins of Table 5

Mean scalar sum-$p_{T}$ of $K^{0}_{S}$ per unit $(\eta, \phi)$ in the transverse region vs. leading-jet $p_{T}$

Statistical covariance between bins of Table 6

Ratio of the multiplicity of $K^{0}_{S}$ to prompt charged particles in the away region vs. leading-jet $p_{T}$

Statistical covariance between bins of Table 7

Ratio of the multiplicity of $K^{0}_{S}$ to prompt charged particles in the towards region vs. leading-jet $p_{T}$

Statistical covariance between bins of Table 8

Ratio of the multiplicity of $K^{0}_{S}$ to prompt charged particles in the transverse region vs. leading-jet $p_{T}$

Statistical covariance between bins of Table 9

Ratio of the scalar sum-pt of $K^{0}_{S}$ to prompt charged particles in the away region vs. leading-jet $p_{T}$

Statistical covariance between bins of Table 10

Ratio of the scalar sum-pt of $K^{0}_{S}$ to prompt charged particles in the towards region vs. leading-jet $p_{T}$

Statistical covariance between bins of Table 11

Ratio of the scalar sum-pt of $K^{0}_{S}$ to prompt charged particles in the transverse region vs. leading-jet $p_{T}$

Statistical covariance between bins of Table 12

Mean-$p_{T}$ of $K^{0}_{S}$ in the away region vs. leading-jet $p_{T}$

Statistical covariance between bins of Table 13

Mean-$p_{T}$ of $K^{0}_{S}$ in the towards region vs. leading-jet $p_{T}$

Statistical covariance between bins of Table 14

Mean-$p_{T}$ of $K^{0}_{S}$ in the transverse region vs. leading-jet $p_{T}$

Statistical covariance between bins of Table 15

Mean multiplicity of $\Lambda$ and $\bar{\Lambda}$ per unit $(\eta, \phi)$ in the away region vs. leading-jet $p_{T}$

Statistical covariance between bins of Table 16

Mean multiplicity of $\Lambda$ and $\bar{\Lambda}$ per unit $(\eta, \phi)$ in the towards region vs. leading-jet $p_{T}$

Statistical covariance between bins of Table 17

Mean multiplicity of $\Lambda$ and $\bar{\Lambda}$ per unit $(\eta, \phi)$ in the transverse region vs. leading-jet $p_{T}$

Statistical covariance between bins of Table 18

Mean scalar sum-$p_{T}$ of $\Lambda$ and $\bar{\Lambda}$ per unit $(\eta, \phi)$ in the away region vs. leading-jet $p_{T}$

Statistical covariance between bins of Table 19

Mean scalar sum-$p_{T}$ of $\Lambda$ and $\bar{\Lambda}$ per unit $(\eta, \phi)$ in the towards region vs. leading-jet $p_{T}$

Statistical covariance between bins of Table 20

Mean scalar sum-$p_{T}$ of $\Lambda$ and $\bar{\Lambda}$ per unit $(\eta, \phi)$ in the transverse region vs. leading-jet $p_{T}$

Statistical covariance between bins of Table 21

Ratio of the multiplicity of $\Lambda$ and $\bar{\Lambda}$ to prompt charged particles in the away region vs. leading-jet $p_{T}$

Statistical covariance between bins of Table 22

Ratio of the multiplicity of $\Lambda$ and $\bar{\Lambda}$ to prompt charged particles in the towards region vs. leading-jet $p_{T}$

Statistical covariance between bins of Table 23

Ratio of the multiplicity of $\Lambda$ and $\bar{\Lambda}$ to prompt charged particles in the transverse region vs. leading-jet $p_{T}$

Statistical covariance between bins of Table 24

Ratio of the scalar sum-pt of $\Lambda$ and $\bar{\Lambda}$ to prompt charged particles in the away region vs. leading-jet $p_{T}$

Statistical covariance between bins of Table 25

Ratio of the scalar sum-pt of $\Lambda$ and $\bar{\Lambda}$ to prompt charged particles in the towards region vs. leading-jet $p_{T}$

Statistical covariance between bins of Table 26

Ratio of the scalar sum-pt of $\Lambda$ and $\bar{\Lambda}$ to prompt charged particles in the transverse region vs. leading-jet $p_{T}$

Statistical covariance between bins of Table 27

Ratio of the multiplicity of $\Lambda$ and $\bar{\Lambda}$ to $K^{0}_{S}$ in the away region vs. leading-jet $p_{T}$

Statistical covariance between bins of Table 28

Ratio of the multiplicity of $\Lambda$ and $\bar{\Lambda}$ to $K^{0}_{S}$ in the towards region vs. leading-jet $p_{T}$

Statistical covariance between bins of Table 29

Ratio of the multiplicity of $\Lambda$ and $\bar{\Lambda}$ to $K^{0}_{S}$ in the transverse region vs. leading-jet $p_{T}$

Statistical covariance between bins of Table 30

Ratio of the scalar sum-pt of $\Lambda$ and $\bar{\Lambda}$ to $K^{0}_{S}$ in the away region vs. leading-jet $p_{T}$

Statistical covariance between bins of Table 31

Ratio of the scalar sum-pt of $\Lambda$ and $\bar{\Lambda}$ to $K^{0}_{S}$ in the towards region vs. leading-jet $p_{T}$

Statistical covariance between bins of Table 32

Ratio of the scalar sum-pt of $\Lambda$ and $\bar{\Lambda}$ to $K^{0}_{S}$ in the transverse region vs. leading-jet $p_{T}$

Statistical covariance between bins of Table 33

Mean-$p_{T}$ of $\Lambda$ and $\bar{\Lambda}$ in the away region vs. leading-jet $p_{T}$

Statistical covariance between bins of Table 34

Mean-$p_{T}$ of $\Lambda$ and $\bar{\Lambda}$ in the towards region vs. leading-jet $p_{T}$

Statistical covariance between bins of Table 35

Mean-$p_{T}$ of $\Lambda$ and $\bar{\Lambda}$ in the transverse region vs. leading-jet $p_{T}$

Statistical covariance between bins of Table 36

Mean multiplicity of $K^{0}_{S}$ per unit $(\eta, \phi)$ in the away region vs. $N_\textrm{ch,trans}$

Statistical covariance between bins of Table 37

Mean multiplicity of $K^{0}_{S}$ per unit $(\eta, \phi)$ in the towards region vs. $N_\textrm{ch,trans}$

Statistical covariance between bins of Table 38

Mean multiplicity of $K^{0}_{S}$ per unit $(\eta, \phi)$ in the transverse region vs. $N_\textrm{ch,trans}$

Statistical covariance between bins of Table 39

Mean scalar sum-$p_{T}$ of $K^{0}_{S}$ per unit $(\eta, \phi)$ in the away region vs. $N_\textrm{ch,trans}$

Statistical covariance between bins of Table 40

Mean scalar sum-$p_{T}$ of $K^{0}_{S}$ per unit $(\eta, \phi)$ in the towards region vs. $N_\textrm{ch,trans}$

Statistical covariance between bins of Table 41

Mean scalar sum-$p_{T}$ of $K^{0}_{S}$ per unit $(\eta, \phi)$ in the transverse region vs. $N_\textrm{ch,trans}$

Statistical covariance between bins of Table 42

Ratio of the multiplicity of $K^{0}_{S}$ to prompt charged particles in the away region vs. $N_\textrm{ch,trans}$

Statistical covariance between bins of Table 43

Ratio of the multiplicity of $K^{0}_{S}$ to prompt charged particles in the towards region vs. $N_\textrm{ch,trans}$

Statistical covariance between bins of Table 44

Ratio of the multiplicity of $K^{0}_{S}$ to prompt charged particles in the transverse region vs. $N_\textrm{ch,trans}$

Statistical covariance between bins of Table 45

Ratio of the scalar sum-pt of $K^{0}_{S}$ to prompt charged particles in the away region vs. $N_\textrm{ch,trans}$

Statistical covariance between bins of Table 46

Ratio of the scalar sum-pt of $K^{0}_{S}$ to prompt charged particles in the towards region vs. $N_\textrm{ch,trans}$

Statistical covariance between bins of Table 47

Ratio of the scalar sum-pt of $K^{0}_{S}$ to prompt charged particles in the transverse region vs. $N_\textrm{ch,trans}$

Statistical covariance between bins of Table 48

Mean-$p_{T}$ of $K^{0}_{S}$ in the away region vs. $N_\textrm{ch,trans}$

Statistical covariance between bins of Table 49

Mean-$p_{T}$ of $K^{0}_{S}$ in the towards region vs. $N_\textrm{ch,trans}$

Statistical covariance between bins of Table 50

Mean-$p_{T}$ of $K^{0}_{S}$ in the transverse region vs. $N_\textrm{ch,trans}$

Statistical covariance between bins of Table 51

Mean multiplicity of $\Lambda$ and $\bar{\Lambda}$ per unit $(\eta, \phi)$ in the away region vs. $N_\textrm{ch,trans}$

Statistical covariance between bins of Table 52

Mean multiplicity of $\Lambda$ and $\bar{\Lambda}$ per unit $(\eta, \phi)$ in the towards region vs. $N_\textrm{ch,trans}$

Statistical covariance between bins of Table 53

Mean multiplicity of $\Lambda$ and $\bar{\Lambda}$ per unit $(\eta, \phi)$ in the transverse region vs. $N_\textrm{ch,trans}$

Statistical covariance between bins of Table 54

Mean scalar sum-$p_{T}$ of $\Lambda$ and $\bar{\Lambda}$ per unit $(\eta, \phi)$ in the away region vs. $N_\textrm{ch,trans}$

Statistical covariance between bins of Table 55

Mean scalar sum-$p_{T}$ of $\Lambda$ and $\bar{\Lambda}$ per unit $(\eta, \phi)$ in the towards region vs. $N_\textrm{ch,trans}$

Statistical covariance between bins of Table 56

Mean scalar sum-$p_{T}$ of $\Lambda$ and $\bar{\Lambda}$ per unit $(\eta, \phi)$ in the transverse region vs. $N_\textrm{ch,trans}$

Statistical covariance between bins of Table 57

Ratio of the multiplicity of $\Lambda$ and $\bar{\Lambda}$ to prompt charged particles in the away region vs. $N_\textrm{ch,trans}$

Statistical covariance between bins of Table 58

Ratio of the multiplicity of $\Lambda$ and $\bar{\Lambda}$ to prompt charged particles in the towards region vs. $N_\textrm{ch,trans}$

Statistical covariance between bins of Table 59

Ratio of the multiplicity of $\Lambda$ and $\bar{\Lambda}$ to prompt charged particles in the transverse region vs. $N_\textrm{ch,trans}$

Statistical covariance between bins of Table 60

Ratio of the scalar sum-pt of $\Lambda$ and $\bar{\Lambda}$ to prompt charged particles in the away region vs. $N_\textrm{ch,trans}$

Statistical covariance between bins of Table 61

Ratio of the scalar sum-pt of $\Lambda$ and $\bar{\Lambda}$ to prompt charged particles in the towards region vs. $N_\textrm{ch,trans}$

Statistical covariance between bins of Table 62

Ratio of the scalar sum-pt of $\Lambda$ and $\bar{\Lambda}$ to prompt charged particles in the transverse region vs. $N_\textrm{ch,trans}$

Statistical covariance between bins of Table 63

Ratio of the multiplicity of $\Lambda$ and $\bar{\Lambda}$ to $K^{0}_{S}$ in the away region vs. $N_\textrm{ch,trans}$

Statistical covariance between bins of Table 64

Ratio of the multiplicity of $\Lambda$ and $\bar{\Lambda}$ to $K^{0}_{S}$ in the towards region vs. $N_\textrm{ch,trans}$

Statistical covariance between bins of Table 65

Ratio of the multiplicity of $\Lambda$ and $\bar{\Lambda}$ to $K^{0}_{S}$ in the transverse region vs. $N_\textrm{ch,trans}$

Statistical covariance between bins of Table 66

Ratio of the scalar sum-pt of $\Lambda$ and $\bar{\Lambda}$ to $K^{0}_{S}$ in the away region vs. $N_\textrm{ch,trans}$

Statistical covariance between bins of Table 67

Ratio of the scalar sum-pt of $\Lambda$ and $\bar{\Lambda}$ to $K^{0}_{S}$ in the towards region vs. $N_\textrm{ch,trans}$

Statistical covariance between bins of Table 68

Ratio of the scalar sum-pt of $\Lambda$ and $\bar{\Lambda}$ to $K^{0}_{S}$ in the transverse region vs. $N_\textrm{ch,trans}$

Statistical covariance between bins of Table 69

Mean-$p_{T}$ of $\Lambda$ and $\bar{\Lambda}$ in the away region vs. $N_\textrm{ch,trans}$

Statistical covariance between bins of Table 70

Mean-$p_{T}$ of $\Lambda$ and $\bar{\Lambda}$ in the towards region vs. $N_\textrm{ch,trans}$

Statistical covariance between bins of Table 71

Mean-$p_{T}$ of $\Lambda$ and $\bar{\Lambda}$ in the transverse region vs. $N_\textrm{ch,trans}$

Statistical covariance between bins of Table 72

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Summary of the total uncertainty in the background prediction for each SR of the tt0L-low, tt0L-high, tt1L and tt2L analysis channels in the statistical combination. Their dominant contributions are indicated by individual lines. Individual uncertainties can be correlated, and do not necessarily add up in quadrature to the total background uncertainty.

Exclusion limits for colour-neutral scalar mediator dark matter models as a function of the mediator mass $m(\phi)$ for a DM mass $m_{\chi} = 1$ GeV. Associated production of DM with both single top quarks ($tW$ and $tj$ channels) and top quark pairs is considered. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for each individual channel and their statistical combination.

Exclusion limits for colour-neutral pseudoscalar mediator dark matter models as a function of the mediator mass $m(a)$ for a DM mass $m_{\chi} = 1$ GeV. Associated production of DM with both single top quarks ($tW$ and $tj$ channels) and top quark pairs is considered. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for each individual channel and their statistical combination.

$E_{\text{T}}^{\text{miss}}$ distribution in SR0X for the tt0L-low analysis. The contributions from all SM backgrounds are shown after the profile likelihood simultaneous fit to all tt0L-low CRs, with the hatched bands representing the total uncertainty. The category '$t\bar{t}$ (other)' represents $t\bar{t}$ events without extra jets or events with extra light-flavour jets. 'Other' includes contributions from $t\bar{t}W$, $tZ$ and $tWZ$ processes. The expected distributions for selected signal models are shown as dashed lines. The overflow events are included in the last bin. The bottom panels show the ratio of the observed data to the total SM background prediction, with the hatched area representing the total uncertainty in the background prediction and the red arrows marking data outside the vertical-axis range.

$E_{\text{T}}^{\text{miss}}$ distribution in SRWX for the tt0L-low analysis. The contributions from all SM backgrounds are shown after the profile likelihood simultaneous fit to all tt0L-low CRs, with the hatched bands representing the total uncertainty. The category '$t\bar{t}$ (other)' represents $t\bar{t}$ events without extra jets or events with extra light-flavour jets. 'Other' includes contributions from $t\bar{t}W$, $tZ$ and $tWZ$ processes. The expected distributions for selected signal models are shown as dashed lines. The overflow events are included in the last bin. The bottom panels show the ratio of the observed data to the total SM background prediction, with the hatched area representing the total uncertainty in the background prediction and the red arrows marking data outside the vertical-axis range.

$E_{\text{T}}^{\text{miss}}$ distribution in SRTX for the tt0L-low analysis. The contributions from all SM backgrounds are shown after the profile likelihood simultaneous fit to all tt0L-low CRs, with the hatched bands representing the total uncertainty. The category '$t\bar{t}$ (other)' represents $t\bar{t}$ events without extra jets or events with extra light-flavour jets. 'Other' includes contributions from $t\bar{t}W$, $tZ$ and $tWZ$ processes. The expected distributions for selected signal models are shown as dashed lines. The overflow events are included in the last bin. The bottom panels show the ratio of the observed data to the total SM background prediction, with the hatched area representing the total uncertainty in the background prediction and the red arrows marking data outside the vertical-axis range.

Exclusion limits for colour-neutral scalar mediator dark matter models as a function of the mediator mass $m(\phi)$ for a DM mass $m_{\chi} = 1$ GeV. Associated production of DM with both single top quarks ($tW$ and $tj$ channels) and top quark pairs is considered. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the nominal cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for the tt0L-high and tt0L-low analyses and their statistical combination.

Exclusion limits for colour-neutral pseudoscalar mediator dark matter models as a function of the mediator mass $m(a)$ for a DM mass $m_{\chi} = 1$ GeV. Associated production of DM with both single top quarks ($tW$ and $tj$ channels) and top quark pairs is considered. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the nominal cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for the tt0L-high and tt0L-low analyses and their statistical combination.

Exclusion limits for colour-neutral scalar mediator dark matter models as a function of the mediator mass $m(\phi)$ for a DM mass $m_{\chi} = 1$ GeV. Only associated production of DM with top quark pairs is considered for this interpretation. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for each individual channel and their statistical combination.

Exclusion limits for colour-neutral pseudoscalar mediator dark matter models as a function of the mediator mass $m(a)$ for a DM mass $m_{\chi} = 1$ GeV. Only associated production of DM with top quark pairs is considered for this interpretation. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for each individual channel and their statistical combination.

Exclusion limits for colour-neutral scalar mediator dark matter models as a function of the mediator mass $m(\phi)$ for a DM mass $m_{\chi} = 1$ GeV. Only associated production of DM with top quark pairs is considered for this interpretation. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the nominal cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for the tt0L-high and tt0L-low analyses and their statistical combination.

Exclusion limits for colour-neutral pseudoscalar mediator dark matter models as a function of the mediator mass $m(a)$ for a DM mass $m_{\chi} = 1$ GeV. Only associated production of DM with top quark pairs is considered for this interpretation. The limits are calculated at 95% CL and are expressed in terms of the ratio of the excluded cross section to the nominal cross section for a coupling assumption of $g = g_q = g_{\chi} = 1$. The solid (dashed) lines show the observed (expected) exclusion limits for the tt0L-high and tt0L-low analyses and their statistical combination.

Representative fit distribution in the different flavour leptons signal region for the tt2L analysis: each bin of such distribution, starting from the red arrow, corresponds to a single SR included in the fit. 'FNP' includes the contribution from fake/non-prompt lepton background arising from jets (mainly $\pi/K$, heavy-flavour hadron decays and photon conversion) misidentified as leptons, estimated in a purely data-driven way. 'Other' includes contributions from $t\bar{t}W$, $tZ$ and $tWZ$ processes. The total uncertainty in the SM expectation is represented with hatched bands and the expected distributions for selected signal models are shown as dashed lines.

Signal acceptance in SR0X, SRWX and SRTX for simplified DM+$t\bar{t}$ model, defined as the number of accepted events at generator level in signal Monte Carlo simulation divided by the total number of events in the sample.

Signal acceptance in SR0X, SRWX and SRTX for simplified DM+$tW$ model, defined as the number of accepted events at generator level in signal Monte Carlo simulation divided by the total number of events in the sample.

Signal acceptance in SR0X, SRWX and SRTX for simplified DM+$tj$ model, defined as the number of accepted events at generator level in signal Monte Carlo simulation divided by the total number of events in the sample.

Signal efficiency in SR0X, SRWX and SRTX for simplified DM+$t\bar{t}$ model, defined as the number of selected reconstructed events divided by the acceptance.

Signal efficiency in SR0X, SRWX and SRTX for simplified DM+$tW$ model, defined as the number of selected reconstructed events divided by the acceptance.

Signal efficiency in SR0X, SRWX and SRTX for simplified DM+$tj$ model, defined as the number of selected reconstructed events divided by the acceptance.

Cutflow for the reference point DM+$t\bar{t}$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SR0X. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 2045000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$t\bar{t}$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SRWX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 2045000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$t\bar{t}$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SRTX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 2045000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$t\bar{t}$ $m(a, \chi) = (10, 1)$ GeV in signal region SR0X. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 400000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$t\bar{t}$ $m(a, \chi) = (10, 1)$ GeV in signal region SRWX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 400000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$t\bar{t}$ $m(a, \chi) = (10, 1)$ GeV in signal region SRTX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 400000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tW$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SR0X. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 120000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tW$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SRWX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 120000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tW$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SRTX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 120000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tW$ $m(a, \chi) = (10, 1)$ GeV in signal region SR0X. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 100000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tW$ $m(a, \chi) = (10, 1)$ GeV in signal region SRWX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 100000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tW$ $m(a, \chi) = (10, 1)$ GeV in signal region SRTX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 100000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tj$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SR0X. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 169000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tj$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SRWX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 169000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tj$ $m(\phi, \chi) = (10, 1)$ GeV in signal region SRTX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 169000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tj$ $m(a, \chi) = (10, 1)$ GeV in signal region SR0X. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 140000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tj$ $m(a, \chi) = (10, 1)$ GeV in signal region SRWX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 140000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Cutflow for the reference point DM+$tj$ $m(a, \chi) = (10, 1)$ GeV in signal region SRTX. The column labelled 'weighted' shows the event yield including all correction factors applied to simulation, and is normalised to 139 fb$^{-1}$. A notable exception concerns the 'weighted' numbers in the first and the second row, labelled 'Total' and 'Filtered', which correspond to $\mathcal{L}\cdot\sigma$ and $\mathcal{L}\cdot\sigma\cdot\epsilon$ expected, respectively. The 'Skim' selection requires the $p_{\text{T}}$ of the leading four jets to be above (80, 60, 40, 40) GeV, the missing transverse momentum $E_{\text{T}}^{\text{miss}} > 140$ GeV, the missing momentum significance $\mathcal{S} > 8$, $\Delta\phi_{\min}(\vec{p}_{\text{T,1-4}},\vec{p}_{\text{T}}^{\text{miss}}) > 0.4$ and a lepton veto. The 'Orthogonalisation' selection is defined in the main body. In total 140000 raw MC events were generated prior to the specified cuts, with the column 'Unweighted yield' collecting the numbers after each cut.

Distribution of ${E}_{T}^{miss}$ after the fit of the multijet backgrounds, in the electron channel, in the $t\bar{t}$ 3-jets VR, without applying the cut on ${E}_{T}^{miss}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Distribution of ${E}_{T}^{miss}$ after the fit of the multijet backgrounds, in the electron channel, in the $t\bar{t}$ 4-jets VR, without applying the cut on ${E}_{T}^{miss}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Distribution of $m_{T}^{W}$ after the fit of the multijet backgrounds, in the muon channel, in the signal region, without applying the cut on $m_{T}^{W}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Distribution of $m_{T}^{W}$ after the fit of the multijet backgrounds, in the muon channel, in the $W$+jets VR, without applying the cut on $m_{T}^{W}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Distribution of $m_{T}^{W}$ after the fit of the multijet backgrounds, in the muon channel, in the $t\bar{t}$ 3-jets VR, without applying the cut on $m_{T}^{W}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Distribution of $m_{T}^{W}$ after the fit of the multijet backgrounds, in the muon channel, in the $t\bar{t}$ 4-jets VR, without applying the cut on $m_{T}^{W}$. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit. The last bin includes the overflow. The uncertainty band indicates the simulation's statistical uncertainty, the normalisation uncertainties for different processes ($40$ % for $W$+jets production, $30$ % for multijet background and $6$ % for top-quark processes) and the multijet background shape uncertainty in each bin, summed in quadrature. The lower panel of the figure shows the ratio of the data to the prediction.

Expected distributions of the MEM discriminant $P(S|X)$ in the SR, for the s-channel single-top signal, and for the $t\bar{t}$ and $W$+jets backgrounds, for MEM discriminant values larger than $2.0\times10^{-4}$. Each distribution is normalised to unity. The binning is the same as the optimised binning used in the signal extraction fit, resulting in a non-linear horizontal scale.

Distribution of the MEM discriminant $P(S|X)$ in the $W$+jets VR. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit presented in Section 5 in the paper. The uncertainty band indicates the simulation's statistical uncertainty and the normalisation uncertainties for the various processes in each bin, summed in quadrature. The ratio of the observed number to the predicted number of events in each bin is shown in the lower panel of the figure, with different vertical axis ranges. The binning is the same as the optimised binning used in the signal extraction fit described in Section 8 in the paper, resulting in a non-linear horizontal scale.

Distribution of the MEM discriminant $P(S|X)$ in the $t\bar{t}$ 3-jets VR. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit presented in Section 5 in the paper. The uncertainty band indicates the simulation's statistical uncertainty and the normalisation uncertainties for the various processes in each bin, summed in quadrature. The ratio of the observed number to the predicted number of events in each bin is shown in the lower panel of the figure, with different vertical axis ranges. The binning is the same as the optimised binning used in the signal extraction fit described in Section 8 in the paper, resulting in a non-linear horizontal scale.

Distribution of the MEM discriminant $P(S|X)$ in the $t\bar{t}$ 4-jets VR. Simulated events are normalised to the expected number of events given the integrated luminosity, after applying the normalisation factors obtained in the multijet fit presented in Section 5 in the paper. The uncertainty bands indicate the simulation's statistical uncertainty and the normalisation uncertainties for the various processes in each bin, summed in quadrature. The ratio of the observed number to the predicted number of events in each bin is shown in the lower panel of the figure, with different vertical axis ranges. The binning is the same as the optimised binning used in the signal extraction fit described in Section 8 in the paper, resulting in a non-linear horizontal scale.

Distribution of the MEM discriminant $P(S|X)$ in the SR before the fit to data, for MEM discriminant values larger than $2.0\times10^{-4}$. The lower panel of the figure shows the ratio of the data to the prediction, with different vertical axis ranges. The uncertainty band indicates the total uncertainties and their correlations in each bin. The uncertainties in the $t\bar{t}$ and $W$+jets normalisation factors, as well as in the s-channel signal cross-section, are not defined pre-fit and therefore not included. The binning is the same as the optimised binning used in the fit, resulting in a non-linear horizontal scale.

Distribution of the MEM discriminant $P(S|X)$ in the SR after the fit to data, for MEM discriminant values larger than $2.0\times10^{-4}$. The lower panel of the figure shows the ratio of the data to the prediction, with different vertical axis ranges. The uncertainty band indicates the total uncertainties and their correlations in each bin. The binning is the same as the optimised binning used in the fit, resulting in a non-linear horizontal scale.

Distribution of the MEM discriminant $P(S|X)$ in the SR after the fit to data, for MEM discriminant values larger than $2.0\times10^{-4}$, after subtraction of all backgrounds. The fitted distribution for the simulation of the signal is shown together with the post-fit uncertainty in the backgrounds. The binning is the same as the optimised binning used in the fit, resulting in a non-linear horizontal scale.

Pre-fit and post-fit event yields in the SR, for MEM discriminant values larger than $2.0\times10^{-4}$. The central value of the event yield for each process is calculated by summing the values of the discriminant bin contents, using the nominal expected yield for the pre-fit value, and the best-fit estimate for the post-fit value. The error includes statistical and systematic uncertainties summed in quadrature. All sources of systematic uncertainties are included, taking into account correlations and anti-correlations in the post-fit case. The uncertainties in the $t\bar{t}$ and $W$+jets normalisation factors, as well as in the s-channel signal cross-section, are not defined pre-fit and therefore only included in the post-fit uncertainties.

Observed impact of the different sources of uncertainty on the measured s-channel signal cross-section, grouped by categories. The impact of each category is obtained by repeating the fit after having fixed the set of nuisance parameters corresponding to that category, subtracting the square of the resulting uncertainty from the square of the uncertainty found in the full fit, and calculating the square root. The 'Systematic uncertainties' category combines all sources of systematic uncertainties. The statistical uncertainty is obtained by repeating the fit after having fixed all nuisance parameters, including the $t\bar{t}$ and $W$+jets normalisation factors. 'Total' gives the total uncertainty on the measurement.

Observed impact of the different sources of $t\bar{t}$ modelling uncertainty on the measured s-channel signal cross-section. The impact of each category is obtained by repeating the fit after having fixed the set of nuisance parameters corresponding to that category, subtracting the square of the resulting uncertainty from the square of the uncertainty found in the full fit, and calculating the square root. 'PS & had.' refers to the parton shower and hadronisation model, and 'ME/PS matching' to the matching of the ME to the parton shower.

Observed impact of the different sources of s-channel modelling uncertainty on the measured s-channel signal cross-section. The impact of each category is obtained by repeating the fit after having fixed the set of nuisance parameters corresponding to that category, subtracting the square of the resulting uncertainty from the square of the uncertainty found in the full fit, and calculating the square root. 'PS & had.' refers to the parton shower and hadronisation model, as described in Section 7 in the paper.

Observed impact of the different sources of t-channel modelling uncertainty on the measured s-channel signal cross-section. The impact of each category is obtained by repeating the fit after having fixed the set of nuisance parameters corresponding to that category, subtracting the square of the resulting uncertainty from the square of the uncertainty found in the full fit, and calculating the square root. 'PS & had.' refers to the parton shower and hadronisation model, as described in Section 7 in the paper.

Observed impact of the different sources of $tW$ modelling uncertainty on the measured s-channel signal cross-section, grouped by categories. The impact of each category is obtained by repeating the fit after having fixed the set of nuisance parameters corresponding to that category, subtracting the square of the resulting uncertainty from the square of the uncertainty found in the full fit, and calculating the square root. 'PS & had.' refers to the parton shower and hadronisation model, and '$t\bar{t}$ overlap' to the algorithm removing the overlap between $tW$ and $t\bar{t}$ production at NLO, as described in Section 7 in the paper.

Observed impact of the different sources of PDF uncertainties on the measured s-channel signal cross-section, grouped by categories. The impact of each category is obtained by repeating the fit after having fixed the set of nuisance parameters corresponding to that category, subtracting the square of the resulting uncertainty from the square of the uncertainty found in the full fit, and calculating the square root.

Comparison between data and prediction after the fit to data in the signal region for the leading-jet $p_{T}$. The last bin includes the overflow. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the leading-jet $\eta$. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the subleading-jet $p_{T}$. The last bin includes the overflow. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the subleading-jet $\eta$. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the lepton $p_{T}$. The last bin includes the overflow. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the lepton $\eta$. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the ${E}_{T}^{miss}$. The last bin includes the overflow. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Comparison between data and prediction after the fit to data in the signal region for the $m_{T}^{W}$. The last bin includes the overflow. The uncertainty band includes all uncertainties and their correlations. The lower panel of the figure shows the ratio of the data to the prediction.

Nuisance parameters ranked according to their post-fit impacts on the best-fit value of the ratio $\mu$ of the measured cross-section to the predicted cross-section. In the figure, only the 20 nuisance parameters with the largest post-fit impacts are shown. The empty (solid) blue rectangles illustrate the pre-fit (post-fit) impact on $\mu$, corresponding to the upper axis. The pre-fit (post-fit) impact of each nuisance parameter, $\Delta\mu$, is calculated as the difference in the fitted value of $\mu$ between the nominal fit and the fit when fixing the corresponding nuisance parameter to $\hat{\theta}\pm\Delta\theta$ ($\hat{\theta}\pm\Delta\hat{\theta}$), where $\hat{\theta}$ is the best-fit value of the nuisance parameter and $\Delta\theta$ ($\Delta\hat{\theta}$) is its pre-fit (post-fit) uncertainty. Several systematic uncertainties are split into different nuisance parameters, which are indicated by NP. JES (JER) indicates jet energy scale (resolution), and $\gamma$ indicates a nuisance parameter associated to the MC statistics in one of the 18 bins numbered from 0 to 17. The black points show the best-fit values of the nuisance parameters, with the error bars representing the post-fit uncertainties. Each nuisance parameter is shown wrt. its nominal value, $\theta_0$, and in units of its pre-fit uncertainty, except the free-floating normalisation factors of the $t\bar{t}$ and $W$+jets backgrounds, and the parameters associated to the MC statistics in each bin, for which the post-fit values and uncertainties are shown.

Correlation matrix of the nuisance parameters and of the ratio $\mu$ of the measured cross-section to the predicted cross-section. The correlations are given after the fit to data. In the figure, only the parameters which have a correlation of at least 0.2 with any other parameter are shown.

Distribution of the MEM discriminant $P(S|X)$ in the SR for MEM discriminant values larger than $2.0\times10^{-4}$, for the collision data used for the measurement, and for 1000 pseudo-data replicas, generated using a bootstrapping technique, in order to assess the statistical correlations between this measurement and others, for the purpose of combinations. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the <a href="https://zenodo.org/record/5361038">BootstrapGenerator</a> software package , which implements a technique described in <a href="https://cds.cern.ch/record/2759945/">ATL-PHYS-PUB-2021-011</a>. The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. Each pseudo-data replica is assigned an index, ranging from 0 to 999, corresponding to the random number index used consistently for each observed data event.

Measured values of the signal cross-section and of the $t\bar{t}$ and $W$+jets normalisation factors, obtained by statistical-only fits to the collision data used for the measurement, and to 1000 pseudo-data replicas, generated using a bootstrapping technique, in order to assess the statistical correlations between this measurement and others, for the purpose of combinations. The central values and their statistical uncertainties are obtained by repeating the fit after having fixed all nuisance parameters, except the $t\bar{t}$ and $W$+jets normalisation factors, which are let free-floating (unlike for the statistical uncertainty on the cross-section quoted in the paper). The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the <a href="https://zenodo.org/record/5361038">BootstrapGenerator</a> software package , which implements a technique described in <a href="https://cds.cern.ch/record/2759945/">ATL-PHYS-PUB-2021-011</a>. The ATLAS event number and run number of each event are used as seed to uniquely but reproducibly initialise the random number generator for each event. Each pseudo-data replica is assigned an index, ranging from 0 to 999, corresponding to the random number index used consistently for each observed data event.

Predicted ratio of cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $0.6<|\eta^{\gamma}|<1.37$.

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

Predicted ratio of cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.56<|\eta^{\gamma}|<1.81$.

Measured ratio of cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.81<|\eta^{\gamma}|<2.37$.

Predicted ratio of cross sections for inclusive isolated-photon production as a function of $E_{\rm T}^{\gamma}$ for $1.81<|\eta^{\gamma}|<2.37$.

Measured double ratio of cross sections for inclusive isolated-photon production and Z-boson production as a function of $E_{\rm T}^{\gamma}$ for $|\eta^{\gamma}|<0.6$.

Predicted double ratio of cross sections for inclusive isolated-photon production and Z-boson production as a function of $E_{\rm T}^{\gamma}$ for $|\eta^{\gamma}|<0.6$.

Measured double ratio of cross sections for inclusive isolated-photon production and Z-boson production as a function of $E_{\rm T}^{\gamma}$ for $0.6<|\eta^{\gamma}|<1.37$.

Predicted double ratio of cross sections for inclusive isolated-photon production and Z-boson production as a function of $E_{\rm T}^{\gamma}$ for $0.6<|\eta^{\gamma}|<1.37$.

Measured double ratio of cross sections for inclusive isolated-photon production and Z-boson production as a function of $E_{\rm T}^{\gamma}$ for $1.56<|\eta^{\gamma}|<1.81$.

Predicted double ratio of cross sections for inclusive isolated-photon production and Z-boson production as a function of $E_{\rm T}^{\gamma}$ for $1.56<|\eta^{\gamma}|<1.81$.

Measured double ratio of cross sections for inclusive isolated-photon production and Z-boson production as a function of $E_{\rm T}^{\gamma}$ for $1.81<|\eta^{\gamma}|<2.37$.

Predicted double ratio of cross sections for inclusive isolated-photon production and Z-boson production as a function of $E_{\rm T}^{\gamma}$ for $1.81<|\eta^{\gamma}|<2.37$.

Acceptance x efficiency versus resonance mass for both lephad and hadhad channels in the RS bulk model with k/MPl = 1

Acceptance x efficiency versus resonance mass for both lephad and hadhad channels in the RS bulk model with k/MPl = 2

Acceptance x efficiency versus resonance mass for both lephad and hadhad channels in the scalar model

Upper limits on the production cross-section times the HH to bbtautau branching ratio for non-resonant HH at 95% CLS and their interpretation as multiples of the SM prediction

Upper limits on the production cross-section times the HH to bbtautau branching ratio divided by the SM prediction for non-resonant HH at 95% CL

Post-fit expected number of signal and background events and observed number of data events after applying the selection criteria and requiring exactly 2 b-tagged jets and assuming a background-only hypothesis

Post-fit expected number of signal and background events and observed number of data events in the last two bins of the non-resonant BDT score distribution of the SM signal after applying the selection criteria and requiring exactly 2 b-tagged jets and assuming a background-only hypothesis

Dijet mass spectra after the background only fit with the background prediction in the high-mass region with two b-tags.

Dijet mass spectra after the background only fit with the background prediction in the low-mass region with two b-tags.

The online b-tagging efficiency with respect to the offline b-tagging efficiency as a function of pT. The b-tagging online and offline working points correspond to an efficiency of 60% and 70%, respectively.

Observed and expected 95% credibility-level upper limits on the cross-section for the b* model in the high-mass region with inclusive b-jet selection.

Observed and expected 95% credibility-level upper limits on the cross-section times branching ratio for the SSM and leptophobic Z' models in the low- and high-mass region with two b-tags selection.

Observed and expected 95% credibility-level upper limits on the cross-section for DM Z' models in the low-mass region with two b-tags selection. The Z' is expected to decay to all five quark flavors other than the top quark and the mediator to SM quark coupling (gSM) equal to 0.1 is assumed.

Observed and expected 95% credibility-level upper limits on the cross-section times branching ratio for DM Z'->bb models in the high-mass region with two b-tags selection. The Z' is expected to decay to bb only and the mediator to SM quark coupling (gSM) equal to 0.25 is assumed.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances with intrinsic width smaller than the detector resolution. The width of the Gaussian reconstructed shape is dominated by the dijet mass resolution. The table shows the limits obtained from the high-mass inclusive one b-tag selection.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances with intrinsic width smaller than the detector resolution. The width of the Gaussian reconstructed shape is dominated by the dijet mass resolution. The table shows the limits obtained from the combined low- and high-mass two b-tags selection.

The mass distributions for the inclusive one b-tagged selection and two b-tagged selection using an integrated luminosity of 36.1 $fb^{-1}$. The inclusive one b-tagged Pythia8 MC distribution is normalized to the inclusive one b-tagged data. The two b-tagged Pythia8 MC distribution is normalized to the two b-tagged data. The systematic uncertainty band is dominated by the b-tagging scale factor and the b-jet energy scale.

Signal acceptance times efficiency for inclusive 1 b-tag and 2 b-tag categories as a function of the simulated b* and the Z' masses.

Signal acceptance times efficiency for 2 b-tag categories as a function of the simulated Z' masses.

The flavor composition of the simulated dijet background as a function of dijet mass before tagging. The fraction of the six combinations of the b-jet , c-jet and light-flavor jet are shown. All offline selections are applied.

The flavor composition of the simulated dijet background as a function of dijet mass with inclusive one b-tag. The fraction of the six combinations of the b-jet , c-jet and light-flavor jet are shown. All offline selections are applied.

The flavor composition of the simulated dijet background as a function of dijet mass with two b-tags. The fraction of the six combinations of the b-jet , c-jet and light-flavor jet are shown. All offline selections are applied.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances exhibiting a generic Gaussian shape at particle level. The table shows the limits obtained from the inclusive b-jet selection. The limits corresponding to Gaussian-shaped resonances with width of Γ(X)/m(X) = 3%.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances exhibiting a generic Gaussian shape at particle level. The table shows the limits obtained from the inclusive b-jet selection. The limits corresponding to Gaussian-shaped resonances with width of Γ(X)/m(X) = 7%.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances exhibiting a generic Gaussian shape at particle level. The table shows the limits obtained from the inclusive b-jet selection. The limits corresponding to Gaussian-shaped resonances with width of Γ(X)/m(X) = 10%.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances exhibiting a generic Gaussian shape at particle level. The table shows the limits obtained from the inclusive b-jet selection. The limits corresponding to Gaussian-shaped resonances with width of Γ(X)/m(X) = 15%.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances exhibiting a generic Gaussian shape at particle level. The table shows the limits obtained from the combined low- and high-mass two b-tags selection. The limits corresponding to Gaussian-shaped resonances with width of Γ(X)/m(X) = 3%.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances exhibiting a generic Gaussian shape at particle level. The table shows the limits obtained from the combined low- and high-mass two b-tags selection. The limits corresponding to Gaussian-shaped resonances with width of Γ(X)/m(X) = 7%.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances exhibiting a generic Gaussian shape at particle level. The table shows the limits obtained from the combined low- and high-mass two b-tags selection. The limits corresponding to Gaussian-shaped resonances with width of Γ(X)/m(X) = 10%.

Observed and expected 95% credibility-level upper limits on cross section times acceptance times branching ratio of X --> bb, including kinematic acceptance and b-tagging efficiencies, for resonances exhibiting a generic Gaussian shape at particle level. The table shows the limits obtained from the combined low- and high-mass two b-tags selection. The limits corresponding to Gaussian-shaped resonances with width of Γ(X)/m(X) = 15%.

Event yields for the different processes estimated with the fit to the $O_\mathrm{NN}$ distribution compared to the numbers of observed events. Only the statistical uncertainties are quoted. The $Z,VV+\mathrm{jets}$ contributions and the multijet background are fixed in the fit; therefore no uncertainty is quoted for these processes.

Detailed list of the contribution from each source of uncertainty to the total uncertainty in the measured values of $\sigma_{\mathrm{fid}}(tq)$ and $\sigma_{\mathrm{fid}}(\bar tq)$. The estimation of the systematic uncertainties has a statistical uncertainty of $0.3\%$. Uncertainties contributing less than $0.5\%$ are marked with ‘<0.5’.

Significant contributions to the total relative uncertainty in the measured value of $R_{t}$. The estimation of the systematic uncertainties has a statistical uncertainty of $0.3~\%$. Uncertainties contributing less than $0.5~\%$ are not shown.

Slopes $a$ of the mass dependence of the measured cross$-$sections.

Predicted (post-fit) and observed event yields for the signal region (SR), after the requirement on the neural network discriminant, $O_{\mathrm{NN}}~>~0.8$. The multijet background prediction is obtained from the fit to the $E_{\mathrm{T}}^{\mathrm{miss}}$ distribution described in Section 6, while all the other predictions and uncertainties are derived from the total cross$-$section measurement. In some cases there is no uncertainty quoted. In these cases the uncertainty is < 0.5.

Predicted (post-fit) and observed event yields for the signal region (SR), after the requirement on the second neural network discriminant, $O_{\mathrm{NN2}}~>~0.8$. The multijet background prediction is obtained from the fit to the $E_{\mathrm{T}}^{\mathrm{miss}}$ distribution described in Section 6, while all the other predictions and uncertainties are derived from the total cross$-$section measurement. In some cases there is no uncertainty quoted. In these cases the uncertainty is < 0.5.

Migration matrix for $p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})$ at the particle level. The pseudo top quark is shown on the $y$-axis and the reconstructed variable is shown on the $x$-axis.

Migration matrix for $p_{\mathrm{T}}(t)$ at the parton level. The parton-level quark is shown on the $y$-axis and the reconstructed variable is shown on the $x$-axis.

Migration matrix for $|y(\hat{t\hspace{-0.2mm}})|$ at the particle level. The pseudo top quark is shown on the $y$-axis and the reconstructed variable is shown on the $x$-axis.

Migration matrix for $|y(t)|$ at the parton level. The parton-level quark is shown on the $y$-axis and the reconstructed variable is shown on the $x$-axis.

Uncertainties in the normalisations of the different backgrounds for all processes, as derived from the total cross-section measurement.

Absolute and normalised unfolded differential $tq$ production cross$-$section as a function of $p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})$ at particle level.

Absolute and normalised unfolded differential $\bar tq$ production cross$-$section as a function of $p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})$ at particle level.

Absolute and normalised unfolded differential $tq$ production cross$-$section as a function of $|y(\hat{t\hspace{-0.2mm}})|$ at particle level.

Absolute and normalised unfolded differential $\bar tq$ production cross$-$section as a function of $|y(\hat{t\hspace{-0.2mm}})|$ at particle level.

Absolute and normalised unfolded differential $tq$ production cross$-$section as a function of $p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})$ at particle level.

Absolute and normalised unfolded differential $\bar tq$ production cross$-$section as a function of $p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})$ at particle level.

Absolute and normalised unfolded differential $tq$ production cross$-$section as a function of $|y(\hat{j\hspace{-0.2mm}})|$ at particle level.

Absolute and normalised unfolded differential $\bar tq$ production cross$-$section as a function of $|y(\hat{j\hspace{-0.2mm}})|$ at particle level.

Absolute and normalised unfolded differential $tq$ production cross$-$section as a function of $p_{\mathrm{T}}(t)$ at parton level.

Absolute and normalised unfolded differential $\bar tq$ production cross$-$section as a function of $p_{\mathrm{T}}(t)$ at parton level.

Absolute and normalised unfolded differential $tq$ production cross$-$section as a function of $|y(t)|$ at parton level.

Absolute and normalised unfolded differential $\bar tq$ production cross$-$section as a function of $|y(t)|$ at parton level.

Statistical correlation matrix for the absolute differential cross-section as a function of $p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})$ for $tq$ events(at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the absolute differential cross-section as a function of $p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})$ for $ \bar tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the normalised differential cross-section as a function of $p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})$ for $tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the normalised differential cross-section as a function of $p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})$ for $\bar tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the absolute differential cross-section as a function of $|y(\hat{t\hspace{-0.2mm}})|$ for $tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the absolute differential cross-section as a function of $|y(\hat{t\hspace{-0.2mm}})|$ for $ \bar tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the normalised differential cross-section as a function of $|y(\hat{t\hspace{-0.2mm}})|$ for $tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the normalised differential cross-section as a function of $|y(\hat{t\hspace{-0.2mm}})|$ for $\bar tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the absolute differential cross-section as a function of $p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})$ for $tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the absolute differential cross-section as a function of $p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})$ for $\bar tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the normalised differential cross-section as a function of $p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})$ for $tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the normalised differential cross-section as a function of $p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})$ for $\bar tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the absolute differential cross-section as a function of $|y(\hat{j\hspace{-0.2mm}})|$ for $tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the absolute differential cross-section as a function of $|y(\hat{j\hspace{-0.2mm}})|$ for $\bar tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the normalised differential cross-section as a function of $|y(\hat{j\hspace{-0.2mm}})|$ for $tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the normalised differential cross-section as a function of $|y(\hat{j\hspace{-0.2mm}})|$ for $\bar tq$ events (at the particle level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the absolute differential cross-section as a function of $p_{\mathrm{T}}(t)$ for $tq$ events (at the parton level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the absolute differential cross-section as a function of $p_{\mathrm{T}}(t)$ for $ \bar tq$ events (at the parton level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the normalised differential cross-section as a function of $p_{\mathrm{T}}(t)$ for $tq$ events (at the parton level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the normalised differential cross-section as a function of $p_{\mathrm{T}}(t)$ for $ \bar tq$ events (at the parton level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the absolute differential cross-section as a function of $|y(t)|$ for $tq$ events (at the parton level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the absolute differential cross-section as a function of $|y(t)|$ for $\bar tq$ events (at the parton level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the normalised differential cross-section as a function of $|y(t)|$ for $tq$ events (at the parton level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Statistical correlation matrix for the normalised differential cross-section as a function of $|y(t)|$ for $\bar tq$ events (at the parton level). It includes the statistical uncertainty due to the number of data events and MC statistics.

Fiducial acceptance $A_{\mathrm{fid}}$ for different $t$-channel single top-quark MC samples. $^{\mathrm{(a)}}$ Calculation taken from AcerMC $+$ $\mathrm{P{\scriptsize YTHIA}6}$. $^{\mathrm{(b)}}$ Calculation taken from $\mathrm{P{\scriptsize OWHEG}}$-$\mathrm{B{\scriptsize OX}}$ $+$ $\mathrm{P{\scriptsize YTHIA}6}$.

Fiducial acceptance $A_{\mathrm{fid}}$ for different $t$-channel single top-antiquark MC samples. $^{\mathrm{(a)}}$ Calculation taken from AcerMC $+$ $\mathrm{P{\scriptsize YTHIA}6}$. $^{\mathrm{(b)}}$ Calculation taken from $\mathrm{P{\scriptsize OWHEG}}$-$\mathrm{B{\scriptsize OX}}$ $+$ $\mathrm{P{\scriptsize YTHIA}6}$.

Uncertainties for the absolute differential $tq$ cross-section as a function of $p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})$ at particle level per bin ([0,35,50,75,100,150,200,300] GeV) in percent of $\dfrac{\mathrm{d}\sigma(tq)}{\mathrm{d}p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the normalised differential $tq$ cross-section as a function of $p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})$ at particle level per bin ([0,35,50,75,100,150,200,300] GeV) in percent of $\left( \dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(tq)}{\mathrm{d}p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the absolute differential $\bar tq$ cross-section as a function of $p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})$ at particle level per bin ([0,35,50,75,100,150,300] GeV) in percent of $\dfrac{\mathrm{d}\sigma(\bar tq)}{\mathrm{d}p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})}$.

Uncertainties for the normalised differential $\bar tq$ cross-section as a function of $p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})$ at particle level per bin ([0,35,50,75,100,150,300] GeV) in percent of $\left( \dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(\bar tq)}{\mathrm{d}p_{\mathrm{T}}(\hat{t\hspace{-0.2mm}})}$.

Uncertainties for the absolute differential $tq$ cross-section as a function of $|y(\hat{t\hspace{-0.2mm}})|$ at particle level per bin ([0,0.15,0.3,0.45,0.7,1.0,1.3,2.2]) in percent of $\dfrac{\mathrm{d}\sigma(tq)}{\mathrm{d}|y(\hat{t\hspace{-0.2mm}})|}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the normalised differential $tq$ cross-section as a function of $|y(\hat{t\hspace{-0.2mm}})|$ at particle level per bin ([0,0.15,0.3,0.45,0.7,1.0,1.3,2.2]) in percent of $\left(\dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(tq)}{\mathrm{d}|y(\hat{t\hspace{-0.2mm}})|}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the absolute differential $\bar tq$ cross-section as a function of $|y(\hat{t\hspace{-0.2mm}})|$ at particle level per bin ([0,0.15,0.3,0.45,0.7,1.0,1.3,2.2]) in percent of $\dfrac{\mathrm{d}\sigma(\bar tq)}{\mathrm{d}|y(\hat{t\hspace{-0.2mm}})|}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the normalised differential $\bar tq$ cross-section as a function of $|y(\hat{t\hspace{-0.2mm}})|$ at particle level per bin ([0,0.15,0.3,0.45,0.7,1.0,1.3,2.2]) in percent of $\left(\dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(\bar tq)}{\mathrm{d}|y(\hat{t\hspace{-0.2mm}})|}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the absolute differential $tq$ cross-section as a function of $p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})$ at particle level per bin ([30,45,60,75,100,150,300] GeV) in percent of $\dfrac{\mathrm{d}\sigma(tq)}{\mathrm{d}p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the normalised differential $tq$ cross-section as a function of $p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})$ at particle level per bin ([30,45,60,75,100,150,300] GeV) in percent of $\left(\dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(tq)}{\mathrm{d}p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the absolute differential $\bar tq$ cross-section as a function of $p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})$ at particle level per bin ([30,45,60,75,100,150,300] GeV) in percent of $\dfrac{\mathrm{d}\sigma(\bar tq)}{\mathrm{d}p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the normalised differential $\bar tq$ cross-section as a function of $p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})$ at particle level per bin ([30,45,60,75,100,150,300] GeV) in percent of $\left(\dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(\bar tq)}{\mathrm{d}p_{\mathrm{T}}(\hat{j\hspace{-0.2mm}})}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the absolute differential $tq$ cross-section as a function of $|y(\hat{j\hspace{-0.2mm}})|$ at particle level per bin ([0.0, 1.2, 1.7, 2.2, 2.7, 3.3, 4.5]) in percent of $\dfrac{\mathrm{d}\sigma(tq)}{\mathrm{d}|y(\hat{j\hspace{-0.2mm}})|}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the normalised differential $tq$ cross-section as a function of $|y(\hat{j\hspace{-0.2mm}})|$ at particle level per bin ([0.0, 1.2, 1.7, 2.2, 2.7, 3.3, 4.5]) in percent of $\left(\dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(tq)}{\mathrm{d}|y(\hat{j\hspace{-0.2mm}})|}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the absolute differential $\bar tq$ cross-section as a function of $|y(\hat{j\hspace{-0.2mm}})|$ at particle level per bin ([0.0, 1.2, 1.7, 2.2, 2.7, 3.3, 4.5]) in percent of $\dfrac{\mathrm{d}\sigma(\bar tq)}{\mathrm{d}|y(\hat{j\hspace{-0.2mm}})|}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the normalised differential $\bar tq$ cross-section as a function of $|y(\hat{j\hspace{-0.2mm}})|$ at particle level per bin ([0.0, 1.2, 1.7, 2.2, 2.7, 3.3, 4.5]) in percent of $\left(\dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(\bar tq)}{\mathrm{d}|y(\hat{j\hspace{-0.2mm}})|}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the absolute differential $tq$ cross-section as a function of $p_{\mathrm{T}}(t)$ at parton level per bin ([0,50,100,150,200,300] GeV) in percent of $\dfrac{\mathrm{d}\sigma(tq)}{\mathrm{d}p_{\mathrm{T}}(t)}$.

Uncertainties for the normalised differential $tq$ cross-section as a function of $p_{\mathrm{T}}(t)$ at parton level per bin ([0,50,100,150,200,300] GeV) in percent of $\left(\dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(tq)}{\mathrm{d}p_{\mathrm{T}}(t)}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the absolute differential $\bar tq $ cross-section as a function of $p_{\mathrm{T}}(t)$ at parton level per bin ([0,50,100,150,300] GeV) in percent of $\dfrac{\mathrm{d}\sigma(\bar tq)}{\mathrm{d}p_{\mathrm{T}}(t)}$.

Uncertainties for the normalised differential $\bar tq $ cross-section as a function of $p_{\mathrm{T}}(t)$ at parton level per bin ([0,50,100,150,300] GeV) in percent of $\left(\dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(\bar tq)}{\mathrm{d}p_{\mathrm{T}}(t)}$.

Uncertainties for the absolute differential $ tq $ cross-section as a function of $|y(t)|$ at parton level per bin ([0,0.3,0.7,1.3,2.2]) in percent of $\dfrac{\mathrm{d}\sigma(tq)}{\mathrm{d}|y(t)|}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the normalised differential $ tq $ cross-section as a function of $|y(t)|$ at parton level per bin ([0,0.3,0.7,1.3,2.2]) in percent of $\left(\dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(tq)}{\mathrm{d}|y(t)|}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.

Uncertainties for the absolute differential $ \bar tq $ cross-section as a function of $|y(t)|$ at parton level per bin ([0,0.3,0.7,1.3,2.2]) in percent of $\dfrac{\mathrm{d}\sigma(\bar tq)}{\mathrm{d}|y(t)|}$.

Uncertainties for the normalised differential $ \bar tq $ cross-section as a function of $|y(t)|$ at parton level per bin ([0,0.3,0.7,1.3,2.2]) in percent of $\left(\dfrac{1}{\sigma}\right)\dfrac{\mathrm{d}\sigma(\bar tq)}{\mathrm{d}|y(t)|}$. If the uncertainty reported in the paper is "0.0" for both the $\textit{plus}$ and $\textit{minus}$ variation, the value "+0.01" is assigned to the $\textit{plus}$ variation for technical reasons.