Searches for pair production of third-generation squarks in \(\sqrt{s}=13\) \(\,\text {TeV}\) pp collisions

The European Physical Journal C, May 2017

Searches are presented for direct production of top or bottom squark pairs in proton–proton collisions at the CERN LHC. Two searches, based on complementary techniques, are performed in all-jet final states that are characterized by a significant imbalance in transverse momentum. An additional search requires the presence of a charged lepton isolated from other activity in the event. The data were collected in 2015 at a centre-of-mass energy of 13\(\,\text {TeV}\) with the CMS detector and correspond to an integrated luminosity of 2.3\(\,\text {fb}^{-1}\). No statistically significant excess of events is found beyond the expected contribution from standard model processes. Exclusion limits are set in the context of simplified models of top or bottom squark pair production. Models with top and bottom squark masses up to 830 and 890\(\,\text {GeV}\), respectively, are probed for light neutralinos. For models with top squark masses of 675\(\,\text {GeV}\), neutralino masses up to 260\(\,\text {GeV}\) are excluded at 95% confidence level.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

https://link.springer.com/content/pdf/10.1140%2Fepjc%2Fs10052-017-4853-2.pdf

Searches for pair production of third-generation squarks in \(\sqrt{s}=13\) \(\,\text {TeV}\) pp collisions

