This paper presents the observation of four-top-quark ($t\bar{t}t\bar{t}$) production in proton-proton collisions at the LHC. The analysis is performed using an integrated luminosity of 140 fb$^{-1}$ at a centre-of-mass energy of 13 TeV collected using the ATLAS detector. Events containing two leptons with the same electric charge or at least three leptons (electrons or muons) are selected. Event kinematics are used to separate signal from background through a multivariate discriminant, and dedicated control regions are used to constrain the dominant backgrounds. The observed (expected) significance of the measured $t\bar{t}t\bar{t}$ signal with respect to the standard model (SM) background-only hypothesis is 6.1 (4.3) standard deviations. The $t\bar{t}t\bar{t}$ production cross section is measured to be $22.5^{+6.6}_{-5.5}$ fb, consistent with the SM prediction of $12.0 \pm 2.4$ fb within 1.8 standard deviations. Data are also used to set limits on the three-top-quark production cross section, being an irreducible background not measured previously, and to constrain the top-Higgs Yukawa coupling and effective field theory operator coefficients that affect $t\bar{t}t\bar{t}$ production.
Post-fit distributions for the number of jets ($N_{j}$) in CR $t\bar{t}W^{+}$+jets. The QmisID represents the backgrounds with a mis-assigned charge. HF e and HF $\mu$ are the backgrounds with fake/non-prompt leptons. Mat. Conv. and Low $m_{\gamma*}$ are the material and virtual photon conversions.
Post-fit distributions for the number of jets ($N_{j}$) in CR $t\bar{t}W^{-}$+jets. The QmisID represents the backgrounds with a mis-assigned charge. HF e and HF $\mu$ are the backgrounds with fake/non-prompt leptons. Mat. Conv. and Low $m_{\gamma*}$ are the material and virtual photon conversions.
Post-fit distributions for the number of jets ($N_{j}$) in CR 1b(+). The QmisID represents the backgrounds with a mis-assigned charge. HF e and HF $\mu$ are the backgrounds with fake/non-prompt leptons. Mat. Conv. and Low $m_{\gamma*}$ are the material and virtual photon conversions.
Differential and double-differential distributions of kinematic variables of leptons from decays of top-quark pairs ($t\bar{t}$) are measured using the full LHC Run 2 data sample collected with the ATLAS detector. The data were collected at a $pp$ collision energy of $\sqrt{s}=13$ TeV and correspond to an integrated luminosity of 140 fb$^{-1}$. The measurements use events containing an oppositely charged $e\mu$ pair and $b$-tagged jets. The results are compared with predictions from several Monte Carlo generators. While no prediction is found to be consistent with all distributions, a better agreement with measurements of the lepton $p_{\text{T}}$ distributions is obtained by reweighting the $t\bar{t}$ sample so as to reproduce the top-quark $p_{\text{T}}$ distribution from an NNLO calculation. The inclusive top-quark pair production cross-section is measured as well, both in a fiducial region and in the full phase-space. The total inclusive cross-section is found to be \[ \sigma_{t\bar{t}} = 829 \pm 1\;(\textrm{stat}) \pm 13\;(\textrm{syst}) \pm 8\;(\textrm{lumi}) \pm 2\; (\textrm{beam})\ \textrm{pb}, \] where the uncertainties are due to statistics, systematic effects, the integrated luminosity and the beam energy. This is in excellent agreement with the theoretical expectation.
Definition of the fiducial phase space with the lepton candidate, electron $e$ and muon $\mu$, and jets.
Breakdown of systematic uncertainties in the measured fiducial cross-section. The impact of the top-quark mass on the cross-section is included in the table and not counted in the total uncertainty entry in the paper.
Data bootstrap post unfolding for the fiducial cross-section. The replicas are obtained by reweighting each observed data event by a random integer generated according to Poisson statistics, using the BootstrapGenerator software package (https://gitlab.cern.ch/atlas-physics/sm/StandardModelTools_BootstrapGenerator/BootstrapGenerator), which implements a technique described in ATL-PHYS-PUB-2021-011 (https://cds.cern.ch/record/2759945). 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. All the provided numbers originate from pseudo-data, including the 0th entry, and are in units of [fb].
The first measurement of the top quark pair ($\mathrm{t\bar{t}}$) production cross section in proton-proton collisions at $\sqrt{s}$ = 13.6 TeV is presented. Data recorded with the CMS detector at the CERN LHC in Summer 2022, corresponding to an integrated luminosity of 1.21 fb$^{-1}$, are analyzed. Events are selected with one or two charged leptons (electrons or muons) and additional jets. A maximum likelihood fit is performed in event categories defined by the number and flavors of the leptons, the number of jets, and the number of jets identified as originating from b quarks. An inclusive $\mathrm{t\bar{t}}$ production cross section of 881 $\pm$ 23 (stat+syst) $\pm$ 20 (lumi) pb is measured, in agreement with the standard model prediction of 924 $^{+32}_{-40}$ pb.
Comparison of the number of observed (points) and predicted (filled histograms) events in the final analysis binning. The predictions are shown before fitting the model to the data. The lower panel of the plot displays the ratio of the event yields in data to the sum of predicted signal and background yields. The vertical bars on the points represent the statistical uncertainties in the data, while the hatched bands represent systematic uncertainty in the predictions, excluding the integrated luminosity. No b jet efficiency scale factors are applied in the plot, and no systematic uncertainty entering into the hatched bands is intended to cover thesefactors, which are free parameters in the fit.
Comparison of the number of observed (points) and predicted (filled histograms) events in the final analysis binning. The predictions are shown after fitting the model to the data. The lower panel of each plot displays the ratio of the event yields in data to the sum of predicted signal and background yields. The vertical bars on the points represent the statistical uncertainties in the data, while the hatched bands represent systematic uncertainty in the predictions, excluding the integrated luminosity. The hatched bands are greatly reduced due to additional constraint of the nuisances parameters as well as correlations between them.
Result for the inclusive ttbar production cross section
A search is reported for excited $\tau$-leptons and leptoquarks in events with two hadronically decaying $\tau$-leptons and two or more jets. The search uses proton-proton (pp) collision data at $\sqrt{s} = 13$ TeV recorded by the ATLAS experiment during the Run 2 of the Large Hadron Collider in 2015-2018. The total integrated luminosity is 139 fb$^{-1}$. The excited $\tau$-lepton is assumed to be produced and to decay via a four-fermion contact interaction into an ordinary $\tau$-lepton and a quark-antiquark pair. The leptoquarks are assumed to be produced in pairs via the strong interaction, and each leptoquark is assumed to couple to a charm or lighter quark and a $\tau$-lepton. No excess over the background prediction is observed. Excited $\tau$-leptons with masses below 2.8 TeV are excluded at 95% CL in scenarios with the contact interaction scale $\Lambda$ set to 10 TeV. At the extreme limit of model validity where $\Lambda$ is set equal to the excited $\tau$-lepton mass, excited $\tau$-leptons with masses below 4.6 TeV are excluded. Leptoquarks with masses below 1.3 TeV are excluded at 95% CL if their branching ratio to a charm quark and a $\tau$-lepton equals 1. The analysis does not exploit flavour-tagging in the signal region.
Observed and expected upper 95% CL limit on the $\tau^\ast$ production cross-section as a function of $m_{\tau^\ast}$ for a fixed value of the contact interaction scale, $\Lambda = 10$ TeV.
Observed and expected lower 95% CL limit on the contact interaction scale $\Lambda$ as a function of $m_{\tau^\ast}$.
Observed and expected upper 95% CL limit on the LQ production cross-section as a function of $m_\mathrm{LQ}$. The LQ couples to a tau lepton and a c-quark. The limits are also valid for scenarios in which the LQ couples to lighter quarks.
A generic search is presented for the associated production of a Z boson or a photon with an additional unspecified massive particle X, pp $\to$ pp + Z/$\gamma$ + X, in proton-tagged events from proton-proton collisions at $\sqrt{s}$ = 13 TeV, recorded in 2017 with the CMS detector and the CMS-TOTEM precision proton spectrometer. The missing mass spectrum is analysed in the 600-1600 GeV range and a fit is performed to search for possible deviations from the background expectation. No significant excess in data with respect to the background predictions has been observed. Model-independent upper limits on the visible production cross section of pp $\to$ pp + Z/$\gamma$ + X are set.
Comparison of the $m_{miss}$ shapes for the simulated signal events within the fiducial region and those outside it, after including the effect of PU protons as describe in the text, for a generated $m_{X}$ mass of 1000 GeV. The distributions are shown for multi(+z)-multi(−z) proton reconstruction categories.
Comparison of the $m_{miss}$ shapes for the simulated signal events within the fiducial region and those outside it, after including the effect of PU protons as describe in the text, for a generated $m_{X}$ mass of 1000 GeV. The distributions are shown for multi(+z)-single(−z) proton reconstruction categories.
Comparison of the $m_{miss}$ shapes for the simulated signal events within the fiducial region and those outside it, after including the effect of PU protons as describe in the text, for a generated $m_{X}$ mass of 1000 GeV. The distributions are shown for single(+z)-multi(−z) proton reconstruction categories.
The production of four top quarks ($\mathrm{t\bar{t}t\bar{t}}$) is studied with LHC proton-proton collision data samples collected by the CMS experiment at a center-of-mass energy of 13 TeV, and corresponding to integrated luminosities of up to 138 fb$^{-1}$. Events that have no leptons (all-hadronic), one lepton, or two opposite-sign leptons (where lepton refers only to prompt electrons or prompt muons) are considered. This is the first $\mathrm{t\bar{t}t\bar{t}}$ measurement that includes the all-hadronic final state. The observed significance of the $\mathrm{t\bar{t}t\bar{t}}$ signal in these final states of 3.9 standard deviations (1.5 expected) provides evidence for $\mathrm{t\bar{t}t\bar{t}}$ production, with a measured cross section of 36 $^{+12}_{-11}$ fb. Combined with earlier CMS results in other final states, the signal significance is 4.0 standard deviations (3.2 expected). The combination returns an observed cross section of 17 $\pm$ 4 (stat) $\pm$ 3 (syst) fb, which is consistent with the standard model prediction.
The jet multiplicity for $N_\textrm{b} \geq 4$ in the opposite-sign dilepton channel for the combined 2017--2018 dataset with dilepton decay categories combined. Here, $\textrm{t}\bar{\textrm{t}} + \geq 1 \textrm{b}$ refers to $\textrm{t}\bar{\textrm{t}}$ events with at least one additional b jet, $\textrm{t}\bar{\textrm{t}} + 0 \textrm{b}$ includes all other $\textrm{t}\bar{\textrm{t}}$ events not produced in association with a boson, and EW refers to events that contain W and Z bosons but no top quarks. The backgrounds and $\textrm{t}\bar{\textrm{t}}\textrm{t}\bar{\textrm{t}}$ signal (derived from the fit) are shown as a stacked histogram. The hatched bands correspond to the estimated total uncertainty after the fit.
The distribution of the BDT discriminants for the 2016--2018 data set for three different categories in the combined single-electron and single-muon channels. The three categories are defined by the number of resolved t tags ($N_\textrm{RT}$), b tags ($N_\textrm{b}$), and jets ($N_\textrm{j}$), selected as representative based on their sensitivity to signal. Here, $\textrm{t}\bar{\textrm{t}} + \geq 1 \textrm{b}$ refers to $\textrm{t}\bar{\textrm{t}}$ events with at least one additional b jet, while $\textrm{t}\bar{\textrm{t}} + 0 \textrm{b}$ includes all other $\textrm{t}\bar{\textrm{t}}$ events not produced in association with a boson. The TOP grouping contains single top quark production along with the other $\textrm{t}\bar{\textrm{t}}$ processes not explicitly shown, and EW refers to events that contain W and Z bosons but no top quarks. The backgrounds and $\textrm{t}\bar{\textrm{t}}\textrm{t}\bar{\textrm{t}}$ signal (derived from the fit) are shown as a stacked histogram. The hatched bands correspond to the estimated total uncertainty after the fit. While the bins are shown to be equal width, they do not correspond to equal width in BDT value.
The distribution of the BDT discriminants for the full 2016--2018 data set in the all-hadronic channel. The sample VR category shown is defined by $N_\textrm{RT}=1$, $N_\textrm{BT} \geq 1$, $H_T > 1400$ GeV. The background from QCD multijet and $\textrm{t}\bar{\textrm{t}}$ production is derived from control regions in the data. Estimates for the $\textrm{t}\bar{\textrm{t}}\textrm{t}\bar{\textrm{t}}$ signal and other backgrounds are shown using simulated samples. The hatched bands correspond to the estimated total uncertainty.
A search for decays to invisible particles of Higgs bosons produced in association with a top-antitop quark pair or a vector boson, which both decay to a fully hadronic final state, has been performed using proton-proton collision data collected at $\sqrt{s}$ = 13 TeV by the CMS experiment at the LHC, corresponding to an integrated luminosity of 138 fb$^{-1}$. The 95% confidence level upper limit set on the branching fraction of the 125 GeV Higgs boson to invisible particles, $\mathcal{B}$(H $\to$ inv), is 0.54 (0.39 expected), assuming standard model production cross sections. The results of this analysis are combined with previous $\mathcal{B}$(H $\to$ inv) searches carried out at $\sqrt{s}$ = 7, 8, and 13 TeV in complementary production modes. The combined upper limit at 95% confidence level on $\mathcal{B}$(H $\to$ inv) is 0.15 (0.08 expected).
Observed and expected 95% CL upper limits on ${{(\sigma_{\text{H}}/\sigma_{\text{H}}^{\mathrm{SM}}) \times {{\mathcal{B}(\text{H} \to \text{inv})}}}}$ for hadronic final states of ttH and resolved VH channels, and their combination, using data from 2016--2018 and assuming a SM Higgs boson with a mass of 125 GeV.
Observed and expected 95% CL upper limits on ${{(\sigma_{\text{H}}/\sigma_{\text{H}}^{\mathrm{SM}}) \times {{\mathcal{B}(\text{H} \to \text{inv})}}}}$ for the VBF, ttH, VH and ggH channels using all available CMS data, and their combination, assuming a SM Higgs boson with a mass of 125 GeV.
Observed and expected 95% CL upper limits on ${{(\sigma_{\text{H}}/\sigma_{\text{H}}^{\mathrm{SM}}) \times {{\mathcal{B}(\text{H} \to \text{inv})}}}}$ for the VBF, ttH, VH and ggH channels using Run2 CMS data, and their combination, assuming a SM Higgs boson with a mass of 125 GeV.
A search for pair-produced scalar or vector leptoquarks decaying into a $b$-quark and a $\tau$-lepton is presented using the full LHC Run 2 (2015-2018) data sample of 139 fb$^{-1}$ collected with the ATLAS detector in proton-proton collisions at a centre-of-mass energy of $\sqrt{s} =13$ TeV. Events in which at least one $\tau$-lepton decays hadronically are considered, and multivariate discriminants are used to extract the signals. No significant deviations from the Standard Model expectation are observed and 95% confidence-level upper limits on the production cross-section are derived as a function of leptoquark mass and branching ratio $B$ into a $\tau$-lepton and $b$-quark. For scalar leptoquarks, masses below 1460 GeV are excluded assuming $B=100$%, while for vector leptoquarks the corresponding limit is 1650 GeV (1910 GeV) in the minimal-coupling (Yang-Mills) scenario.
Acceptance $\times$ efficiency for the $\tau_\text{lep}\tau_\text{had}$ signal region assuming $\beta$ = 0.5 as a function of m$_\text{LQ}$.
Acceptance $\times$ efficiency for the $\tau_\text{had}\tau_\text{had}$ signal region assuming $\beta$ = 0.5 as a function of m$_\text{LQ}$.
The observed and expected 95% CL upper limits on the scalar LQ pair production cross-sections assuming B = 1 as a function of m$_\text{LQ}$.
A search for the electroweak production of a vector-like quark T$'$, decaying to a top quark and a Higgs boson is presented. The search is based on a sample of proton-proton collision events recorded at the LHC at $\sqrt{s}$ = 13 TeV, corresponding to an integrated luminosity of 138 fb$^{-1}$. This is the first T$'$ search that exploits the Higgs boson decay to a pair of photons. For narrow isospin singlet T$'$ states with masses up to 1.1 TeV, the excellent diphoton invariant mass resolution of 1-2% results in an increased sensitivity compared to previous searches based on the same production mechanism. The electroweak production of a T$'$ quark with mass up to 960 GeV is excluded at 95% confidence level, assuming a coupling strength $\kappa_\mathrm{T}$ = 0.25 and a relative decay width $\Gamma/M_{\mathrm{T}'}$ $\lt$ 5%.
The combined, leptonic plus hadronic, expected (dotted black) and observed (solid black) upper limits at 95% CL on $\sigma_{T^\prime bq} \times \mathcal{B}_{T^\prime \rightarrow tH}$ displayed as a function of $\rm{M}_{T^\prime}$. The green (yellow) band represents the 68% (95%) of the limit values expected under the background-only hypothesis. The theoretical cross sections for the singlet $T^\prime$ production with representative $\kappa_T$-values fixed at 0.1, 0.15, 0.2 and 0.25 (for $\Gamma/\rm{M}_{T^\prime} < 5\%$) are shown as solid red lines. The theoretical cross sections for the singlet $T^\prime$ production with $\Gamma/\rm{M}_{T^\prime}$ = 1 and 5% are also shown in table.
The combined, leptonic plus hadronic, expected (dotted black) and observed (solid black) upper limits at 95% CL on the $T^\prime$ coupling to third-generation quarks, $\kappa_{T}$, under the narrow width approximation (NWA) displayed as a function of $\rm{M}_{T^\prime}$. The green (yellow) band represents the 68% (95%) of the limit values expected under the background-only hypothesis. The theoretical $\kappa_{T}$ values corresponding to the $\Gamma/\rm{M}_{T^\prime}$-values fixed at 1, 2, 3, 4, and 5% are shown as red dashed lines.
The expected (dotted black) and observed (solid black) upper limits at 95% CL on $\sigma_{T^\prime bq(T^\prime \rightarrow tH)}$ in the leptonic channel displayed as a function of $\rm{M}_{T^\prime}$. The green (yellow) band represents the 68% (95%) of the limit values expected under the background-only hypothesis. The theoretical cross sections for the singlet $T^\prime$ production with representative $\kappa_T$-values fixed at 0.1, 0.15, 0.2 and 0.25 (for $\Gamma/\rm{M}_{T^\prime} < 5\%$) are shown as solid red lines.
A search for a charged Higgs boson, $H^{\pm}$, produced in top-quark decays, $t \rightarrow H^{\pm}b$, is presented. The search targets $H^{\pm}$ decays into a bottom and a charm quark, $H^{\pm} \rightarrow cb$. The analysis focuses on a selection enriched in top-quark pair production, where one top quark decays into a leptonically decaying $W$ boson and a bottom quark, and the other top quark decays into a charged Higgs boson and a bottom quark. This topology leads to a lepton-plus-jets final state, characterised by an isolated electron or muon and at least four jets. The search exploits the high multiplicity of jets containing $b$-hadrons, and deploys a neural network classifier that uses the kinematic differences between the signal and the background. The search uses a dataset of proton-proton collisions collected at a centre-of-mass energy $\sqrt{s}=13$ TeV between 2015 and 2018 with the ATLAS detector at CERN's Large Hadron Collider, amounting to an integrated luminosity of 139 fb$^{-1}$. Observed (expected) 95% confidence-level upper limits between 0.15% (0.09%) and 0.42% (0.25%) are derived for the product of branching fractions $\mathscr{B}(t\rightarrow H^{\pm}b) \times \mathscr{B}(H^{\pm}\rightarrow cb)$ for charged Higgs boson masses between 60 and 160 GeV, assuming the SM production of the top-quark pairs.
The observed 95% CL upper limits on $\mathscr{B}=\mathscr{B}(t\rightarrow H^{\pm}b) \times \mathscr{B}(H^{\pm}\rightarrow cb)$ as a function of $m_{H^{\pm}}$ and the expectation (dashed) under the background-only hypothesis. The inner green and outer yellow shaded bands show the $\pm 1\sigma$ and $\pm 2\sigma$ uncertainties of the expected limits. The exclusion limits are presented for $m_{H^{\pm}}$ between 60 and 160 GeV with 10 GeV $m_{H^{\pm}}$ spacing and linear interpolation between adjacent mass points. Superimposed on the upper limits, the predictions from the 3HDM are shown, corresponding to three benchmark values for the parameters $X$, $Y$, and $Z$
Pre-fit event yields in each of the nine analysis regions. The $H^{\pm}$ signal yields for $m_{H^{\pm}}=130$ GeV and $m_{H^{\pm}}=70$ GeV are normalised to $\mathscr{B}_{\mathrm{ref}}=1\%$. The quoted uncertainties are the sum in quadrature of statistical and systematic uncertainties of the yields, computed taking into account correlations among processes resulting from the data-based $t\bar{t}$ correction procedure.
Post-fit yields in each of the nine analysis regions considered. The total prediction is shown after the fit to data under the signal-plus-background hypothesis assuming $H^{\pm}$ signal with $m_{H^{\pm}}=130$ GeV. The predicted yileds for the $H^{\pm}$ signal with $m_{H^{\pm}}=70$ GeV are also shown for reference. The best fit-values of $\mathscr{B}$ for $H^{\pm}$ signal with $m_{H^{\pm}}=130$ GeV and $m_{H^{\pm}}=70$ GeV are 0.16% and 0.07% respectively. The quoted uncertainties are the sum in quadrature of statistical and systematic uncertainties of the yields, computed taking into account correlations among nuisance parameters and among processes.