Searches for scalar leptoquarks in pp collisions at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\sqrt{s}}$$\end{document}s = 8 TeV with the ATLAS detector

Searches for pair-produced scalar leptoquarks are performed using 20 fb\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{-1}$$\end{document}-1 of proton–proton collision data provided by the LHC and recorded by the ATLAS detector at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s}=8$$\end{document}s=8 TeV. Events with two electrons (muons) and two or more jets in the final state are used to search for first (second)-generation leptoquarks. The results from two previously published ATLAS analyses are interpreted in terms of third-generation leptoquarks decaying to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$b\nu _{\tau }\bar{b}\bar{\nu _{\tau }}$$\end{document}bντb¯ντ¯ and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$t\nu _{\tau }\bar{t}\bar{\nu _{\tau }}$$\end{document}tντt¯ντ¯ final states. No statistically significant excess above the Standard Model expectation is observed in any channel and scalar leptoquarks are excluded at 95 % CL with masses up to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_{\mathrm {LQ1}}<$$\end{document}mLQ1< 1050 GeV for first-generation leptoquarks, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_{\mathrm {LQ2}}<$$\end{document}mLQ2< 1000 GeV for second-generation leptoquarks, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_{\mathrm {LQ3}}<$$\end{document}mLQ3< 625 GeV for third-generation leptoquarks in the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$b\nu _{\tau }\bar{b}\bar{\nu _{\tau }}$$\end{document}bντb¯ντ¯ channel, and 200 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$< m_{\mathrm {LQ3}}< $$\end{document}


Introduction
Leptoquarks (LQ) are predicted by many extensions of the Standard Model (SM) [1][2][3][4][5][6][7] and may provide an explanation for the many observed similarities between the quark and lepton sectors of the SM. LQs are colour-triplet bosons with fractional electric charge. They carry non-zero values of both baryon and lepton number [8]. They can be either scalar or vector bosons and are expected to decay directly to lepton-quark pairs (where the lepton can be either charged or neutral).
The coupling strength between scalar LQs and the leptonquark pairs depends on a single Yukawa coupling, termed λ LQ→ q . The additional magnetic moment and electric quadrupole moment interactions of vector LQs are governed by two coupling constants [9]. The coupling constants for both the scalar and vector LQs, and the branching fraction of the LQ decay into a quark and a charged lepton, β, are model dependent. The production cross-section and couplings of vector LQs are enhanced relative to the contribution of scalar LQs, although the acceptance is expected to be similar in both cases. This analysis considers the simpler scenario of scalar LQ pair-production, for which the form of the interaction is known and which provides more conservative limits on LQ pair-production than for vector LQ pair-production.
In proton-proton collisions, LQs can be produced singly and in pairs. The production of single LQs, which happens at hadron colliders in association with a lepton, depends directly on the unknown Yukawa coupling λ LQ→ q . However, LQ pair-production is not sensitive to the value of the coupling. In pp collisions with a centre-of-mass energy √ s = 8 TeV, the dominant pair-production mechanism for LQ masses below ∼1 TeV is gluon fusion, while the qq-annihilation production process becomes increasingly important with increasing LQ mass.
The minimal Buchmüller-Rückl-Wyler model (mBRW) [10] is used as a benchmark model for scalar LQ production. It postulates additional constraints on the LQ properties, namely that the couplings have to be purely chiral, and makes the assumption that LQs are grouped into three families (first, second and third-generation) that couple only to leptons and quarks within the same generation. The latter requirement excludes the possibility of flavour-changing neutral currents (FCNC) [11], which have not been observed to date.
Previous searches for pair-produced LQs have been performed by the ATLAS Collaboration with 1.03 fb −1 of data collected at √ s = 7 TeV, excluding at 95 % confidence level (CL) the existence of scalar LQs with masses below 660 (607) GeV for first-generation LQs at β = 1 (0.5) [12] and 685 (594) GeV for second-generation LQs at β = 1 (0.5) [13]. The CMS Collaboration excluded at 95 % CL the existence of scalar LQs with masses below 830 (640) GeV for first-generation LQs at β = 1 (0.5) and 840 (650) GeV for second-generation LQs at β = 1 (0.5) with 5.0 fb −1 of data collected at √ s = 7 TeV [14]. Pair-produced third-generation scalar LQs decaying to bν τbντ have been excluded by the CMS Collaboration for masses below 700 GeV at β = 0, and for masses below 560 GeV over the full β range using 19.7 fb −1 of data collected at √ s = 8 TeV [15]. Third-generation scalar LQs have been excluded in the bτ +b τ − channel at β = 1 for masses up to 740 GeV by the CMS Collaboration using 19.7 fb −1 of data collected at √ s = 8 TeV [16], and by the ATLAS Collaboration at β = 1 for masses up to 534 GeV using 4.7 fb −1 of data collected at √ s = 7 TeV [17]. The CMS Collaboration also excluded third-generation scalar LQs in the tτ −t τ + channel at β = 1 for masses up to 685 GeV using 19.7 fb −1 of data collected at √ s = 8 TeV [15].
In this paper, searches for pair-produced first-and secondgeneration scalar LQs (LQ1 and LQ2, respectively) are performed by selecting events with two electrons or muons plus two jets in the final state (denoted by eejj and μμjj, respectively). In addition, limits are placed on pair-produced third-generation scalar LQs (LQ3) by reinterpreting ATLAS searches for supersymmetry (SUSY) in two different channels [18,19]. LQ production and decay mechanisms can be similar to those of stop quarks (t) and sbottom quarks (b). For example,tt → ttχ 0χ 0 gives the same event topology as LQ3 LQ3 → tν τtντ in the limit where the neutralino (χ 0 ) is massless. Two ATLAS analyses optimised for these SUSY processes are therefore used to set limits on the equivalent LQ decay processes: LQ3 LQ3 → bν τbντ and LQ3 LQ3 → tν τtντ .
The results for each LQ3 channel cannot be combined since the parent LQs have different electric charges in the two cases (− 1 3 e for the LQ3 LQ3 → bν τbντ channel and 2 3 e for the LQ3 LQ3 → tν τtντ channel, where e is the elementary electric charge). The branching fractions of LQ3 decays to bν τ and tν τ are assumed to be equal to 100 % in each case. Although complementary decays of a charge − 1 3 e ( 2 3 e) LQ into a tτ −t τ + (bτ +b τ − ) final state are also allowed, kinematic suppression factors which favour LQ decays to bquarks over t-quarks and the relative strengths of the Yukawa couplings would have to be considered. Since these suppression factors are model dependent, limits are not provided as a function of β for the LQ3 channels.

