Measurement of the $t\bar{t}$ production cross-section using $e\mu$ events with $b$-tagged jets in $pp$ collisions at $\sqrt{s}=7$ and 8 TeV with the ATLAS detector

The inclusive top quark pair ($t\bar{t}$) production cross-section $\sigma_{t\bar{t}}$ has been measured in $pp$ collisions at $\sqrt{s}=7$ TeV and $\sqrt{s}=8$ TeV with the ATLAS experiment at the LHC, using $t\bar{t}$ events with an opposite-charge $e\mu$ pair in the final state. The measurement was performed with the 2011 7 TeV dataset corresponding to an integrated luminosity of 4.6 fb$^{-1}$ and the 2012 8 TeV dataset of 20.3 fb$^{-1}$. The cross-section was measured to be: $\sigma_{t\bar{t}}=182.9\pm 3.1\pm 4.2\pm 3.6 \pm 3.3$ pb ($\sqrt{s}=7$ TeV) and $\sigma_{t\bar{t}}=242.9\pm 1.7\pm 5.5\pm 5.1\pm 4.2$ pb ($\sqrt{s}=8$ TeV, updated as described in the Addendum), where the four uncertainties arise from data statistics, experimental and theoretical systematic effects, the knowledge of the integrated luminosity and of the LHC beam energy. The results are consistent with recent theoretical QCD calculations at next-to-next-to-leading order. Fiducial measurements corresponding to the experimental acceptance of the leptons are also reported, together with the ratio of cross-sections measured at the two centre-of-mass energies. The inclusive cross-section results were used to determine the top quark pole mass via the dependence of the theoretically-predicted cross-section on $m_t^{\rm pole}$, giving a result of $m_t^{\rm pole}=172.9^{+2.5}_{-2.6}$ GeV. By looking for an excess of $t\bar{t}$ production with respect to the QCD prediction, the results were also used to place limits on the pair-production of supersymmetric top squarks $\tilde{t}_1$ with masses close to the top quark mass decaying via $\tilde{t}_1\rightarrow t\tilde{\chi}^0_1$ to predominantly right-handed top quarks and a light neutralino $\tilde{\chi}_0^1$, the lightest supersymmetric particle. Top squarks with masses between the top quark mass and 177 GeV are excluded at the 95% confidence level.

where the four uncertainties arise from data statistics, experimental and theoretical systematic effects, knowledge of the integrated luminosity and of the LHC beam energy.The results are consistent with recent theoretical QCD calculations at next-to-next-to-leading order.Fiducial measurements corresponding to the experimental acceptance of the leptons are also reported, together with the ratio of cross-sections measured at the two centre-of-mass energies.The inclusive cross-section results were used to determine the top quark pole mass via the dependence of the theoretically predicted crosssection on m pole t giving a result of m pole t = 172.9+2.5 −2.6 GeV.By looking for an excess of t t production with respect to the QCD prediction, the results were also used to place limits on the pair-production of supersymmetric top squarks t1 with masses close to the top quark mass, decaying via t1 → t χ0

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 t), and the precise prediction of the corresponding inclusive cross-section (σ t t) is a substantial challenge for quantum chromodynamics (QCD) calculational techniques.Precise measurements of σ t 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 t production are governed by the decay modes of the two W bosons.This paper describes a measurement in the dileptonic eµ channel, t t → W + bW −b → e ± µ ∓ ννb b, selecting events with an eµ pair with opposite-sign electric charges, 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µ pair and one or two tagged b-jets were used to measure simultaneously the t 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 τ decays t → W b → τ νb → e/µνννb, were included as part of the t t signal.
The main background is W t, the associated production of a W boson and a single top quark.Other background contributions arise from Z → τ τ → eµ+jets (+4ν) 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 σ t 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 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.

Theoretical cross-section predictions
Calculations of σ t t for hadron collisions are now available at full next-to-next-to-leading-order (NNLO) accuracy in the strong coupling constant α s , including the resummation of next-to-next-to-leading logarithmic (NNLL) soft gluon terms [1,2].At a centre-of-mass energy of √ s = 7 TeV and assuming m t = 172.5 GeV, these calculations give a prediction of 177.3 ± 9.0 +4.6 −6.0 pb, where the first uncertainty is due to PDF and α s uncertainties, and the second to QCD scale uncertainties.The corresponding prediction at √ s = 8 TeV is 252.9 ± 11.7 +6.4 −8.6 pb.These values were calculated using the top++ 2.0 program [3].The PDF and α s uncertainties were calculated using the PDF4LHC prescription [4] with the MSTW2008 68 % CL NNLO [5], CT10 NNLO [6,7] and NNPDF2.35f FFN [8] 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 crosssections at √ s = 8 TeV and √ s = 7 TeV is predicted to be 1.430 ± 0.013 (PDF+α s ) ±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 [9].For comparison, the corresponding nextto-leading-order (NLO) predictions, also calculated using top++ 2.0 with the same set of PDFs, are 157±12±24 pb at √ s = 7 TeV and 225 ± 16 ± 29 pb at √ s = 8 TeV, where again the first quoted uncertainties are due to PDF and α 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.

