1 Introduction

The top quark is the heaviest known fundamental particle, with a mass (\({m}_{t}\)) that is much larger than any of the other quarks, and close to the scale of electroweak symmetry breaking. The study of its production and decay properties forms a core part of the ATLAS physics programme at the CERN Large Hadron Collider (LHC). At the LHC, top quarks are primarily produced in quark–antiquark pairs (\(t\overline{t}\)), and the precise prediction of the corresponding inclusive cross-section (\({\sigma _{t\overline{t}}}\)) is a substantial challenge for quantum chromodynamics (QCD) calculational techniques. Precise measurements of \({\sigma _{t\overline{t}}}\) are sensitive to the gluon parton distribution function (PDF), the top quark mass, and potential enhancements of the cross-section due to physics beyond the Standard Model.

Within the Standard Model (SM), the top quark decays almost exclusively to a \(W\) boson and a \(b\) quark, so the final-state topologies in \(t\overline{t}\) production are governed by the decay modes of the two \(W\) bosons. This paper describes a measurement in the dileptonic \(e\mu \) channel, \(t\overline{t}\rightarrow W^{+}bW^{-}\bar{b}\rightarrow {e}^{\pm }\mu ^{\mp }{\nu }{{\overline{\nu }}}{{b\overline{b}}}\), selecting events with an \(e\mu \) pair with opposite-sign electric charges,Footnote 1 and one or two hadronic jets from the \(b\) quarks. Jets originating from \(b\) quarks were identified (‘tagged’) using a \(b\)-tagging algorithm exploiting the long lifetime, high decay multiplicity, hard fragmentation and high mass of \(B\) hadrons. The rates of events with an \(e\mu \) pair and one or two tagged \(b\)-jets were used to measure simultaneously the \(t\overline{t}\) production cross-section and the combined probability to reconstruct and \(b\)-tag a jet from a top quark decay. Events with electrons or muons produced via leptonic \(\tau \) decays \(t\rightarrow Wb\rightarrow \tau \nu b\rightarrow e/\mu \nu \nu \nu b\), were included as part of the \(t\overline{t}\) signal.

The main background is \(Wt\), the associated production of a \(W\) boson and a single top quark. Other background contributions arise from \(Z\rightarrow \tau \tau \rightarrow e\mu \)+jets (\(+4\nu \)) production, diboson+jets production and events where at least one reconstructed lepton does not arise from a \(W\) or \(Z\) boson decay.

Theoretical predictions for \({\sigma _{t\overline{t}}}\) are described in Sect. 2, followed by the data and Monte Carlo (MC) simulation samples in Sect. 3, the object and event selection in Sect. 4, and the extraction of the \(t\overline{t}\) cross-section in Sect. 5. Systematic uncertainties are discussed in Sect. 6, the results, including fiducial cross-section measurements, the extraction of the top quark mass from the measured cross-section and a limit on the production of supersymmetric top squarks, are given in Sect. 7, and conclusions are drawn in Sect. 8.

2 Theoretical cross-section predictions

Calculations of \({\sigma _{t\overline{t}}}\) for hadron collisions are now available at full next-to-next-to-leading-order (NNLO) accuracy in the strong coupling constant \(\alpha _\mathrm{s}\), including the resummation of next-to-next-to-leading logarithmic (NNLL) soft gluon terms [16]. At a centre-of-mass energy of \({\sqrt{s}\,{=}\,7~\mathrm{TeV}}\) and assuming \({m}_{t}\) \(=172.5\) GeV, these calculations give a prediction of \(177.3\pm 9.0\,^{+4.6}_{-6.0}\) pb, where the first uncertainty is due to PDF and \(\alpha _\mathrm{s}\) uncertainties, and the second to QCD scale uncertainties. The corresponding prediction at \({\sqrt{s}=8~\mathrm{TeV}}\) is \(252.9\pm 11.7\,^{+6.4}_{-8.6}\) pb. These values were calculated using the top++ 2.0 program [7]. The PDF and \(\alpha _\mathrm{s}\) uncertainties were calculated using the PDF4LHC prescription [8] with the MSTW2008 68 % CL NNLO [9, 10], CT10 NNLO [11, 12] and NNPDF2.3 5f FFN [13] PDF sets, and added in quadrature to the QCD scale uncertainty. The latter was obtained from the envelope of predictions with the renormalisation and factorisation scales varied independently by factors of two up and down from their default values of \({m}_{t}\), whilst never letting them differ by more than a factor of two. The ratio of cross-sections at \(\sqrt{s}=8~\) TeV and \(\sqrt{s}=7~\) TeV is predicted to be \(1.430\pm 0.013\) (PDF+\(\alpha _\mathrm{s}\)) \(\pm 0.001\) (QCD scale). The total relative uncertainty is only 0.9 %, as the cross-section uncertainties at the two centre-of-mass energies are highly correlated.

The NNLO+NNLL cross-section values are about 3 % larger than the exact NNLO predictions, as implemented in Hathor 1.5 [14]. For comparison, the corresponding next-to-leading-order (NLO) predictions, also calculated using top++ 2.0 with the same set of PDFs, are \(157\pm 12\pm 24\) pb at \(\sqrt{s}=7~\) TeV and \(225\pm 16\pm 29\) pb at \(\sqrt{s}=8~\) TeV, where again the first quoted uncertainties are due to PDF and \(\alpha _\mathrm{s}\) uncertainties, and the second to QCD scale uncertainties. The total uncertainties of the NLO predictions are approximately 15 %, about three times larger than the NNLO+NNLL calculation uncertainties quoted above.

3 Data and simulated samples

The ATLAS detector [15] at the LHC covers nearly the entire solid angle around the collision point, and consists of an inner tracking detector surrounded by a thin superconducting solenoid magnet producing a 2 T axial magnetic field, electromagnetic and hadronic calorimeters, and an external muon spectrometer incorporating three large toroid magnet assemblies. The inner detector consists of a high-granularity silicon pixel detector and a silicon microstrip tracker, together providing precision tracking in the pseudorapidityFootnote 2 range \(|\eta |<2.5\), complemented by a transition radiation tracker providing tracking and electron identification information for \(|\eta |<2.0\). A lead/liquid-argon (LAr) electromagnetic calorimeter covers the region \(|\eta |<3.2\), and hadronic calorimetry is provided by steel/scintillator tile calorimeters for \(|\eta |<1.7\) and copper/LAr hadronic endcap calorimeters. The forward region is covered by additional LAr calorimeters with copper and tungsten absorbers. The muon spectrometer consists of precision tracking chambers covering the region \(|\eta |<2.7\), and separate trigger chambers covering \(|\eta |<2.4\). A three-level trigger system, using custom hardware followed by two software-based levels, is used to reduce the event rate to about 400 Hz for offline storage.

The analysis was performed on the ATLAS 2011–2012 proton–proton collision data sample, corresponding to integrated luminosities of 4.6 \({\mathrm{fb}^{-1}}\) at \(\sqrt{s}=7~\) TeV and 20.3 \({\mathrm{fb}^{-1}}\) at \(\sqrt{s}=8~\) TeV after the application of detector status and data quality requirements. Events were required to pass either a single-electron or single-muon trigger, with thresholds chosen in each case such that the efficiency plateau is reached for leptons with \(p_\mathrm{T}>25\) GeV passing offline selections. Due to the high instantaneous luminosities achieved by the LHC, each triggered event also includes the signals from on average about 9 (\(\sqrt{s}=7~\) TeV) or 20 (\(\sqrt{s}=8~\) TeV) additional inelastic \(pp\) collisions in the same bunch crossing (known as pileup).

Monte Carlo simulated event samples were used to develop the analysis, to compare to the data and to evaluate signal and background efficiencies and uncertainties. Samples were processed either through the full ATLAS detector simulation [16] based on GEANT4 [17], or through a faster simulation making use of parameterised showers in the calorimeters [18]. Additional simulated \(pp\) collisions generated either with Pythia6 [19] (for \({\sqrt{s}\,{=}\,7~\mathrm{TeV}}\) simulation) or Pythia8 [20] (for \({\sqrt{s}=8~\mathrm{TeV}}\)) were overlaid to simulate the effects of both in- and out-of-time pileup, from additional \(pp\) collisions in the same and nearby bunch crossings. All simulated events were then processed using the same reconstruction algorithms and analysis chain as the data. Small corrections were applied to lepton trigger and selection efficiencies to better model the performance seen in data, as discussed further in Sect. 6.

The baseline \(t\overline{t}\) full simulation sample was produced using the NLO matrix element generator Powheg [2123] interfaced to Pythia6 [19] with the Perugia 2011C tune (P2011C) [24] for parton shower, fragmentation and underlying event modelling, and CT10 PDFs [11], and included all \(t\overline{t}\) final states involving at least one lepton. The \(W\rightarrow \ell \nu \) branching ratio was set to the SM expectation of \(0.1082\) [25], and \({m}_{t}\) was set to 172.5 GeV. Alternative \(t\overline{t}\) samples were produced with the NLO generator MC@NLO [26, 27] interfaced to Herwig [28] with Jimmy [29] for the underlying event modelling, with the ATLAS AUET2 [30] tune and CT10 PDFs; and with the leading-order (LO) multileg generator Alpgen [31] interfaced to either Pythia6 or Herwig and Jimmy, with the CTEQ6L1 PDFs [32]. These samples were all normalised to the NNLO+NNLL cross-section predictions given in Sect. 2 when comparing simulation with data.

