Prepared for submission to JHEP

The Light Gluino Gap

Jared A. Evansa,b and David McKeenc aDepartment of Physics, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA bDepartment of Physics, University of Cincinnati, Cincinnati, Ohio 45221, USA cPittsburgh , Astrophysics, and Cosmology Center, Department of Physics and Astronomy, University of Pittsburgh, Pittsburgh, USA E-mail: [email protected], [email protected]

Abstract: A variety of experimental searches and theoretical efforts have constrained gluinos that undergo an all-hadronic decay with no missing energy, as can arise in R-parity violating or stealth . Although the g˜ → jj decay is robustly excluded, there is gap in current experimental coverage for gluinos with masses between 51–76 GeV that decay into three light-flavor quarks. In this work, we probe this gap with published multi- data from the UA2 experiment. Despite setting the strongest current limit on this region, we find that UA2 data is unable to close this gap for g˜ → jjj decays. In addition to this three-jet gap, we note an additional gap for all- hadronic g˜ → n parton decays with n ≥ 4 for light gluinos (51 GeV < mg˜ . 300 GeV) not covered by the current search program. arXiv:1803.01880v1 [hep-ph] 5 Mar 2018 Contents

1 Introduction1

2 The UA2 Experiment4

3 Limits on gluinos from multi-jets6

4 Discussion7

1 Introduction

The large collider (LHC) has pushed the lower bound on the gluino mass to above 1 TeV for nearly any decay path [1], and for particular decay paths, limits are above 2 TeV; see, e.g., [2]. However, there are three exceptions to this 1 TeV bound. First, any gluino process that involves some exotic detector object requires a dedicated study and has the potential to evade the current search program if that object is not searched for explicitly (e.g., non-isolated [3]). Second, gluinos that decay into all-hadronic final states exhibiting a jet pT hierarchy [1], where two jets are hard and the rest are very soft, appear more QCD-like and can have significantly weaker constraints [4]. Lastly, the topic of this work, very light gluinos with all-hadronic decays, which would evade constraints by being buried within the QCD background.

Independent of any production or decay modes, gluinos with mg˜ < 51 GeV were robustly excluded by Kaplan and Schwartz using jet thrust data from LEP [5], essen- tially a measurement of αs at the Z pole. In a similar spirit, LHC measurements of the electroweak running have been used to place relatively robust constraints on new physics [6]. Although efforts have been made to extend this procedure to αs informa- tion contained within the 3-jet to 2-jet ratio at the LHC [7], these results currently do not place robust constraints on new physics. In the CMS measurement of this ratio [8], a value for αs(mZ ) is determined by varying the PDF input αs(mZ ) parameter and minimizing this χ2 against the data. This relies on SM QCD running for the profile. Although the authors of ref. [7] argue that the PDF profile is minimally impacted by including additional states, it is not obvious that these leading order conclusions can be extended to higher order PDFs. While CMS uses a fixed choice of factorization and renormalization scale, the QCD hard process should in truth be sensitive to multiple

– 1 – scales, and the results cannot be directly used [7]. Moreover, very light gluinos have the potential to contribute directly to the 3-jet to 2-jet ratio, skewing any interpretation that does not contain their production and decays directly. Of course, complementary direct limits would still be valuable even if there were a study that included all aspects necessary to indirectly exclude light gluinos.

