CMS-PAS-FTR-16-006

CMS-PAS-FTR-16-006
ECFA 2016: Prospects for selected standard model measurements with the CMS experiment at the High-Luminosity LHC
CMS Collaboration
May 2017

Abstract: The prospects for selected standard model measurements at the High-Luminosity LHC presented at ECFA 2016 workshop are summarized. The extrapolations assume proton-proton collision data at a centre-of-mass energy of 14 TeV corresponding to an integrated luminosity of up to 3 ab $^{-1}$ . The achievable precision for top quark mass measurements based on different analysis strategies is estimated. Searches for flavour-changing neutral currents in top quark decays are studied and expected limits are set, based on different scenarios for the extrapolation of systematic uncertainties to the High-Luminosity LHC run conditions. The feasibility of a dedicated track trigger for the $B_s \rightarrow \phi \phi$ decay studies is discussed.
Links: CDS record (PDF) ; inSPIRE record ; CADI line (restricted) ;

Figures & Tables	Summary	References	CMS Publications

Figures
png pdf	Figure 1: Total uncertainty on top quark mass ( $m_{\mathrm{ t } }$ ) obtained with different measurement methods and their projections to the HL-LHC for running conditions foreseen after the phase II upgrade. The projections for $\sqrt{s} =$ 14 TeV, with 0.3 ab $^{-1}$ or 3 ab $^{-1}$ of data, are based on $m_{\mathrm{ t } }$ measurements performed at the LHC Run-1, assuming that an upgraded detector will maintain the same physics performance despite a severe pileup.
png pdf	Figure 2: Transverse momentum (left) and pseudorapidity (right) of the photon candidates from t+ $\gamma$ production due to tu $\gamma$ FCNC interaction and various background processes with V = $\gamma$ , Z, W. The distributions are obtained using DELPHES simulation for the upgraded CMS detector at $\sqrt{s} =$ 14 TeV and on average 200 interactions per bunch crossing.
png pdf	Figure 2-a: Transverse momentum of the photon candidates from t+ $\gamma$ production due to tu $\gamma$ FCNC interaction and various background processes with V = $\gamma$ , Z, W. The distributions are obtained using DELPHES simulation for the upgraded CMS detector at $\sqrt{s} =$ 14 TeV and on average 200 interactions per bunch crossing.
png pdf	Figure 2-b: Pseudorapidity of the photon candidates from t+ $\gamma$ production due to tu $\gamma$ FCNC interaction and various background processes with V = $\gamma$ , Z, W. The distributions are obtained using DELPHES simulation for the upgraded CMS detector at $\sqrt{s} =$ 14 TeV and on average 200 interactions per bunch crossing.
png pdf	Figure 3: Upper limits at 95% CL on the branching fractions of t $\rightarrow$ u+ $\gamma$ (left) and t $\rightarrow$ c+ $\gamma$ (right) for an integrated luminosity up to 3 ab $^{-1}$ at $\sqrt{s} =$ 14 TeV with 200 interactions per bunch crossing on average. The black curve is the expected upper limit at 95% CL and green and yellow bands show the $\pm$ 1 and $\pm$ 2 standard deviations from the expected limits. The results are obtained for the scenario 2 that is described in the text.
png pdf	Figure 3-a: Upper limits at 95% CL on the branching fractions of t $\rightarrow$ u+ $\gamma$ for an integrated luminosity up to 3 ab $^{-1}$ at $\sqrt{s} =$ 14 TeV with 200 interactions per bunch crossing on average. The black curve is the expected upper limit at 95% CL and green and yellow bands show the $\pm$ 1 and $\pm$ 2 standard deviations from the expected limits. The results are obtained for the scenario 2 that is described in the text.
png pdf	Figure 3-b: Upper limits at 95% CL on the branching fractions of t $\rightarrow$ c+ $\gamma$ for an integrated luminosity up to 3 ab $^{-1}$ at $\sqrt{s} =$ 14 TeV with 200 interactions per bunch crossing on average. The black curve is the expected upper limit at 95% CL and green and yellow bands show the $\pm$ 1 and $\pm$ 2 standard deviations from the expected limits. The results are obtained for the scenario 2 that is described in the text.
png pdf	Figure 4: Expected upper limits at 95% CL on B(t $\rightarrow$ q+Z) and B(t $\rightarrow$ q+ $\gamma$ ) obtained from preliminary projections based on a DELPHES simulation. The horizontal dashed line corresponds to upper limit on B(t $\rightarrow$ q+Z) at 14 TeV with 3 ab $^{-1}$ [56]. The two vertical dashed and dashed-dotted lines show the results of this analysis. The two vertical solid lines are the observed CMS results on B(t $\rightarrow$ u+ $\gamma$ ) and B(t $\rightarrow$ c+ $\gamma$ ) at 95% CL [44] and the two solid horizontal lines are the current observed 95% CL upper limits on B(t $\rightarrow$ u+Z) and B(t $\rightarrow$ c+Z) from 8 TeV CMS data [57].
png pdf	Figure 5: Feynman graph of the dominant amplitude contributing to the decay ${ {B^0_s} \to \phi \phi }$ .
png pdf	Figure 6: Invariant mass distribution of all track pairs with opposite charges, ${ {\| d_{z} \| } }<$ 1 cm, ${ {d_{xy}} }<$ 1 cm, track ${p_{\mathrm{T}}} >$ 2 GeV, and assuming that the tracks are arising from kaons. The event sample does not have a preliminary selection on the ${B^0_s}$ mass window. The distributions are normalized to unit area. The blue solid histogram corresponds to the signal events reconstructed with offline tracks, the red dashed one with tracks from L1 trigger system and the green filled area represents the background events.
png pdf	Figure 7: ${ {\Delta R} }$ ( $\phi$ -pair) distribution for all $\phi$ -pairs with 0.99 $< M_{K^+K^-} <$ 1.04 GeV, ${ {\| d_{z} \| } }<$ 1 cm, ${ {d_{xy}} }<$ 1 cm. The event sample does not have a preliminary selection on the ${B^0_s}$ mass window. The distributions are normalized to unit area. The blue solid histogram corresponds to the signal events reconstructed with offline tracks, the red dashed one with tracks from L1 trigger system and the green filled area represents the background events.
png pdf	Figure 8: Invariant mass distribution of all the $\phi$ -pairs with ${ {\| d_{z} \| } } (\phi \text {-pair}) <$ 1 cm, ${ {d_{xy}} } (\phi \text {-pair})<$ 1 cm, 0.2 $< { {\Delta R} } (\phi \text {-pair}) <$ 1, ${ {\Delta R} }(K^{+}, K^{-}) <$ 0.12. The distributions are normalized to unit area. The blue solid histogram corresponds to the signal events reconstructed with offline tracks, the red dashed one with tracks from L1 trigger system and the green filled area represents the background events.
png pdf	Figure 9: Efficiency and rate for different selection baselines and for different pileup scenarios. Uncertainties are statistical only.

