CMS logoCMS event Hgg
Compact Muon Solenoid
LHC, CERN

CMS-PAS-EXO-19-012
A search for dijet resonances in proton-proton collisions at $\sqrt{s}= $ 13 TeV with a new background prediction method
Abstract: Searches for narrow and broad resonances with a mass greater than 1.8 TeV decaying to a pair of jets are presented using proton-proton collision data at $\sqrt{s}= $ 13 TeV corresponding to an integrated luminosity of 137 fb$^{-1}$. The background arising from standard model processes is predicted using two complementary data-driven methods. The dijet invariant mass spectrum is well described by both methods, and no significant evidence for the production of new particles is observed. Model independent upper limits at 95% confidence level are reported on the production cross section of narrow resonances, and broad resonances with widths up to 30% of the resonance mass, extending previous searches. The limits are applied to various models of narrow resonances and exclude, at 95% confidence level: string resonances with masses below 7.9 TeV, scalar diquarks below 7.5 TeV, axigluons and colorons below 6.6 TeV, excited quarks below 6.3 TeV, color-octet scalars below 3.7 TeV, W' bosons below 3.6 TeV, Z' bosons with SM-like couplings below 2.9 TeV and between 3.1 TeV and 3.3 TeV, Randall-Sundrum Gravitons below 2.6 TeV TeV, and dark matter mediators below 2.8 TeV. Limits on dark matter mediators are presented as a function of resonance mass and width, and equivalently coupling, which exclude at 95% confidence level a dark matter mediator with mass less than 4.7 TeV for a width equal to 25% of the mass.
Figures & Tables Summary Additional Figures References CMS Publications
Figures

png pdf
Figure 1:
The pseudorapidity separation between the two wide jets for the signal and control regions. Data (black points) are compared to QCD predictions from the PYTHIA -8 MC with detector simulation (red histogram) normalized to the data. A signal from an RS Graviton decaying into a quark and an anti-quark is also shown (blue histogram) normalized to the same number of events as the data.

png pdf
Figure 2:
The dijet mass spectra of the data and PYTHIA simulation in the SR (black points and red histogram), $ {\mathrm {CR}_{\mathrm {middle}}} $ (triangles and blue line), and $ {\mathrm {CR}_{\mathrm {high}}} $ (crosses and magenta line).

png pdf
Figure 3:
Three dimensional display of the event with the second highest dijet invariant mass at 8 TeV. The display shows the energy deposited in the electromagnetic (red) and hadronic (blue) calorimeters and the reconstructed tracks of charged particles (green). The grouping of four observed jets into two wide jets (purple) is discussed in the text.

png pdf
Figure 4:
Left: $ {R_{\mathrm {aux.}}} $ auxiliary transfer factor for data (black points), PYTHIA (blue line), and POWHEG with electroweak corrections added on the top, along with their ratio fitted with the correction function on the bottom (magenta line with 95% CL error band). Right: $ {R} $ transfer factor for data (black points), PYTHIA (blue line), POWHEG with electroweak corrections added (red line) and corrected PYTHIA (magenta line) shown on the top. The ratio of data to PYTHIA (black points) fitted with the correction function (blue line with 95% CL error band), along with the ratio of POWHEG with electroweak corrections added to PYTHIA (red line), and the ratio of corrected PYTHIA to uncorrected PYTHIA (magenta line) are shown on the bottom. The corrected PYTHIA transfer factor (magenta line) using $ {\mathrm {CR}_{\mathrm {middle}}} $ agrees with the data one (blue line).

png pdf
Figure 4-a:
$ {R_{\mathrm {aux.}}} $ auxiliary transfer factor for data (black points), PYTHIA (blue line), and POWHEG with electroweak corrections added on the top, along with their ratio fitted with the correction function on the bottom (magenta line with 95% CL error band).

png pdf
Figure 4-b:
$ {R} $ transfer factor for data (black points), PYTHIA (blue line), POWHEG with electroweak corrections added (red line) and corrected PYTHIA (magenta line) shown on the top. The ratio of data to PYTHIA (black points) fitted with the correction function (blue line with 95% CL error band), along with the ratio of POWHEG with electroweak corrections added to PYTHIA (red line), and the ratio of corrected PYTHIA to uncorrected PYTHIA (magenta line) are shown on the bottom. The corrected PYTHIA transfer factor (magenta line) using $ {\mathrm {CR}_{\mathrm {middle}}} $ agrees with the data one (blue line).