All-hadronic decays of gluinos with no appreciable missing energy (E/T ) are achiev- able in R-parity violating (RPV) models [9] or stealth SUSY models [10]. Within these models, the parameters can be chosen so that the gluino would decay consistently to n ≥ 2 colored partons. While the RPV couplings require some structure to evade indi- 00 c c c rect constraints [9, 11], with all other couplings small, the λcdsC D S coupling in the superpotential may be O (1) without conflicting with any current precision constraints, allowing for prompt g˜ → jjj decays (although a complicated story for baryogenesis may be required [12, 13]). Stealth SUSY, on the other hand, introduces a nearly de- generate hidden sector supersymmetric singlet, (S, S˜), that couples to QCD allowing for g˜ → gS˜ → gSa˜ → g(gg)˜a decays, where a˜ is the lightest supersymmetric particle, e.g. the gravitino or axino. If ma˜ < mS˜ − mS  mS, the neutral, detector-stable a˜ will generically contribute missing energy far below the detector resolution, giving again a prompt g˜ → jjj decay. For direct limits on all-hadronic decays, a lot of effort, from both theory and experiment, has been devoted to probing the possible all-hadronic signatures from a variety of RPV and stealth SUSY models [1, 14–28]. Although a pure g˜ → jj decay is not possible, if mg˜ − mS˜ . O (GeV), stealth SUSY would allow for a compressed decay of g˜ → S˜ + {soft} → Sa˜ → (gg)˜a to mimic a 2-jet signal. Fortunately, experimental results at CDF [29], ATLAS [30], and CMS [31, 32] reliably close this possibility up to high masses. The combined exclusion from these searches is summarized in figure1, where the ratio of the cross-section limit to the gluino production cross-section is shown as a function of mg˜. Although this interpretation is simply using the efficiencies from t˜ → jj as a function of resonance mass, the deviations introduced from switching the triplet to an octet should be O (1), at most. Many limits have been placed on pair-produced gluinos that each decay to three- jets by CDF [33], CMS [34–36], and ATLAS [37–39]. The LHC bounds do not begin until 100 GeV [37], and the CDF bound begins at 76 GeV [33]. Between the robust lower bound from LEP of 51 GeV and 76 GeV, where CDF becomes sensitive, there is currently a gap in the g˜ → jjj coverage. More complicated, higher-multiplicity, all-hadronic gluino decays, such as the sim- ple five-parton RPV decay, g˜ → jjχ˜0 → jj(jjj), or, more generally, any n ≥ 4- parton decays that can be achieved in stealth SUSY or hidden valley models [1, 40, 41], have been excluded by black hole [42] and all-hadronic [38] LHC searches for mg˜ &

– 2 –  � →{����}+�→ �� 1

0.100 σexcl

∼ σg 0.010

CDF 0.001 LEP ATLAS (Kaplan, Schwartz) CMS 8 (13) 0 50 100 150 200

∼ mg

Figure 1. Limits on stealth decays of light gluinos into (effectively) two jets, g˜ → jj + soft. The robust lower bound of 51 GeV from LEP jet data is in gold [5]. The CDF limit is in red [29], the ATLAS limit is in green [30], and the CMS limits from their 8 TeV direct search [31] and 13 TeV substructure search [32] are in dark and light blue, respectively.

300 GeV [1]. Below ∼ 300 GeV limits weaken, as the significant boost required for the event to pass the trigger causes the decay products to be too collimated to be resolved as separate jets.

In this work, we confront one facet of the light gluino gap, namely the RPV decay of g˜ → jjj within the 51–76 GeV band, by examining multi-jet data from UA2 [43]. In the next section, we will describe the UA2 detector and our UA2 simulator. In order to validate this digital simulacrum, we reproduce the efficiencies from the only hadronic SUSY search at UA2 [44]. We then describe how we reinterpret the multi-jet measurement to place limits on light, R-parity violating gluinos. While these results are close to exclusion, they are unable to close this gap for gluinos. We conclude by discussing how the light gluino gap could be probed with targeted substructure-based LHC searches, likely closing the gap not only for three-jet decays, but also for n ≥ 4- parton decays.

– 3 – 2 The UA2 Experiment