Tables
png pdf	Table 1: Summary of the systematic uncertainties on $m_{\mathrm{ t } }$ for the reference measurement in lepton+jets channel. Experimental uncertainties are separated from theoretical ones.
png pdf	Table 2: Summary of the systematic uncertainties on $m_{\mathrm{ t } }$ for the measurements in the single-top quark $t$ -channel. Experimental uncertainties are separated from theoretical ones.
png pdf	Table 3: Summary of the systematic uncertainties on $m_{\mathrm{ t } }$ for the measurement from $m_{\text {sv}\ell }$ . Experimental uncertainties are separated from theoretical ones.
png pdf	Table 4: Summary of the systematic uncertainties on $m_{\mathrm{ t } }$ for the measurement from $m_{\mathrm{J}/\psi +\ell }$ . Experimental uncertainties are separated from theoretical ones.
png pdf	Table 5: Summary of the systematic uncertainties on $m_{\mathrm{ t } }$ for the measurement from $\sigma _{{\mathrm{ t } {}\mathrm{ \bar{t} } } }$ . Experimental uncertainties are separated from theoretical ones.
png pdf	Table 6: Upper limits at 95% CL for B $($ t $\rightarrow$ u+ $\gamma )$ and B $($ t $\rightarrow$ c+ $\gamma )$ , obtained with the 8 TeV data and the projections for 14 TeV with an integrated luminosity of 3 ab $^{-1}$ using CMS DELPHES simulation for two scenarios presented in the text.
png pdf	Table 7: Baseline event selection conditions. The variable ${ {d_{z}} }$ represents distance between a pair of tracks or trajectories of a pair of reconstructed particles along the beam axis ( $z$ ) while ${ {d_{xy}} }$ represents that in the plane perpendicular to the beam axis ( $xy$ ).
png pdf	Table 8: Efficiency and rate for loose, medium and tight baselines respectively. Pileup dependence of event rate is also presented for $< \text{PU } > =$ 70, 140 and 200. Uncertainties are statistical only.

