UvA-DARE (Digital Academic Repository)
Search for chargino and neutralino production in final states with a Higgs boson and missing transverse momentum at √s = 13 TeV with the ATLAS detector
Aaboud, M.; The ATLAS Collaboration DOI 10.1103/PhysRevD.100.012006 Publication date 2019 Document Version Final published version Published in Physical Review D. Particles and Fields License CC BY Link to publication
Citation for published version (APA): Aaboud, M., & The ATLAS Collaboration (2019). Search for chargino and neutralino production in final states with a Higgs boson and missing transverse momentum at √s = 13 TeV with the ATLAS detector. Physical Review D. Particles and Fields, 100(1), [012006]. https://doi.org/10.1103/PhysRevD.100.012006
General rights It is not permitted to download or to forward/distribute the text or part of it without the consent of the author(s) and/or copyright holder(s), other than for strictly personal, individual use, unless the work is under an open content license (like Creative Commons).
Disclaimer/Complaints regulations If you believe that digital publication of certain material infringes any of your rights or (privacy) interests, please let the Library know, stating your reasons. In case of a legitimate complaint, the Library will make the material inaccessible and/or remove it from the website. Please Ask the Library: https://uba.uva.nl/en/contact, or a letter to: Library of the University of Amsterdam, Secretariat, Singel 425, 1012 WP Amsterdam, The Netherlands. You will be contacted as soon as possible.
UvA-DARE is a service provided by the library of the University of Amsterdam (https://dare.uva.nl) Download date:23 Sep 2021 PHYSICAL REVIEW D 100, 012006 (2019)
Search for chargino and neutralino production in finalpffiffi states with a Higgs boson and missing transverse momentum at s =13 TeV with the ATLAS detector
M. Aaboud et al.* (ATLAS Collaboration)
(Received 27 December 2018; published 30 July 2019)
A search is conducted for the electroweak pair production of a chargino and a neutralino 0 0 pp → χ˜1 χ˜2, where the chargino decays into the lightest neutralino and a W boson, χ˜1 → χ˜1W , while the neutralino decays into the lightest neutralino and a Standard Model-like 125 GeV Higgs 0 0 boson, χ˜2 → χ˜1h. Fully hadronic, semileptonic, diphoton, and multilepton (electrons, muons) final states with missing transverse momentum are considered in this search. Higgs bosons in the final state are identified by either two jets originating from bottom quarks (h → bb¯), two photons (h → γγ), → → → ττ 36 1 −1 orpffiffiffi leptons from the decay modes h WW, h ZZ or h . The analysis is based on . fb of s ¼ 13 TeV proton-proton collision data recorded by the ATLAS detector at the Large Hadron Collider. Observations are consistent with the Standard Model expectations, and 95% con- 0 fidence-level limits of up to 680 GeV in χ˜1 =χ˜2 mass are set in the context of a simplified supersymmetric model.
DOI: 10.1103/PhysRevD.100.012006
I. INTRODUCTION conservation hypothesis, the lightest supersymmetric Theoretical and experimental arguments suggest that particle (LSP) is stable. If the LSP is weakly interacting, the Standard Model (SM) is an effective theory valid up it may provide a viable DM candidate. The nature of the to a certain energy scale. The observation by the ATLAS LSP is defined by the mechanism that spontaneously and CMS collaborations of a particle consistent with the breaks supersymmetry and the parameters of the chosen SM Higgs boson [1–4] has brought renewed attention to theoretical framework. the mechanism of electroweak symmetry breaking and In the SUSY scenarios considered as benchmarks in this χ˜0 the hierarchy problem [5–8]: the Higgs boson mass is paper, the LSP is the lightest of the neutralinos ( ) which, χ˜ strongly sensitive to quantum corrections from physics at together with the charginos ( ), represent the mass γ very high energy scales and demands a high level of eigenstates formed from the mixture of the , W, Z and Higgs bosons’ superpartners (the higgsinos, winos and fine-tuning. Supersymmetry (SUSY) [9–14] resolves the binos). The neutralinos and charginos are collectively hierarchy problem by introducing for each known boson referred to as electroweakinos. Specifically, the electro- or fermion a new partner (superpartner) that shares the weakino mass eigenstates are designated in order of same mass and internal quantum numbers if supersym- increasing mass as χ˜ (i ¼ 1, 2) (charginos) and χ˜0 metry is unbroken. However, these superpartners have i j (j ¼ 1, 2, 3, 4) (neutralinos). In the models considered not been observed, so SUSY must be a broken symmetry in this paper, the compositions of the lightest chargino (χ˜ ) and the mass scale of the supersymmetric particles is as 1 χ˜0 yet undetermined. The possibility of a supersymmetric and next-to-lightest neutralino ( 2) are wino-like and the dark matter (DM) candidate [15,16] is related closely to two particles are nearly mass degenerate, while the lightest χ˜0 the conservation of R-parity [17]. Under the R-parity neutralino ( 1) is assumed to be bino-like. Naturalness considerations [18,19] suggest that the lightest of the charginos and neutralinos have masses near *Full author list given at the end of the article. the electroweak scale. Their direct production may be the dominant mechanism at the Large Hadron Collider (LHC) Published by the American Physical Society under the terms of if the superpartners of the gluon and quarks are heavier than the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to a few TeV. In SUSY models where the masses of the the author(s) and the published article’s title, journal citation, heaviest (pseudoscalar, charged) MSSM Higgs boson and and DOI. Funded by SCOAP3. the superpartners of the leptons have masses larger than
2470-0010=2019=100(1)=012006(37) 012006-1 © 2019 CERN, for the ATLAS Collaboration M. AABOUD et al. PHYS. REV. D 100, 012006 (2019)
FIG. 1. Diagrams illustrating the signal scenarios considered for the pair production of chargino and next-to-lightest neutralino targeted by the (a) hadronic (0lbb¯) and (b) 1lbb¯, (c) 1lγγ, (d) l l , 3l leptonic channel selections. In (a) and (b) the Higgs boson decays into two b-quarks. In (c), the diphoton channel is shown with h → γγ. In (d), the visible multilepton final state of the Higgs boson is shown. Leptons are either electrons or muons. those of the lightest chargino and next-to-lightest neutra- set to 125 GeV and its branching ratios are assumed to be 0 lino, the former might decay into the χ˜1 and a W boson the same as in the SM. The Higgs boson candidate can be 0 0 ¯ ¯ (χ˜1 → Wχ˜1), while the latter could decay into the χ˜1 and fully reconstructed with 0lbb, 1lbb and 1lγγ signatures, the lightest MSSM Higgs boson (h, SM-like), or Z boson while l l and 3l final states are sensitive to decays 0 0 (χ˜2 → h=Zχ˜1) [17,20,21]. The decay via the Higgs boson is h → WW, h → ZZ and h → ττ. Previous searches for dominant for many choices of the parameters as long as the charginos and neutralinos at the LHC targeting decays via mass-splitting between the two lightest neutralinos is larger the Higgs boson into leptonic final states have been than the Higgs boson mass and the higgsinos are heavier reported by the ATLAS [28] and CMS [29] collabora- than the winos. SUSY models of this kind, where sleptons tions; a search in the hadronic channel is also reported in are not too heavy although with masses above that of χ˜1 this paper. 0 and χ˜2, could provide a possible explanation for the discrepancy between measurements of the muon’s anoma- II. ATLAS DETECTOR lous magnetic moment g − 2 and SM predictions [22–25]. The ATLAS detector [30] is a multipurpose particle This paper presents a search in proton-proton collision pffiffiffi detector with a forward-backward symmetric cylindrical produced at the LHC at a center-of-mass energy s ¼ geometry and nearly 4π coverage in solid angle.1 The inner 13 TeV for the direct pair production of mass-degenerate tracking detector consists of pixel and microstrip silicon charginos and next-to-lightest neutralinos that promptly detectors covering the pseudorapidity region jηj < 2.5, χ˜ → χ˜0 χ˜0 → χ˜0 decay as 1 W 1 and 2 h 1. The search targets surrounded by a transition radiation tracker which enhances hadronic and leptonic decays of both the W and Higgs electron identification in the region jηj < 2.0. A new inner bosons. Three Higgs decay modes are considered: decays pixel layer, the insertable B-layer [31,32], was added at a into a pair of b-quarks, a pair of photons, or a pair of W or mean radius of 3.3 cm during the period between Run 1 and Z bosons or τ-leptons, where at least one of the W=Z=τ Run 2 of the LHC. The inner detector is surrounded by a decays leptonically. Four signatures are considered, illus- thin superconducting solenoid providing an axial 2 T trated in Fig. 1. All final states contain missing transverse magnetic field and by a fine-granularity lead/liquid-argon ⃗miss miss jηj 3 2 momentum (pT , with magnitude ET ) from neutrali- (LAr) electromagnetic calorimeter covering < . . nos, and in some cases neutrinos. Events are characterized A steel/scintillator-tile calorimeter provides hadronic cov- by the various decay modes of the W and Higgs bosons. erage in the central pseudorapidity range (jηj < 1.7). The The signatures considered have the following: jets, with end cap and forward regions (1.5 < jηj < 4.9) of the two of them originating from the fragmentation of ¯ b-quarks, called b-jets, and either no leptons [0lbb, 1ATLAS uses a right-handed coordinate system with its origin Fig. 1(a)], or exactly one lepton (l ¼ e, μ)[1lbb¯, at the nominal interaction point in the center of the detector. The Fig. 1(b)]; two photons and one lepton [1lγγ,Fig.1(c)]; positive x axis is defined by the direction from the interaction point to the center of the LHC ring, with the positive y axis only leptons [Fig. 1(d)] such that the final state contains pointing upwards, while the beam direction defines the z axis. either two leptons with the same electric charge, l l ,or Cylindrical coordinates ðr; ϕÞ are used in the transverse plane, ϕ three leptons, 3l. being the azimuthal angle around the z axis. The pseudorapidity η θ η ¼ − ðθ 2Þ A simplified SUSY model [26,27] is considered for the is defined in terms of the polar angle by ln tan = . Rapidity is defined as y ¼ 0.5 ln½ðE þ pzÞ=ðE − pzÞ where E optimization of the search and the interpretation of denotes the energy and p is the component of the momentum 0 0 0 z results. The χ˜1 → Wχ˜1 and χ˜2 → hχ˜1 decays are assumed palongffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi the beam direction. The angular distance ΔR is defined as to have 100% branching ratio. The Higgs boson mass is ðΔyÞ2 þðΔϕÞ2.
012006-2 SEARCH FOR CHARGINO AND NEUTRALINO PRODUCTION IN … PHYS. REV. D 100, 012006 (2019)
TABLE I. List of generators used for the different processes. Information is given about the underlying-event tunes, the PDF sets and the perturbative QCD highest-order accuracy (LO, NLO, next-to-next-to-leading order, NNLO, and next-to-next-to-leading-log, NNLL) used for the normalization of the different samples.
Process Generator þ fragmentation/hadronization Tune PDF set Cross-section
W=Z þ jets SHERPA-2.2.1 [47] Default NNPDF3.0NNLO [42] NNLO tt¯ POWHEG-BOX v2 [49,50] PERUGIA2012 [51] CT10 [52] NNLO þ PYTHIA-6.428 [53] þNNLL Single top POWHEG-BOX v1 or v2 PERUGIA2012 CT10 NNLO þ PYTHIA-6.428 þNNLL Diboson WW, WZ, ZZ SHERPA-2.2.1 Default NNPDF3.0NNLO NLO tt¯ þ X ttW=Z¯ MADGRAPH-2.2.2 [38] A14 [40] NNPDF2.3 NLO 4 top quarks þ PYTHIA-8.186 [39] tth¯ [email protected] UEEE5 [54] CT10 NLO þ HERWIG++-2.7.1 Wh, Zh PYTHIA-8.186 [48] A14 NNPDF2.3 NLO
hadronic calorimeter are made of LAr active layers with Monte Carlo (MC) samples of simulated events are used either copper or tungsten as the absorber material. A muon to model the signal and to aid in the estimation of SM spectrometer with an air-core toroid magnet system sur- background processes, with the exception of multijet rounds the calorimeters. Three layers of high-precision processes, which are estimated from data. All simulated tracking chambers provide coverage in the range jηj < 2.7, samples were produced using the ATLAS simulation while dedicated fast chambers allow triggering in the infrastructure [36] and GEANT4 [37], or a faster simulation region jηj < 2.4. The ATLAS trigger system consists of a based on a parameterization of the calorimeter response and hardware-based level-1 trigger followed by a software- GEANT4 for the other detector systems. The simulated based high-level trigger [33]. events were reconstructed with the same algorithm as that used for data. III. DATA AND MONTE CARLO SIMULATION SUSY signal samples were generated with ¯ MADGRAPH5_aMC@NLO v2.2.3 [38] (v2.3.3 for 0lbb) The data used in this analysis were collected in pp at leading order (LO) and interfaced to PYTHIA v8.186 [39] collisions at the LHC with a center-of-mass energy of 0l ¯ 13 TeV and a 25 ns proton bunch crossing interval during (v8.212 for bb) with the A14 [40] set of tuned param- 2015 and 2016. The full dataset corresponds to an eters (tune) for the modeling of the parton showering (PS), integrated luminosity of 36.1 fb−1 after requiring that all hadronization and underlying event. The matrix element detector subsystems were operational during data record- (ME) calculation was performed at tree level and includes ing. The uncertainty in the combined 2015 þ 2016 inte- the emission of up to two additional partons. The ME-PS grated luminosity is 2.1%. It is derived, following a matching was done using the CKKW-L [41] prescription, methodology similar to that detailed in Ref. [34], and with a matching scale set to one quarter of the chargino and using the LUCID-2 detector for the baseline luminosity next-to-lightest neutralino mass. The NNPDF23LO [42] measurements [35], from calibration of the luminosity scale parton distribution function (PDF) set was used. The using x − y beam-separation scans. Each event includes on cross sections used to evaluate the signal yields are average 13.7 and 24.9 inelastic pp collisions in the same calculated to next-to-leading-order (NLO) accuracy in bunch crossing (pileup) in the 2015 and 2016 datasets, the strong coupling constant, adding the resummation of respectively. In the 0lbb¯ and 1lbb¯ channels, events are soft gluon emission at next-to-leading-logarithm accuracy miss (NLO þ NLL) [43–45]. The nominal cross section and its required to pass ET triggers with period-dependent thresholds. These triggers are fully efficient for events uncertainty are taken as the midpoint and half-width of an miss 200 envelope of cross-section predictions using different PDF with ET > GeV reconstructed offline. Data for the 1lγγ signature were collected with a diphoton trigger sets and factorization and renormalization scales, as which selects events with at least two photons, with described in Ref. [46]. transverse momentum thresholds on the highest- and Background samples were simulated using different MC second-highest pT photons of 35 GeV and 25 GeV, event generators depending on the process. All background respectively. A combined set of dilepton and single-lepton processes are normalized to the best available theoretical triggers was used for event selection in the l l and 3l calculation of their respective cross sections. The event channels. generators, the accuracy of theoretical cross sections, the
012006-3 M. AABOUD et al. PHYS. REV. D 100, 012006 (2019) underlying-event parameter tunes, and the PDF sets used in neutral hadrons decaying into two photons [59]. Strict simulating the SM background processes are summarized identification requirements based on calorimeter shower in Table I. For all samples, except those generated using shapes are used to identify the so-called tight photons, which SHERPA [47], the EVTGEN v1.2.0 [48] program was used to are used in the 1lγγ analysis channel. In this case, photons simulate the properties of the bottom- and charm-hadron are required to satisfy an isolation criterion based on the sum decays. Several samples produced without detector simu- of the calorimeter energy in a cone of ΔR ¼ 0.4 centered on lation are employed to estimate systematic uncertainties the direction of the candidate photon, minus the energy of associated with the specific configuration of the MC the photon candidate itself and energy expected from pileup generators used for the nominal SM background samples. interactions. The resulting isolation transverse energy is They include variations of the renormalization and factori- 2 45 þ 0 022 γ γ required to be less than . GeV . × ET, where ET zation scales, the CKKW-L matching scale, as well as is the candidate photon’s transverse energy. Photons are different PDF sets and fragmentation/hadronization mod- calibrated using comparisons of data with MC simulation 25 els. Details of the MC modeling uncertainties are discussed [57] and required to have ET > GeV. For both the in Sec. VII. electrons and photons, additional criteria are applied to remove poor quality or fake electromagnetic clusters result- IV. EVENT RECONSTRUCTION AND ing from instrumental problems. OBJECT DEFINITIONS Muon candidates are reconstructed from matching tracks Common event-quality criteria and object reconstruction in the inner detector and muon spectrometer. They are definitions are applied for all analysis channels, including required to meet medium quality and identification criteria standard data-quality requirements to select events taken and to be isolated, as described in Ref. [61], and to have 20 10 1lγγ during optimal detector operation. In addition, each analy- pT > GeV (pT > GeV for the analysis) and jηj 2 5 sis channel applies selection criteria that are specific to the < . . These muons are used in the overlap removal objects and kinematics of interest in those final states, procedure and to apply lepton selections and vetoes in the which are described in Sec. VI. various analysis channels, in some cases with additional Events are required to have at least one primary vertex, selections applied. Events containing muons from calo- defined as the vertex associated with at least two tracks with rimeter punch-through or poorly measured tracks are p > 0.4 GeV and with the highest sum of squared trans- rejected if any muon has a large relative q=p error, or T σð Þ j j 0 2 verse momenta of associated tracks [55]. Quality criteria q=p = q=p > . , where q is the charge of the track and are imposed to reject events that contain at least one jet p is the momentum. Cosmic-ray muons are rejected after arising from noncollision sources or detector noise [56]. the muon-jet overlap removal by requiring the transverse Electron candidates are reconstructed from energy clus- and longitudinal impact parameters to be jd0j < 0.25 mm ters in the electromagnetic calorimeter and inner-detector and jz0 sin θj < 0.5 mm, respectively. tracks. They are required to satisfy the loose likelihood Jets are reconstructed from three-dimensional topologi- identification criteria, have B-layer hits (the loose require- cal energy clusters [62] in the calorimeter using the anti-kt ment), and be isolated [57,58]. These identification jet algorithm [63] with a radius parameter of 0.4. Each criteria are based on several properties of the electron topological cluster is calibrated to the electromagnetic scale candidates, including calorimeter-based shower shapes, prior to jet reconstruction. The reconstructed jets are then 2 inner-detector track hits and impact parameters, and com- calibrated to the energy scale of stable final state particles parisons of calorimeter cluster energy to track momentum. in the MC simulationpffiffiffi by a jet energy scale (JES) correction Corrections for energy contributions due to pileup are derived from s ¼ 13 TeV data and simulations [64]. included. For all but the 1lγγ channel, electrons are also Further selections are applied to reject jets within jηj < 2.4 20 jηj 2 47 1lγγ required to have pT > GeV and < . ; for the that originate from pileup interactions by means of a 15 channel they are required to have pT > GeV and multivariate algorithm using information about the tracks jηj < 2.37. These electrons are used in the overlap removal matched to each jet [64,65]. Candidate jets are required to 20 jηj 2 8 procedure that is described below, and to apply lepton have pT > GeV and < . . selections and vetoes in the various analysis channels, in A jet is tagged as a b-jet by means of a multivariate some cases with additional selections applied. algorithm called MV2c10 using information about the Photon candidates are reconstructed from energy impact parameters of inner-detector tracks matched to clusters in the electromagnetic calorimeter [59] in the the jet, the presence of displaced secondary vertices, and region jηj < 2.37, after removing the transition region the reconstructed flight paths of b-andc-hadrons inside the between barrel and end cap calorimeters, 1.37
012006-4 SEARCH FOR CHARGINO AND NEUTRALINO PRODUCTION IN … PHYS. REV. D 100, 012006 (2019) operating points are available, corresponding to various event is useful for rejecting events with mismeasured ¯ miss efficiencies obtained in tt simulated events. The 77% jet energies leading to ET in the event, and is efficiency point was found to be optimal for most SUSY defined as models considered in this paper and is used in the analysis. This configuration corresponds to a background rejection of 4j miss jet Δϕ ¼ min ≤4Δϕðp⃗ ; p⃗ Þ 6forc-jets, 22 for τ-leptons and 134 for light-quark and min i T T;i – ¯ gluon jets [66 68], estimated using tt simulated events. Δϕ miss where mini≤4 selects the jet the minimizes . The ET in the event is defined as the magnitude of the (e) The effective mass meff is defined as the scalar sum negative vector sum of the pT of all selected and calibrated miss physics objects in the event, with an extra term added to of the pT of jets, leptons and ET , which aids in account for soft energy in the event that is not associated establishing the mass scale of the processes being with any of the selected objects. This soft term is calculated probed, and is defined as from inner-detector tracks matched to the primary vertex to XNjet NXlepton ¼ jet þ l þ miss make it more resilient to pileup contamination [69]. meff pT;i pT;j ET : Overlaps between reconstructed objects are accounted i j for and removed in a prioritized sequence. If a recon- structed muon shares an inner-detector track with an (f) mbb¯ is the invariant mass of the two leading b-jets in electron, the electron is removed. Jets within ΔR ¼ 0.2 the event, and serves as a selection criterion for jet of an electron are removed. Electrons that are reconstructed pairs to be considered as Higgs boson candidates. within ΔR ¼ 0.4 of any surviving jet are then removed, (g) mqq¯ corresponds to the invariant mass of the two 0l ¯ Δ ¼ highest-pT jets in the event not identified as b-jets. except in the case of the bb channel, where R ¯ ð0 4 0 04 þ 10 e Þ This observable, used in the 0lbb channel, serves as min . ; . GeV=pT , thereby allowing a high-pT electron to be slightly closer to a jet than ΔR ¼ 0.4. If a jet a selection criterion for jet pairs to be considered as is reconstructed within ΔR ¼ 0.2 of a muon and the jet has W boson candidates. fewer than three associated tracks or the muon energy (h) mγγ is the invariant mass of the two leading photons constitutes a large fraction (>50%) of the jet energy, then in the event, and serves as a selection criterion for the jet is removed. Muons reconstructed within a cone of photon pairs to be considered as Higgs boson Δ ¼ ð0 4 0 04 þ 10 μ Þ candidates. size R min . ; . GeV=pT around the axis of any surviving jet are removed. If an electron (muon) and (i) mljðjÞ is the invariant mass of the jet (when requiring a photon are found within ΔR ¼ 0.4, the object is inter- exactly one jet), or the two leading jets system (when preted as electron (muon) and the overlapping photon is requiring two or more jets), and the closest lepton. Δ removed from the event. Finally, if a jet and a photon are The angular distance R is used as the distance found within ΔR<0.2, the object is interpreted as photon measure between the lepton and the jet. and the overlapping jet is removed from the event; other- (j) mlll is the invariant mass of the three selected Δ 0 4 leptons. wise, if R< . , the object is interpreted as a jet and the miss overlapping photon is discarded. (k) mT is the transverse mass formed by the ET and the leading lepton in the event. It is defined as V. KINEMATIC REQUIREMENTS AND EVENT qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi VARIABLES ¼ 2 l missð1 − Δϕðl ⃗missÞÞ mT pTET cos ; pT Different analysis channels’ signal regions are optimized to target different mass hierarchies of the particles involved. and is used to reduce the W þ jets and tt¯ back- The event selection criteria are defined on the basis of grounds. b;min kinematic requirements for the objects described in the (l) mT is the minimum transverse mass formed by miss previous section and event variables are presented below. In ET and up to two of the highest-pT b-jets in the the following, jets are ordered according to decreasing pT, event, defined as and pT thresholds depend on the analysis channel. jηj 2 8 (a) Njet is the number of jets with < . and pT b;min mT above an analysis-dependent pT threshold. qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi jηj 2 5 (b) Nb-jet is the number of b-jets with < . with pT ¼ 2 b-jeti missð1− Δϕð ⃗miss ⃗b-jeti ÞÞ mini≤2 pT ET cos pT ;pT : above an analysis-dependent pT threshold. (c) Δηll is the pseudorapidity difference between the two leading leptons in the case of multilepton where mini≤2 selects the b-jet the minimizes the channels. transverse mass. 4 γ Δϕ j miss γ W i (d) The minimum azimuthal angle min between the (m) The lepton-ET - transverse mass mT is calcu- ⃗miss ⃗ th γ pT and the pT of each of the four leading jets in the lated with respect to the i photon i, ordered in
012006-5 M. AABOUD et al. PHYS. REV. D 100, 012006 (2019)
miss terms of decreasing ET, the ET , and the identified with a profile likelihood fit [74], referred to as a back- lepton l. It is defined as ground-only fit. The background-only fit uses the observed
γ 2 γ event yield in the associated CRs as a constraint to adjust ð W i Þ ¼ 2 i missð1 − Δϕðγ ⃗missÞÞ mT ET ET cos i; pT the normalization of the dominant background processes þ 2 l missð1 − Δϕðl ⃗missÞÞ assuming that no signal is present. The CRs are designed to pTET cos ; pT γ be enriched in specific background contributions relevant þ 2 i lð1 − Δϕðγ lÞÞ ET pT cos i; : to the analysis, while minimizing the signal contamination, and they are orthogonal to the SRs. The inputs to the (n) m is the contransverse mass variable [70,71] and CT background-only fit for each SR include the number of is defined for the bb¯ system as events observed in the associated CR and the number of qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi events predicted by simulation in each region for all ¼ 2 b1 b2 ð1 þ Δϕ Þ mCT pT pT cos bb ; background processes. They are both described by Poisson statistics. The systematic uncertainties, described b1 b2 where pT and pT are transverse momenta of the two in Sec. VII, are included in the fit as nuisance parameters. Δϕ leading b-jets and bb is the azimuthal angle They are constrained by Gaussian distributions with widths between them. It is one of the main discriminating corresponding to the sizes of the uncertainties and are variables in selections targeting Higgs bosons treated as correlated, when appropriate, between the various decaying into b-quarks and is effective in suppressing regions. The product of the various probability density the background from top-quark pair production. functions forms the likelihood, which the fit maximizes by (o) mT2 is referred to as the stransverse mass and is adjusting the background normalization and the nuisance closely related to mT. It is used to bound the masses parameters. Finally, the reliability of the MC extrapolation of particles produced in pairs and each decaying into of the SM background estimates outside of the control one particle that is detected and another particle that regions is evaluated in validation regions orthogonal to CRs contributes to the missing transverse momentum and SRs. [72,73]. In the case of a dilepton final state, it is defined by A. Fully hadronic signature (0lbb¯) ¼ ½ ð ð ⃗l;1 ⃗Þ ð ⃗l;2 ⃗miss − ⃗ÞÞ mT2 min max mT pT ;qT ;mT pT ;pT qT ; qT 1. Event selection The fully hadronic analysis channel exploits the large where q⃗ is the transverse vector that minimizes the T branching ratios for both W → qq¯ and h → bb¯. Missing ⃗l;1 larger of the two transverse masses mT, and pT and transverse momentum triggers are used for the trigger ⃗l;2 pT are the leading and subleading transverse selection for the analysis, with an offline requirement of momenta of the two leptons in the pair. Emiss > 200 GeV. Stringent event selections based on the 1lγγ Δϕ T (p) The variable W;h is the azimuthal angle masses of both the W and Higgs boson candidates, the between the W boson and Higgs boson candidates. presence of exactly two b-jets, and the kinematic relation- ⃗l The W boson is defined by the sum of the lepton pT ships of the final-state jets and Emiss, are required in order to ⃗miss T and pT vectors, and the Higgs boson by the sum of reduce the significant backgrounds from tt¯, Z þ jets, W þ the transverse momentum vectors of the two photons. jets and single-top Wt production. Events are characterized 30 by having four or five jets with pT > GeV, exactly two VI. ANALYSIS STRATEGY of which are identified as b-jets, and large meff, mCT, and b;min The hadronic and leptonic decay modes of the W and mT . Two signal regions are defined, specifically target- 0 Higgs bosons are divided into four independent and ing either high (HM) or low (LM) χ˜2 and χ˜1 masses mutually exclusive analysis channels according to key (SRHad-High and SRHad-Low, respectively). The selec- b;min features of the visible final state: hadronic decays of both tions used are shown in Table II. The meff and mT the W and h (0lbb¯, Sec. VI A); hadronic h decays with selections are particularly effective in reducing the tt¯ leptonic W decays (1lbb¯, Sec. VI B); diphoton h decays contributions, which is the dominant background for both with leptonic W decays (1lγγ, Sec. VI C); multilepton h signal regions. The Z þ jets and single-top contributions decays via W, Z, τ and leptonic W decays (l l and 3l, are also significant, whereas the contribution from multijet Sec. VI D). Event selections and background estimation production is found to be negligible and is not included. methods specific to each analysis channel are described Control regions are used to constrain the normalizations of here, as well as the signal, control, and validation region the tt¯, Z þ jets, and Wt backgrounds with the data, while definitions (SR, CR, and VR, respectively). other processes are estimated using simulation. The bb¯ The expected SM backgrounds are determined sepa- invariant mass is required to be consistent with the Higgs 105 135 rately for each SR, and independently for each channel, boson mass, 012006-6 SEARCH FOR CHARGINO AND NEUTRALINO PRODUCTION IN … PHYS. REV. D 100, 012006 (2019) TABLE II. Signal region definitions for the fully hadronic 0lbb¯ analysis channel. Variable SRHad-High SRHad-Low ¼ 0 ¼ 0 Nlepton Njet (pT>30 GeV) ∈ ½4; 5 ∈ ½4; 5 Nb-jet ¼ 2 ¼ 2 Δϕ4j >0.4 >0.4 min miss 250 200 ET [GeV] > > 900 700 meff [GeV] > > ∈ ½105 135 ∈ ½105 135 mbb¯ [GeV] ; ; mqq¯ [GeV] ∈ ½75; 90 ∈ ½75; 90 140 190 mCT [GeV] > > b;min >160 >180 mT [GeV] All CRs and VRs select sidebands in the mbb¯ spectrum in order to remain orthogonal to the two SRs. These are further described in Sec. VI A 2. 2. Background estimation The background contributions to SRHad-High and SRHad-Low are estimated using fits to the data for tt¯, Z þ jets, and single-top production in specially designed control regions. The three control regions used for estimating the tt¯ (CRHad-TT), Z þ jets (CRHad-Zj), and Wt (CRHad-ST) contributions are further divided into high-mass (HM) and low-mass (LM) categories in order to follow the design of the SRs. These control regions are defined primarily by b;min ¯ inverting the selections on mbb, mCT, mT , and by requiring the presence of a lepton in some cases. The tt¯ b;min 140 background is estimated using mCT;mT < GeV and FIG. 2. Comparisons of data with SM predictions in tt¯ and 135 mbb¯ > GeV selections, while retaining the other SR Z þ jets control regions for representative kinematic distribu- requirements. This approach isolates the tt¯ contribution miss ¯ tions: (a) ET for the tt high-mass control region and (b) mbb¯ for while suppressing single-top and Z þ jets events, yielding a the Z þ jets low-mass control region. Predictions from simulation sample estimated to be 94% pure in tt¯ events with are shown after the background-only fit. The arrow indicates negligible signal contamination. Background events from the selection on that variable used to define the corresponding Wt are estimated by requiring exactly one lepton and CRs. The uncertainty bands include statistical and systematic b;min 200 180 ¯ 195 uncertainties. mCT > GeV, mT > GeV, and mbb > GeV, and relaxing the meff requirement for HM to 700 þ meff > GeV. The Z jets contribution is isolated for tt¯, Z þ jets, and Wt in the high-mass (low-mass) 2l using an opposite-sign, same-flavor, high-pT require- signal region, respectively. The errors include statistical l 140 75 105 ment with pT;1 > GeV and