Backgrounds were classified into two types: those with two real prompt leptons from \(W\) or \(Z\) boson decays (including those produced via leptonic \(\tau \) decays), and those where at least one of the reconstructed lepton candidates is misidentified, i.e. a non-prompt lepton from the decay of a bottom or charm hadron, an electron from a photon conversion, hadronic jet activity misidentified as an electron, or a muon produced from an in-flight decay of a pion or kaon. The first category with two prompt leptons includes \(Wt\) single top production, modelled using Powheg + Pythia6 [33] with the CT10 PDFs and the P2011C tune; \(Z\rightarrow \tau \tau \)+jets modelled using Alpgen + Herwig + Jimmy (\(\sqrt{s}=7~\) TeV) or Alpgen + Pythia6 including LO matrix elements for \(Z{b\overline{b}}\) production, with CTEQ6L1 PDFs; and diboson (\(WW\), \(WZ\), \(ZZ\)) production in association with jets, modelled using Alpgen + Herwig + Jimmy. The \(Wt\) background was normalised to approximate NNLO cross-sections of \(15.7\pm 1.2\) pb at \({\sqrt{s}\,{=}\,7~\mathrm{TeV}}\) and \(22.4\pm 1.5\) pb at \({\sqrt{s}=8~\mathrm{TeV}}\), determined as in Ref. [34]. The inclusive \(Z\) cross-sections were set to the NNLO predictions from FEWZ [35], but the normalisation of \(Z\rightarrow \tau \tau \rightarrow e\mu 4\nu \) backgrounds with \(b\)-tagged jets were determined from data as described in Sect. 5.1. The diboson background was normalised to the NLO QCD inclusive cross-section predictions calculated with MCFM [36]. Production of \(t\overline{t}\) in association with a \(W\) or \(Z\) boson, which contributes to the sample with same-sign leptons, was simulated with Madgraph [37] interfaced to Pythia with CTEQ6L1 PDFs, and normalised to NLO cross-section predictions [38, 39].

Backgrounds with one real and one misidentified lepton include \(t\overline{t}\) events with one hadronically decaying \(W\); \(W\)+jets production, modelled as for \(Z\)+jets; \(W\gamma \)+jets, modelled with Sherpa [40] with CT10 PDFs; and \(t\)-channel single top production, modelled using AcerMC [41] interfaced to Pythia6 with CTEQ6L1 PDFs. Other backgrounds, including processes with two misidentified leptons, are negligible after the event selections used in this analysis.

4 Object and event selection

The analysis makes use of reconstructed electrons, muons and \(b\)-tagged jets. Electron candidates were reconstructed from an isolated electromagnetic calorimeter energy deposit matched to an inner detector track and passing tight identification requirements [42], with transverse energy \(E_\mathrm{T}>25\) GeV and pseudorapidity \(|\eta |<2.47\). Electron candidates within the transition region between the barrel and endcap electromagnetic calorimeters, \(1.37<|\eta |<1.52\), were removed. Isolation requirements were used to reduce background from non-prompt electrons. The calorimeter transverse energy within a cone of size \(\varDelta R=0.2\) and the scalar sum of track \(p_\mathrm{T}\) within a cone of size \(\varDelta R=0.3\), in each case excluding the contribution from the electron itself, were each required to be smaller than \(E_\mathrm{T}\) and \(\eta \)-dependent thresholds calibrated to separately give nominal selection efficiencies of 98 % for prompt electrons from \(Z\rightarrow ee\) decays.

Muon candidates were reconstructed by combining matching tracks reconstructed in both the inner detector and muon spectrometer [43], and were required to satisfy \(p_\mathrm{T}>25\) GeV and \(|\eta |<2.5\). In the \(\sqrt{s}=7~\) TeV dataset, the calorimeter transverse energy within a cone of size \(\varDelta R=0.2\), excluding the energy deposited by the muon, was required to be less than 4 GeV, and the scalar sum of track \(p_\mathrm{T}\) within a cone of size \(\varDelta R=0.3\), excluding the muon track, was required to be less than 2.5 GeV. In the \(\sqrt{s}=8~\) TeV dataset, these isolation requirements were replaced by a cut \(I<0.05\), where \(I\) is the ratio of the sum of track \(p_\mathrm{T}\) in a variable-sized cone of radius \(\varDelta R=10~\mathrm{~{\hbox {GeV}}}/p_\mathrm{T}^\mu \) to the transverse momentum \(p_\mathrm{T}^\mu \) of the muon [44]. Both sets of isolation requirements have efficiencies of about 97 % for prompt muons from \(Z\rightarrow \mu \mu \) decays.

Jets were reconstructed using the anti-\(k_t\) algorithm [45, 46] with radius parameter \(R=0.4\), starting from calorimeter energy clusters calibrated at the electromagnetic energy scale for the \(\sqrt{s}=7~\) TeV dataset, or using the local cluster weighting method for \({\sqrt{s}=8~\mathrm{TeV}}\) [47]. Jets were calibrated using an energy- and \(\eta \)-dependent simulation-based calibration scheme, with in-situ corrections based on data, and were required to satisfy \(p_\mathrm{T}>25\) GeV and \(|\eta |<2.5\). To suppress the contribution from low-\(p_\mathrm{T}\) jets originating from pileup interactions, a jet vertex fraction requirement was applied: at \({\sqrt{s}\,{=}\,7~\mathrm{TeV}}\) jets were required to have at least 75 % of the scalar sum of the \(p_\mathrm{T}\) of tracks associated with the jet coming from tracks associated with the event primary vertex. The latter was defined as the reconstructed vertex with the highest sum of associated track \(p_\mathrm{T}^2\). Motivated by the higher pileup background, in the \(\sqrt{s}=8~\) TeV dataset this requirement was loosened to 50 %, only applied to jets with \(p_\mathrm{T}<50\) GeV and \(|\eta |<2.4\), and the effects of pileup on the jet energy calibration were further reduced using the jet-area method as described in Ref. [48]. Finally, to further suppress non-isolated leptons likely to have come from heavy-flavour decays inside jets, electrons and muons within \(\varDelta R=0.4\) of selected jets were also discarded.

Jets were \(b\)-tagged as likely to have originated from \(b\) quarks using the MV1 algorithm, a multivariate discriminant making use of track impact parameters and reconstructed secondary vertices [49, 50]. Jets were defined to be \(b\)-tagged if the MV1 discriminant value was larger than a threshold (working point) corresponding approximately to a 70 % efficiency for tagging \(b\)-quark jets from top decays in \(t\overline{t}\) events, with a rejection factor of about 140 against light-quark and gluon jets, and about five against jets originating from charm quarks.

Events were required to have at least one reconstructed primary vertex with at least five associated tracks, and no jets failing jet quality and timing requirements. Events with muons compatible with cosmic-ray interactions and muons losing substantial fractions of their energy through bremsstrahlung in the detector material were also removed. A preselection requiring exactly one electron and one muon selected as described above was then applied, with at least one of the leptons being matched to an electron or muon object triggering the event. Events with an opposite-sign \(e\mu \) pair constituted the main analysis sample, whilst events with a same-sign \(e\mu \) pair were used in the estimation of the background from misidentified leptons.

5 Extraction of the \(t\overline{t}\) cross-section

The \(t\overline{t}\) production cross-section \({\sigma _{t\overline{t}}}\) was determined by counting the numbers of opposite-sign \(e\mu \) events with exactly one (\(N_1\)) and exactly two (\(N_2\)) \(b\)-tagged jets. No requirements were made on the number of untagged jets; such jets originate from \(b\)-jets from top decays which were not tagged, and light-quark, charm-quark or gluon jets from QCD radiation. The two event counts can be expressed as:

$$\begin{aligned}&\!\!\!N_1= L {{\sigma _{t\overline{t}}}}\, {{\epsilon _{e\mu }}} {2{\epsilon _{b}}} (1-{{{C}_{b}}}{{\epsilon _{b}}}) + {N_1^\mathrm{bkg}}\nonumber \\&\!\!\!N_2= L {{\sigma _{t\overline{t}}}}\, {{\epsilon _{e\mu }}} {{{C}_{b}}} {{\epsilon _{b}}}^{2} + {N_2^\mathrm{bkg}}\end{aligned}$$
(1)

where \(L\) is the integrated luminosity of the sample, \({\epsilon _{e\mu }}\) is the efficiency for a \(t\overline{t}\) event to pass the opposite-sign \(e\mu \) preselection and \({{C}_{b}}\) is a tagging correlation coefficient close to unity. The combined probability for a jet from the quark \(q\) in the \(t\rightarrow Wq\) decay to fall within the acceptance of the detector, be reconstructed as a jet with transverse momentum above the selection threshold, and be tagged as a \(b\)-jet, is denoted by \({\epsilon _{b}}\). Although this quark is almost always a \(b\) quark, \({\epsilon _{b}}\) thus also accounts for the approximately \(0.2~\%\) of top quarks that decay to \(Ws\) or \(Wd\) rather than \(Wb\), slightly reducing the effective \(b\)-tagging efficiency. Furthermore, the value of \({\epsilon _{b}}\) is slightly increased by the small contributions to \(N_1\) and \(N_2\) from mistagged light-quark, charm-quark or gluon jets from radiation in \(t\overline{t}\) events, although more than 98 % of the tagged jets are expected to contain particles from \(B\)-hadron decays in both the one and two \(b\)-tag samples.

Table 1 Observed numbers of opposite-sign \(e\mu \) events with one and two \(b\)-tagged jets (\(N_1\) and \(N_2\)) for each data sample, together with the estimates of backgrounds and associated total uncertainties described in Sect. 6

If the decays of the two top quarks and the subsequent reconstruction of the two \(b\)-tagged jets are completely independent, the probability to tag both \(b\)-jets \({\epsilon _{bb}}\) is given by \({\epsilon _{bb}}={\epsilon _{b}}^2\). In practice, small correlations are present for both kinematic and instrumental reasons, and these are taken into account via the tagging correlation \({{C}_{b}}\), defined as \({{C}_{b}}={\epsilon _{bb}}/{\epsilon _{b}}^2\) or equivalently \({{C}_{b}}=4 N^{t\overline{t}}_{e\mu } N_2^{t\overline{t}}/(N^{t\overline{t}}_1+2 N^{t\overline{t}}_2)^2\), where \(N^{t\overline{t}}_{e\mu }\) is the number of preselected \(e\mu \) \(t\overline{t}\) events and \(N^{t\overline{t}}_1\) and \(N^{t\overline{t}}_2\) are the numbers of \(t\overline{t}\) events with one and two \(b\)-tagged jets. Values of \({{C}_{b}}\) greater than one correspond to a positive correlation, where a second jet is more likely to be selected if the first one is already selected, whilst \({{C}_{b}}=1\) corresponds to no correlation. This correlation term also compensates for the effect on \({\epsilon _{b}}\), \(N_1\) and \(N_2\) of the small number of mistagged charm-quark or gluon jets from radiation in the \(t\overline{t}\) events.

Fig. 1
figure 1