The ATLAS detector
The ATLAS experiment [20] is a multi-purpose detector with a forward-backward symmetric cylindrical geometry and nearly 4π coverage in solid angle. The three major sub-components of ATLAS are the tracking detector, the calorimeter and the muon spectrometer. Charged-particle tracks and vertices are reconstructed by the inner detector (ID) tracking system, comprising silicon pixel and microstrip detectors covering the pseudorapidity 1 range |η| < 2.5, and a straw tube tracker that covers |η| < 2.0. The ID is immersed in a homogeneous 2 T magnetic field provided by a solenoid. Electron, photon, jet and tau energies are measured with sampling calorimeters. The ATLAS calorimeter system covers a pseudorapidity range of |η| < 4.9. Within the region |η| < 3.2, electromagnetic calorimetry is provided by barrel and endcap high-granularity lead/liquid argon (LAr) calorimeters, with an additional thin LAr presampler covering |η| < 1.8, to correct for energy loss in material upstream of the calorimeters. Hadronic calorimetry is provided by a steel/scintillatortile calorimeter, segmented into three barrel structures within |η| < 1.7, and two copper/LAr hadronic endcap calorimeters. The forward region (3.1 < |η| < 4.9) is instrumented by a LAr calorimeter with copper (electromagnetic) and tungsten (hadronic) absorbers. Surrounding the calorimeters is a muon spectrometer (MS) with air-core toroids, a system of precision tracking chambers providing coverage over |η| < 2.7, and detectors with triggering capabilities over |η| < 2.4 to provide precise muon identification and momentum measurements.

Data and Monte Carlo samples
The results presented here are based on proton-proton collision data at a centre-of-mass energy of √ s = 8 TeV, collected by the ATLAS detector at the LHC during 2012. Data samples corresponding to an integrated luminosity of 20.3fb −1 are used by all channels except for the LQ3 LQ3 → bν τbντ analysis which uses 20.1 fb −1 because of requirements made by the trigger used in the analysis.
Simulated signal events of pair-produced scalar LQs decaying to e + e − qq, μ + μ − qq, tν τtντ , and bν τbντ final states are produced using the Pythia 8.160 [21] event generator with CTEQ6L1 [22] parton distribution functions (PDFs). The coupling λ LQ→ q which determines the LQ lifetime and width [23] is set to √ 0.01 × 4πα, where α is the finestructure constant. This value gives the LQ a full width of less than 100 MeV, which is much smaller than the detector resolution. For LQ masses in the ranges considered here (200 GeV ≤ m LQ ≤ 1300 GeV, in steps of 50 GeV), the value of the coupling used is such that the LQs can be considered to decay promptly. The production cross-section of pairproduced LQs is assumed to be independent of the coupling strength. The signal process is normalised to the expected next-to-leading-order (NLO) cross-sections for scalar LQ pair-production [24]. The signal production cross-section is 23.5 fb for a LQ mass of 600 GeV, and 0.40 fb for a 1 TeV LQ and is the same for each generation.

