Search for pair production of first- and second-generation scalar leptoquarks in pp collisions at sqrt(s) = 7 TeV

Results are presented from a search for the pair production of first- and second-generation scalar leptoquarks in proton-proton collisions at sqrt(s) = 7 TeV. The data sample corresponds to an integrated luminosity of 5.0 inverse femtobarns, collected by the CMS detector at the LHC. The search signatures involve either two charged leptons of the same-flavor (electrons or muons) and at least two jets, or a single charged lepton (electron or muon), missing transverse energy, and at least two jets. If the branching fraction of the leptoquark decay into a charged lepton and a quark is assumed to be beta = 1, leptoquark pair production is excluded at the 95% confidence level for masses below 830 GeV and 840 GeV for the first and second generations, respectively. For beta = 0.5, masses below 640 GeV and 650 GeV are excluded. These limits are the most stringent to date.


Introduction
The structure of the standard model (SM) of particle physics suggests a fundamental relationship between quarks and leptons. There are many models beyond the SM that predict the existence of leptoquarks (LQ), hypothetical particles that carry both baryon number and lepton number and couple to both quarks and leptons. Among these scenarios are grand unified theories [1,2], composite models [3], extended technicolor models [4][5][6], and superstring-inspired models [7]. LQs are color triplets with fractional electric charge, and can be either scalar or vector particles. A leptoquark couples to a lepton and a quark with a coupling strength λ, and it decays to a charged lepton and a quark with an unknown branching fraction β or to a neutrino and a quark with branching fraction 1 − β. In order to satisfy constraints from bounds on flavor-changing-neutral-currents and from rare pion and kaon decays [3,8], it is assumed that LQs couple only to quarks and leptons of a single generation. LQs are classified as first-, second-, or third-generation, depending on the generation of leptons they couple to. The dominant mechanisms for the production of LQ pairs at the Large Hadron Collider (LHC) are gluon-gluon (gg) fusion and quark-antiquark (qq) annihilation, shown in Fig. 1. The dominant processes only depend on the strong coupling constant and have been calculated at next-toleading order (NLO) [9]. The cross section for production via the unknown Yukawa coupling λ of an LQ to a lepton and a quark is typically smaller.
This paper reports on a search for pair production of scalar leptoquarks. Several experiments have searched for pair-produced scalar LQs but none has obtained evidence for them [10][11][12][13][14][15]. This search uses a data sample corresponding to an integrated luminosity of 5.0 fb −1 recorded with the Compact Muon Solenoid (CMS) detector during the 2011 proton-proton run of the LHC at √ s = 7 TeV. The analysis performed in this paper considers the decay of leptoquark pairs into two charged leptons of the same flavor (either electrons or muons) and two quarks; or into a charged lepton, a neutrino, and two quarks. As a result, two distinct classes of events are selected: one with two high-transverse-momentum (p T ) electrons or muons and at least two high-p T jets (lljj) and the other with one high-p T electron or muon, large missing transverse energy (E miss T ), and at least two high-p T jets (lνjj). The CMS detector, described in detail elsewhere [16], uses a cylindrical coordinate system with the z axis along the counterclockwise beam axis. The detector consists of an inner tracking system and electromagnetic (ECAL) and hadron (HCAL) calorimeters surrounded by a 3.8 T solenoid. The inner tracking system consists of a silicon pixel and strip tracker, providing the required granularity and precision for the reconstruction of vertices of charged particles in the range 0 ≤ φ ≤ 2π in azimuth and |η| < 2.5, where the pseudorapidity η is defined as η = − ln[tan(θ/2)], and θ is the polar angle measured with respect to the z axis. The crystal 2 Dataset and object reconstruction ECAL and the brass/scintillator sampling HCAL are used to measure with high resolution the energies of photons, electrons, and hadrons for |η| < 3.0. The three muon systems surrounding the solenoid cover a region |η| < 2.4 and are composed of drift tubes in the barrel region (|η| < 1.2), of cathode strip chambers in the endcaps (0.9 < |η| < 2.4), and of resistive plate chambers in both the barrel region and the endcaps (|η| < 1.6). Events are recorded based on a trigger decision using information from either the calorimeter or muon systems. The final trigger decision is based on the information from all subsystems, which is passed on to the highlevel trigger (HLT), consisting of a farm of computers running a version of the reconstruction software optimized for fast processing.
The lljj and lνjj analyses are performed separately and the results are combined as a function of the branching fraction β and the leptoquark mass M LQ for first and second generations independently. The analysis in all four decay channels searches for leptoquarks in an excess of events characteristic of the decay of heavy objects. Various triggers are used to collect events depending on the decay channel and the data taking periods as described in Section 2. An initial selection isolates events with high-p T final-state particles (two or more isolated leptons and two or more jets; or one isolated lepton, two or more jets, and large E miss T indicative of the emission of a neutrino). Kinematic variables are then identified to further separate a possible leptoquark signal from the expected backgrounds, and optimized thresholds on the values of these variables are derived in order to maximize the sensitivity to the possible presence of a signal in each decay mode. The variables used in the optimization are the invariant mass of jet-lepton pairs (M lj ), the scalar sum (S T ) of the p T of each of the final state objects, and either the invariant mass of the dilepton pair (M ll ) in the lljj channels or E miss T in the lνjj channels.
Major sources of SM background are Z/γ * +jets, W+jets processes, and tt. Smaller contributions arise from single top production, diboson processes, and QCD multijet processes. The major backgrounds are determined either from control samples in data or from Monte Carlo (MC) simulated samples normalized to data in selected control regions.
After final selection, the data are well described by the SM background predictions, and upper limits on the LQ cross section are set using a CL S modified frequentist approach [17,18]. Using Poisson statistics, 95% confidence level (CL) upper limits are obtained on the leptoquark pairproduction cross-section times branching fraction as a function of leptoquark mass (M LQ ). This is compared with the NLO predictions [9] to determine lower limits on M LQ for β = 1 and β = 0.5. The lljj and lνjj channels are combined to further maximize the exclusion in β and M LQ , especially for the case β ∼ 0.5, where combining the two channels increases the sensitivity of the search.
The data and the initial event selection are detailed in Section 2 of this paper, followed by a description of the signal modeling and background estimates in Section 3 and Section 4, respectively. Section 5 contains the final event selection, and Section 6 describes the systematic uncertainties. The results of the search are presented in Section 7 and summarized in Section 8.

