Search for standard model production of four top quarks with same-sign and multilepton final states in proton–proton collisions at s=13TeV\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s} = 13\,\text {TeV} $$\end{document}

A search for standard model production of four top quarks (tt¯tt¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm {t}\overline{\mathrm {t}} \mathrm {t}\overline{\mathrm {t}} $$\end{document}) is reported using events containing at least three leptons (e,μ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm {e}, \mathrm {\mu }$$\end{document}) or a same-sign lepton pair. The events are produced in proton–proton collisions at a center-of-mass energy of 13TeV\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\,\text {TeV}$$\end{document} at the LHC, and the data sample, recorded in 2016, corresponds to an integrated luminosity of 35.9fb-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\,\text {fb}^{-1}$$\end{document}. Jet multiplicity and flavor are used to enhance signal sensitivity, and dedicated control regions are used to constrain the dominant backgrounds. The observed and expected signal significances are, respectively, 1.6 and 1.0 standard deviations, and the tt¯tt¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm {t}\overline{\mathrm {t}} \mathrm {t}\overline{\mathrm {t}} $$\end{document} cross section is measured to be 16.9-11.4+13.8\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$16.9^{+13.8}_{-11.4}$$\end{document}fb\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\,\text {fb}$$\end{document}, in agreement with next-to-leading-order standard model predictions. These results are also used to constrain the Yukawa coupling between the top quark and the Higgs boson to be less than 2.1 times its expected standard model value at 95% confidence level.

After the decays of the top quarks, the final state contains several jets resulting from the hadronization of light quarks and b quarks (b jets), and may contain isolated leptons and missing transverse momentum depending on the decays of the W bosons [19]. Among these final states, the same-sign dilepton and the three-(or more) lepton final states, considering = e, μ, correspond to branching fractions in tttt events of 8 and 1%, respectively, excluding the small contribution from W → τ ν, which is included in selected events. However, due to the low level of backgrounds, these channels are the most sensitive to tttt production in the regime with SM-like kinematic properties. The ATLAS and CMS Collaborations at the CERN LHC have previously searched for SM tttt production in √ s = 8 and 13 TeV pp collisions [20][21][22][23][24]. The most sensitive of these results is a re-interpretation of the CMS same-sign dilepton search for BSM physics at 13 TeV [23], with an observed (expected) tttt cross section upper limit (assuming no SM tttt signal) of 42 (27 +13 −8 ) fb at the 95% confidence level (CL).
The previous search is inclusive, exploring the final state with two same-sign leptons and at least two jets, using an integrated luminosity of 35.9 fb −1 [23]. The analysis described in this paper is based on the same data set and improves on the previous search by optimizing the signal selection for sensitivity to SM tttt production, by using an improved b jet identification algorithm, and by employing background estimation techniques that are adapted to take into account the higher jet and b jet multiplicity requirements in the signal regions.

Background and signal simulation
Monte Carlo (MC) simulations at NLO are used to evaluate the tttt signal acceptance and to estimate the back- ground from diboson (WZ, ZZ, Zγ , W ± W ± ) and triboson (WWW, WWZ, WZZ, ZZZ, WWγ , WZγ ) processes, as well as from production of single top quarks (tZq, tγ ), or tt produced in association with a boson (ttW, ttZ/γ * , ttH). These samples are generated using the NLO Mad- Graph5_amc@nlo 2.2.2 [17] program with up to one additional parton in the matrix-element calculation, except for W ± W ± which is generated with up to two additional partons, and the WZ, ZZ and ttH samples, which are generated with no additional partons with the powheg box v2 [25,26] program. The ttZ/γ * sample with Z/γ * → is generated with a dilepton invariant mass greater than 1 GeV. The LO MadGraph5_amc@nlo generator, scaled to NLO cross sections, is used to estimate the Wγ and ttγ processes with up to three additional partons. Other rare backgrounds, such as tt production in association with dibosons (ttWW, ttWZ, ttZZ, ttWH, ttZW, ttHH), triple top quark production (ttt, tttW), and tWZ are generated using LO MadGraph5_amc@nlo without additional partons, and scaled to NLO cross sections [27]. The NNPDF3.0LO (NNPDF3.0NLO) [28] parton distribution functions (PDFs) are used to generate all LO (NLO) samples. Parton showering and hadronization, as well as W ± W ± from double-parton-scattering, are modeled by the pythia 8.205 [29] program, while the MLM [30] and FxFx [31] prescriptions are employed in matching additional partons in the matrix-element calculations to parton showers in the LO and NLO samples, respectively. The top quark mass in the generators is set to 172.5 GeV. The Geant4 package [32] is used to model the response of the CMS detector. Additional proton-proton interactions (pileup) within the same or nearby bunch crossings are also included in the simulated events.
To improve the MC modeling of the multiplicity of additional jets from initial-state radiation (ISR), simulated ttW and ttZ/γ * events are reweighted based on the number of ISR jets (N ISR jets ). The reweighting is based on a comparison of the light-flavor jet multiplicity in dilepton tt events in data and simulation. The method requires exactly two jets identified as originating from b quarks in dilepton tt events, and assumes that all other jets are from ISR. Weighting factors are obtained as a function of N ISR jets to bring data and MC into agreement. These weights are then applied, keeping the total cross section constant, to ttW and ttZ/γ * MC as a function of the number of jets not originating from top quark, W, or Z decays. To improve the modeling of the flavor of additional jets, the simulation is also corrected to account for the measured ratio of ttbb/ttjj cross sections reported in Ref. [33]. More details on these corrections and their uncertainties are provided in Sect. 6.