png pdf
Figure 5:
Dijet mass spectrum in the signal region (points) compared to a fitted parameterization of the background (solid curve) and the one obtained from the control region (green squares). For the displayed signal a cross section at the 95% CL observed exclusion limit is being used. The lower panel shows the difference between the data and the fitted parametrization (red), and the data and the prediction obtained from the control region (green), divided by the statistical uncertainty of the data, which for the "ratio method" includes the one in $ {\mathrm {CR}_{\mathrm {high}}} $ as well. The ratio of the expected signal showed in the upper panel to the statistical uncertainty of the data is also shown for three different resonance masses and signals models.

png pdf
Figure 6:
The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for dijet resonances decaying to quark-quark (top left), quark-gluon (top right), gluon-gluon (bottom left), and for RS Gravitons (bottom right). The corresponding expected limits (dashed) and their variations at the 1 and 2 standard deviation levels (shaded bands) are also shown. Limits are compared to predicted cross sections for string resonances [1,2], excited quarks [4,5], axigluons [6], colorons [8], scalar diquarks [3], color-octet scalars [9], new gauge bosons W' and Z' with SM-like couplings [10], dark matter mediators for $ {m_{\text {DM}}} =$ 1 GeV [15,14], and RS Gravitons [11].

png pdf
Figure 6-a:
The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for dijet resonances decaying to quark-quark. The corresponding expected limits (dashed) and their variations at the 1 and 2 standard deviation levels (shaded bands) are also shown. Limits are compared to predicted cross sections for axigluons [6], scalar diquarks [3], new gauge bosons W' and Z' with SM-like couplings [10], and dark matter mediators for $ {m_{\text {DM}}} =$ 1 GeV [15,14].

png pdf
Figure 6-b:
The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for dijet resonances decaying to quark-gluon. The corresponding expected limits (dashed) and their variations at the 1 and 2 standard deviation levels (shaded bands) are also shown. Limits are compared to predicted cross sections for string resonances [1,2], and excited quarks [4,5].

png pdf
Figure 6-c:
The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for dijet resonances decaying to gluon-gluon. The corresponding expected limits (dashed) and their variations at the 1 and 2 standard deviation levels (shaded bands) are also shown. Limits are compared to predicted cross sections for color-octet scalars [9].

png pdf
Figure 6-d:
The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for RS Gravitons. The corresponding expected limits (dashed) and their variations at the 1 and 2 standard deviation levels (shaded bands) are also shown. Limits are compared to predicted cross sections for RS Gravitons [11].

png pdf
Figure 7:
The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for quark-quark, quark-gluon, and gluon-gluon type dijet resonances. Limits are compared to predicted cross sections for string resonances [1,2], excited quarks [4,5], axigluons [6], colorons [8], scalar diquarks [3], color-octet scalars [9], new gauge bosons W' and Z' with SM-like couplings [10], dark matter mediators for $ {m_{\text {DM}}} =$ 1 GeV [15,14], and RS Gravitons [11].

png pdf
Figure 8:
Local significance for a $\mathrm{q} \mathrm{q} $ resonance with the "ratio method" (blue) and the "fit method" (red).

png pdf
Figure 9:
The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for spin 2 resonances produced and decaying in the quark-quark (upper left) and gluon-gluon (upper right) channels, as well as for spin 1 resonances decaying in the quark-quark channel (bottom), shown for various values of intrinsic width and resonance mass.

png pdf
Figure 9-a:
The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for spin 2 resonances produced and decaying in the quark-quark channel, shown for various values of intrinsic width and resonance mass.

png pdf
Figure 9-b:
The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for

png pdf
Figure 9-c:
The observed 95% CL upper limits on the product of the cross section, branching fraction, and acceptance for spin 1 resonances decaying in the quark-quark channel, shown for various values of intrinsic width and resonance mass.