Data and simulated samples
The ATLAS detector [10] 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 pseudorapidity2 range |η| < 2.5, complemented by a transition radiation tracker providing tracking and electron identification information for |η| < 2.0.A lead/liquidargon (LAr) electromagnetic calorimeter covers the region |η| < 3.2, and hadronic calorimetry is provided by steel/scintillator tile calorimeters for |η| < 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 |η| < 2.7, and separate trigger chambers covering |η| < 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 fb −1 at √ s = 7 TeV and 20.3 fb −1 at √ 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 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 ( √ s = 7 TeV) or 20 ( √ 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 [11] based on GEANT4 [12], or through a faster simulation making use of parameterised showers in the calorimeters [13].Additional simulated pp collisions generated either with Pythia6 [14] (for √ s = 7 TeV simulation) or Pythia8 [15] (for √ s = 8 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 t full simulation sample was produced using the NLO matrix element generator Powheg [16] interfaced to Pythia6 [14] with the Perugia 2011C tune (P2011C) [17] for parton shower, fragmentation and underlying event modelling, and CT10 PDFs [6], and included all t t final states involving at least one lepton.The W → ν branching ratio was set to the SM expectation of 0.1082 [18], and m t was set to 172.5 GeV.Alternative t t samples were produced with the NLO generator MC@NLO [19] interfaced to Herwig [20] with Jimmy [21] for the underlying event modelling, with the ATLAS AUET2 [22] tune and CT10 PDFs; and with the leading-order (LO) multileg generator Alpgen [23] interfaced to either Pythia6 or Herwig and Jimmy, with the CTEQ6L1 PDFs [24].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 τ 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 W t single top production, modelled using Powheg + Pythia6 [25] with the CT10 PDFs and the P2011C tune; Z → τ τ +jets modelled using Alpgen + Herwig + Jimmy ( √ s = 7 TeV) or Alpgen + Pythia6 including LO matrix elements for Zb b production, with CTEQ6L1 PDFs; and diboson (W W , W Z, ZZ) production in association with jets, modelled using Alpgen + Herwig + Jimmy.The W t background was normalised to approximate NNLO cross-sections of 15.7 ± 1.2 pb at √ s = 7 TeV and 22.4 ± 1.5 pb at √ s = 8 TeV, determined as in Ref. [26].The inclusive Z cross-sections were set to the NNLO predictions from FEWZ [27], but the normalisation of Z → τ τ → eµ4ν 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 [28].Production of t t in association with a W or Z boson, which contributes to the sample with same-sign leptons, was simulated with Madgraph [29] interfaced to Pythia with CTEQ6L1 PDFs, and normalised to NLO cross-section predictions [30].
Backgrounds with one real and one misidentified lepton include t t events with one hadronically decaying W ; W +jets production, modelled as for Z+jets; W γ+jets, modelled with Sherpa [31] with CT10 PDFs; and t-channel single top production, modelled using AcerMC [32] interfaced to Pythia6 with CTEQ6L1 PDFs.Other backgrounds, including processes with two misidentified leptons, are negligible after the event selections used in this analysis.

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 [33], with transverse energy E T > 25 GeV and pseudorapidity |η| < 2.47.Electron candidates within the transition region between the barrel and endcap electromagnetic calorimeters, 1.37 < |η| < 1.52, were removed.Isolation requirements were used to reduce background from non-prompt electrons.The calorimeter transverse energy within a cone of size ∆R = 0.2 and the scalar sum of track p T within a cone of size ∆R = 0.3, in each case excluding the contribution from the electron itself, were each required to be smaller than E T and η-dependent thresholds calibrated to separately give nominal selection efficiencies of 98 % for prompt electrons from Z → ee decays.
Muon candidates were reconstructed by combining matching tracks reconstructed in both the inner detector and muon spectrometer [34], and were required to satisfy p T > 25 GeV and |η| < 2.5.In the √ s = 7 TeV dataset, the calorimeter transverse energy within a cone of size ∆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 T within a cone of size ∆R = 0.3, excluding the muon track, was required to be less than 2.5 GeV.In the √ 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 T in a variable-sized cone of radius ∆R = 10 GeV/p µ T to the transverse momentum p µ T of the muon.This isolation criterion performs better for highly boosted top quark decays [35].Both sets of isolation requirements have efficiencies of about 97 % for prompt muons from Z → µµ decays.
Jets were reconstructed using the anti-k t algorithm [36] with radius parameter R = 0.4, starting from calorimeter energy clusters calibrated at the electromagnetic energy scale for the √ s = 7 TeV dataset, or using the local cluster weighting method for √ s = 8 TeV [37].Jets were calibrated using an energy-and η-dependent simulationbased calibration scheme, with in-situ corrections based on data, and were required to satisfy p T > 25 GeV and |η| < 2.5.To suppress the contribution from low-p T jets originating from pileup interactions, a jet vertex fraction requirement was applied: at √ s = 7 TeV jets were required to have at least 75 % of the scalar sum of the p T of tracks associated with the jet come from tracks associated with the event primary vertex.The latter was defined as the reconstructed vertex with the highest sum of associated track p 2 T .Motivated by the higher pileup background, in the √ s = 8 TeV dataset this requirement was loosened to 50 %, only applied to jets with p T < 50 GeV and |η| < 2.4, and the effects of pileup on the jet energy calibration were further reduced using the jet-area method as described in Ref. [38].Finally, to further suppress non-isolated leptons likely to have come from heavy-flavour decays inside jets, electrons and muons within ∆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 [39,40].Jets were defined to be b-tagged if the MV1 discriminant value was larger than a threshold corresponding approximately to a 70 % efficiency for tagging b-quark jets from top decays in t t events, with a rejection factor of about 140 against lightquark 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 oppositesign eµ pair constituted the main analysis sample, whilst events with a same-sign eµ pair were used in the estimation of the background from misidentified leptons.

Extraction of the t t cross-section
The t t production cross-section σ t t was determined by counting the numbers of opposite-sign eµ 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: where L is the integrated luminosity of the sample, eµ is the efficiency for a t t event to pass the opposite-sign eµ 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 → W q 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 b .Although this quark is almost always a b quark, b thus also accounts for the approximately 0.2 % of top quarks that decay to W s or W d rather than W b, slightly reducing the effective b-tagging efficiency.Furthermore, the value of 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 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.
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 The data are shown compared to the expectation from simulation, broken down into contributions from t t, W t 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 t signal samples generated with Powheg + Pythia6 (PY), MC@NLO + Herwig (HW) and Alpgen + Herwig, and with the cyan band indicating the statistical uncertainty.
bb is given by bb = 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 2 , where N t t eµ is the number of preselected eµ t t events and N t t 1 and N t t 2 are the numbers of t 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 b , N 1 and N 2 of the small number of mistagged charm-quark or gluon jets from radiation in the t t events.
Background from sources other than t t → eµννb b also contributes to the event counts N 1 and N 2 , and is given by the terms N bkg 1 and N bkg 2 .The preselection efficiency eµ and tagging correlation C b were taken from t t event simulation, and the background contributions N bkg 1 and N bkg 2 were estimated using a combination of simulationand data-based methods, allowing the two equations in Eq. ( 1) to be solved numerically yielding σ t t and b .
A total of 11796 events passed the eµ opposite-sign preselection in √ s = 7 TeV data, and 66453 in √ s = 8 TeV data.Table 1 shows the number of events with one and two b-tagged jets, together with the estimates of non-t 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 t events, with the dominant background coming from W t 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 t events, with W t production again being the dominant background.
Distributions of the number of b-tagged jets in oppositesign eµ events are shown in Fig. 1, and compared to the expectations with several t t simulation samples.Distributions of the number of jets, the b-tagged jet p T , and the electron and muon |η| and p T are shown for oppositesign eµ events with at least one b-tagged jet in Fig. 2 ( √ s = 7 TeV) and Fig. 3 ( √ s = 8 TeV), with the simulation normalised to the same number of entries as the data.In general, the agreement between data and simulation is good, within the range of predictions from the different t t simulation samples.
The value of σ t t extracted from Eq. ( 1) is inversely proportional to the assumed value of eµ , with (dσ t t/d eµ ) / (σ t t/ eµ ) = −1.Uncertainties on eµ therefore translate directly into uncertainties on σ t t.The value of eµ was determined from simulation to be about 0.8 % for both centre-of-mass energies, and includes the t t → eµννb b branching ratio of about 3.2 % including W → τ → e/µ decays.Similarly, σ t t is proportional to the value of C b , also determined from simulation, giving a dependence with the opposite sign, (dσ t t /dC b )/(σ t t/C b ) = 1.The systematic uncertainties on eµ and C b are discussed in Sect.6.
With the kinematic cuts and b-tagging working point chosen for this analysis, the sensitivities of σ t t to knowledge of the backgrounds N bkg 1 and N bkg 2 are given by (dσ t t/dN bkg 1 ) / (σ t t/N bkg 1 ) = −0.12 and (dσ t t/dN bkg 2 ) / (σ t t/N bkg 2 ) = −0.004.The fitted cross-sections are therefore most sensitive to the systematic uncertainties on N bkg 1 , whilst for the chosen b-tagging working point, the measurements of N 2 serve mainly to constrain 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.

Background estimation
The W t single top and diboson backgrounds were estimated from simulation as discussed in Sect.3. The Z+jets background (with Z → τ τ → eµ4ν) at √ s = 8 TeV was estimated from simulation using Alpgen+Pythia, scaled by the ratios of Z → ee or Z → µµ accompanied by btagged 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-tosimulation 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 t production and misidentified leptons.Combining the results from both dilepton channels, the scale factors were determined to be 1.43 ± 0.03 and 1.13 ± 0.08 for the one and two b-tag backgrounds, after normalising the simulation to the inclusive Z cross-section prediction from FEWZ [27].The uncertainties include systematic components derived from a comparison of results from the ee and µµ channels, and from studying the variation of scale factors with Z boson p T .The average p T is higher in selected Z → τ τ → eµ4ν events than in Z → ee/µµ events due to the momentum lost to the undetected neutrinos from the τ decays.The same procedure was used for the √ s = 7 TeV dataset, resulting in scale factors of 1.23 ± 0.07 (one btag) and 1.14 ± 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µ 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 samesign contributions, multiplied by the opposite-to samesign misidentified-lepton ratios R j = N mis,OS j /N mis,SS j estimated from simulation for events with j = 1 and 2 btagged 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 TeV data are compared to the expectation from simulation, broken down into contributions from t 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 t signal samples and with the cyan band indicating the statistical uncertainty.The last bin includes the overflow.TeV data are compared to the expectation from simulation, broken down into contributions from t 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 t signal samples and with the cyan band indicating the statistical uncertainty.The last bin includes the overflow.Table 1.Observed numbers of opposite-sign eµ events with one and two b-tagged jets (N1 and N2) for each data sample, together with the estimates of backgrounds and associated total uncertainties described in Sect.6 from prompt electrons produced from t → W q → eνq in signal t t → eµννb 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 ± 0.5 and R 2 = 1.2 ± 0.5 at √ s = 7 TeV, and R 1 = 1.2 ± 0.3 and R 2 = 1.6 ± 0.5 at √ s = 8 TeV.The uncertainties were derived by considering the range of R j values for different components of the misidentifiedlepton 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 t events with an additional leptonically decaying W or Z boson, diboson decays producing two samesign leptons, and wrong-sign t t → eµννb 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 [33] 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 |η| and p T distributions of the leptons in √ s = 8 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.

Systematic uncertainties
The systematic uncertainties on the measured cross-sections σ t t are shown in detail in Table 3 together with the individual uncertainties on eµ and C b .A summary of the uncertainties on σ t 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 ±1 standard deviation.Systematic correlations between 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 pT distributions, the last bin includes the overflows.
input parameters (in particular significant anti-correlations between eµ and C b which contribute with opposite signs to σ t t) were thus taken into account.The total uncertainties on σ t t and 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.
t t modelling: Uncertainties on eµ and C b due to the simulation of t t events were assessed by studying the predictions of different t t generators and hadronisation models as detailed in Sect.3. The prediction for eµ 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 dis-cussed below.The remaining uncertainties on eµ relating to lepton reconstruction, identification and leptonjet overlap removal, were evaluated from the differences between the predictions from the baseline Powheg + Pythia t 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 eµ 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 eµ up to 2 % higher than that of Powheg+Pythia, due largely to a more central predicted η distribution for the leptons.How- ever, this sample uses a leading-order generator and PDFs, and gives an inferior description of the electron and muon η distributions (see Figs. 3(c) and 3(e)), so was not used to set the systematic uncertainty on eµ .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 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 t events at √ s = 7 TeV [41], generated using both Ac-erMC + 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 eµ , C b and the W t single top background due to uncertainties on the proton PDFs were evaluated using the error sets of the CT10 NLO [6], MSTW 2008 68 % CL NLO [5] and NNPDF 2.3 NLO [8] 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 [4].QCD scale choices: The lepton p T and η distributions, and hence eµ , are sensitive to the choices of QCD renormalisation and factorisation scales.This effect was investigated using √ s = 8 TeV generator-level Powheg + Pythia t 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 2 T,t .The systematic uncertainty for each scale was taken as half the difference in eµ 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 eµ , assumed to be valid for both centre-of-mass energies.Single top modelling: Uncertainties related to W t 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 W t with an additional bjet, a NLO contribution to W t which can interfere with the t t final state.The sensitivity to this interference was studied by comparing the predictions of Powheg with the diagram-removal (baseline) and diagramsubtraction schemes [25,42], giving additional singletop/t 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 W t single top cross-section were taken to be 7.6 % and 6.8% at √ s = 7 TeV and √ s = 8 TeV, based on Ref. [26].The uncertainties on the diboson cross-sections were set to 5 % [28].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 [31] 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 W H and ZH associated production modes, with H → W W and H → τ τ , 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 → ee/µµ to Z → τ τ events result from statistical uncertainties, comparing the results from ee and µµ, which have different background compositions, and considering the dependence of the scale factors on Z boson p 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 → ee/µµ, J/ψ → ee/µµ and W → eν events in data and simulation, using the techniques described in Refs.[33,34,43].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µ 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µ samples as described in Sect.5, or by using the distributions of lepton impact parameter significance |d 0 |/σ d0 , where d 0 is the distance of closest approach of the lepton track to the event primary vertex in the transverse plane, and σ d0 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 eµ , with uncertainties dominated by the limited sizes of the samesign and high impact-parameter significance samples used for background estimation.Similar results were found from studies in Z → ee and Z → µµ events, after correcting the results for the larger average amount of hadronic activity near the leptons in t t → eµννb b events.Jet-related uncertainties: Although the efficiency to reconstruct and b-tag jets from t 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 eµ via the lepton-jet ∆R separation cuts.The jet energy scale was varied in simulation according to the uncertainties derived from simulation and in-situ calibration measurements [37,44], using a model with 21 ( √ s = 7 TeV) or 22 ( √ s = 8 TeV) separate orthogonal uncertainty components which were then added in quadrature.The jet energy resolution was found to be well modelled by simulation [45], and remaining uncertainties were assessed by applying additional smearing, which reduces eµ .The calorimeter jet reconstruction efficiency was measured in data us-ing 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 → ee/µµ+jets events.b-tagging uncertainties: The efficiency for b-tagging jets from t 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 W t single top and diboson backgrounds.The tagging correlation C b is also slightly sensitive to the efficiencies for tagging heavy-and lightflavour jets.The uncertainties in the simulation modelling of the b-tagging performance were assessed using studies of b-jets containing muons [40,46], jets containing D * + mesons [47] and inclusive jet events [48].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 √ s = 7 TeV dataset is 1.8 % [49].Using beam-separation scans performed in November 2012, the same methodology was applied to determine the √ s = 8 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-ofmass energies.The relative uncertainties on the crosssection measurements are slightly larger than those on the luminosity measurements because the W t 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 ± 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 [50].Since this calibration is compatible with the nominal √ s of 8 TeV, no correction was applied to the measured σ t t value.However, an uncertainty of 1.72 %, corresponding to the expected change in σ t t for a 0.66 % change in √ s is quoted separately on the final result.Following Ref. [50], the same relative uncertainty on the LHC beam energy is applied for the √ s = 7 TeV dataset, giving a slightly larger uncertainty of 1.79 % due to the steeper relative dependence of σ t t on √ s in this region.These uncertainties are much larger than those corresponding to the very small dependence of eµ on √ 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 t and W t events, and the W t 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 dσ t t/dm t = −0.28%/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.As shown in Tables 3 and 4, the largest systematic uncertainties on σ t t come from t t modelling and PDFs, and knowledge of the integrated luminosities and LHC beam energy.

Additional correlation studies
The tagging correlation C b was determined from simulation to be 1.009 ± 0.002 ± 0.007 ( √ s = 7 TeV) and 1.007 ± 0.002 ± 0.006 ( √ s = 8 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 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 b or cc pairs in t 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 σ t 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 ± 0.  prediction of 2.4 ± 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 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 T and η.For example, t 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 T , and longitudinal boosts of the t t system along the beamline give rise to η correlations between the two jets.These effects were probed by increasing the jet p 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 ≈ 1.2 for a 75 GeV p T cut.As shown for the √ 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 T cut range, suggesting that any such kinematic correlations are well modelled by the simulation.The results were also found to be stable within the uncorrelated components of the statistical and systematic uncertainties when tightening the jet and lepton η cuts, raising the lepton p 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.

Results
Combining the estimates of eµ and C b from simulation samples, the estimates of the background N bkg 1 and N bkg 2 shown in Table 1 and the data integrated luminosities, the √ s = 7 TeV using the ee, µµ and eµ channels [51] and using a fit to jet multiplicities and missing transverse momentum in the eµ channel [52].The √ s = 7 TeV measurements are displaced slightly for clarity.The NNLO+NNLL prediction [2,3] described in Sect. 2 is also shown as a function of √ s, for fixed mt = 172.5 GeV and with the uncertainties from PDFs, αs and QCD scale choices indicated by the green band.t t cross-section was determined by solving Eq. ( 1) to be: 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 √ s = 7 TeV and 10.3 pb (4.3 %) at √ 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 dσ t t/dm t = −0.28%/GeV.The product of jet reconstruction and b-tagging efficiencies b was measured to be 0.557 ± 0.009 at √ s = 7 TeV and 0.540 ± 0.006 at √ s = 8 TeV, in both cases consistent with the values in simulation.
The results are shown graphically as a function of √ s in Fig. 6, together with previous ATLAS measurements of σ t t at √ s = 7 TeV in the ee, µµ and eµ dilepton channels using a count of the number of events with two leptons and at least two jets in an 0.7 fb −1 dataset [51], and using a fit of jet multiplicities and missing transverse momentum in the eµ dilepton channel alone with the full 4.6 fb −1 dataset [52].The √ 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 centreof-mass energies are consistent with the NNLO+NNLL QCD calculations discussed in Sect. 2. The √ s = 7 TeV result is 13 % higher than a previous measurement by the CMS collaboration [53], whilst the √ s = 8 TeV result is consistent with that from CMS [54].
From the present analysis, the ratio of cross-sections R t t =σ t t(8 TeV)/σ t t(7 TeV) was determined to be: R t t = 1.326 ± 0.024 ± 0.015 ± 0.049 ± 0.001 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±0.013(see Sect. 2), which in addition to the quoted PDF, α s and QCD scale uncertainties varies by only ±0.001 for a ±1 GeV variation of m t .

Fiducial cross-sections
The preselection efficiency eµ can be written as the product of two terms eµ = A eµ G eµ , where the acceptance A eµ represents the fraction of t t events which have a true opposite-sign eµ pair from t → W → decays (including via W → τ → ), each with p T > 25 GeV and within |η| < 2.5, and G eµ 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 σ fid t t can then be defined as σ fid t t = A eµ σ t t, and measured by replacing σ t t eµ with σ fid t t G eµ in Eq. ( 1), leaving the background terms unchanged.Measurement of the fiducial cross-section avoids the systematic uncertainties associated with A eµ , i.e. the extrapolation from the measured lepton phase space to the full phase space populated by inclusive t 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 b , no restrictions on jet kinematics are imposed in the definition of σ fid t t .In calculating A eµ and G eµ from the various t 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 ∆R = 0.1 around the lepton direction, excluding photons from hadron decays or produced in interactions with detector material.The values of A eµ are about 1.4 % (including the t t → eµννb b branching ratio), and those of G eµ about 55 %, at both centre-of-mass energies.
The measured fiducial cross-sections at √ s = 7 TeV and √ s = 8 TeV, for leptons with p T > 25 GeV and |η| < 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 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±0.1 % of t t events in the fiducial region have at least one lepton produced via W → τ → 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 T > 30 GeV and |η| < 2.4, a common definition agreed with the CMS collaboration.

Top quark mass determination
The strong dependence of the theoretical prediction for σ t t on m t offers the possibility of interpreting measurements of σ t t as measurements of m t .The theoretical calculations use the pole mass m pole t , 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 [55] may differ from the pole mass by O(1 GeV) [56].It is therefore interesting to compare the values of m t determined from the two approaches, as explored previously by the D0 [57] and CMS [58] collaborations.
The dependence of the cross-section predictions (calculated as described in Sect.2) on m pole t is shown in Fig. 7  ), a 1 and a 2 are free parameters.This function was used to parameterise the dependence of σ t t on m t separately for each of the NNLO PDF sets CT10, MSTW and NNPDF2.3,together with their uncertainty envelopes.
Figure 7 also shows the small dependence of the experimental measurement of σ t t on the assumed value of m t , arising from variations in the acceptance and W t single top background, as discussed in Sect.6.This dependence was also parameterised using Eq. ( 2), giving a derivative of dσ t t/dm t = −0.28 ± 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 TeV, for different requirements on the minimum lepton pT and maximum lepton |η|, and with or without the inclusion of leptons from W → τ → 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.
Here, G(x|µ, ρ) represents a Gaussian probability density in the variable x with mean µ and standard deviation ρ.The first Gaussian term represents the experimental measurement σ t t with its dependence on m pole t and uncertainty ρ exp , and the second Gaussian term represents the theoretical prediction given by Eq. ( 2) with its asymmetric uncertainty ρ ± theo from PDFs, α s and QCD scale choices as discussed in Sect. 2. The likelihood in Eq. ( 3) was maximised separately for each PDF set and centreof-mass energy to give the m pole t values shown in Table 6.The theoretical contributions to the total uncertainty are slightly larger than the experimental uncertainties from the measurement of σ t t; the latter are about 1.4 GeV at √ s = 7 TeV and 1.6 GeV at √ s = 8 TeV.A single m pole t value was derived for each centre-of-mass energy by defin- 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 pole t value as the Gaussian compatibility between the theoretically predicted and experimentally measured values, and fixing the theory uncertainties to those at m pole t = 172.5 GeV.The results differ from those of the Bayesian approach by at most 0.2 GeV.
Finally, m pole t was extracted from the combined √ s = 7 TeV and √ 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 σ theo t t was considered fully correlated between the two datasets.The resulting value using the envelope of all three considered PDF sets is m pole t = 172.9+2.5 −2.6 GeV 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 [57] [GeV]  [57,58].The average of top mass measurements from direct reconstruction [59] is also shown.
and CMS [58].All extracted values are consistent with the average of measurements from kinematic reconstruction of t t events of 173.34 ± 0.76 GeV [59], showing good compatibility of top quark masses extracted using very different techniques and assumptions.

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 [60], 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 [61].In this scenario, the lighter top squark mass eigenstate t1 would be produced in pairs, and could decay via t1 → t χ0 1 if m t1 > m t + m χ0 1 , where χ0 1 , the lightest neutralino, is the lightest supersymmetric particle and is therefore stable.Stop-pair production could therefore give rise to t t χ0 1 χ0 1 intermediate states, appearing like t t production with additional missing transverse momentum carried away by the escaping neutralinos.The predicted cross-sections at √ s = 8 TeV are about 40 pb for m t1 = 175 GeV, falling to 20 pb for 200 GeV.If the top squark mass m t1 is smaller than about 200 GeV, such events would look very similar to SM QCD t t production, making traditional searches exploiting kinematic differences very difficult, but producing a small excess in the measured t t cross-section.
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 χ0 1 < m t1 − m t using simulated samples generated with Herwig++ [62] with the CTEQ6L1 PDFs [24], and NLO+NLL production crosssections [63].The mixing matrices for the top squarks and the neutralinos were chosen such that the top quark produced in the t1 → t χ0 1 decay has a right-handed polarisation in 95 % of the decays.Due to the slightly more central |η| distribution of the leptons from the subsequent t → W q, W → ν decay, the preselection efficiency eµ for these events is typically 10-20 % higher than for SM QCD t t, increasing with m t1 .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 t production is therefore to increase the measured cross-section by R t1 t1 σ t1 t1 , where R t1 t1 is the ratio of eµ values for stoppair and SM t t production, and σ t1 t1 is the stop-pair production cross-section.
Limits were set on stop-pair production by fitting the effective production cross-section R t1 t1 σ t1 t1 multiplied by a signal strength µ to the difference between the measured cross-sections (σ t t) and the theoretically predicted SM QCD production cross-sections (σ theo t t ).The two datasets were fitted simultaneously, assuming values of σ theo t t = 177.3+11.5 −12.0 pb for √ s = 7 TeV and 252.9 +15.3 −16.3 pb for √ s = 8 TeV, including the uncertainty due to a ±1 GeV variation in the top quark mass.The limits were determined using a profile likelihood ratio in the asymptotic limit [64], using nuisance parameters to account for correlated theoretical and experimental uncertainties.
The observed and expected limits on µ at the 95 % confidence level (CL) were extracted using the CLs prescription [65] and are shown in Fig. 9. Due to the rapidly decreasing stop-pair production cross-section with increasing m t1 , 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 [66]), stop masses between the top mass threshold and 177 GeV are excluded, assuming 100 % branching ratio for t1 → t χ0 1 and m χ0 1 = 1 GeV.The limits from considering the √ s = 7 TeV and √ s = 8 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 t1 = 200 GeV for m χ0 1 = 20 GeV.However, they depend more strongly on the assumed top quark polarisation; in a scenario with m t1 = 175 GeV and m χ0 1 = 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 t1 = 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 t1 → bW χ0 1 for [GeV]

9.
Expected and observed limits at 95 % CL on the signal strength µ as a function of mt 1 , for pair produced top squarks t1 decaying with 100 % branching ratio via t1 → t χ0 1 to predominantly right-handed top quarks, assuming m χ0 1 = 1 GeV.The black dotted line shows the expected limit with ±1σ 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.m t1 < m t , the b-jets become softer than those from SM t t production, affecting the determination of b .Therefore, no limits are set for scenarios with m t1 < m t .

Conclusions
The inclusive t t production cross-section has been measured at the LHC using the full ATLAS 2011-2012 pp collision data sample of 4.6 fb −1 at √ s = 7 TeV and 20.3 fb −1 at √ s = 8 TeV, in the dilepton t t → eµννb b decay channel.The numbers of opposite-sign eµ events with one and two b-tagged jets were counted, allowing a simultaneous determination of the t t cross-section σ t t and the probability to reconstruct and b-tag a jet from a t t decay.Assuming a top quark mass of m t = 172.5 GeV, the results are: σ t t = 182.9± 3.1 ± 4. 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 √ s = 7 TeV and √ s = 8 TeV.The dependence of the results on the assumed value of m t is dσ t t/dm 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 t cross-sections have been used to determine the top quark pole mass via the dependence of the predicted cross-section on m pole t , giving a value of m pole t = 172.9+2.5 −2.6 GeV, compatible with the mass measured from kinematic reconstruction of t t events.
The results have also been used to search for pairproduced supersymmetric top squarks decaying to top quarks and light neutralinos.Assuming 100 % branching ratio for the decay t1 → t χ0 1 , 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.

Fig. 1 .
Fig. 1.Distributions of the number of b-tagged jets in preselected opposite-sign eµ events in (a) √ s = 7 TeV and (b) √ s = 8 TeV data.The data are shown compared to the expectation from simulation, broken down into contributions from t t, W t 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 t signal samples generated with Powheg + Pythia6 (PY), MC@NLO + Herwig (HW) and Alpgen + Herwig, and with the cyan band indicating the statistical uncertainty.

Fig. 2 .
Fig.2.Distributions of (a) the number of jets, (b) the transverse momentum pT of the b-tagged jets, (c) the |η| of the electron, (d) the pT of the electron, (e) the |η| of the muon and (f) the pT of the muon, in events with an opposite-sign eµ pair and at least one b-tagged jet.The √ s = 7 TeV data are compared to the expectation from simulation, broken down into contributions from t 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 t signal samples and with the cyan band indicating the statistical uncertainty.The last bin includes the overflow.

Fig. 3 .
Fig.3.Distributions of (a) the number of jets, (b) the transverse momentum pT of the b-tagged jets, (c) the |η| of the electron, (d) the pT of the electron, (e) the |η| of the muon and (f) the pT of the muon, in events with an opposite-sign eµ pair and at least one b-tagged jet.The √ s = 8 TeV data are compared to the expectation from simulation, broken down into contributions from t 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 t signal samples and with the cyan band indicating the statistical uncertainty.The last bin includes the overflow.

Fig. 4 .
Fig.4.Distributions of electron and muon |η| and pT in same-sign eµ events at √ s = 8 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 pT distributions, the last bin includes the overflows.

Fig. 5 .
Fig. 5. Measured t t cross-section at √ s = 8 TeV as a function of the b-tagged jet pT cut.The error bars show the uncorrelated part of the statistical uncertainty with respect to the baseline measurement with jet pT > 25 GeV.

Fig. 6 .
Fig.6.Measurements of the t t cross-section at √ s = 7 TeV and √ s = 8 TeV from this analysis (eµ b-tag) together with previous ATLAS results at √ s = 7 TeV using the ee, µµ and eµ channels[51] and using a fit to jet multiplicities and missing transverse momentum in the eµ channel[52].The √ s = 7 TeV measurements are displaced slightly for clarity.The NNLO+NNLL prediction[2,3] described in Sect. 2 is also shown as a function of √ s, for fixed mt = 172.5 GeV and with the uncertainties from PDFs, αs and QCD scale choices indicated by the green band.

Fig. 7 .
Fig. 7. Predicted NNLO+NNLL t t production cross-sections at √ s = 7 TeV and √ s = 8 TeV as a function of m pole t , showing the central values and total uncertainty bands with several PDF sets.The yellow band shows the QCD scale uncertainty.The measurements of σ t t are also shown, with their dependence on the assumed value of mt through acceptance and background corrections parameterised using Eq.(2).

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µ event samples at √ s = 7 TeV and √ s = 8 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.

Table 3 .
Detailed breakdown of the symmetrised relative statistical, systematic and total uncertainties on the measurements of the t t production cross-section σ t t at √ s = 7 TeV and √ s = 8 TeV.Uncertainties quoted as '0.00' are smaller than 0.005, whilst '-' indicates the corresponding uncertainty is not applicable.The uncertainties on eµ and C b are also shown, with their relative signs indicated where relevant.They contribute with opposite signs to the uncertainties on σ t t, which also include uncertainties from estimates of the background terms N bkg 1 and N bkg 2 .The lower part of the table gives the systematic uncertainties that are different for the measurement of the fiducial cross-section σ fid t t , together with the total analysis systematic and total uncertainties on σ fid t t .

Table 4 .
Summary of the relative statistical, systematic and total uncertainties on the measurements of the t t production cross-section σ t t at

Table 5 .
Fiducial cross-section measurement results at

Table 6 .
Measurements of the top quark pole mass determined from the t t cross-section measurements at Comparison of top quark pole mass values determined from this and previous cross-section measurements