The CMS detector and event reconstruction
The central feature of the CMS detector is a superconducting solenoid of 6 m internal diameter, providing a magnetic field of 3.8 T. Within the solenoid volume are a silicon pixel and strip tracker, a lead tungstate crystal electromagnetic calorimeter (ECAL), and a brass and scintillator hadron calorimeter (HCAL), each composed of a barrel and two endcap sections. Forward calorimeters extend the pseudorapidity (η) coverage provided by the barrel and endcap detectors. Muons are measured in gas-ionization detectors embedded in the steel flux-return yoke outside the solenoid. A more detailed description of the CMS detector, together with a definition of the coordinate system used and the relevant kinematic variables, can be found in Ref. [34].
Events of interest are selected using a two-tiered trigger system [35]. The first level (L1), composed of custom hardware processors, uses information from the calorimeters and muon detectors to select events at a rate of around 100 kHz within a time interval of less than 4 μs. The second level, known as the high-level trigger (HLT), consists of a farm of processors running a version of the full event reconstruction software optimized for fast processing, and reduces the event rate to less than 1 kHz before data storage.
Events are processed using the particle-flow (PF) algorithm [36], which reconstructs and identifies each individual particle with an optimized combination of information from the various elements of the CMS detector. The energy of photons is directly obtained from the ECAL measurement. The energy of electrons is determined from a combination of the electron momentum at the primary interaction vertex as determined by the tracker, the energy of the corresponding ECAL cluster, and the energy sum of all bremsstrahlung photons spatially compatible with the electron track [37]. The momentum of muons is obtained from the curvature of the corresponding track, combining information from the silicon tracker and the muon system [38]. The energy of charged hadrons is determined from a combination of their momentum measured in the tracker and the matching ECAL and HCAL energy deposits, corrected for the response function of the calorimeters to hadronic showers. The energy of neutral hadrons is obtained from the corresponding corrected ECAL and HCAL energies.
Hadronic jets are clustered from neutral PF candidates and charged PF candidates associated with the primary vertex, using the anti-k T algorithm [39,40] with a distance parameter of 0.4. The jet momentum is determined as the vectorial sum of all PF candidate momenta in the jet. An offset correction is applied to jet energies to take into account the contribution from pileup. Jet energy corrections are derived from simulation, and are improved with in situ measurements of the energy balance in dijet, multijet, γ +jet and leptonically decaying Z+jet events [41,42]. Additional selection criteria are applied to each event to remove spurious jet-like features originating from isolated noise patterns in certain HCAL regions. Jets originating from b quarks are identified as b-tagged jets using a deep neural network algorithm [43], with a working point chosen such that the efficiency to identify a b jet is 55-70% for a jet transverse momentum ( p T ) between 20 and 400 GeV. The misidentification rate for a light-flavor jet is 1-2% in the same jet p T range. The vector p miss T is defined as the projection on the plane perpendicular to the beams of the negative vector sum of the momenta of all reconstructed PF candidates in an event [44]. Its magnitude, called missing transverse momentum, is referred to as p miss T . The scalar p T sum of all jets in an event is referred to as H T .