UA2 was one of the two primary experiments at the CERN pp¯ SPS collider. Together with UA1, the collaboration co-discovered the W [45] and Z bosons [46] (see [47] for more history). Years later, when the two studies [43, 44] that are the focus of this work √ were performed, the UA2 detector had collected ∼ 7.5pb−1 of s = 630 GeV data, covered nearly 4π of solid angle, displayed exceptional calorimetry, and possessed no magnetic field. UA2 consisted of a central calorimeter covering |η| < 1 with cells of polar and az- imuthal widths of ∆θ = 10◦ and ∆φ = 15◦ , respectively. The two end cap calorimeters covered 1 < |η| < 3, with cell widths of ∆η = 0.2 and ∆φ = 15◦ , out to |η| = 2.2 where the final cells had ∆φ = 30◦ and 2.2 < |η| < 2.5 and 2.5 < |η| < 3.0. The cells contained multilayer lead- and iron-scintillator sandwiches for electro- magnetic and hadronic calorimetry, respectively. The calorimeter exhibited an energy 3 E  4 resolution σE ∼ 0.32 GeV GeV, which we apply to the cells with a truncated gaus- sian about the truth-level value determined by the showering in Pythia 8 [48]. A would deposit up to 1.5 GeV of energy during its minimally ionizing traverse, with any remaining energy appearing as E/T . As our signal of interest does not contain electrons at appreciable rates, energy from these particles is lumped into jets for simplicity. After cell contents are recorded by the the detector, jets are formed by dropping all cells with E < 400 MeV and combining adjacent cells above this energy threshold until no cells remain. Unlike in some other UA2 results [47], cell clusters are not split in either search due to the presence of energetic valleys [49]. The energy and size of Pcells the jet are defined to be E = i Ei, and

cells q X Ei R = R2 + R2 where R = R , (2.1) φ θ x x,i E i where Rφ(θ),i is the angle corresponding to the center of the cell. In standard fash- ion, missing transverse energy is determined by performing a vector-like sum of the calorimeter cells, and projecting onto the x − y plane. In order to validate our simulation of the UA2 detector, we recast the UA2 search for squarks and gluinos [44] using their benchmark models. The search requires at least two central jets with ET,1(2) > 25(15) GeV and |η| < 0.85 that do not behave like QCD dijet events, and E/T > 40 GeV not aligned with jets in the event. Event samples are generated and showered in Pythia 8 [48]. As many of the simulation details provided by UA2 are antiquated, we focus on the efficiency between our recast and their model. These results are presented in figure2, using data from figure 4 of ref. [44] for three 5 benchmark mass curves: mg˜ = mq˜, mq˜ = 200 GeV, and mg˜ = 10 GeV (in all cases

– 4 – mg∼ = mq∼ mq∼ = 200 GeV mg∼ = 100 TeV 2 5 mq˜ = mg˜ mq˜ = 200 GeV mg˜ = 10 GeV

1 ) pb

( 0.5 σ⨯ϵ

0.2

(a) (b) (c) 0.1 90 100 110 120 60 70 80 90 60 70 80

mg∼ (GeV) mg∼ (GeV) mq∼ (GeV)

Figure 2. Recast of the UA2 SUSY search [44]. The excluded effective cross-sections are shown as functions of (a) mg˜ = mq˜, (b) mg˜ with mq˜ = 200 GeV, and (c) mq˜ with mg˜ = 105 GeV. In blue are the curves from UA2. Our simulation rescaled to their cross-section is shown in orange. In green, we show our simulation with an FSR enhancement (obtained by increasing αs(mZ ) to 0.25) to better mimic the aggressive FSR profile used in the UA2 study.

mχ˜0 ≈ 0). Although we were unable to reproduce the quoted values for their SUSY cross-sections using any tools at our disposal, we rescale our cross-section values to match the values presented in Table 1 of [44] via

σUA2(200, 80) σuse(200, mg˜) = σsim(200, mg˜), (2.2) σsim(200, 80) and similarly for the two other benchmark lines in the study. While we do not expect this approximation to be reliable far from the central point, we believe it is sufficient for determining the reliability of our simulation. In Pythia, we disable multiple parton in- teractions, initial state radiation, and final state radiation (FSR) in process, but enable FSR in resonances to reflect their treatment of the Monte Carlo SUSY model. The UA2 result is shown in blue, with our simulation in orange. Overall, the agreement is rather good, with the exceptions of the 50 GeV squark point (in figure2c), where the efficiency is very low, and in the entirely gluino sample (figure2b) where disagreement is about 40% for all mass choices. In this latter case, we believe this is due to a disagreement

– 5 – between FSR modeling in modern Monte Carlo tools and the fragmentation scheme used by UA2 [50], with the latter being overly aggressive. To support our suspicions, we reran our simulation with αs(mZ ) set to the highest value allowed in Pythia (0.25) modifying the efficiencies resulting in the green curve in figure2. This value more accurately reflects the gluino case, with only small changes to the other cases. Further refinements to our archaeological efforts would not be valuable, as our simulation is mimicking the UA2 detector to the precision we require. In the main study of this paper, we will use the default showering settings, as the modern understanding of pp¯ interactions will better reflect the actual environment in SPS collisions.