Summary

The three physics proposals for the upgrade studies for the HL-LHC CMS detector discussed in this note were prepared for and presented at the ECFA 2016 workshop.

It is demonstrated that with 3 ab

$^{-1}$ of data the top quark mass analyses will be limited by systematic uncertainties, and especially by theoretical modeling uncertainties. The reference method, which is the most precise one, is expected to yield an ultimate relative precision below 0.1%. The other techniques, with alternative systematic sensitivity, are expected to reach a precision good enough to carry weight in a combination with the reference method. This would make it possible to further reduce the systematic uncertainties, which are related mostly to the JES and hard process modeling.

According to the projections for a search for the FCNC process in the top quark production associated with a photon at a luminosity of 3 ab

$^{-1}$ upper limits at 95% CL on the branching fractions B(t

$\rightarrow$ u+

$\gamma$ )

$<$ 0.0027% and B( t

$\rightarrow$ c+

$\gamma$ )

$<$ 0.020% are expected.

The

$\mathrm{B}_s \to 4 \mathrm{K}$ channel is used to investigate capabilities of the HL-LHC CMS detector to trigger events in the low-

$p_{\mathrm{T}}$ region for fully-hadronic final states. The study uses the track trigger to estimate the efficiency for selecting the signal events and the trigger rate of the background events. It is demonstrated that with the track trigger sufficient efficiency can be achieved, while the trigger rate requires further improvement e.g. by including displaced vertex finding tool for low

$p_{\mathrm{T}}$ tracks and a mitigation of pileup effects.