Event selection and search strategy
The definitions of objects and the baseline event selection follow closely those of Refs. [23,45]. Electron identification is based on a multivariate discriminant using shower shape and track quality variables, while for muons it is based on the quality of the geometrical matching between the tracker and muon system measurements. Isolation and impact parameter requirements are applied to both lepton flavors, as well as specific selections designed to improve the accuracy of the charge reconstruction. The combined reconstruction and identification efficiency is in the range of 45-70% (70-90%) for electrons (muons), increasing as a function of p T and converging to the maximum value for p T > 60 GeV. The number of leptons (N ), the number of jets (N jets ), and the number of b-tagged jets (N b ) are counted after the application of the basic kinematic requirements summarized in Table 1.
Signal events are selected using triggers that require two leptons with p T > 8 GeV and H T > 300 GeV. The trigger efficiency is greater than 95% for di-electron (ee) and electron-muon (eμ) events and about 92% for di-muon (μμ) events. The baseline selections require H T > 300 GeV and p miss T > 50 GeV, at least two jets (N jets ≥ 2), at least two b-tagged jets (N b ≥ 2), a leading lepton with p T > 25 GeV, and a second lepton of the same charge with p T > 20 GeV. To reduce the background from Drell-Yan with a chargemisidentified electron, events with same-sign electron pairs with mass below 12 GeV are rejected. Events where a third lepton with p T larger than 5 (7) GeV for muons (electrons) forms an opposite-sign (OS) same-flavor pair with mass below 12 GeV or between 76 and 106 GeV are also rejected. If the third lepton has p T > 20 GeV and the invariant mass of the pair is between 76 and 106 GeV, these rejected events are used to populate a ttZ background control region (CRZ). The signal acceptance in the baseline region, including the leptonic W boson branching fraction, is approximately 1.5%. After these requirements, we define eight mutually exclusive signal regions (SRs) and a control region for the ttW background (CRW), based on N jets , N b , and N , as detailed in Table 2. The observed and predicted yields in the control and  Table 2 Definitions of the eight SRs and the two control regions for ttW (CRW) and ttZ (CRZ) Inverted Z veto CRZ signal regions are used to measure σ (pp → tttt), following the procedure described in Sect. 7.

Backgrounds
The main backgrounds to the tttt process in the same-sign dilepton and three-(or more) lepton final states arise from rare multilepton processes, such as ttW, ttZ/γ * , and ttH (H → WW), and single-lepton or OS dilepton processes with an additional "nonprompt lepton". Nonprompt leptons consist of electrons from conversions of photons in jets and leptons from the decays of heavy-or light-flavor hadrons. In this category we include also hadrons misidentified as leptons. The minor background from OS dilepton events with a charge-misidentified lepton is also taken into account. Rare multilepton processes are estimated using simulated events. Control regions are used to constrain the normalization of the ttW and ttZ backgrounds, as described in Sect. 7, while for other processes the normalization is based on the NLO cross sections referenced in Sect. 2. Processes such as the associated production of a tt pair with a pair of bosons (W, Z, H) are grouped into a "ttVV" category. Associated photon production processes such as Wγ , Zγ , ttγ , and tγ , where an electron is produced in an unidentified photon conversion, are grouped into a "Xγ " category. All residual processes with very small contributions, including diboson (WZ, ZZ, W ± W ± from single-and double-parton scattering), triboson (WWW, WWZ, WZZ, ZZZ, WWγ , WZγ ), and rare single top quark (tZq, tWZ) and triple top quark processes (ttt and tttW), are grouped into a "Rare" category.
The nonprompt lepton and charge-misidentified lepton backgrounds are estimated following the methods described in Ref. [23]. For nonprompt leptons, an estimate referred to as the "tight-to-loose" method defines two control regions by modifying the lepton identification (including isolation) and event kinematic requirements, respectively. An "application region" is defined for every SR by requiring at least one lepton to fail the standard identification ("tight") while satisfying a more relaxed one ("loose"). To obtain the nonprompt lepton background estimate in the corresponding SR, the event yield in each application region is weighted by a factor of TL /(1 − TL ) for each lepton failing the tight requirement. The TL parameter is the probability that a nonprompt lepton that satisfies a loose lepton selection also satisfies the tight selection. It is extracted as a function of lepton flavor and kinematic properties from a "measurement region" that consists of a single-lepton events with event kinematic properties designed to suppress the W → ν contribution.
For charge-misidentified leptons, an OS dilepton control region is defined for each same-sign dilepton signal region. Its yield is then weighted by the charge misidentification probability estimated in simulation, which ranges between 10 −5 and 10 −3 for electrons and is negligible for muons.