Monte Carlo for background predictions
The Monte Carlo (MC) samples used to estimate the contributions from SM backgrounds to the LQ1 and LQ2 searches are discussed here. Details about the MC models used for estimating backgrounds in the LQ3 searches are available in Refs. [18] (for the bν τbντ channel) and [19] (for the tν τtντ channel).
The MC samples used to model the Z /γ * +jets background with a dilepton invariant mass (m ) less than 120 GeV and high-mass Drell-Yan backgrounds (m ≥ 120 GeV) are generated with SHERPA 1.4.1 [25]. The high-mass Drell-Yan samples are generated assuming massive c-and b-quarks instead of the conventional massless treatment.
Samples of tt events are produced with POWHEG box [26,27] interfaced with PYTHIA 6. MC samples representing the W W , W Z, and Z Z diboson decays are generated with HERWIG 6.52 [28] and use the AUET2 [29] values for the tunable parameters (the 'AUET2 MC tune'). Samples of single-top-quark events in the W t and s-channel are gener-ated with MC@NLO 4.01 [30,31] and the AUET2 MC tune, while the t-channel samples are generated with AcerMC 3.8 [32] interfaced with PYTHIA 8 and use the AUET2B [33] MC tune. The hadronisation and parton showering of the samples produced with MC@NLO are done using HERWIG 6.52 coupled to JIMMY 4.31 [34]. The W +jets samples are produced with ALPGEN 2.14 interfaced with JIMMY 4.31, also with the AUET2 MC tune applied. The choice of PDFs used to produce the MC simulated samples is generator dependent: AcerMC, PYTHIA, HERWIG and ALPGEN use CTEQ6L1, while MC@NLO uses CT10 [35]. For all samples, the detector response is modelled [36] using GEANT4 [37], except for the Drell-Yan background samples, which use a fast detector simulation where the calorimeter response is parameterised. The differences between fast and full simulation in terms of kinematic spectra and modelling of relevant objects are evaluated to be negligible.
The cross-sections of background processes used in the analysis are taken from theoretical predictions. Single-top production cross-sections in the s-channel [38], t-channel [39], and in associated production with a W boson [40], are calculated to NLO+NNLL accuracy. W +jets and Z → τ τ cross-sections with NNLO accuracy are used [41]. The crosssections for W W , W Z, and Z Z processes are calculated at NLO [42,43]. The theoretical cross-section for W W production is scaled by a factor 1.2 and the uncertainty is increased by an extra 20 %, in order to take into account the ATLAS [44] and CMS measurements [45], which showed an excess in data at the level of 20 % (see Refs. [46,47] for more discussion about possible causes of the excess).
For the Z /γ * +jets and tt backgrounds, LO and NLO cross-sections, respectively are used. These backgrounds are constrained using two control regions (CRs), as described in Sect. 4.5.

Searches for first-and second-generation LQs
The first-and second-generation analyses exploit similarities in the final states and use common search strategies to select dilepton plus dijet final states. Control regions are used to constrain estimates of the dominant backgrounds to the data. A set of discriminating variables is used to define signal regions (SRs) that are used for a counting analysis.

Trigger and data collection
Selected data events are required to have all relevant components of the ATLAS detector in good working condition. For the LQ1 (eejj) analysis, the trigger requires at least two electromagnetic calorimeter clusters, defined as energy deposits in the cells of the electromagnetic calorimeter. The leading cluster is required to have transverse momentum p T > 35 GeV and the sub-leading one p T > 25 GeV. This trigger selects electrons without imposing any requirement on the isolation and this allows a data-driven estimate of the background contribution from jets in the final state that pass the electron selection, as described in detail in Ref. [48]. The trigger is 98 % efficient with respect to the offline selection, which requires p T above 40 (30) GeV for the leading (sub-leading) electron.
For the LQ2 (μμjj) analysis, events are selected from data using a trigger which requires the presence of at least one muon candidate in the event with p T above 36 GeV. This trigger is fully efficient relative to the offline selection for muons with p T above 40 GeV [49].

Object selection
Electrons are selected and identified by imposing requirements on the shape of the cluster of energy deposits in the calorimeter, as well as on the quality of the track, and on the track-to-cluster matching. The identification efficiency is on average 85 % [50]. Electron candidates must have transverse energy E T > 30 GeV and |η| < 2.47. Electron candidates associated with clusters in the transition region between the barrel and endcap calorimeters (1.37 < |η| < 1.52) are excluded. All electrons are required to be reconstructed with cluster-based or combined cluster-and track-based algorithms and to satisfy calorimeter quality criteria. Requirements are made on the transverse (|d 0 |) and longitudinal (|z 0 |) impact parameters of the electron relative to the primary vertex and must satisfy |d 0 | < 1 mm and |z 0 | < 5 mm. In addition, electrons are required to be isolated by imposing requirements on the E R<0.2 T measured in the calorimeter within a cone of size R = ( η) 2 + ( φ) 2 = 0.2 around the electron cluster excluding the electron cluster energy, and corrected to account for leakage (i.e. energy deposited by the electron outside of the cluster) and the average number of proton-proton interactions per bunch-crossing. The isolation requirements are optimised for highp T electrons following the strategy in Ref. [48]. The leading electron is required to have E R<0.2 T < 0.007 × E T + 5 GeV, and the sub-leading electron is required to have E R<0.2 Muon tracks are reconstructed independently in the ID and the MS. Tracks are required to have a minimum number of hits in each system, and must be compatible in terms of geometrical and momentum matching. In particular, in order to prevent mis-measurements at high p T , muons are required to have hits in all three MS stations, as described in Ref. [48]. In order to increase the muon identification efficiency, when one muon in the event satisfies the three-stations requirement, the criteria for the second muon in the event are relaxed to require hits in only two MS stations. Information from both the ID and MS is used in a combined fit to refine the measure-ment of the momentum of each muon [51]. Muon candidates are required to have p T > 40 GeV, |η| < 2.4, |d 0 | < 0.2 mm and |z 0 | < 1.0 mm. Muons must also pass a relative-isolation requirement T is the sum of the transverse momenta of all the tracks with p T above 1 GeV (except for the muon track) within a cone of R < 0.2 around the muon track, and p T is the transverse momentum of the muon.
Jets are reconstructed from clusters of energy deposits detected in the calorimeter using the anti-k t algorithm [52] with a radius parameter R = 0.4 [53]. They are calibrated using energy-and η-dependent correction factors derived from simulation and with residual corrections from in-situ measurements. The jets used in the analysis must satisfy p T > 30 GeV and |η| < 2.8. Jets reconstructed within a cone of R = 0.4 around a selected electron or muon are removed. Additional jet quality criteria are also applied to remove fake jets caused by detector effects. A detailed description of the jet energy scale measurement and its systematic uncertainties is given in Ref. [54].