Dataset and object reconstruction
For the first-generation eejj analysis, events are required to pass a double-electron trigger or a double-photon trigger, with an electron or photon p T > 33 GeV. For the first-generation eνjj analysis, events are required to pass either a single-electron trigger or a trigger based on the requirement of one electron with p T threshold between 17 and 30 GeV, missing transverse energy threshold between 15 and 20 GeV, and two jets with p T threshold between 25 and 30 GeV. The trigger thresholds vary according to the run period. For the second-generation leptoquark analyses, events are required to pass a single-muon trigger without isolation requirements and with a p T threshold of 40 GeV. For the eejj channel, the trigger efficiency is greater than 99%. For the eνjj channel, the electron trigger efficiency is measured to be 95%. For the µµjj and µνjj channels, the single muon trigger efficiency is measured to be 92%.
Electron candidates [19] are required to have an electromagnetic cluster with p T > 40 GeV and pseudorapidity |η| < 2.5 (2.2) for the eejj (eνjj) analysis, excluding the transition region between the barrel and the endcap detectors, 1.44 < |η| < 1.57. The eνjj analysis requires lower electron |η| to reduce the QCD multijet background, with negligible reduction of the signal acceptance. Electron candidates are required to have an electromagnetic cluster in the ECAL that is spatially matched to a reconstructed track in the central tracking system in both η and the azimuthal angle φ, and to have a shower shape consistent with that of an electromagnetic shower. Electron candidates are further required to be isolated from additional energy deposits in the calorimeter and from reconstructed tracks beyond the matched track in the central tracking system. In addition, to reject electrons coming from photon conversion in the tracker material, the track associated with the reconstructed electron is required to have hits in all inner tracker layers.
Muons are reconstructed as tracks in the muon system that are matched to the tracks reconstructed in the inner tracking system [20]. Muons are required to have p T > 40 GeV, and to be reconstructed in the HLT fiducial volume, i.e. with |η| < 2.1. In addition, muons must be isolated by requiring that the tracker-only relative isolation be less than 0.1. Here, the relative isolation is defined as a sum of the transverse momenta of all tracks in the tracker in a cone of ∆R = (∆φ) 2 + (∆η) 2 = 0.3 around the muon track (excluding the muon track), divided by the muon p T . To have a precise measurement of the transverse impact parameter of the muon track relative to the beam spot, only muons with tracks containing more than 10 hits in the silicon tracker and at least one hit in the pixel detector are considered. To reject muons from cosmic rays, the transverse impact parameter with respect to the primary vertex is required to be less than 2 mm.
Jets and E miss T are reconstructed using a particle-flow algorithm [21], which identifies and measures stable particles by combining information from all CMS sub-detectors. The E miss T calculation uses calorimeter estimates improved by high precision inner tracking information as well as corrections based on particle-level information in the event. Jets are reconstructed using the anti-k T [22] algorithm with a distance parameter of R = 0.5. The jet energy is calibrated using p T balance of dijet and γ+jet events [23]. In the eejj and µµjj (eνjj and µνjj) channels, jets are required to have p T > 30 (40) GeV, and |η| < 2.4. Furthermore, jets are required to have a spatial separation from electron or muon candidates of ∆R > 0.3.
The initial selection of eejj or µµjj events requires two electrons or two muons and at least two jets satisfying the conditions described above. The two leptons and the two highest-p T jets are selected as the decay products from a pair of leptoquarks. The invariant mass of the two electrons (muons) is required to be M ll > 60 (50) GeV. To reduce the combinatorial background, events with a scalar transverse energy S ll T = p T (l 1 ) + p T (l 2 ) + p T (j 1 ) + p T (j 2 ) below 250 GeV are rejected.
For the eνjj (µνjj) initial selection, events are required to contain one electron (muon) satisfying the conditions described above and at least two jets with p T > 40 GeV and E miss T > 55 GeV. The jet p T threshold is higher than that in the dilepton channels to account for jet p T thresholds in triggers used in the eνjj channel. A veto on the presence of extra muons (electrons) is also applied. The angle in the transverse plane between the leading p T jet and the E miss T vector is required to be ∆φ(E miss T , j 1 ) > 0.5 to reject events with misreconstructed E miss T . In order to 4 4 Background estimate reduce the contribution from QCD multijet events, the lepton and the E miss T are required to be separated by ∆φ(E miss T , l) > 0.8. In addition, events are rejected if the scalar transverse energy S lν The initial selection criteria are summarized in Table 1. Table 1: Initial selection criteria in the eejj, µµjj, eνjj, and µνjj channels.