Distributions of the number of \(b\)-tagged jets in preselected opposite-sign \(e\mu \) events in a \({\sqrt{s}\,{=}\,7~\mathrm{TeV}}\) and b \({\sqrt{s}=8~\mathrm{TeV}}\) data. The data are shown compared to the expectation from simulation, broken down into contributions from \(t\overline{t}\), \(Wt\) single top, \(Z\)+jets, dibosons, and events with misidentified electrons or muons, normalised to the same integrated luminosity as the data. The lower parts of the figure show the ratios of simulation to data, using various \(t\overline{t}\) signal samples generated with Powheg + Pythia6 (PY), MC@NLO + Herwig (HW) and Alpgen + Herwig, and with the cyan band indicating the statistical uncertainty

Background from sources other than \(t\overline{t}\rightarrow e\mu \nu {\overline{\nu }}{b\overline{b}}\) also contributes to the event counts \(N_1\) and \(N_2\), and is given by the terms \({N_1^\mathrm{bkg}}\) and \({N_2^\mathrm{bkg}}\). The preselection efficiency \({\epsilon _{e\mu }}\) and tagging correlation \({{C}_{b}}\) were taken from \(t\overline{t}\) event simulation, and the background contributions \({N_1^\mathrm{bkg}}\) and \({N_2^\mathrm{bkg}}\) were estimated using a combination of simulation- and data-based methods, allowing the two equations in Eq. (1) to be solved numerically yielding \({\sigma _{t\overline{t}}}\) and \({\epsilon _{b}}\).

A total of 11796 events passed the \(e\mu \) opposite-sign preselection in \({\sqrt{s}\,{=}\,7~\mathrm{TeV}}\) data, and 66453 in \({\sqrt{s}=8~\mathrm{TeV}}\) data. Table 1 shows the number of events with one and two \(b\)-tagged jets, together with the estimates of non-\(t\overline{t}\) background and their systematic uncertainties discussed in detail in Sect. 5.1 below. The samples with one \(b\)-tagged jet are expected to be about 89 % pure in \(t\overline{t}\) events, with the dominant background coming from \(Wt\) single top production, and smaller contributions from events with misidentified leptons, \(Z\)+jets and dibosons. The samples with two \(b\)-tagged jets are expected to be about 96 % pure in \(t\overline{t}\) events, with \(Wt\) production again being the dominant background.

Distributions of the number of \(b\)-tagged jets in opposite-sign \(e\mu \) events are shown in Fig. 1, and compared to the expectations with several \(t\overline{t}\) simulation samples. The histogram bins with one and two \(b\)-tagged jets correspond to the data event counts shown in Table 1. Distributions of the number of jets, the \(b\)-tagged jet \(p_\mathrm{T}\), and the electron and muon \(|\eta |\) and \(p_\mathrm{T}\) are shown for opposite-sign \(e\mu \) events with at least one \(b\)-tagged jet in Fig. 2 (\(\sqrt{s}=7~\) TeV) and Fig. 3 (\(\sqrt{s}=8~\) TeV), with the simulation normalised to the same number of entries as the data. The lepton \(|\eta |\) distributions reflect the differing acceptances and efficiencies for electrons and muons, in particular the calorimeter transition region at \(1.37<|\eta |<1.52\). In general, the agreement between data and simulation is good, within the range of predictions from the different \(t\overline{t}\) simulation samples.

The value of \({\sigma _{t\overline{t}}}\) extracted from Eq. (1) is inversely proportional to the assumed value of \({\epsilon _{e\mu }}\), with \((\mathrm{d}{\sigma _{t\overline{t}}}/\mathrm{d}{\epsilon _{e\mu }})\)/\(({\sigma _{t\overline{t}}}/{\epsilon _{e\mu }})=-1\). Uncertainties on \({\epsilon _{e\mu }}\) therefore translate directly into uncertainties on \({\sigma _{t\overline{t}}}\). The value of \({\epsilon _{e\mu }}\) was determined from simulation to be about 0.8 % for both centre-of-mass energies, and includes the \(t\overline{t}\rightarrow e\mu \nu {\overline{\nu }}{b\overline{b}}\) branching ratio of about 3.2 % including \(W\rightarrow \tau \rightarrow e/\mu \) decays. Similarly, \({\sigma _{t\overline{t}}}\) is proportional to the value of \({{C}_{b}}\), also determined from simulation, giving a dependence with the opposite sign, \((\mathrm{d}{\sigma _{t\overline{t}}}\) /\( \mathrm{d}{{C}_{b}})/({\sigma _{t\overline{t}}}/{{C}_{b}})=1\). The systematic uncertainties on \({\epsilon _{e\mu }}\) and \({{C}_{b}}\) are discussed in Sect. 6.

With the kinematic cuts and \(b\)-tagging working point chosen for this analysis, the sensitivities of \({\sigma _{t\overline{t}}}\) to knowledge of the backgrounds \({N_1^\mathrm{bkg}}\) and \({N_2^\mathrm{bkg}}\) are given by \((\mathrm{d}{\sigma _{t\overline{t}}}/\mathrm{d}{N_1^\mathrm{bkg}})\)/\(({\sigma _{t\overline{t}}}/{N_1^\mathrm{bkg}})\,{=}-0.12\) and \((\mathrm{d}{\sigma _{t\overline{t}}}/\mathrm{d}{N_2^\mathrm{bkg}})\)/\(({\sigma _{t\overline{t}}}/{N_2^\mathrm{bkg}})=-0.004\). The fitted cross-sections are therefore most sensitive to the systematic uncertainties on \({N_1^\mathrm{bkg}}\), whilst for the chosen \(b\)-tagging working point, the measurements of \(N_2\) serve mainly to constrain \({\epsilon _{b}}\). As discussed in Sect. 6.1, consistent results were also obtained at different \(b\)-tagging efficiency working points that induce greater sensitivity to the background estimate in the two \(b\)-tag sample.

5.1 Background estimation

The \(Wt\) single top and diboson backgrounds were estimated from simulation as discussed in Sect. 3. The \(Z\)+jets background (with \(Z\rightarrow \tau \tau \rightarrow e\mu 4\nu \)) at \({\sqrt{s}=8~\mathrm{TeV}}\) was estimated from simulation using Alpgen+Pythia, scaled by the ratios of \(Z\rightarrow ee\) or \(Z\rightarrow \mu \mu \) accompanied by \(b\)-tagged jets measured in data and simulation. The ratio was evaluated separately in the one and two \(b\)-tag event samples. This scaling eliminates uncertainties due to the simulation modelling of jets (especially heavy-flavour jets) produced in association with the \(Z\) bosons. The data-to-simulation ratios were measured in events with exactly two opposite-sign electrons or muons passing the selections given in Sect. 4 and one or two \(b\)-tagged jets, by fitting the dilepton invariant mass distributions in the range 60–120 GeV, accounting for the backgrounds from \(t\overline{t}\) production and misidentified leptons. Combining the results from both dilepton channels, the scale factors were determined to be \(1.43\pm 0.03\) and \(1.13\pm 0.08\) for the one and two \(b\)-tag backgrounds, after normalising the simulation to the inclusive \(Z\) cross-section prediction from FEWZ [35]. The uncertainties include systematic components derived from a comparison of results from the \(ee\) and \(\mu \mu \) channels, and from studying the variation of scale factors with \(Z\) boson \(p_\mathrm{T}\). The average \(p_\mathrm{T}\) is higher in selected \(Z\rightarrow \tau \tau \rightarrow e\mu 4\nu \) events than in \(Z\rightarrow ee/\mu \mu \) events due to the momentum lost to the undetected neutrinos from the \(\tau \) decays. The same procedure was used for the \({\sqrt{s}\,{=}\,7~\mathrm{TeV}}\) dataset, resulting in scale factors of \(1.23\pm 0.07\) (one \(b\)-tag) and \(1.14\pm 0.18\) (two \(b\)-tags) for the Alpgen + Herwig \(Z\)+jets simulation, which predicts different numbers of events with heavy-flavour jets than Alpgen + Pythia.

The background from events with one real and one misidentified lepton was estimated using a combination of data and simulation. Simulation studies show that the samples with a same-sign \(e\mu \) pair and one or two \(b\)-tagged jets are dominated by events with misidentified leptons, with rates comparable to those in the opposite-sign sample. The contributions of events with misidentified leptons were therefore estimated using the same-sign event counts in data after subtraction of the estimated prompt same-sign contributions, multiplied by the opposite- to same-sign misidentified-lepton ratios \({{R_j} = N_j^{\mathrm{mis,OS}}}/ {N_j^{\mathrm{mis,SS}}}\) estimated from simulation for events with \(j=1\) and 2 \(b\)-tagged jets. The procedure is illustrated by Table 2, which shows the expected numbers of events with misidentified leptons in opposite- and same-sign samples. The contributions where the electron is misidentified, coming from a photon conversion, the decay of a heavy-flavour hadron or other sources (such as a misidentified hadron within a jet), and where the muon is misidentified, coming either from heavy-flavour decay or other sources (e.g. decay in flight of a pion or kaon) are shown separately. The largest contributions come from photon conversions giving electron candidates, and most of these come from photons radiated from prompt electrons produced from \(t\rightarrow Wq\rightarrow e\nu q\) in signal \(t\overline{t}\rightarrow e\mu \nu {\overline{\nu }}{b\overline{b}}\) events. Such electrons populate both the opposite- and same-sign samples, and are treated as misidentified-lepton background.

The ratios \(R_j\) were estimated from simulation to be \(R_1=1.4\pm 0.5\) and \(R_2=1.1\pm 0.5\) at \({\sqrt{s}\,{=}\,7~\mathrm{TeV}}\), and \(R_1=1.2\pm 0.3\) and \(R_2=1.6\pm 0.5\) at \({\sqrt{s}=8~\mathrm{TeV}}\). The uncertainties were derived by considering the range of \(R_j\) values for different components of the misidentified-lepton background, including the small contributions from sources other than photon conversions and heavy-flavour decays, which do not significantly populate the same-sign samples. As shown in Table 2, about 25 % of the same-sign events have two prompt leptons, which come mainly from semileptonic \(t\overline{t}\) events with an additional leptonically decaying \(W\) or \(Z\) boson, diboson decays producing two same-sign leptons, and wrong-sign \(t\overline{t}\rightarrow e\mu \nu {\overline{\nu }}{b\overline{b}}\) events where the electron charge was misreconstructed. A conservative uncertainty of 50 % was assigned to this background, based on studies of the simulation modelling of electron charge misidentification [42] and uncertainties in the rates of contributing physics processes.