Event pre-selection
Multiple pp interactions during bunch-crossings (pile-up) can give rise to multiple reconstructed vertices in events. The primary vertex of the event, from which the leptons are required to originate, is defined as the one with the largest sum of squared transverse momenta of its associated tracks. Events are selected if they contain a primary vertex with at least three associated tracks satisfying p T,track > 0.4 GeV.
MC events are corrected to better describe the data by applying a per-event weight to match the distribution of the average number of primary vertices observed in data. A weighting factor is also applied in order to improve the modelling of the vertex position in z. Scale factors are applied to account for differences in lepton identification and selection efficiency between data and MC simulation. The scale factors depend on the lepton kinematics and are described in detail in Ref. [51] for muons, and in Ref. [55] for electrons. The energy and momentum of the selected physics objects are corrected to account for the resolution and scale measured in data, as described in Ref. [51] for muons, in Ref. [55] for electrons and in Ref. [54] for jets.
Events are selected in the eejj channel if they contain exactly two electrons with p T > 40 (30) GeV for the leading (sub-leading) electron and at least two jets with p T > 30 GeV. For the μμjj channel, events are selected if they contain exactly two muons with p T > 40 GeV and opposite-sign charge, and at least two jets with p T > 30 GeV. No requirements are placed on the charges of the electron candidates due to inefficiencies in determining the charge of highp T tracks associated with electrons. These sets of requirements form the basic event 'pre-selection' for the analyses, which is used to build the control and signal regions discussed in the following sections.

Signal regions
After applying the event pre-selection requirements, a set of signal regions is defined using additional kinematic variables in order to discriminate LQ signals from SM background processes and to enhance the signal-to-background ratio. The variables used are: Signal regions are determined by optimising the statistical significance as defined in Ref. [56]. The optimisation procedure is performed in a three-dimensional space constructed by m , S T and m min LQ for each of the signal mass points. Several adjacent mass points may be grouped into a single SR. The signal acceptance of the selection requirements is estimated to be ≈50 % in the μμjj channel and between 65 and 80 % in the eejj channel (assuming β = 1.0). The difference is due to tighter quality selection requirements in the μμjj channel used to prevent muon mis-measurements in MS regions with poor alignment or missing chambers. The optimised signal regions are presented in Table 1 together with the mass of the corresponding LQ hypothesis. Each LQ mass hypothesis is tested in only one signal region, where limits on σ × β are extracted. Table 1 The minimum values of m , S T , and m min LQ used to define each of the signal regions targeting different LQ masses in the eejj and μμjj channels. Each signal region is valid for one or more mass hypotheses, as shown in the second column

Background estimation
The main SM background processes to the LQ1 and LQ2 searches are the production of Z /γ * +jets events, tt events where both top quarks decay leptonically, and diboson events. Additional small contributions are expected from Z → τ τ and single-top processes. Multi-jets, W +jets, tt (where one or more top quarks decays hadronically), and single-top events with mis-identified or non-prompt leptons arising from hadron decays or photon conversions can also contribute. These fake lepton backgrounds are estimated separately in the eejj and μμjj channels using the same datadriven techniques as described in Ref. [48] and are found to be negligible for the μμjj channel. Normalisation factors, derived using background-enriched control regions, are applied to the MC predictions for Z /γ * +jets and tt backgrounds to predict as accurately as possible the background in the signal regions. These control regions are constructed to be mutually exclusive to the signal region and the assumption is made that normalisation factors and their associated uncertainties in the signal region are the same as in the backgroundenriched control regions.