Eur. Phys. J. C Searches for pair production of third-generation squarks in √ s = 13 TeV pp collisions CMS Collaboration 0 1 2 3 4 5 8 9 11 12 13 0 CERN , 1211 Geneva 23 , Switzerland 1 Department of Physics, University of Helsinki , Helsinki , Finland P. Eerola, J. Pekkanen, M. Voutilainen 2 University of Split, Faculty of Science , Split , Croatia Z. Antunovic, M. Kovac 3 State Key Laboratory of Nuclear Physics and Technology, Peking University , Beijing, China Y. Ban, G. Chen, Q. Li, S. Liu, Y. Mao, S. J. Qian, D. Wang, Z. Xu 4 University of Hamburg , Hamburg, Germany V. Blobel, M. Centis Vignali, A. R. Draeger, T. Dreyer, E. Garutti, D. Gonzalez, J. Haller, M. Hoffmann, A. Junkes, R. Klanner, R. Kogler, N. Kovalchuk, T. Lapsien, T. Lenz, I. Marchesini, D. Marconi, M. Meyer, M. Niedziela, D. Nowatschin, F. Pantaleo 5 RWTH Aachen University, I. Physikalisches Institut , Aachen, Germany C. Autermann, S. Beranek, L. Feld, M. K. Kiesel, K. Klein, M. Lipinski, M. Preuten, C. Schomakers, J. Schulz, T. Verlage, V. Zhukov 6 , T. Peiffer , A. Perieanu, J. Poehlsen, C. Sander, C. Scharf, P. Schleper, A. Schmidt, S. Schumann, J. Schwandt, H. Stadie, G. Steinbrück, F. M. Stober, M. Stöver, H. Tholen, D. Troendle, E. Usai, L. Vanelderen, A. Vanhoefer, B. Vormwald 7 , P. Raics, Z. L. Trocsanyi, B. Ujvari 8 Institute of Physics, University of Debrecen , Debrecen , Hungary M. Bartók 9 Institute of Experimental Physics, Faculty of Physics, University of Warsaw , Warsaw , Poland K. Bunkowski, A. Byszuk 10 , K. Doroba , A. Kalinowski, M. Konecki, J. Krolikowski, M. Misiura, M. Olszewski, M. Walczak 11 Istanbul Technical University , Istanbul , Turkey A. Cakir, K. Cankocak, S. Sen 12 Bogazici University , Istanbul , Turkey E. Gülmez, M. Kaya 13 Middle East Technical University, Physics Department , Ankara , Turkey B. Bilin, S. Bilmis, B. Isildak 14 , M. Yalvac , M. Zeyrek Searches are presented for direct production of top or bottom squark pairs in proton-proton collisions at the CERN LHC. Two searches, based on complementary techniques, are performed in all-jet final states that are characterized by a significant imbalance in transverse momentum. An additional search requires the presence of a charged lepton isolated from other activity in the event. The data were collected in 2015 at a centre-of-mass energy of 13 TeV with the CMS detector and correspond to an integrated luminosity of 2.3 fb−1. No statistically significant excess of events is found beyond the expected contribution from standard model processes. Exclusion limits are set in the context of simplified models of top or bottom squark pair production. Models with top and bottom squark masses up to 830 and 890 GeV, respectively, are probed for light neutralinos. For models with top squark masses of 675 GeV, neutralino masses up to 260 GeV are excluded at 95% confidence level. - The standard model (SM) has been extremely successful at describing particle physics phenomena. Nevertheless, it suffers from shortcomings such as the hierarchy problem [1– 6], the need for fine-tuned cancellations of large quantum corrections to keep the Higgs boson mass near the electroweak scale. Supersymmetry (SUSY), based on a symmetry between bosons and fermions, is an attractive extension of the SM. A key feature of SUSY is the existence of a superpartner for every SM particle with the same quantum numbers, except for spin, which differs by one half unit. In R-parity conserving SUSY models [7,8], supersymmetric particles are created in pairs, and the lightest supersymmetric particle (LSP) is stable [9,10] and considered to be a candidate for dark matter [11]. Supersymmetry can potentially provide a “natural”, i.e. not fine-tuned, solution to the hierarchy problem through the cancellation of quadratic divergences in particle and sparticle loop corrections to the Higgs boson mass. In natural SUSY models light top and bottom squarks with masses close to the electroweak scale are preferred. This paper presents three complementary searches for direct production of a pair of top (t1t1) or bottom squarks (b1b1), where the subscript here denotes the less massive partner of the corresponding SM fermion’s chirality states. The first search targets top squark pair production in the alljet final state, while the second focuses on the single-lepton final state. These two analyses were explicitly designed for complementarity, allowing for a combination of the results to enhance the sensitivity. The third search targets bottom squark pair production in the all-jet final state. The searches are performed using the data collected in proton–proton collisions at a centre-of-mass energy of 13 TeV with the CMS detector at the CERN LHC in 2015, corresponding to an integrated luminosity of 2.3 fb−1. The results of similar searches were previously reported by the ATLAS and CMS collaborations using proton–proton collisions at 7 and 8 TeV [12–25] and by the CDF and D0 collaborations in pp¯ collisions at 1.96 TeV at the Fermilab Tevatron [26–30]. With the increase in LHC collision energy from 8 to 13 TeV, the cross section to produce signal events is enhanced by a factor of 8–12 for a top or bottom squark mass in the range 700–1000 GeV [31,32]. Therefore, new territory can be explored even with the relatively small amount of data collected in 2015. The CMS and ATLAS collaborations have already provided first exclusion results for these models in the all-jet and single-lepton final states [33–36]. Unlike the more generic searches for new phenomena presented by the CMS collaboration in Refs. [33– 35], the searches described in this paper directly target top and bottom squark production through the design of search regions that exploit the specific characteristics of these signal models, for instance through the use of a top quark tagging algorithm in the top squark search in the all-jet final state to identify boosted hadronically decaying top quarks originating from top squark decays. The decay modes of top squarks depend on the sparticle mass spectrum. Figure 1 illustrates the top and bottom squark W− ¯ ~ t1 ¯ ~ t1 ~ ¯ b1 ~ ¯ t1 ~¯ t1 Fig. 1 Feynman diagrams for pair production of top and bottom squarks via the decay modes considered in this paper. The model with 50% branching fractions for t1 → t(∗)χ10 and t1 → bχ1± → bW±∗χ10 decays leads to the final states in a–c decay modes explored in this paper. The simplest top squark decay modes are t1 → t(∗)χ10 and t1 → bχ1± → bW±(∗)χ10, with χ1± representing the lightest chargino, and with intermediate particles that can be virtual marked by asterisks. In these decay modes, the neutralino and charginos are mixtures of the superpartners of electroweak gauge and Higgs bosons, and χ10 is considered to be an LSP that escapes detection, leading to a potentially large transverse momentum imbalance in the detector. The two analyses of top squark pair production in the all-jet and single-lepton final states probe both of these t1 decay modes. In the t1 → t(∗)χ10 decay mode, the top quark is produced off-shell when m ≡ mt1 − mχ10 < mt, while in the t1 → bχ1± decay mode, the experimental signature is affected by the mass of the chargino. We consider a model in which both top squarks decay via the t1 → t(∗)χ 0 1 decay mode. A second model in which the branching fraction for each of the two top squark decay modes is 50% is also considered, under the assumption of a compressed mass spectrum in which the mass of χ1± is only 5 GeV greater than that of χ 0, with the W bosons resulting from 1 chargino decays consequently being produced heavily offshell. If m < mW, t1 can decay through a four-body decay involving an SM fermion pair ff as t1 → bffχ10, or through a flavour changing neutral current decay t1 → cχ10. The analysis of bottom squark pair production considers the decay mode b1 → bχ10 within the allowed phase space, and also probes top squark pair production in the t1 → cχ10 decay scenario. This paper is organized as follows. Section 2 contains a brief description of the CMS detector, while Sect. 3 discusses the event reconstruction and simulation. Sections 4, 5, and 6 present details for the all-jet top squark search, the singlelepton top squark search, and the all-jet bottom squark search, respectively. Section 7 describes the systematic uncertainties affecting the results of the three analyses. The interpretation of the results in the form of exclusion limits on models of top or bottom squark pair production is discussed in Sect. 8, followed by a summary in Sect. 9. 2 The CMS detector The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diameter, providing a magnetic field of 3.8 T. Within the solenoid volume are an all-silicon pixel and strip tracker, a lead tungstate crystal electromagnetic calorimeter (ECAL), and a brass and scintillator hadron calorimeter (HCAL), each composed of a barrel and two endcap sections. Forward calorimeters extend the pseudorapidity (η) coverage provided by the barrel and endcap detectors. Muons are measured in gas-ionization detectors embedded in the steel flux-return yoke outside the solenoid. The first level of the CMS trigger system, composed of custom hardware processors, uses information from the calorimeters and muon detectors to select the most interesting events in a fixed time interval of less than 4 µs. The high-level trigger processor farm further decreases the event rate from around 100 kHz to around 1 kHz, before data storage. A more detailed description of the CMS detector, together with a definition of the coordinate system used and the relevant kinematic variables, can be found in Ref. [37]. 3 Reconstruction algorithms and simulation Event reconstruction uses the particle-flow (PF) algorithm [38,39], combining information from the tracker, calorimeter, and muon systems to identify charged hadrons, neutral hadrons, photons, electrons, and muons in an event. The missing transverse momentum, pTmiss, is computed as the negative vector sum of the transverse momenta ( pT) of all PF candidates reconstructed in an event, and its magnitude E miss is an T important discriminator between signal and SM background. Events selected for the searches are required to pass filters designed to remove detector- and beam-related noise and must have at least one reconstructed vertex. Usually more than one such vertex is reconstructed, due to pileup, i.e. multiple pp collisions within the same or neighbouring bunch crossings. The reconstructed vertex with the largest pT2 of associated tracks is designated as the primary vertex. Charged particles originating from the primary vertex, photons, and neutral hadrons are clustered into jets using the anti-kT algorithm [40] implemented in FastJet [41] with a distance parameter of 0.4. The jet energy is corrected to account for the contribution of additional pileup interactions in an event and to compensate for variations in detector response [41,42]. Jets considered in the searches are required to have their axes within the tracker volume, within the range |η| < 2.4. Jets originating from b quarks are identified with the combined secondary vertex (CSV) algorithm [43,44] using two different working points, referred to as “loose” and “medium”. The b tagging efficiency for jets originating from b quarks is about 80 and 60% for the loose and medium working point, respectively, while the misidentification rates for jets from charm quarks, and from light quarks or gluons are about 45 and 12%, and 10 and 2%, respectively. The “CMS top (quark) tagging” (CTT) algorithm [45– 47] is used to identify highly energetic top quarks decaying to jets with the help of observables related to jet substructure [48,49] and mass. For a relativistic top quark with a Lorentz boost γ = E /m, the W boson and b quark produced in the top quark decay are expected to be separated √ by a distance R ≡ ( η)2 + ( φ)2 ≈ 2/γ (where φ is the azimuthal angle in radians). In cases where the W boson subsequently decays hadronically, the three resulting jets from the W boson decay and the hadronization of the b quark are likely to be merged into a single jet by a clustering algorithm with a distance parameter larger than 2/γ . To identify hadronically decaying top quarks with pT > 400 GeV, we therefore use jets reconstructed using the anti-kT algorithm with a distance parameter of 0.8 to try to cluster the top quark decay products into a single jet. The next step of top quark reconstruction is an attempt to decompose the candidate jet into at least three subjets with the help of the CambridgeAachen jet clustering algorithm [50,51], the invariant mass of which is required to be consistent with the top quark mass (140–250 GeV). The final requirement of top quark identification is that the minimum invariant mass of any pair of the three subjets with the highest pT must exceed 50 GeV. The efficiency of the CTT algorithm to identify jets originating from top quark decays is measured to be about 30–40% while the misidentification rate is found to be about 4–6%, depending on the pT of the top quark candidates. No disambiguation is performed between top quark candidates and jets reconstructed with a distance parameter of 0.4. Electron candidates are reconstructed by first matching clusters of energy deposited in the ECAL to reconstructed tracks. Selection criteria based on the distribution of the shower shape, track–cluster matching, and consistency between the cluster energy and track momentum are then used in the identification of electron candidates [52]. Muon candidates are reconstructed by requiring consistent hit patterns in the tracker and muon systems [53]. Electron and muon candidates are required to be consistent with originating from the primary vertex by imposing restrictions on the size of their impact parameters in the transverse plane and longitudinal direction with respect to the beam axis. The relative isolation variable Irel for these candidates is defined as the scalar sum of the transverse momenta of all PF candidates, excluding the lepton, within a pT-dependent cone size of radius R around the candidate’s trajectory, divided by the pT ≤ 50 GeV, R = ⎪ 10 GeV/ pT, 50 < pT < 200 GeV, ⎩ pT ≥ 200 GeV. The shrinking cone radius for higher- pT leptons allows us to maintain high efficiency for the collimated decay products of boosted heavy objects. The isolation sum is corrected for contributions originating from pileup interactions through an area-based estimate [42] of the pileup energy deposited in the cone. Hadronically decaying τ lepton (τh) candidates are reconstructed using the CMS hadron-plus-strips (HPS) algorithm [54]. The constituents of the reconstructed jets are used to identify individual τ lepton decay modes with one charged hadron and up to two neutral pions, or three charged hadrons. The presence of extra particles within the jet, not compatible with the reconstructed decay mode, is used as a criterion to discriminate τh decays from other jets. Photon candidates are reconstructed from energy deposited in the ECAL, and selected using the distribution of the shower shape variable, the photon isolation, and the amount of leakage of the photon shower into the HCAL [55]. Monte Carlo (MC) simulations of events are used to study the properties of SM backgrounds and signal models. The MadGraph 5_amc@nlo 2.2.2 generator [56] is used in leading-order (LO) mode to simulate events originating from tt, W+jets, Z+jets, γ +jets, and quantum chromodynamics multijet processes (QCD), as well as signal events, based on LO NNPDF3.0 [57] parton distribution functions (PDFs). Single top quark events produced in the t W channel and tt samples used in the single-lepton analysis are generated at next-to-leading order (NLO) with Powheg v2 [58–61], while rare SM processes such as ttZ and ttW are generated at NLO using the MadGraph 5_amc@nlo 2.2.2 program, using NLO NNPDF3.0 PDFs. Parton showering and hadronization is generated using Pythia8.205 [62]. The response of the CMS detector for the SM backgrounds is simulated via the Geant4 [63] package. The CMS fast simulation package [64] is used to simulate all signal samples, and is verified to provide results that are consistent with those obtained from the full Geant4-based simulation. Event reconstruction is performed in the same manner as for collision data. A nominal distribution of pileup interactions is used when producing the simulated samples. The samples are then reweighted to match the pileup profile observed in the collected data. The signal production cross sections are calculated using NLO with next-to-leading logarithm (NLL) soft-gluon resummation calculations [31]. The most precise cross section calculations are used to normalize the SM simulated samples, corresponding most often to next-to-next-to-leading order (NNLO) accuracy. 4 Search for top squarks in the fully-hadronic final state The top squark search in the all-jet final state is characterized by the categorization of events into exclusive search regions based on selection criteria applied to kinematic variables related to jets and E miss, and the use of the CTT algo T rithm to identify boosted top quark candidates. The main backgrounds in the search regions are estimated from dedicated data control samples. 4.1 Analysis strategy The events in this analysis are recorded using a trigger that requires the presence of two or more energetic jets within the tracker acceptance and large E miss. To be efficient, events T selected offline are therefore required to have at least two jets with pT > 75 GeV, |η| < 2.4, and E miss > 250 GeV. In T order to reduce SM backgrounds with intrinsic E miss such as T leptonic tt and W+jets events, we reject events with isolated electrons or muons that have pT > 5 GeV, |η| < 2.4, and Irel less than 0.1 or 0.2, respectively. The contribution from events in which a W boson decays to a τ lepton is reduced by rejecting events containing isolated charged-hadron PF candidates with pT > 10 GeV and |η| < 2.5 that are consistent with τh decays. The isolation requirement applied is based on a discriminant obtained from a multivariate boosted decision tree (BDT) trained to distinguish the characteristics of charged hadrons originating from τh decays. The transverse mass MT of the system comprising the charged-hadron PF candidate and pmiss is required to be less than 100 GeV assur T ing consistency with τh originating from a W boson decay, to minimize loss of signal at high jet multiplicity. The transverse mass for a particle q (in this case, the τh candidate) is defined as: MT(q, pTmiss) ≡ with qT denoting the particle transverse momentum, and φ the azimuthal separation between the particle and pmiss. T Events selected for the search sample must also have at least five jets with pT > 20 GeV, at least two of which must be b-tagged satisfying the loose working point of the CSV algorithm, with one or more of them required to additionally satisfy the medium working point. In addition, the absolute value of the azimuthal angle between pmiss and the closest of T the four highest- pT (leading) jets, φ1234, must be at least 0.5. An imbalance in event pT is produced in QCD events through a mismeasurement of jet pT, and is often aligned with one of the leading jets in the event. The requirement on φ1234 therefore greatly reduces the contribution of the QCD background. The set of selection criteria defined above will be referred to as the “baseline selection” for this search. After imposing the baseline selection, we subdivide the event sample into categories based on kinematic observables related to jets and E miss to improve the power of the anal T ysis to discriminate between signal and the remaining SM background. The dominant sources of SM background are tt, W+jets, and Z+jets events. The contribution from tt and W+jets processes arises from events with W bosons decaying leptonically, in which the charged lepton either falls outside of the kinematic acceptance, or, in most cases, evades identification, and may be reconstructed as a jet. Large E miss T can be generated by the associated neutrino, allowing such events to satisfy the baseline selection criteria. This background is collectively referred to as the “lost-lepton background”. Contributions arising from ttW and single top quark processes also enter this category, but with lesser importance. The contributions from Z+jets and ttZ events arise when the Z boson decays to neutrinos, producing thereby a significant amount of E miss. The QCD background is reduced to a sub T dominant level by the requirements of the baseline selection. In tt events with a lost lepton, the transverse mass of the b quark pmiss system arising from the same top quark decay T as the lost lepton has a kinematic endpoint at the mass of the top quark. The observable MT(b1,2, pTmiss) is defined as MT(b1,2, pTmiss) ≡ min[MT(b1, pTmiss), MT(b2, pTmiss)], (3) where b1, b2 are the two selected b-tagged jets with highest values in the CSV discriminant. Imposing a minimum requirement of 175 GeV on MT(b1,2, pmiss) reduces a sigT nificant portion of the tt background, but also results in a loss in signal efficiency for models with small m, as seen in Fig. 2, in which signal models with different top squark and neutralino mass hypotheses are shown, with the first number indicating the assumed top squark mass in units of GeV and the second the neutralino mass. To benefit from the separation power provided by this variable, we define two search categories, one with MT(b1,2, pmiss) ≥ 175 GeV, T taking advantage of the corresponding reduction in tt background for signal models with large m, and another with MT(b1,2, pmiss) < 175 GeV to retain the statistical power of T events populating the low-MT(b1,2, pTmiss) region for models with small m. Signal events with all-jet top quark decays should have at least six jets in the final state, although in the case of signals with compressed mass spectra these jets can be too soft in pT to satisfy the jet selection threshold. Additional jets may be produced through initial-state radiation (ISR). The jet multiplicity is lower for the semileptonic tt background, as well as for the other backgrounds remaining after the baseline selection. A requirement of higher reconstructed jet multiplicity therefore improves the discrimination of signal events from the SM background. We consider two regions in jet multiplicity for the analysis, a high-Nj region (≥7 jets) that benefits Fig. 2 The MT(b1,2, pmiss) distribution after the baseline selection of T the top squark search in the all-jet final state (top), and the number of reconstructed top quarks for events in the high-MT(b1,2, pmiss) category T (bottom). Signal models with different top squark and neutralino mass hypotheses are shown, with the first number indicating the assumed top squark mass in units of GeV and the second the neutralino mass. The expected signal yields are scaled up by a factor of 10 to facilitate comparison of the distributions with expectations from SM backgrounds. In this and subsequent figures, the last bin shown includes the overflow events from this improved discrimination, and a medium-Nj region (5–6 jets) to preserve signal events with fewer reconstructed jets. The high-Nj region in conjunction with the low threshold on the pT of selected jets improves sensitivity for signal models with soft decay products in the final state. 2.3 fb-1 (13 TeV) 2.3 fb-1 (13 TeV) 50 100 150 200 250 300 350 400 450 500 In the high-MT(b1,2, pmiss) category, requiring the pres T ence of at least one top quark reconstructed by the CTT algorithm (Nt ≥ 1) ensures a high-purity selection of signal events with highly boosted top quarks, at the sacrifice of some loss in signal efficiency. To benefit from this highpurity region, without giving up signal events that would enter the Nt = 0 region, we use both regions to extract the final signal. Figure 2 shows the Nt distribution for events in the high-MT(b1,2, pmiss) category. Subdividing each Nt T region by the number of b-tagged jets (Nb) that satisfy the medium working point of the CSV algorithm provides even greater discrimination of signal from background. Since there are relatively few events in the Nt ≥ 1 category, the subcategorization in Nj is not performed for these events because it provides no additional gain after the Nb subdivision. The event categorization according to MT(b1,2, pTmiss), Nj, Nb, and Nt is summarized in Table 1. In each of these categories, we use E miss as the final discriminant to characterize T and distinguish potential signal from the SM background by defining five E miss regions. The analysis is therefore carried T out in a total of 50 disjoint search regions (SRs). 4.2 Background estimation 4.2.1 Estimation of the lost-lepton background The lost-lepton background is estimated from a single-lepton control sample, selected using the same trigger as the search sample, and consisting of events that have at least one lepton ( ) obtained by inverting the electron and muon rejection criteria. Studies in simulation indicate that the event kinematics for different lepton flavours are similar enough to estimate them collectively from the same control sample. Potential signal contamination is suppressed by requiring MT( , pTmiss) < 100 GeV. If there is more than one lepton satisfying the selection criteria, the lepton used to determine MT( , pTmiss) is chosen randomly. The events selected in the lepton control sample are further subdivided into control regions (CRs) using the same selection criteria as in the search sample, according to MT(b1,2, pTmiss), Nj, Nt, and E miss. However with the requirement Nb ≥ 1 the distribution T in E miss originating from lost-lepton processes is indepen T dent of Nb, and therefore the CRs are not subdivided according to the number of b-tagged jets. These CRs generally have a factor of 2–4 more events than the corresponding SRs. The estimation of the lost-lepton background in each SR is based on the event count in data in the corresponding singlelepton CR (N1data). We translate this event count to the SR by means of a lost-lepton transfer factor TLL obtained from simulation. The lost-lepton background prediction can therefore be extracted as NLpLred = N1data TLL, TLL = where N sim and N sim are the simulated lost-lepton back0 1 ground yields in the corresponding zero- and single-lepton regions, respectively, taking into account contributions from tt and W+jets events, with smaller contributions from single top quark and ttW processes. The contamination from other SM processes in the single-lepton CRs is found to be negligible in studies of simulated events. Monte Carlo simulated samples are used to estimate the small component of the lost-lepton background that originates from leptons falling outside the kinematic acceptance, since this component is not accounted for in the CRs. To improve the statistical power of the estimation, CRs with Nt ≥ 1 are summed over E miss bins as well as over T Nb. We rely on the simulation through N sim to provide the 0 E miss-dependence and to predict the yield in each of the SRs T with Nt ≥ 1. We check this procedure by computing the data-to-simulation ratios N data/N sim in the higher-statistics 1 1 region of MT(b1,2, pmiss) ≥ 175 GeV with Nt = 0, and T find no evidence of a dependence on E miss. We assign the T relative statistical uncertainties of these ratios as systematic uncertainties in the SRs. The dominant uncertainty in the lost-lepton prediction is due to the limited number of events in the CRs, and can be as large as 100%. The statistical uncertainties in the simulated samples also affect the uncertainty in the prediction via the transfer factors. The effect in the uncertainty ranges between 3 and 50%. A source of bias in the prediction can arise from a possible difference between data and simulation in the background composition, which is assessed by independently changing the cross sections of the W+jets and tt processes by ±20% based on CMS differential cross section measurements [65, 66]. The effect of these changes is as large as 11% for the transfer factors. The uncertainties in the measurements of correction factors in lepton efficiency that are applied to the simulation to reduce discrepancies with the data lead to a systematic uncertainty of up to 7% in TLL. All other sources of systematic uncertainty, to be discussed in Sect. 7, have a negligible effect on the prediction. 4.2.2 Estimation of the Z → νν background Two methods are traditionally used to estimate the Z → νν background in searches involving all-jet final states with large E miss. The first method relies on a sample dominated by T Z → +jets events, which has the advantage of accessing very similar kinematics to the Z → νν process, after correcting for the difference in acceptance between charged-lepton pairs and pairs of neutrinos, but is statistically limited in regions defined with stringent requirements on jets and E miss. T The second method utilizes γ +jets events that have a significantly larger production cross section than the Z → +jets process, but similar leading-order Feynman diagrams. The two main differences between the processes that must be taken into account, namely, different quark–boson couplings and the massive nature of the Z boson, become less important at large Z boson pT, which is the kinematic region we are probing in this search. We have therefore adopted a hybrid method to estimate the Z → νν background by combining information from Z+jets, with Z → , and γ +jets events. Z → events are used to obtain the normalization for the Z → νν background in different ranges of Nb to account for potential effects related to heavy-flavour production, while the much higher yields from the γ +jets sample are exploited to extract corrections to distributions of variables used to characterize the SRs. The Z → events are obtained from dielectron and dimuon triggers, with the leading lepton required to have pT > 20 GeV, and the trailing lepton pT > 15 and > 10 GeV for electrons and muons, respectively. Both leptons must also have |η| < 2.4. The γ +jets sample is collected through a singlephoton trigger, and consists of events containing photons with pT > 180 GeV and |η| < 2.5. The transverse momentum of the dilepton or photon system is added vectorially to pmiss T in each event of the corresponding data samples to emulate the kinematics of the Z → νν process. The modified E miss, T denoted by ETmiss, and ETmiss,γ for the Z → and γ +jets processes, respectively, is used to calculate related kinematic variables. The prediction for the Z → νν background is given by: ETmiss,γ distributions for Z → +jets and γ +jets samples, respectively, and found them to agree. The factor RZ is calculated by comparing the observed and expected Z → yields for a relaxed version of the baseline selection. In particular, we remove the requirements on φ1234 after confirming that this does not bias the result, and relax the requirements on ETmiss, from a threshold of 250 GeV to a threshold of 100 GeV. To increase the purity of the Z → events, we require the dilepton invariant mass to lie within the Z boson mass window of 80 < M < 100 GeV. The normalization of the nonnegligible tt contamination is estimated in the region outside the Z boson mass window (20 < M < 80 or M > 100 GeV) and taken into account. Small contributions from tZ and ttZ production, estimated from simulation, are included in the Z → sample when measuring RZ. Contributions from tW and ttW are included in the simulation sample used to obtain the normalization factor for the tt contamination. As discussed previously, we calculate RZ separately for different Nb requirements. The values obtained are 0.94 ± 0.13 and 0.84 ± 0.19 for Nb = 1 and ≥2, respectively. The uncertainty in RZ originates from the limited event counts in data and simulation, and from the extrapolation in E miss. T The quantity Sγ is the correction factor related to the modelling of the distributions in the kinematic variables of Z νν events. It is calculated via a comparison of the E →miss,γ distributions of γ +jets events in simulation T and data. The simulation is normalized to the number of events in data after applying the baseline selection. To suppress potential contamination from signal and avoid overlap with the search sample, we only consider events with E miss < 200 GeV. The Sγ factor is estimated separately for T each SR to account for any potential mismodelling of the observables MT(b1,2, E miss,γ ), Nj, ETmiss,γ , and Nt in simula T tion. Since no statistically significant dependence of E miss,γ T on Nb is observed, we improve the statistical power of the correction by combining the Nb = 1 and Nb ≥ 2 subsets of the γ +jets sample to extract the Sγ corrections. The correction factors range between 0.3 and 2, with uncertainties of up to 100% due to the limited number of events in the data sample. The γ +jets control data have contributions from three main components: prompt photons produced directly or via fragmentation, and other objects misidentified as photons. The prompt photon purity measured in Ref. [33] shows good agreement between data and simulation. In addition, the impact of varying the fraction of misidentified photons, or those produced via fragmentation, by 50% in simulated events results in a bias of less than 5% in the E miss distribuT tion from the predicted Z → νν background. We therefore rely on simulation to estimate the relative contributions of the three different components. The statistical uncertainty in the γ +jets control data and the uncertainty in RZ are the main sources of uncertainty in the Z → νν prediction. The statistical uncertainties in the simulated samples, ranging up to 50% in both the SRs and in the γ +jets CRs, also makes sizeable contributions. 4.2.3 Estimation of the QCD background The QCD background is estimated using a data CR selected with the same trigger as the SR and enriched in QCD events by imposing a threshold on the azimuthal separation between pmiss and the closest of the three leading jets, T namely φ123 < 0.1. After correcting for the contribution from other SM processes (i.e. tt and W+jets), estimated by applying the normalization factor obtained in the corresponding single-lepton control sample to simulation, we translate the observation in this CR to a prediction in the SR by means of transfer factors obtained from simulation. Each transfer factor is defined as the ratio of the expected QCD events satisfying φ1234 > 0.5 to the expected QCD events with φ123 < 0.1. The estimation is carried out in each search category. Since the distributions in key observables show little dependence on Nb, the QCD CR is summed over Nb to improve the statistical precision of the estimation. The main source of QCD events populating the SR is from severe mismeasurement of the pT of one or more jets in the event. Correct modelling of jet mismeasurement in simulation is therefore an important part of the QCD prediction. The level of mismeasurement of a simulated event is parameterized by the jet response of the most mismeasured jet, which is the jet with the greatest absolute difference between the reconstructed and generated pT. The jet response, rjet, is defined as the ratio of the reconstructed pT of a jet to its generated pT, computed without including the loss of visible momentum due to neutrinos. We use the observable rjpesteudo, defined as the ratio of the pT of a jet to the magnitude of the vector sum of its transverse momentum and pmiss, as T an approximate measure of the true jet response in data, and extract mismeasurement correction factors for the simulation by comparing rjpesteudo of the jet closest in φ to pmiss between T data and simulation. The correction factors extracted from simulation are parameterized by rjet and the flavour of the most mismeasured jet. The correction factors range between 0.44 and 1.13, and are applied in the simulation on an eventby-event basis. The largest sources of uncertainty in the QCD prediction originate from the limited event counts in data and simulated samples surviving the selection, giving rise to uncertainties of up to 100% in the estimated QCD background contribution in some SRs. The uncertainty due to jet response corrections is up to 15%, while the uncertainty due to contributions from non-QCD processes in the data CR ranges from 7 to 35%. 4.2.4 Estimation of the ttZ background Contributions from the ttZ process are generally small since this is a relatively rare process. However, it has a final state very similar to signal when the Z boson decays to neutrinos and both top quarks decay only into jets, which can constitute up to 25% of the total SM background in some SRs with large E miss and Nt ≥ 1. The ttZ prediction is obtained from T simulation. We assign a 30% uncertainty to the ttZ cross section, based on the 8 TeV CMS measurement [67]. Additional theoretical and experimental uncertainties in the prediction are evaluated as will be discussed in Sect. 7, and range up to 25 and 20%, respectively, depending on the SR. We also take into consideration the statistical uncertainty in the simulation, which ranges from 5 to 100% for regions with small ttZ contributions. 4.3 Results Figure 3 shows the yields in each of the SR bins, as well as the predicted SM backgrounds based on the background estimation methods discussed in Sect. 4.2. The results are also summarized in Table 2. Expected yields are also shown for two benchmark models for the pure t1 → t(∗)χ10 decay and one for the mixed (t1 → tχ10 or t1 → bχ1±) decay. No statistically significant deviation from the SM prediction is observed in the data. 5 Search for top squarks in the single-lepton final state We also perform a search for top squarks in events with exactly one isolated electron or muon and considerable E miss. The main SM backgrounds originating from tt and T W+jets processes are suppressed using dedicated kinematic variables. The dominant remaining backgrounds arise from lost-lepton processes and the surviving W+jets background, both of which are estimated from control samples in data. 5.1 Analysis strategy The search sample is selected using triggers that require either large E miss or the presence of an isolated electron T or muon. The combined trigger efficiency for a selection of E miss > 250 GeV and at least one lepton, as measured T in a data sample with large HT, is found to be 99% with an asymmetric uncertainty of +1%. Selected events are required −3 to have at least two jets with pT > 30 GeV, at least one of which must be b-tagged using the medium working point. We require exactly one well-identified and isolated electron or muon with pT > 20 GeV, |η| < 1.442 or < 2.4, respectively, and Irel < 0.1. Electrons in the forward region of the detector are not considered in this search due to a significant 10−1 10−1 ts 104 n e v E 103 10−1 10−1 Nb = 1 Nb ≥ 2 Nb ≥ 2 Nb = 1 Nb ≥ 2 Nb = 1 Nb ≥ 2 0 250-300 300-400 400-500 500-600 > 600 250-300 300-400 400-500 500-600 > 600 Emiss [GeV] T 0 250-300 300-400 400-500 500-600 > 600 250-300 300-400 400-500 500-600 > 600 Emiss [GeV] T 10−1 Nb = 1 Nb ≥ 2 400–500, 500–600, > 600 GeV for Nb = 1, and the second five bins correspond to the same E miss binning for Nb ≥ 2. The SM background T predictions shown do not include the effects of the maximum likelihood fit to the data. The ratio of the data to the SM prediction extracted from CRs is shown in the lower panel of each plot. The shaded black band represents the statistical and systematic uncertainty in the background prediction Table 2 Observed and predicted background yields in the different search regions for the top squark search in the all-jet final state. The total uncertainty is given for each background prediction MT(b1,2, pTmiss) < 175 GeV, 5 ≤ Nj ≤ 6, Nb = 1 250–300 60 ± 6 14 ± 3 4.1 ± 1.7 300–400 23 ± 3 7.4 ± 1.9 1.5 ± 0.8 400–500 2.5 ± 1.0 1.6 ± 0.8 0.21 ± 0.15 500–600 1.9 ± 1.0 0.25 +−00..2275 0.14 +−00..1154 >600 0.28 +−00..3218 0.13 −+00..1135 0.01 ± 0.01 MT(b1,2, pTmiss) < 175 GeV, 5 ≤ Nj ≤ 6, Nb ≥ 2 250–300 61 ± 6 4.7 ± 1.4 1.1 ± 0.5 300–400 24 ± 3 3.0 ± 1.0 0.44 ± 0.23 400–500 2.8 ± 1.2 0.61 ± 0.33 0.16 ± 0.13 500–600 1.7 ± 0.9 0.13 −+00..1135 0.05 −+00..0056 >600 0.38 +−00..4318 0.04 +−00..0064 <0.01 MT(b1,2, pTmiss) < 175 GeV, Nj ≥ 7, Nb = 1 250–300 30 ± 4 3.0 ± 1.0 1.8 ± 0.6 300–400 17 ± 3 4.6 ± 1.6 1.1 ± 0.5 400–500 2.9 ± 0.9 0.82 ± 0.64 0.40 ± 0.27 500–600 1.3 ± 0.7 0.09 −+00..0191 0.05 ± 0.05 >600 <0.56 0.39 +−00..4369 0.02 ± 0.02 MT(b1,2, pTmiss) < 175 GeV, Nj ≥ 7, Nb ≥ 2 250–300 36 ± 4 0.96 ± 0.38 1.1 ± 0.5 300–400 20 ± 3 2.1 ± 0.9 0.34 ± 0.19 400–500 4.5 ± 1.4 0.15 ± 0.13 0.07 ± 0.05 500–600 1.5 ± 0.8 0.09 +−00..1019 0.01 ± 0.01 >600 <0.59 0.10 −+00..1102 0.01 ± 0.01 MT(b1,2, pTmiss) ≥ 175 GeV, 5 ≤ Nj ≤ 6, Nt = 0, Nb = 1 250–300 20 ± 3 12 ± 3 0.66 ± 0.37 300–400 9.6 ± 2.3 17 ± 4 0.63 ± 0.32 400–500 4.4 ± 1.9 8.6 ± 2.6 0.52 ± 0.35 500–600 0.82 ± 0.63 3.8 ± 1.8 0.40 ± 0.35 >600 <0.4 1.2 ± 0.7 0.05 ± 0.05 MT(b1,2, pTmiss) ≥ 175 GeV, 5 ≤ Nj ≤ 6, Nt = 0, Nb ≥ 2 250–300 11 ± 2 4.5 ± 1.4 0.45 ± 0.27 300–400 4.9 ± 1.2 6.3 ± 1.8 0.37 ± 0.23 400–500 1.6 ± 0.7 3.1 ± 1.1 0.18 ± 0.17 500–600 0.29 ± 0.24 1.4 ± 0.8 0.01 ± 0.01 >600 <0.49 0.32 ± 0.20 0.01 +−00..0021 MT(b1,2, pTmiss) ≥ 175 GeV, Nj ≥ 7, Nt = 0, Nb = 1 250–300 8.8 ± 1.9 2.5 ± 1.0 1.2 ± 0.6 300–400 7.1 ± 1.8 3.9 ± 1.5 0.76 ± 0.46 400–500 2.0 ± 0.8 1.3 ± 0.7 0.08 ± 0.07 500–600 0.38 +−00..4308 0.40 +−00..4430 0.02 ± 0.02 >600 0.28 +−00..3238 2.2 ± 1.2 0.02 +−00..0032 MT(b1,2, pTmiss) ≥ 175 GeV, Nj ≥ 7, Nt = 0, Nb ≥ 2 250–300 5.9 ± 1.3 1.2 ± 0.5 0.46 ± 0.24 300–400 3.8 ± 1.0 1.6 ± 0.7 0.08 ± 0.06 Table 2 continued >600 0.06 −+00..0067 0.73 ± 0.44 MT(b1,2, pTmiss) ≥ 175 GeV, Nj ≥ 5, Nt ≥ 1, Nb = 1 250–300 1.2 ± 0.5 0.30 ± 0.25 300–400 1.5 ± 0.7 0.34 ± 0.26 MT(b1,2, pTmiss) ≥ 175 GeV, Nj ≥ 5, Nt ≥ 1, Nb ≥ 2 250–300 0.66 ± 0.26 0.11 ± 0.09 300–400 0.92 ± 0.39 0.12 ± 0.10 0.11 −+00..1112 0.20 −+00..2202 0.22 −+00..2225 0.15 −+00..1353 0.13 −+00..1239 rate for a jet to be misidentified as an electron. To reduce the dilepton background originating from tt and tW production, events are rejected if they contain a second electron or muon with pT > 5 GeV and Irel < 0.2. A significant fraction of the remaining SM background originates from events with τh decays. This contribution is reduced by rejecting events that have an isolated τh candidate reconstructed using the HPS algorithm with pT > 20 GeV and |η| < 2.4. A further veto is placed on events containing isolated charged-hadron PF candidates with pT > 10 GeV and |η| < 2.5. Candidates are categorized as being isolated if their isolation sum, i.e. the scalar sum of the pT of charged PF candidates within a fixed cone of R = 0.3 around the candidate, is less than 6 GeV and smaller than 10% of the candidate pT. Single-lepton backgrounds originating from semileptonic tt, W+jets, and single top quark processes are suppressed through the MT of the lepton–neutrino system. Background processes containing a single lepton from W boson decay have a kinematic endpoint for MT at the W boson mass, modulo detector resolution and off-shell W boson mass effects. In this analysis we require MT > 150 GeV, which significantly reduces single-lepton backgrounds. To further reduce the tt background, we require the absolute value of the azimuthal angle between pmiss and the closest of the two highest- pT T jets, φ12, to be larger than 0.8, since the events that satisfy the E miss and MT requirements tend to have higher- pT T top quarks, and therefore smaller values of φ12 than signal events. The remaining background after the preselection is dominated by dilepton events from tt and tW production, where one of the leptons is not reconstructed or identified, and the presence of the additional neutrino from the second leptoni0.02 +−00..0032 0.13 +−00..1173 0.12 +−00..1162 0.02 +−00..0032 0.04 +−00..0054 cally decaying W boson makes it possible to satisfy the MT requirement. Kinematic properties of signal events such as E miss, MT, T and jet multiplicity depend on the decay modes of top squarks, as well as on the mass splittings ( m) between the top squark, neutralino, and chargino (if present). As a basis for the search strategy in the topologies shown in Fig. 1a, b, we require the presence of at least four jets. Events are then categorized based on the value of the MTW2 variable [68], which is calculated for each event under the assumption that it originates from the dilepton tt process with a lost lepton: MTW2 ≡ Min{my, consistent with: × [ p12 = 0, ( p1 + p )2 = p22 = M W2, pT1 + pT2 = ETmiss, ( p1 + p + pb1 )2 = ( p2 + pb2 )2 = my2]}, (6) where my is the fitted parent particle mass, and p1, p , p2, pb1 , and pb2 are the four momenta of the neutrino corresponding to the visible W boson decay, the lepton from the same decay, the W boson whose decay gives rise to the undetected lepton, and the two b jet candidates, respectively. To select the b jet candidates, we examine all possible pairings with the three jets that have the highest CSV discriminator values. The pairing that gives the lowest value of MTW2 defines the final estimate. The reconstruction of an event using the MTW2 variable helps discriminate signal from the dominant dilepton tt background. For large mass differences between the top squark and the neutralino, the MTW2 > 200 GeV requirement significantly reduces the background while maintaining reasonable signal efficiency. In contrast, for small- m models, such a requirement results in a significant loss in signal efficiency. To preserve sensitivity to both high- and low- m scenarios, we subdivide the search sample into two event categories with MTW2 > 200 GeV and ≤ 200 GeV. The MTW2 distribution for events with at least four jets is shown in Fig. 4 (top). In signals with a large difference in mass between the top squark and the neutralino, a significant fraction of events can contain two quarks that merge into a single jet as a result of the large boost of the top quark or W boson that decay into jets. These events would fail the four-jet requirement. To recover acceptance for such topologies, we define an additional SR in events with three jets. Since this region targets large m signal scenarios, only events with MTW2 > 200 GeV are considered. To increase the sensitivity of this analysis to a mixed decay scenario (Fig. 1c) when the chargino and neutralino are nearly degenerate in mass, SRs with exactly two jets are added. In events with low jet multiplicity the modified topness variable (tmod) [69] provides improved dilepton tt rejection: This equation uses the mass constraints for the particles and also the assumption that pTmiss = pT,W + pT,ν . The first term constrains the W boson whose lepton decay product is the detected lepton, while the second term constrains the top quark for which the lepton from the W boson decay is lost in the reconstruction. Once again, we consider all possible pairings of b jet candidates with up to three jets with highest CSV discriminator values. The calculation of modified topness uses the resolution parameters aW = 5 GeV and at = 15 GeV, which determine the relative weighting of the mass shell conditions. We select events with tmod > 6.4. The definition of topness used in this analysis is modified from the one originally proposed in Ref. [69]: namely, the terms corresponding to the detected leptonic top quark decay and the centre-of-mass energy are dropped since in events with low jet multiplicity the second b jet is often not identified. In these cases, the discriminating power of the topness variable is reduced when a light-flavour jet is used instead in the calculation. The modified topness is more robust against such effects and provides better signal sensitivity in these SRs than the MTW2 variable. The distribution of modified topness for events with at least two jets is shown in Fig. 4 (bottom). Finally, events in each of the categories described above are further classified into different SRs based on the value of E miss. This results in a total of nine exclusive SRs as T summarized in Table 3. Fig. 4 The MTW2 (top) and tmod (bottom) distributions for signal and backgrounds after the preselection are shown. The MTW2 variable is shown for events with four or more jets, while tmod is shown for events with at least two jets. Signal models with different top squark and neutralino mass hypotheses are shown for comparison 5.2 Background estimation Three categories of backgrounds originating from SM processes remain after the preselection described in Sect. 5.1. The dominant contribution arises from backgrounds with a lost lepton, primarily from the dilepton tt process. A second class of background events originates from SM processes with a single leptonically decaying W boson. Preselection requirements of E miss > 250 GeV and MT > 150 GeV T strongly suppress this background. The suppression is much stronger for events with a W boson originating from the decay 10−1 10−2 ts 104 n e v E 103 10−1 2.3 fb-1 (13 TeV) 2.3 fb-1 (13 TeV) of a top quark than for direct W boson production, as the mass of the top quark imposes a constraint on MW. As a result, large values of MT in semileptonic tt events are dominated by E miss resolution effects, while for events in which T the W boson is produced directly (W+jets) they are mainly a function of the width of the W boson. The third class of background events includes rare SM processes such as WZ and ttZ (where the Z boson decays to neutrinos), with smaller contributions from ttW, ttγ , and processes with two or three electroweak vector bosons. The QCD background is negligible in this search due to requirements on the presence of a high- pT isolated lepton, large E miss, and large MT. T 5.2.1 Lost-lepton background The lost-lepton background is estimated from data in dilepton CRs, where we require the presence of a second lepton passing the rejection requirements but with pT > 10 GeV, an isolated track, or a τh candidate. This is done again by extrapolating the data in the dilepton CRs to the SRs using transfer factors obtained from simulation. We use the same preselection requirements on E miss and MT as in the search regions. T We remove the subdivision in E miss and the separation of the T three and at least four jet regions to increase the statistical power of the CRs, and arrive at three CRs: exactly two jets and tmod > 6.4, at least three jets and MTW2 ≤ 200 GeV, and at least three jets and MTW2 > 200 GeV. These control regions have a purity in dilepton events of >97%. Additional transfer factors are therefore needed to account for the extrapolation in jet multiplicity and E miss requirements; these are derived T from simulation. The background estimate can be written as follows: where N data is the number of events observed in data in the 2 dilepton CR. The largest systematic uncertainty in the background estimate is due to the statistical uncertainties of the event yields in data CRs and the estimates from simulated samples (10–30%). The signal contamination in this CR is around 10% for the bulk of the studied parameter space and is taken into account in the final interpretation. The transfer factor TLL is obtained from simulation, and estimates the probability that a lepton is not identified in the detector, accounting for the kinematic acceptance and the efficiency of the lepton selection criteria. The second transfer factor, TEmiss,Nj , extrapolates the inclusive estimate to individual SR T bins. This transfer factor, also obtained from simulation, is validated by checking the modelling of the jet multiplicity and of the E miss spectrum in dedicated data CRs, which will T be described in the following paragraphs. The dilepton tt background contributes to the SRs with three or more jets only if jets from ISR or final-state radiation (FSR) are also present, or when a τh decay is misidentified as an additional jet. The modelling of jet multiplicity is checked in a high-purity dedicated dilepton data control sample with one electron and one muon, at least two b-tagged jets, and E miss > 250 GeV. The differences between data and simula T tion are used to estimate scale factors relative to the baseline selection of events with at least two jets. The scale factors are 1.10 ± 0.06 for three-jet events and 0.94 ± 0.06 for events with at least four jets. Within statistical uncertainties, these factors display no E miss dependence. The scale factors are T applied to the dilepton tt simulation when extrapolating the inclusive background prediction into the specified jet multiplicity bins. The statistical uncertainties in these scale factors are also propagated to the predictions in the SRs. The uncertainty in the modelling of the jet multiplicity ranges up to 3%. The extrapolation in E miss is carried out through simula T tion, and it must be verified that its resolution is accurately modelled. Changing the resolution can lead to a different E miss spectrum. In this analysis we are interested in the effect T of the E miss resolution in events containing intrinsic E miss T T because of the presence of neutrinos in the events. This effect is estimated by comparing a γ +jets sample in data with simulation. The events are selected using a single-photon trigger with pT > 165 GeV and |η| < 2.4. Photons are required to pass stringent identification criteria. We use the photons to mimic the neutrinos in the event, with the photon momentum serving as an estimate of the sum of the neutrino momenta. The photon pT spectrum in data and in simulation is reweighted to match that of the neutrinos in the backgroundsimulation sample. For dilepton tt events, this corresponds to the νν- pT spectrum. To model the E miss resolution, the T transverse momentum of the photon system is added vectorially to the pmiss and the resulting E miss spectrum is compared T T between data and simulation. We use this modified E miss defT inition to calculate our discriminants. For this CR, we then apply selection criteria close to the SR criteria, except that selections related to the lepton are dropped, the presence of a well-identified photon is required, and the requirement of a b-tagged jet is reversed so as to suppress effects related to semileptonic heavy-flavour decays. Corrections for the observed differences, which can go up to 15%, are applied to events in the simulated samples and the uncertainties propagated to the final background estimate, resulting in an uncertainty of 1–4% in the lost-lepton background prediction. 5.2.2 One-lepton background In SRs with a high MTW2 or modified topness requirement, the W+jets background is estimated using a data control sample containing no b-tagged jets. For SRs with a low-MTW2 requirement, this background constitutes less than 10% of the total SM background. In these SRs we do not employ an estimate based on data, but instead use the W+jets background estimate directly from simulation. The semileptonic tt background is also estimated from simulation. The CRs used to extract the W+jets background in the SRs with a high MTW2 or modified topness requirement are again not subdivided in E miss to have a sufficient number of T events to carry out the prediction. We therefore use three CRs for this background estimate: exactly two jets with tmod > 6.4, exactly three jets with MTW2 > 200 GeV, and at least four jets with MTW2 > 200 GeV. We extrapolate the yields from the CRs to the SRs by applying transfer factors from simulation for the extrapolation in E miss and number of b-tagged jets: T N Wpr+edjets = (N Ndabt=a0 − N Nnobn=-W0+jets) TEmiss TNb , T with N Ndabt=a0 − N Nnobn=-W0+jets representing the event yield in the CR after subtracting the estimated contribution from other SM background processes. The non-1 contribution in the CRs, N non-W+jets, is estimated from simulation and amounts Nb=0 to roughly 25–35%. A 50% uncertainty is assigned to the subtraction. The largest source of uncertainty is again the limited size of the data and simulation samples. The statistical uncertainty of these samples results in an uncertainty of 20– 40% in the W+jets background estimate. The transfer factor TEmiss extrapolates the yields from the T inclusive CR with E miss > 50 GeV to the exclusive E miss T T regions. The main uncertainties in this extrapolation factor can be attributed to the modelling of the neutrino pT spectrum, the W boson width, and the E miss resolution. The neuT trino pT spectrum is checked in a data sample enriched in W+jets, with no b-tagged jets and 60 < MT < 120 GeV. No large mismodelling of E miss is observed. Therefore, we T do not apply any corrections to the neutrino pT spectrum but only propagate the statistical limitation of this study as the uncertainty (6–22%) in the modelling of the neutrino pT spectrum. The uncertainty in the W boson width (3% [70]) is estimated by scaling the four-vectors of the W boson decay products appropriately. The E miss resolution effects on this T background are studied using the same method as described in Sect. 5.2.1, giving rise to a 1–3% uncertainty. The other transfer factor, TNb , performs the extrapolation in the number of b-tagged jets for each E miss bin. Scale facT tors are applied to the simulation to match the b tagging efficiency in data. The largest uncertainty in this transfer factor is the fraction of the heavy-flavour component in the W+jets sample; we assign a 50% uncertainty to this component. We performed a dedicated cross-check in a CR with one or two jets and at least 50 GeV of E miss. Data and simulation were T found to be in agreement in the b jet multiplicity within uncertainties. After taking into consideration the additional sources of systematic uncertainty described in Sect. 7, the total uncertainty in the W+jets estimate varies from 50 to 70%. The semileptonic tt background is never larger than 10% of the total background estimate. We rely on simulation to estimate it. The main source of uncertainty in this estimate is the modelling of the E miss resolution because poor resolu T tion can enhance the contributions at large MT. The studies of E miss resolution presented in Sect. 5.2.1 indicate that it T could be mismodelled by about 10% in simulation. Changes in the simulated E miss resolution by a corresponding amount T provide an uncertainty of 100% in the semileptonic tt estimate. 5.2.3 Rare standard model backgrounds The “rare” background category includes tt production in association with a vector boson (W, Z, or γ ), diboson, and triboson events. Within this category, WZ events dominate the SRs with two jets, and ttZ events with the Z boson decaying into a pair of neutrinos (Z → νν) dominate regions of higher jet multiplicity. The expected contributions from these backgrounds are small, and the simulation is expected to model the kinematics of these processes well in the regions of phase space relevant to the SRs. The rare backgrounds are therefore estimated using simulation. We assess the theoretical and experimental uncertainties affecting the estimates as described in Sect. 7, resulting in a total uncertainty of 15– 26%, depending on the SR. 5.3 Results The background expectations and the corresponding yields for each SR are summarized in Table 4 and in Fig. 5. Overall, the observed and predicted yields agree within two standard deviations (SD) in all SRs. For signals of top squark pair production for different mass hypotheses, the maximum observed significance obtained by combining the results in different SRs is 1.2 SD for a top squark mass of ≈400 GeV 1 (not from top) ts 40 en 35 v E30 2.3 fb-1 (13 TeV) 0 250 −350 > 350 250 −350 > 350 250 −325 > 325 250 −350 350 −450 > 450 Emiss [GeV] T Fig. 5 Background estimates from data and simulation, together with the observed yields in the SRs of the single-lepton analysis, described in Table 3. The uncertainties, which are the quadratic sums of statistical and systematic uncertainties, are indicated by the cross-hatched areas. The SM background predictions shown do not include the effects of the maximum likelihood fit to the data. Three signal hypotheses are overlaid. The hypothesis t1 → tχ10/t1 → bχ1± has branching fractions B(t1 → tχ10) = B(t1 → bχ1±) = 0.5 and a massless LSP hypothesis. We therefore find no evidence for top squark pair production. 6 Search for pair production of bottom squarks or of top squarks decaying to charm quarks This search is motivated by the production of pairs of bottom or top squarks, in which each b1 or t1 decays, respectively, into a bottom or a charm quark and a neutralino. In the latter search, the difference between the t1 and χ10 masses is assumed to be less than 80 GeV, and the only top squark Table 4 Background estimates from data and simulation, and observed data yields for the single-lepton top squark analysis using 2.3 fb−1 of data collected during 2015 pp collisions. The uncertainties are the quadratic sums of statistical and systematic uncertainties decay mode considered is through a flavour changing neutnreamultcr=uarlirmneonbt1lte−oacvmχeχ10l10i.tStblmeetwvalielseimbnlaetshseesnpteolripgttyoinrignbsot htteommde=steqcmutaot1rrk−,mamnadχk10itnhoger signal events difficult to distinguish from SM background. However, events with an energetic ISR jet recoiling against the pmiss originating from the neutralino can provide a dis T tinct topology for signals with compressed mass spectra, i.e. with small m. We thus perform a search for events with an ISR jet and significant E miss. T 6.1 Analysis strategy Events in the search sample are recorded using the same trigger as that for the top squark search in the all-jet final state, requiring the presence of large E miss and at least two enerT getic jets within the tracker acceptance. After applying an offline selection requiring E miss > 250 GeV and at least two T jets with pT > 60 GeV, we find the trigger efficiency to be greater than 97%. We veto events that have at least four jets with pT above 50 GeV. The veto and its threshold are motivated by the harder pT spectrum of the fourth jet in semileptonic tt events compared to the signal, in which extra jets originate from ISR or FSR. To reduce the SM background from processes with a leptonically decaying W boson, we reject events containing isolated electrons or muons with Irel < 0.1 and |η| < 2.5, or Irel < 0.2 and |η| < 2.4, respectively, and with pT > 10 GeV. The contribution containing τh decays is reduced by placing a veto on events containing chargedhadron PF candidates with pT > 10 GeV, |η| < 2.5, and an isolation sum smaller than 10% of the candidate pT. The dominant SM background sources are Z+jets production with Z → νν, and the lost-lepton background originating from W+jets, tt, and single top quark processes with leptonic W boson decays. A smaller background contribution comes from QCD events in which large E miss originates T from jet mismeasurements and the direction of pmiss is often T aligned with one of the jets. To suppress this background we require that the absolute difference in azimuthal angle between the pmiss and the closest of the three leading jets T ( φ123) is greater than 0.4. Two sets of SRs are defined to optimize the sensitivity for signal models with either compressed or noncompressed mass spectra. In addition to the criteria discussed above, for regions targeting noncompressed scenarios we require that the pT of the leading jet be above 100 GeV and that the event contain at least one additional jet with pT above 75 GeV. We also require that the two highest- pT jets be identified as b jets. These requirements suppress events originating from W and Z boson production, for which the leading jets have a softer pT spectrum since they are produced by ISR or FSR. To maintain a stable b tagging efficiency as a function of jet pT, both the loose and medium working points of the b tagging algorithm are used to identify b jets. The b tagging efficiency of the medium working point depends strongly on the jet pT and degrades by about 20–30% for jets with pT above 500 GeV, while the efficiency of the loose working point is more stable with increasing jet pT. Specifically, we use the loose working point to identify b-tagged jets when the leading jet has pT above 500 GeV, and the medium working point otherwise. Since such high- pT jets are less likely to occur in SM processes, the higher misidentification rate of the loose working point results in only a small increase in the SM background. The distribution of MT (j1,2, pmiss) ≡ min[MT(j1, pTmiss), T MT(j2, pmiss)], where j1, j2 are the two highest- pT jets, is T expected to have a kinematic endpoint at the mass of the top quark when pmiss and the closest jet originate from the T semileptonic decay of a top quark. In the noncompressed search sample we require MT (j1,2, pTmiss) to be greater than 250 GeV. Events in this sample are then categorized by HT,12, defined for the purposes of this analysis as the scalar sum of the pT of the two leading jets, and the mCT kinematic variable. The boost-corrected cotransverse mass [71,72], mCT, is defined by: m2CT(j1, j2) = 2 pT(j1) pT(j2)[1 + cos φ (j1, j2)]. For scenarios in which two particles are pair-produced and have the same decay chain, the mCT distribution has an endpoint determined by the masses of the parent and decayproduct particles. For b1 → bχ10 this endpoint is at (m(b1)2− m(χ10)2)/m(b1). For signals with compressed mass spectra, high- pT ISR is required to be able to reconstruct the quarks as jets and obtain a large value of E miss. Compressed SRs require therefore a T leading jet with pT > 250 GeV that is back-to-back relative to the pmiss ( φ(j1, pmiss) > 2.3). Since such ISR jets are T T not expected to originate from b quarks, we require that the leading jet fail the loose b-tagging requirement. We relax the thresholds on the second jet pT and on the MT (j1,2, pmiss) to 60 and 200 GeV, respectively, and cate T gorize events in the search sample according to the number of b-tagged jets. The mCT observable loses its discriminating power for these compressed signal models due to the small mass splitting between the parent particle and χ10. The E miss is therefore used as the main discriminant, with differ T ent E miss thresholds applied to define the final SRs. T The baseline selections for both noncompressed and compressed regions are summarized in Table 5, while the definitions of the two sets of SRs are described in Table 6. 6.2 Background estimation The SM background contributions originating from Z → νν, lost-lepton, and QCD processes are estimated from dedicated data CRs as discussed below. Smaller contributions from other SM processes, such as diboson (VV) processes, are estimated from simulation, and an uncertainty of 50% is assigned to these contributions. 6.2.1 Estimation of the Z → νν background The Z → νν background is estimated from a high-purity data sample of Z → μ+μ− events in which we remove the muons and recalculate the relevant kinematic variables to emulate Z → νν events. The trigger used to collect this CR requires the presence of a high- pT muon with |η| < 2.1. In keeping with the trigger constraints, the sample is selected by requiring the presence of two isolated muons in the event with pT > 50 (10) GeV and |η| < 2.1 (2.4) for the leading (trailing) muon. The invariant mass of the dimuon pair is required to be within 15 GeV of the Z boson mass [70]. Each muon is required to be separated from jets in the event by R > 0.3. Apart from the lepton selection, we apply the same object and event selection criteria as described in Sect. 6.1 to this sample, with the exception that b jets are selected using the loose working point of the b tagging algorithm to improve the statistical power of the data CR. Events in the selected sample are subdivided into CRs corresponding to the noncompressed and compressed SRs. The observed events in these data CRs, Nμdaμta, are translated into an estimation of the Z → νν contribution in the SRs with the help of simulation, as follows: Table 5 A summary of the baseline selections used for the noncompressed and compressed b1 and t1 → cχ10 compressed SRs scenarios, and in Nb and Emiss where Nμnoμn-Z, representing the small contamination in the CRs due to tt, W+jets, single top quark, and diboson processes, is estimated from simulation. The corrected dimuon event yield is scaled by the kinematic and detector acceptance of muons from Z bosons, A, and the muon reconstruction, identification, and isolation efficiency . The acceptance and efficiency are determined from simulation. Efficiency scale factors are applied to correct for differences between data and simulation. These scale factors are determined with a “tag-and-probe” method in Z → μ+μ− events [53]. The product of the muon acceptance and efficiency, A , varies from 0.6 in the low-mCT and low-ETmiss regions to 0.9 in the high-mCT and high-ETmiss regions. The correction factor RZμμ→νν = 5.942 ± 0.019 [70] represents the ratio of the Z boson branching fractions to neutrinos and leptons. The remaining term, κ, accounts for differences in the b tagging efficiency and misidentification rate between the CRs and SRs, resulting from the use of different b tagging working points. These κ factors are determined from Z → simulation and corrected for known differences in the performance of the b tagging algorithm between data and simulation as measured in samples of multijet and tt events [44]. The value of the b tagging κ factor ranges from 0.10 to 0.15 for the noncompressed SRs, and from 0.20 to 0.25 for the Nb = 1 compressed SRs, while it is about 0.15 for the Nb = 2 compressed SR. The largest uncertainty in the Z → νν background estimate arises from the limited event yields in the dimuon CR, corresponding to a 10–100% uncertainty in the Z → νν prediction. We correct for the estimated contributions to the CR from SM processes other than Z → μ+μ− using simulation samples with an assigned uncertainty of 50% in their normalization. This leads to an uncertainty of 2–20% in the background estimate. Other experimental and theoretical sources of uncertainty, to be discussed in Sect. 7, result in an additional 2–8% uncertainty in A , and a 2% uncertainty is assigned in all SRs to account for the uncertainty in the Z boson branching fractions. The uncertainty in the b tagging κ factors is assessed by varying the data-to-simulation b tagging correction factors according to their measured uncertainties. Additionally, the dependence of κ on the heavyflavour content in Z boson events is evaluated by varying the Z+bb and Z+cc fractions in simulation by 20% based on the uncertainty in the CMS Z+bb measurement [73], resulting in an additional uncertainty of 10–20% in the Z → νν estimate. 6.2.2 Estimation of the lost-lepton background The lost-lepton background in each SR is estimated from a single-lepton CR in data selected by inverting the electron and muon vetoes in events collected with the same trigger as used to record the signal sample. We relax the b tagging requirement in the CRs using the loose working point in the noncompressed selection, while keeping the same requirement as in the SRs for the compressed regions. In all other respects, the CRs are defined through the same selection criteria as the corresponding SRs, including requirements on the HT,12, mCT, Nb, and E miss, to remove any dependence T of the prediction on the modelling of these kinematic variables in simulation. The possible contamination from signal in the single-lepton CR is negligible, less than 1%, so no extra requirement on MT( , pTmiss) is made. The lost-lepton component of the SM background in each SR, NLpLred, is estimated once again from the corresponding data via a transfer factor, TLL, determined from simulation: where N data is the observed event yield in the single-lepton 1 CR. The transfer factor TLL accounts for effects related to lepton acceptance and efficiency. The largest uncertainty in the lost-lepton estimate is, as in the previous analyses, due to the statistical uncertainty in the event yields, ranging from 3 to 50%, depending on the SR. Contributions to the CRs from Z → and diboson processes are subtracted using estimates from simulation, and a 50% uncertainty is applied to this subtraction, which leads to an uncertainty of 3–10% in the lost-lepton prediction. The limited event counts in the simulation sample result in a 2– 12% uncertainty, while uncertainties related to discrepancies between the lepton selection efficiency in data and simulation give rise to a 3–4% uncertainty in the final estimate. An additional uncertainty of 7% in the τh component accounts for differences in isolation efficiency between muons and singleprong τh decays, as determined from studies with simulated samples of W+jets and tt events. A systematic uncertainty of 5–10% is found for the uncertainties in the b tagging scale factors that are applied to the simulation for the differences in b tagging performance between data and simulation and the different b tagging working points. Finally, we estimate a systematic uncertainty in the transfer factor to account for differences in the tt and W+jets admixture in the search and control regions. This results in a 20–30% uncertainty in the final prediction. 6.2.3 Estimation of the QCD background The φ123 > 0.4 requirement reduces the QCD contribution to a small fraction of the total background in all SRs for both compressed and noncompressed models. We estimate this contribution for each SR by applying a transfer factor to the number of events observed in a CR enriched in QCD events. The CRs are obtained by inverting the φ123 requirement. The transfer factor, TQCD, is measured in a sideband region in data with ETmiss ∈ [200, 250] GeV and the same requirements on the other variables as in the SRs. This factor is the ratio between the number of QCD events in the φ123 > 0.4 and φ123 < 0.4 subsets of this sideband region. The estimated contribution of other SM processes (tt, W+jets, single top quark, and diboson production) based on simulated samples is subtracted from the event yields in the CR and each subset of the sideband. The transfer factor for the noncompressed regions does not vary significantly as a function of HT,12 and mCT. Therefore, we extract the value of TQCD used for the noncompressed SRs from a sideband selected with an inclusive requirement on HT,12 and mCT to reduce the statistical uncertainty in the transfer factor. The transfer factors for the compressed SRs are obtained from sidebands that are subdivided according to the number of b-tagged jets into Nb = 0 and Nb ≥ 1 regions, with the latter used to extract the QCD predictions for the Nb = 1 and Nb ≥ 2 SRs. The statistical uncertainties due to the limited number of events in the data CRs and the non-QCD simulated samples are propagated to the final QCD estimate, ranging from 10 to 100%. The main uncertainty in TQCD also originates from the statistical uncertainty of the observed and simulated event yields in the sideband region. We assign additional uncertainties for differences in b tagging efficiency between data and simulation and for the subtraction of the non-QCD background contribution in the sideband. The total systematic uncertainty in the QCD prediction varies between 27 and 76% in the compressed SRs, but can be as large as 550% in the noncompressed SRs due to the small event samples in the corresponding sideband in data. 6.3 Results The expected SM background yields and the number of events observed in data are summarized in Table 7 and shown in Fig. 6. The observed yields agree well with the predicted SM background. 7 Systematic uncertainties Several categories of systematic uncertainties apply to all three analyses. These include uncertainties arising from the limited event counts in control samples, uncertainties related to the use of simulation in SM background predictions, and a 2.7% uncertainty in integrated luminosity [74] that applies to the estimated signal yields and contributions from rare background processes that are taken directly from simulation, without the use of data control samples. Nb = 0 Nb = 1 Nb = 2 HT,12 > 500 GeV Table 7 Observed number of events and background prediction in the different SRs for the b1 and t1 → cχ10 searches. The total uncertainty in the background predictions is also shown The limited number of simulated events surviving the stringent requirements on jets and E miss in all three searches T can lead to a significant statistical uncertainty in background predictions. In the case of background predictions that rely on simulation for accurate modelling of the relevant event kinematics, we assess theoretical uncertainties, primarily those associated with missing higher-order corrections, in the simulated samples by varying the renormalization and factorization scales up and down by a factor of two [75, 76] and by variations of PDFs. The PDF uncertainties are defined by the SD obtained from 100 variations of the NNPDF3.0 [57] PDFs. The uncertainties are then propagated to the final background estimates. When the simulation of the detector response does not adequately describe the data, correction factors are applied to account for the observed discrepancies. Differences in the efficiencies for selecting isolated leptons between simulation and data are measured in Z → events in the case of electrons and muons and in a tt-enriched sample for hadronically decaying τ leptons. The observed deviations are accounted for in the form of corrections to the simulation, and the corresponding uncertainties are propagated to the predicted SM yields in the SRs. Correction factors and uncertainties based on measurements of b tagging performance in data and simulation [44] are also applied. They are parameterized by jet kinematics and flavour. We also assess an uncertainty related to the modelling of additional interactions in the simulation. For the rare SM backgrounds with top quarks, predominantly from tt production in association with a Z boson, where the Z boson decays to a pair of neutrinos, an extra uncertainty is estimated to account for the possible mismodelling of the top quark pT spectrum. The efficiency and misidentification rates for the top quark tagging algorithm are compared between data and simulation in CRs as a function of the key kinematic variables. The correction factors are found not to be strongly dependent on the different kinematic variables. The efficiency estimated in simulation agrees with the measured efficiency while the misidentification rate has to be corrected by 30%. Both correction factors have a 10% uncertainty, estimated from the variations of the efficiency measurement. All these uncertainties are propagated to the different signal and background estimates to which they apply. The back2.3 fb-1 (13 TeV) tsn 105 e v E 104 10−1 10−1 Nb=0 Nb=1 Nb=2 200 GeV < HT,12 ≤ 500 GeV Fig. 6 Observed events and estimated SM background and signal yields for the compressed (top) and noncompressed (bottom) SRs for the bottom squark search in the all-jet final state. The observed data yield is shown as black points and the total background predictions are shown in solid area. The SM background predictions shown do not include the effects of the maximum likelihood fit to the data. The bottom panel shows the ratio of data to the total background prediction in each search bin. Only statistical uncertainties are propagated to the ratio ground predictions from control samples in data are affected through the transfer factors that are calculated from simulation corrected to reproduce data. In general these uncertainties are subdominant and the uncertainty in the final background estimate is dominated by the statistical uncertainty of the data control sample. For the signal samples differences between the fast simulation and the full Geant4-based model are also taken into account. Lepton selection efficiencies and b tagging performance are found to be different in the fast simulation. We derive appropriate corrections for the fast simulation and propagate the corresponding uncertainties to the predicted signal yields. We also assess an additional uncertainty for the difference in E miss resolution between the fast simulation T and the full Geant4-based model. This difference in E miss T resolution has the largest impact on signal models with small intrinsic E miss, as is the case for compressed mass spectra. T The modelling of the ISR plays an important role in cases where the top squark and χ10 masses are very similar. The uncertainty is determined by comparing the simulated and observed pT spectra of the system recoiling against the ISR jets in tt events, using the method described in Ref. [21]. The effect is generally found to be small, although in scenarios with a compressed mass spectrum the effect can be as large as 30%. The uncertainties in the signal modelling are determined in each analysis for every SR. The dominant uncertainties in the predicted signal yield arise from the size of the simulated samples in some of the SRs (1–100%), jet energy scale corrections (1–50%), b tagging efficiency corrections used to scale simulation to data (1–35%), and ISR (1–30%). The largest uncertainties are in SRs that have small signal acceptance to a specific model. The statistical uncertainties of the signal samples are uncorrelated, whereas all other signal systematic uncertainties are considered to be fully correlated among the different SRs and analyses. Since the three analyses predict the backgrounds with different CRs, the treatment of systematic uncertainties is mostly uncorrelated among analyses, except for the estimates based on simulation. Here only the statistical component of the uncertainty is treated as uncorrelated. Systematic uncertainties due to jet energy scale corrections, b tagging efficiency and selection efficiencies are treated as correlated among the different background estimates. 8 Interpretation The data in all three searches are consistent with the background expected from SM processes. The results are interpreted as limits on supersymmetric particle masses in the context of simplified models [77–80] of top or bottom squark pair production. Different decay modes are considered for top squark pair production. For mass splittings m larger than the W boson mass, we consider two decay modes for the top squark: to a top quark and a neutralino, or to a bottom quark and a chargino, where the chargino decays to an LSP. Scenarios with t1 → t(∗)χ10 branching fractions of 50 or 100% are considered. The results of the top squark searches in the all-jet and single-lepton final states are combined for these interpretations. For m smaller than the W boson mass, only the decay of top squarks to a charm quark and an LSP is considered in this paper. For the pair production of bottom squarks, all bottom squarks are assumed to decay to a bottom quark and an LSP. The signal yield is corrected for signal contamination of data CRs for each mass hypothesis and each analysis. Typical values are around 5–10%, except for compressed mass spectra, where it can vary between 10 and 50%. The signal contamination is most significant for the top squark production models with a 100% t1 → t(∗)χ10 branching fraction, a light LSP, and m close to the top quark mass. The 95% confidence level (CL) upper limits on SUSY production cross sections are calculated using a modified frequentist approach with the CLS criterion [81,82] and asymptotic results for the test statistic [83,84]. The SRs and CRs for top squark searches in the all-jet and single-lepton final states are mutually exclusive. We combine the results of the two searches, treating the systematic uncertainties assigned to the predicted signal and background yields as correlated or uncorrelated depending on the source, as detailed in Sect. 7. Figure 7 shows 95% CL exclusion limits for pp → t1t1 → t(∗)χ 0t(∗)χ 0, assuming the top quarks in the decay to be 1 1 unpolarized, together with the upper limit at 95% CL on the excluded signal cross section. All top squarks are assumed to decay to a top quark and an LSP. For m < mt the signal samples assume a three-body decay without an off-shell top quark as intermediate particle. The expected exclusion is given by the dashed red line, with the one SD experimental uncertainty. The observed exclusion curve is shown as a solid black line together with the 1 SD uncertainty in the theoretical cross section. We do not interpret in the region near m ≈ mt when χ10 is very light because of the difficulty in modelling rapidly varying kinematics in this region. In this region an indirect search for top squark pair production can be performed by looking for a small excess in the measured tt cross section compared to the SM expectation [20,85]. We exclude top squark masses from 280 to 830 GeV for a massless LSP and LSP masses up to 260 for 675 GeV top squarks. At 8 TeV top squark masses were excluded up to 780 GeV for a massless LSP [25]. For models with heavy top squarks and light LSPs, the sensitivity is driven by the top squark analysis in the all-jet final state of Sect. 4, which is more sensitive than the single-lepton analysis (Sect. 5) because of the larger acceptance for signal. The combination extends the expected reach in top squark mass by about 45 GeV. When the LSP is heavier, the cleaner search in the single-lepton final state becomes more important. Both analyses have similar sensitivity in this area of parameter space, and combining them extends the reach in LSP mass by about 30 GeV. Figure 8 shows the 95% CL exclusion limits for t1t1 production, assuming equal probabilities for the decay modes t1 → t(∗)χ10 and t1 → bχ1±. The chargino in the latter mode decays to a W boson and an LSP. In this model, the chargino is considered to be nearly mass-degenerate with the oLrSigPin(matiχn±g =fromm 10th+e c5hGaregVin).oTdheecaWy baroesovnerdyecsaoyftpbreocdauucstes 1 χ of the small mass splitting, and might not be detectable. For intermediate LSP masses, top squark masses are probed up to 725 GeV. The LSP masses up to 210 GeV are probed for a top squark mass of around 500 GeV. Here, the single-lepton Fig. 7 Exclusion limits at 95% CL for direct top squark pair production for the decay mode t1 → t(∗)χ0. The interpretation is performed in the 1 two-dimensional space of mt1 vs. mχ10 . The color indicates the 95% CL upper limit on the product of cross section and branching fraction abtlaecakchcuprovienst rineptrheesemntt1 -thmeχ1o0bpslearnvee.dTehxecrluegsiioonnsaetn9c5l%oseCdLb,ywthheilethtihcek dashed red lines indicate the expected limits at 95% CL and their ±1 SD experimental uncertainties. The thin black lines show the impact of the ±1 SD theoretical uncertainties in the signal cross section. The magenta short-dashed curve and the blue dotted curve show the expected limits for the analysis in the all-jet (Sect. 4) and single-lepton (Sect. 5) final states, respectively. The limits in the region near m ≈ mt and low χ10 mass are not shown due to the difficulty in modelling rapidly varying kinematics in this region analysis does not contribute much to the combination because of the larger acceptance in the all-jet final state, except at low LSP masses. In most of the mass parameter space the combination reaches ≈15 GeV higher than the analysis in the all-jet final state. The compressed SRs from the bottom squark analysis in the all-jet final state (Sect. 6) are used to set upper limits on the top squark cross sections when the mass splitting between the top squark and the LSP is smaller than the mass of the W boson. Figure 9 shows the expected and observed 95% CL upper limits on the top squark cross sections in the mt1 mχ0 plane assuming the top squark always decays to a charm 1 quark and an LSP. Top squarks with masses below 240 GeV are probed in this model, when the mass splitting between the top squark and the LSP is close to 10 GeV. At 8 TeV top squark masses up to 270 GeV were probed for the same m [20]. Figure 10 shows the expected and observed 95% CL upper limits on the bottom squark cross sections in the mb1 -mχ0 1 plane using both the compressed and noncompressed SRs of the bottom squark analysis. We probe bottom squark masses up to 890 GeV for small LSP masses. With 8 TeV data bottom squark masses below 650 GeV were excluded. [20,24]. 2.3 fb-1 (13 TeV) 2.3 fb-1 (13 TeV) Fig. 9 Exclusion limits at 95% CL for direct top squark pair production with decay t1 → cχ10 using the compressed SRs of the bottom squark analysis (Sect. 6). The interpretation is done in the two-dimensional space of mt1 vs. mχ10 . The caption of Fig. 7 explains the use of lines and colors in detail 9 Summary Results are presented from three complementary searches for top or bottom squark–antisquark pairs in data collected with the CMS detector in proton–proton collisions at a centreof-mass energy of 13 TeV, corresponding to an integrated luminosity of 2.3 fb−1. The search for top squarks is carried out in the all-jet and single-lepton final states, which are combined for the final result. A second search in all-jet events is designed for bottom squark pairs and for top squarks decaying to charm quarks through a flavour changing neutral current process. No statistically significant excess of events 0 200 300 400 500 600 700 800 900 m ~t1 [GeV] Fig. 8 Exclusion limits at 95% CL for direct top squark pair production assuming equal branching fractions for the decays t1 → t(∗)χ10 and t1 → bχ1±. The interpretation is performed in the two-dimensional space of mt1 vs. mχ0 . The chargino is considered to be nearly mass1 degenerate with the LSP (mχ1± = mχ10 + 5 GeV). The caption of Fig. 7 explains the use of lines and colors in detail 2.3 fb-1 (13 TeV) Fig. 10 Exclusion limits at 95% CL for direct bottom squark pair production with decay b1 → bχ10. The interpretation is performed in the two-dimensional space of mb1 vs. mχ10 using the results of the bottom squark analysis (Sect. 6). The caption of Fig. 7 explains the use of lines and colors in detail is observed above the expected standard model background, and exclusion limits are set at 95% confidence level in the context of simplified models of direct top and bottom squark pair production. Limits for top squark masses of 830 GeV are established for a massless LSP, and for LSP masses up to 260 GeV for a 675 GeV top squark mass, when all top squarks are assumed to decay to a top quark and an LSP. When the top squarks can also decay to a bottom quark and a chargino, this reach is reduced. Assuming a mass splitting between the top squark and the LSP close to 10 GeV, and top squarks that decay to a charm quark and an LSP, top squark mass limits up to 240 GeV are established. Finally, bottom squark mass limits up to 890 GeV are established for small LSP masses. The results extend the reach with respect to previous limits obtained from LHC Run 1 data in most of the parameter space. Acknowledgements We congratulate our colleagues in the CERN accelerator departments for the excellent performance of the LHC and thank the technical and administrative staffs at CERN and at other CMS institutes for their contributions to the success of the CMS effort. In addition, we gratefully acknowledge the computing centres and personnel of the Worldwide LHC Computing Grid for delivering so effectively the computing infrastructure essential to our analyses. Finally, we acknowledge the enduring support for the construction and operation of the LHC and the CMS detector provided by the following funding agencies: the Austrian Federal Ministry of Science, Research and Economy and the Austrian Science Fund; the Belgian Fonds de la Recherche Scientifique, and Fonds voor Wetenschappelijk Onderzoek; the Brazilian Funding Agencies (CNPq, CAPES, FAPERJ, and FAPESP); the Bulgarian Ministry of Education and Science; CERN; the Chinese Academy of Sciences, Ministry of Science and Technology, and National Natural Science Foundation of China; the Colombian Funding Agency (COLCIENCIAS); the Croatian Ministry of Science, Education and Sport, and the Croatian Science Foundation; the Research Promotion Foundation, Cyprus; the Secretariat for Higher Education, Science, Technology and Innovation, Ecuador; the Ministry of Education and Research, Estonian Research Council via IUT23-4 and IUT236 and European Regional Development Fund, Estonia; the Academy of Finland, Finnish Ministry of Education and Culture, and Helsinki Institute of Physics; the Institut National de Physique Nucléaire et de Physique des Particules/CNRS, and Commissariat à l’Énergie Atomique et aux Énergies Alternatives/CEA, France; the Bundesministerium für Bildung und Forschung, Deutsche Forschungsgemeinschaft, and Helmholtz-Gemeinschaft Deutscher Forschungszentren, Germany; the General Secretariat for Research and Technology, Greece; the National Scientific Research Foundation, and National Innovation Office, Hungary; the Department of Atomic Energy and the Department of Science and Technology, India; the Institute for Studies in Theoretical Physics and Mathematics, Iran; the Science Foundation, Ireland; the Istituto Nazionale di Fisica Nucleare, Italy; the Ministry of Science, ICT and Future Planning, and National Research Foundation (NRF), Republic of Korea; the Lithuanian Academy of Sciences; the Ministry of Education, and University of Malaya (Malaysia); the Mexican Funding Agencies (BUAP, CINVESTAV, CONACYT, LNS, SEP, and UASLP-FAI); the Ministry of Business, Innovation and Employment, New Zealand; the Pakistan Atomic Energy Commission; the Ministry of Science and Higher Education and the National Science Centre, Poland; the Fundação para a Ciência e a Tecnologia, Portugal; JINR, Dubna; the Ministry of Education and Science of the Russian Federation, the Federal Agency of Atomic Energy of the Russian Federation, Russian Academy of Sciences, the Russian Foundation for Basic Research and the Russian Competitiveness Program of NRNU MEPhI (M.H.U.); the Ministry of Education, Science and Technological Development of Serbia; the Secretaría de Estado de Investigación, Desarrollo e Innovación and Programa Consolider-Ingenio 2010, Spain; the Swiss Funding Agencies (ETH Board, ETH Zurich, PSI, SNF, UniZH, Canton Zurich, and SER); the Ministry of Science and Technology, Taipei; the Thailand Center of Excellence in Physics, the Institute for the Promotion of Teaching Science and Technology of Thailand, Special Task Force for Activating Research and the National Science and Technology Development Agency of Thailand; the Scientific and Technical Research Council of Turkey, and Turkish Atomic Energy Authority; the National Academy of Sciences of Ukraine, and State Fund for Fundamental Researches, Ukraine; the Science and Technology Facilities Council, UK; the US Department of Energy, and the US National Science Foundation. Individuals have received support from the Marie-Curie programme and the European Research Council and EPLANET (European Union); the Leventis Foundation; the A. P. Sloan Foundation; the Alexander von Humboldt Foundation; the Belgian Federal Science Policy Office; the Fonds pour la Formation à la Recherche dans l’Industrie et dans l’Agriculture (FRIABelgium); the Agentschap voor Innovatie door Wetenschap en Technologie (IWT-Belgium); the Ministry of Education, Youth and Sports (MEYS) of the Czech Republic; the Council of Science and Industrial Research, India; the HOMING PLUS programme of the Foundation for Polish Science, cofinanced from European Union, Regional Development Fund, the Mobility Plus programme of the Ministry of Science and Higher Education, the National Science Center (Poland), contracts Harmonia 2014/14/M/ST2/00428, Opus 2014/13/B/ST2/02543, 2014/15/B/ST2/03998, and 2015/19/B/ST2/02861, Sonata-bis 2012/07/E/ST2/01406; the Thalis and Aristeia programmes cofinanced by EU-ESF and the Greek NSRF; the National Priorities Research Program by Qatar National Research Fund; the Programa Clarín-COFUND del Principado de Asturias; the Rachadapisek Sompot Fund for Postdoctoral Fellowship, Chulalongkorn University and the Chulalongkorn Academic into Its 2nd Century Project Advancement Project (Thailand); and the Welch Foundation, contract C-1845. Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP3. CMS Collaboration Yerevan Physics Institute, Yerevan, Armenia A. M. Sirunyan, A. Tumasyan Institut für Hochenergiephysik, Wien, Austria W. Adam, E. Asilar, T. Bergauer, J. Brandstetter, E. Brondolin, M. Dragicevic, J. Erö, M. Flechl, M. Friedl, R. Frühwirth1, V. M. Ghete, C. Hartl, N. Hörmann, J. Hrubec, M. Jeitler1, A. König, I. Krätschmer, D. Liko, T. Matsushita, I. Mikulec, D. Rabady, N. Rad, B. Rahbaran, H. Rohringer, J. Schieck1, J. Strauss, W. Waltenberger, C.-E. Wulz1 Institute for Nuclear Problems, Minsk, Belarus V. Chekhovsky, O. Dvornikov, Y. Dydyshka, I. Emeliantchik, A. Litomin, V. Makarenko, V. Mossolov, R. Stefanovitch, J. Suarez Gonzalez, V. Zykunov National Centre for Particle and High Energy Physics, Minsk, Belarus N. Shumeiko Universiteit Antwerpen, Antwerpen, Belgium S. Alderweireldt, E. A. De Wolf, X. Janssen, J. Lauwers, M. Van De Klundert, H. Van Haevermaet, P. Van Mechelen, N. Van Remortel, A. Van Spilbeeck Vrije Universiteit Brussel, Brussel, Belgium S. Abu Zeid, F. Blekman, J. D’Hondt, N. Daci, I. De Bruyn, K. Deroover, S. Lowette, S. Moortgat, L. Moreels, A. Olbrechts, Q. Python, K. Skovpen, S. Tavernier, W. Van Doninck, P. Van Mulders, I. Van Parijs Université Libre de Bruxelles, Bruxelles, Belgium H. Brun, B. Clerbaux, G. De Lentdecker, H. Delannoy, G. Fasanella, L. Favart, R. Goldouzian, A. Grebenyuk, G. Karapostoli, T. Lenzi, A. Léonard, J. Luetic, T. Maerschalk, A. Marinov, A. Randle-conde, T. Seva, C. Vander Velde, P. Vanlaer, D. Vannerom, R. Yonamine, F. Zenoni, F. Zhang2 Ghent University, Ghent, Belgium A. Cimmino, T. Cornelis, D. Dobur, A. Fagot, M. Gul, I. Khvastunov, D. Poyraz, S. Salva, R. Schöfbeck, M. Tytgat, W. Van Driessche, E. Yazgan, N. Zaganidis Université Catholique de Louvain, Louvain-la-Neuve, Belgium H. Bakhshiansohi, C. Beluffi3, O. Bondu, S. Brochet, G. Bruno, A. Caudron, S. De Visscher, C. Delaere, M. Delcourt, B. Francois, A. Giammanco, A. Jafari, M. Komm, G. Krintiras, V. Lemaitre, A. Magitteri, A. Mertens, M. Musich, K. Piotrzkowski, L. Quertenmont, M. Selvaggi, M. Vidal Marono, S. Wertz Centro Brasileiro de Pesquisas Fisicas, Rio de Janeiro, Brazil W. L. Aldá Júnior, F. L. Alves, G. A. Alves, L. Brito, C. Hensel, A. Moraes, M. E. Pol, P. Rebello Teles Universidade do Estado do Rio de Janeiro, Rio de Janeiro, Brazil E. Belchior Batista Das Chagas, W. Carvalho, J. Chinellato4, A. Custódio, E. M. Da Costa, G. G. Da Silveira5, D. De Jesus Damiao, C. De Oliveira Martins, S. Fonseca De Souza, L. M. Huertas Guativa, H. Malbouisson, D. Matos Figueiredo, C. Mora Herrera, L. Mundim, H. Nogima, W. L. Prado Da Silva, A. Santoro, A. Sznajder, E. J. Tonelli Manganote4, A. Vilela Pereira Universidade Estadual Paulistaa , Universidade Federal do ABCb, São Paulo, Brazil S. Ahujaa , C. A. Bernardesa , S. Dograa , T. R. Fernandez Perez Tomeia , E. M. Gregoresb, P. G. Mercadanteb, C. S. Moona , S. F. Novaesa , Sandra S. Padulaa , D. Romero Abadb, J. C. Ruiz Vargasa Institute for Nuclear Research and Nuclear Energy, Sofia, Bulgaria A. Aleksandrov, R. Hadjiiska, P. Iaydjiev, M. Rodozov, S. Stoykova, G. Sultanov, M. Vutova University of Sofia, Sofia, Bulgaria A. Dimitrov, I. Glushkov, L. Litov, B. Pavlov, P. Petkov Universidad de Los Andes, Bogota, Colombia C. Avila, A. Cabrera, L. F. Chaparro Sierra, C. Florez, J. P. Gomez, C. F. González Hernández, J. D. Ruiz Alvarez, J. C. Sanabria University of Split, Faculty of Electrical Engineering, Mechanical Engineering and Naval Architecture, Split, Croatia N. Godinovic, D. Lelas, I. Puljak, P. M. Ribeiro Cipriano, T. Sculac Institute Rudjer Boskovic, Zagreb, Croatia V. Brigljevic, D. Ferencek, K. Kadija, B. Mesic, T. Susa Charles University, Prague, Czech Republic M. Finger8, M. Finger Jr.8 Universidad San Francisco de Quito, Quito, Ecuador E. Carrera Jarrin University of Cyprus, Nicosia, Cyprus A. Attikis, G. Mavromanolakis, J. Mousa, C. Nicolaou, F. Ptochos, P. A. Razis, H. Rykaczewski, D. Tsiakkouri Academy of Scientific Research and Technology of the Arab Republic of Egypt, Egyptian Network of High Energy Physics, Cairo, Egypt E. El-khateeb9, S. Elgammal10, A. Mohamed11 National Institute of Chemical Physics and Biophysics, Tallinn, Estonia M. Kadastik, L. Perrini, M. Raidal, A. Tiko, C. Veelken Helsinki Institute of Physics, Helsinki, Finland J. Härkönen, T. Järvinen, V. Karimäki, R. Kinnunen, T. Lampén, K. Lassila-Perini, S. Lehti, T. Lindén, P. Luukka, J. Tuominiemi, E. Tuovinen, L. Wendland Lappeenranta University of Technology, Lappeenranta, Finland J. Talvitie, T. Tuuva IRFU, CEA, Université Paris-Saclay, Gif-sur-Yvette, France M. Besancon, F. Couderc, M. Dejardin, D. Denegri, B. Fabbro, J. L. Faure, C. Favaro, F. Ferri, S. Ganjour, S. Ghosh, A. Givernaud, P. Gras, G. Hamel de Monchenault, P. Jarry, I. Kucher, E. Locci, M. Machet, J. Malcles, J. Rander, A. Rosowsky, M. Titov Institut Pluridisciplinaire Hubert Curien (IPHC), Université de Strasbourg, CNRS-IN2P3, Strasbourg, France J.-L. Agram12, J. Andrea, A. Aubin, D. Bloch, J.-M. Brom, M. Buttignol, E. C. Chabert, N. Chanon, C. Collard, E. Conte12, X. Coubez, J.-C. Fontaine12, D. Gelé, U. Goerlach, M. Jansová, A.-C. Le Bihan, P. Van Hove Centre de Calcul de l’Institut National de Physique Nucleaire et de Physique des Particules, CNRS/IN2P3, Villeurbanne, France S. Gadrat Université de Lyon, Université Claude Bernard Lyon 1, CNRS-IN2P3, Institut de Physique Nucléaire de Lyon, Villeurbanne, France S. Beauceron, C. Bernet, G. Boudoul, C. A. Carrillo Montoya, R. Chierici, D. Contardo, B. Courbon, P. Depasse, H. El Mamouni, J. Fan, J. Fay, S. Gascon, M. Gouzevitch, G. Grenier, B. Ille, F. Lagarde, I. B. Laktineh, M. Lethuillier, L. Mirabito, A. L. Pequegnot, S. Perries, A. Popov13, D. Sabes, V. Sordini, M. Vander Donckt, P. Verdier, S. Viret Georgian Technical University, Tbilisi, Georgia T. Toriashvili14 Tbilisi State University, Tbilisi, Georgia Z. Tsamalaidze8 RWTH Aachen University, III. Physikalisches Institut A, Aachen, Germany A. Albert, M. Brodski, E. Dietz-Laursonn, D. Duchardt, M. Endres, M. Erdmann, S. Erdweg, T. Esch, R. Fischer, A. Güth, M. Hamer, T. Hebbeker, C. Heidemann, K. Hoepfner, S. Knutzen, M. Merschmeyer, A. Meyer, P. Millet, S. Mukherjee, M. Olschewski, K. Padeken, T. Pook, M. Radziej, H. Reithler, M. Rieger, F. Scheuch, L. Sonnenschein, D. Teyssier, S. Thüer RWTH Aachen University, III. Physikalisches Institut B, Aachen, Germany V. Cherepanov, G. Flügge, B. Kargoll, T. Kress, A. Künsken, J. Lingemann, T. Müller, A. Nehrkorn, A. Nowack, C. Pistone, O. Pooth, A. Stahl15 Institut für Experimentelle Kernphysik, Karlsruhe, Germany M. Akbiyik, C. Barth, S. Baur, C. Baus, J. Berger, E. Butz, R. Caspart, T. Chwalek, F. Colombo, W. De Boer, A. Dierlamm, S. Fink, B. Freund, R. Friese, M. Giffels, A. Gilbert, P. Goldenzweig, D. Haitz, F. Hartmann15, S. M. Heindl, U. Husemann, I. Katkov13, S. Kudella, H. Mildner, M. U. Mozer, Th. Müller, M. Plagge, G. Quast, K. Rabbertz, S. Röcker, F. Roscher, M. Schröder, I. Shvetsov, G. Sieber, H. J. Simonis, R. Ulrich, S. Wayand, M. Weber, T. Weiler, S. Williamson, C. Wöhrmann, R. Wolf Institute of Nuclear and Particle Physics (INPP), NCSR Demokritos, Aghia Paraskevi, Greece G. Anagnostou, G. Daskalakis, T. Geralis, V. A. Giakoumopoulou, A. Kyriakis, D. Loukas, I. Topsis-Giotis National and Kapodistrian University of Athens, Athens, Greece S. Kesisoglou, A. Panagiotou, N. Saoulidou, E. Tziaferi University of Ioánnina, Ioánnina, Greece I. Evangelou, G. Flouris, C. Foudas, P. Kokkas, N. Loukas, N. Manthos, I. Papadopoulos, E. Paradas MTA-ELTE Lendület CMS Particle and Nuclear Physics Group, Eötvös Loránd University, Budapest, Hungary N. Filipovic Wigner Research Centre for Physics, Budapest, Hungary G. Bencze, C. Hajdu, D. Horvath19, F. Sikler, V. Veszpremi, G. Vesztergombi20, A. J. Zsigmond Institute of Nuclear Research ATOMKI, Debrecen, Hungary N. Beni, S. Czellar, J. Karancsi21, A. Makovec, J. Molnar, Z. Szillasi National Institute of Science Education and Research, Bhubaneswar, India S. Bahinipati, S. Choudhury22, P. Mal, K. Mandal, A. Nayak23, D. K. Sahoo, N. Sahoo, S. K. Swain Panjab University, Chandigarh, India S. Bansal, S. B. Beri, V. Bhatnagar, R. Chawla, U. Bhawandeep, A. K. Kalsi, A. Kaur, M. Kaur, R. Kumar, P. Kumari, A. Mehta, M. Mittal, J. B. Singh, G. Walia University of Delhi, Delhi, India Ashok Kumar, A. Bhardwaj, B. C. Choudhary, R. B. Garg, S. Keshri, S. Malhotra, M. Naimuddin, N. Nishu, K. Ranjan, R. Sharma, V. Sharma Saha Institute of Nuclear Physics, Kolkata, India R. Bhattacharya, S. Bhattacharya, K. Chatterjee, S. Dey, S. Dutt, S. Dutta, S. Ghosh, N. Majumdar, A. Modak, K. Mondal, S. Mukhopadhyay, S. Nandan, A. Purohit, A. Roy, D. Roy, S. Roy Chowdhury, S. Sarkar, M. Sharan, S. Thakur Indian Institute of Technology Madras, Madras, India P. K. Behera Bhabha Atomic Research Centre, Mumbai, India R. Chudasama, D. Dutta, V. Jha, V. Kumar, A. K. Mohanty15, P. K. Netrakanti, L. M. Pant, P. Shukla, A. Topkar Tata Institute of Fundamental Research-A, Mumbai, India T. Aziz, S. Dugad, G. Kole, B. Mahakud, S. Mitra, G. B. Mohanty, B. Parida, N. Sur, B. Sutar Tata Institute of Fundamental Research-B, Mumbai, India S. Banerjee, S. Bhowmik24, R. K. Dewanjee, S. Ganguly, M. Guchait, Sa. Jain, S. Kumar, M. Maity24, G. Majumder, K. Mazumdar, T. Sarkar24, N. Wickramage25 Indian Institute of Science Education and Research (IISER), Pune, India S. Chauhan, S. Dube, V. Hegde, A. Kapoor, K. Kothekar, S. Pandey, A. Rane, S. Sharma Institute for Research in Fundamental Sciences (IPM), Tehran, Iran S. Chenarani26, E. Eskandari Tadavani, S. M. Etesami26, M. Khakzad, M. Mohammadi Najafabadi, M. Naseri, S. Paktinat Mehdiabadi27, F. Rezaei Hosseinabadi, B. Safarzadeh28, M. Zeinali University College Dublin, Dublin, Ireland M. Felcini, M. Grunewald INFN Sezione di Baria , Università di Barib, Politecnico di Baric, Bari, Italy M. Abbresciaa ,b, C. Calabriaa ,b, C. Caputoa ,b, A. Colaleoa , D. Creanzaa ,c, L. Cristellaa ,b, N. De Filippisa ,c, M. De Palmaa ,b, L. Fiorea , G. Iasellia ,c, G. Maggia ,c, M. Maggia , G. Minielloa ,b, S. Mya ,b, S. Nuzzoa ,b, A. Pompilia ,b, G. Pugliesea ,c, R. Radognaa ,b, A. Ranieria , G. Selvaggia ,b, A. Sharmaa , L. Silvestrisa ,15, R. Vendittia ,b, P. Verwilligena Chonnam National University, Institute for Universe and Elementary Particles, Kwangju, Korea H. Kim Hanyang University, Seoul, Korea J. A. Brochero Cifuentes, T. J. Kim Korea University, Seoul, Korea S. Cho, S. Choi, Y. Go, D. Gyun, S. Ha, B. Hong, Y. Jo, Y. Kim, K. Lee, K. S. Lee, S. Lee, J. Lim, S. K. Park, Y. Roh Seoul National University, Seoul, Korea J. Almond, J. Kim, H. Lee, S. B. Oh, B. C. Radburn-Smith, S. h. Seo, U. K. Yang, H. D. Yoo, G. B. Yu University of Seoul, Seoul, Korea M. Choi, H. Kim, J. H. Kim, J. S. H. Lee, I. C. Park, G. Ryu, M. S. Ryu Sungkyunkwan University, Suwon, Korea Y. Choi, J. Goh, C. Hwang, J. Lee, I. Yu Vilnius University, Vilnius, Lithuania V. Dudenas, A. Juodagalvis, J. Vaitkus National Centre for Particle Physics, Universiti Malaya, Kuala Lumpur, Malaysia I. Ahmed, Z. A. Ibrahim, J. R. Komaragiri, M. A. B. Md Ali32, F. Mohamad Idris33, W. A. T. Wan Abdullah, M. N. Yusli, Z. Zolkapli Centro de Investigacion y de Estudios Avanzados del IPN, Mexico City, Mexico H. Castilla-Valdez, E. De La Cruz-Burelo, I. Heredia-De La Cruz34, A. Hernandez-Almada, R. Lopez-Fernandez, R. Magaña Villalba, J. Mejia Guisao, A. Sanchez-Hernandez Universidad Iberoamericana, Mexico City, Mexico S. Carrillo Moreno, C. Oropeza Barrera, F. Vazquez Valencia Benemerita Universidad Autonoma de Puebla, Puebla, Mexico S. Carpinteyro, I. Pedraza, H. A. Salazar Ibarguen, C. Uribe Estrada Universidad Autónoma de San Luis Potosí, San Luis Potosí, Mexico A. Morelos Pineda National Centre for Physics, Quaid-I-Azam University, Islamabad, Pakistan A. Ahmad, M. Ahmad, Q. Hassan, H. R. Hoorani, W. A. Khan, A. Saddique, M. A. Shah, M. Shoaib, M. Waqas National Centre for Nuclear Research, Swierk, Poland H. Bialkowska, M. Bluj, B. Boimska, T. Frueboes, M. Górski, M. Kazana, K. Nawrocki, K. Romanowska-Rybinska, M. Szleper, P. Zalewski Laboratório de Instrumentação e Física Experimental de Partículas, Lisbon, Portugal P. Bargassa, C. Beirão Da Cruz E. Silva, B. Calpas, A. Di Francesco, P. Faccioli, P. G. Ferreira Parracho, M. Gallinaro, J. Hollar, N. Leonardo, L. Lloret Iglesias, M. V. Nemallapudi, J. Rodrigues Antunes, J. Seixas, O. Toldaiev, D. Vadruccio, J. Varela, P. Vischia Joint Institute for Nuclear Research, Dubna, Russia V. Alexakhin, P. Bunin, M. Gavrilenko, I. Golutvin, I. Gorbunov, A. Kamenev, V. Karjavin, A. Lanev, A. Malakhov, V. Matveev36,37, V. Palichik, V. Perelygin, M. Savina, S. Shmatov, N. Skatchkov, V. Smirnov, N. Voytishin, A. Zarubin Petersburg Nuclear Physics Institute, Gatchina, St. Petersburg, Russia L. Chtchipounov, V. Golovtsov, Y. Ivanov, V. Kim38, E. Kuznetsova39, V. Murzin, V. Oreshkin, V. Sulimov, A. Vorobyev Institute for Nuclear Research, Moscow, Russia Yu. Andreev, A. Dermenev, S. Gninenko, N. Golubev, A. Karneyeu, M. Kirsanov, N. Krasnikov, A. Pashenkov, D. Tlisov, A. Toropin Institute for Theoretical and Experimental Physics, Moscow, Russia V. Epshteyn, V. Gavrilov, N. Lychkovskaya, V. Popov, I. Pozdnyakov, G. Safronov, A. Spiridonov, M. Toms, E. Vlasov, A. Zhokin Moscow Institute of Physics and Technology, Moscow, Russia A. Bylinkin37 National Research Nuclear University ‘Moscow Engineering Physics Institute’ (MEPhI), Moscow, Russia O. Markin, E. Popova, E. Zhemchugov P.N. Lebedev Physical Institute, Moscow, Russia V. Andreev, M. Azarkin37, I. Dremin37, M. Kirakosyan, A. Leonidov37, A. Terkulov Skobeltsyn Institute of Nuclear Physics, Lomonosov Moscow State University, Moscow, Russia A. Baskakov, A. Belyaev, E. Boos, V. Bunichev, M. Dubinin40, L. Dudko, A. Ershov, A. Gribushin, V. Klyukhin, O. Kodolova, I. Lokhtin, I. Miagkov, S. Obraztsov, S. Petrushanko, V. Savrin Novosibirsk State University (NSU), Novosibirsk, Russia V. Blinov41, Y. Skovpen41, D. Shtol41 State Research Center of Russian Federation, Institute for High Energy Physics, Protvino, Russia I. Azhgirey, I. Bayshev, S. Bitioukov, D. Elumakhov, V. Kachanov, A. Kalinin, D. Konstantinov, V. Krychkine, V. Petrov, R. Ryutin, A. Sobol, S. Troshin, N. Tyurin, A. Uzunian, A. Volkov University of Belgrade, Faculty of Physics and Vinca Institute of Nuclear Sciences, Belgrade, Serbia P. Adzic42, P. Cirkovic, D. Devetak, M. Dordevic, J. Milosevic, V. Rekovic Centro de Investigaciones Energéticas Medioambientales y Tecnológicas (CIEMAT), Madrid, Spain J. Alcaraz Maestre, M. Barrio Luna, E. Calvo, M. Cerrada, M. Chamizo Llatas, N. Colino, B. De La Cruz, A. Delgado Peris, A. Escalante Del Valle, C. Fernandez Bedoya, J. P. Fernández Ramos, J. Flix, M. C. Fouz, P. Garcia-Abia, O. Gonzalez Lopez, S. Goy Lopez, J. M. Hernandez, M. I. Josa, E. Navarro De Martino, A. Pérez-Calero Yzquierdo, J. Puerta Pelayo, A. Quintario Olmeda, I. Redondo, L. Romero, M. S. Soares Universidad Autónoma de Madrid, Madrid, Spain J. F. de Trocóniz, M. Missiroli, D. Moran Universidad de Oviedo, Oviedo, Spain J. Cuevas, J. Fernandez Menendez, I. Gonzalez Caballero, J. R. González Fernández, E. Palencia Cortezon, S. Sanchez Cruz, I. Suárez Andrés, J. M. Vizan Garcia Institute for Scintillation Materials of National Academy of Science of Ukraine, Kharkov, Ukraine B. Grynyov National Scientific Center, Kharkov Institute of Physics and Technology, Kharkov, Ukraine L. Levchuk, P. Sorokin University of Bristol, Bristol, UK R. Aggleton, F. Ball, L. Beck, J. J. Brooke, D. Burns, E. Clement, D. Cussans, H. Flacher, J. Goldstein, M. Grimes, G. P. Heath, H. F. Heath, J. Jacob, L. Kreczko, C. Lucas, D. M. Newbold63, S. Paramesvaran, A. Poll, T. Sakuma, S. Seif El Nasr-storey, D. Smith, V. J. Smith Rutherford Appleton Laboratory, Didcot, UK K. W. Bell, A. Belyaev64, C. Brew, R. M. Brown, L. Calligaris, D. Cieri, D. J. A. Cockerill, J. A. Coughlan, K. Harder, S. Harper, E. Olaiya, D. Petyt, C. H. Shepherd-Themistocleous, A. Thea, I. R. Tomalin, T. Williams Baylor University, Waco, USA A. Borzou, K. Call, J. Dittmann, K. Hatakeyama, H. Liu, N. Pastika Catholic University of America, Washington, DC, USA R. Bartek, A. Dominguez The University of Alabama, Tuscaloosa, USA S. I. Cooper, C. Henderson, P. Rumerio, C. West Boston University, Boston, USA D. Arcaro, A. Avetisyan, T. Bose, D. Gastler, D. Rankin, C. Richardson, J. Rohlf, L. Sulak, D. Zou Brown University, Providence, USA G. Benelli, D. Cutts, A. Garabedian, J. Hakala, U. Heintz, J. M. Hogan, O. Jesus, K. H. M. Kwok, E. Laird, G. Landsberg, Z. Mao, M. Narain, S. Piperov, S. Sagir, E. Spencer, R. Syarif University of California, Davis, Davis, USA R. Breedon, D. Burns, M. Calderon De La Barca Sanchez, S. Chauhan, M. Chertok, J. Conway, R. Conway, P. T. Cox, R. Erbacher, C. Flores, G. Funk, M. Gardner, W. Ko, R. Lander, C. Mclean, M. Mulhearn, D. Pellett, J. Pilot, S. Shalhout, M. Shi, J. Smith, M. Squires, D. Stolp, K. Tos, M. Tripathi University of California, Los Angeles, USA C. Bravo, R. Cousins, A. Dasgupta, A. Florent, J. Hauser, M. Ignatenko, N. Mccoll, D. Saltzberg, C. Schnaible, V. Valuev, M. Weber University of California, Riverside, Riverside, USA E. Bouvier, K. Burt, R. Clare, J. Ellison, J. W. Gary, S. M. A. Ghiasi Shirazi, G. Hanson, J. Heilman, P. Jandir, E. Kennedy, F. Lacroix, O. R. Long, M. Olmedo Negrete, M. I. Paneva, A. Shrinivas, W. Si, H. Wei, S. Wimpenny, B. R. Yates University of California, San Diego, La Jolla, USA J. G. Branson, G. B. Cerati, S. Cittolin, M. Derdzinski, R. Gerosa, A. Holzner, D. Klein, V. Krutelyov, J. Letts, I. Macneill, D. Olivito, S. Padhi, M. Pieri, M. Sani, V. Sharma, S. Simon, M. Tadel, A. Vartak, S. Wasserbaech66, C. Welke, J. Wood, F. Würthwein, A. Yagil, G. Zevi Della Porta University of California, Santa Barbara-Department of Physics, Santa Barbara, USA N. Amin, R. Bhandari, J. Bradmiller-Feld, C. Campagnari, A. Dishaw, V. Dutta, M. Franco Sevilla, C. George, F. Golf, L. Gouskos, J. Gran, R. Heller, J. Incandela, S. D. Mullin, A. Ovcharova, H. Qu, J. Richman, D. Stuart, I. Suarez, J. Yoo California Institute of Technology, Pasadena, USA D. Anderson, J. Bendavid, A. Bornheim, J. Bunn, J. Duarte, J. M. Lawhorn, A. Mott, H. B. Newman, C. Pena, M. Spiropulu, J. R. Vlimant, S. Xie, R. Y. Zhu Carnegie Mellon University, Pittsburgh, USA M. B. Andrews, T. Ferguson, M. Paulini, J. Russ, M. Sun, H. Vogel, I. Vorobiev, M. Weinberg University of Colorado Boulder, Boulder, USA J. P. Cumalat, W. T. Ford, F. Jensen, A. Johnson, M. Krohn, T. Mulholland, K. Stenson, S. R. Wagner Cornell University, Ithaca, USA J. Alexander, J. Chaves, J. Chu, S. Dittmer, K. Mcdermott, N. Mirman, G. Nicolas Kaufman, J. R. Patterson, A. Rinkevicius, A. Ryd, L. Skinnari, L. Soffi, S. M. Tan, Z. Tao, J. Thom, J. Tucker, P. Wittich, M. Zientek Fermi National Accelerator Laboratory, Batavia, USA S. Abdullin, M. Albrow, G. Apollinari, A. Apresyan, S. Banerjee, L. A. T. Bauerdick, A. Beretvas, J. Berryhill, P. C. Bhat, G. Bolla, K. Burkett, J. N. Butler, H. W. K. Cheung, F. Chlebana, S. Cihangir†, M. Cremonesi, V. D. Elvira, I. Fisk, J. Freeman, E. Gottschalk, L. Gray, D. Green, S. Grünendahl, O. Gutsche, D. Hare, R. M. Harris, S. Hasegawa, J. Hirschauer, Z. Hu, B. Jayatilaka, S. Jindariani, M. Johnson, U. Joshi, B. Klima, B. Kreis, S. Lammel, J. Linacre, D. Lincoln, R. Lipton, M. Liu, T. Liu, R. Lopes De Sá, J. Lykken, K. Maeshima, N. Magini, J. M. Marraffino, S. Maruyama, D. Mason, P. McBride, P. Merkel, S. Mrenna, S. Nahn, V. O’Dell, K. Pedro, O. Prokofyev, G. Rakness, L. Ristori, E. Sexton-Kennedy, A. Soha, W. J. Spalding, L. Spiegel, S. Stoynev, J. Strait, N. Strobbe, L. Taylor, S. Tkaczyk, N. V. Tran, L. Uplegger, E. W. Vaandering, C. Vernieri, M. Verzocchi, R. Vidal, M. Wang, H. A. Weber, A. Whitbeck, Y. Wu University of Florida, Gainesville, USA D. Acosta, P. Avery, P. Bortignon, D. Bourilkov, A. Brinkerhoff, A. Carnes, M. Carver, D. Curry, S. Das, R. D. Field, I. K. Furic, J. Konigsberg, A. Korytov, J. F. Low, P. Ma, K. Matchev, H. Mei, G. Mitselmakher, D. Rank, L. Shchutska, D. Sperka, L. Thomas, J. Wang, S. Wang, J. Yelton Florida International University, Miami, USA S. Linn, P. Markowitz, G. Martinez, J. L. Rodriguez Florida State University, Tallahassee, USA A. Ackert, T. Adams, A. Askew, S. Bein, S. Hagopian, V. Hagopian, K. F. Johnson, H. Prosper, A. Santra, R. Yohay Florida Institute of Technology, Melbourne, USA M. M. Baarmand, V. Bhopatkar, S. Colafranceschi, M. Hohlmann, D. Noonan, T. Roy, F. Yumiceva University of Illinois at Chicago (UIC), Chicago, USA M. R. Adams, L. Apanasevich, D. Berry, R. R. Betts, I. Bucinskaite, R. Cavanaugh, O. Evdokimov, L. Gauthier, C. E. Gerber, D. J. Hofman, K. Jung, I. D. Sandoval Gonzalez, N. Varelas, H. Wang, Z. Wu, M. Zakaria, J. Zhang The University of Iowa, Iowa City, USA B. Bilki67, W. Clarida, K. Dilsiz, S. Durgut, R. P. Gandrajula, M. Haytmyradov, V. Khristenko, J.-P. Merlo, H. Mermerkaya68, A. Mestvirishvili, A. Moeller, J. Nachtman, H. Ogul, Y. Onel, F. Ozok69, A. Penzo, C. Snyder, E. Tiras, J. Wetzel, K. Yi Johns Hopkins University, Baltimore, USA I. Anderson, B. Blumenfeld, A. Cocoros, N. Eminizer, D. Fehling, L. Feng, A. V. Gritsan, P. Maksimovic, M. Osherson, J. Roskes, U. Sarica, M. Swartz, M. Xiao, Y. Xin, C. You The University of Kansas, Lawrence, USA A. Al-bataineh, P. Baringer, A. Bean, S. Boren, J. Bowen, J. Castle, L. Forthomme, R. P. Kenny III, S. Khalil, A. Kropivnitskaya, D. Majumder, W. Mcbrayer, M. Murray, S. Sanders, R. Stringer, J. D. Tapia Takaki, Q. Wang Kansas State University, Manhattan, USA A. Ivanov, K. Kaadze, Y. Maravin, A. Mohammadi, L. K. Saini, N. Skhirtladze, S. Toda Lawrence Livermore National Laboratory, Livermore, USA F. Rebassoo, D. Wright University of Maryland, College Park, USA C. Anelli, A. Baden, O. Baron, A. Belloni, B. Calvert, S. C. Eno, C. Ferraioli, J. A. Gomez, N. J. Hadley, S. Jabeen, R. G. Kellogg, T. Kolberg, J. Kunkle, Y. Lu, A. C. Mignerey, F. Ricci-Tam, Y. H. Shin, A. Skuja, M. B. Tonjes, S. C. Tonwar Massachusetts Institute of Technology, Cambridge, USA D. Abercrombie, B. Allen, A. Apyan, V. Azzolini, R. Barbieri, A. Baty, R. Bi, K. Bierwagen, S. Brandt, W. Busza, I. A. Cali, M. D’Alfonso, Z. Demiragli, L. Di Matteo, G. Gomez Ceballos, M. Goncharov, D. Hsu, Y. Iiyama, G. M. Innocenti, M. Klute, D. Kovalskyi, K. Krajczar, Y. S. Lai, Y.-J. Lee, A. Levin, P. D. Luckey, B. Maier, A. C. Marini, C. Mcginn, C. Mironov, S. Narayanan, X. Niu, C. Paus, C. Roland, G. Roland, J. Salfeld-Nebgen, G. S. F. Stephans, K. Tatar, M. Varma, D. Velicanu, J. Veverka, J. Wang, T. W. Wang, B. Wyslouch, M. Yang University of Minnesota, Minneapolis, USA A. C. Benvenuti, R. M. Chatterjee, A. Evans, P. Hansen, S. Kalafut, S. C. Kao, Y. Kubota, Z. Lesko, J. Mans, S. Nourbakhsh, N. Ruckstuhl, R. Rusack, N. Tambe, J. Turkewitz University of Mississippi, Oxford, USA J. G. Acosta, S. Oliveros University of Nebraska-Lincoln, Lincoln, USA E. Avdeeva, K. Bloom, D. R. Claes, C. Fangmeier, R. Gonzalez Suarez, R. Kamalieddin, I. Kravchenko, A. Malta Rodrigues, F. Meier, J. Monroy, J. E. Siado, G. R. Snow, B. Stieger State University of New York at Buffalo, Buffalo, USA M. Alyari, J. Dolen, A. Godshalk, C. Harrington, I. Iashvili, J. Kaisen, A. Kharchilava, A. Parker, S. Rappoccio, B. Roozbahani Northeastern University, Boston, USA G. Alverson, E. Barberis, A. Hortiangtham, A. Massironi, D. M. Morse, D. Nash, T. Orimoto, R. Teixeira De Lima, D. Trocino, R.-J. Wang, D. Wood Northwestern University, Evanston, USA S. Bhattacharya, O. Charaf, K. A. Hahn, A. Kumar, N. Mucia, N. Odell, B. Pollack, M. H. Schmitt, K. Sung, M. Trovato, M. Velasco University of Notre Dame, Notre Dame, USA N. Dev, M. Hildreth, K. Hurtado Anampa, C. Jessop, D. J. Karmgard, N. Kellams, K. Lannon, N. Marinelli, F. Meng, C. Mueller, Y. Musienko36, M. Planer, A. Reinsvold, R. Ruchti, G. Smith, S. Taroni, M. Wayne, M. Wolf, A. Woodard The Ohio State University, Columbus, USA J. Alimena, L. Antonelli, B. Bylsma, L. S. Durkin, S. Flowers, B. Francis, A. Hart, C. Hill, R. Hughes, W. Ji, B. Liu, W. Luo, D. Puigh, B. L. Winer, H. W. Wulsin Princeton University, Princeton, USA S. Cooperstein, O. Driga, P. Elmer, J. Hardenbrook, P. Hebda, D. Lange, J. Luo, D. Marlow, T. Medvedeva, K. Mei, J. Olsen, C. Palmer, P. Piroué, D. Stickland, A. Svyatkovskiy, C. Tully Purdue University Calumet, Hammond, USA N. Parashar, J. Stupak Rice University, Houston, USA A. Adair, B. Akgun, Z. Chen, K. M. Ecklund, F. J. M. Geurts, M. Guilbaud, W. Li, B. Michlin, M. Northup, B. P. Padley, J. Roberts, J. Rorie, Z. Tu, J. Zabel University of Rochester, Rochester, USA B. Betchart, A. Bodek, P. de Barbaro, R. Demina, Y. t. Duh, T. Ferbel, M. Galanti, A. Garcia-Bellido, J. Han, O. Hindrichs, A. Khukhunaishvili, K. H. Lo, P. Tan, M. Verzetti Rutgers, The State University of New Jersey, Piscataway, USA A. Agapitos, J. P. Chou, Y. Gershtein, T. A. Gómez Espinosa, E. Halkiadakis, M. Heindl, E. Hughes, S. Kaplan, R. Kunnawalkam Elayavalli, S. Kyriacou, A. Lath, K. Nash, H. Saka, S. Salur, S. Schnetzer, D. Sheffield, S. Somalwar, R. Stone, S. Thomas, P. Thomassen, M. Walker University of Tennessee, Knoxville, USA A. G. Delannoy, M. Foerster, J. Heideman, G. Riley, K. Rose, S. Spanier, K. Thapa Texas A & M University, College Station, USA O. Bouhali70, A. Celik, M. Dalchenko, M. De Mattia, A. Delgado, S. Dildick, R. Eusebi, J. Gilmore, T. Huang, E. Juska, T. Kamon71, R. Mueller, Y. Pakhotin, R. Patel, A. Perloff, L. Perniè, D. Rathjens, A. Safonov, A. Tatarinov, K. A. Ulmer Texas Tech University, Lubbock, USA N. Akchurin, C. Cowden, J. Damgov, F. De Guio, C. Dragoiu, P. R. Dudero, J. Faulkner, E. Gurpinar, S. Kunori, K. Lamichhane, S. W. Lee, T. Libeiro, T. Peltola, S. Undleeb, I. Volobouev, Z. Wang Vanderbilt University, Nashville, USA S. Greene, A. Gurrola, R. Janjam, W. Johns, C. Maguire, A. Melo, H. Ni, P. Sheldon, S. Tuo, J. Velkovska, Q. Xu University of Virginia, Charlottesville, USA M. W. Arenton, P. Barria, B. Cox, J. Goodell, R. Hirosky, A. Ledovskoy, H. Li, C. Neu, T. Sinthuprasith, X. Sun, Y. Wang, E. Wolfe, F. Xia Wayne State University, Detroit, USA C. Clarke, R. Harr, P. E. Karchin, J. Sturdy University of Wisconsin-Madison, Madison, WI, USA D. A. Belknap, J. Buchanan, C. Caillol, S. Dasu, L. Dodd, S. Duric, B. Gomber, M. Grothe, M. Herndon, A. Hervé, P. Klabbers, A. Lanaro, A. Levine, K. Long, R. Loveless, I. Ojalvo, T. Perry, G. A. Pierro, G. Polese, T. Ruggles, A. Savin, N. Smith, W. H. Smith, D. Taylor, N. Woods † Deceased 12: Also at Université de Haute Alsace, Mulhouse, France 13: Also at Skobeltsyn Institute of Nuclear Physics, Lomonosov Moscow State University, Moscow, Russia 14: Also at Tbilisi State University, Tbilisi, Georgia 15: Also at CERN, European Organization for Nuclear Research, Geneva, Switzerland 16: Also at RWTH Aachen University, III. Physikalisches Institut A, Aachen, Germany 17: Also at University of Hamburg, Hamburg, Germany 18: Also at Brandenburg University of Technology, Cottbus, Germany 19: Also at Institute of Nuclear Research ATOMKI, Debrecen, Hungary 20: Also at MTA-ELTE Lendület CMS Particle and Nuclear Physics Group, Eötvös Loránd University, Budapest, Hungary 21: Also at Institute of Physics, University of Debrecen, Debrecen, Hungary 22: Also at Indian Institute of Science Education and Research, Bhopal, India 23: Also at Institute of Physics, Bhubaneswar, India 24: Also at University of Visva-Bharati, Santiniketan, India 25: Also at University of Ruhuna, Matara, Sri Lanka 26: Also at Isfahan University of Technology, Isfahan, Iran 27: Also at Yazd University, Yazd, Iran 28: Also at Plasma Physics Research Center, Science and Research Branch, Islamic Azad University, Tehran, Iran 29: Also at Università degli Studi di Siena, Siena, Italy 30: Also at Laboratori Nazionali di Legnaro dell’INFN, Legnaro, Italy 31: Also at Purdue University, West Lafayette, USA 32: Also at International Islamic University of Malaysia, Kuala Lumpur, Malaysia 33: Also at Malaysian Nuclear Agency, MOSTI, Kajang, Malaysia 34: Also at Consejo Nacional de Ciencia y Tecnología, Mexico city, Mexico 35: Also at Warsaw University of Technology, Institute of Electronic Systems, Warsaw, Poland 36: Also at Institute for Nuclear Research, Moscow, Russia 37: Now at National Research Nuclear University ’Moscow Engineering Physics Institute’ (MEPhI), Moscow, Russia 38: Also at St. Petersburg State Polytechnical University, St. Petersburg, Russia 39: Also at University of Florida, Gainesville, USA 40: Also at California Institute of Technology, Pasadena, USA 41: Also at Budker Institute of Nuclear Physics, Novosibirsk, Russia 42: Also at Faculty of Physics, University of Belgrade, Belgrade, Serbia 43: Also at INFN Sezione di Roma; Università di Roma, Rome, Italy 44: Also at University of Belgrade, Faculty of Physics and Vinca Institute of Nuclear Sciences, Belgrade, Serbia 45: Also at Scuola Normale e Sezione dell’INFN, Pisa, Italy 46: Also at National and Kapodistrian University of Athens, Athens, Greece 47: Also at Riga Technical University, Riga, Latvia 48: Also at Institute for Theoretical and Experimental Physics, Moscow, Russia 49: Also at Albert Einstein Center for Fundamental Physics, Bern, Switzerland 50: Also at Gaziosmanpasa University, Tokat, Turkey 51: Also at Adiyaman University, Adiyaman, Turkey 52: Also at Istanbul Aydin University, Istanbul, Turkey 53: Also at Mersin University, Mersin, Turkey 54: Also at Cag University, Mersin, Turkey 55: Also at Piri Reis University, Istanbul, Turkey 56: Also at Ozyegin University, Istanbul, Turkey 57: Also at Izmir Institute of Technology, Izmir, Turkey 58: Also at Marmara University, Istanbul, Turkey 59: Also at Kafkas University, Kars, Turkey 60: Also at Istanbul Bilgi University, Istanbul, Turkey 61: Also at Yildiz Technical University, Istanbul, Turkey 62: Also at Hacettepe University, Ankara, Turkey 63: Also at Rutherford Appleton Laboratory, Didcot, UK 64: Also at School of Physics and Astronomy, University of Southampton, Southampton, UK 65: Also at Instituto de Astrofísica de Canarias, La Laguna, Spain 66: Also at Utah Valley University, Orem, USA 67: Also at Argonne National Laboratory, Argonne, USA 68: Also at Erzincan University, Erzincan, Turkey 69: Also at Mimar Sinan University, Istanbul, Istanbul, Turkey 70: Also at Texas A&M University at Qatar, Doha, Qatar 71: Also at Kyungpook National University, Taegu, Korea 1. G. t 'Hooft, Naturalness, chiral symmetry, and spontaneous chiral symmetry breaking . NATO Sci. Ser. B 59 , 135 ( 1980 ) 2. E. Witten , Dynamical breaking of supersymmetry . Nucl. Phys. B 188 , 513 ( 1981 ). doi:10.1016/ 0550 -3213(81)90006- 7 3. M. Dine , W. Fischler , M. Srednicki , Supersymmetric technicolor. Nucl. Phys. B 189 , 575 ( 1981 ). doi:10.1016/ 0550 -3213(81)90582- 4 4. S. Dimopoulos , S. Raby , Supercolor. Nucl. Phys . B 192 , 353 ( 1981 ). doi:10.1016/ 0550 -3 2 5. S. Dimopoulos , H. Georgi , Softly broken supersymmetry and SU ( 5 ). Nucl. Phys . B 193 , 150 ( 1981 ). doi:10.1016/ 0550 -3213( 8 6. R.K. Kaul , P. Majumdar , Cancellation of quadratically divergent mass corrections in globally supersymmetric spontaneously broken gauge theories . Nucl. Phys. B 199 , 36 ( 1982 ). doi:10.1016/ 0550 -3213(82) 90565 - X 7. J. Wess , B. Zumino , Supergauge transformations in fourdimensions . Nucl. Phys. B 70 , 39 ( 1974 ). doi:10.1016/ 0550 -32 1 8. G.R. Farrar , P. Fayet , Phenomenology of the production, decay, and detection of new hadronic states associated with supersymmetry . Phys. Lett. B 76 , 575 ( 1978 ). doi:10.1016/ 0370 -2693(78)90858- 4 9. C. Boehm , A. Djouadi , M. Drees , Light scalar top quarks and supersymmetric dark matter . Phys. Rev. D 62 , 035012 ( 2000 ). doi:10. 1103/PhysRevD.62.035012. arXiv:hep-ph/9911496 10. C. Balázs , M. Carena , C.E.M. Wagner , Dark matter, light stops and electroweak baryogenesis . Phys. Rev. D 70 , 015007 ( 2004 ). doi:10. 1103/PhysRevD.70.015007. arXiv:hep-ph/403224 11. G. Jungman , M. Kamionkowski , K. Griest , Supersymmetric dark matter . Phys. Rept . 267 , 195 ( 1996 ). doi:10.1016/ 0370-1573(95)00058-5. arXiv:hep-ph/9506380 12. ATLAS Collaboration, Search for a supersymmetric partner to the top quark in final states with jets and missing transverse momentum at √s = 7 TeV with the ATLAS detector . Phys. Rev. Lett . 109 , 211802 ( 2012 ). doi:10.1103/PhysRevLett.109.211802. arXiv:1208.1447 13. ATLAS Collaboration, Search for direct top squark pair production in final states with one isolated lepton, jets, and missing transverse momentum in √s = 7 TeV pp collisions using 4.7 fb−1 of ATLAS data . Phys. Rev. Lett . 109 , 211803 ( 2012 ). doi:10.1103/ PhysRevLett.109.211803. arXiv:1208.2590 14. ATLAS Collaboration, Search for a heavy top-quark partner in final states with two leptons with the ATLAS detector at the LHC . JHEP 11 , 094 ( 2012 ). doi:10.1007/ JHEP11(2012)094 . arXiv:1209.4186 15. ATLAS Collaboration, Search for direct top-squark pair production in final states with two leptons in pp collisions at √s = 8 TeV with the ATLAS detector . JHEP 06 , 124 ( 2014 ). doi:10.1007/ JHEP06(2014)124 . arXiv:1403.4853 16. ATLAS Collaboration, Search for direct third-generation squark pair production in final states with missing transverse momentum and two b-jets in √s = 8 TeV pp collisions with the ATLAS detector . JHEP 10 , 189 ( 2013 ). doi:10.1007/ JHEP10(2013)189 . arXiv:1308.2631 17. ATLAS Collaboration, Measurement of spin correlation in topantitop quark events and search for top squark pair production in pp collisions at √s = 8 TeV using the ATLAS detector . Phys. Rev. Lett . 114 , 142001 ( 2015 ). doi:10.1103/PhysRevLett.114.142001. arXiv:1412.4742 18. ATLAS Collaboration, Search for pair-produced third-generation squarks decaying via charm quarks or in compressed supersymmetric scenarios in pp collisions at √s = 8 TeV with the ATLAS detector . Phys. Rev. D 90 , 052008 ( 2014 ). doi:10.1103/PhysRevD. 90.052008. arXiv:1407.0608 19. ATLAS Collaboration, Search for direct pair production of the top squark in all-hadronic final states in proton-proton collisions at √s = 8 TeV with the ATLAS detector . JHEP 09 , 015 ( 2014 ). doi:10.1007/ JHEP09(2014)015 . arXiv:1406.1122 20. ATLAS Collaboration, ATLAS Run 1 searches for direct pair production of third-generation squarks at the large hadron collider . Eur. Phys. J. C 75 , 510 ( 2015 ). doi:10.1140/epjc/s10052- 015 - 3726 -9. arXiv: 1506 .08616 21. CMS Collaboration, Search for top-squark pair production in the single-lepton final state in pp collisions at √s = 8 TeV . Eur. Phys. J. C 73 , 2677 ( 2013 ). doi:10.1140/epjc/s10052- 013 - 2677 -2. arXiv: 1308 .1586 22. CMS Collaboration, Search for supersymmetry in hadronic final states with missing transverse energy using the variables αT and b-quark multiplicity in pp collisions at √s = 8 T eV . Eur . Phys. J. C 73 , 2568 ( 2013 ). doi:10.1140/epjc/s10052- 013 - 2568 -6. arXiv: 1303 .2985 23. CMS Collaboration, Search for supersymmetry using razor variables in events with b-tagged jets in pp collisions at √s = 8 TeV. Phys. Rev. D 91 , 052018 ( 2015 ). doi:10.1103/PhysRevD.91. 052018. arXiv:1502.00300 24. CMS Collaboration, Searches for third-generation squark production in fully hadronic final states in proton-proton collisions at √s = 8 TeV . JHEP 06 , 116 ( 2015 ). doi:10.1007/ JHEP06(2015)116 . arXiv:1503.08037 25. CMS Collaboration, Search for direct pair production of supersymmetric top quarks decaying to all-hadronic final states in pp collisions at √s = 8 TeV . Eur. Phys. J. C 76 , 460 ( 2016 ). doi:10. 1140/epjc/s10052- 016 - 4292 -5. arXiv: 1603 .00765 26. D0 Collaboration, Search for 3- and 4-body decays of the scalar top quark in pp¯ collisions at √s = 1.8 TeV . Phys. Lett . B 581 , 147 ( 2004 ). doi:10.1016/j.physletb. 2003 .12.001 27. D0 Collaboration, Search for pair production of the scalar top quark in muon+tau final states . Phys. Lett. B 710 , 578 ( 2012 ). doi:10. 1016/j.physletb. 2012 .03.028. arXiv:1202.1978 28. D0 Collaboration, Search for the lightest scalar top quark in events with two leptons in pp¯ collisions at √s = 1.96 TeV . Phys. Lett . B 659 , 500 ( 2008 ). doi:10.1016/j.physletb. 2007 .11.086. arXiv:0707.2864 29. CDF Collaboration, Search for the supersymmetric partner of the top quark in pp¯ collisions at √ s = 1.96 TeV. Phys. Rev. D 82 , 092001 ( 2010 ). doi:10.1103/PhysRevD.82.092001. arXiv:1009.0266 30. CDF Collaboration, Search for the supersymmetric partner of the top quark in dilepton events from pp¯ collisions at √ s = 1.8 TeV. Phys. Rev. Lett . 90 , 251801 ( 2003 ). doi:10.1103/PhysRevLett.90. 251801. arXiv:hep-ex/0302009 31. C. Borschensky et al., Squark and gluino production cross sections in pp collisions at √s = 13, 14, 33 and 100 TeV. Eur. Phys. J. C 74 , 3174 ( 2014 ). doi:10.1140/epjc/s10052- 014 - 3174 -y. arXiv:1407.5066 32. LHC SUSY Cross Section Working Group, SUSY cross sections twiki ( 2015 ). https://twiki.cern.ch/twiki/bin/view/LHCPhysics/ SUSYCrossSections 33. CMS Collaboration, Search for new physics in the all-hadronic final state with the MT 2 variable . JHEP 10 , 006 ( 2016 ). doi:10. 1007/ JHEP10(2016)006 . arXiv:1603.04053 34. CMS Collaboration, Inclusive search for supersymmetry using razor variables in pp collisions at √s = 13 TeV. Phys. Rev. D 95 , 012003 ( 2017 ). doi:10.1103/PhysRevD.95.012003. arXiv:1609.07658 35. CMS Collaboration, A search for new phenomena in pp collisions at √s = 13 TeV in final states with missing transverse momentum and at least one jet using the αT variable . Eur. Phys. J. C 77 , 294 ( 2017 ). doi:10.1140/epjc/s10052- 017 - 4787 -8. arXiv: 1611 .00338 36. ATLAS Collaboration, Search for top squarks in final states with one isolated lepton, jets, and missing transverse momentum in √s = 13 TeV collisions with the ATLAS detector . Phys. Rev. D 94 , 052009 ( 2016 ). doi:10.1103/PhysRevD.94.052009. arXiv:1606.03903 37. CMS Collaboration, The CMS experiment at the CERN LHC . JINST 3, S08004 ( 2008 ). doi:10.1088/ 1748 - 0221 /3/08/S08004 38. CMS Collaboration, Particle-flow event reconstruction in CMS and performance for jets , taus, and E miss. CMS Physics AnalyT sis Summary CMS-PAS-PFT-09-001 ( 2009 ). http://cdsweb.cern. ch/record/1194487CMS 39. CMS Collaboration, Commissioning of the particle-flow event reconstruction with the first LHC collisions recorded in the CMS detector . CMS Physics Analysis Summary CMS-PAS-PFT-10-001 ( 2010 ). http://cdsweb.cern. ch/record/1247373 40. M. Cacciari , G.P. Salam , G. Soyez , The anti-kt jet clustering algorithm . JHEP 04 , 063 ( 2008 ). doi:10.1088/ 1126 - 6708 /2008/04/063. arXiv: 0802 .1189 41. M. Cacciari , G.P. Salam , G. Soyez , FastJet user manual . Eur. Phys. J. C 72 , 1896 ( 2012 ). doi:10.1140/epjc/s10052- 012 - 1896 -2. arXiv: 1111 .6097 42. CMS Collaboration, Study of pileup removal algorithms for jets . CMS Physics Analysis Summary CMS-PAS-JME-14-001 ( 2014 ). https://cds.cern. ch/record/1751454 43. CMS Collaboration, Identification of b-quark jets with the CMS experiment . JINST 8, P04013 ( 2013 ). doi:10.1088/ 1748 - 0221 /8/ 04/P04013. arXiv: 1211 .4462 44. CMS Collaboration, Identification of b quark jets at the CMS experiment in the LHC Run 2 . CMS Physics Analysis Summary CMSPAS-BTV-15-001 ( 2016 ). http://cds.cern. ch/record/2138504 45. CMS Collaboration, Boosted top jet tagging at CMS . CMS Physics Analysis Summary CMS-PAS-JME-13-007 ( 2014 ). http:// cds.cern. ch/record/1647419 46. CMS Collaboration, Top tagging with new approaches . CMS Physics Analysis Summary CMS-PAS-JME-15-002 ( 2016 ). http:// cds.cern. ch/record/2126325 47. D.E. Kaplan , K. Rehermann , M.D. Schwartz , B. Tweedie , Top tagging: a method for identifying boosted hadronically decaying top quarks . Phys. Rev. Lett . 101 , 142001 ( 2008 ). doi:10.1103/ PhysRevLett.101.142001. arXiv:0806.0848 48. M. Dasgupta , A. Fregoso , S. Marzani , G.P. Salam , Towards an understanding of jet substructure . JHEP 09 , 029 ( 2013 ). doi:10. 1007/ JHEP09(2013)029 . arXiv:1307.0007 49. M. Cacciari , G.P. Salam , Dispelling the N 3 myth for the kt jetfinder . Phys. Lett. B 641 , 57 ( 2006 ). doi:10.1016/j.physletb. 2006 . 08.037. arXiv:hep-ph/0512210 50. Y.L. Dokshitzer , G.D. Leder , S. Moretti , B.R. Webber , Better jet clustering algorithms . JHEP 08 , 001 ( 1997 ). doi:10.1088/ 1126 - 6708 /1997/08/001. arXiv:hep-ph/ 9707323 51. M. Wobisch , T. Wengler , Hadronization corrections to jet cross sections in deep- inelastic scattering (1998). arXiv:hep-ph/9907280 52. CMS Collaboration, Performance of electron reconstruction and selection with the CMS detector in proton-proton collisions at √s = 8 TeV . JINST 10 , P06005 ( 2015 ). doi:10.1088/ 1748 - 0221 / 10/06/P06005. arXiv: 1502 .02701 53. CMS Collaboration, Performance of CMS muon reconstruction in pp collision events at √s = 7 TeV . JINST 7, P10002 ( 2012 ). doi:10.1088/ 1748 - 0221 /7/10/P10002. arXiv: 1206 .4071 54. CMS Collaboration, Reconstruction and identification of τ lepton decays to hadrons and ντ at CMS . JINST 11 , P01019 ( 2016 ). doi:10.1088/ 1748 - 0221 /11/01/P01019. arXiv: 1510 .07488 55. CMS Collaboration, Performance of photon reconstruction and identification with the CMS detector in proton-proton collisions at √s = 8TeV . JINST 10, P08010 ( 2015 ). doi:10.1088/ 1748 - 0221 / 10/08/P08010 56. J. Alwall et al., The automated computation of tree-level and nextto-leading order differential cross sections, and their matching to parton shower simulations . JHEP 07 , 079 ( 2014 ). doi:10.1007/ JHEP07(2014)079 . arXiv:1405.0301 57. NNPDF Collaboration, Parton distributions for the LHC Run II . JHEP 04 , 040 ( 2015 ). doi:10.1007/ JHEP04(2015)040 . arXiv:1410.8849 58. P. Nason , A new method for combining NLO QCD with shower Monte Carlo algorithms . JHEP 11 , 040 ( 2004 ). doi:10.1088/ 1126 - 6708 /2004/11/040. arXiv:hep-ph/ 0409146 59. S. Frixione , P. Nason , C. Oleari , Matching NLO QCD computations with Parton Shower simulations: the POWHEG method . JHEP 11 , 070 ( 2007 ). doi:10.1088/ 1126 - 6708 /2007/11/070. arXiv: 0709 . 2092 60. S. Alioli , P. Nason , C. Oleari , E. Re , A general framework for implementing NLO calculations in shower Monte Carlo programs: the POWHEG BOX . JHEP 06 , 043 ( 2010 ). doi:10.1007/ JHEP06(2010)043 . arXiv:1002.2581 61. E. Re , Single-top W t -channel production matched with parton showers using the POWHEG method . Eur. Phys. J. C 71 , 1547 ( 2011 ). doi:10.1140/epjc/s10052- 011 - 1547 -z. arXiv:1009.2450 62. T. Sjöstrand et al., An introduction to PYTHIA 8.2. Comput. Phys. Commun . 191 , 159 ( 2015 ). doi:10.1016/j.cpc. 2015 .01.024. arXiv:1410.3012 63. GEANT4 Collaboration, GEANT4-a simulation toolkit . Nucl. Instrum. Meth. A 506 , 250 ( 2003 ). doi:10.1016/ S0168-9002(03)01368- 8 64. S. Abdullin et al., The fast simulation of the CMS detector at LHC . J. Phys. Conf. Ser . 331 , 032049 ( 2011 ). doi:10.1088/ 1742 - 6596 / 331/3/032049 65. CMS Collaboration, Differential cross section measurements for the production of a W boson in association with jets in protonproton collisions at √s = 7 TeV . Phys. Lett . B 741 , 12 ( 2015 ). doi:10.1016/j.physletb. 2014 .12.003. arXiv:1406.7533 66. CMS Collaboration, Measurement of differential cross sections for top quark pair production using the lepton+jets final state in proton-proton collisions at 13 TeV . Phys. Rev . D ( 2016 ). arXiv:1610.04191. (Submitted) 67. CMS Collaboration, Observation of top quark pairs produced in association with a vector boson in pp collisions at √s = 8 TeV . JHEP 01 , 096 ( 2016 ). doi:10.1007/ JHEP01(2016)096 . arXiv:1510.01131 68. Y. Bai , H.-C. Cheng , J. Gallicchio , J. Gu , Stop the top background of the stop search . JHEP 07 , 110 ( 2012 ). doi:10.1007/ JHEP07(2012)110 . arXiv:1203.4813 69. M.L. Graesser , J. Shelton , Hunting mixed top squark decays . Phys. Rev. Lett . 111 , 121802 ( 2013 ). doi:10.1103/PhysRevLett. 111.121802. arXiv:1212.4495 70. Particle Data Group , K. A. Olive et al., Review of particle physics. Chin. Phys. C 38 , 090001 ( 2014 ). doi:10.1088/ 1674 - 1137 /38/9/ 090001 71. G. Polesello , D. Tovey , Supersymmetric particle mass measurement with boost-corrected cotransverse mass . JHEP 03 , 030 ( 2010 ). doi:10.1007/ JHEP03(2010)030 . arXiv:0910.0174 72. D. Tovey , On measuring the masses of pair-produced semi-invisibly decaying particles at hadron colliders . JHEP 04 , 034 ( 2008 ). doi:10. 1088/ 1126 - 6708 /2008/04/034. arXiv: 0802 .2879 73. CMS Collaboration, Measurement of the production cross sections for a Z boson and one or more b jets in pp collisions at √s = 7T eV . JHEP 06 , 120 ( 2014 ). doi:10.1007/ JHEP06(2014)120 . arXiv:1402.1521 74. CMS Collaboration, CMS luminosity measurement for the 2015 data taking period . CMS Physics Analysis Summary CMS-PASLUM-15-001 ( 2016 ). https://cds.cern. ch/record/2138682 75. S. Catani , D. de Florian, M. Grazzini , P. Nason , Soft gluon resummation for Higgs boson production at hadron colliders . JHEP 07 , 028 ( 2003 ). doi:10.1088/ 1126 - 6708 /2003/07/028. arXiv:hep-ph/ 0306211 76. M. Cacciari et al., The tt cross-section at 1.8 TeV and 1.96 TeV: a study of the systematics due to parton densities and scale dependence . JHEP 04 , 068 ( 2004 ). doi:10.1088/ 1126 - 6708 /2004/04/ 068. arXiv:hep-ph/ 0303085 77. CMS Collaboration, Interpretation of searches for supersymmetry with simplified models . Phys. Rev. D 88 , 052017 ( 2013 ). doi:10. 1103/PhysRevD.88.052017. arXiv:1301.2175 78. J. Alwall , P. Schuster , N. Toro , Simplified models for a first characterization of new physics at the LHC . Phys. Rev. D 79 , 075020 ( 2009 ). doi:10.1103/PhysRevD.79.075020 79. J. Alwall , M.-P. Le , M. Lisanti , J. Wacker , Model-independent jets plus missing energy searches . Phys. Rev. D 79 , 015005 ( 2009 ). doi:10.1103/PhysRevD.79.015005 80. LHC New Physics Working Group, D. Alves et al ., Simplified models for LHC new physics searches. J. Phys. G 39 , 105005 ( 2012 ). doi:10.1088/ 0954 - 3899 /39/10/105005. arXiv: 1105 .2838 81. T. Junk , Confidence level computation for combining searches with small statistics . Nucl. Instrum. Meth. A 434 , 435 ( 1999 ). doi:10. 1016/ S0168-9002(99)00498-2. arXiv:hep-ex/9902006 82. A.L. Read , Presentation of search results: the C L S technique. J. Phys. G 28 , 2693 ( 2002 ). doi:10.1088/ 0954 - 3899 /28/10/313 83. ATLAS and CMS Collaborations, LHC Higgs Combination Group, Procedure for the LHC Higgs boson search combination in Summer 2011 . Technical Report ATL-PHYS-PUB 2011-11, CMS NOTE 2011 / 005 ( 2011 ). http://cdsweb.cern. ch/record/1379837 84. G. Cowan , K. Cranmer , E. Gross , O. Vitells , Asymptotic formulae for likelihood-based tests of new physics . Eur. Phys. J. C 71 , 1554 ( 2011 ). doi:10.1140/epjc/s10052- 011 - 1554 -0. arXiv: 1007 .1727. (Erratum: doi:10.1140/epjc/s10052- 013 - 2501 -z) 85. CMS Collaboration, Measurement of the tt production cross section in the eμ channel in proton-proton collisions at √s = 7 and 8 TeV. JHEP 08, 029 (2016). doi:10.1007/JHEP08(2016)029. arXiv:1603.02303 Institute of High Energy Physics, Beijing, China M. Ahmad, J. G. Bian, G. M. Chen, H. S. Chen, M. Chen, Y. Chen7, T. Cheng, C. H. Jiang, D. Leggat, Z. Liu, F. Romeo, M. Ruan, S. M. Shaheen, A. Spiezia, J. Tao, C. Wang, Z. Wang, H. Zhang, J. Zhao Laboratoire Leprince-Ringuet, Ecole Polytechnique, IN2P3-CNRS, Palaiseau, France A. Abdulsalam, I. Antropov, S. Baffioni, F. Beaudette, P. Busson, L. Cadamuro, E. Chapon, C. Charlot, O. Davignon, R. Granier de Cassagnac, M. Jo, S. Lisniak, P. Miné, M. Nguyen, C. Ochando, G. Ortona, P. Paganini, P. Pigard, S. Regnard, R. Salerno, Y. Sirois, T. Strebler, Y. Yilmaz, A. Zabi, A. Zghiche Deutsches Elektronen-Synchrotron, Hamburg, Germany M. Aldaya Martin, T. Arndt, C. Asawatangtrakuldee, K. Beernaert, O. Behnke, U. Behrens, A. A. Bin Anuar, K. Borras16, A. Campbell, P. Connor, C. Contreras-Campana, F. Costanza, C. Diez Pardos, G. Dolinska, G. Eckerlin, D. Eckstein, T. Eichhorn, E. Eren, E. Gallo17, J. Garay Garcia, A. Geiser, A. Gizhko, J. M. Grados Luyando, A. Grohsjean, P. Gunnellini, A. Harb, J. Hauk, M. Hempel18, H. Jung, A. Kalogeropoulos, O. Karacheban18, M. Kasemann, J. Keaveney, C. Kleinwort, I. Korol, D. Krücker, W. Lange, A. Lelek, J. Leonard, K. Lipka, A. Lobanov, W. Lohmann18, R. Mankel, I.-A. Melzer-Pellmann, A. B. Meyer, G. Mittag, J. Mnich, A. Mussgiller, E. Ntomari, D. Pitzl, R. Placakyte, A. Raspereza, B. Roland, M. Ö. Sahin, P. Saxena, T. Schoerner-Sadenius, C. Seitz, S. Spannagel, N. Stefaniuk, G. P. Van Onsem, R. Walsh, C. Wissing INFN Sezione di Bolognaa , Università di Bolognab, Bologna, Italy G. Abbiendia , C. Battilana, D. Bonacorsia ,b, S. Braibant-Giacomellia ,b, L. Brigliadoria ,b, R. Campaninia ,b, P. Capiluppia ,b, A. Castroa ,b, F. R. Cavalloa , S. S. Chhibraa ,b, G. Codispotia ,b, M. Cuffiania ,b, G. M. Dallavallea , F. Fabbria , A. Fanfania ,b, D. Fasanellaa ,b, P. Giacomellia , C. Grandia , L. Guiduccia ,b, S. Marcellinia , G. Masettia , A. Montanaria , F. L. Navarriaa ,b, A. Perrottaa , A. M. Rossia ,b, T. Rovellia ,b, G. P. Sirolia ,b, N. Tosia ,b,15 INFN Sezione di Cataniaa , Università di Cataniab, Catania, Italy S. Albergoa ,b, S. Costaa ,b, A. Di Mattiaa , F. Giordanoa ,b, R. Potenzaa ,b, A. Tricomia ,b, C. Tuvea ,b INFN Sezione di Firenzea , Università di Firenzeb, Firenze, Italy G. Barbaglia , V. Ciullia ,b, C. Civininia , R. D'Alessandroa ,b, E. Focardia ,b, P. Lenzia ,b, M. Meschinia , S. Paolettia , L. Russoa ,29, G. Sguazzonia , D. Stroma , L. Viliania ,b,15 INFN Laboratori Nazionali di Frascati, Frascati, Italy L. Benussi, S. Bianco, F. Fabbri, D. Piccolo, F. Primavera15 INFN Sezione di Genovaa , Università di Genovab, Genoa, Italy V. Calvellia ,b, F. Ferroa , M. R. Mongea ,b, E. Robuttia , S. Tosia ,b INFN Sezione di Milano-Bicoccaa , Università di Milano-Bicoccab, Milan, Italy L. Brianzaa ,b,15, F. Brivioa ,b, V. Ciriolo, M. E. Dinardoa ,b, S. Fiorendia ,b,15, S. Gennaia , A. Ghezzia ,b, P. Govonia ,b, M. Malbertia ,b, S. Malvezzia , R. A. Manzonia ,b, D. Menascea , L. Moronia , M. Paganonia ,b, D. Pedrinia , S. Pigazzinia ,b, S. Ragazzia ,b, T. Tabarelli de Fatisa ,b INFN Sezione di Napolia , Università di Napoli 'Federico II' b, Napoli, Italy, Università della Basilicatac, Potenza, Italy, Università G. Marconid , Rome, Italy S. Buontempoa , N. Cavalloa ,c, G. De Nardo, S. Di Guidaa ,d ,15, M. Espositoa ,b, F. Fabozzia ,c, F. Fiengaa ,b, A. O. M. Iorioa ,b, G. Lanzaa , L. Listaa , S. Meolaa ,d ,15, P. Paoluccia ,15, C. Sciaccaa ,b, F. Thyssena INFN Sezione di Padovaa , Università di Padovab, Padova, Italy, Università di Trentoc, Trento, Italy P. Azzia ,15, N. Bacchettaa , M. Bellatoa , L. Benatoa ,b, M. Biasottoa ,30, D. Biselloa ,b, A. Bolettia ,b, P. Checchiaa , M. Dall'Ossoa ,b, P. De Castro Manzanoa , U. Dossellia , F. Gasparinia ,b, U. Gasparinia ,b, A. Gozzelinoa , S. Lacapraraa , M. Margonia ,b, A. T. Meneguzzoa ,b, J. Pazzinia ,b, N. Pozzobona ,b, P. Ronchesea ,b, F. Simonettoa ,b, E. Torassaa , S. Venturaa , M. Zanettia ,b, P. Zottoa ,b, G. Zumerlea ,b INFN Sezione di Paviaa , Università di Paviab, Pavia, Italy A. Braghieria , F. Fallavollitaa ,b, A. Magnania ,b, P. Montagnaa ,b, S. P. Rattia ,b, V. Rea , C. Riccardia ,b, P. Salvinia , I. Vaia ,b, P. Vituloa ,b INFN Sezione di Perugiaa , Università di Perugiab, Perugia, Italy L. Alunni Solestizia ,b, G. M. Bileia , D. Ciangottinia ,b, L. Fanòa ,b, P. Laricciaa ,b, R. Leonardia ,b, G. Mantovania ,b, M. Menichellia , A. Sahaa , A. Santocchiaa ,b INFN Sezione di Pisaa , Università di Pisab, Scuola Normale Superiore di Pisac, Pisa, Italy K. Androsova ,29, P. Azzurria ,15, G. Bagliesia , J. Bernardinia , T. Boccalia , R. Castaldia , M. A. Cioccia ,29, R. Dell'Orsoa , S. Donatoa ,c, G. Fedi, A. Giassia , M. T. Grippoa ,29, F. Ligabuea ,c, T. Lomtadzea , L. Martinia ,b, A. Messineoa ,b, F. Pallaa , A. Rizzia ,b, A. Savoy-Navarroa ,31, P. Spagnoloa , R. Tenchinia , G. Tonellia ,b, A. Venturia , P. G. Verdinia INFN Sezione di Romaa , Università di Romab, Rome, Italy L. Baronea ,b, F. Cavallaria , M. Cipriania ,b, D. Del Rea ,b,15, M. Diemoza , S. Gellia ,b, E. Longoa ,b, F. Margarolia ,b, B. Marzocchia ,b, P. Meridiania , G. Organtinia ,b, R. Paramattia , F. Preiatoa ,b, S. Rahatloua ,b, C. Rovellia , F. Santanastasioa ,b INFN Sezione di Torinoa , Università di Torinob, Turin, Italy, Università del Piemonte Orientalec, Novara, Italy N. Amapanea ,b, R. Arcidiaconoa ,c,15, S. Argiroa ,b, M. Arneodoa ,c, N. Bartosika , R. Bellana ,b, C. Biinoa , N. Cartigliaa , F. Cennaa ,b, M. Costaa ,b, R. Covarellia ,b, A. Deganoa ,b, N. Demariaa , L. Fincoa ,b, B. Kiania ,b, C. Mariottia , S. Masellia , E. Migliorea ,b, V. Monacoa ,b, E. Monteila ,b, M. Montenoa , M. M. Obertinoa ,b, L. Pachera ,b, N. Pastronea , M. Pelliccionia , G. L. Pinna Angionia ,b, F. Raveraa ,b, A. Romeroa ,b, M. Ruspaa ,c, R. Sacchia ,b, K. Shchelinaa ,b, V. Solaa , A. Solanoa ,b, A. Staianoa , P. Traczyka ,b INFN Sezione di Triestea , Università di Triesteb, Trieste, Italy S. Belfortea , M. Casarsaa , F. Cossuttia , G. Della Riccaa ,b, A. Zanettia Kyungpook National University, Taegu, Korea D. H. Kim, G. N. Kim, M. S. Kim, S. Lee, S. W. Lee, Y. D. Oh, S. Sekmen, D. C. Son, Y. C. Yang Instituto de Física de Cantabria (IFCA), CSIC-Universidad de Cantabria, Santander, Spain I. J. Cabrillo, A. Calderon, E. Curras, M. Fernandez, J. Garcia-Ferrero, G. Gomez, A. Lopez Virto, J. Marco, C. Martinez Rivero, F. Matorras, J. Piedra Gomez, T. Rodrigo, A. Ruiz-Jimeno, L. Scodellaro, N. Trevisani, I. Vila, R. Vilar Cortabitarte CERN, European Organization for Nuclear Research, Geneva, Switzerland D. Abbaneo, E. Auffray, G. Auzinger, M. Bachtis, P. Baillon, A. H. Ball, D. Barney, P. Bloch, A. Bocci, C. Botta, T. Camporesi, R. Castello, M. Cepeda, G. Cerminara, Y. Chen, D. d'Enterria, A. Dabrowski, V. Daponte, A. David, M. De Gruttola, A. De Roeck, E. Di Marco43, M. Dobson, B. Dorney, T. du Pree, D. Duggan, M. Dünser, N. Dupont, A. Elliott-Peisert, P. Everaerts, S. Fartoukh, G. Franzoni, J. Fulcher, W. Funk, D. Gigi, K. Gill, M. Girone, F. Glege, D. Gulhan, S. Gundacker, M. Guthoff, P. Harris, J. Hegeman, V. Innocente, P. Janot, J. Kieseler, H. Kirschenmann, V. Knünz, A. Kornmayer15, M. J. Kortelainen, K. Kousouris, M. Krammer1, C. Lange, P. Lecoq, C. Lourenço, M. T. Lucchini, L. Malgeri, M. Mannelli, A. Martelli, F. Meijers, J. A. Merlin, S. Mersi, E. Meschi, P. Milenovic44, F. Moortgat, S. Morovic, M. Mulders, H. Neugebauer, S. Orfanelli, L. Orsini, L. Pape, E. Perez, M. Peruzzi, A. Petrilli, G. Petrucciani, A. Pfeiffer, M. Pierini, A. Racz, T. Reis, G. Rolandi45, M. Rovere, H. Sakulin, J. B. Sauvan, C. Schäfer, C. Schwick, M. Seidel, A. Sharma, P. Silva, P. Sphicas46, J. Steggemann, M. Stoye, Y. Takahashi, M. Tosi, D. Treille, A. Triossi, A. Tsirou, V. Veckalns47, G. I. Veres20, M. Verweij, N. Wardle, H. K. Wöhri, A. Zagozdzinska35, W. D. Zeuner Paul Scherrer Institut, Villigen, Switzerland W. Bertl, K. Deiters, W. Erdmann, R. Horisberger, Q. Ingram, H. C. Kaestli, D. Kotlinski, U. Langenegger, T. Rohe Institute for Particle Physics, ETH Zurich, Zurich, Switzerland F. Bachmair, L. Bäni, L. Bianchini, B. Casal, G. Dissertori, M. Dittmar, M. Donegà, C. Grab, C. Heidegger, D. Hits, J. Hoss, G. Kasieczka, W. Lustermann, B. Mangano, M. Marionneau, P. Martinez Ruiz del Arbol, M. Masciovecchio, M. T. Meinhard, D. Meister, F. Micheli, P. Musella, F. Nessi-Tedaldi, F. Pandolfi, J. Pata, F. Pauss, G. Perrin, L. Perrozzi, M. Quittnat, M. Rossini, M. Schönenberger, A. Starodumov48, V. R. Tavolaro, K. Theofilatos, R. Wallny Universität Zürich, Zurich, Switzerland T. K. Aarrestad, C. Amsler49, L. Caminada, M. F. Canelli, A. De Cosa, C. Galloni, A. Hinzmann, T. Hreus, B. Kilminster, J. Ngadiuba, D. Pinna, G. Rauco, P. Robmann, D. Salerno, Y. Yang, A. Zucchetta National Central University, Chung-Li, Taiwan V. Candelise, T. H. Doan, Sh. Jain, R. Khurana, M. Konyushikhin, C. M. Kuo, W. Lin, Y. J. Lu, A. Pozdnyakov, S. S. Yu National Taiwan University (NTU), Taipei, Taiwan Arun Kumar, P. Chang, Y. H. Chang, Y. Chao, K. F. Chen, P. H. Chen, F. Fiori, W.-S. Hou, Y. Hsiung, Y. F. Liu, R.-S. Lu, M. Miñano Moya, E. Paganis, A. Psallidas, J. f. Tsai Imperial College, London, UK M. Baber, R. Bainbridge, O. Buchmuller, A. Bundock, D. Burton, S. Casasso, M. Citron, D. Colling, L. Corpe, P. Dauncey, G. Davies, A. De Wit, M. Della Negra, R. Di Maria, P. Dunne, A. Elwood, D. Futyan, Y. Haddad, G. Hall, G. Iles, T. James, R. Lane, C. Laner, R. Lucas63, L. Lyons, A.-M. Magnan, S. Malik, L. Mastrolorenzo, J. Nash, A. Nikitenko48, J. Pela, B. Penning, M. Pesaresi, D. M. Raymond, A. Richards, A. Rose, C. Seez, S. Summers, A. Tapper, K. Uchida, M. Vazquez Acosta65, T. Virdee15, J. Wright, S. C. Zenz Brunel University, Uxbridge, UK J. E. Cole, P. R. Hobson, A. Khan, P. Kyberd, I. D. Reid, P. Symonds, L. Teodorescu, M. Turner