3 Limits on gluinos from multi-jets

The UA2 search of interest is a measurement of multi-jet events and search for double- parton scattering [43]. The search region of interest in this study requires six jet clusters

(as defined in the previous section), each with ET,j > 15 GeV and |η| < 2.0, at least 40% of the total transverse energy in the calorimeter to be contained within these six jets, the average cluster radius of the jets, as defined in (2.1), to be less than Rav < 0.349, no additional jet clusters with ET,j > 10 GeV, and small missing energy, E/T < 20 GeV. Within the search region, UA2 observed only 7 of these six-jet events with their 7.6 pb−1 of data. Although this study is intended as a measurement of standard model properties, the signature from prompt, all-hadronic, R-parity violating decays of gluino pairs into three quarks, g˜g˜ → (jjj)(jjj), can actually overwhelm the QCD background so much that pp¯ → g˜g˜ → (jjj)(jjj) events alone would produce many more than the seven events UA2 observed in their signal region. We produce our pp¯ → g˜g˜ → (jjj)(jjj) signal events in MadGraph 5 [51], shower in Pythia 8 [48], and process through our custom UA2 simulator (detailed in the previous section) before applying the specific cuts used in this search. The SPS gluino cross-sections were determined with Prospino2.1 [52, 53] (modified slightly to treat SPS energies). We do not match additional jets, which is justified as hard radiation from the parton shower providing more than 15 GeV, or ∼ 2% of the total collision energy would be extremely rare. While there are several known sources of uncertainty, including the luminosity, our Monte Carlo statistics, and the PDF, the largest by far is our UA2 modeling uncertainty. To this end, we conservatively assign a 33% uncertainty on our signal for all choice of gluino masses, implying that signals yielding more than N95,excl = 19.6 events (σexcl = 2.58 pb) are excluded. In figure3, we show the cross-section as a function of gluino mass excluded by this search in blue, normalized to the gluino production cross-section value at that mass. We find that UA2 constrains the gluino production

– 6 – �� →(���)(���) 5

1 σexcl 0.50 ∼ σg

0.10 UA2 LEP CDF 0.05 (Kaplan, Schwartz) ATLAS 0 50 100 150 200

∼ mg

Figure 3. Limits on RPV decays of light gluinos, g˜ → jjj. The recast performed in this work of the UA2 study [43] is in blue, the CDF limit is in red [33], the lowest LHC constraint from ATLAS is in green [37], and the robust lower bound of 51 GeV from LEP jet data is in gold [5]. cross-section for R-parity violating decays to within 2-3 times the expected values across the entire light gluino to three-jet gap.

4 Discussion

In this work, we illustrated that six-jet event data from UA2 constrains the allowed cross-section for all-hadronic, three-body decays of gluinos in R-parity violating SUSY with masses from 51–76 GeV, between the sensitivity of LEP and CDF. While unable to close the light gluino to three-jet gap, the UA2 data does exlcude gluino production cross-sections that are a factor of 2–3 larger than the expected values, setting what is presently the strongest bound within this region. As discussed in the introduction, high-multiplicity, all-hadronic gluino decays, g˜ → n partons with n ≥ 4, have been excluded for mg˜ & 300 GeV by LHC searches [1, 38, 42]. However, these states are not currently robustly constrained at lower masses, as the decay products become too collimated for the existing multi-jet searches to have

– 7 – sensitivity. Further LHC studies, such as dedicated substructure-based searches within dijet events, have the potential to close this gap for all n-parton decays. Jet substructure techniques (for a recent review and original references, see [54]) have been successfully used to reduce large QCD backgrounds in a variety of searches. These method have, for instance, allowed CMS to search for new bosons with masses as low as 50 GeV that decay into dijets [55], and to search for paired dijet resonances with masses as low at 80 GeV [32] (displayed in figure1). In addition to probing g˜ → n partons, such substructure-based searches would have the potential to constrain light squark multi-jet signatures, such as q˜ → jχ˜0 → j(jjj) in RPV or q˜ → jS˜ → j(gg)˜a in stealth SUSY, or all-hadronic decays of non-SUSY particles such as colorons or axigluons, see e.g. [56– 58]. Direct probes of light hadronically-decaying particles at the LHC are desirable, due the relatively small theoretical and modeling uncertainties, and timely, as rising trigger thresholds will make accessing these regions more difficult.