Control regions for Z /γ * +jets and tt backgrounds
Two control regions with negligible signal contributions are defined to validate the modelling accuracy of the MC simulated background events and to derive normalisation scale factors. The Z /γ * +jets control region is defined by the pre-selection requirements with an additional requirement of 60 < m ee < 120 GeV in the eejj channel and 70 < m μμ < 110 GeV in the μμjj channel. These control regions define a pure sample of Z /γ * +jets events. The tt control region is defined in both channels by applying the preselection requirements, but demanding exactly one muon and one electron (both with p T above 40 GeV) in the offline selection instead of two same-flavour leptons. In the case of the tt control region for the eejj channel, the trigger requirement is modified by requiring a single isolated electron with p T above 24 GeV, which is fully efficient relative to the offline selection for electrons with p T above 30 GeV. In both cases, the same selection criteria are applied to data and MC events. Normalisation factors are applied to the MC predictions for the Z /γ * +jets and tt background processes. They are obtained by performing a combined maximum likelihood fit to the observed yields in the control regions and signal region under consideration. Systematic uncertainties on the predicted MC yields related to the uncertainty on the crosssections are taken into account by the fit through the use of dedicated nuisance parameters. The fit procedure is performed using the HistFitter package [57], which is a tool based on the RooStats framework [58]. The normalisation scale factor obtained from a background-only fit for the Z /γ * +jets background in the eejj (μμjj) channel is 1.1±0.2 (0.97 ± 0.15), while the normalisation scale factor for tt is 1.10±0.05 (1.01±0.05). The fitted background scale factors have little sensitivity to the inclusion of signal regions and the eventual presence of a signal.

Kinematic distributions
The distributions of the kinematic variables after performing the background-only fits in the control regions, and applying the event pre-selection requirements are shown in Figs. 1, 2 and 3 for the data, background estimates, and for three LQ masses of 300, 600 and 1000 GeV (with β = 1.0).

Systematic uncertainties
The theoretical uncertainty on the NLO cross-section is taken into account for diboson, single-top, W +jets, and Z → τ τ processes. For the two dominant backgrounds (tt and Z /γ * +jets) the modelling uncertainties are estimated using the symmetrised deviation from unity of the ratio of data to MC events in the tt and Z /γ * +jets control regions, which is fitted with a linear function for S T > 400 GeV. The modelling systematic uncertainty is then applied as a function of S T , in the form of a weighting factor. The choice of S T for such a purpose is motivated by its sensitivity to mismodelling of the kinematics of jets and leptons. It varies in the eejj (μμjj) channel between 8 % (10 %) and 25 % (30 %) for the Z /γ * +jets background and between 6 % (10 %) and 24 % (40 %) for the tt background. It increases for signal regions targeting higher m LQ .
The jet energy scale (JES) uncertainty depends on p T and η and contains additional factors, which are used to correct for pile-up effects. They are derived as a function of the number of primary vertices in the event to take into account additional pp collisions in a recorded event (in-time pile-up), or as a function of the expected number of interactions per bunchcrossing to constrain past and future collisions affecting the measurement of energies in the current bunch-crossing (outof-time pile-up). An additional uncertainty on the jet energy resolution (JER) is taken into account. The relative impact on the background event yields from the JES (JER) uncertainty is between 8 % (1 %) in SR1 and 26 % (1 %) in SR8. The signal selection efficiency change due to the JES uncertainties ranges between 3 % in SR1 and 1 % in SR8, while the effect of the JER is negligible.
The electron energy scale and resolution are corrected to provide better agreement between MC predictions and data. The uncertainties on these corrections are propagated through the analysis as sources of systematic uncertainty. Uncertainties are taken into account for the electron trigger (∼0.1 %), identification (∼1 %) and reconstruction (∼1 %) efficiencies, and for uncertainties associated with the isolation requirements (∼0.1 %).
Scaling and smearing corrections are applied to the p T of the muons in order to minimise the differences in resolution between data and MC simulated events. The uncertainty on these corrections is below 1 %. Differences in the identification efficiency and in the efficiency of the trigger selection are taken into account and are less than 1 %.
QCD renormalisation and factorisation scales are varied by a factor of two to estimate the impact of higher orders on Events / 20 GeV The ratio of the number of data events to the number of background events (and its statistical uncertainty) is also shown. The hashed bands represent all sources of statistical and systematic uncertainty on the background prediction the signal production cross-section. The variation is found to be approximately 14 % for all mass points. The uncertainty on the signal cross-section related to the choice of PDF set is evaluated as the envelope of the prediction of 40 different CTEQ6.6 NLO error sets [24]. The uncertainty ranges from 18 % at m LQ = 300 GeV to 56 % at m LQ = 1300 GeV. These uncertainties are the same for all LQ generations. The effect on the choice of PDF set on the signal acceptance times reconstruction efficiency is estimated using the Hes-sian method [59]. The final PDF uncertainties on the signal samples are approximately 1 % for most mass points, rising to 4 % for some higher LQ masses. The impact of the choice of PDF set on the acceptance times reconstruction efficiency for each background process is estimated using the Hessian method (using the same method as for signals). The uncertainties range from 4 % in the low-mass signal regions to 17 % in the high-mass signal regions.