Signal and background modeling
The MC samples for the signal processes are generated for a range of leptoquark mass hypotheses between 250 and 900 GeV, with a renormalization and factorization scale µ ≈ M LQ . The MC generation uses the PYTHIA generator [24] (version 6.422) and CTEQ6.6 parton distribution functions (PDF). The MC samples used to estimate the contribution from SM background processes are tt+jets events, generated with MADGRAPH [25,26]; single-top events (s, t, and tW channels), generated with MADGRAPH; Z/γ * +jets events and W+jets events, generated with SHERPA [27]; VV events, where V either represents a W or a Z boson, generated with PYTHIA; QCD muon-enriched multijet events, generated with PYTHIA in bins of transverse momentum of the hard-scattering process from 15 GeV to the kinematic limit. The simulation of the CMS detector is based on GEANT4 [28], and includes multiple collisions in a single bunch crossing corresponding to the luminosity profile of the LHC during the data taking periods of interest.

Background estimate
The main processes that can mimic the signature of a leptoquark signal in the lljj channels are Z/γ * +jets, tt, VV+jets, W+jets, and QCD multijets. The Z/γ * +jets background is determined by comparing events from data and MC samples in two different regions: in the region of low (L) dilepton invariant mass around the Z boson mass (70 < M ll < 100 GeV for electrons and 80 < M ll < 100 GeV for muons) and in the region of high (H) mass M ll > 100 GeV. The low mass scaling factor R Z = N L /N MC L is measured to be 1.27 ± 0.02 for both eejj and µµjj channels, where N L and N MC L are the number of data and MC events, respectively, in the Z mass window.
The number of Z/γ * +jets events above 100 GeV is then estimated as: where N MC H is the number of MC events with M ll > 100 GeV. The estimated number of Z/γ * +jets events is obtained with the selection criteria optimized for different LQ mass hypotheses and it is used in the limit setting procedure.
The number and kinematic distributions for the tt events with two leptons of the same flavor is estimated from the number of data events that contain one electron and one muon. This type of background is expected to produce ee plus µµ and eµ final state events with equal probability. In the data the number of ee or µµ events is estimated to be: where µ and e are the muon and electron reconstruction and identification efficiencies and trig are the HLT efficiencies to select ee, µµ, and eµ events.
No QCD multijet MC events pass the µµjj final selection. A crosscheck made using a data control sample containing same sign muons confirms that the QCD multijet background is negligible in this channel.
For the first-generation LQ analyses the multijet background contribution is estimated from a data control sample as follows. The probability that an electron candidate passing loose electron requirements additionally passes all electron requirements is measured as a function of p T and η in a data sample with one and only one electron candidate, two or more jets and low E miss T . This sample is dominated by QCD multijet events and is similar in terms of jet activity to the eejj and eνjj analysis samples. A correction for a small contamination of genuine electrons passing all electron requirements is derived from MC simulations. The QCD multijet background in the final eejj (eνjj) selection is predicted by applying twice (once) the above probability to a sample with two electron candidates (one electron candidate and large E miss T ), and two or more jets, which satisfy all the requirements of the signal selections. The resulting estimate is ∼1% (∼8%) of the total background for the selections corresponding to the region of leptoquark masses where the exclusion limits are placed.
Contributions to the lljj background from VV+jets processes and single-top production are small and they are estimated using MC simulation.
In the lνjj channel, the main backgrounds come from three sources: processes that lead to the production of a genuine W bosons such as W+jets, tt, single-top production, diboson processes (WW, WZ); instrumental background, mostly caused by the misidentification of jets as leptons in multijet processes, thus creating misidentified electrons or muons and misreconstructed E miss T in the final state; and Z boson production, such as Z/γ * +jets and ZZ processes. The contribution from the principal background, W+jets, is estimated with MC simulation normalized to data at the initial selection level in the region 50 < M T < 110 GeV, where M T is the transverse mass calculated from the lepton p T and the E miss T .
In the eνjj analysis, the region 50 < M T < 110 GeV is used to determine both the W+jets and the tt normalization factors using two mutually exclusive selections (less than four jets or at least four jets with p T > 40 GeV and |η| < 2.4) that separately enhance the samples with W+jets and with tt events. The results of these two selections are used to form a system of equations: where N i , N i,W , N i,O , N i,tt , and N i,QCD are the number of events in data, W+jets, other MC backgrounds, tt-MC, and QCD multijet events obtained from data, passing selection i. The solution of the system yields the following normalization factors: R tt = 0.82 ± 0.04 (stat.)±0.02 (syst.) and R W = 1.21 ± 0.03 (stat.)±0.02 (syst.), where the systematic errors reflect the uncertainties on the QCD multijet background estimate.
In the µνjj analysis the W+jets normalization factor is calculated in the region 50 < M T < 110 GeV as: where N, N tt , N Z , N O , and N W are the number of events in data, tt, Z/γ * +jets, other backgrounds (QCD, diboson, single top production), and W+jets MC samples, normalized to the integrated luminosity of the data sample. R W is measured to be 1.21 ± 0.01 (stat.). To estimate the total tt background in the µνjj analysis, the MC prediction is compared with data in a control sample enriched with tt events (50 < M T < 110 GeV, E miss T > 55 GeV, and at least four jets). The normalization factor of the tt background, is calculated to be 0.88 ± 0.02 (stat.). The calculation of the normalization factors R W and R tt is repeated using an iterative procedure.
The contribution from QCD multijet processes after all selection criteria are applied to the µνjj sample is estimated to be negligible. No QCD MC events survive the full selection criteria optimized for any leptoquark mass hypothesis. However, as multijet processes are difficult to accurately model by MC simulation, several crosschecks are made with data control samples to ensure that the QCD multijet background in the µνjj analysis is negligible. The method used to determine the QCD multijet background in the eνjj analysis is similar to the one used for the eejj channel.