Acknowledgments

We thank J. P. Chou, S. D. Ellis, G. Gollin, D. Goncalves, E. Halkaidakis, J. Incandela, A. Larkoski, S. Mukhopadhyay, A. Nelson, H. Russell, D. Shih, and S. Thomas for useful discussions. We thank J. Shelton for collaboration in the early stages of this work. We are especially grateful to Karlheinz Meier for providing specific details concerning the jet algorithms at UA2. The work of JAE is supported in part by DE-SC0015655. The work of DM is supported by PITT PACC through the Samuel P. Langley Fellowship.

References

[1] J. A. Evans, Y. Kats, D. Shih and M. J. Strassler, Toward Full LHC Coverage of Natural Supersymmetry, JHEP 07 (2014) 101,[1310.5758]. [2] ATLAS collaboration, M. Aaboud et al., Search for squarks and gluinos in final states √ with jets and missing transverse momentum using 36 fb−1 of s=13 TeV pp collision data with the ATLAS detector, 1712.02332. [3] C. Brust, P. Maksimovic, A. Sady, P. Saraswat, M. T. Walters and Y. Xin, Identifying boosted new physics with non-isolated leptons, JHEP 04 (2015) 079,[1410.0362]. [4] CMS collaboration, V. Khachatryan et al., Search for new physics in events with high jet multiplicity and low missing transverse momentum in -proton collisions at √ s = 8 TeV, Submitted to: Phys. Lett. (2016) , [1608.01224]. [5] D. E. Kaplan and M. D. Schwartz, Constraining Light Colored Particles with Event Shapes, Phys. Rev. Lett. 101 (2008) 022002,[0804.2477].

– 8 – [6] D. S. M. Alves, J. Galloway, J. T. Ruderman and J. R. Walsh, Running Electroweak Couplings as a Probe of New Physics, JHEP 02 (2015) 007,[1410.6810]. [7] D. Becciolini, M. Gillioz, M. Nardecchia, F. Sannino and M. Spannowsky, Constraining new colored matter from the ratio of 3 to 2 jets cross sections at the LHC, Phys. Rev. D91 (2015) 015010,[1403.7411]. [8] CMS collaboration, S. Chatrchyan et al., Measurement of the ratio of the inclusive 3-jet √ cross section to the inclusive 2-jet cross section in pp collisions at s = 7 TeV and first determination of the strong coupling constant in the TeV range, Eur. Phys. J. C73 (2013) 2604,[1304.7498]. [9] R. Barbier, C. Berat, M. Besancon, M. Chemtob, A. Deandrea et al., R-parity violating supersymmetry, Phys.Rept. 420 (2005) 1,[hep-ph/0406039]. [10] J. Fan, M. Reece and J. T. Ruderman, Stealth Supersymmetry, JHEP 1111 (2011) 012, [1105.5135]. [11] Y. Kao and T. Takeuchi, Single-Coupling Bounds on R-parity violating Supersymmetry, an update, 0910.4980. [12] P. W. Graham, D. E. Kaplan, S. Rajendran and P. Saraswat, Displaced Supersymmetry, JHEP 07 (2012) 149,[1204.6038]. [13] Y. Cui, Natural Baryogenesis from Unnatural Supersymmetry, JHEP 12 (2013) 067, [1309.2952]. [14] J. Fan, M. Reece and J. T. Ruderman, A Stealth Supersymmetry Sampler, JHEP 1207 (2012) 196,[1201.4875]. [15] C. Brust, A. Katz and R. Sundrum, SUSY Stops at a Bump, JHEP 08 (2012) 059, [1206.2353]. [16] J. A. Evans and Y. Kats, LHC Coverage of RPV MSSM with Light Stops, JHEP 1304 (2013) 028,[1209.0764]. [17] D. Curtin, R. Essig and B. Shuve, Boosted Multijet Resonances and New Color-Flow Variables, Phys. Rev. D88 (2013) 034019,[1210.5523]. [18] D. Duggan, J. A. Evans, J. Hirschauer, K. Kaadze, D. Kolchmeyer, A. Lath et al., Sensitivity of an Upgraded LHC to R-Parity Violating Signatures of the MSSM, 1308.3903. [19] Y. Bai, A. Katz and B. Tweedie, Pulling Out All the Stops: Searching for RPV SUSY with Stop-Jets, JHEP 01 (2014) 040,[1309.6631]. [20] J. A. Evans and Y. Kats, LHC searches examined via the RPV MSSM, PoS EPS-HEP2013 (2013) 287, [1311.0890]. [21] J. A. Evans, A Swarm of Bs, JHEP 08 (2014) 073,[1402.4481].