Results
The observed and expected yields in three representative signal regions for the eejj and the μμjj channels after the combined maximum likelihood fits are shown in Tables 2 and  3, respectively. The fit maximizes the likelihood constructed using the two CRs and the SR under study. When contructing the likelihood, the signal stregth and the background scale factors are treated as free parameters, the systematic uncertainties are treated as nuisance parameters. No significant excess above the SM expectation is observed in any of the signal regions and a modified frequentist CL s method [60] is used to set limits on the strength of the LQ signal, by constructing a profile likelihood ratio. Pseudo-experiments are used to determine the limits.
The cross-section limits on scalar LQ pair-production are presented as a function of β for both channels in Fig. 4. Also shown are the results of the ATLAS searches for firstand second-generation LQs using 1.03 fb −1 data at √ s= 7 TeV which also included searches in the eνjj and μνjj decay channels and therefore provide better sensitivity at low values of β. First (second)-generation scalar LQs are excluded for β = 1 at 95 % CL for m LQ1 < 1050 GeV (m LQ2 < 1000 GeV). The expected exclusion ranges are the same as the observed ones. First (second)-generation scalar LQs are excluded for m LQ1 < 650 GeV (m LQ2 < 650 GeV) at β = 0.2 and m LQ1 < 900 GeV (m LQ2 < 850 GeV) at β = 0.5.

Search for third-generation LQs in the bν τbντ channel
The ATLAS search for pair-production of third-generation supersymmetric partners of bottom quarks (sbottom,b) [18] is reinterpreted in terms of the LQ model, in the case where each LQ decays to a b-quark and a ν τ neutrino. In the original analysis, theb is assumed to decay viab→ bχ 0 , and

2-muons + 2-jets
All limits at 95% CL jj ν μ jj+ μ μ expected limit observed limit σ 1 ± expected σ 2 ± expected -1 = 7 TeV, 1.03 fb s (d) Fig. 4 The cross-section limits on scalar LQ pair-production times the square of the branching ratio as a function of mass (left) and the excluded branching ratio as a function of the LQ mass (right) to eq for the eejj channel (top) and to μq for the μμjj channel (bottom). The ±1(2)σ uncertainty bands on the expected limit represent all sources of systematic and statistical uncertainty. The expected NLO production cross-section (β = 1.0) for scalar LQ pair-production and its corresponding theoretical uncertainty due to the choice of PDF set and renormalisation/factorisation scale are also included. The exclusion limits on LQ1 [12] and LQ2 [13] set by ATLAS in the eejj +eνjj and μμjj +μνjj search channels using 1.03 fb −1 of data collected at √ s = 7 TeV are also showñ t viat→ bχ ± in the case where mχ± − mχ0 is small and theχ ± decay products are undetectable. The search is performed for final states with large missing transverse momentum (p miss T , with magnitude E miss T ) and two jets identified as originating from b-quarks. The full analysis strategy is covered in Ref. [18]. A complete description of the analysis, including treatment of systematic uncertainties on background processes can be found there, but the event selection and background estimation methods used are summarised here for clarity.

Object and event selection
Events are required to have exactly two b-tagged [61] jets with p T > 20 GeV and |η| < 2.5, and E miss T > 150 GeV. Additional jets in the event are accepted if they have p T > 20 GeV and |η| < 4.9. Events with one or more electrons (muons) with p T > 7 (6) GeV are vetoed. Candidate signal events are selected from data using a E miss T trigger which is 99 % efficient for events passing the offline selection. Several variables are defined and used to optimise the event selection:    m CT : The contransverse mass [62], used to measure the masses of pair-produced heavy particles that decay semiinvisibly (i.e. decays where one of the decay products can be detected, but the other cannot).
In the original analysis, different signal regions were optimised according to the masses of the third-generation squark and the lightest supersymmetric particle (LSP). In the case of the LQ model reinterpretation the signal regions corresponding to the case where the mass of the LSP is approximately zero have best sensitivity, but all the signal regions are retained for coherence with the original analysis. The different signal region definitions are given in Table 4. Signal region A (SRA) has five different m CT thresholds. Signal region B (SRB) is optimised towards the region where the squark and LSP masses are approximately equal. The signal region with the best expected limit is used for each point in the exclusion plots.

Background estimation
The dominant background process is the production of Z bosons in association with heavy-flavour jets where the Z boson subsequently decays to two neutrinos [Z (→ νν) + bb]. Its contribution is estimated from data in an oppositesign dilepton control region. Top quark pair-production (tt) and W bosons produced in association with heavy flavour quarks also contribute significantly and are normalised in dedicated control regions before being extrapolated to the signal regions using MC simulation. Different control regions are defined for each signal region, requiring one or two leptons plus additional requirements similar to the corresponding signal region. The contributions from Z +jets, W +jets, and top quark production are estimated simultaneously with a profile likelihood fit to the three control regions. Contributions from diboson and tt+W/Z processes are estimated from MC simulation in all regions. The contribution from multi-jet events is estimated from data by taking well-measured multijet events from data and smearing the jets with jet response functions taken from MC simulation and validated in data. This procedure is described in detail in Ref. [63]. The contribution from multi-jet events in signal regions is found to be negligible.