Event selection optimization
After the initial selection, the sensitivity of the search is optimized by maximizing the Gaussian signal significance S/ √ S + B in all channels. Optimized thresholds on the following variables are applied for each leptoquark mass hypothesis in the lljj channels: M ll , S ll T , and M min lj . The invariant mass of the dilepton pair, M ll , is used to remove the majority of the contribution from the Z/γ * +jets background. The variable M min lj is defined as the smaller lepton-jet invariant mass for the assignment of jets and leptons to LQs which minimizes the LQ -LQ invariant mass difference.
Thresholds on the following variables are optimized for each leptoquark mass hypothesis in the lνjj channels: E miss T , S lν T , and M lj . A minimum threshold on E miss T is used, primarily to reduce the dominant W+jets background. The variable M lj is defined as the invariant mass of the lepton-jet combination which minimizes the LQ -LQ transverse mass difference. In addition, a lower threshold is applied on the transverse mass of the lepton (electron or muon) and E miss T in the event, M T > 120 GeV.
The resulting optimized thresholds are summarized in Tables 2 and 3 for the lljj and the lνjj channels, respectively.         Data and background predictions after final selection are also shown in Figs. 6-9, which compare S T and the best combination for the lepton-jet invariant mass for a signal leptoquark mass of 600 GeV in the four decay channels considered.