Systematic uncertainties
The sources of experimental and theoretical uncertainty for the data and simulations are summarized in Table 3. The uncertainty in the integrated luminosity is 2.5% [46]. The simulation is reweighted to match the distribution in the number of pileup collisions per event in data. The uncertainty in the inelastic cross section propagated to the final yields provides an uncertainty of at most 6%. Trigger efficiencies are measured with an uncertainty of 2% in an independent data sample selected using singlelepton triggers. Lepton-efficiency scale factors, used to account for differences in the reconstruction and identification efficiencies between data and simulation, are measured using a "tag-and-probe" method in data enriched in Z → events [37,38]. The scale factors are applied to all simulated processes with an uncertainty per lepton of approximately 3% for muons and 4% for electrons.
The uncertainty in the calibration of the jet energy scale depends on the p T and η of the jet and results in a 1-15% variation in the event yield in a given SR. The uncertainty due to the jet energy resolution is estimated by broadening the resolution in simulation [42], and the resulting effect is a change of 1-5% in the SR yields. The b tagging efficiency in simulation is corrected using scale factors determined from efficiencies measured in data and simulation [47]. The uncertainty in the measured scale factors results in an overall effect between 1 and 15%, again depending on the SR.
As mentioned in Sect. 2, ttW and ttZ/γ * simulated events are reweighted to match the number of additional jets observed in data. The reweighting factors vary between 0.92 for N ISR jets = 1 and 0.77 for N ISR jets ≥ 4. Half of the difference from unity is taken as a systematic uncertainty in these reweighting factors to cover differences observed between data and simulation when the factors are used to reweight simulation in a control sample enriched in single-lepton tt events. Uncertainties in the reweighting factors are treated as correlated among regions. Simulated ttW and ttZ/γ * events with two b quarks not originating from top quark decay are also weighted to account for the CMS measurement of the ratio of cross sections σ (ttbb)/σ (ttjj), which was found to be a factor of 1.7 ± 0.6 larger than the MC prediction [33]. In signal regions requiring four b-tagged jets, where the effect is dominant, this results in a systematic uncertainty of up to 15% on the total background prediction. In signal regions requiring three b-tagged jets, the dominant origin of the additional b-tagged jet is a charm quark from a W decay, so the effect is negligible.
Uncertainties in the renormalization and factorization scales (varied by a factor of two) and from the choice of PDF [48,49] affect the number of events expected (normalization) in the simulated background processes, as well as the acceptance for the tttt signal. The effects of these uncertainties on the relative distribution of events in the signal regions (shape) are also considered. For the ttW and ttZ/γ * backgrounds, the normalization uncertainty is 40%, while for ttH a 50% normalization uncertainty reflects the signal strength of 1.5 ± 0.5 measured by CMS [50]. The processes in the Rare category along with Xγ and ttVV, many of which have never been observed, are expected to give small contributions to the event yields in the signal regions. We assign separate 50% normalization uncertainties to each of these three categories. The shape uncertainty resulting from variations of the renormalization and factorization scales is as large as 15% for the ttW, ttZ/γ * , and ttH backgrounds, and 10% for the tttt signal, while the effect from the PDF is only 1%. For the signal, the uncertainty in the acceptance from variations of the scales (PDFs) is 2% (1%). In addition, for the tttt signal, the scales that determine ISR and final-state radiation (FSR) in the parton shower are also varied, resulting in a 6% change in the acceptance and shape variations as large as 15%.
For nonprompt and charge-misidentified lepton backgrounds, the statistical uncertainty from the application region depends on the SR considered. The background from misidentified charge is assigned a systematic uncertainty of 20%, based on comparisons of the expected number of samesign events estimated from an OS control sample and the observed same-sign yield in a control sample enriched in Z → e + e − events with one electron or positron having a misidentified charge.
In addition to the statistical uncertainty, the nonprompt lepton background is assigned an overall normalization uncertainty of 30% to cover variations observed in closure tests performed with simulated multijet and tt events. This uncertainty is increased to 60% for electrons with p T > 50 GeV, to account for trends observed at high p T in the closure tests. We also include an uncertainty related to the subtraction of events with prompt leptons (from electroweak processes with a W or Z boson) in the measurement region, which has an effect between 1% and 50%, depending on the SR. The prompt lepton contamination was also checked in the application region, where it was found to be below 1%.
Experimental uncertainties are treated as correlated among signal regions for all signal and background processes. Systematic uncertainties in data-driven estimates and theoretical uncertainties are treated as uncorrelated between processes, but correlated among signal regions. Statistical uncertainties from the limited number of simulated events or in the number of events in data control regions are considered uncorrelated.