The simulation modelling of the different components of the misidentified-lepton background was checked by studying kinematic distributions of same-sign events, as illustrated for the \(|\eta |\) and \(p_\mathrm{T}\) distributions of the leptons in \({\sqrt{s}=8~\mathrm{TeV}}\) data in Fig. 4. The simulation generally models the normalisation and shapes of distributions well in both the one and two \(b\)-tag event samples. The simulation modelling was further tested in control samples with relaxed electron or muon isolation requirements to enhance the relative contributions of electrons or muons from heavy-flavour decays, and similar levels of agreement were observed.

Fig. 2
figure 2

Distributions of a the number of jets, b the transverse momentum \(p_\mathrm{T}\) of the \(b\)-tagged jets, c the \(|\eta |\) of the electron, d the \(p_\mathrm{T}\) of the electron, e the \(|\eta |\) of the muon and f the \(p_\mathrm{T}\) of the muon, in events with an opposite-sign \(e\mu \) pair and at least one \(b\)-tagged jet. The \({\sqrt{s}\,{=}\,7~\mathrm{TeV}}\) data are compared to the expectation from simulation, broken down into contributions from \(t\overline{t}\), single top, \(Z\)+jets, dibosons, and events with misidentified electrons or muons, normalised to the same number of entries as the data. The lower parts of the figure show the ratios of simulation to data, using various \(t\overline{t}\) signal samples and with the cyan band indicating the statistical uncertainty. The last bin includes the overflow

Fig. 3
figure 3

Distributions of a the number of jets, b the transverse momentum \(p_\mathrm{T}\) of the \(b\)-tagged jets, c the \(|\eta |\) of the electron, d the \(p_\mathrm{T}\) of the electron, e the \(|\eta |\) of the muon and f the \(p_\mathrm{T}\) of the muon, in events with an opposite-sign \(e\mu \) pair and at least one \(b\)-tagged jet. The \({\sqrt{s}=8~\mathrm{TeV}}\) data are compared to the expectation from simulation, broken down into contributions from \(t\overline{t}\), single top, \(Z\)+jets, dibosons, and events with misidentified electrons or muons, normalised to the same number of entries as the data. The lower parts of the figure show the ratios of simulation to data, using various \(t\overline{t}\) signal samples and with the cyan band indicating the statistical uncertainty. The last bin includes the overflow

Table 2 Breakdown of estimated misidentified-lepton contributions to the one (\(1b\)) and two (\(2b\)) \(b\)-tag opposite- and same-sign (OS and SS) \(e\mu \) event samples at \({\sqrt{s}\,{=}\,7~\mathrm{TeV}}\) and \({\sqrt{s}=8~\mathrm{TeV}}\). The different misidentified-lepton categories are described in the text. For the same-sign samples, the contributions from wrong-sign (where the electron charge sign is misreconstructed) and right-sign prompt lepton events are also shown, and the total expectations are compared to the data. The uncertainties shown are due to the limited size of the simulated samples, and values and uncertainties quoted as ‘0.0’ are smaller than 0.05
Fig. 4
figure 4

Distributions of electron and muon \(|\eta |\) and \(p_\mathrm{T}\) in same-sign \(e\mu \) events at \({\sqrt{s}=8~\mathrm{TeV}}\) with at least one \(b\)-tagged jet. The simulation prediction is normalised to the same integrated luminosity as the data, and broken down into contributions where both leptons are prompt, or one is a misidentified lepton from a photon conversion originating from a top quark decay or from background, or from heavy-flavour decay. In the \(p_\mathrm{T}\) distributions, the last bin includes the overflows

Table 3 Detailed breakdown of the symmetrised relative statistical, systematic and total uncertainties on the measurements of the \(t\overline{t}\) production cross-section \({\sigma _{t\overline{t}}}\) at \({\sqrt{s}\,{=}\,7~\mathrm{TeV}}\) and \({\sqrt{s}=8~\mathrm{TeV}}\). Uncertainties quoted as ‘0.00’ are smaller than 0.005, whilst ‘–’ indicates the corresponding uncertainty is not applicable. The uncertainties on \({\epsilon _{e\mu }}\) and \({{C}_{b}}\) are also shown, with their relative signs indicated where relevant. They contribute with opposite signs to the uncertainties on \({\sigma _{t\overline{t}}}\), which also include uncertainties from estimates of the background terms \({N_1^\mathrm{bkg}}\) and \({N_2^\mathrm{bkg}}\). The lower part of the table gives the systematic uncertainties that are different for the measurement of the fiducial cross-section \(\sigma _{t\overline{t}}^\mathrm{fid}\), together with the total analysis systematic and total uncertainties on \(\sigma _{t\overline{t}}^\mathrm{fid}\)

6 Systematic uncertainties