Systematics Uncertainties
The uncertainty on the integrated luminosity is taken as 2.2% [29].
The statistical uncertainties on the values of R Z (R W ) after initial selection requirements are used as an estimate of the uncertainty on the normalization of the Z/γ * +jets (W+jets) backgrounds. The uncertainty on the shape of the Z/γ * +jets and W+jets distributions is calculated to be 15% (10%) and 20% (11%) for first (second) generation, respectively, by comparing the predictions of MADGRAPH samples produced with factorization or renormalization scales and matrix element-parton shower matching threshold varied up and down by a factor of two.
The uncertainty on the estimate of the tt background in the eejj and µµjj channels is derived from the statistical uncertainty of the eµjj data sample and the ratio of electron and muon       reconstruction uncertainties, which is calculated to be 2% and 3%, respectively. A 5.5% (5%) uncertainty on the normalization of the estimated tt background in the eνjj (µνjj) channel is given by the statistical uncertainty on the value of R tt after the initial selection requirements. A 10% (10%) uncertainty on the shape of the tt background distribution in the eνjj (µνjj) channel is estimated by comparing the predictions of MADGRAPH samples produced with factorization or renormalization scales and matrix element-parton shower matching thresholds varied up and down by a factor of two.
A systematic uncertainty of 50% (25%) on the QCD multijet background estimate for the eejj (eνjj) channel is estimated from the difference between the number of observed data and the background prediction in a QCD multijet-enriched data sample with lower jet multiplicity.
PDF uncertainties on the theoretical cross section of LQ production and on the final selection acceptance have been calculated using the PDF4LHC [30] prescriptions. Uncertainties on the cross section vary from 10 to 30 % for leptoquarks in the mass range of 200-900 GeV, while the effect of the PDF uncertainties on signal acceptance varies from 1 to 3%. The PDF uncertainties are not considered for background sources with uncertainties determined from data. An uncertainty on the modeling of pileup interactions in the MC simulation is determined by varying the mean of the distribution of pileup interactions by 8%.
Energy and momentum scale uncertainties are estimated by assigning a 4% uncertainty on the jet energy scale, a 1% (3%) uncertainty on the electron energy scale for the barrel (endcap) region of ECAL, and a 1% uncertainty on the muon momentum scale. The effect of electron energy, muon momentum, and jet energy resolution on expected signal and backgrounds is assessed by smearing the electron energy by 1% and 3% in the barrel and endcaps, respectively, the muon momentum by 4%, and by varying the jet energy resolution by an η-dependent value in the range 5-14%. In the lνjj analyses, the uncertainty on the energy and momentum scales and resolutions are propagated to the measurement of E miss T . The effect of these uncertainties is calculated for the (minor) background sources for which no data rescaling is applied. For the background sources for which data rescaling is applied, residual uncertainties are calculated (i.e. relative to the initial selection used to derive the rescaling factor).
Recent measurements of the muon reconstruction, identification, trigger, and isolation efficiencies using Z → µµ events show very good agreement between data and MC events [31]. A ∼ 1% discrepancy is observed in the data-to-MC comparison of the muon trigger efficiency. This discrepancy is taken as a systematic uncertainty per muon, assigned to both signal and estimated background. The electron trigger and reconstruction and identification uncertainties contribute 3% (4%) to the uncertainty in both signal and estimated background for the eejj (eνjj) channel.
The systematic uncertainties, and their effects on signal and background are summarized in Table 8 for all channels, corresponding to the final selection optimized for M LQ = 600 GeV.

