Search for top and bottom squarks from gluino pair production in final states with missing transverse energy and at least three b-jets with the ATLAS detector

This letter reports the results of a search for top and bottom squarks from gluino pair production in 4.7 fb^-1 of pp collisions at sqrt(s) = 7 TeV using the ATLAS detector at the LHC. The search is performed in events with large missing transverse momentum and at least three jets identified as originating from a b-quark. Exclusion limits are presented for a variety of gluino-mediated models with gluino masses up to 1 TeV excluded.

Eur. Phys. J. C manuscript No. (will be inserted by the editor) Search for top and bottom squarks from gluino pair production in final states with missing transverse energy and at least three b-jets with the ATLAS detector The ATLAS Collaboration July 19, 2012 Abstract This letter reports the results of a search for top and bottom squarks from gluino pair production in 4.7 fb −1 of pp collisions at √ s = 7 TeV using the ATLAS detector at the LHC. The search is performed in events with large missing transverse momentum and at least three jets identified as originating from a bquark. Exclusion limits are presented for a variety of gluino-mediated models with gluino masses up to 1 TeV excluded.
Supersymmetry (SUSY) [1][2][3][4][5][6][7][8][9] provides an extension of the Standard Model (SM) which resolves the hierarchy problem [10][11][12][13] by introducing supersymmetric partners of the known bosons and fermions. In the framework of the R-parity conserving minimal supersymmetric extension of the SM (MSSM) [14][15][16][17][18], SUSY particles are produced in pairs and the lightest supersymmetric particle (LSP) is stable, providing a possible candidate for dark matter. In a large variety of models, the LSP is the lightest neutralino (χ 0 1 ). The coloured superpartners of quarks and gluons, the squarks (q) and gluinos (g), if not too heavy, would be produced in strong interaction processes at the Large Hadron Collider (LHC) and decay via cascades ending with the LSP. The undetected LSP results in missing transverse momentum -whose magnitude is referred to as E miss T -while the rest of the cascade yields final states with multiple jets and possibly leptons. In the MSSM, the right-handed and left-handed squarks,q R andq L , can mix to form two mass eigenstatesq 1 andq 2 . The mixing effect is proportional to the masses of the SM fermion partners and can therefore be large for the third generation. This may lead to the lightest sbottom (b 1 ) and stop (t 1 ) mass eigenstates being much lighter than the other squarks. As a consequence,b 1 andt 1 could be produced with relatively large cross-sections at the LHC, either directly in pairs, or throughgg production followed byg →b 1 b org →t 1 t decays.
This letter extends the search for gluino-mediated b 1 andt 1 production at ATLAS reported earlier [19], which used 2.05 fb −1 of data collected in the first half of 2011 at a centre-of-mass energy of 7 TeV. The present analysis comprises the full 2011 dataset of 4.7 fb −1 and adopts an improved selection that requires large E miss T , no electron or muon and at least three jets identified as originating from b-quarks (b-jets) in the final state. Results are interpreted in four simplified models where sbottoms or stops are the only squarks produced in the gluino decays, leading to final states with four b-quarks. Searches in similar scenarios have also been reported by the CMS Collaboration [20].
The ATLAS detector [21] consists of inner tracking devices surrounded by a superconducting solenoid, electromagnetic and hadronic calorimeters and a muon spectrometer with a toroidal magnetic field. The inner detector provides precision tracking of charged particles for |η| < 2.5 1 . It is immersed in a 2 T magnetic field from the solenoid and consists of a silicon pixel detector, a silicon strip detector and a straw tube tracker that also provides transition radiation measurements for electron identification. The calorimeter system covers the pseudorapidity range |η| < 4.9. It is composed of sampling calorimeters with either liquid argon (LAr) or scintillating tiles as the active medium. The muon 2 The ATLAS Collaboration spectrometer has separate trigger and high-precision tracking chambers which provide muon identification and momentum measurement for |η| < 2.7.
Samples of simulated events are used for the description of the background and to model the SUSY signal. The dominant sources of background come from events with b-quarks in the final state. Monte Carlo (MC) samples of tt, W /Z and diboson events produced in association with light-and heavy-flavour jets are generated with ALPGEN [22] and the parton distribution function (PDF) set CTEQ6L1 [23]. These samples are generated with different maximum numbers of additional partons at the matrix-element level. Diboson samples are generated with up to three additional partons, tt+jet and Z(→ ℓ + ℓ − )+jets (ℓ = e,µ,τ ) samples with up to five additional partons, and W (→ ℓν)+jets and Z(→νν)+jet samples with up to six additional partons. Single top quark production is simulated with MC@NLO [24] using the next-to-leading-order (NLO) PDF set CTEQ6.6 [25]. The fragmentation and hadronisation for the ALPGEN and MC@NLO samples are performed with HERWIG [26,27], using JIMMY [28] for the underlying event. Samples of tt+W , tt+Z and tt+W W events are generated with MADGRAPH [29] interfaced to PYTHIA [30]. The signal samples are generated using Herwig++ [31]. The MC samples are processed through the ATLAS detector simulation [32] based on GEANT4 [33] taking into account the effect of multiple pp interactions per bunch crossing. For the comparison with data, all SM background cross-sections are normalised to the results of higherorder calculations when available, using the same values as in Ref. [19].
Jets are reconstructed from three-dimensional calorimeter energy clusters using the anti-k t jet algorithm [34,35] with a radius parameter of 0.4. The measured jet energy is corrected for inhomogeneities and for the noncompensating nature of the calorimeter by using p Tand η-dependent correction factors, and additional corrections for pile-up are applied [36]. Jets are required to have p T > 20 GeV, and are reconstructed in the range |η| < 4.9. Events are rejected if they include jets failing the quality criteria described in Ref. [36], or if there is any selected jet with |η| < 2 for which the scalar sum of the transverse momenta of its associated tracks is less than 5% of the jet p T . A neural-network-based algorithm [37] is used to identify jets containing a b-hadron decay. This uses as inputs the output weights of different algorithms exploiting the impact parameter of the inner detector tracks, the secondary vertex reconstruction and the topology of b-and c-hadron decays inside the jet. Three operating points are used, corresponding to efficiencies of 60%, 70% and 75% for tagging b-jets in a MC sample of tt events. In all cases the tagging rate is less than 2% for light-quark and gluon jets, 10% for τ leptons decaying hadronically and 25% for c-quark jets. The b-jets are identified within the nominal acceptance of the inner detector (|η| < 2.5) and are required to have p T > 30 GeV. To compensate for the differences between the b-tagging efficiency and the mistag rates in data and MC simulation, b-tagging scale factors are applied to each jet in the simulations, as described in Refs. [37][38][39].
Electrons are reconstructed from energy clusters in the electromagnetic calorimeter matched to a track in the inner detector. Electron candidates are required to have p T > 20 GeV and |η| < 2.47 and must satisfy the "medium" selection criteria described in Ref. [40]. Muons candidates are identified using a match between an extrapolated inner detector track and one or more track segments in the muon spectrometer, and are required to have p T > 10 GeV and |η| < 2.4.
To resolve overlaps between reconstructed jets and leptons, jets within a distance of ∆R = 0.2 of an electron candidate are rejected. Furthermore, any lepton candidate with a distance ∆R < 0.4 to the closest remaining jet is discarded. Events containing any remaining electrons and muons are vetoed in the control and signal regions defined in Tables 1 and 2.
The measurement of the missing transverse momentum two-dimensional vector (and its magnitude E miss T ) is based on the transverse momenta of all remaining jets with |η| < 4.9, all electron and muon candidates and all calorimeter clusters not associated to such objects.
Events are selected using triggers requiring one high p T jet and E miss T . Different trigger thresholds were used to cope with the increasing luminosity. These triggers are fully efficient for this analysis, which requires one jet with p T > 130 GeV and E miss T > 160 GeV at the offline reconstruction stage. Events must pass basic quality criteria to reject detector noise and non-collision backgrounds. They are also required to have a reconstructed primary vertex associated with five or more tracks with p T > 0.4 GeV; when more than one such vertex is found, the vertex with the largest summed p 2 T of the associated tracks is chosen as the hard interaction vertex. Events are required to have at least three b-tagged jets, and two jet-multiplicity regions (N J ≥ 4 and N J ≥ 6) are considered by selecting jets with |η| < 2.8 and p T > 50 GeV.
Two variables are calculated from the reconstructed objects to further select the events: m eff and ∆φ min . The effective mass m eff is defined as the scalar sum of the E miss T and the p T of all selected jets in a given jetmultiplicity region. The ∆φ min is defined as the minimum azimuthal separation between the selected jets and the missing transverse momentum direction. Plac-Search for top and bottom squarks in final states with missing transverse energy and at least three b-jets 3 Common criteria: lepton veto, p j 1 T > 130 GeV, Table 1 Definition of the five signal regions based on the number of jets (N J ), the E miss T , m eff requirements and the b-tagging operating point (OP).
ing the requirements ∆φ min > 0.4 and E miss T /m eff > 0.2 reduces the amount of multi-jet background, where E miss T results from mis-reconstructed jets or from neutrinos emitted close to the direction of the jet axis.
Two sets of signal regions are defined which yield good signal sensitivity for the various models and parameter values studied here. They are characterised by having at least four (SR4) or six (SR6) jet candidates, no electron or muon, and are further classified as loose (L), medium (M) or tight (T) depending on the E miss T and m eff thresholds and on the b-tagging operating point. The requirements that characterise each signal region are summarised in Table 1.
The main source of reducible background is the production of tt events in association with additional jets followed by the leptonic decay of one W boson, where the lepton is not reconstructed or is misidentified as a jet (mainly through the hadronic decays of a τ lepton). This background is estimated by normalising the MC event yield in the signal region to the extrapolated event yield observed in a tt-dominated control region. Systematic uncertainties that are correlated between the control and the signal regions largely cancel out in this procedure. Additional sources of reducible background are single top, tt+W /Z and W /Z+heavy-flavour jets. Their contributions are taken from MC simulation and account for 10% to 20% of the total background depending on the signal region. The irreducible background tt+bb is also estimated from MC simulation and accounts for about 10% of the total background in all signal regions. The reducible contribution from multi-jet events is estimated with a data-driven method, based on a jet response smearing technique [41], and is found to account for less than 5% of the total background in all signal regions.
Four control regions where the tt+jets background accounts for more than 80% of the total yield are defined by applying the same jet requirements and lepton veto as in the signal regions, but requiring exactly two b-jets instead of three or more. The requirements m eff > 500 GeV and E miss T > 160 GeV are applied to all control regions to make them kinematically similar Common criteria: lepton veto, p j 1 T > 130 GeV,  Table 3 Expected numbers of SM events and observed data events in the four tt control regions. The contribution from tt+jets events is taken directly from MC simulation. The column "others" includes the contributions from single top, tt+bb, tt+W /Z and W /Z+jets processes, also estimated from MC simulation, and the contribution from multi-jet events which is estimated with the jet smearing technique and accounts for less than 5% of the total background. The column "SM" shows the total expected background and is the sum of the columns "tt+jets" and "others". The uncertainties include all detector-related systematic uncertainties.
The dominant detector-related systematic effects are due to the jet energy scale (JES) and resolution (JER) uncertainties, and the uncertainty on the b-tagging efficiency and mistag rates. The JES uncertainty is derived from a combination of simulations, test beam data and in-situ measurements [36], and includes additional uncertainties due to the jet flavour and nearby jets. Uncertainties on the JER are obtained with an in-situ measurement of the jet response asymmetry in di-jet events. These uncertainties on jets are propagated to the E miss T measurement, and additional uncertainties on E miss T arising from energy deposits not associated with any reconstructed objects are also included. The b-tagging uncertainty is evaluated by varying the η-, p Tand flavour-dependent scale factors applied to each jet in the simulation within a range that reflects the systematic uncertainty on the measured tagging efficiency and mistag rates. The systematic uncertainties in the modelling of the tt+jets background are assessed as follows: the uncertainty due to the choice of the MC generator is estimated by comparing the leading-order ALPGEN generator to the MC@NLO generator; the uncertainty due to the factorisation and matching scale ambiguities in ALPGEN are estimated by independently varying their nominal settings by factors of one half and two; the uncertainty due to the finite number of additional partons at the matrix-element level is assessed by comparing inclusive ALPGEN samples generated with up to three and up to five extra partons. Finally the PDF uncertainties are estimated using the MSTW2008NNLO PDF set. Uncertainties of 100% are assumed for the multi-jet prediction and for the cross-section of tt and W /Z production in association with a bb pair. For tt+W /Z production, an uncertainty of approximatively 70% has been derived from the variations of the factorisation and renormalisation scales and from the PDF uncertainties [42].
The tt+jets yield in each signal region is extrapolated from the measured number of events in the corresponding control region (as per Table 2) using a fit based on the profile likelihood method [43]. Each pair of control and signal regions is fitted separately, assuming no signal events. The free parameter in each fit is the tt+jets overall normalisation scale, while the contributions from subdominant background processes are fixed at the expected value. Systematic uncertainties are treated as nuisance parameters constrained with a Gaussian function and correlations are taken into account where appropriate. The results of the fits and the numbers of observed events for each signal region are summarised in Table 4. The fitted values of the normalisation factors for tt+jets are compatible with one and the main impact of the data-driven estimate is a reduction in the uncertainty by approximately a factor of two. Figure 1 shows the measured m eff distributions and the MC predictions for the SM backgrounds in each signal region. Also shown are the prediction of two benchmark signal models described below.
The reliability of the MC extrapolation of the tt background to larger b-jet multiplicities has been checked in validation regions defined with kinematic cuts similar to those used in the control and signal regions, except that exactly one isolated electron or muon is required. The transverse mass of the lepton and the E miss T is required to be less than 100 GeV in all validation regions to minimise the possible contamination from stop production. The extrapolated event yield in the validation regions with at least three b-jets from the validation regions with exactly two b-jets is found to be consistent with the number of observed events for all b-tagging operating points.  Table 4 Comparison between the results of the fits and the numbers of observed events in the five signal regions. The tt+jets event yield predicted by the MC simulation is quoted in parentheses. The column "others" includes the contributions from single top, tt+bb, tt+W /Z, W /Z+jets and multijet processes. Multi-jet events contribute less than 5% of the total background. The column "SM" shows the total expected background and is the sum of the columns "tt+jets" and "others". The uncertainties include the statistical and systematic uncertainties.
The background predictions have been further validated using a data-driven method that simultaneously estimates all SM background contributions with at least one misidentified b-jet. This method consists of predicting the number of jets originating from b-quarks in each event by solving a system of equations based on the number of b-tagged and non b-tagged jets in the event, along with the b-tagging efficiency and mistag rates. Consistent background predictions with respect to the fit results have been found in all signal regions.
Limits for non-SM signal at 95% confidence level (CL) are derived by testing the signal plus background hypothesis in each signal region with the CL s prescription [43]. These limits are obtained with fits similar to those used to estimate the background in each signal region, except that the number of observed events in the signal region is added as an input to the fit and a second free parameter for the non-SM signal strength, constrained to be non-negative, is adjusted in the likelihood maximisation. This additional free parameter ensures a proper treatment of the expected signal contamination in the control regions when the results are interpreted in the framework of specific SUSY scenarios. Modelindependent upper limits at 95% CL on the number of signal events and on the visible cross-section (defined as the cross-section times kinematic acceptance times experimental efficiency) for non-SM contributions derived for each signal region are given in Table 5.
These data have been used to derive limits in the parameter space of the following SUSY models: Gluino-sbottom model: MSSM scenarios where theb 1 is the lightest squark, all other squarks are heavier than the gluino, and mg > mb branching ratio forg →b 1 b decays is 100%. Sbottoms are produced viagg or byb 1b1 direct pair production and are assumed to decay exclusively viab 1 → bχ 0 1 , where mχ 0 1 is set to 60 GeV. Exclusion limits are presented in the (mg, mb 1 ) plane.
Gbb model: Simplified scenarios, whereb 1 is the lightest squark but mg < mb 1 . Pair production of gluinos is the only process taken into account since the masses of all other sparticles apart from theχ 0 1 are set above the TeV scale. A three-body decay via an off-shell sbottom is assumed for the gluino, yielding a 100% BR for the decayg → bbχ 0 1 . The sbottom mass has no impact on the kinematics of the decay and the exclusion limits are presented in the (mg, mχ 0 1 ) plane. Gluino-stop model: MSSM scenarios where thẽ t 1 is the lightest squark, all other squarks are heavier than the gluino, and mg > mt 1 + m t , so the branching ratio forg →t 1 t decays is 100%. Stops are produced viagg andt 1t1 and are assumed to decay exclusively viat 1 → bχ ± 1 . The neutralino mass is set to 60 GeV, the chargino mass to 120 GeV and the latter is assumed to decay through a virtual W boson. Exclusion limits are presented in the (mg, mt 1 ) plane.
Gtt model: Simplified scenarios, wheret 1 is the lightest squark but mg < mt 1 . Pair production of gluinos is the only process taken into account since the mass of all other sparticles apart from theχ 0 1 are above the TeV scale. A three-body decay via off-shell stop is assumed for the gluino, yielding a 100% BR for the decaỹ g → ttχ 0 1 . The stop mass has no impact on the kinemat- The ATLAS Collaboration    [44,45], D0 [46] and ATLAS [19,42,47,48] analyses.
ics of the decay and the exclusion limits are presented in the (mg, mχ 0 1 ) plane.
The SR4 regions are mostly sensitive to the SUSY models where sbottom production dominates, whilst the SR6 regions are used to set exclusion limits in models characterised by on-shell or off-shell stop production, where top-enriched final states are expected. The signal region with the best expected sensitivity at each point in the parameter space is used to derive the limits at 95% CL. Signal cross-sections are calculated to nextto-leading order in the strong coupling constant, including the resummation of soft gluon emission at next-toleading-logarithmic accuracy (NLO+NLL) [49][50][51][52][53]. The nominal cross-section and the uncertainty σ SUSY Theory are taken from an envelope of cross-section predictions using different PDF sets and factorisation and renormalisation scales, as described in Ref. [54]. All detectorrelated systematic uncertainties are treated as fully correlated between signal and backgrounds. In the Gbb scenario, the impact of initial-state radiation (ISR) is expected to be large in the region with low mg − mχ0 1 due to the small signal acceptance. Therefore, an uncertainty on the modelling of ISR is assessed by comparing the signal acceptance obtained with the Herwig++ samples to the one obtained with dedicated MADGRAPH samples generated with additional jets. This uncertainty varies from 4% to 35% as a function of mg − mχ0 1 and is included in the ±1σ SUSY Theory band.
The expected and observed 95% CL exclusion limits in the four models considered above are shown in Figure 2. In the gluino-sbottom model, gluino masses below 1000 GeV are excluded for sbottom masses up to about 870 GeV using the most conservative −1σ SUSY Theory hypothesis. This extends by approximatively 100 GeV the limits derived in the same scenario by the previous ATLAS analysis performed with 2 fb −1 [19] and is complementary to the ATLAS search for direct sbottom pair production, also carried out with 2 fb −1 [47]. The exclusion is less stringent in the region with low mg − mb 1 , where softer jets are expected. Because of the kinematic cuts applied, the limits depend on the neutralino mass assumption for low mass splitting be-tween the sbottom and the neutralino as shown for the Gbb model where gluino masses below 1020 GeV are excluded for neutralino masses up to about 400 GeV, improving the previous ATLAS limits [19] by approximatively 100 GeV. In the gluino-stop model, gluino masses below 820 GeV are excluded for stop masses up to 640 GeV, extending the previous ATLAS limits [19,42] by approximatively 150 GeV. In the Gtt model, gluino masses below 940 GeV are excluded for mχ0 1 < 50 GeV while neutralino masses below 320 GeV are excluded for mg = 800 GeV. This search extends the exclusion limits on the gluino mass from the AT-LAS multi-jet analysis carried out with the same data set [48] and from the CMS same-sign dilepton analysis performed with 5 fb −1 [20] by approximatively 60 GeV and 130 GeV, respectively, for neutralino masses below 100 GeV. In the region with low mg − mχ0 1 , the limits obtained with the CMS analysis are most stringent due to the softer kinematic cuts.
In summary, this letter presents results from a search for top and bottom squarks in the decay of gluino pairs produced in pp collisions at √ s = 7 TeV, based on 4.7 fb −1 of ATLAS data. The events are selected with large E miss T , four or six jets and at least three jets originating from b-quarks in the final state. The results are in agreement with the SM background prediction and translate into 95% CL upper limits on excluded masses for a variety of SUSY benchmark scenarios. Gluino masses up to 1 TeV are excluded, depending on the model, which significantly extends the previous results.
We thank CERN for the very successful operation of the LHC, as well as the support staff from our institutions without whom ATLAS could not be oper-