References
1	G. Apollinari et al.	High-Luminosity Large Hadron Collider (HL-LHC) : Preliminary Design Report
2	D. Contardo et al.	Technical Proposal for the Phase-II Upgrade of the CMS Detector	Technical Report CERN-LHCC-2015-010. LHCC-P-008. CMS-TDR-15-02, Geneva, Jun
3	CMS Collaboration	Technical Proposal for the Phase-II Upgrade of the CMS Detector	CMS-PAS-TDR-15-002	CMS-PAS-TDR-15-002
4	ECFA	3rd ECFA High Luminosity LHC Workshop	2016
5	G. Degrassi et al.	Higgs mass and vacuum stability in the Standard Model at NNLO	JHEP 08 (2012) 1	1205.6497
6	F. Bezrukov, M. Y. Kalmykov, B. A. Kniehl, and M. Shaposhnikov	Higgs boson mass and new physics	JHEP 10 (2012) 140	1205.2893
7	F. Bezrukov and M. Shaposhnikov	The Standard Model Higgs boson as the inflaton	PLB 659 (2008) 703	0710.3755
8	A. D. Simone, M. P. Herzberg, and F. Wilczek	Running inflation in the Standard Model	PLB 678 (2009) 1	0812.4946
9	CMS Collaboration	Projected improvement of the accuracy of top-quark mass measurements at the upgraded LHC	CDS
10	CMS Collaboration	Measurement of the top quark mass using proton-proton data at $\sqrt{s} =$ 7 and 8 TeV	PRD 93 (2016) 072004	CMS-TOP-14-022 1509.04044
11	A. H. Hoang	The Top Mass: Interpretation and Theoretical Uncertainties	in Proceedings, 7th International Workshop on Top Quark Physics (TOP2014): Cannes, 2014	1412.3649
12	CMS Collaboration	Measurement of the $\mathrm{ t \bar{t} }$ production cross section in the e $\mu$ channel in proton-proton collisions at $\sqrt{s} =$ 7 and 8 TeV	JHEP 08 (2016) 029	CMS-TOP-13-004 1603.02303
13	ATLAS Collaboration	Measurement of the Top Quark Mass from $\sqrt{s} =$ 7 TeV ATLAS Data using a 3-dimensional Template Fit	ATLAS Conference Note ATLAS-CONF-2013-046, CERN
14	CMS Collaboration	Calculation of residual energy correction for b jets using Z+b events in 8 TeV pp Collisions	CMS-PAS-JME-13-001	CMS-PAS-JME-13-001
15	CMS Collaboration	Measurement of the top quark mass in the dileptonic $\text{t}\bar{\text{t}}$ decay channel using the $\text{M}_{\text{b}\ell}$ , $\text{M}_{\text{T}2}$ , and MAOS $\text{M}_{\text{b}\ell\nu}$ observables	CMS-PAS-TOP-15-008	CMS-PAS-TOP-15-008
16	CMS Collaboration Collaboration	Measurement of the top quark mass using proton-proton data at $\sqrt{s} =$ 7 and 8 TeV	PRD 93 (2015) 072004
17	CMS Collaboration	Pileup Removal Algorithms	CMS-PAS-JME-14-001	CMS-PAS-JME-14-001
18	M. Czakon, P. Fiedler, and A. Mitov	Total Top-Quark Pair-Production Cross Section at Hadron Colliders Through $\mathcal{O}(\alpha^4_S)$	PRL 110 (2013) 252004	1303.6254
19	CMS Collaboration	Study of the underlying event, b-quark fragmentation and hadronization properties in $\mathrm{ t \bar{t} }$ events	CDS
20	CMS Collaboration	Measurement of the top quark mass using charged particles in pp collisions at $\sqrt{s} =$ 8 TeV	PRD 93 (2016) 2006	CMS-TOP-12-030 1603.06536
21	CMS Collaboration	Underlying event measurement with $\mathrm{ t \bar{t} }+X$ events with p-p collision data at $\sqrt{s} =$ 13 TeV	CMS-PAS-TOP-15-017	CMS-PAS-TOP-15-017
22	CMS Collaboration	Comparisons of Theory Predictions for the $\mathrm{ t \bar{t} }$ Process with Data from pp Collisions at $\sqrt{s} =$ 8 TeV	CMS-PAS-TOP-15-011	CMS-PAS-TOP-15-011
23	CMS Collaboration	Investigations of the impact of the parton shower tuning in Pythia 8 in the modelling of $\mathrm{ t \bar{t} }$ at $\sqrt{s}=$ 8 and 13 TeV	CMS-PAS-TOP-16-021	CMS-PAS-TOP-16-021
24	CMS Collaboration	Measurement of the differential cross section for top quark pair production in pp collisions at $\sqrt{s} =$ 8 TeV	EPJC 75 (2015) 542	CMS-TOP-12-028 1505.04480
25	M. Czakon, D. Heymes, and A. Mitov	High-precision differential predictions for top-quark pairs at the LHC	PRL 116 (2016) 2003	1511.00549
26	CMS Collaboration Collaboration	Measurement of the top quark mass with muon+jets final states in pp collisions at $\sqrt{s}=$ 13 TeV	Technical Report CMS-PAS-TOP-16-022, CERN, Geneva
27	CMS Collaboration	Measurement of the top quark mass from single-top production events	CMS-PAS-TOP-15-001	CMS-PAS-TOP-15-001
28	CMS Collaboration	Measurement of the top quark mass using charged particles in pp collisions at $\sqrt {s} =$ 8 TeV	PRD 93 (2016) 2006	CMS-TOP-12-030 1603.06536
29	CMS Collaboration	Measurement of the mass of the top quark in decays with a $\mathrm{J}/\psi$ meson in pp collisions at $8 TeV$	JHEP 12 (2016) 123	CMS-TOP-15-014 1608.03560
30	M. Butenschoen et al.	Top Quark Mass Calibration for Monte Carlo Event Generators	PRL 117 (2016), no. 23, 232001	1608.01318
31	CMS Collaboration	Measurement of the $\mathrm{ t \bar{t} }$ production cross section using events with one lepton and at least one jet in pp collisions at $\sqrt{s} =$ 13 TeV		CMS-TOP-16-006 1701.06228
32	J. Kieseler, K. Lipka, and S.-O. Moch	Calibration of the Top-Quark Monte Carlo Mass	PRL 116 (2016) 2001	1511.00841
33	J. Wenninger	Energy Calibration of the LHC Beams at 4 TeV	Technical Report CERN-ATS-2013-040, CERN, Geneva
34	J. Wenninger and E. Todesco	Large Hadron Collider momentum calibration and accuracy	Technical Report CERN-ACC-2017-0007, CERN, Geneva, Feb
35	S. L. Glashow, J. Iliopoulos, and L. Maiani	Weak Interactions with Lepton Hadron Symmetry	PRD 2 (1970) 1285
36	J. A. Aguilar-Saavedra and B. M. Nobre	Rare top decays $t \rightarrow c \gamma, t \rightarrow cg$ and CKM unitarity	PLB 553 (2003) 251	hep-ph/0210360
37	G. Couture, M. Frank, and H. Konig	Supersymmetric QCD flavor changing top quark decay	PRD 56 (1997) 4213	hep-ph/9704305
38	G. R. Lu, F. R. Yin, X. L. Wang, and L. D. Wan	Rare top quark decays t $\to$ cV in the topcolor assisted technicolor model	PRD 68 (2003) 015002	hep-ph/0303122
39	L3 Collaboration	Search for Single Top Production at LEP	PLB 549 (2002) 290	hep-ex/0210041
40	H1 Collaboration	Search for Single Top Quark Production at HERA	PLB 678 (2009) 450	0904.3876
41	ZEUS Collaboration	Search for single-top production in ep collisions at HERA	PLB 708 (2012) 27	1111.3901
42	CDF Collaboration	Search for Flavor-Changing Neutral Current Decays of the Top Quark in $\rm p\bar{p}$ Collisions at $\sqrt{s} =$ 1.8 TeV	PRL 80 (1998) 2525
43	S. Khatibi and M. Mohammadi Najafabadi	Constraints on top quark flavor changing neutral currents using diphoton events at the LHC	Nucl. Phys. B 909 (2016) 607	1511.00220
44	CMS Collaboration	Search for anomalous single top quark production in association with a photon in pp collisions at $\sqrt{s}=$ 8 TeV	JHEP 04 (2016) 035	CMS-TOP-14-003 1511.03951
45	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
46	T. Sj\"ostrand, S. Mrenna, and P. Skands	PYTHIA 6.4 physics and manual	JHEP 05 (2006) 026	hep-ph/0603175
47	R. D. Ball et al.	Parton distributions with LHC data	Nucl. Phys. B867 (2013) 244--289	1207.1303
48	A. Avetisyan et al.	Methods and Results for Standard Model Event Generation at $\sqrt{s} =$ 14 TeV, 33 TeV and 100 TeV Proton Colliders (A Snowmass Whitepaper)		1308.1636
49	J. Anderson et al.	Snowmass Energy Frontier Simulations		1309.1057
50	J. de Favereau et al.	DELPHES 3, A modular framework for fast simulation of a generic collider experiment	JHEP 1402 (2013) 057	1307.6346
51	CMS Collaboration	Performance of the b-jet identification in CMS	CDS
52	CMS Collaboration	Measurement of the $t$ -channel single top quark production cross section in $pp$ collisions at $\sqrt{s} =$ 7 TeV	PRL 107 (2011) 091802
53	A. L. Read	Presentation of search results: the CL $_\mathrm{s}$ technique	JPG 28 (2002) 2693
54	T. Junk	Confidence level computation for combining searches with small statistics	NIMA 434 (1999) 435
55	L. Moneta et al.	The RooStats Project	in 13 $^th$ International Workshop on Advanced Computing and Analysis Techniques in Physics Research (ACAT2010) SISSA, 2010	1009.1003
56	CMS Collaboration	Projections for Top FCNC Searches in 3000 fb $^{-1}$ at the LHC	CDS
57	CMS Collaboration	Search for associated production of a Z boson with a single top quark and for tZ flavour-changing interactions in pp collisions at $\sqrt{s}=$ 8 TeV		CMS-TOP-12-039 1702.01404
58	M. Kobayashi and T. Maskawa	CP Violation in the Renormalizable Theory of Weak Interaction	Prog. Theor. Phys. 49 (1973) 652--657
59	N. Cabibbo	Unitary Symmetry and Leptonic Decays	PRL 10 (1963) 531--533, .[,648(1963)]
60	J. Charles et al.	Predictions of selected flavour observables within the Standard Model	PRD84 (2011) 033005	1106.4041
61	D0 Collaboration	Measurement of the CP-violating phase $\phi_s^{J/\psi \phi}$ using the flavor-tagged decay $B_s^0 \rightarrow J/\psi \phi$ in 8 fb $^{-1}$ of $p \bar p$ collisions	PRD85 (2012) 032006	1109.3166
62	CDF Collaboration	Measurement of the CP-Violating Phase $\beta_s^{J/\psi\phi}$ in $B^0_s \to J/\psi \phi$ Decays with the CDF II Detector	PRD85 (2012) 072002	1112.1726
63	LHCb Collaboration	Measurement of the CP-violating phase $\phi_s$ in the decay $B^0_s\to J/\psi \phi$	PRL 108 (2012) 101803	1112.3183
64	ATLAS Collaboration	Time-dependent angular analysis of the decay $B_{s}^{0} \to J/{\psi\phi}$ and extraction of $\Delta\Gamma_{s}$ and the CP-violating weak phase $\phi_s$ by ATLAS	JHEP 12 (2012) 072	1208.0572
65	CMS Collaboration	Measurement of the CP-violating weak phase $\phi_s$ and the decay width difference $\Delta \Gamma_s$ using the $B_s^0 \to \jpsi\phi(1020)$ decay channel in pp collisions at $\sqrt{s}=$ 8 TeV	PLB757 (2016) 97--120	CMS-BPH-13-012 1507.07527
66	LHCb Collaboration	Measurement of CP violation and the $B^0_s$ meson decay width difference with $B^0_s \to\jpsi K^+K^-$ and $B^0_s\to\jpsi\pip\pim$ decays	PRD87 (2013) 2010	1304.2600
67	LHCb Collaboration	Measurement of the CP-violating phase $\phi_s$ in $\overline{B}^0_s\rightarrow J/\psi \pi^+\pi^-$ decays	PLB 736 (2014) 186--195	1405.4140
68	LHCb Collaboration	Precision measurement of $CP$ violation in $B_s^0 \to J/\psi K^+K^-$ decays	PRL 114 (2015) 1801	1411.3104
69	CDF Collaboration	First evidence for $B_s^0 \to \phi \phi$ decay and measurements of branching ratio and $A_{CP}$ for $B^+ \to \phi K^+$	PRL 95 (2005) 031801	hep-ex/0502044
70	CDF Collaboration	Measurement of Polarization and Search for CP-Violation in $B_s^0 \to \phi\phi$ Decays	PRL 107 (2011) 261802	1107.4999
71	LHCb Collaboration	First measurement of the CP-violating phase in $B_s^0 \to \phi \phi$ decays	PRL 110 (2013) 1802	1303.7125
72	LHCb Collaboration	Measurement of the polarization amplitudes and triple product asymmetries in the $B_s^0 \to \phi\phi$ decay	PLB713 (2012) 369--377	1204.2813
73	LHCb Collaboration	Measurement of CP violation in $B_s^0 \to \phi \phi$ decays	PRD90 (2014) 2011	1407.2222
74	LHCb Collaboration	Measurement of the $B_s^0 \to \phi \phi$ branching fraction and search for the decay $B^0 \to \phi \phi$	JHEP 10 (2015) 053	1508.00788
75	M. Bartsch, G. Buchalla, and C. Kraus	$B\to V(L) V(L)$ Decays at Next-to-Leading Order in QCD		0810.0249
76	M. Beneke, J. Rohrer, and D. Yang	Branching fractions, polarisation and asymmetries of $B\to VV$ decays	Nucl. Phys. B 774 (2007) 64--101	hep-ph/0612290
77	H.-Y. Cheng and C.-K. Chua	QCD Factorization for Charmless Hadronic $B_s$ Decays Revisited	PRD 80 (2009) 114026	0910.5237
78	CMS Collaboration	Measurement of the Strange $B$ Meson Production Cross Section with $J/\psi\phi$ Decays in $pp$ Collisions at $\sqrt{s}=$ 7 TeV	PRD84 (2011) 052008	CMS-BPH-10-013 1106.4048
79	CMS Collaboration	Measurement of the $B_s \to \mu^+ \mu^-$ branching fraction and search for $B^0 \to \mu^+ \mu^-$ with the CMS Experiment	PRL 111 (2013) 101804	CMS-BPH-13-004 1307.5025
80	D. J. Lange	The EvtGen particle decay simulation package	NIMA462 (2001) 152--155