png pdf
Figure 10:
The 95% CL upper limits on the universal quark coupling $ {g_{\mathrm {q}}} $ as a function of resonance mass for a vector mediator of interactions between quarks and DM particles (left), and quarks only (right). The dotted horizontal lines on the right plot show the coupling strength for which the cross section for dijet production in this model is the same as for a DM mediator for $ {g_{\mathrm {q}}} =$ 0.25. The right vertical axis shows the natural width of the mediator divided by its mass. The expected limits (dashed) and their variation at the 1 and 2 standard deviation levels (shaded bands) are also shown.

png pdf
Figure 10-a:
The 95% CL upper limits on the universal quark coupling $ {g_{\mathrm {q}}} $ as a function of resonance mass for a vector mediator of interactions between quarks and DM particles. The right vertical axis shows the natural width of the mediator divided by its mass. The expected limits (dashed) and their variation at the 1 and 2 standard deviation levels (shaded bands) are also shown.

png pdf
Figure 10-b:
The 95% CL upper limits on the universal quark coupling $ {g_{\mathrm {q}}} $ as a function of resonance mass for a vector mediator of interactions between quarks only. The dotted horizontal lines show the coupling strength for which the cross section for dijet production in this model is the same as for a DM mediator for $ {g_{\mathrm {q}}} =$ 0.25. The right vertical axis shows the natural width of the mediator divided by its mass. The expected limits (dashed) and their variation at the 1 and 2 standard deviation levels (shaded bands) are also shown.
Tables

png pdf
Table 1:
Observed and expected mass limits at 95% CL from this analysis with 137 fb$^{-1}$ at $\sqrt {s}=$ 13 TeV compared to previously published limits on narrow resonances with 36 fb$^{-1}$ at $\sqrt {s}=$ 13 TeV [16]. The listed models are excluded between 1.8 TeV and the indicated mass limit by this analysis. The Z' model within the mass interval between 3.1 and 3.3 TeV is also excluded.
Summary
Searches for resonances decaying into a pair of jets have been performed using proton-proton collisions at $\sqrt{s}=$ 13 TeV corresponding to an integrated luminosity of 137 fb$^{-1}$.
The dijet mass spectra are observed to be smoothly falling distributions, and the QCD background is predicted using using two methods. The fit method uses an empirical functional form to fit the background in the signal region, dijet $|\Delta\eta| < $ 1.1, while the ratio method uses two control regions at higher values of $|\Delta\eta|$ to predict the background in the signal region. The ratio method is a new background prediction method which is independent of and complementary to the fit method. No evidence for resonant particle production is observed. Generic upper limits are presented on the product of the cross section, the branching fraction, and the acceptance for narrow and broad quark-quark, quark-gluon, and gluon-gluon resonances. The limits are applied to various models of new resonances and exclude at 95% confidence level, string resonances with masses below 7.9 TeV, scalar diquarks below 7.5 TeV, axigluons and colorons below 6.6 TeV, excited quarks below 6.3 TeV, color-octet scalars below 3.7 TeV, W' bosons below 3.6 TeV, Z' bosons with SM-like couplings below 2.9 TeV and between 3.1 and 3.3 TeV, Randall-Sundrum Gravitons below 2.6 TeV, and dark matter mediators below 2.8 TeV. This search extends previously reported limits on narrow resonances by between 200 and 800 GeV. Limits are also presented for resonances with intrinsic width as large as 30% of the resonance mass, and are used to improve and extend the exclusions of a dark matter mediator to larger values of the resonance mass and coupling to quarks. In the search for broad resonances, the ratio method provides significantly enhanced sensitivity compared to the fit method, resulting in the exclusion at 95% confidence level of a dark matter mediator with mass less than 4.7 TeV for a width equal to 25% of the mass, which corresponds to a coupling to quarks $g_q=$ 0.68.
Additional Figures

png pdf
Additional Figure 1:
The 95% CL upper limits on the universal quark coupling $g_{q}$ as a function of resonance mass for a vector mediator of interactions between quarks and dark matter, from the narrow resonance search, which are only valid for a width up to approximately 10% of the resonance mass. The right vertical axis shows the natural width of the mediator divided by its mass. The observed, expected limits (dashed) and their variation at the 1 and 2 standard deviation levels (shaded bands) are shown. The exclusions are computed for a spin-1 mediator and Dirac DM with a mass $m_{DM}=1$ GeV and a coupling $g_{DM}=$ 1.0.