Results
The number of data events observed in each signal region is reported in Table 5, together with the SM background expectation after the background-only fit, and the expected number of signal events for different LQ masses. The signal acceptance efficiency is around 2 % for all but the lowest LQ masses targeted (dropping to 0.27 % efficiency for m LQ = 200 GeV). All sources of systematic and statistical uncertainty are taken into account. The dominant systematic uncertainties on the background prediction are the jet energy scale (JES 1-5 %) and resolution (JER 1-8 %), and the b-tagging uncertainty (2-10 %). Detector-related systematic uncertainties on the signal prediction are dominated by uncertainties on the btagging efficiency (∼30 %). The second-largest source of uncertainty is due to the JES and is around 3 %. The expected (dashed) and observed (solid) 95 % CL upper limits on third-generation scalar LQ pair-production cross-section times the square of the branching ratio to bν τ as a function of LQ mass, for the bν τbντ channel. The ±1(2)σ uncertainty bands on the expected limit represent all sources of systematic and statistical uncertainty. The expected NLO production cross-section (β = 0.0) for scalar LQ pair-production and its corresponding theoretical uncertainty due to the choice of PDF set and renormalisation/factorisation scale are also included The uncertainties on the signal production cross-section are estimated using the methods described in Sect. 4.6. These uncertainties are the same for all LQ genertions but the uncertainty due to the choice of PDF set varies with m LQ . Since the third-generation analyses consider a different mass range to the first-and second-generation analyses, in this case the uncertainty due to the choice of PDF set ranges from 7.1 % at m LQ = 200 GeV to 30 % at m LQ = 800 GeV. Effects on the acceptance due to the choice of PDF set are negligible.
No significant excess above the SM expectation is observed in any of the signal regions. Figure 5 shows the observed and expected exclusion limits for the scalar LQ3 pair-production scenario obtained by taking, for each signal mass configuration, the signal region with the best expected limit. These limits are obtained using the methods described in Sect. 4.7. These methods compare the observed numbers of events in the signal regions with the fitted background expectation and accounting for signal contamination in the corresponding CRs for a given model. Pair-produced third-generation scalar LQs decaying to bν τbντ are excluded at 95 % CL for m LQ3 < 625 GeV. The expected excluded range is m LQ3 < 640 GeV.

Search for third-generation LQs in the tν τtντ channel
The ATLAS search for pair-production of the supersymmetric partner of the top quark (stop quark,t) [19] is reinterpreted in terms of the LQ model, in the case where each LQ decays τ -veto loose τ particle ID --modified, see [19].
Exclusion setup shape-fit in m T and cut-and-count to a top quark and a ν τ neutrino. The original analysis has dedicated signal regions targetingt decays into tχ 0 and the subsequent semileptonic decay of the tt pair. Events compatible with tt plus extra E miss T are selected with final states containing one isolated lepton, jets, and E miss T . A complete description of the analysis strategy, including the treatment of systematic uncertainties on background processes can be found in Ref. [19]. The event selection and the background estimation methods are summarised here for clarity.

Object and event selection
Events are required to contain exactly one electron with p T > 25 GeV and |η| < 2.47, or muon with p T > 25 GeV and |η| < 2.4. Events containing more than one electron or muon with looser identification and p T requirements (10 GeV for both) are vetoed. In some signal regions, events are vetoed if they are consistent with containing a hadronically decaying τ lepton. Events are required to have a minimum of four jets with p T > 20 GeV and |η| < 2.5, with at least one of these passing b-tagging requirements [61]. In addition, selected events must have E miss T > 100 GeV. Several variables are used to further select signal events and reject background processes: m T : The transverse mass of the electron or muon and the E miss T . am T and m τ T,2 : These are two variants on the stransverse mass (m T,2 ) [64][65][66] which is a generalisation of the transverse mass when applied to signatures with two invisible particles in the final state. The asymmetric stransverse mass am T , aims to reject dileptonic tt events where one of the leptons is not reconstructed or is outside the acceptance (and therefore adds to the E miss T of the event). The second implementation of this variable, the τ stransverse mass m τ T,2 , targets tt events where one top decays leptonically and the other top decays into a τ that subsequently decays hadronically.
topness: This variable is designed to reject dileptonic tt events where one lepton is assumed to be lost, as detailed in Ref. [67]. The topness variable is based on the minimisation of a χ 2 -type function. m had−top : This quantity is used to reject dileptonic tt events but retain signal events that contain a hadronically decaying on-shell top quark, as in the LQ → t + ν τ andt 1 → tχ 0 1 scenarios. -φ(jet 1,2 , p miss T ): The azimuthal opening angle between the leading or sub-leading jet and p miss T used to suppress multi-jet events where p miss -E miss T / √ H T : An approximation of the E miss T significance, where H T is defined as the scalar p T sum of the leading four jets.
-H miss T,sig : An object-based missing transverse momentum, divided by the per-event resolution of the jets, and shifted to the scale of the background [68].
The variables listed above are used to define three cutand-count SRs and one shape-fit SR. Table 6 details the event selections for these signal regions. The SR labelled tN_boost targets LQ/stop masses of 700 GeV and takes advantage of the 'boosted' topology of such a heavy parent particle. The selection assumes that either all decay products of the hadronically decaying top quark, or at least the decay products of the hadronically decaying W boson, collimate into a jet reconstructed with a radius parameter R = 1.0 [69,70].