Results and interpretation
The properties of events in the signal regions (SR 1-8 as defined in Table 2) are shown in Fig. 2, where distributions of the main kinematic variables in the data (N jets , N b , H T , and p miss T ) are compared to SM background predictions. The N jets and N b distributions for CRW and CRZ are shown in Fig. 3. In both figures we overlay the expected SM tttt signal, scaled by a factor of 5. The SM predictions are generally consistent with the observations, with some possible underestimation in CRW and CRZ.
The yields from SR 1-8, CRW, and CRZ are combined in a maximum-likelihood fit, following the procedures described in Ref.
[51], to estimate a best-fit cross section for tttt, the significance of the observation relative to the background-only hypothesis, and the upper limit on σ (pp → tttt). The experimental and theoretical uncertainties described in Sect. 6 are incorporated in the likelihood as "nuisance" parameters and are profiled in the fit. Nuisance parameters corresponding to systematic uncertainties are parameterized as log-normal distributions. The fitted values of the nuisance parameters are found to be consistent with their initial values within uncertainties. The nuisance parameters corresponding to the ttW and ttZ/γ * normalizations are scaled by 1.2 ± 0.3 and 1.3±0.3, respectively, while other background contributions including ttH are scaled up by 1.1 or less. The signal and control region results after the maximum-likelihood fit (postfit) are shown in Fig. 4, with the fitted tttt signal contribution added to the background predictions, which are given in 18 x 5 t t t t x 5 t t t t Data x 5 t t t t 10 x 5 t t t t  120 x 5 t t t t x 5 t t t t x 5 t t t t Data    [54,55], is found to be 41.7 fb. The corresponding expected upper limit, assuming no SM tttt contribution to the data, is 20.8 +11.2 −6.9 fb, showing a significant improvement relative to the value of 27 fb of Ref. [23].
The pp → tttt process has contributions from diagrams with virtual Higgs bosons, as shown in Fig. 1. Experimental information on σ (pp → tttt) can therefore be used to constrain the Yukawa coupling, y t , between the top quark and the Higgs boson. We constrain y t assuming that the signal acceptance is not affected by the relative contribution of the virtual Higgs boson diagrams. As the cross section for the ttH background also depends on the top quark Yukawa coupling, for the purpose of constraining y t the fit described above is repeated with the ttH contribution scaled by the square of the absolute value of the ratio of the top quark Yukawa coupling to its SM value (|y t /y SM t | 2 ), where This results in a dependence of the measured σ (pp → tttt) on |y t /y SM t | which is shown in Fig. 5 and is compared to its theoretical prediction. The prediction is obtained from the LO calculation of Ref. [16], with an NLO/LO K -factor of 1.27 [18]. The LO calculation is used instead of the NLO one, as Ref.
[16] provides a breakdown of the contributions to the cross section according to powers of y t . The prediction also includes the uncertainty associated with varying the renormalization and factorization scales in the LO calculation by a factor of 2. The central, upper and lower values of the theoretical cross section provide respective 95% CL limits for |y t /y SM t | < 2.1, < 1.9 and < 2.4.

Summary
The results of a search for standard model (SM) production of tttt at the LHC have been presented, using data from √ s = 13 TeV proton-proton collisions corresponding to an integrated luminosity of 35.9 fb −1 , collected with the CMS detector in 2016. The analysis strategy uses same-sign dilepton as well as three-(or more) lepton events, relying on jet multiplicity and jet flavor to define search regions that are used to probe the tttt process. Combining these regions yields a significance of 1.6 standard deviations relative to the background-only hypothesis, and a measured value for the tttt cross section of 16.9 +13.8 −11.4 fb, in agreement with the standard model predictions. The results are also re-interpreted to constrain the ratio of the top quark Yukawa coupling to its SM value, yielding |y t /y SM t | < 2.1 at 95% confidence level. straint on the top quark Yukawa coupling. We congratulate our colleagues in the CERN accelerator departments for the excellent performance of the LHC and thank the technical and administrative staffs at CERN and at other CMS institutes for their contributions to the success of the CMS effort. In addition, we gratefully acknowledge the computing centers and personnel of the Worldwide LHC Com-puting Grid for delivering so effectively the computing infrastructure essential to our analyses. Finally, we acknowledge the enduring support for the construction and operation of the LHC and the CMS detector provided by the following funding agencies: BMWFW and FWF (