JHEP10(2018)064 , miss T E Springer , October 9, 2018 August 14, 2018 : : September 12, 2018 : Received Published and , Accepted b,c Published for SISSA by S. Moretti https://doi.org/10.1007/JHEP10(2018)064 a [email protected] , [email protected] , final state we explore the sensitivity of a current LHC H.U. Flaecher, c miss T d E analysis to such scenarios. + b ¯ bb ¯ miss T b . 3 E U. Ellwanger, 1807.10672 The Authors. c Phenomenology a,b,c

We examine scenarios in the Next-to-Minimal Supersymmetric , [email protected] decay giving a b ¯ b → Universit´eParis-Saclay, 91405 Orsay, France E-mail: [email protected] [email protected] Bristol, BS8 1TL, United Kingdom School of Physics andHighfield, Astronomy, Southampton, University SO17 of 1BJ, Southampton, UnitedParticle Kingdom Physics Department, Rutherford AppletonChilton, Laboratory, Didcot, Oxon, OX11 0QX,Laboratoire de United Physique Kingdom Th´eorique,UMR 8627, CNRS, Universit´ede Paris-Sud, H.H. Wills Physics Laboratory, University of Bristol, b c d a Open Access Article funded by SCOAP ArXiv ePrint: is very low.H Performing two-dimensional parameter scansgeneral-purpose jets+ and focusing on the hadronic Keywords: Abstract: (NMSSM), where pair-produced squarksin and a stable decay via (SM)-like two as Higgs cascades, Lightest , each Supersymmetric ending with mass (LSP) spectra and such a that Standard the Model missing transverse energy, A. Titterton, C.H. Shepherd-Themistocleous Exploring sensitivity to NMSSMmissing signatures transverse with energy low at the LHC JHEP10(2018)064 3 18 12 13 20 20 ], have pushed 7 , 6 well in excess of 1 TeV and ˜ g M kinematic variable [ T 10 10 α 4 23 – 1 – 16 and mass ˜ q 7 M ]. These experimental limits are of course dependent 7 miss T also as high as 1 TeV for certain regions of parameter 13 ], utilising the E – 0 5 1 8 ˜ 14 X M 4 8 ], this has allowed lower bounds to be placed on the masses of T 11 4 ∗ – H -tagged hadronic jets φ 1 b 1 24 miss T H Recent LHC analyses such as [ = 7–8 TeV and a further three with 13 TeV, searches for physics Beyond the SM (BSM) 6.1 MSSM-like scenarios with light LSP 5.3 Number of5.4 hadronic jets Number of 5.5 Minimum ∆ 5.6 Angular separation between bottom jets from decays 4.3 Signal, background4.4 and data event Validation yields of cut and count analysis tools 5.1 Total scalar 5.2 4.1 Two-dimensional mass4.2 scans Event selections 2.1 Other factors affecting signal processes s the mass of aspace neutralino of LSP simplified SUSY modelsupon [ various properties ofdecay the branching fractions decay and cascade thestill such kinematic ruled as distributions out of the a the masses large decay area of products, of but other parameter have sparticles, space the in the MSSM. have so farsymmetry not (SUSY) observed [ anysupersymmetric significant excesses. such as squarks, Particularly, gluinos, in and searches . the for lower bounds Super- on the squark mass 1 Introduction After four√ years of -proton collisions at the LHC with a centre-of-mass energy of 6 Results 7 Conclusion 5 Signal properties 3 Simulation techniques 4 Mass scans and event selections Contents 1 Introduction 2 Testing scenarios with low JHEP10(2018)064 , in but ˜ g scalar miss T M T E p or ˜ is the gauge q ] is examined, . In a scenario M S 5 , is just slightly 0 2 H ˜ χ M generally have weaker throughout this paper, as well as low , where ˆ ], having only weak cou- S T of CMS data from Run II range. In order to gain un- miss T 16 miss H T in an event within various 1 2 E E ], we consider the case where c − 15 9 fb . scenarios. In this SUSY scenario, 5 TeV/ . , meaning events whose jet not much smaller than miss T is reduced considerably. An example of T p 0 2 E ˜ χ hadronic jets miss T M = 13 TeV. E is defined as the norm of the two-dimensional s – 2 – √ miss T H . ] scenarios proposed in [ and various isolation requirements. On the other hand miss ], this will be used in lieu of jets, implying events with low η T 14 -based general-purpose BSM analysis in [ . In addition, long SUSY decay cascades often imply many which is frequently used in searches involving hadronic final 5 , the now heavy LSP will gain fairly little extra momentum , T H 2 T T c p α p miss T miss T ], parameter scans are developed which envelop these points E E 15 ], these NMSSM scenarios with low 15 , denoted T throughout [ + 125 GeV/ H , since the much heavier Higgs boson would inherit most of the momentum 0 1 H ˜ χ miss T ]. M H M 17 + . , taken from [ 0 1 , is also very high. ˜ χ 1 T miss T M includes all objects within these acceptance regions, such as . Following the H E Starting from the eight Benchmark Points (BPs), denoted BP1, . . . , BP8 and presented A similar quantity to As discussed in [ Furthermore, in the case where we have If, in contrast, the LSP were the singlino of the NMSSM [ Starting from the NMSSM [ The majority of the SUSY search effort so far has relied upon the notion that an miss -parity-conserving supersymmetric model is expected to generate events featuring large T -problem [ approach of noting that these quantities should betight essentially veto equivalent on in events the containing presence of final a state sufficiently leptons. in table of the LHC at a centre-of-mass energy of states is missing- vector sum ofacceptance the regions, transverse such momenta as E of the constraints imposed by experimentalRun searches, I however, of the the LHC, focusderstanding with of therein all the sparticle was sensitivity masses mainly to inefforts the on these 1–1 models in which Run might be II,which attained focuses the with on current CMS an search all-hadronic final state and uses 35 compared to the initially producedthe sparticles, now even though (in it comparison)decay will cascade lighter inherit will more Higgs mean of boson. low- thiscontrast than to Additionally, the the typical signature small of mass a gaps more minimal in SUSY the model. the singlino is thesinglet fermionic field counterpart to ofµ which a the singlet superfield coupling of the twostill Higgs close to doublets of the MSSM solves the high plings to sparticles, thenplays the the role heavier sparticles ofand will the a decay MSSM-like Higgs into LSP. boson, a thus This bino-like allowing NLSP NLSP, for which then these decays low- into a soft, “true” LSP the Next-to-LSP (NLSP) decays intowhere a the LSP LSP plus is a verysuch SM-like soft a Higgs in scenario boson, this would decay beabove the for a very lightfrom LSP the where NLSP. the NLSPskip However, mass, this in NLSP the decay MSSM, step most and of decay the straight into heavier the sparticles LSP would instead, prefer thus to still generating R missing transverse energy, hadronic jets with largesum, transverse momentum, JHEP10(2018)064 miss T final E and the miss T 2 c E + b final state will leptons, produced miss T ] τ 2 E c for each of the eight BPs = 750 = 750 = 750 = 750 [GeV/ ˜ ˜ t t b b ˜ ˜ 1 b ˜ , weakening the lower bounds ˜ M M t, M M decoupled decoupled decoupled decoupled 2 c ]. M 15 ] 2 c ] concentrates specifically on all-hadronic 5 with the squark and gluino masses around 900 1010 1410 1300 1410 1110 1300 1200 [GeV/ 2 ˜ g c is quite low for these NMSSM scenarios, more – 3 – M leptons produced, thus seeking a 4 Original BPs in [ miss T τ ] miss T 2 E c E LSP. 2 c Table 1. 1000 1400 1100 1500 1400 1100 1500 1400 [GeV/ ˜ q distributions are shown in figure M miss T E BP1 BP2 BP3 BP4 BP5 BP6 BP7 BP8 Point ] the Higgs boson in question is SM-like with a mass of around 125 GeV/ , featuring a LSP mass of 3 GeV/ squarks and a 3 GeV/ , along with a simplified MSSM-like scenario which exhibits typically larger values . This MSSM-like scenario features the same sparticle masses as BP1, such as 2 2 1 15 c c miss T However, since the analysis considered in [ In [ Here, it is clear that the mean For example, the -antilepton pair. This was chosen in order to suppress SM background processes E τ final states, this paperbottom will quark-antiquark pair, focus with onstate. no the This scenario retains where thefraction each maximum for Higgs the proportion Higgs boson of boson signal decays decay events, into is since that a the to a largest bottom branching quark-antiquark pair. Likewise, the final state under consideration contains twowhen bottom one quark of jets the and two Higgs two a decays to a bottomcompared quark-antiquark pair with and an the all-hadronic other final to a state bottom such as quark-antiquark that pair. where both Higgs bosons decay into 1 TeV/ akin to that fromcurrent SM experimental processes searches such concentrating aslikely on top not quark a be pair hadronic optimally production. jets tuned plus In to this turn, type this of suggests SUSY signature. on the masses of the squarks and gluino. in table 1 TeV/ of Current searches forbounds SUSY on have sparticle masses yet for simplifiedspace, to models. it find is However, for possible substantial certain that regionssignal evidence, a of whilst supersymmetric parameter increasing still decay cascade having the may a in lower LSP fact give mass a of very only low a few GeV/ in order to determinesimulated events the are current then sensitivitydetector generated of response at calculated this parton-level, for CMS each with analysis. point the in decays, these Monte hadronisation scans. Carlo2 and (MC) Testing scenarios with low JHEP10(2018)064 (2.1) scenar- . miss T T E p ] at NLO. 19 ] along with an MSSM- 15 ] decays the particles according to 20 . and larger Higgs boson T , our final state decays will appear as H H (jets) + LSP b ¯ H b jets from the SUSY decay cascade as well as T – 4 – → → p H NLSP distribution for the eight benchmark points defined in [ ], which has built-in support for NMSSM processes at Leading Order (LO), 18 miss T E The initial squarks, antisquarks and gluinos, plus up to two additional hard jets, are As the squark and gluino masses increase, so does the mass gap between these and Thus denoting this SM-like Higgs boson as is used. However MadGraph cannotfore, generate the NMSSM events cross-section at value Next-to-LOproduction stated (NLO), there- by cross-sections MadGraph are calculated is (inclusively) not using used. Prospino Instead, [ squarkgenerated and by gluino MadGraph at LO, after which Pythia8 [ 3 Simulation techniques Firstly, the mass spectra areios. chosen In in order order to toMadGraph generate [ simulated exploit events these for light-LSP, processes low- calculated at matrix element level, cross-section, which is toand a gluino. large extent dependent only upon thethe NLSP. masses of This in theincreased turn squarks NLSP will momentum, lead resulting to in high- higher 2.1 Other factorsWhilst affecting the signal masses processes ofof the the kinematical and Higgs variables, boson these will have dictate no the real behaviour effect of on some the initial sparticle production contribution towards the background fromjets SM formed events through involving Quantum large Chromo-Dynamics numberswill (QCD) of also interactions, hadronic i.e., increase. multijet events, Figure 1. like simplified scenario. JHEP10(2018)064 , we see 1 lighter than the squarks. In 2 c heavier than squarks. The gluinos thus ], denoted BP1 and BP2 in table 2 c 15 – 5 – . show examples of the processes by which we may q 2 + 0 2 ˜ χ ]. → ]. q 23 , ; ˜ 21 q 22 + ˜ q → g shows example decay cascades for each of these BPs. BP2 and BP4 have been 2 The Feynman diagrams in figure Figure BP7 and BP8 also involve the stop (BP7) or sbottom (BP8) squark, but with the For BP5 and BP6 the squarks are lighter than the gluino, as in points BP1 and BP3 and BP4 differ in that the gluino is 200 GeV/ Taking, for example the, first two BPs in [ produced along with the initial sparticles or typically softer radiated jets. squark/quark. produce a final stateproduce with one two LSPs squark and andhowever, two one both Higgs gluino squark-squark bosons. and directly, gluino-gluino In withditionally pair each each in production of decaying the processes these via are event diagrams an included generation. we example ad- cascade, These processes may also involve extra hadronic jets omitted, since the possibleextra diagrams diagrams do are not included for differright-handed BP3 from squarks and those may BP7 decay. for to Additionally, illustrate BP1 theto the and alternative possible BP8 BP7 BP3, routes by decay however, by chain simply which may the switching apply each of the stop/top for the corresponding sbottom/bottom flavoured quark, whilst the right-handed squarkgluino decays and with a 70% quark and or 30%the directly BR into gluino into a always either NLSP a decays and intowith quark, the respectively. the respective In third stop/sbottom both of squark generationtop/bottom these and quark. squark points top/bottom decaying quark, furthermore into a NLSP and corresponding quark. The first- andcorresponding second-generation quark squarks as decay in with BP1 100% and BR BP2. into thegluino NLSP and lighter than theBP3 first and BP4, two the generation left-handed squarks, squarks as always in decay into BP3 a and gluino BP4. and a Much correspondingly like in BP2. However, thesufficiently respective so stop/sbottom-type that squark gluino isassumed two-body now to decays lighter always are than decaying possible. the into top/bottom gluino, quark, In a which these stop in two (BP5) turn BPs or decays the sbottom into gluino a (BP6) is NLSP squark and and correspondingly correspond- flavoured these scenarios the left-handedingly squarks flavoured quark, always decay whilstcombination into or the the skip right-handed this gluino squarkscorresponding step and decay BRs entirely a either of and correspond- 70% into decaymated and directly a by 30%, into gluino-quark NMSDECAY respectively, a [ with NLSP these and fractions a close quark, to with those esti- performed using Delphes [ that gluinos are indecay both first cases into around squarks, 10flavoured GeV/ quark: with each ˜ squark decaying into a NLSP and a correspondingly their respective Branching Ratios (BRs). A fast simulation of the CMS detector is then JHEP10(2018)064 BP5 possible decay cascade. BP3 possible decay cascade. BP7 possible decay cascade. (f) (d) (b) – 6 – BP1 possible decay cascade. BP6 possible decay cascade. BP3 alternative decay cascade. (e) (a) (c) JHEP10(2018)064 , ˜ g 2 ˜ q, c are M 0 1 ˜ χ M − 0 2 = 3 GeV/ ˜ χ 0 1 ˜ χ M M and and ˜ g 2 c M − ˜ q M BP8 possible decay cascade. = 130 GeV/ 0 2 ˜ χ (h) M particles by MadGraph, they must be on-shell. – 7 – ], where , the initial gluino and squark may be produced 15 2a now treated as two independent parameters. 0 2 ˜ χ , 0 1 final state ˜ χ . However, if one were to direct MadGraph to produce two . Thus there is an overlap between the events generated M 3a 3b and ˜ g ˜ q, M heavier in the latter. Therefore, one would find that a mass scan about 2 c ] which contain unique mass hierarchies. Here, BP1, BP3, BP5, BP6, BP7 and BP8 BP7 alternative decay cascade. 15 Feynman diagrams showing example processes by which we may produce two singlino (g) However, it may be noted that BP1 and BP2 are essentially the same, but for In order to remove the possibility of subprocesses being counted more than once in the Considering the diagram in figure being 400 GeV/ BP1 would encapsulate BP2however, anyway. BP5 This to is BP8 alsostop transform the squark into case four is for independent involved BP3BP8 in scans, the with the since converse respect decay for applies. to cascades BP5 BP4, but and BP7 not the the sbottom squark. For BP6 and 4.1 Two-dimensional mass scans Considering each of the existing BPsa in two-dimensional [ mass scan iskept constructed. constant, with The mass gaps Thus, by summing the complementaryover on- which and the off-shell terms, squarks the and entire gluinos momentum may space decay is obtained. 4 Mass scans and event selections which are in turn decayed by Pythia into squarks. overall calculation, it isparton-level required must that be any off-shell. squarkPythia Since will or have any gluino been squarks whose treated or as decay gluinos is whose performed decay at is performed by of the gluino, shown insquarks, figure allowing the productioninvolve of a extra jets, squark some stemmingat of from the parton-level, the subprocesses as decay considered inby of would figure MadGraph a when gluino, asked all for of two which squarks would and be those calculated generated when producing gluinos, LSP along with twothe Standard BPs Model-like in Higgs [ bosons.are shown An in example (a), diagram (b,c), (d), is (e), given (f,g) for and each (h), of respectively. at parton-level in MadGraph, with the second squark produced by Pythia from the decay Figure 2. JHEP10(2018)064 4 ] . 2 250 250 c 250 250 2 − − − − 0 1 ˜ ˜ ˜ q g ˜ 0 1 q g SM SM M M M M [GeV/ q H H q ect. The initial squarks, b ˜ ↵ = = = = ˜ t, q decoupled decoupled Pythia8 ˜ t ˜ t b b ˜ ˜ 0 2 0 2 M M M M M ] ˜ q 2 c ˜ q ˜ g + 127 + 127 + 127 + 127 + 127 + 127 [GeV/ 0 1 0 1 0 1 0 1 0 1 0 1 ˜ ˜ ˜ ˜ ˜ ˜ χ χ χ χ χ χ 0 2 ˜ χ M M M M M M M MadGraph } } } } } } ] Gluino decays in MadGraph. 2 c 20 20 20 20 200 200 − (b) − − − − − ˜ b b ˜ q ˜ ˜ g P P ˜ ˜ t t 4 M M M M [GeV/ The above Feynman diagram shows an example of one of the processes by which Considering the diagram in Figure 2, the initial gluino and squark may be produced Therefore we ask MadGraph to include additional hard jets at parton level, which The mass hierarchy in a particular point in mass-space dictates the possible decay M M in Table 1, we see that gluinos are around 10GeV heavier than squarks. The gluinos 0 1 ˜ we may produce thisproduce final one state squark of directly, two whilst LSPcorresponding the and to second two a squark Higgs stems mass bosons. fromtwo hierarchy generations. the where In decay the Such this of gluino processes example a isalong may we gluino, heavier with also than the involve extra squarks initial jets sparticles of such or the as typically first hard softer jets jetsat produced parton-level radiated in afterwards. MadGraph, with theof second the squark gluino. produced by However, Pythia if from onethe the were production decay to of direct extra MadGraph jets, to some produce of two squarks, the allowing subprocesses considered would involve a squark CONTENTS then allows a wider spread of sparticleantisquarks momenta, and negating gluinos, this plus e up toat two additional leading hard order, jets, after are which generatedCMS by Pythia8 MadGraph detector [7] fast performs simulation showering, decays is and then hadronisation. performed using Delphescascades we [8]. must consider. Taking forP1 example the first benchmark pointthus in decay [3], denoted first intodecays further, squarks, producing with anjets LSP each are and squark produced a decaying SM-like at Higgs each into decay boson. an step, Appropriately-flavoured as NLSP, shown which below. then Figure 2: Feynman diagram showingsinglino an LSP example along process with by two which Standard we may Model-like produce Higgs two boson. χ → { → { → { → { distribution, the masses of the NLSP and LSP → { → { – 8 – M 3 3 3 3 3 3 ] contains many measurement bins for various ob- 0 1 0 1 5 miss ] T 2 SM SM c H q H H q ect. The initial squarks, ↵ 200 200 200 q Pythia8 0 2 0 2 + 10 + 10 + 10 − − − ˜ ˜ ˜ q q q ˜ ˜ ˜ [GeV/ q q q M M M ˜ g ˜ q M M M Table showing various mass ranges in the scans. ˜ M q ] ˜ g 2 c 3000 3000 3000 3000 3000 3000 Table 2. → → → → → → [GeV/ MadGraph ˜ q Gluino decays in Pythia. 1200 1200 1200 1200 1200 1200 M lower than the lightest of the squark and gluino, such that the gluino may still MadGraph event generation of two squarks and one jet (left) and one squark and one 2 (a) c P P Considering the diagram in Figure 2, the initial gluino and squark may be produced The above Feynman diagram shows an example of one of the processes by which The mass hierarchy in a particular point in mass-space dictates the possible decay Therefore we ask MadGraph to include additional hard jets at parton level, which For the first two scans the NLSP mass is increased up to just below the lighter of the Whilst it is expected that the sensitivity to these NMSSM scenarios be lowest for the in Table 1, we see that gluinos are around 10GeV heavier than squarks. The gluinos BP5 BP6 BP7 BP8 two generations. Such processesalong may with also the involve extra initial jets sparticles such or as typically hard softer jets jetsat produced parton-level radiated in afterwards. MadGraph, with theof second the squark gluino. produced by However, Pythia if from onethe the were production decay to of direct extra MadGraph jets, to some produce of two squarks, the allowing subprocesses considered would involve a squark we may produce thisproduce final one state squark of directly, two whilst LSPcorresponding the and to second two a squark Higgs stems mass bosons. from hierarchy the where In decay the this of gluino example a is we gluino, heavier than squarks of the first Figure 2: Feynman diagram showingsinglino an LSP example along process with by two which Standard we may Model-like produce Higgs two boson. CONTENTS then allows a wider spread of sparticleantisquarks momenta, and negating gluinos, this plus e up toat two additional leading hard order, jets, after are which generatedCMS by Pythia8 MadGraph detector [7] fast performs simulation showering, decays is and then hadronisation. performed using Delphescascades we [8]. must consider. Taking forP1 example the first benchmark pointthus in decay [3], denoted first intodecays further, squarks, producing with anjets LSP each are and squark produced a decaying SM-like at Higgs each into decay boson. an step, Appropriately-flavoured as NLSP, shown which below. then BP1/BP2 BP3/BP4 a NLSP. 4.2 Event selections The considered experimental analysis [ servables in the data, background and signal channels, in particular for the number of scans, where the respective stop250 GeV/ or sbottom squark isdecay into non-decoupled, the its relevant mass third iscases, generation set squark, the to along NLSP with be mass ansquark appropriate may may quark. still In still be these decay increased, in so long an that on-shell the fashion involved into third its generation respectively flavoured quark and lighter of the squark andsparticle gluino. mass This range, is such donea in that range order it to of is cover LSP possible the masses. entirety to of Thus, compare the the the available six same independentsquark NMSSM mass and scenario scans gluino with are masses defined whilst as still in table allowing for on-shell decay. For the remaining four lightest LSP, due to thein suppressed each scan are increased as high as is possible whilst remaining below the mass of the Figure 3. gluino (right), sharing the same Feynman diagram after Pythia has performed the gluino decay. JHEP10(2018)064 T is 4. H T R > ≥ ]. α . Fo- 7 -tagged , noting -jets b miss b T 4 N H pseudo-jets and where T c ) is applied to re- . H ? c . requirement. φ 2 ), 6. Anticipating a high T ! i ≥ H j y -jets b p N jets 400 GeV/ 40 GeV/ 2 =1 i N 4 bins. Additionally, since we X 5, ≥ . . ≥ c

2 ≥ T T , we focus on the upper-most p T − < H | 2 . H η -jets | ! b 3 ] regarding vetoing events containing i j x 5 N 3, or , cut due to the p 400 GeV/ 2 T and ≤ j T T 2 5, =1 . i > α c X M E 2

T -jets = < b 4, – 9 – is implemented as described in table H . − | T = 3 and N 2 η 2 T | α between these two pseudo-jets is minimised [ α ! -tagged hadronic jets ( > ] is 4, i.e., separate bins. 10 GeV/ i b T | 5 j T -jets p b η ≥ | > and E N c T 2 =1 p i -jets X and b

c N = 2, 4 bin there is no binning defined in table v u u t 4 for , and jets. . ≥ 0 20 GeV/ = -jets for events where = 3, b final state, with plenty of jets from both cascades, we consider the miss T T > c -jets N H b 40 GeV/ T ], in the case where an event contains just two hadronic jets, R > M -jets N p ), the number of miss 6 T b -dependent cut on > E T N and jets T + H bin, with the exception of the case where we have four or more p T b N ¯ 5. . c H 0 bb ¯ = 2, 1200 GeV/ b , is the transverse energy of the second-leading hadronic jet and ≥ ≥ 2 ? j -jets -jets T φ T b b E with electrons and muons with jets with N H ∆ At least 6 hadronic jets, where anyN jet must have Finally an As defined in [ In addition, event vetoes are defined such that events will fail should they contain any Furthermore, in this analysis, a cut on ‘biased Delta-phi’ (∆ • • • • • • • • 1200 GeV/ 3 for photons or ∆ . For events with more thanin two hadronic such jets, a the way latter that are the combined difference to in create calculated as follows: where where isolation is defined as0 being separated from other objects by angular distance ∆ that for all but the of the following isolated objects: as follows: > hadronic jets, where the only bin in [ duce the QCDisolated background. leptons Additional and cutstain photons in forward/backward-oriented are hadronic [ also jets. performed, The as event well selection is as therefore failing detailed events which con- cusing on a bin with the highestchance number of of jets, bottom i.e., quarkswe the being consider one mis-identified, for both especially which the ifare the primarily interested Higgs in bosons topologies are generating boosted, high hadronic jets ( JHEP10(2018)064 . 2 c miss miss T T E H ]. 5 ]. How- 70 02 04 08 16 18 32 44 02 5 ...... 0 0 0 0 0 0 0 0 1 and ∞ 0          > -jets 900– b N 46 02 09 25 40 19 62 65 51 ...... 52 2 0 0 0 0 0 0 1 2 . Background Yield 0 ] for the same example 5 > cuts. and the LSP 100 GeV/ 400–900 2 T c α 53 . 4 0 1 0 1 0 2 0 0 0 ] for each of the 5 > 350–400 Data Yield -dependent T 55 ] H . c 0 > 300–350 – 10 – [GeV/ ∞ ∞ ∞ 6 . miss 0 T 200– 900– 600–900 400–600 200–400 900– 600–900 400–600 200–400 H > 250–300 Table detailing the 65 . ] 0 c it is clear that the respective efficiencies of each of the event 5 > 200–250 ∞ ∞ ∞ ∞ ∞ ∞ ∞ ∞ ∞ ]. Table 4. [GeV/ ] 24 c [ T 400– 1200– 1200– 1200– 1200– 1200– 1200– 1200– 1200– H Data and background yields for the bins used in this analysis, taken from [ T bbbb α 1 [GeV/ is a simplified Supersymmetric model whereby pair produced gluinos each un- T . In this example the gluino has a mass of 1900 GeV/ T 4 H contains the data and background yields from [ bbbb 1 4 4b Table 3. As shown in table T These data yields and background estimations are then used to calculate lower bounds Bin ≥ = 3b = 3b = 3b = 3b = 2b = 2b = 2b = 2b in the event generation, as was done with theselections NMSSM are signal all mass within points. abenchmark few percent model. of those Therefore, takensuch in directly as from order Initial [ to State take Radiation (ISR), into the account choice other of sources 25% of appears uncertainty to be appropriate. is chosen, dergo a three-body decay into ain bottom figure quark-antiquark pair and an LSPIn neutralino, addition shown to the pair produced gluinos, up to two hard jets are considered at parton level In order to checkused that one to may implement rely theuncertainty on is event the appropriate, selection signal it and event isever indeed yields important of to that calculated course compare the by the experimental these the estimationscenario analysis yields software in such of with question as those the does in those not systematic [ under feature an consideration NMSSM in low- this work, so a reference benchmark model on the sparticle masses giventhe the systematic signal uncertainty yields on for the each signal mass yields point assumed in to each be mass4.4 scan, 25%. with Validation of cut and count analysis tools 4.3 Signal, background andTable data event yields bins, satisfying the remaining event selection criteria. JHEP10(2018)064 bench- . Plotting bbbb 2 1 c T = 100 GeV LSP more typ- 98.2 98.1 98.1 98.1 92.2 74.4 73.7 73.7 23.7 100.0 0 1 2 χ c M Delphes & own software ] 5 ]. = 1900 GeV, Benchmark Model ˜ q 24 [ M : bbbb 99.4 98.7 93.9 93.9 88.5 69.9 69.3 69.2 25.1 100.0 1 T bbbb T1 Efficiencies from [ – 11 – jets stemming from the decay of each Higgs boson. b 4) . 2 > jet angular separation, since this quantity taking a value | b η LSP and the other a mid-range 953 GeV/ | thresholds 2 c T α values, i.e., the generated values without detector simulation being 25 . truth 1 95 . -dependent < 0 T 5 . Feynman diagram showing gluino pair production and decay in the < Cumulative percentages of events passing the event selections compared with those H 0 200 GeV miss T  1 j h > E 200 GeV > 100 / ] for a standard reference benchmark model, 5 > < f - and > min ∗ First consider the BP1-BP8-type mass scans, taking from each two mass points where miss miss φ T T T 1 1 jet j T . Event Selection Before selection Isolated , , vetos p 0 H H Event veto for forward jets ( H n ∆ Figure 4. mark model. tagging and the ability to resolve both one has the lightestical 3 GeV/ of simplified SUSYquantities models, of choosing interest for an all example events, squark before mass any of event 2 selection TeV/ is applied, and normalising We examine the observablelections properties have of the these greatestconsider BPs, the impact in MC on order the toapplied, experimental explore for quantities which sensitivity. such eventbelow as Additionally se- the resolution we of also the detector will have a large effect on the efficiency of Table 5. from [ 5 Signal properties JHEP10(2018)064 T H BP5 BP8 (f) (c) ] will likely thwart most 5 , due to the large number of jets near the observed limit in the BP1- T H . 2 LSP c M regions. These are then overlayed with the BP3 0 1 BP7 for both QCD multijet and top pair produc- ˜ χ T – 12 – M . (e) (b) H c background processes and an MSSM-like scenario with a ¯ t t distributions for two example points from each of the six mass T H T H BP6 BP1 distributions for low and mid-range T shows the and Higgs boson enriched scenario. In this model the squark mass is set at (a) (d) H 5 and the effective “LSP” has mass 3 GeV/ miss T 2 c E Figure Additionally, a simplified MSSM-like model is presented for comparison. This scenario scans. Whilst it istion background clear processes that is thecross-sections far mean are lower much than higher. thatseparated, Therefore, for it whilst the would the be signal peaks expected processes,icant for of compared their the these with tails respective distributions of the are both signal well background processes processes when to considering still the be event signif- yields in an analysis. 5.1 Total scalar A dominant feature ofproduced many in SUSY the cascades many isdistribution decays. high with We mean note well that over these 2 TeV/ low-LSP mass scenarios deliver an with the squarks and gluinosare decaying produced instead and directly the to LSPa the momentum low is LSP. no Thus longer no2 suppressed, TeV/ Higgs removing bosons the possibility of of the signal eventscorresponding generated normalised in distributions the for low- backgroundevents processes and containing background QCD contributions multijet from pair production. is derived from the BP1-type scan, with the difference being that the NLSP is dropped, Figure 5. type scan, compared withlight LSP. QCD and to unity for comparison, it can be shown how many of the cuts in [ JHEP10(2018)064 ¯ ν re- ¯ νν T ν H → ∗ ZZ is quite high. → minimum miss T H c H , it becomes clear that BP8-type mass scan BP5-type Mass Scan is much lower where the 6 (f) T (c) H values due to the presence of miss T H of the jet(s) emitted as the squark or T of any emitted quark is suppressed. p 1200 GeV for the BP1-BP8-type mass scans. T p , that the mean > 6 T – 13 – H BP7-type mass scan BP3-type Mass Scan (e) (b) can arise via Higgs bosons decaying to final states other than a bottom rather distinct distributions between the two mass points may be observed. miss T Fraction of events with total 7 H BP6-type mass scan miss T BP1-type Mass Scan H (d) (a) Even in these light-LSP scenarios, there are still events whose In figure However, it is also apparent, in figure In fact, considering the fraction of events passing the 1200 GeV/ Figure 6. NLSP decay, leaving a soft LSP. This larger quark-antiquark pair, though the branching fractionsare for very decays small. such as A more likely phenomenon is for one or more of the bottom stem- back-to-back. The heavy LSP scenariossuch give a a heavy, fairly boosted widesince and spread invisible the of LSP. considerably heavier The Higgs light-LSP scenarios boson however will suppress inherit this, most of the momentum from the in the decay cascade are reduced and so the 5.2 Whilst of course thealso existence be of some these missing LSP net with non-zero transverse momenta energy, dictates so there long must that the two LSPs are exactly not is not the caseof are the generally squarks limited and togluino decays. gluinos, those where thus the reducing LSP the mass is closeLSP mass to approaches the that masses of the squarks and gluinos. This is the case since the mass gaps quirement across the BP1-typefor mass many scan range, of shown the in mass figure points essentially all events pass this cut. The regions where this JHEP10(2018)064 and c ]. It is clear 5 40 GeV/ BP5 BP8 > T , rises considerably p (f) (c) c cuts, thus decreasing miss T ], with H 5 , the fraction of events with 8 bin, 200 GeV/ near the observed limit in the BP1-type requirement used from [ miss T LSP H miss T M H BP3 BP7 decay channel, meaning a larger contribution – 14 – (e) (b) γγ → . H T . Additionally, hadronic Higgs boson decay final states H background processes. , the minimum c miss ¯ t T t and H . miss miss T T H H decreases drastically in the limit of a light-LSP in the NMSSM, with distributions for low and mid-range miss BP6 T BP1 H shows the number of hadronic jets as defined in [ miss T (a) (d) are difficult to reconstruct owing to resolution and, given the higher number of 9 H b ¯ b 4, for the usual signal and background processes. . greater than the lower edge of the minimum 2 Figure Additionally, as seen in the colour map plots in figure Here the < miss T | η | Typical SUSY decay cascades often involveThe a large model number of considered hadronic jetsdecays, in being and produced. this some paper includingmore features squarks. stop two and such sbottom-type cascades squarks ending whose in decays Higgs produce boson even as the LSP mass increases.since, Above a as certain the threshold, LSP however, thisdecay mass fraction cascade approaches begins shrinks, that to drop meaning ofenergy, that the and we squarks thus are low and left gluinos, with the heavy mass5.3 LSPs gaps with in very the Number little of kinetic hadronic jets the peak well below 200for GeV/ these areas of massexperimental space sensitivity many of to the this events type are of lost model. due to H contributing to the overall such as jets compared with, say, ato the cleaner uncertainty in Figure 7. scan, compared with QCD and ming from the decaysor of muon, the along Higgs with bosons the to appropriately-flavoured decay , into with a the charm neutrino quark momentum and an electron JHEP10(2018)064 BP5 BP8 (f) (c) BP8-type mass scan BP5-type Mass Scan (f) (c) near the observed limit in the LSP M for the BP1-BP8-type mass scans. c BP3 BP7 200 GeV/ > – 15 – (e) (b) background processes and an MSSM-like scenario with ¯ t miss t T BP7-type mass scan BP3-type Mass Scan H (e) (b) BP6 BP1 Fraction of events with (a) (d) Number of hadronic jets for low and mid-range BP6-type mass scan BP1-type Mass Scan (d) (a) Figure 8. Figure 9. BP1-type scan, compared with QCDa and light LSP. JHEP10(2018)064 of some of these jets T 5 for the BP1-BP8-type BP8-type mass scan BP5-type Mass Scan p > (f) (c) – 16 – BP7-type mass scan BP3-type Mass Scan (e) (b) minimum threshold and so will not be considered. c -tagged hadronic jets b Fraction of events with total number of hadronic jets that whilst in general most signal events pass this selection, the three scans BP6-type mass scan 10 BP1-type Mass Scan (d) (a) Additionally in each mass scan the fraction of events passing this selection is gen- This three-body gluino decay means that each cascade, of which these scenarios include The primary reason for this behaviour lies in the decay products of the sparticles in Considering the fraction of events containing at least six hadronic jets, it can be seen will fall below the 40 GeV/ 5.4 Number of In this NMSSM scenarioproduced. where The the BRs LSP for is SM-like a Higgs singlino, boson two decay SM-like dictate Higgs that bosons the will most be likely decay is event in the detector. erally highest for asquark/gluino lighter mass. LSP, This dropping dropin considerably in the as efficiency decay the for cascades LSP meaning heavy softer LSP mass hadronic is approaches jets, due the such to that the the small mass gaps and a neutralino. Conversely, scenariosthe where gluino the decaying squark into is aand lighter squark only than and one the a quark. gluino quark, involve with the squark decaying intotwo, a produces neutralino an extra quark, thus increasing the expected number of hadronic jets per where the gluino isthis lighter cut. than the squark have a much higherthe efficiency various with mass respect hierarchies. to involve squarks The decaying to decay a cascade gluino for and a the quark, BP3, where BP6 the gluino and decays BP8-type into scans two quarks Figure 10. mass scans. in figure JHEP10(2018)064 . 11 -tagging b . -tagged jets for BP5 BP8 b tagged - b (f) (c) near the observed limit LSP M -tagged hadronic jets, as shown in b -tagged jets. However, the efficiency b -tagged hadronic jets is not as clear as b background processes and an MSSM-like scenario BP3 BP7 ¯ t . t as being a bottom quark, or – 17 – -tagged jets are necessarily hadronic jets contain- (e) (b) b 11f and tagged 11e , say, there is generally a decrease in efficiency as the LSP miss T H -tagged hadronic jets for low and mid-range b for 8 BP6 BP1 -tag multiplicities in these scans stem from the decay cascades, shown in b -tagged hadronic jets, it may be noted that in general at least around half of -tagging a bottom quark is around 70%, whereas the likelihood of b b (a) (d) Number of . However it is clear that in the BP7 and BP8-type scans this efficiency increases . In the case of the BP7 scan, up to four top quarks are produced, whose decays 2 12 Additionally, it may be noted that whilst the correlations between sparticle masses The high Considering the efficiency of the binning imposed, taking the fraction of events with Due to uncertainty, of course, not all scan, it is possible toof obtain each up Higgs to boson. four bottom quarks without even consideringand the the decay fraction ofwas events containing seen at in least figure two to almost 100%, shownshown also for by these the scans large in number figures of such jets infigure the example mass points may lead to the production of bottom quarks. More simply, in the case of the BP8-type each of the example mass points in each mass scanat is least two quite large, asthe shown events in in figure thesefigure signal points contain at least two hadronic jets in each event will be ing bottom quarks, norof will correctly all bottom quarksa form lighter flavour quark is only 1% or so. As such the average number of Figure 11. in the BP1-type scan, comparedwith with QCD a and light LSP. that to a bottom quark-antiquark pair, therefore it is expected that a large number of the JHEP10(2018)064 is φ , since this would be ∗ . φ BP8-type mass scan BP5-type Mass Scan (f) (c) such as erroneous measure- miss T H ] in order to reduce the background 5 values whereas SM processes should threshold necessary to be considered. without that jet ∗ c φ -tagged hadronic jets greater than or equal b distribution is softer as the mass gaps in . Considering the fraction of events with . T p 13 miss T – 18 – calculated E BP7-type mass scan BP3-type Mass Scan miss T H (e) (b) , is a variable used in [ values is then taken as minimum ∆ ∗ , it becomes clear that a large fraction of events are rejected φ φ 4. For each of these jets the difference in azimuthal angle values produced in these signal mass points indicates this partic- . 14 ], despite this variable being designed to reject QCD background 2 ∗ 5 φ < | ∗ , as shown in figure η ∗ | φ φ and would be expected to have large ∆ c Fraction of events with total number of miss T 5, shown in figure H BP6-type mass scan . BP1-type Mass Scan 0 or 40 GeV/ > (d) (a) The relatively low ∆ The minimum of these ∆ However, the NMSSM scenarios considered in this paper do not generate many events This quantity is defined as follows. For each event, first consider the hadronic jets with > ∗ φ miss T T in this paper. ∆ by the event selection in [ events and to allow events with genuine ular background-reduction variable is not tuned to the type of scenario under consideration expected to be mostment robust of against the sources momentum of of spurious hadronic jets. containing large ∆ E generate small values, typically less than 0.5. p calculated between the jet and the transverse momentum required to pass the 40 GeV/ 5.5 Minimum ∆ Minimum biased Delta-phi, ∆ contribution from QCD multijet events, designed in such a way that events with genuine Figure 12. to 2 for the BP1-BP8-type mass scans. mass approaches the masses ofof the squark hadronic and jets, gluino. this Similarlythe to is SUSY the likely case cascade since for decrease, the the number meaning jet some of the bottom quarks will not have enough JHEP10(2018)064 BP5 BP8 (f) (c) BP8-type mass scan BP5-type Mass Scan (f) (c) near the observed limit in the BP1- LSP M 5 for the BP1-BP8-type mass scans. . 0 > BP3 BP7 ∗ φ – 19 – (e) (b) BP7-type mass scan BP3-type Mass Scan background processes and an MSSM-like scenario with a ¯ t t (e) (b) Fraction of events with ∆ distributions for low and mid-range BP6 BP1 ∗ φ ∆ (a) (d) BP6-type mass scan BP1-type Mass Scan Figure 14. (d) (a) Figure 13. type scan, compared withlight LSP. QCD and JHEP10(2018)064 5 as . 5, since the . how for ex- 0 , defined as the 15 > µ ∗ φ , of which there are c 40 GeV/ cut would indeed kill much of is less than the jet cone radius ∗ > φ R T p bottom quark jet. 4. If two such jets from a boosted Higgs . ] are standard AK4 jets, formed using the 5 fat = 0 – 20 – R lies in the number of hadronic jets produced. Without of value. As the LSP mass increases, however, the bot- ∗ φ R value computed without that jet to exceed the cut of 0 miss T H only involves hadronic jets with ∗ φ 4. . -tagging both hadronic jets from each Higgs boson will become very difficult. b is performed in order to determine the strength parameter algorithm with a jet cone radius 3 T k between the bottom quarks from each Higgs boson decay, at MC generation (“truth”) Here, it may be noted that the behaviour is essentially the same for each of these Thus, in the extreme light LSP limit, it is expected that a large fraction of the bot- Considering the BP1-BP8-type mass scans it may be seen in figure In order to examine the angular separation of the bottom quark jets stemming from Given the very large average number of jets per event in the light-LSP mass points Accordingly, the regions of parameter space for which the LSP mass is close to that A likely reason for small ∆ R upper limit placed on the fraction of the signal cross-section which could not quite be ruled 6 Results A fit across the signal,in background and table data yields in the nine measurement bins described resolving and six mass scans,boosted since Higgs in boson all and,antiquark cases pair. as heavy such, quarks small angular combined separation with between a the light bottom LSP quark- leads to a thus decreasing thetom mean quark ∆ jets becomethat of more the separated, associatedHiggs since Higgs boson. the boson, LSP resulting momentum in increases a relative more to boostedtom quarks LSP stemming and from a the less decay boosted of each Higgs boson will overlap to the extent that of all individual particles inthe the hadronised event, jets, rather allowing thancone for the radius measurements emulated of of detector 0 angular measurement separation of smaller than theample jet a heavy squark and a light LSP will correspond to a more boosted Higgs boson, anti- boson are close enoughthen that it their might angular only separation be possible ∆ to resolve one the Higgs bosons, the MC truth information is examined. This contains the four-momenta 5.6 Angular separationOne between such quantity bottom which characterises quark the BPs jets∆ in from this signal Higgs model is boson thelevel. decays angular separation The bottom quark jets considered in [ of the squarks andcalculation gluino of display ∆ afewer larger in these fraction regions. of events with ∆ presented in this paper,the it signal would events in be these expected regions that of a parameter space. ∆ course assuming a uniformbe distribution in expected the to become direction increasinglya of unlikely the hadronic for hadronic jet the jets, and azimuthalthe it the angular number would separation of still between jets in each event increases. JHEP10(2018)064 to 2 = 1 c µ BP5 BP8 (f) (c) . The black contour at near the observed limit in the BP1- µ . 6 LSP M 1 the point may not be ruled out. µ > BP3 BP7 – 21 – (e) (b) -axis represent the masses of the squark and LSP, respectively, , with the systematic uncertainty on the signal yields assumed Y 3 in some cases, as shown in table ] decreases dramatically in the limit of a light LSP in these NMSSM 2 5 c - and distributions for low and mid-range it can be seen that both the observed and expected limits are far weaker BP6 BP1 X R 16 ∆ (a) (d) 1 this point may be excluded, whereas if Additionally whilst the black observed limit contour is generally further to the left In figure The expected limit, shown as a red contour, is defined as the upper limit we would These limits are calculated given the data and background yields as well as background Here we present the contour limit plots for the six scan types performed within the much lower squark masses areof no the longer excluded. analysis in Asscenarios, such [ with it the would appear loweras the bound sensitivity little on as the 1 squark TeV/ mass decreasing from aroundthan 2 the TeV/ red, expected limit contour, indicating a slight excess in some of the data yields observe if the dataquantify yields any matched excess the or deficit backgroundthe in expectation background-only the hypothesis. exactly, observation and compared allows to us what we to wouldfor expect given a very light LSP mass, with both contours bending to the left such that points with while the colourthen scale identifies shows the the areasand strength in outside parameter, parameter which space maysignal inside estimations. not which be all excluded, mass given points the are data ruled yield out, and the background and uncertainties in table as 25%. NMSSM. The Figure 15. BP8-type mass scans. out at 95% Confidenceµ Level < (CL). As such, if for a particular point in parameter space JHEP10(2018)064 -axes Y - and X BP6-type mass scan BP8-type mass scan BP3-type Mass Scan (f) (d) (b) – 22 – . µ BP7-type mass scan BP1-type Mass Scan BP5-type Mass Scan (e) (c) (a) Observed and expected limits for the BP1-BP8-type mass scans. The represent the squark andlimit on LSP the masses, strength respectively, parameter whilst the colour scale represents the upper Figure 16. JHEP10(2018)064 2 c , to the miss T H . Thus, in order 6a BP8 1200 1000 and low T H BP7 1250 1050 . the converse becomes true, ], the observed and expected 2 c 5 16 BP6 1000 1010 lower than the squark mass, rather 2 c scenarios under consideration relies on BP5 1250 1260 range would be required than is considered miss T E T scenarios compared with the simplified models H – 23 – BP3 1200 1000 in these heavier LSP regions. miss T E miss T BP1 H 1000 1010 ] is not optimised, decreases dramatically. ] ] 5 2 2 shows the sensitivity to the simplified MSSM-like scenario c c 17 . 6 [GeV/ [GeV/ Scan min min LSP. ˜ ˜ q, g, higher, and the BP5/BP6- and BP7/BP8-type scans are the same as the 2 ]. c M M 2 5 c lighter than the squark/gluino, rather than being decoupled. 2 Approximate lower bounds on the squark mass and corresponding gluino mass at 95% c For comparison, figure The observed and expected limits for these remaining mass scans exhibit a similar Similar experimental limits can be placed on the other types of scan. We recall here However, as the LSP mass is increased closer towards the masses of the squarks and It may be noted that the lower bounds on the squark and gluinoAs masses are the consider- LSP mass is increased above around 100 GeV/ These weaker lower bounds on the squark masses for the lightest LSP mass of 3 GeV/ cut, since few events in this region pass this cut as shown in figure T the LSP being singlino.in In an this NLSP case it decaying is to possible an for LSP the and decay aseen cascades Higgs to previously. boson, end which Applying exclusively is thelimits not same are the event calculated case selections at in 95% as the CL, in MSSM. akin [ to the limits in figure sensitivity for regions with lowerlatter LSP of masses which and the featuring analysis high in [ 6.1 MSSM-like scenarios withThe light main LSP feature of the light LSP and low BP1- and BP3-types, respectively, but with250 GeV/ the appropriate stop or sbottom squark masses behaviour, that is, themid-range LSP cross-section mass, appears where to the dominate contours the are sensitivity closer for to vertical. points with However, in all cases the to explore this area ofin mass this space paper. a wider that the BP3-type scan has thethan gluino 10 mass GeV/ 200 GeV/ This is expected due to the larger gluinos, the sensitivity appearslight once LSP again region, to however, decrease thisH lack for of heavier sensitivity neutralinos. for heavy Unlike LSP the likely arises from the high ably weaker for these light-LSP,considered low- in [ with the limits being more harsh for these NMSSM scenarios than for the simplified models. compared with the backgroundsonably estimation, strong. the agreement between the two limits isare summarised rea- in table Table 6. CL for a 3 GeV/ JHEP10(2018)064 0 1 ], ˜ χ 15 M , with 6 data from 1 − 9 fb . or so. However, below 2 c ] with 35 5 mass scenarios. Conversely, , however, this would require . The eight original BPs in [ miss T miss T 2 c H E scenarios present further challenges for miss T it may be observed that the experimental E -tagged jets per event are lower. Therefore, it b – 24 – 17 ] would be expected to increase the sensitivity to 5 have been placed on the squark/gluino masses for a 2 c ] at the LHC. In order to develop a search for these and 5 16 LSP, despite having a large direct production cross-section, 2 c , contrary to the NMSSM-specific low- 2 c is weaker for this model compared with the NMSSM scenarios considered. 2 c Observed (black) and expected (red) limits for an MSSM-like scenario demonstrating based searches akin to [ miss T E It is clear that these light LSP and low In the limit plots in figures 200 GeV/ -tagged jets, so they will not pass the event selections which are imposed in this paper. the CMS detector at the LHC, thus they cannotjets+ be completely excluded atstealthy this scenarios stage. one might wish tocareful access regions techniques of so low as to not allow yields from background processes to dominate. singlino-like LSP in the NMSSMthis for mass a the limits LSP weaken mass bylimits above a for 100 considerable some GeV/ amount scenarios in decreasing all byall as cases, featuring much as as a shown 1 in light TeV/ table are 3 GeV/ still on or around the limit of sensitivity of the analysis in [ this MSSM-like scenario. 7 Conclusion Fairly strong limits of around 2 TeV/ Since the decay cascadethe is expected truncated numbers and of there hadronicis are jets unlikely and no for Higgs as bosonb many decays events in to this containHowever, more model, exploiting than the full five 254 hadronic bins jets in or [ as many as two sensitivity to this simplifiedmass of MSSM-like 3 GeV/ modelit is is in also fact clear that> strongest the for overall sensitivity the for lightest the regions LSP of parameter space in which Figure 17. the higher sensitivity to regions with low LSP mass. JHEP10(2018)064 T ]. k 27 B 70 13 -tagging (2018) b jets. One -tagged jet (1984) 1 ) sensitivity 05 b 11 110 tagged Nucl. Phys. ] whereby, in the - , b JHEP ]. This might lead 26 , 16 . Phys. Rept. 1 , H ]. Using the double- Pisma Zh. Eksp. Teor. Fiz. 16 [ ]. 1 H SPIRE IN -tagger, one of which being developed by (1971) 323 [ b [ 13 – 25 – Gauge Models with Spontaneously Broken Local 8. This tool thus allows for a way to distinguish -tagged jet mass can be studied. In any scenario . Extension of the Algebra of Poincar´eGroup Generators b (1982) 343 ]. JETP Lett. ), which permits any use, distribution and reproduction in , B 119 SPIRE -tag multiplicities in these scenarios (cf. figure b IN Supergauge Transformations in Four-Dimensions ][ ] and another by the ATLAS collaboration [ 25 ]. Search for natural and split supersymmetry in proton-proton collisions at ]. CC-BY 4.0 Phys. Lett. , This article is distributed under the terms of the Creative Commons SPIRE SPIRE IN Supersymmetry, Supergravity and IN [ ]. TeV in final states with jets and missing transverse momentum -tag multiplicities occur also in scenarios where the rˆoleof the SM Higgs boson collaboration, b arXiv:1802.02110 = 13 [ SPIRE s IN Supersymmetry [ CMS √ 025 and Violation of p(1971) Invariance 452] [ (1974) 39 High Given the expected high H.P. Nilles, J. Wess and B. Zumino, R. Barbieri, S. Ferrara and C.A. Savoy, Yu. A. Golfand and E.P. Likhtman, [4] [5] [2] [3] [1] Attribution License ( any medium, provided the original author(s) and source areReferences credited. and U.E. acknowledge support fromMinimalHiggs). the H2020-MSCA-RISE-2014 U.E. grant also no.agreement acknowledges 645722 no. support (Non- 690575 from (InvisiblesPlus). the Marie Sklodowska-CurieOpen grant Access. Acknowledgments SM, CHS-T and ATported are by supported an in STFC partST/L000296/1. studentship. through H.F SM the acknowledges acknowledges NExT financial financial Institute, support support from from AT is the STFC CG also STFC CG ST/N000250/1. sup- SM method, distributions ofinvolving the SM double- or NMSSMmass specific should exhibit Higgs excesses bosons,quark over distributions pair a production of smoothly around the decreasingto the double- background a corresponding from simultaneous Higgs discovery QCD of boson and an mass top extra [ lighter Higgs boson two bottom quarks originating from thealgorithm same with decay vertex. a Here jetbetween AK8 single cone refers bottom to radius the quark anti- of jetsfound and originating 0 two from overlapping the bottom decay quark of jets, a such boosted as Higgsis may boson, played be a by technique a employed lighter in NMSSM [ specific Higgs boson may be boosted, though, ifmay one employ concentrates on new signal techniques regionsthe with such CMS more as collaboration a [ double former, AK8 jets areis formed assigned and, by depending a on Boosted the Decision substructure, Tree a (BDT), discriminator describing output how likely the jet is to contain JHEP10(2018)064 ]. ] ]. (2017) (2014) (2017) SPIRE B 138 12 IN 07 ]. ][ Phys. Rev. Comput. SPIRE , (2017) 005 , C 77 IN ][ JHEP 10 SPIRE JHEP , , IN Comput. Phys. [ Phys. Lett. , ]. ]. , arXiv:1307.6346 TeV with the ATLAS [ JHEP , (2017) 032003 SPIRE SPIRE Eur. Phys. J. = 13 TeV IN , IN arXiv:1406.7221 s variable in the all-hadronic [ ]. ][ D 96 ][ √ ]. TeV T2 arXiv:1412.6394 (2014) 057 = 13 [ M s hep-ph/9611232 02 ]. SPIRE , √ SPIRE ]. IN = 13 IN s ][ (2014) 113 DELPHES 3, A modular framework for ][ Phys. Rev. √ SPIRE , JHEP 10 TeV pp collision data with the ATLAS SPIRE IN , (2015) 172 The Next-to-Minimal Supersymmetric IN TeV with the ATLAS detector ][ 04 arXiv:0910.1785 ][ = 13 A Brief Introduction to PYTHIA 8.1 arXiv:1712.02332 [ NMSDECAY: A Fortran Code for Supersymmetric s PROSPINO: A program for the production of JHEP [ – 26 – = 13 , √ ]. s JHEP of √ , 1 ]. ]. (2010) 1 Dijet Searches for Supersymmetry at the LHC NMSSM with a singlino LSP: possible challenges for Excessive Higgs pair production with little MET from arXiv:0710.3820 − ]. [ SPIRE arXiv:1711.01901 [ IN ]. ]. 496 [ SPIRE SPIRE arXiv:1106.5633 (2018) 112001 SPIRE The mu Problem and the Strong CP Problem -jets in proton-proton collisions at IN IN [ arXiv:0806.1049 b IN ][ ][ [ SPIRE SPIRE Search for squarks and gluinos in final states with jets and missing Search for supersymmetry in final states with missing transverse Search for a scalar partner of the top quark in the jets plus missing ][ IN IN (2008) 852 D 97 Search for supersymmetry in multijet events with missing transverse Search for direct production of supersymmetric partners of the top quark Search for new phenomena with the Search strategy for exclusive multi-jet events from supersymmetry at ]. ][ ][ (2018) 107 Phys. Rept. 178 06 The automated computation of tree-level and next-to-leading order , (2012) 774 collaboration, J. de Favereau et al., SPIRE IN [ 183 (2008) 221803 Phys. Rev. JHEP collaboration, collaboration, collaboration, ]. , , collaboration, collaboration, collaboration, collaboration, CMS-PAS-SUS-09-001 101 , arXiv:1405.0301 arXiv:1709.04183 arXiv:1705.04650 [ [ [ SPIRE IN arXiv:1704.07781 arXiv:1707.03316 [ Lett. CMS CMS ATLAS transverse momentum using 36 fb DELPHES 3 fast simulation of a[ generic collider experiment Particle Decays in theCommun. Next-to-Minimal Supersymmetric Standard Model 079 supersymmetric particles in next-to-leading order QCD Phys. Commun. searches for supersymmetry at the LHC (1984) 150 differential cross sections and their matching to parton shower simulations Standard Model squarks and gluinos in the NMSSM momentum in proton-proton collisions[ at 13 TeV CMS in the all-jets final state in proton-proton collisions at 085 CMS final state produced in proton-proton710 collisions at CMS ATLAS momentum and multiple detector ATLAS transverse momentum final state at detector L. Randall and D. Tucker-Smith, D. Das, U. Ellwanger and A.M. Teixeira, W. Beenakker, R. Hopker and M. Spira, T. Sj¨ostrand,S. Mrenna and P.Z. Skands, J.E. Kim and H.P. Nilles, J. Alwall et al., U. Ellwanger, C. Hugonie and A.M. Teixeira, U. Ellwanger and A.M. Teixeira, U. Ellwanger and A.M. Teixeira, [9] [7] [8] [6] [21] [22] [19] [20] [17] [18] [14] [15] [16] [13] [11] [12] [10] JHEP10(2018)064 , ]. , SPIRE . IN ][ (2005) 46 168 ]. PoS(ICHEP2016)1129 arXiv:1712.08501 , [ SPIRE IN ][ SDECAY: A Fortran code for the decays of Comput. Phys. Commun. (2018) 241801 , – 27 – Tagger in Run II b 120 ¯ b → arXiv:1301.2175 ]. H [ ]. SPIRE IN Boosted [ Phys. Rev. Lett. SPIRE Identification of double-b quark jets in boosted event topologies Search for Physics Beyond the Standard Model in Events with Interpretation of Searches for Supersymmetry with simplified Models , IN (2013) 052017 ][ D 88 collaboration, collaboration, collaboration, collaboration, hep-ph/0311167 ATLAS CMS High-Momentum Higgs Bosons andCollisions Missing at Transverse Momentum 13 TeV in Proton-Proton the supersymmetric particles in[ the MSSM CMS Phys. Rev. CMS CMS-PAS-BTV-15-002 M. Muhlleitner, A. Djouadi and Y. Mambrini, [26] [27] [24] [25] [23]