Background estimation
The dominant sources of background are the production of tt events and W +jets where the W boson decays leptonically. Other background processes considered are single top, dibosons, Z +jets, tt produced in association with a vector boson (tt V ), and multi-jets.
The predicted numbers of tt and W +jets background events in the SRs are estimated from data using a fit to the number of observed events in dedicated control regions. Each SR has an associated CR for each of the tt and W +jets backgrounds. The CRs are designed to select events as similar as possible to those selected by the corresponding SR while keeping the contamination from other backgrounds and potential signal low. This is achieved by e.g. requiring that 60 < m T < 90 GeV and in the case of the W +jets CR, inverting the b-jet requirement so that it becomes a b-jet veto. The simulation is used to extrapolate the background predictions into the signal region. The background fit predictions are validated using dedicated event samples, referred to as validation regions (VRs), and one or more VR is defined for each of these. Most VRs are defined by changing the m T windows to 90 < m T < 120 GeV. The VRs are designed to be kinematically close to the associated SRs to test the background estimates in regions of phase space as similar as possible to the SRs.
The multi-jet background is estimated from data using a matrix method described in Refs. [71,72]. The contribution is found to be negligible. All other (small) backgrounds are determined entirely from simulation and normalised to the most accurate theoretical cross-sections available. Table 7 The number of observed events in the three cut-and-count signal regions, together with the expected number of background events and signal events for different LQ masses (assuming β = 0.0) in the tν τtντ channel [19] tN_med tN_high tN_boost

Results
The number of events observed in each signal region is reported in Tables 7 and 8, together with the SM background expectation and the expected number of signal events for different LQ masses. The signal acceptance is between 1.5 and 3 % depending on the LQ mass. All sources of systematic uncertainty and statistical uncertainty are taken into account. The dominant sources of uncertainty on the background prediction come from uncertainties related to the JES, JER, tt background modelling, the b-tagging efficiency, and statistical uncertainties. Detector-related systematic effects are evaluated for signal using the same methods used for the backgrounds (see Ref. [19] for details). The dominant detector-related systematic effects are the uncertainties on the JES (4 %) and the b-tagging efficiency (3 %).
The uncertainties on the signal production cross-section are estimated using the methods described in Sect. 4.6. The effect on the choice of PDF set on the signal acceptance is less than 1 % for most mass points, but increases to 1.7 % for m LQ = 800 GeV.
Similar methods as described in Sect. 4.7 are used to assess the compatibility of the SM background-only hypothesis with the observations in the signal regions. The observed number of events is found to agree well with the expected number of background events in all signal regions. No significant excess over the expected background from SM processes is observed and the data are used to derive one-sided limits at 95 % CL. The results are obtained from a profile likelihoodratio test following the CL s prescription [60]. The likelihood of the simultaneous fit is configured to include all CRs and one SR or shape-fit bin. The 'exclusion setup' event selection is applied (see Table 6), and all uncertainties except the theoretical signal uncertainty are included in the fit. Table 8 The number of observed events in the shape-fit signal region, together with the expected number of background events and signal events for different LQ masses (assuming β = 0.0) in the tν τtντ channel [19]   The expected (dashed) and observed (solid) 95 % CL upper limits on the third-generation scalar LQ pair-production cross-section times the square of the branching ratio to tν τ as a function of LQ mass, for the tν τtντ channel. The ±1(2)σ uncertainty bands on the expected limit represent all sources of systematic and statistical uncertainty. The expected NLO production cross-section (β = 0.0) for third-generation scalar LQ pair-production and its corresponding theoretical uncertainty due to the choice of PDF set and renormalisation/factorisation scale are also included Exclusion limits are obtained by selecting a priori the signal region with the lowest expected CL s value for each signal grid point. The expected and observed limits on the LQ3 LQ3 → tν τtντ process are shown in Fig. 6. Thirdgeneration scalar LQs decaying to tν τtντ are excluded at 95 % CL in the mass range 210 < m LQ3 < 640 GeV. The expected exclusion range is 200 < m LQ3 < 685 GeV. The limits for stop production in the case where the neutralino is massless are slightly stronger than the limits set on LQ3 pro- Table 9 Expected and observed exclusion ranges at 95 % CL for each of the four LQ decay channels considered This leads to the top quarks being polarised in such a way that the acceptance increases. The limit worsens at low mass, due to the effect of greater contamination from top backgrounds.

Summary and conclusions
Searches for pair-production of first-, second-and thirdgeneration scalar leptoquarks have been performed with the ATLAS detector at the LHC using an integrated luminosity of 20 fb −1 of data from pp collisions at √ s = 8 TeV. No significant excess above the SM background expectation is observed in any channel. The results are summarised in Table 9.
The results presented here significantly extend the sensitivity in mass compared to previous searches. Low-mass regions are also considered and limits on the cross-sections are provided for the different final states analysed. Since β is not constrained by the theory, searches in the low mass regions are also important in order to extract limits for low-β values for the LQ1 and LQ2 analyses.