png pdf
Additional Figure 2:
Ratio of expected limits between the Fit method and the Ratio method for wide Spin 1 resonances decaying to a pair of quarks, and for differentwidths.

png pdf
Additional Figure 3:
Ratio of expected limits between the Fit method and the Ratio method for wide Spin 2 resonances decaying to a pair of gluons, and for different widths.

png pdf
Additional Figure 4:
Three dimensional display of the event with the second highest dijet invariant mass at 8 TeV. The display shows the energy deposited in the electromagnetic (green) and hadronic (blue) calorimeters and the reconstructed tracks of charged particles (yellow).

png pdf
Additional Figure 5:
Three dimensional display of the event with the second highest dijet invariant mass at 8 TeV. The display shows the energy deposited in the electromagnetic (green) and hadronic (blue) calorimeters and the reconstructed tracks of charged particles (yellow).

png pdf
Additional Figure 6:
Three dimensional display of the event with the second highest dijet invariant mass at 8 TeV. The display shows the energy deposited in the electromagnetic (green) and hadronic (blue) calorimeters and the reconstructed tracks of charged particles (yellow).

png pdf
Additional Figure 7:
Three dimensional display of the event with the second highest dijet invariant mass at 8 TeV. The display shows the energy deposited in the electromagnetic (green) and hadronic (blue) calorimeters and the reconstructed tracks of charged particles (yellow).

png pdf
Additional Figure 8:
Three dimensional display of the event with the second highest dijet invariant mass at 8 TeV. The display shows the energy deposited in the electromagnetic (green) and hadronic (blue) calorimeters and the reconstructed tracks of charged particles (yellow).