The systematic uncertainties on the measured cross-sections \({\sigma _{t\overline{t}}}\) are shown in detail in Table 3 together with the individual uncertainties on \({\epsilon _{e\mu }}\) and \({{C}_{b}}\). A summary of the uncertainties on \({\sigma _{t\overline{t}}}\) is shown in Table 4. Each source of uncertainty was evaluated by repeatedly solving Eq. (1) with all relevant input parameters simultaneously changed by \(\pm 1\) standard deviation. Systematic correlations between input parameters (in particular significant anti-correlations between \({\epsilon _{e\mu }}\) and \({{C}_{b}}\) which contribute with opposite signs to \({\sigma _{t\overline{t}}}\)) were thus taken into account. The total uncertainties on \({\sigma _{t\overline{t}}}\) and \({\epsilon _{b}}\) were calculated by adding the effects of all the individual systematic components in quadrature, assuming them to be independent. The sources of systematic uncertainty are discussed in more detail below; unless otherwise stated, the same methodology was used for both \(\sqrt{s}=7~\) TeV and \(\sqrt{s}=8~\) TeV datasets.

  • \(t\overline{t}\) modelling: Uncertainties on \({\epsilon _{e\mu }}\) and \({{C}_{b}}\) due to the simulation of \(t\overline{t}\) events were assessed by studying the predictions of different \(t\overline{t}\) generators and hadronisation models as detailed in Sect. 3. The prediction for \({\epsilon _{e\mu }}\) was found to be particularly sensitive to the amount of hadronic activity near the leptons, which strongly affects the efficiency of the lepton isolation requirements described in Sect. 4. These isolation efficiencies were therefore measured directly from data, as discussed below. The remaining uncertainties on \({\epsilon _{e\mu }}\) relating to lepton reconstruction, identification and lepton–jet overlap removal, were evaluated from the differences between the predictions from the baseline Powheg + Pythia \(t\overline{t}\) sample and a sample generated using MC@NLO + Herwig, thus varying both the hard-scattering event generator and the fragmentation and hadronisation model. The MC@NLO + Herwig sample gave a larger value of \({\epsilon _{e\mu }}\) but a smaller value of \({{C}_{b}}\). Additional comparisons of Powheg + Pythia samples with the AUET2 rather than P2011C tune and with Powheg + Herwig, i.e. changing only the fragmentation/hadronisation model, gave smaller uncertainties. The Alpgen + Herwig and Alpgen + Pythia samples gave values of \({\epsilon _{e\mu }}\) up to 2 % higher than that of Powheg+Pythia, due largely to a more central predicted \(\eta \) distribution for the leptons. However, this sample uses a leading-order generator and PDFs, and gives an inferior description of the electron and muon \(\eta \) distributions (see Fig. 3c, e), so was not used to set the systematic uncertainty on \({\epsilon _{e\mu }}\). In contrast, the Alpgen samples were considered in setting the uncertainty on \({{C}_{b}}\), taken as the largest difference between the predictions of Powheg + Pythia and any of the other generators. The effect of extra radiation in \(t\overline{t}\) events was also considered explicitly by using pairs of simulation samples with different Pythia tunes whose parameters span the variations compatible with ATLAS studies of additional jet activity in \(t\overline{t}\) events at \({\sqrt{s}\,{=}\,7~\mathrm{TeV}}\) [51], generated using both AcerMC + Pythia and Alpgen + Pythia. These samples predicted large variations in the lepton isolation efficiencies (which were instead measured from data), but residual variations in other lepton-related uncertainties and \({{C}_{b}}\) within the uncertainties set from other simulation samples.

  • Parton distribution functions: The uncertainties on \({\epsilon _{e\mu }}\), \({{C}_{b}}\) and the \(Wt\) single top background due to uncertainties on the proton PDFs were evaluated using the error sets of the CT10 NLO [11], MSTW 2008 68 % CL NLO [9, 10] and NNPDF 2.3 NLO [13] sets. The final uncertainty was calculated as half the envelope encompassing the predictions from all three PDF sets along with their associated uncertainties, following the PDF4LHC recommendations [8].

  • QCD scale choices: The lepton \(p_\mathrm{T}\) and \(\eta \) distributions, and hence \({\epsilon _{e\mu }}\), are sensitive to the choices of QCD renormalisation and factorisation scales. This effect was investigated using \(\sqrt{s}=8~\) TeV generator-level Powheg + Pythia \(t\overline{t}\) samples where the two scales were separately varied up and down by a factor of two from their default values of \(Q^2=\) \({m}_{t}\) \(^2+p_{\mathrm{T},t}^2\). The systematic uncertainty for each scale was taken as half the difference in \({\epsilon _{e\mu }}\) values between the samples with increased and decreased QCD scale, and the uncertainties for the renormalisation and factorisation scales were then added linearly to give a total scale uncertainty of 0.30 % on \({\epsilon _{e\mu }}\), assumed to be valid for both centre-of-mass energies.

  • Single top modelling: Uncertainties related to \(Wt\) single top modelling were assessed by comparing the predictions from Powheg + Pythia, Powheg + Herwig, MC@NLO + Herwig, and AcerMC + Pythia with two tunes producing different amounts of additional radiation, in all cases normalising the total production rate to the approximate NNLO cross-section prediction. The resulting uncertainties are about 5 % and 20 % on the one and two \(b\)-tag background contributions. The background in the two \(b\)-tag sample is sensitive to the production of \(Wt\) with an additional \(b\)-jet, a NLO contribution to \(Wt\) which can interfere with the \(t\overline{t}\) final state. The sensitivity to this interference was studied by comparing the predictions of Powheg with the diagram-removal (baseline) and diagram-subtraction schemes [33, 52], giving additional single-top/\(t\overline{t}\) interference uncertainties of 1–2 % and 20 % for the one and two \(b\)-tag samples. The production of single top quarks in association with a \(Z\) boson gives contributions which are negligible compared to the above uncertainties. Production of single top quarks via the \(t\)- and \(s\)-channels gives rise to final states with only one prompt lepton, and is accounted for as part of the misidentified-lepton background.

  • Background cross-sections: The uncertainties on the \(Wt\) single top cross-section were taken to be 7.6 % and 6.8 % at \({\sqrt{s}\,{=}\,7~\mathrm{TeV}}\) and \({\sqrt{s}=8~\mathrm{TeV}}\), based on Ref. [34]. The uncertainties on the diboson cross-sections were set to 5 % [36].

  • Diboson modelling: Uncertainties in the backgrounds from dibosons with one or two additional \(b\)-tagged jets were assessed by comparing the baseline prediction from Alpgen + Herwig with that of Sherpa [40] including massive \(b\) and \(c\) quarks, and found to be about 20 %. The background from 125 GeV SM Higgs production in the gluon fusion, vector-boson fusion, and \(WH\) and \(ZH\) associated production modes, with \(H\rightarrow WW\) and \(H\rightarrow \tau \tau \), was evaluated to be smaller than the diboson modelling uncertainties, and was neglected.

  • Z \(+\) jets extrapolation: The uncertainties on the extrapolation of the \(Z\)+jets background from \(Z\rightarrow ee/\mu \mu \) to \(Z\rightarrow \tau \tau \) events result from statistical uncertainties, comparing the results from \(ee\) and \(\mu \mu \), which have different background compositions, and considering the dependence of the scale factors on \(Z\) boson \(p_\mathrm{T}\).

  • Lepton identification and measurement: The modelling of the electron and muon identification efficiencies, energy scales and resolutions (including the effects of pileup) were studied using \(Z\rightarrow ee/\mu \mu \), \(J/\psi \rightarrow ee/\mu \mu \) and \(W\rightarrow e\nu \) events in data and simulation, using the techniques described in Refs. [42, 43, 53]. Small corrections were applied to the simulation to better model the performance seen in data, and the associated systematic uncertainties were propagated to the cross-section measurement.

  • Lepton isolation: The efficiency of the lepton isolation requirements was measured directly in data, from the fraction of selected opposite-sign \(e\mu \) events with one or two \(b\)-tags where either the electron or muon fails the isolation cut. The results were corrected for the contamination from misidentified leptons, estimated using the same-sign \(e\mu \) samples as described in Sect. 5, or by using the distributions of lepton impact parameter significance \(|d_0|/\sigma _{d_0}\), where \(d_0\) is the distance of closest approach of the lepton track to the event primary vertex in the transverse plane, and \(\sigma _{d_0}\) its uncertainty. Consistent results were obtained from both methods, and showed that the baseline Powheg+Pythia simulation overestimates the efficiencies of the isolation requirements by about 0.5 % for both the electrons and muons. These corrections were applied to \({\epsilon _{e\mu }}\), with uncertainties dominated by the limited sizes of the same-sign and high impact-parameter significance samples used for background estimation. Similar results were found from studies in \(Z\rightarrow ee\) and \(Z\rightarrow \mu \mu \) events, after correcting the results for the larger average amount of hadronic activity near the leptons in \(t\overline{t}\rightarrow e\mu \nu {\overline{\nu }}{b\overline{b}}\) events.

  • Jet-related uncertainties: Although the efficiency to reconstruct and \(b\)-tag jets from \(t\overline{t}\) events is extracted from the data, uncertainties in the jet energy scale, energy resolution and reconstruction efficiency affect the backgrounds estimated from simulation and the estimate of the tagging correlation \({{C}_{b}}\). They also have a small effect on \({\epsilon _{e\mu }}\) via the lepton–jet \(\varDelta R\) separation cuts. The jet energy scale was varied in simulation according to the uncertainties derived from simulation and in-situ calibration measurements [47, 54], using a model with 21 (\({\sqrt{s}\,{=}\,7~\mathrm{TeV}}\)) or 22 (\({\sqrt{s}=8~\mathrm{TeV}}\)) separate orthogonal uncertainty components which were then added in quadrature. The jet energy resolution was found to be well modelled by simulation [55], and remaining uncertainties were assessed by applying additional smearing, which reduces \({\epsilon _{e\mu }}\). The calorimeter jet reconstruction efficiency was measured in data using track-based jets, and is also well described by the simulation; the impact of residual uncertainties was assessed by randomly discarding jets. The uncertainty associated with the jet vertex fraction requirement was assessed from studies of \(Z\rightarrow ee/\mu \mu \)+jets events.

  • b -tagging uncertainties:   The efficiency for \(b\)-tagging jets from \(t\overline{t}\) events was extracted from the data via Eq. (1), but simulation was used to predict the number of \(b\)-tagged jets and mistagged light-quark, gluon and charm jets in the \(Wt\) single top and diboson backgrounds. The tagging correlation \({{C}_{b}}\) is also slightly sensitive to the efficiencies for tagging heavy- and light-flavour jets. The uncertainties in the simulation modelling of the \(b\)-tagging performance were assessed using studies of \(b\)-jets containing muons [50, 56], jets containing \(D^{*+}\) mesons [57] and inclusive jet events [58].

  • Misidentified leptons: The uncertainties on the number of events with misidentified leptons in the one and two \(b\)-tagged samples were derived from the statistical uncertainties on the numbers of same-sign lepton events, the systematic uncertainties on the opposite- to same-sign ratios \(R_j\), and the uncertainties on the numbers of prompt same-sign events, as discussed in detail in Sect. 5.1. The overall uncertainties on the numbers of misidentified leptons vary from 30 to 50 %, dominated by the uncertainties on the ratios \(R_j\).

  • Integrated luminosity: The uncertainty on the integrated luminosity of the \({\sqrt{s}\,{=}\,7~\mathrm{TeV}}\) dataset is 1.8 % [59]. Using beam-separation scans performed in November 2012, the same methodology was applied to determine the \({\sqrt{s}=8~\mathrm{TeV}}\) luminosity scale, resulting in an uncertainty of 2.8 %. These uncertainties are dominated by effects specific to each dataset, and so are considered to be uncorrelated between the two centre-of-mass energies. The relative uncertainties on the cross-section measurements are slightly larger than those on the luminosity measurements because the \(Wt\) single top and diboson backgrounds are evaluated from simulation, so are also sensitive to the assumed integrated luminosity.

  • LHC beam energy: The LHC beam energy during the 2012 \(pp\) run was calibrated to be \(0.30\pm 0.66\) % smaller than the nominal value of 4 TeV per beam, using the revolution frequency difference of protons and lead ions during \(p\)+Pb runs in early 2013 [60]. Since this calibration is compatible with the nominal \(\sqrt{s}\) of 8 TeV, no correction was applied to the measured \({\sigma _{t\overline{t}}}\) value. However, an uncertainty of 1.72 %, corresponding to the expected change in \({\sigma _{t\overline{t}}}\) for a 0.66 % change in \(\sqrt{s}\) is quoted separately on the final result. This uncertainty was calculated using top++ 2.0, assuming that the relative change of \({\sigma _{t\overline{t}}}\) for a 0.66 % change in \(\sqrt{s}\) is as predicted by the NNLLO+NNLL calculation. Following Ref. [60], the same relative uncertainty on the LHC beam energy is applied for the \({\sqrt{s}\,{=}\,7~\mathrm{TeV}}\) dataset, giving a slightly larger uncertainty of 1.79 % due to the steeper relative dependence of \({\sigma _{t\overline{t}}}\) on \(\sqrt{s}\) in this region. These uncertainties are much larger than those corresponding to the very small dependence of \({\epsilon _{e\mu }}\) on \(\sqrt{s}\), which changes by only 0.5 % between 7 and 8 TeV.

  • Top quark mass: The simulation samples used in this analysis were generated with \({m}_{t}\) \(=172.5\) GeV, but the acceptance for \(t\overline{t}\) and \(Wt\) events, and the \(Wt\) background cross-section itself, depend on the assumed \({m}_{t}\) value. Alternative samples generated with \({m}_{t}\) varied in the range 165–180 GeV were used to quantify these effects. The acceptance and background effects partially cancel, and the final dependence of the result on the assumed \({m}_{t}\) value was determined to be \({\mathrm{d}{\sigma _{t\overline{t}}}/\mathrm{d}{{m}_{t}}}=-0.28~\%/\hbox {GeV}.\) The result of the analysis is reported assuming a fixed top mass of 172.5 GeV, and the small dependence of the cross-section on the assumed mass is not included as a systematic uncertainty.

Table 4 Summary of the relative statistical, systematic and total uncertainties on the measurements of the \(t\overline{t}\) production cross-section \({\sigma _{t\overline{t}}}\) at \({\sqrt{s}\,{=}\,7~\mathrm{TeV}}\) and \({\sqrt{s}=8~\mathrm{TeV}}\)

As shown in Tables 3 and 4, the largest systematic uncertainties on \({\sigma _{t\overline{t}}}\) come from \(t\overline{t}\) modelling and PDFs, and knowledge of the integrated luminosities and LHC beam energy.

Fig. 5
figure 5

Measured \(t\overline{t}\) cross-section at \({\sqrt{s}=8~\mathrm{TeV}}\) as a function of the \(b\)-tagged jet \(p_\mathrm{T}\) cut. The error bars show the uncorrelated part of the statistical uncertainty with respect to the baseline measurement with jet \(p_\mathrm{T}>25\) GeV

6.1 Additional correlation studies

The tagging correlation \({{C}_{b}}\) was determined from simulation to be \(1.009\pm 0.002\pm 0.007\) (\({\sqrt{s}\,{=}\,7~\mathrm{TeV}}\)) and \(1.007\pm 0.002\pm 0.006\) (\({\sqrt{s}=8~\mathrm{TeV}}\)), where the first uncertainty is due to limited sizes of the simulated samples, and the second is dominated by the comparison of predictions from different \(t\overline{t}\) generators. Additional studies were carried out to probe the modelling of possible sources of correlation. One possible source is the production of additional \({b\overline{b}}\) or \({c\overline{c}}\) pairs in \(t\overline{t}\) production, which tends to increase both \({{C}_{b}}\) and the number of events with three or more \(b\)-tagged jets, which are not used in the measurement of \({\sigma _{t\overline{t}}}\). The ratio \(R_{32}\) of events with at least three \(b\)-tagged jets to events with at least two \(b\)-tagged jets was used to quantify this extra heavy-flavour production in data. It was measured to be \(R_{32}\) \(=2.7\pm 0.4\) % (\({\sqrt{s}\,{=}\,7~\mathrm{TeV}}\)) and \(2.8\pm 0.2\) % (\({\sqrt{s}=8~\mathrm{TeV}}\)), where the uncertainties are statistical. These values are close to the Powheg + Pythia prediction of \(2.4\pm 0.1\) % (see Fig. 1), and well within the spread of \(R_{32}\) values seen in the alternative simulation samples.