This is a preview of a remote PDF: https://link.springer.com/content/pdf/10.1140%2Fepjc%2Fs10052-017-4853-2.pdf

A. M. Sirunyan, A. Tumasyan, W. Adam, E. Asilar, T. Bergauer, J. Brandstetter, E. Brondolin, M. Dragicevic, J. Erö, M. Flechl, M. Friedl, R. Frühwirth, V. M. Ghete, C. Hartl, N. Hörmann, J. Hrubec, M. Jeitler, A. König, I. Krätschmer, D. Liko, T. Matsushita, I. Mikulec, D. Rabady, N. Rad, B. Rahbaran, H. Rohringer, J. Schieck, J. Strauss, W. Waltenberger, C.-E. Wulz, V. Chekhovsky, O. Dvornikov, Y. Dydyshka, I. Emeliantchik, A. Litomin, V. Makarenko, V. Mossolov, R. Stefanovitch, J. Suarez Gonzalez, V. Zykunov, N. Shumeiko, S. Alderweireldt, E. A. De Wolf, X. Janssen, J. Lauwers, M. Van De Klundert, H. Van Haevermaet, P. Van Mechelen, N. Van Remortel, A. Van Spilbeeck, S. Abu Zeid, F. Blekman, J. D’Hondt, N. Daci, I. De Bruyn, K. Deroover, S. Lowette, S. Moortgat, L. Moreels, A. Olbrechts, Q. Python, K. Skovpen, S. Tavernier, W. Van Doninck, P. Van Mulders, I. Van Parijs, H. Brun, B. Clerbaux, G. De Lentdecker, H. Delannoy, G. Fasanella, L. Favart, R. Goldouzian, A. Grebenyuk, G. Karapostoli, T. Lenzi, A. Léonard, J. Luetic, T. Maerschalk, A. Marinov, A. Randle-conde, T. Seva, C. Vander Velde, P. Vanlaer, D. Vannerom, R. Yonamine, F. Zenoni, F. Zhang, A. Cimmino, T. Cornelis, D. Dobur, A. Fagot, M. Gul, I. Khvastunov, D. Poyraz, S. Salva, R. Schöfbeck, M. Tytgat, W. Van Driessche, E. Yazgan, N. Zaganidis, H. Bakhshiansohi, C. Beluffi, O. Bondu, S. Brochet, G. Bruno, A. Caudron, S. De Visscher, C. Delaere, M. Delcourt, B. Francois, A. Giammanco, A. Jafari, M. Komm, G. Krintiras, V. Lemaitre, A. Magitteri, A. Mertens, M. Musich, K. Piotrzkowski, L. Quertenmont, M. Selvaggi, M. Vidal Marono, S. Wertz, N. Beliy, W. L. Aldá Júnior, F. L. Alves, G. A. Alves, L. Brito, C. Hensel, A. Moraes, M. E. Pol, P. Rebello Teles, E. Belchior Batista Das Chagas, W. Carvalho, J. Chinellato, A. Custódio, E. M. Da Costa, G. G. Da Silveira, D. De Jesus Damiao, C. De Oliveira Martins, S. Fonseca De Souza, L. M. Huertas Guativa, H. Malbouisson, D. Matos Figueiredo, C. Mora Herrera, L. Mundim, H. Nogima, W. L. Prado Da Silva, A. Santoro, A. Sznajder, E. J. Tonelli Manganote, A. Vilela Pereira, S. Ahuja, C. A. Bernardes, S. Dogra, T. R. Fernandez Perez Tomei, E. M. Gregores, P. G. Mercadante, C. S. Moon, S. F. Novaes, Sandra S. Padula, D. Romero Abad, J. C. Ruiz Vargas, A. Aleksandrov, R. Hadjiiska, P. Iaydjiev, M. Rodozov, S. Stoykova, G. Sultanov, M. Vutova, A. Dimitrov, I. Glushkov, L. Litov, B. Pavlov, P. Petkov, W. Fang, M. Ahmad, J. G. Bian, G. M. Chen, H. S. Chen, M. Chen, Y. Chen, T. Cheng, C. H. Jiang, D. Leggat, Z. Liu, F. Romeo, M. Ruan, S. M. Shaheen, A. Spiezia, J. Tao, C. Wang, Z. Wang, H. Zhang, J. Zhao, Y. Ban, G. Chen, Q. Li, S. Liu, Y. Mao, S. J. Qian, D. Wang, Z. Xu, C. Avila, A. Cabrera, L. F. Chaparro Sierra, C. Florez, J. P. Gomez, C. F. González Hernández, J. D. Ruiz Alvarez, J. C. Sanabria, N. Godinovic, D. Lelas, I. Puljak, P. M. Ribeiro Cipriano, T. Sculac, Z. Antunovic, M. Kovac, V. Brigljevic, D. Ferencek, K. Kadija, B. Mesic, T. Susa, A. Attikis, G. Mavromanolakis, J. Mousa, C. Nicolaou, F. Ptochos, P. A. Razis, H. Rykaczewski, D. Tsiakkouri, M. Finger, M. Finger Jr., E. Carrera Jarrin, E. El-khateeb, S. Elgammal, A. Mohamed, M. Kadastik, L. Perrini, M. Raidal, A. Tiko, C. Veelken, P. Eerola, J. Pekkanen, M. Voutilainen, J. Härkönen, T. Järvinen, V. Karimäki, R. Kinnunen, T. Lampén, K. Lassila-Perini, S. Lehti, T. Lindén, P. Luukka, J. Tuominiemi, E. Tuovinen, L. Wendland, J. Talvitie, T. Tuuva, M. Besancon, F. Couderc, M. Dejardin, D. Denegri, B. Fabbro, J. L. Faure, C. Favaro, F. Ferri, S. Ganjour, S. Ghosh, A. Givernaud, P. Gras, G. Hamel de Monchenault, P. Jarry, I. Kucher, E. Locci, M. Machet, J. Malcles, J. Rander, A. Rosowsky, M. Titov, A. Abdulsalam, I. Antropov, S. Baffioni, F. Beaudette, P. Busson, L. Cadamuro, E. Chapon, C. Charlot, O. Davignon, R. Granier de Cassagnac, M. Jo, S. Lisniak, P. Miné, M. Nguyen, C. Ochando, G. Ortona, P. Paganini, P. Pigard, S. Regnard, R. Salerno, Y. Sirois, T. Strebler, Y. Yilmaz, A. Zabi, A. Zghiche, J.-L. Agram, J. Andrea, A. Aubin, D. Bloch, J.-M. Brom, M. Buttignol, E. C. Chabert, N. Chanon, C. Collard, E. Conte, X. Coubez, J.-C. Fontaine, D. Gelé, U. Goerlach, M. Jansová, A.-C. Le Bihan, P. Van Hove, S. Gadrat, S. Beauceron, C. Bernet, G. Boudoul, C. A. Carrillo Montoya, R. Chierici, D. Contardo, B. Courbon, P. Depasse, H. El Mamouni, J. Fan, J. Fay, S. Gascon, M. Gouzevitch, G. Grenier, B. Ille, F. Lagarde, I. B. Laktineh, M. Lethuillier, L. Mirabito, A. L. Pequegnot, S. Perries, A. Popov, D. Sabes, V. Sordini, M. Vander Donckt, P. Verdier, S. Viret, T. Toriashvili, Z. Tsamalaidze, C. Autermann, S. Beranek, L. Feld, M. K. Kiesel, K. Klein, M. Lipinski, M. Preuten, C. Schomakers, J. Schulz, T. Verlage, V. Zhukov, A. Albert, M. Brodski, E. Dietz-Laursonn, D. Duchardt, M. Endres, M. Erdmann, S. Erdweg, T. Esch, R. Fischer, A. Güth, M. Hamer, T. Hebbeker, C. Heidemann, K. Hoepfner, S. Knutzen, M. Merschmeyer, A. Meyer, P. Millet, S. Mukherjee, M. Olschewski. Searches for pair production of third-generation squarks in \(\sqrt{s}=13\) \(\,\text {TeV}\) pp collisions, The European Physical Journal C, 2017, 327, DOI: 10.1140/epjc/s10052-017-4853-2