Results
The number of observed events in data passing the full selection criteria is consistent with the SM background prediction in all decay channels. An upper limit on the leptoquark production cross section is therefore set using the CL S modified frequentist approach [17,18]. A log-normal probability function is used to integrate over the systematic uncertainties. Uncertainties of statistical nature are described with Γ distributions with widths determined by the number of events simulated in MC samples or observed in data control regions.   The 95% CL upper limits on σ × β 2 or σ × 2β(1 − β) as a function of leptoquark mass are shown together with the NLO predictions for the scalar leptoquark pair production cross section in Figs. 10 and 11. The theoretical cross sections are represented for different values of the renormalization and factorization scale, µ, varied between half and twice the LQ mass (blue shaded region). The PDF uncertainties are taken into account in the theoretical cross section values.
By comparing the observed upper limit with the theoretical cross section values, first-and second-generation scalar leptoquarks with masses less than 830 (640) GeV and 840 (620) GeV, respectively, are excluded with the assumption that β = 1 (0.5). This is to be compared with median expected limits of 790 (640) GeV and 780 (610) GeV.
The observed and expected limits on the branching fraction β as a function of leptoquark mass can be further improved using the combination of the lljj and lνjj channels, as shown in Figure 12. These combinations lead to the exclusion of first-and second-generation scalar leptoquarks with masses less than 640 and 650 GeV for β = 0.5, compared with median expected limits of 680 and 660 GeV.

Summary
In summary, a search for pair production of first-and second-generation scalar leptoquarks has been performed in decay channels with either two charged leptons of the same-flavour (electrons or muons) and at least two jets, or a single charged lepton (electron or muon), missing transverse energy, and at least two jets, using 7 TeV proton-proton collisions data corresponding to an integrated luminosity of 5 fb −1 . The selection criteria have been optimized for each leptoquark signal mass hypothesis. The number of observed candidates for each hypothesis agree with the estimated number of background events. The CL S modified frequentist approach has been used to set limits on the leptoquark cross section times the branching fraction for the decay of a leptoquark pair. At 95% confidence level, the pair-production of first-and second-generation leptoquarks is excluded with masses below 830 (640) GeV and 840 (650) GeV for β = 1 (0.5). These are the most stringent limits to date.   Figure 10: Left (right) frame: the expected and observed upper limits at 95% CL on the leptoquark pair production cross section times β (2β(1 − β)) as a function of the first-generation leptoquark mass obtained with the eejj (eνjj) analysis. The expected limits and uncertainty bands represent the median expected limits and the 68% and 95% confidence intervals. The systematic uncertainties reported in Table 8 are included in the calculation. The left and middle shaded regions are excluded by the current ATLAS limit [15] and CMS limits [12, 13] for β = 1(0.5) in the eejj (eνjj) channel only. The right shaded regions are excluded by these results. The σ theory curves and their bands represent, respectively, the theoretical scalar leptoquark pair production cross section and the uncertainties due to the choice of PDF and renormalization/factorization scales.  Figure 11: Left (right) frame: the expected and observed upper limits at 95% CL on the leptoquark pair production cross section times β (2β(1 − β)) as a function of the second-generation leptoquark mass obtained with the µµjj (µνjj) analysis. The expected limits and uncertainty bands represent the median expected limits and the 68% and 95% confidence intervals. The systematic uncertainties reported in Table 8 are included in the calculation. The left and middle shaded regions in the left frame are excluded by the current ATLAS limit [14] and CMS limit [11] for β = 1 in the µµjj channel only. The left shaded region in the right frame is excluded by the current ATLAS limit [14] for β = 0.5 in the µνjj channel only. The right shaded regions are excluded by these results. The σ theory curves and their bands represent, respectively, the theoretical scalar leptoquark pair production cross section and the uncertainties due to the choice of PDF and renormalization/factorization scales.  Figure 12: Left (right) frame: the expected and observed exclusion limits at 95% CL on the first-(second-)generation leptoquark hypothesis in the β versus mass plane using the central value of signal cross section for the individual eejj and eνjj (µµjj and µνjj) channels and their combination. The dark green and light yellow expected limit uncertainty bands represent the 68% and 95% confidence intervals. Solid lines represent the observed limits in each channel, and dashed lines represent the expected limits. The systematic uncertainties reported in Table 8 are included in the calculation. The shaded region is excluded by the current ATLAS limits [14,15]. financed from European Union, Regional Development Fund.