Kinematic correlations between the two \(b\)-jets produced in the \(t\overline{t}\) decay could also produce a positive tagging correlation, as the efficiency to reconstruct and tag \(b\)-jets is not uniform as a function of \(p_\mathrm{T}\) and \(\eta \). For example, \(t\overline{t}\) pairs produced with high invariant mass tend to give rise to two back-to-back collimated top quark decay systems where both \(b\)-jets have higher than average \(p_\mathrm{T}\), and longitudinal boosts of the \(t\overline{t}\) system along the beamline give rise to \(\eta \) correlations between the two jets. These effects were probed by increasing the jet \(p_\mathrm{T}\) cut in steps from the default of 25 GeV up to 75 GeV; above about 50 GeV, the simulation predicts strong positive correlations of up to \({{C}_{b}}\approx 1.2\) for a 75 GeV \(p_\mathrm{T}\) cut. As shown for the \(\sqrt{s}=8~\) TeV dataset in Fig. 5, the cross-sections fitted in data after taking these correlations into account remain stable across the full \(p_\mathrm{T}\) cut range, suggesting that any such kinematic correlations are well modelled by the simulation. Similar results were seen at \(\sqrt{s}=7~\) TeV. The results were also found to be stable within the uncorrelated components of the statistical and systematic uncertainties when tightening the jet and lepton \(\eta \) cuts, raising the lepton \(p_\mathrm{T}\) cut up to 55 GeV and changing the \(b\)-tagging working point between efficiencies of 60 % and 80 %. No additional uncertainties were assigned as a result of these studies.

7 Results

Combining the estimates of \({\epsilon _{e\mu }}\) and \({{C}_{b}}\) from simulation samples, the estimates of the background \({N_1^\mathrm{bkg}}\) and \({N_2^\mathrm{bkg}}\) shown in Table 1 and the data integrated luminosities, the \(t\overline{t}\) cross-section was determined by solving Eq. (1) to be:

$$\begin{aligned} {\sigma _{t\overline{t}}}&= 182.9\pm 3.1\pm 4.2\pm 3.6\pm 3.3~\mathrm{pb}~ ({\sqrt{s}\,{=}\,7~\mathrm{TeV}})\quad \mathrm{and} \\ {\sigma _{t\overline{t}}}&= 242.4\pm 1.7\pm 5.5\pm 7.5\pm 4.2~\mathrm{pb}~ ({\sqrt{s}=8~\mathrm{TeV}}), \end{aligned}$$

where the four uncertainties arise from data statistics, experimental and theoretical systematic effects related to the analysis, knowledge of the integrated luminosity and of the LHC beam energy. The total uncertainties are 7.1 pb (3.9 %) at \(\sqrt{s}=7~\) TeV and 10.3 pb (4.3 %) at \(\sqrt{s}=8~\) TeV. A detailed breakdown of the different components is given in Table 3. The results are reported for a fixed top quark mass of \({m}_{t}\) \(=172.5\) GeV, and have a dependence on this assumed value of \({\mathrm{d}{\sigma _{t\overline{t}}}/\mathrm{d}{{m}_{t}}}=-0.28\) %/GeV. The product of jet reconstruction and \(b\)-tagging efficiencies \({\epsilon _{b}}\) was measured to be \(0.557\pm 0.009\) at \({\sqrt{s}\,{=}\,7~\mathrm{TeV}}\) and \(0.540\pm 0.006\) at \({\sqrt{s}=8~\mathrm{TeV}}\), in both cases consistent with the values in simulation.

The results are shown graphically as a function of \(\sqrt{s}\) in Fig. 6, together with previous ATLAS measurements of \({\sigma _{t\overline{t}}}\) at \(\sqrt{s}=7~\) TeV in the \(ee\), \(\mu \mu \) and \(e\mu \) dilepton channels using a count of the number of events with two leptons and at least two jets in an 0.7 \({\mathrm{fb}^{-1}}\) dataset [61], and using a fit of jet multiplicities and missing transverse momentum in the \(e\mu \) dilepton channel alone with the full 4.6 \({\mathrm{fb}^{-1}}\) dataset [62]. The \(\sqrt{s}=7~\) TeV results are all consistent, but cannot be combined as they are not based on independent datasets. The measurements from this analysis at both centre-of-mass energies are consistent with the NNLO+NNLL QCD calculations discussed in Sect. 2. The \({\sqrt{s}\,{=}\,7~\mathrm{TeV}}\) result is 13 % higher than a previous measurement by the CMS collaboration [63], whilst the \({\sqrt{s}=8~\mathrm{TeV}}\) result is consistent with that from CMS [64].

Fig. 6
figure 6

Measurements of the \(t\overline{t}\) cross-section at \({\sqrt{s}\,{=}\,7~\mathrm{TeV}}\) and \({\sqrt{s}=8~\mathrm{TeV}}\) from this analysis (\(e\mu \) \(b\)-tag) together with previous ATLAS results at \({\sqrt{s}\,{=}\,7~\mathrm{TeV}}\) using the \(ee\), \(\mu \mu \) and \(e\mu \) channels [61] and using a fit to jet multiplicities and missing transverse momentum in the \(e\mu \) channel [62]. The uncertainties in \(\sqrt{s}\) due to the LHC beam energy uncertainty are displayed as horizontal error bars, and the vertical error bars do not include the corresponding cross-section uncertainties. The three \({\sqrt{s}\,{=}\,7~\mathrm{TeV}}\) measurements are displaced horizontally slightly for clarity. The NNLO+NNLL prediction [6, 7] described in Sect. 2 is also shown as a function of \(\sqrt{s}\), for fixed \({m}_{t}\) \(=172.5~{\hbox {GeV}}\) and with the uncertainties from PDFs, \(\alpha _\mathrm{s}\) and QCD scale choices indicated by the green band

From the present analysis, the ratio of cross-sections \({R_{t\overline{t}}}={\sigma _{t\overline{t}}}\)(8 TeV)/\({\sigma _{t\overline{t}}}\)(7 TeV) was determined to be:

$$\begin{aligned} {R_{t\overline{t}}}=1.326\pm 0.024 \pm 0.015 \pm 0.049 \pm 0.001 \end{aligned}$$

with uncertainties defined as above, adding in quadrature to a total of 0.056. The experimental systematic uncertainties (apart from the statistical components of the lepton isolation and misidentified lepton uncertainties, which were evaluated independently from data in each dataset) and the LHC beam energy uncertainty are correlated between the two centre-of-mass energies. The luminosity uncertainties were taken to be uncorrelated between energies. The result is consistent with the QCD NNLO+NNLL predicted ratio of \(1.430\pm 0.013\) (see Sect. 2), which in addition to the quoted PDF, \(\alpha _\mathrm{s}\) and QCD scale uncertainties varies by only \(\pm 0.001\) for a \(\pm 1\) GeV variation of \({m}_{t}\).

Table 5 Fiducial cross-section measurement results at \({\sqrt{s}\,{=}\,7~\mathrm{TeV}}\) and \({\sqrt{s}=8~\mathrm{TeV}}\), for different requirements on the minimum lepton \(p_\mathrm{T}\) and maximum lepton \(|\eta |\), and with or without the inclusion of leptons from \(W\rightarrow \tau \rightarrow \ell \) decays. In each case, the first uncertainty is statistical, the second due to analysis systematic effects, the third due to the integrated luminosity and the fourth due to the LHC beam energy

7.1 Fiducial cross-sections

The preselection efficiency \({\epsilon _{e\mu }}\) can be written as the product of two terms \({\epsilon _{e\mu }}={{A}_{e\mu }}{{G}_{e\mu }}\), where the acceptance \({A}_{e\mu }\) represents the fraction of \(t\overline{t}\) events which have a true opposite-sign \(e\mu \) pair from \(t\rightarrow W\rightarrow \ell \) decays (including via \(W\rightarrow \tau \rightarrow \ell \)), each with \(p_\mathrm{T}>25\) GeV and within \(|\eta |<2.5\), and \({G}_{e\mu }\) represents the reconstruction efficiency, i.e. the probability that the two leptons are reconstructed and pass all the identification and isolation requirements. A fiducial cross-section \(\sigma _{t\overline{t}}^\mathrm{fid}\) can then be defined as \({\sigma _{t\overline{t}}^\mathrm{fid}}={{A}_{e\mu }}{\sigma _{t\overline{t}}}\), and measured by replacing \({\sigma _{t\overline{t}}}{\epsilon _{e\mu }}\) with \({\sigma _{t\overline{t}}^\mathrm{fid}}{{G}_{e\mu }}\) in Eq. (1), leaving the background terms unchanged. Measurement of the fiducial cross-section avoids the systematic uncertainties associated with \({A}_{e\mu }\), i.e. the extrapolation from the measured lepton phase space to the full phase space populated by inclusive \(t\overline{t}\) production. In this analysis, these come mainly from knowledge of the PDFs and the QCD scale uncertainties. Since the analysis technique naturally corrects for the fraction of jets which are outside the kinematic acceptance through the fitted value of \({\epsilon _{b}}\), no restrictions on jet kinematics are imposed in the definition of \(\sigma _{t\overline{t}}^\mathrm{fid}\). In calculating \({A}_{e\mu }\) and \({G}_{e\mu }\) from the various \(t\overline{t}\) simulation samples, the lepton four-momenta were taken after final-state radiation, and including the four-momenta of any photons within a cone of size \(\varDelta R=0.1\) around the lepton direction, excluding photons from hadron decays or produced in interactions with detector material. The values of \({A}_{e\mu }\) are about 1.4 % (including the \(t\overline{t}\rightarrow e\mu \nu {\overline{\nu }}{b\overline{b}}\) branching ratio), and those of \({G}_{e\mu }\) about 55 %, at both centre-of-mass energies.