– 9 – [22] P. W. Graham, S. Rajendran and P. Saraswat, Supersymmetric crevices: Missing signatures of R -parity violation at the LHC, Phys. Rev. D90 (2014) 075005, [1403.7197]. [23] J. Fan, R. Krall, D. Pinner, M. Reece and J. T. Ruderman, Stealth Supersymmetry Simplified, JHEP 07 (2016) 016,[1512.05781]. [24] A. Monteux, New signatures and limits on R-parity violation from resonant squark production, JHEP 03 (2016) 216,[1601.03737]. [25] M. R. Buckley, D. Feld, S. Macaluso, A. Monteux and D. Shih, Cornering Natural SUSY at LHC Run II and Beyond, JHEP 08 (2017) 115,[1610.08059]. [26] S. Diglio, L. Feligioni and G. Moultaka, Stashing the stops in multijet events at the LHC, Phys. Rev. D96 (2017) 055032,[1611.05850]. [27] H. Beauchesne, K. Earl and T. Gregoire, LHC phenomenology and baryogenesis in

supersymmetric models with a U(1)R baryon number, JHEP 06 (2017) 122, [1703.03866]. [28] D. Dercks, H. Dreiner, M. E. Krauss, T. Opferkuch and A. Reinert, R-Parity Violation at the LHC, Eur. Phys. J. C77 (2017) 856,[1706.09418].

[29] CDF collaboration, T. Aaltonen et al., Search for Pair Production of Strongly √ Interacting Particles Decaying to Pairs of Jets in pp¯ Collisions at s = 1.96 TeV, Phys. Rev. Lett. 111 (2013) 031802,[1303.2699].

[30] ATLAS collaboration, M. Aaboud et al., A search for pair-produced resonances in √ four-jet final states at s=13 TeV with the ATLAS detector, 1710.07171.

[31] CMS collaboration, V. Khachatryan et al., Search for pair-produced resonances √ decaying to jet pairs in protonÐproton collisions at s = 8 TeV, Phys. Lett. B747 (2015) 98–119,[1412.7706].

[32] CMS collaboration, Search for low-mass pair-produced dijet resonances using jet substructure techniques in proton-proton collisions at a center-of-mass energy of √ s = 13 TeV, CMS-PAS-EXO-16-029 (2016).

[33] CDF collaboration, T. Aaltonen et al., First Search for Multijet Resonances in √ s = 1.96 TeV pp¯ Collisions, Phys. Rev. Lett. 107 (2011) 042001,[1105.2815].

[34] CMS collaboration, S. Chatrchyan et al., Search for Three-Jet Resonances in pp √ Collisions at s = 7 TeV, Phys.Rev.Lett. 107 (2011) 101801,[1107.3084].

[35] CMS collaboration, S. Chatrchyan et al., Search for three-jet resonances in pp √ collisions at s = 7 TeV, Phys.Lett. B718 (2012) 329–347,[1208.2931].

[36] CMS collaboration, S. Chatrchyan et al., Searches for light- and heavy-flavour three-jet

– 10 – √ resonances in pp collisions at s = 8 TeV, Phys. Lett. B730 (2014) 193–214, [1311.1799].