png pdf
Additional Figure 9:
Three dimensional display of the event with the second highest dijet invariant mass at 8 TeV. The display shows the energy deposited in the electromagnetic (green) and hadronic (blue) calorimeters and the reconstructed tracks of charged particles (yellow).
References
1 L. A. Anchordoqui et al. Dijet signals for low mass strings at the LHC PRL 101 (2008) 241803 0808.0497
2 S. Cullen, M. Perelstein, and M. E. Peskin TeV strings and collider probes of large extra dimensions PRD 62 (2000) 055012 hep-ph/0001166
3 J. L. Hewett and T. G. Rizzo Low-energy phenomenology of superstring-inspired E(6) models PR 183 (1989) 193
4 U. Baur, I. Hinchliffe, and D. Zeppenfeld Excited quark production at hadron colliders Int. J. Mod. Phys. A 02 (1987) 1285
5 U. Baur, M. Spira, and P. M. Zerwas Excited quark and lepton production at hadron colliders PRD 42 (1990) 815
6 P. H. Frampton and S. L. Glashow Chiral color: An alternative to the standard model PLB 190 (1987) 157
7 R. S. Chivukula, E. H. Simmons, A. Farzinnia, and J. Ren Hadron collider production of massive color-octet vector bosons at next-to-leading order PRD 87 (2013) 094011 1303.1120
8 E. H. Simmons Coloron phenomenology PRD 55 (1997) 1678 hep-ph/9608269
9 T. Han, I. Lewis, and Z. Liu Colored resonant signals at the LHC: largest rate and simplest topology JHEP 12 (2010) 085 1010.4309
10 E. Eichten, I. Hinchliffe, K. D. Lane, and C. Quigg Supercollider physics Rev. Mod. Phys. 56 (1984) 579
11 L. Randall and R. Sundrum An alternative to compactification PRL 83 (1999) 4690 hep-th/9906064
12 M. Chala et al. Constraining dark sectors with monojets and dijets JHEP 07 (2015) 089 1503.05916
13 D. Abercrombie et al. Dark Matter Benchmark Models for Early LHC Run-2 Searches: Report of the ATLAS/CMS Dark Matter Forum 1507.00966
14 J. Abdallah et al. Simplified models for dark matter searches at the LHC Phys. Dark Univ. 9-10 (2015) 8 1506.03116
15 A. Boveia et al. Recommendations on presenting LHC searches for missing transverse energy signals using simplified $ s $-channel models of dark matter 1603.04156
16 CMS Collaboration Search for narrow and broad dijet resonances in proton-proton collisions at $ \sqrt{s}= $ 13 TeV and constraints on dark matter mediators and other new particles JHEP 08 (2018) 130 CMS-EXO-16-056
1806.00843
17 R. S. Chivukula, E. H. Simmons, and N. Vignaroli Distinguishing dijet resonances at the LHC PRD 91 (2015) 055019 1412.3094
18 ATLAS Collaboration Search for new phenomena in dijet events using 37 fb$ ^{-1} $ of $ pp $ collision data collected at $ \sqrt{s}= $ 13 TeV with the ATLAS detector PRD96 (2017) 052004 1703.09127
19 CMS Collaboration Search for dijet resonances in proton-proton collisions at $ \sqrt{s} = $ 13 TeV and constraints on dark matter and other models PLB 769 (2017) 520 CMS-EXO-16-032
1611.03568
20 CMS Collaboration Search for narrow resonances decaying to dijets in proton-proton collisions at $ \sqrt{s} = $ 13 TeV PRL 116 (2016) 071801 CMS-EXO-15-001
1512.01224
21 ATLAS Collaboration Search for new phenomena in dijet mass and angular distributions from pp collisions at $ \sqrt{s}= $ 13 TeV with the ATLAS detector PLB 754 (2016) 302 1512.01530
22 CMS Collaboration Search for narrow resonances in dijet final states at $ \sqrt{s}= $ 8 TeV with the novel CMS technique of data scouting PRL 117 (2016) 031802 CMS-EXO-14-005
1604.08907
23 CMS Collaboration Search for resonances and quantum black holes using dijet mass spectra in proton-proton collisions at $ \sqrt{s} = $ 8 TeV PRD 91 (2015) 052009 CMS-EXO-12-059
1501.04198
24 ATLAS Collaboration Search for new phenomena in the dijet mass distribution using pp collision data at $ \sqrt{s}= $ 8 TeV with the ATLAS detector PRD 91 (2015) 052007 1407.1376
25 CMS Collaboration Search for narrow resonances using the dijet mass spectrum in pp collisions at $ \sqrt{s} = $ 8 TeV PRD 87 (2013) 114015 CMS-EXO-12-016
1302.4794
26 CMS Collaboration Search for narrow resonances and quantum black holes in inclusive and b-tagged dijet mass spectra from pp collisions at $ \sqrt{s}= $ 7 TeV JHEP 01 (2013) 013 CMS-EXO-11-094
1210.2387
27 ATLAS Collaboration Search for new physics in the dijet mass distribution using 1 fb$ ^{-1} $ of $ pp $ collision data at $ \sqrt{s} = $ 7 TeV collected by the ATLAS detector PLB 708 (2012) 37 1108.6311
28 ATLAS Collaboration ATLAS search for new phenomena in dijet mass and angular distributions using pp collisions at $ \sqrt{s}= $ 7 TeV JHEP 01 (2013) 029 1210.1718
29 CMS Collaboration Search for resonances in the dijet mass spectrum from 7 TeV pp collisions at CMS PLB 704 (2011) 123 CMS-EXO-11-015
1107.4771
30 ATLAS Collaboration Search for new physics in dijet mass and angular distributions in $ pp $ collisions at $ \sqrt{s} = $ 7 TeV measured with the ATLAS detector New J. Phys. 13 (2011) 053044 1103.3864
31 CMS Collaboration Search for Dijet Resonances in 7 TeV pp Collisions at CMS PRL 105 (2010) 211801 CMS-EXO-10-010
1010.0203
32 ATLAS Collaboration Search for new particles in two-jet final states in 7 tev proton-proton collisions with the ATLAS detector at the LHC PRL 105 (2010) 161801 1008.2461
33 R. M. Harris and K. Kousouris Searches for dijet resonances at hadron colliders Int. J. Mod. Phys. A 26 (2011) 5005 1110.5302
34 CMS Collaboration The CMS experiment at the CERN LHC JINST 3 (2008) S08004 CMS-00-001
35 CMS Collaboration Particle-flow reconstruction and global event description with the cms detector JINST 12 (2017) P10003 CMS-PRF-14-001
1706.04965
36 M. Cacciari and G. P. Salam Dispelling the $ N^{3} $ myth for the $ k_\mathrm{t} $ jet-finder PLB 641 (2006) 57 hep-ph/0512210
37 M. Cacciari, G. P. Salam, and G. Soyez The anti-$ k_t $ jet clustering algorithm JHEP 04 (2008) 063 0802.1189
38 M. Cacciari, G. P. Salam, and G. Soyez FastJet user manual EPJC 72 (2012) 1896 1111.6097
39 M. Cacciari and G. P. Salam Pileup subtraction using jet areas PLB 659 (2008) 119 0707.1378
40 CMS Collaboration Jet energy scale and resolution in the CMS experiment in pp collisions at 8 TeV JINST 12 (2017), no. 02, P02014 CMS-JME-13-004
1607.03663
41 CMS Collaboration The CMS trigger system JINST 12 (2017), no. 01, P01020 CMS-TRG-12-001
1609.02366
42 CMS Collaboration Jet algorithms performance in 13 TeV data CMS-PAS-JME-16-003 CMS-PAS-JME-16-003
43 CMS Collaboration Searches for dijet resonances in pp collisions at $ \sqrt{s}= $ 13 TeV using data collected in 2016
44 T. Sjostrand, S. Mrenna, and P. Skands A brief introduction to PYTHIA 8.1 Comp. Phys. Comm. 178 (2008) 852 0710.3820
45 CMS Collaboration Event generator tunes obtained from underlying event and multiparton scattering measurements EPJC 76 (2016) 155 CMS-GEN-14-001
1512.00815
46 P. Skands, S. Carrazza, and J. Rojo Tuning PYTHIA 8.1: the Monash 2013 tune EPJC 74 (2014) 3024 1404.5630
47 GEANT4 Collaboration GEANT4 --- a simulation toolkit NIMA 506 (2003) 250
48 B. A. Dobrescu, R. M. Harris, and J. Isaacson Ultraheavy resonances at the LHC: beyond the QCD background 1810.09429
49 CDF Collaboration Search for new particles decaying into dijets in proton-antiproton collisions at $ \sqrt{s} = $ 1.96 ~TeV PRD 79 (2009) 112002 0812.4036
50 S. Frixione, P. Nason, and C. Oleari Matching NLO QCD computations with Parton Shower simulations: the POWHEG method JHEP 11 (2007) 070 0709.2092
51 S. Dittmaier, A. Huss, and C. Speckner Weak radiative corrections to dijet production at hadron colliders JHEP 11 (2012) 095 1210.0438
52 CMS Collaboration CMS luminosity measurements for the 2016 data taking period CMS-PAS-LUM-17-001 CMS-PAS-LUM-17-001
53 CMS Collaboration CMS luminosity measurements for the 2017 data taking period CMS-PAS-LUM-17-004 CMS-PAS-LUM-17-004
54 CMS Collaboration CMS luminosity measurements for the 2018 data taking period CMS-PAS-LUM-18-002 CMS-PAS-LUM-18-002
55 T. Junk Confidence level computation for combining searches with small statistics Nucl. Instr. Meth. A 434 (1999) 435 hep-ex/9902006
56 A. L. Read Presentation of search results: the $ CL_s $ technique JPG 28 (2002) 2693
57 LHC Higgs Combination Group Procedure for the LHC Higgs boson search combination in Summer 2011 CMS-NOTE-2011-005
58 G. Cowan, K. Cranmer, E. Gross, and O. Vitells Asymptotic formulae for likelihood-based tests of new physics EPJC 71 (2011) 1554 1007.1727
59 J. Pumplin et al. New generation of parton distributions with uncertainties from global QCD analysis JHEP 07 (2002) 012 hep-ph/0201195
60 V. D. Barger and R. J. N. Phillips Collider Physics 1987ISBN 9780201149456
61 B. A. Dobrescu and F. Yu Coupling-mass mapping of dijet peak searches PRD 88 (2013) 035021 1306.2629
62 J. Alwall 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
Compact Muon Solenoid
LHC, CERN