The measured fiducial cross-sections at \(\sqrt{s}=7~\) TeV and \(\sqrt{s}=8~\) TeV, for leptons with \(p_\mathrm{T}>25\) GeV and \(|\eta |<2.5\), are shown in the first row of Table 5. The relative uncertainties are shown in the lower part of Table 3; the PDF uncertainties are substantially reduced compared to the inclusive cross-section measurement, and the QCD scale uncertainties are reduced to a negligible level. The \(t\overline{t}\) modelling uncertainties, evaluated from the difference between Powheg+Pythia and MC@NLO+Herwig samples increase slightly, though the differences are not significant given the sizes of the simulated samples. Overall, the analysis systematics on the fiducial cross-sections are 6–11 % smaller than those on the inclusive cross-section measurements.

Simulation studies predict that \(11.9\pm 0.1\) % of \(t\overline{t}\) events in the fiducial region have at least one lepton produced via \(W\rightarrow \tau \rightarrow \ell \) decay. The second row in Table 5 shows the fiducial cross-section measurements scaled down to remove this contribution. The third and fourth rows show the measurements scaled to a different lepton fiducial acceptance of \(p_\mathrm{T}>30\) GeV and \(|\eta |<2.4\), a common phase space accessible to both the ATLAS and CMS experiments.

7.2 Top quark mass determination

The strong dependence of the theoretical prediction for \({\sigma _{t\overline{t}}}\) on \({m}_{t}\) offers the possibility of interpreting measurements of \({\sigma _{t\overline{t}}}\) as measurements of \({m}_{t}\). The theoretical calculations use the pole mass \({m}_{t}^\mathrm{pole}\), corresponding to the definition of the mass of a free particle, whereas the top quark mass measured through direct reconstruction of the top decay products [6568] may differ from the pole mass by \(O(1\) GeV) [69, 70]. It is therefore interesting to compare the values of \({m}_{t}\) determined from the two approaches, as explored previously by the D0 [71, 72] and CMS [73] collaborations.

The dependence of the cross-section predictions (calculated as described in Sect. 2) on \({m}_{t}^\mathrm{pole}\) is shown in Fig. 7 at both \({\sqrt{s}\,{=}\,7~\mathrm{TeV}}\) and \({\sqrt{s}=8~\mathrm{TeV}}\). The calculations were fitted to the parameterisation proposed in Ref. [6], namely:

$$\begin{aligned} {\sigma ^\mathrm{theo}_{t\overline{t}}}({{m}_{t}^\mathrm{pole}})=\sigma ({{m}_{t}^\mathrm{ref}})\left( \frac{{{m}_{t}^\mathrm{ref}}}{{{m}_{t}^\mathrm{pole}}}\right) ^4 (1+a_1x+a_2x^2) \end{aligned}$$
(2)

where the parameterisation constant \({{m}_{t}^\mathrm{ref}}=172.5\) GeV, \(x=({{m}_{t}^\mathrm{pole}}-{{m}_{t}^\mathrm{ref}})/{{m}_{t}^\mathrm{ref}}\), and \(\sigma ({{m}_{t}^\mathrm{ref}})\), \(a_1\) and \(a_2\) are free parameters. This function was used to parameterise the dependence of \({\sigma _{t\overline{t}}}\) on \({m}_{t}\) separately for each of the NNLO PDF sets CT10, MSTW and NNPDF2.3, together with their uncertainty envelopes.

Fig. 7
figure 7

Predicted NNLO+NNLL \(t\overline{t}\) production cross-sections at \({\sqrt{s}\,{=}\,7~\mathrm{TeV}}\) and \({\sqrt{s}=8~\mathrm{TeV}}\) as a function of \({m}_{t}^\mathrm{pole}\), showing the central values (solid lines) and total uncertainties (dashed lines) with several PDF sets. The yellow band shows the QCD scale uncertainty. The measurements of \({\sigma _{t\overline{t}}}\) are also shown, with their dependence on the assumed value of \({m}_{t}\) through acceptance and background corrections parameterised using Eq. (2)

Figure 7 also shows the small dependence of the experimental measurement of \({\sigma _{t\overline{t}}}\) on the assumed value of \({m}_{t}\), arising from variations in the acceptance and \(Wt\) single top background, as discussed in Sect. 6. This dependence was also parameterised using Eq. (2), giving a derivative of \({\mathrm{d}{\sigma _{t\overline{t}}}/\mathrm{d}{{m}_{t}}}=-0.28\pm 0.03\) %/GeV at 172.5 GeV for both centre-of-mass energies, where the uncertainty is due to the limited size of the simulated samples. Here, \({m}_{t}\) represents the top quark mass used in the Monte Carlo generators, corresponding to that measured in direct reconstruction, rather than the pole mass. However, since this experimental dependence is small, differences between the two masses of up to 2 GeV have a negligible effect (\({<}0.2\) GeV) on the pole mass determination. A comparison of the theoretical and experimental curves shown in Fig. 7 therefore allows an unambiguous extraction of the top quark pole mass.

The extraction is performed by maximising the following Bayesian likelihood as a function of the top quark pole mass \({m}_{t}^\mathrm{pole}\):

$$\begin{aligned} {\mathcal L}({{m}_{t}^\mathrm{pole}})&= \int G({\sigma ^{\prime }_{t\overline{t}}}|{\sigma _{t\overline{t}}}({{m}_{t}^\mathrm{pole}}),{\rho _\mathrm{exp}}) \nonumber \\&\quad \cdot G({\sigma ^{\prime }_{t\overline{t}}}|{\sigma ^\mathrm{theo}_{t\overline{t}}}({{m}_{t}^\mathrm{pole}}),\rho ^{\pm }_\mathrm{theo}~) \mathrm{d}{\sigma ^{\prime }_{t\overline{t}}}. \end{aligned}$$
(3)

Here, \(G(x|\mu ,\rho )\) represents a Gaussian probability density in the variable \(x\) with mean \(\mu \) and standard deviation \(\rho \). The first Gaussian term represents the experimental measurement \({\sigma _{t\overline{t}}}\) with its dependence on \({m}_{t}^\mathrm{pole}\) and uncertainty \(\rho _\mathrm{exp}\). The second Gaussian term represents the theoretical prediction given by Eq. (2) with its asymmetric uncertainty \(\rho ^{\pm }_\mathrm{theo}\) obtained from the quadrature sum of PDF+\(\alpha _\mathrm{s}\) and QCD scale uncertainties evaluated as discussed in Sect. 2. The likelihood in Eq. (3) was maximised separately for each PDF set and centre-of-mass energy to give the \({m}_{t}^\mathrm{pole}\) values shown in Table 6. A breakdown of the contributions to the total uncertainties is given for the CT10 PDF results in Table 7; it can be seen that the theoretical contributions are larger than those from the experimental measurement of \({\sigma _{t\overline{t}}}\). A single \({m}_{t}^\mathrm{pole}\) value was derived for each centre-of-mass energy by defining an asymmetric Gaussian theoretical probability density in Eq. (3) with mean equal to the CT10 prediction, and a \(\pm 1\) standard deviation uncertainty envelope which encompasses the \(\pm 1\) standard deviation uncertainties from each PDF set following the PDF4LHC prescription [8], giving:

$$\begin{aligned} {{m}_{t}^\mathrm{pole}}&= 171.4\pm 2.6~{\hbox {GeV}}\ ({\sqrt{s}\,{=}\,7~\mathrm{TeV}})\quad \mathrm{and} \\ {{m}_{t}^\mathrm{pole}}&= 174.1\pm 2.6~{\hbox {GeV}}\ ({\sqrt{s}=8~\mathrm{TeV}}). \end{aligned}$$

Considering only uncorrelated experimental uncertainties, the two values are consistent at the level of 1.7 standard deviations. The top pole mass was also extracted using a frequentist approach, evaluating the likelihood for each \({m}_{t}^\mathrm{pole}\) value as the Gaussian compatibility between the theoretically predicted and experimentally measured values, and fixing the theory uncertainties to those at \({{m}_{t}^\mathrm{pole}}=172.5\) GeV. The results differ from those of the Bayesian approach by at most 0.2 GeV.

Table 6 Measurements of the top quark pole mass determined from the \(t\overline{t}\)  cross-section measurements at \({\sqrt{s}\,{=}\,7~\mathrm{TeV}}\) and \({\sqrt{s}=8~\mathrm{TeV}}\) using various PDF sets
Table 7 Summary of experimental and theoretical uncertainty contributions to the top quark pole mass determination at \({\sqrt{s}\,{=}\,7~\mathrm{TeV}}\) and \({\sqrt{s}=8~\mathrm{TeV}}\) with the CT10 PDF set

Finally, \({m}_{t}^\mathrm{pole}\) was extracted from the combined \(\sqrt{s}=7~\) TeV and \(\sqrt{s}=8~\) TeV dataset using the product of likelihoods (Eq. (3)) for each centre-of-mass energy and accounting for correlations via nuisance parameters. The same set of experimental uncertainties was considered correlated as for the cross-section ratio measurement, and the uncertainty on \(\sigma ^\mathrm{theo}_{t\overline{t}}\) was considered fully correlated between the two datasets. The resulting value using the envelope of all three considered PDF sets is

$$\begin{aligned} {{m}_{t}^\mathrm{pole}}= 172.9^{+2.5}_{-2.6} ~{\hbox {GeV}}\end{aligned}$$

and has only a slightly smaller uncertainty than the individual results at each centre-of-mass energy, due to the large correlations, particularly for the theoretical predictions. The results are shown in Fig. 8, together with previous determinations using similar techniques from D0 [71, 72] and CMS [73]. All extracted values are consistent with the average of measurements from kinematic reconstruction of \(t\overline{t}\) events of \(173.34\pm 0.76\) GeV [74], showing good compatibility of top quark masses extracted using very different techniques and assumptions.

Fig. 8
figure 8

Comparison of top quark pole mass values determined from this and previous cross-section measurements [7173]. The average of top mass measurements from direct reconstruction [74] is also shown

7.3 Constraints on stop-pair production