[37] ATLAS collaboration, G. Aad et al., Search for pair production of massive particles √ decaying into three quarks with the ATLAS detector in s = 7 TeV pp collisions at the LHC, JHEP 12 (2012) 086,[1210.4813].

[38] ATLAS collaboration, G. Aad et al., Search for massive supersymmetric particles √ decaying to many jets using the ATLAS detector in pp collisions at s = 8 TeV, Phys. Rev. D91 (2015) 112016,[1502.05686].

[39] ATLAS collaboration, Search for massive supersymmetric particles in multi-jet final √ states produced in pp collisions at s = 13 TeV using the ATLAS detector at the LHC, ATLAS-CONF-2016-057 (2016). [40] M. J. Strassler and K. M. Zurek, Echoes of a hidden valley at hadron colliders, Phys.Lett. B651 (2007) 374–379,[hep-ph/0604261]. [41] M. J. Strassler, Possible effects of a hidden valley on supersymmetric phenomenology, hep-ph/0607160.

[42] CMS collaboration, S. Chatrchyan et al., Search for microscopic black holes in pp √ collisions at s = 8 TeV, JHEP 1307 (2013) 178,[1303.5338].

[43] UA2 collaboration, J. Alitti et al., A Study of multi-jet events at the CERN pp¯ collider and a search for double parton scattering, Phys. Lett. B268 (1991) 145–154.

[44] UA2 collaboration, J. Alitti et al., A Search for Squark and Gluino Production at the CERN pp¯ Collider, Phys. Lett. B235 (1990) 363–372.

[45] UA2 collaboration, M. Banner et al., Observation of Single Isolated Electrons of High Transverse Momentum in Events with Missing Transverse Energy at the CERN pp¯ Collider, Phys. Lett. 122B (1983) 476–485.

[46] UA2 collaboration, P. Bagnaia et al., Evidence for Z0 → e+e− at the CERN pp¯ Collider, Phys. Lett. B129 (1983) 130–140. [47] K. Jakobs, The Physics results of the UA2 experiment at the Cern pp¯ collider, Int. J. Mod. Phys. A9 (1994) 2903–2978. [48] T. Sjöstrand, S. Ask, J. R. Christiansen, R. Corke, N. Desai, P. Ilten et al., An Introduction to PYTHIA 8.2, Comput. Phys. Commun. 191 (2015) 159–177, [1410.3012]. [49] K. Meier. Personal communication. [50] R. D. Field and R. P. Feynman, A Parametrization of the Properties of Quark Jets, Nucl. Phys. B136 (1978) 1.

– 11 – [51] J. Alwall, R. Frederix, S. Frixione, V. Hirschi, F. Maltoni, O. Mattelaer et al., The automated computation of tree-level and next-to-leading order differential cross sections, and their matching to parton shower simulations, JHEP 07 (2014) 079,[1405.0301]. [52] W. Beenakker, R. Hopker, M. Spira and P. Zerwas, Squark and gluino production at hadron colliders, Nucl.Phys. B492 (1997) 51,[hep-ph/9610490]. [53] W. Beenakker, R. Hopker and M. Spira, PROSPINO: A Program for the production of supersymmetric particles in next-to-leading order QCD, hep-ph/9611232. [54] A. J. Larkoski, I. Moult and B. Nachman, Jet Substructure at the : A Review of Recent Advances in Theory and Machine Learning, 1709.04464. [55] CMS collaboration, A. M. Sirunyan et al., Search for low mass vector resonances √ decaying into quark-antiquark pairs in proton-proton collisions at s = 13 TeV, JHEP 01 (2018) 097,[1710.00159]. [56] T. Han, I. Lewis and Z. Liu, Colored Resonant Signals at the LHC: Largest Rate and Simplest Topology, JHEP 12 (2010) 085,[1010.4309]. [57] Y. Bai and J. Shelton, Composite Octet Searches with Jet Substructure, JHEP 07 (2012) 067,[1107.3563]. [58] M. Gresham, J. Shelton and K. M. Zurek, Open windows for a light axigluon explanation of the top forward-backward asymmetry, JHEP 03 (2013) 008,[1212.1718].

– 12 –