Supersymmetry (SUSY) theories predict new bosonic partners for the Standard Model fermions and fermionic partners for the bosons. In the framework of a generic \(R\)-parity conserving minimal supersymmetric extension of the SM [7579], SUSY particles are produced in pairs and the lightest supersymmetric particle is stable. If SUSY is realised in nature and responsible for the solution to the hierarchy problem, naturalness arguments suggest that the supersymmetric partners of the top quark—the top squarks or stops—should have mass close to \({m}_{t}\) in order to effectively cancel the top quark loop contributions to the Higgs mass [80, 81]. In this scenario, the lighter top squark mass eigenstate \({\tilde{t}_{1}}\) would be produced in pairs, and could decay via \({\tilde{t}_{1}}\rightarrow t{\tilde{\chi }_{1}^{0}}\) if \({m_{\tilde{t}_1}}>{{m}_{t}}+{m_{\tilde{\chi }_1^0}}\), where \(\tilde{\chi }_{1}^{0}\), the lightest neutralino, is the lightest supersymmetric particle and is therefore stable. Stop-pair production could therefore give rise to \(t\overline{t}{\tilde{\chi }_{1}^{0}}{\tilde{\chi }_{1}^{0}}\) intermediate states, appearing like \(t\overline{t}\) production with additional missing transverse momentum carried away by the escaping neutralinos. The predicted cross-sections at \(\sqrt{s}=8~\) TeV are about 40 pb for \({m_{\tilde{t}_1}}=175\) GeV, falling to 20 pb for 200 GeV. If the top squark mass \(m_{\tilde{t}_1}\) is smaller than about 200 GeV, such events would look very similar to SM QCD \(t\overline{t}\) production, making traditional searches exploiting kinematic differences very difficult, but producing a small excess in the measured \(t\overline{t}\) cross-section, as discussed e.g. in Refs. [82, 83].

Fig. 9
figure 9

Expected and observed limits at 95 % CL on the signal strength \(\mu \) as a function of \(m_{\tilde{t}_1}\), for pair produced top squarks \({\tilde{t}_{1}}\) decaying with 100 % branching ratio via \({\tilde{t}_{1}}\rightarrow t{\tilde{\chi }_{1}^{0}}\) to predominantly right-handed top quarks, assuming \({m_{\tilde{\chi }_1^0}}=1\) GeV. The black dotted line shows the expected limit with \(\pm 1\sigma \) contours, taking into account all uncertainties except the theoretical cross-section uncertainties on the signal. The red solid line shows the observed limit, with dotted lines indicating the changes as the nominal signal cross-section is scaled up and down by its theoretical uncertainty

The potential stop-pair signal yield was studied for top squark masses in the range 175–225 GeV and neutralino masses in the range 1 GeV\(<{m_{\tilde{\chi }_1^0}}<{m_{\tilde{t}_1}}-{{m}_{t}}\) using simulated samples generated with Herwig++ [84] with the CTEQ6L1 PDFs [32], and NLO+NLL production cross-sections [8587]. The mixing matrices for the top squarks and the neutralinos were chosen such that the top quark produced in the \({\tilde{t}_{1}}\rightarrow t{\tilde{\chi }_{1}^{0}}\) decay has a right-handed polarisation in 95 % of the decays. Due to the slightly more central \(|\eta |\) distribution of the leptons from the subsequent \(t\rightarrow Wq\), \(W\rightarrow \ell \nu \) decay, the preselection efficiency \({\epsilon _{e\mu }}\) for these events is typically 10–20 % higher than for SM QCD \(t\overline{t}\), increasing with \(m_{\tilde{t}_1}\). However, the fraction of preselected events with one or two \(b\)-tagged jets is very similar to the SM case. The effect of a small admixture of stop-pair production in addition to the SM \(t\overline{t}\) production is therefore to increase the measured cross-section by \({R_{\tilde{t}_1\tilde{t}_1}}{\sigma _{\tilde{t}_1\tilde{t}_1}}\), where \(R_{\tilde{t}_1\tilde{t}_1}\)is the ratio of \({\epsilon _{e\mu }}\) values for stop-pair and SM \(t\overline{t}\) production, and \(\sigma _{\tilde{t}_1\tilde{t}_1}\) is the stop-pair production cross-section.

Limits were set on stop-pair production by fitting the effective production cross-section \({R_{\tilde{t}_1\tilde{t}_1}}{\sigma _{\tilde{t}_1\tilde{t}_1}}\) multiplied by a signal strength \(\mu \) to the difference between the measured cross-sections (\({\sigma _{t\overline{t}}}\)) and the theoretically predicted SM QCD production cross-sections (\(\sigma ^\mathrm{theo}_{t\overline{t}}\)). The two datasets were fitted simultaneously, assuming values of \({\sigma ^\mathrm{theo}_{t\overline{t}}}=177.3^{+11.5}_{-12.0}\) pb for \(\sqrt{s}=7~\) TeV and \(252.9^{+15.3}_{-16.3}\) pb for \(\sqrt{s}=8~\) TeV, including the uncertainty due to a \(\pm 1\) GeV variation in the top quark mass. The limits were determined using a profile likelihood ratio in the asymptotic limit [88], using nuisance parameters to account for correlated theoretical and experimental uncertainties.

The observed and expected limits on \(\mu \) at the 95 % confidence level (CL) were extracted using the CLs prescription [89] and are shown in Fig. 9. Due to the rapidly decreasing stop-pair production cross-section with increasing \(m_{\tilde{t}_1}\), the analysis is most sensitive below 180 GeV. Adopting the convention of reducing the estimated SUSY production cross-section by one standard deviation of its theoretical uncertainty (15 %, coming from PDFs and QCD scale uncertainties [90]), stop masses between the top mass threshold and 177 GeV are excluded, assuming 100 % branching ratio for \({\tilde{t}_{1}}\rightarrow t{\tilde{\chi }_{1}^{0}}\) and \({m_{\tilde{\chi }_1^0}}=1\) GeV. The limits from considering the \({\sqrt{s}\,{=}\,7~\mathrm{TeV}}\) and \({\sqrt{s}=8~\mathrm{TeV}}\) datasets separately are only slightly weaker, due to the large correlations in the systematic uncertainties between beam energies, particularly for the theoretical predictions. At each energy, they correspond to excluded stop-pair production cross-sections of 25–27 pb at 95 % CL.

The combined cross-section limits depend only slightly on the neutralino mass, becoming e.g. about 3 % weaker at \({m_{\tilde{t}_1}}=200\) GeV for \({m_{\tilde{\chi }_1^0}}=20\) GeV. However, they depend more strongly on the assumed top quark polarisation; in a scenario with \({m_{\tilde{t}_1}}=175\) GeV and \({m_{\tilde{\chi }_1^0}}=1\) GeV, and squark mixing matrices chosen such that the top quarks are produced with full left-handed polarisation, the limits are 4 % weaker than the predominantly right-handed case, rising to 14 % weaker for \({m_{\tilde{t}_1}}=200\) GeV. In this scenario, top squarks with masses from \({m}_{t}\) to 175 GeV can be excluded. Although the analysis has some sensitivity to three-body top squark decays of the form \({\tilde{t}_{1}}\rightarrow bW{\tilde{\chi }_{1}^{0}}\) for \({m_{\tilde{t}_1}}<{{m}_{t}}\), the \(b\)-jets become softer than those from SM \(t\overline{t}\) production, affecting the determination of \({\epsilon _{b}}\). Therefore, no limits are set for scenarios with \({m_{\tilde{t}_1}}<{{m}_{t}}\).

8 Conclusions

The inclusive \(t\overline{t}\) production cross-section has been measured at the LHC using the full ATLAS 2011–2012 \(pp\) collision data sample of 4.6 \({\mathrm{fb}^{-1}}\) at \(\sqrt{s}=7~\) TeV and 20.3 \({\mathrm{fb}^{-1}}\) at \(\sqrt{s}=8~\) TeV, in the dilepton \(t\overline{t}\rightarrow e\mu \nu {\overline{\nu }}{b\overline{b}}\) decay channel. The numbers of opposite-sign \(e\mu \) events with one and two \(b\)-tagged jets were counted, allowing a simultaneous determination of the \(t\overline{t}\) cross-section \({\sigma _{t\overline{t}}}\) and the probability to reconstruct and \(b\)-tag a jet from a \(t\overline{t}\) decay. Assuming a top quark mass of \({m}_{t}\) \(=172.5\) GeV, the results are:

$$\begin{aligned} {\sigma _{t\overline{t}}}&= 182.9\,{\pm }\,3.1\,{\pm }\,4.2\,{\pm }\,3.6\,{\pm }\,3.3~\mathrm{pb}~ ({\sqrt{s}\,{=}\,7~\mathrm{TeV}})\quad \mathrm{and} \\ {\sigma _{t\overline{t}}}&= 242.4\,{\pm }\,1.7\,{\pm }\,5.5\,{\pm }\,7.5\,{\pm }\,4.2~\mathrm{pb}~ ({\sqrt{s}=8~\mathrm{TeV}}), \end{aligned}$$

where the four uncertainties arise from data statistics, experimental and theoretical systematic effects, knowledge of the integrated luminosity, and of the LHC beam energy, giving total uncertainties of 7.1 pb (3.9 %) and 10.3 pb (4.3 %) at \(\sqrt{s}=7~\) TeV and \(\sqrt{s}=8~\) TeV. The dependence of the results on the assumed value of \({m}_{t}\) is \({\mathrm{d}{\sigma _{t\overline{t}}}/\mathrm{d}{{m}_{t}}}=-0.28\) %/GeV, and the associated uncertainty is not included in the totals given above. The results are consistent with recent NNLO+NNLL QCD calculations, and have slightly smaller uncertainties than the theoretical predictions. The ratio of the two cross-sections, and measurements in fiducial ranges corresponding to the experimental acceptance, have also been reported.

The measured \(t\overline{t}\) cross-sections have been used to determine the top quark pole mass via the dependence of the predicted cross-section on \({m}_{t}^\mathrm{pole}\), giving a value of \({{m}_{t}^\mathrm{pole}}=172.9^{+2.5}_{-2.6}\) GeV, compatible with the mass measured from kinematic reconstruction of \(t\overline{t}\) events.

The results have also been used to search for pair-produced supersymmetric top squarks decaying to top quarks and light neutralinos. Assuming 100 % branching ratio for the decay \({\tilde{t}_{1}}\rightarrow t{\tilde{\chi }_{1}^{0}}\), and the production of predominantly right-handed top quarks, top squark masses between the top quark mass and 177 GeV are excluded at 95 % CL.