EUROPEAN ORGANIZATION FOR NUCLEAR RESEARCH (CERN)

​​​ CERN-EP-2024-158 LHCb-PAPER-2024-005 July 16, 2024

Search for the rare ${{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}{\mu^{+}\mu^{-}}$ decay

LHCb collaboration

A search for the nonresonant ${{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}{\mu^{+}\mu^{-}}$ decay is performed using proton-proton collision data recorded at a centre-of-mass energy of $13\,\text{\,Te\kern-1.00006ptV}$ by the LHCb experiment, corresponding to an integrated luminosity of $5.4\text{\,fb}^{-1}$. No evidence for the decay is found in the dimuon invariant-mass regions where the expected contributions of resonances is subdominant. The upper limit on the branching fraction of the ${{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}{\mu^{+}\mu^{-}}$ decay is determined to be $2.9~{}(3.2)\times 10^{-8}~{}{\rm at}~{}90\%~{}(95\%)$ confidence level. The branching fractions in the dimuon invariant-mass regions dominated by the $\eta$, $\rho$ and $\omega$ resonances are also determined.

Submitted to Phys. Rev. D

© 2024 CERN for the benefit of the LHCb collaboration. CC BY 4.0 licence.

 

1 Introduction

The flavour-changing neutral-current decay ${{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}{\mu^{+}}{\mu^{-}}$ is heavily suppressed in the Standard Model (SM) as it is forbidden at tree level and is affected by the Glashow–Iliopoulos–Maiani (GIM) mechanism [1]. The phenomenology of the decay is discussed in Ref. [2] and the SM predictions for $q^{2}$ differential branching fractions, where $q^{2}$ is the dimuon invariant-mass squared, are calculated in Ref. [3]. The branching fraction ($\mathcal{B}$) associated with the short-distance contributions to the $c\rightarrow u\mu^{+}\mu^{-}$ transition is expected to be below $10^{-8}$ in the SM, with an enhancement up to $10^{-6}$ due to long-distance processes involving intermediate vector-meson resonances that decay to a lepton pair [4, 5]. The estimation of resonant SM contributions relies on the measured values for the corresponding branching fractions and the isospin relation between $\rho$ and $\omega$ states, to overcome the lack of experimental information on $\mathcal{B}$(${{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}\rho$). Strong phases between the resonances, and thus the effect of interference, are a priori unknown.

Interesting deviations from the SM predictions have been observed in the equivalent type of transition for the ${b}$ quark, ${b}\!\rightarrow{s}\ell\ell$, in particular in angular observables [6]. These anomalies motivate the extension of this physics program into the charm sector, which is far less explored both experimentally and theoretically. The first results, such as an angular analysis of charm mesons [7], are already available. A discussion of the discovery potential of effects beyond the SM in rare charm-baryon decays can be found in Ref. [8].

The LHCb collaboration previously searched for the ${{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}{\mu^{+}}{\mu^{-}}$ decay using the proton-proton ($pp$) data collected in 2011–2012 [9] and placed an upper limit of $7.7\times 10^{-8}$ at 90% confidence level (CL), thereby improving the result from the BaBar experiment [10] by two orders of magnitude. The analysis of data collected in 2016–2018 presented in this paper follows the methodology of the previous analysis, and relies on improved techniques with the use of the CL_s method [11] implementation based on unbinned likelihoods, the incorporation of background due to particle misidentification into the data model, a different definition of the short-distance signal and a more detailed treatment of the entire range of dimuon invariant mass where resonances are expected to dominate. Accordingly, the new results for the short distance contribution should not be combined with the results based on the 2011–2012 data.

The measurement is carried out in different regions of dimuon invariant mass $m_{\mu\mu}$ and the ${{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}\phi(\rightarrow{\mu^{+}\mu^{-}})$ decay is used for normalisation since it has the same final state as the signal and is relatively abundant.

2 Detector and simulation

The LHCb detector [12, 13] is a single-arm forward spectrometer covering the pseudorapidity range $2<\eta<5$, designed for the study of particles containing $b$ or $c$ quarks. The detector includes a high-precision tracking system consisting of a silicon-strip vertex detector surrounding the $pp$ interaction region [14], a large-area silicon-strip detector located upstream of a dipole magnet with a bending power of about $4{\mathrm{\,T\,m}}$, and three stations of silicon-strip detectors and straw drift tubes [15] placed downstream of the magnet. The tracking system provides a measurement of the momentum, $p$, of charged particles with a relative uncertainty that varies from 0.5% at low momentum to 1.0% at 200$\text{\,Ge\kern-1.00006ptV\!/}c$. The minimum distance of a track to a primary $pp$ collision vertex (PV), the impact parameter (IP), is measured with a resolution of $(15+29/p_{\mathrm{T}})\,\upmu\text{m}$, where $p_{\mathrm{T}}$ is the component of the momentum transverse to the beam, in $\text{\,Ge\kern-1.00006ptV\!/}c$. Different types of charged hadrons are distinguished using information from two ring-imaging Cherenkov detectors [16]. Photons, electrons and hadrons are identified by a calorimeter system consisting of scintillating-pad and preshower detectors, an electromagnetic and a hadronic calorimeter. Muons are identified by a system composed of alternating layers of iron and multiwire proportional chambers [17]. The online event selection is performed by a trigger system [18], which consists of a hardware stage, based on information from the calorimeter and muon systems, followed by a software stage, which applies a full event reconstruction.

Simulated samples of ${\mathchar 28931\relax}^{+}_{c}$ baryons are produced to model the effects of the detector acceptance and the imposed selection requirements. The $pp$ collisions are generated using Pythia [19] with a specific LHCb configuration [21]. Decays of unstable particles are described by EvtGen [22], in which final state radiation is generated using Photos [23]. In the case of ${{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}{\mu^{+}\mu^{-}}$ decays, a model with final state particles equally distributed in the phase space is used. The interaction of the generated particles with the detector, and its response, are implemented using the Geant4 toolkit [24] as described in Ref. [26]. The underlying $pp$ interaction and ${\mathchar 28931\relax}^{+}_{c}$ baryon are reused multiple times to speed up the simulation process, with independently generated ${\mathchar 28931\relax}^{+}_{c}$ decay products [27].

3 Candidate selection

The selection criteria largely follow those described in Ref. [9]. The ${\mathchar 28931\relax}^{+}_{c}$ candidates are reconstructed by combining a pair of oppositely charged tracks identified as muons with an additional track identified as a proton. The three tracks are required to be compatible with originating from a common vertex and have significant $\chi^{2}_{\text{IP}}$ with respect to all PVs. Here, $\chi^{2}_{\text{IP}}$ is defined as the difference in the vertex-fit $\chi^{2}$ of a given PV, reconstructed with and without the track under consideration. The ${\mathchar 28931\relax}^{+}_{c}$ candidates must have a decay vertex displaced from any PV and be compatible with originating from one of the PVs.

The background formed by random combinations of tracks, referred to as combinatorial background, is suppressed using a multivariate classifier. A boosted decision tree (BDT) algorithm [28], as implemented in the TMVA toolkit [29, 30], is trained using supervised learning with ten-fold cross-validation [31] to achieve an unbiased classifier response. The BDT discrimination is done in two steps. A first BDT classifier employs topological and kinematic variables of the ${\mathchar 28931\relax}^{+}_{c}$ candidates. A relatively loose requirement on the response of that BDT is applied as the last preselection step to enhance the purity of the normalisation data sample. After preselection, simulated samples of ${{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}\phi(\rightarrow{\mu^{+}\mu^{-}})$ decays are compared with data, and weights are determined to improve the agreement between data and simulation for the ${\mathchar 28931\relax}^{+}_{c}$ transverse momentum and the number of tracks in the event. A second BDT classifier is trained using a set of variables that includes those of the first classifier and, in addition, kinematic, vertex and isolation variables related to the decay products. For the two classifiers, simulated events containing ${{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}{\mu^{+}\mu^{-}}$ decays are used as a proxy for the signal, while the background is modelled using sideband events with a $p{\mu^{+}}{\mu^{-}}$ invariant-mass within [2140, 2251] and [2321, 2430] $\text{\,Me\kern-1.00006ptV\!/}c^{2}$. The optimisation of final requirements is performed in three dimensions: the response value of the second BDT classifier and two variables characterising the proton and muon particle identification (PID). In this optimisation, the simulation is corrected for imperfections in modeling of the PID performance [32]. The final requirements are determined using pseudoexperiments by finding the best expected upper limit (using the CL_s method) on the ratio of branching fractions of the signal relative to the normalisation channel. At this stage, only two background sources are expected to significantly contribute: combinatorial background and ${{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}\pi^{+}\pi^{-}$ decays where both pions are misidentified as muons. In each pseudoexperiment, pseudodata for the signal and these two background components are drawn from the three-dimensional probability density functions determined after the preselection stage.

Six dimuon invariant-mass ranges are considered as presented in Table 1. The nonresonant signal decay is searched for in two $m_{\mu\mu}$ regions, $low\text{-}m$ and $high\text{-}m$, where the expected contributions of resonances are subdominant. In the following, the combination of these two ranges is referred to as the signal region. The non signal range is divided into regions where either the $\eta$, $\rho$, $\omega$ or $\phi$ resonance dominates. The branching fraction results quoted for a given resonant region are labeled by the resonant decay expected to dominate in that region. However, due to the relatively low sample size, no attempt to separate the resonant contributions within a given dimuon invariant-mass range is made. Therefore it has to be emphasised that the results must be interpreted as corresponding to that particular $m_{\mu\mu}$ region rather than to a specific resonant decay.

Table 1: Definitions of the dimuon invariant-mass ranges. The lower (upper) bound of the $low\text{-}m$ ($high\text{-}m$) region corresponds to the phase-space limits.

Region	Range [$\text{\,Me\kern-1.00006ptV\!/}c^{2}$ ]
$low\text{-}m$ region	$~{}~{}~{}211.32<m_{{\mu^{+}\mu^{-}}}<~{}507.86$
$high\text{-}m$ region	$~{}~{}1059.46<m_{{\mu^{+}\mu^{-}}}<1348.13$
$\eta$ region	$~{}~{}~{}507.86<m_{{\mu^{+}\mu^{-}}}<~{}587.86$
$\omega$ region	$~{}~{}~{}742.65<m_{{\mu^{+}\mu^{-}}}<~{}822.65$
$\rho$ region	$~{}~{}~{}587.86<m_{{\mu^{+}\mu^{-}}}<~{}742.65$
	or $822.65~{}<m_{{\mu^{+}\mu^{-}}}<~{}965.20~{}$
$\phi$ region	$~{}~{}~{}~{}979.46~{}<m_{{\mu^{+}\mu^{-}}}<1059.46$

4 Normalisation and invariant-mass fit

The ${{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}{\mu^{+}\mu^{-}}$ branching fraction is determined from the measured signal yield $N_{\rm sig}$ according to

$${{\cal B}({{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}{\mu^{+}\mu^{-}})=\alpha\times N_{\rm sig}},$$ (1)

where the normalisation factor $\alpha$ is given by

$${\alpha=\frac{\rm{\epsilon\mathstrut_{norm}}}{\rm{\epsilon\mathstrut_{\rm sig}}\times\it{N}_{\rm norm}}}\times{{\cal B}({{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}\phi)\times{\cal B}(\phi\!\rightarrow{\mu^{+}\mu^{-}})}.$$ (2)

Here, $N_{\rm norm}$ is the yield of ${{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}\phi(\rightarrow{\mu^{+}\mu^{-}})$ decays, while $\rm\epsilon\mathstrut_{sig}$ and $\rm\epsilon\mathstrut_{norm}$ represent the corresponding total efficiencies for signal and normalisation channels, respectively. In the case of the ratio of branching fractions ${\cal B}({{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}{\mu^{+}\mu^{-}})$/${\cal B}({{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}\phi(\rightarrow{\mu^{+}\mu^{-}}))$, the corresponding normalisation factor $\alpha$ takes the form ${\rm\epsilon\mathstrut_{norm}}/({\rm\epsilon\mathstrut_{sig}\times\it{N}_{\rm norm})}$. The ratio of efficiencies ${\rm\epsilon\mathstrut_{norm}}/{\rm\epsilon\mathstrut_{sig}}$ for the signal regions is compatible with unity within statistical uncertainty.

An extended unbinned maximum-likelihood fit to the ${p}{\mu^{+}}{\mu^{-}}$ invariant mass is performed in each $m_{\mu\mu}$ region, where the signal shape of ${\mathchar 28931\relax}^{+}_{c}$ candidates is described by a Crystal Ball function [33]. The background is described using two distinct components: the combinatorial background, modeled by an exponential shape, and the ${{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}\pi^{+}\pi^{-}$ peaking background, represented by the sum of a Crystal Ball function and a Gaussian function The shapes of both the signal and the ${{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}\pi^{+}\pi^{-}$ background are determined from simulation for each $m_{\mu\mu}$ region and fixed in the fit to data. Weights are applied to simulated ${{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}\pi^{+}\pi^{-}$ events to match the phase-space distribution observed in Ref. [34].

In each $m_{\mu\mu}$ region, there are four parameters in the fit. Three of them are free parameters: the signal yield, the combinatorial background yield and the exponent describing the shape of the combinatorial background. The fourth one corresponds to the ${{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}\pi^{+}\pi^{-}$ background. Its total yield in all regions is a free parameter, and the fraction in each region is constrained based on simulation.

The results of the fit to the data for the signal region and for the $low\text{-}m$ and $high\text{-}m$ regions are presented in Fig. 1. The results of the fit for the normalisation mode and other regions dominated by resonances are presented in Fig. 2. The yields and the ${\mathchar 28931\relax}^{+}_{c}$ signal significance for the six $m_{\mu\mu}$ regions are reported in Table 2, where the significance is determined using Wilks’ theorem [35]. For the signal regions the upper limits for the corresponding branching fractions are determined. The residual contribution from resonant decays into signal regions is estimated from simulation neglecting interference effects. For the $high\text{-}m$ region, where the significance is the highest, only 2.8 events are attributed to ${{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}\phi(\rightarrow{\mu^{+}\mu^{-}})$ decays while for the $low\text{-}m$ region the expected contribution is negligible.

Figure 1: Distributions of ${p}{\mu^{+}\mu^{-}}$ invariant mass for (top) signal, (bottom left) ${low\text{-}m}$ and (bottom right) ${high\text{-}m}$ regions. The results of the fit are overlaid.

Figure 2: Distributions of ${p}{\mu^{+}\mu^{-}}$ invariant mass for (top left) $\phi$, (top right) $\omega$, (bottom left) $\rho$ and (bottom right) $\eta$ regions. The results of the fit are overlaid.

Table 2: Results of the fit to the ${p}{\mu^{+}}{\mu^{-}}$ invariant-mass distribution for various $m_{\mu\mu}$ regions. The significance of the ${{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}{\mu^{+}\mu^{-}}$ component is also given.

Region	${{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}{\mu^{+}\mu^{-}}$	${{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}\pi^{+}\pi^{-}$	Combinatorial	Significance
	yield	yield	yield	${{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}{\mu^{+}\mu^{-}}$
$signal$	$\phantom{0}18\pm 10$	$\phantom{0,}3\pm 7\phantom{0,}$	$681\pm 28$	$2.0\sigma$
${low\text{-}m}$	$\phantom{00}1\pm 5\phantom{0}$	$\phantom{0,}4\pm 4\phantom{0,}$	$241\pm 17$	$0.3\sigma$
${high\text{-}m}$	$\phantom{0}21\pm 8\phantom{0}$	$\phantom{0,}4\pm 4\phantom{0,}$	$432\pm 22$	$2.8\sigma$
$\eta$	$\phantom{0}12\pm 5\phantom{0}$	$2.2\pm 1.6$	$\phantom{0}84\pm 10$	$3.0\sigma$
$\rho$	$\phantom{0}43\pm 10$	$\phantom{,}20\pm 6\phantom{0,}$	$382\pm 22$	$5.6\sigma$
$\omega$	$\phantom{0}81\pm 10$	$4.8\pm 2.1$	$101\pm 11$	$>7\sigma$
$\phi$	$423\pm 22$	$3.8\pm 2.4$	$173\pm 15$	$>7\sigma$

5 Systematic uncertainties

The systematic uncertainties taken into account in the determination of the upper limit on the branching fraction are summarised in Table 3. The two main contributions are the statistical uncertainty of the simulation samples used to determine the efficiency ratio and the uncertainty on the yield of the normalisation mode. The latter is dominated by the statistical uncertainty on the yield returned from the fit. Other sources of systematic uncertainty on the efficiency ratio are associated with:

Table 3: Summary of systematic uncertainties for the signal regions and for the $\eta$, $\rho$ and $\omega$ resonant regions.

Uncertainty source [%]	signal	$low\text{-}m$	$high\text{-}m$	$\eta$	$\rho$	$\omega$
Normalisation channel	5.3	5.3	5.3	5.3	5.3	5.3
Efficiency ratio (stat.)	5.7	8.2	6.8	12.9	5.3	4.6
Efficiency ratio (syst.)	1.7	1.7	1.7	1.7	1.7	1.7
Shape of signal	1.6	3.4	1.1	0.6	1.7	1.2
Shape of ${{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}\pi^{+}\pi^{-}$	0.1	3.2	0.2	0.4	5.3	1.2
Shape of combinatorial	0.1	0.3	0.2	0.1	0.0	0.0
${{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}\pi^{+}\pi^{-}$ decay model	0.1	0.1	0.1	0.1	0.1	0.1
Fit bias	0.1	0.1	0.2	0.3	0.1	0.1
Total	8.1	11.0	8.9	14.1	9.5	7.4

• the level of agreement between data and simulation. The maximum difference between ratios of efficiencies determined for a few variants of the baseline procedure of weighting simulation to match the data is taken as a systematic uncertainty;

• the efficiency determination from simulation. The difference in the efficiency obtained including or excluding ${\mathchar 28931\relax}^{+}_{c}$ candidates for which one or more decay products are not correctly matched to the generated particles, is taken as a systematic uncertainty. The procedure is applied consistently for the signal and normalisation modes;

• the correction of PID variables in simulation. The difference of the efficiency ratio for the baseline and an alternative parameterisation of the PID calibration procedure [32] is taken as a systematic uncertainty.

These three sources added in quadrature amount to a relative uncertainty of 1.7%.

The systematic uncertainties related to the fit model are estimated using pseudoexperiments. An alternative model for each component is used to generate pseudodata, which are then fitted using the baseline model described in Sec. 4. A Johnson function [36] is used as an alternative shape for the signal, while the ${{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}\pi^{+}\pi^{-}$ background is parameterised as the sum of a Johnson and a Gaussian function, and a linear function is used as an alternative model for the combinatorial background. The difference between the generated signal yield and the average fitted yield is taken as systematic uncertainty. In addition, the uncertainty related to the ${{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}\pi^{+}\pi^{-}$ decay model is determined by applying different corrections to the ${{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}\pi^{+}\pi^{-}$ simulated sample and determining an alternative invariant-mass shape. The bias due to maximum-likelihood fit is determined from the same data model used for pseudodata generation and fitting. The signal yield difference between the generated value and the average of fitted values from the pseudoexperiments is taken as a systematic uncertainty.

The systematic uncertainties associated with the ratio of branching fractions in the resonant regions are determined in similar way as for the signal regions and are presented in Table 3.

6 Results

6.1 Signal regions

No significant contribution of ${{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}{\mu^{+}\mu^{-}}$ decays is observed in either of the two signal regions. Upper limits on the branching fractions and their ratios are determined using the CL_s method implemented in the GammaCombo package [37]. The systematic uncertainties are included in the determination of CL_s. The following upper limit in the signal region is determined to be

$$\frac{{\cal B}({{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}{\mu^{+}\mu^{-}})}{{\cal B}({{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}\phi){\cal B}(\phi\!\rightarrow{\mu^{+}\mu^{-}})}<0.09~{}(0.10)~{}~{}{\rm at}~{}~{}90\%~{}(95\%)~{}{\rm CL}.$$

Using the values of the branching fractions for ${{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}\phi$ and $\phi\!\rightarrow{\mu^{+}\mu^{-}}$ decays from Ref. [38] and including their uncertainties in the CL_s determination, an upper limit on the branching fraction is determined to be

$${\cal B}({{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}{\mu^{+}\mu^{-}})<2.9~{}(3.2)\times 10^{-8}~{}~{}~{}~{}{\rm at}~{}~{}90\%~{}(95\%)~{}{\rm CL}.$$

The CL_s as a function of branching fraction is shown in Fig. 3.

Figure 3: Observed and expected (background only hypothesis) limits for ${{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}{\mu^{+}\mu^{-}}$ branching fraction.

The limits are determined for $low\text{-}m$ and $high\text{-}m$ regions separately with the same assumptions as for the signal region. The following upper limits are determined:

$$\frac{{\cal B}({{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}{\mu^{+}\mu^{-}})}{{\cal B}({{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}\phi){\cal B}(\phi\!\rightarrow{\mu^{+}\mu^{-}})}<0.029~{}(0.034)~{}~{}{\rm at}~{}~{}90\%~{}(95\%)~{}{\rm CL}~{}({low\text{-}m}),$$

$$\frac{{\cal B}({{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}{\mu^{+}\mu^{-}})}{{\cal B}({{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}\phi){\cal B}(\phi\!\rightarrow{\mu^{+}\mu^{-}})}<0.094~{}(0.10)~{}~{}{\rm at}~{}~{}90\%~{}(95\%)~{}{\rm CL}~{}({high\text{-}m}),$$

$$\vspace{0.5cm}{\cal B}({{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}{\mu^{+}\mu^{-}})<0.93~{}(1.1)\times 10^{-8}~{}~{}~{}~{}{\rm at}~{}~{}90\%~{}(95\%)~{}{\rm CL}~{}({low\text{-}m}),$$

$${\cal B}({{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}{\mu^{+}\mu^{-}})<3.0~{}(3.3)\times 10^{-8}~{}~{}~{}~{}{\rm at}~{}~{}90\%~{}(95\%)~{}{\rm CL}~{}({high\text{-}m}).$$

The upper limit placed using data from the $high\text{-}m$ region is driven by the signal yield whose significance is close to but lower than 3$\sigma$, while in the $low\text{-}m$ region the invariant-mass distribution is fully compatible with no signal.

Assuming a decay model with final state particles equally distributed in the phase space, the values of upper limits extrapolated from the signal region to the full phase space are determined to be

$$\frac{{\cal B}({{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}{\mu^{+}\mu^{-}})}{{\cal B}({{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}\phi){\cal B}(\phi\!\rightarrow{\mu^{+}\mu^{-}})}<0.23~{}(0.25)~{}~{}{\rm at}~{}~{}90\%~{}(95\%)~{}{\rm CL},$$

$${\cal B}({{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}{\mu^{+}\mu^{-}})<7.3~{}(8.2)\times 10^{-8}~{}~{}~{}~{}{\rm at}~{}~{}90\%~{}(95\%)~{}{\rm CL}.$$

In spite of the enlarged data sample, the extrapolated upper limit on the branching fraction is only slightly better than that determined for the 2011–2012 data [9], $7.7~{}(9.6)\times 10^{-8}~{}{\rm at}~{}~{}90\%~{}(95\%)~{}{\rm CL}$. This relatively small improvement is driven mainly by the 2.8$\sigma$ effect observed in the $high\text{-}m$ region in this analysis.

6.2 Resonant regions

Excluding the normalisation channel, the significance of the ${\mathchar 28931\relax}^{+}_{c}$ signal exceeds the 5$\sigma$ threshold in the $\rho$ and $\omega$ regions while the significance for the ${{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}\eta(\rightarrow{\mu^{+}\mu^{-}})$ decay is at $3\sigma$. The corresponding branching fractions are determined using Eqs. 1,2 where the decay parameters for the nonresonant decays are replaced by those for the resonant decays. It has to be emphasised that disentangling the resonant components to account for interference effects has not been attempted. The efficiency for a given $m_{\mu\mu}$ region is determined using simulation samples of the corresponding resonant decay, thus neglecting possible contributions from other resonances. Under this assumption, the ratio of branching fractions with respect to the normalisation channel are determined for the $\omega$, $\rho$ and $\eta$ regions to be

$$\frac{{\cal B}({{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}\omega(\rightarrow{\mu^{+}\mu^{-}}))}{{\cal B}({{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}\phi(\rightarrow{\mu^{+}\mu^{-}}))}=0.240\pm 0.030~{}(\rm stat.)\pm 0.018~{}(\rm syst.),$$

$$\frac{{\cal B}({{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}\rho(\rightarrow{\mu^{+}\mu^{-}}))}{{\cal B}({{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}\phi(\rightarrow{\mu^{+}\mu^{-}}))}=0.229\pm 0.051~{}(\rm stat.)\pm 0.022~{}(\rm syst.),$$

$$\frac{{\cal B}({{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}\eta(\rightarrow{\mu^{+}\mu^{-}}))}{{\cal B}({{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}\phi(\rightarrow{\mu^{+}\mu^{-}}))}=0.032\pm 0.013~{}(\rm stat.)\pm 0.004~{}(\rm syst.).$$

The corresponding branching fractions are determined to be

$${\cal B}({{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}\omega)=(9.82\pm 1.23~{}(\rm stat.)\pm 0.73~{}(\rm syst.)\pm 2.79~{}(ext.))\times 10^{-4},$$

$${\cal B}({{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}\rho)=(1.52\pm 0.34~{}(\rm stat.)\pm 0.14~{}(\rm syst.)\pm 0.24~{}(ext.))\times 10^{-3},$$

$${\cal B}({{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}\eta)=(1.67\pm 0.69~{}(\rm stat.)\pm 0.23~{}(\rm syst.)\pm 0.34~{}(ext.))\times 10^{-3},$$

where the second uncertainties are induced by the ratio of branching fractions presented in Table 3 and the third uncertainties are due to the following branching fractions: ${\cal B}({{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}\phi)$, ${\cal B}(\phi\!\rightarrow{\mu^{+}\mu^{-}})$, ${\cal B}(\omega\!\rightarrow{\mu^{+}\mu^{-}})$, ${\cal B}(\rho\!\rightarrow{\mu^{+}}{\mu^{-}})$ and ${\cal B}(\eta\!\rightarrow{\mu^{+}}{\mu^{-}})$ [38], used as external inputs.

7 Summary

A search for the nonresonant ${{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}{\mu^{+}\mu^{-}}$ decay is presented using data samples of proton-proton collisions collected by the LHCb experiment in 2016–2018, corresponding to an integrated luminosity of $5.4\text{\,fb}^{-1}$. In the dimuon invariant-mass regions where the contributions of resonant decays is subdominant, i.e. below the $\eta$ and above the $\phi$ resonances, the following upper limits for the branching fractions are derived: ${\cal B}({{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}{\mu^{+}\mu^{-}})<0.93~{}(1.1)\times 10^{-8}$ at 90% (95%) CL and ${\cal B}({{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}{\mu^{+}\mu^{-}})<3.0~{}(3.3)\times 10^{-8}$ at 90% (95%) CL, respectively where ${{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}\phi(\rightarrow{\mu^{+}\mu^{-}})$ as a normalisation channel. Combining these two regions, one finds ${\cal B}({{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}{\mu^{+}\mu^{-}})<2.9~{}(3.2)\times 10^{-8}$ at 90% (95%) CL. The results for the short-distance contribution presented in this paper should not be combined with the results published for 2011–2012 data [9].

In a narrow region around the $\omega$ resonance, a ${\mathchar 28931\relax}^{+}_{c}$ signal is observed with a significance of more than 7$\sigma$, corresponding to a branching fraction $(9.82\pm 1.23~{}(\rm stat.)\pm 0.73~{}(\rm syst.)\pm 2.79~{}(ext.))\times 10^{-4}$, compatible with and more precise than the result of the previous analysis for the data collected in 2011–2012 [9].

Excluding the region around the $\omega$ meson, a signal in the $\rho$ region is found with a significance of 5.6$\sigma$. The corresponding branching fraction determined with respect to the normalisation channel is measured to be $(1.52\pm 0.34~{}(\rm stat.)\pm 0.14~{}(\rm syst.)\pm 0.24~{}(ext.))\times 10^{-3}$, giving strong evidence for the ${{\mathchar 28931\relax}^{+}_{c}}\!\rightarrow{p}\rho$ decay, qualified by the fact that interference effects are not accounted for. For the $\eta$ region the branching fraction was measured to be $(1.67\pm 0.69~{}(\rm stat.)\pm 0.23~{}(\rm syst.)\pm 0.34~{}(ext.))\times 10^{-3}$ with a significance of 3.0$\sigma$. The value is compatible with the world average [38] but less precise.

LHCb collaboration

R. Aaij³⁶ , A.S.W. Abdelmotteleb⁵⁵ , C. Abellan Beteta⁴⁹, F. Abudinén⁵⁵ , T. Ackernley⁵⁹ , A. A.  Adefisoye⁶⁷ , B. Adeva⁴⁵ , M. Adinolfi⁵³ , P. Adlarson⁷⁹ , C. Agapopoulou¹³ , C.A. Aidala⁸⁰ , Z. Ajaltouni¹¹, S. Akar⁶⁴ , K. Akiba³⁶ , P. Albicocco²⁶ , J. Albrecht¹⁸ , F. Alessio⁴⁷ , M. Alexander⁵⁸ , Z. Aliouche⁶¹ , P. Alvarez Cartelle⁵⁴ , R. Amalric¹⁵ , S. Amato³ , J.L. Amey⁵³ , Y. Amhis^13,47 , L. An⁶ , L. Anderlini²⁵ , M. Andersson⁴⁹ , A. Andreianov⁴² , P. Andreola⁴⁹ , M. Andreotti²⁴ , D. Andreou⁶⁷ , A. Anelli^29,p , D. Ao⁷ , F. Archilli^35,v , M. Argenton²⁴ , S. Arguedas Cuendis⁹ , A. Artamonov⁴² , M. Artuso⁶⁷ , E. Aslanides¹² , M. Atzeni⁶³ , B. Audurier¹⁴ , D. Bacher⁶² , I. Bachiller Perea¹⁰ , S. Bachmann²⁰ , M. Bachmayer⁴⁸ , J.J. Back⁵⁵ , P. Baladron Rodriguez⁴⁵ , V. Balagura¹⁴ , W. Baldini²⁴ , H.  Bao⁷ , J. Baptista de Souza Leite⁵⁹ , M. Barbetti^25,m , I. R. Barbosa⁶⁸ , R.J. Barlow⁶¹ , M. Barnyakov²³ , S. Barsuk¹³ , W. Barter⁵⁷ , M. Bartolini⁵⁴ , J. Bartz⁶⁷ , F. Baryshnikov⁴² , J.M. Basels¹⁶ , G. Bassi³³ , B. Batsukh⁵ , A. Bay⁴⁸ , A. Beck⁵⁵ , M. Becker¹⁸ , F. Bedeschi³³ , I.B. Bediaga² , S. Belin⁴⁵ , V. Bellee⁴⁹ , K. Belous⁴² , I. Belov²⁷ , I. Belyaev³⁴ , G. Benane¹² , G. Bencivenni²⁶ , E. Ben-Haim¹⁵ , A. Berezhnoy⁴² , R. Bernet⁴⁹ , S. Bernet Andres⁴³ , A. Bertolin³¹ , C. Betancourt⁴⁹ , F. Betti⁵⁷ , J.  Bex⁵⁴ , Ia. Bezshyiko⁴⁹ , J. Bhom³⁹ , M.S. Bieker¹⁸ , N.V. Biesuz²⁴ , P. Billoir¹⁵ , A. Biolchini³⁶ , M. Birch⁶⁰ , F.C.R. Bishop¹⁰ , A. Bitadze⁶¹ , A. Bizzeti , T. Blake⁵⁵ , F. Blanc⁴⁸ , J.E. Blank¹⁸ , S. Blusk⁶⁷ , V. Bocharnikov⁴² , J.A. Boelhauve¹⁸ , O. Boente Garcia¹⁴ , T. Boettcher⁶⁴ , A.  Bohare⁵⁷ , A. Boldyrev⁴² , C.S. Bolognani⁷⁶ , R. Bolzonella^24,l , N. Bondar⁴² , F. Borgato^31,q,47 , S. Borghi⁶¹ , M. Borsato^29,p , J.T. Borsuk³⁹ , S.A. Bouchiba⁴⁸ , T.J.V. Bowcock⁵⁹ , A. Boyer⁴⁷ , C. Bozzi²⁴ , M.J. Bradley⁶⁰, A. Brea Rodriguez⁴⁸ , N. Breer¹⁸ , J. Brodzicka³⁹ , A. Brossa Gonzalo⁴⁵ , J. Brown⁵⁹ , D. Brundu³⁰ , E. Buchanan⁵⁷, A. Buonaura⁴⁹ , L. Buonincontri^31,q , A.T. Burke⁶¹ , C. Burr⁴⁷ , A. Butkevich⁴² , J.S. Butter⁵⁴ , J. Buytaert⁴⁷ , W. Byczynski⁴⁷ , S. Cadeddu³⁰ , H. Cai⁷², R. Calabrese^24,l , S. Calderon Ramirez⁹ , L. Calefice⁴⁴ , S. Cali²⁶ , M. Calvi^29,p , M. Calvo Gomez⁴³ , P. Camargo Magalhaes^2,z , J. I. Cambon Bouzas⁴⁵ , P. Campana²⁶ , D.H. Campora Perez⁷⁶ , A.F. Campoverde Quezada⁷ , S. Capelli²⁹ , L. Capriotti²⁴ , R. Caravaca-Mora⁹ , A. Carbone^23,j , L. Carcedo Salgado⁴⁵ , R. Cardinale^27,n , A. Cardini³⁰ , P. Carniti^29,p , L. Carus²⁰, A. Casais Vidal⁶³ , R. Caspary²⁰ , G. Casse⁵⁹ , J. Castro Godinez⁹ , M. Cattaneo⁴⁷ , G. Cavallero^24,47 , V. Cavallini^24,l , S. Celani²⁰ , D. Cervenkov⁶² , S.  Cesare^28,o , A.J. Chadwick⁵⁹ , I. Chahrour⁸⁰ , M. Charles¹⁵ , Ph. Charpentier⁴⁷ , C.A. Chavez Barajas⁵⁹ , M. Chefdeville¹⁰ , C. Chen¹² , S. Chen⁵ , Z. Chen⁷ , A. Chernov³⁹ , S. Chernyshenko⁵¹ , V. Chobanova⁷⁸ , S. Cholak⁴⁸ , M. Chrzaszcz³⁹ , A. Chubykin⁴² , V. Chulikov⁴² , P. Ciambrone²⁶ , X. Cid Vidal⁴⁵ , G. Ciezarek⁴⁷ , P. Cifra⁴⁷ , P.E.L. Clarke⁵⁷ , M. Clemencic⁴⁷ , H.V. Cliff⁵⁴ , J. Closier⁴⁷ , C. Cocha Toapaxi²⁰ , V. Coco⁴⁷ , J. Cogan¹² , E. Cogneras¹¹ , L. Cojocariu⁴¹ , P. Collins⁴⁷ , T. Colombo⁴⁷ , A. Comerma-Montells⁴⁴ , L. Congedo²² , A. Contu³⁰ , N. Cooke⁵⁸ , I. Corredoira ⁴⁵ , A. Correia¹⁵ , G. Corti⁴⁷ , J.J. Cottee Meldrum⁵³, B. Couturier⁴⁷ , D.C. Craik⁴⁹ , M. Cruz Torres^2,g , E. Curras Rivera⁴⁸ , R. Currie⁵⁷ , C.L. Da Silva⁶⁶ , S. Dadabaev⁴² , L. Dai⁶⁹ , X. Dai⁶ , E. Dall’Occo¹⁸ , J. Dalseno⁴⁵ , C. D’Ambrosio⁴⁷ , J. Daniel¹¹ , A. Danilina⁴² , P. d’Argent²² , A.  Davidson⁵⁵ , J.E. Davies⁶¹ , A. Davis⁶¹ , O. De Aguiar Francisco⁶¹ , C. De Angelis^30,k , F. De Benedetti⁴⁷ , J. de Boer³⁶ , K. De Bruyn⁷⁵ , S. De Capua⁶¹ , M. De Cian^20,47 , U. De Freitas Carneiro Da Graca^2,b , E. De Lucia²⁶ , J.M. De Miranda² , L. De Paula³ , M. De Serio^22,h , P. De Simone²⁶ , F. De Vellis¹⁸ , J.A. de Vries⁷⁶ , F. Debernardis²² , D. Decamp¹⁰ , V. Dedu¹² , L. Del Buono¹⁵ , B. Delaney⁶³ , H.-P. Dembinski¹⁸ , J. Deng⁸ , V. Denysenko⁴⁹ , O. Deschamps¹¹ , F. Dettori^30,k , B. Dey⁷⁴ , P. Di Nezza²⁶ , I. Diachkov⁴² , S. Didenko⁴² , S. Ding⁶⁷ , L. Dittmann²⁰ , V. Dobishuk⁵¹ , A. D.  Docheva⁵⁸ , C. Dong⁴ , A.M. Donohoe²¹ , F. Dordei³⁰ , A.C. dos Reis² , A. D.  Dowling⁶⁷ , W. Duan⁷⁰ , P. Duda⁷⁷ , M.W. Dudek³⁹ , L. Dufour⁴⁷ , V. Duk³² , P. Durante⁴⁷ , M. M. Duras⁷⁷ , J.M. Durham⁶⁶ , O. D.  Durmus⁷⁴ , A. Dziurda³⁹ , A. Dzyuba⁴² , S. Easo⁵⁶ , E. Eckstein¹⁷, U. Egede¹ , A. Egorychev⁴² , V. Egorychev⁴² , S. Eisenhardt⁵⁷ , E. Ejopu⁶¹ , L. Eklund⁷⁹ , M. Elashri⁶⁴ , J. Ellbracht¹⁸ , S. Ely⁶⁰ , A. Ene⁴¹ , E. Epple⁶⁴ , J. Eschle⁶⁷ , S. Esen²⁰ , T. Evans⁶¹ , F. Fabiano^30,k,47 , L.N. Falcao² , Y. Fan⁷ , B. Fang⁷² , L. Fantini^32,r , M. Faria⁴⁸ , K.  Farmer⁵⁷ , D. Fazzini^29,p , L. Felkowski⁷⁷ , M. Feng^5,7 , M. Feo^18,47 , M. Fernandez Gomez⁴⁵ , A.D. Fernez⁶⁵ , F. Ferrari²³ , F. Ferreira Rodrigues³ , M. Ferrillo⁴⁹ , M. Ferro-Luzzi⁴⁷ , S. Filippov⁴² , R.A. Fini²² , M. Fiorini^24,l , K.M. Fischer⁶² , D.S. Fitzgerald⁸⁰ , C. Fitzpatrick⁶¹ , F. Fleuret¹⁴ , M. Fontana²³ , L. F.  Foreman⁶¹ , R. Forty⁴⁷ , D. Foulds-Holt⁵⁴ , M. Franco Sevilla⁶⁵ , M. Frank⁴⁷ , E. Franzoso^24,l , G. Frau²⁰ , C. Frei⁴⁷ , D.A. Friday⁶¹ , J. Fu⁷ , Q. Fuehring¹⁸ , Y. Fujii¹ , T. Fulghesu¹⁵ , E. Gabriel³⁶ , G. Galati²² , M.D. Galati³⁶ , A. Gallas Torreira⁴⁵ , D. Galli^23,j , S. Gambetta⁵⁷ , M. Gandelman³ , P. Gandini²⁸ , B.  Ganie⁶¹ , H. Gao⁷ , R. Gao⁶² , Y. Gao⁸ , Y. Gao⁶ , Y. Gao⁸, M. Garau^30,k , L.M. Garcia Martin⁴⁸ , P. Garcia Moreno⁴⁴ , J. García Pardiñas⁴⁷ , K. G.  Garg⁸ , L. Garrido⁴⁴ , C. Gaspar⁴⁷ , R.E. Geertsema³⁶ , L.L. Gerken¹⁸ , E. Gersabeck⁶¹ , M. Gersabeck⁶¹ , T. Gershon⁵⁵ , Z. Ghorbanimoghaddam⁵³, L. Giambastiani^31,q , F. I. Giasemis^15,e , V. Gibson⁵⁴ , H.K. Giemza⁴⁰ , A.L. Gilman⁶² , M. Giovannetti²⁶ , A. Gioventù⁴⁴ , P. Gironella Gironell⁴⁴ , C. Giugliano^24,l , M.A. Giza³⁹ , E.L. Gkougkousis⁶⁰ , F.C. Glaser^13,20 , V.V. Gligorov^15,47 , C. Göbel⁶⁸ , E. Golobardes⁴³ , D. Golubkov⁴² , A. Golutvin^60,42,47 , A. Gomes^2,a,† , S. Gomez Fernandez⁴⁴ , F. Goncalves Abrantes⁶² , M. Goncerz³⁹ , G. Gong⁴ , J. A. Gooding¹⁸ , I.V. Gorelov⁴² , C. Gotti²⁹ , J.P. Grabowski¹⁷ , L.A. Granado Cardoso⁴⁷ , E. Graugés⁴⁴ , E. Graverini^48,t , L. Grazette⁵⁵ , G. Graziani , A. T. Grecu⁴¹ , L.M. Greeven³⁶ , N.A. Grieser⁶⁴ , L. Grillo⁵⁸ , S. Gromov⁴² , C.  Gu¹⁴ , M. Guarise²⁴ , M. Guittiere¹³ , V. Guliaeva⁴² , P. A. Günther²⁰ , A.-K. Guseinov⁴⁸ , E. Gushchin⁴² , Y. Guz^6,42,47 , T. Gys⁴⁷ , K. Habermann¹⁷ , T. Hadavizadeh¹ , C. Hadjivasiliou⁶⁵ , G. Haefeli⁴⁸ , C. Haen⁴⁷ , J. Haimberger⁴⁷ , M. Hajheidari⁴⁷, M.M. Halvorsen⁴⁷ , P.M. Hamilton⁶⁵ , J. Hammerich⁵⁹ , Q. Han⁸ , X. Han²⁰ , S. Hansmann-Menzemer²⁰ , L. Hao⁷ , N. Harnew⁶² , M. Hartmann¹³ , J. He^7,c , F. Hemmer⁴⁷ , C. Henderson⁶⁴ , R.D.L. Henderson^1,55 , A.M. Hennequin⁴⁷ , K. Hennessy⁵⁹ , L. Henry⁴⁸ , J. Herd⁶⁰ , P. Herrero Gascon²⁰ , J. Heuel¹⁶ , A. Hicheur³ , G. Hijano Mendizabal⁴⁹, D. Hill⁴⁸ , S.E. Hollitt¹⁸ , J. Horswill⁶¹ , R. Hou⁸ , Y. Hou¹¹ , N. Howarth⁵⁹, J. Hu²⁰, J. Hu⁷⁰ , W. Hu⁶ , X. Hu⁴ , W. Huang⁷ , W. Hulsbergen³⁶ , R.J. Hunter⁵⁵ , M. Hushchyn⁴² , D. Hutchcroft⁵⁹ , D. Ilin⁴² , P. Ilten⁶⁴ , A. Inglessi⁴² , A. Iniukhin⁴² , A. Ishteev⁴² , K. Ivshin⁴² , R. Jacobsson⁴⁷ , H. Jage¹⁶ , S.J. Jaimes Elles^46,73 , S. Jakobsen⁴⁷ , E. Jans³⁶ , B.K. Jashal⁴⁶ , A. Jawahery^65,47 , V. Jevtic¹⁸ , E. Jiang⁶⁵ , X. Jiang^5,7 , Y. Jiang⁷ , Y. J.  Jiang⁶ , M. John⁶² , D. Johnson⁵² , C.R. Jones⁵⁴ , T.P. Jones⁵⁵ , S. Joshi⁴⁰ , B. Jost⁴⁷ , N. Jurik⁴⁷ , I. Juszczak³⁹ , D. Kaminaris⁴⁸ , S. Kandybei⁵⁰ , Y. Kang⁴ , C. Kar¹¹ , M. Karacson⁴⁷ , D. Karpenkov⁴² , A. Kauniskangas⁴⁸ , J.W. Kautz⁶⁴ , F. Keizer⁴⁷ , M. Kenzie⁵⁴ , T. Ketel³⁶ , B. Khanji⁶⁷ , A. Kharisova⁴² , S. Kholodenko^33,47 , G. Khreich¹³ , T. Kirn¹⁶ , V.S. Kirsebom^29,p , O. Kitouni⁶³ , S. Klaver³⁷ , N. Kleijne^33,s , K. Klimaszewski⁴⁰ , M.R. Kmiec⁴⁰ , S. Koliiev⁵¹ , L. Kolk¹⁸ , A. Konoplyannikov⁴² , P. Kopciewicz^38,47 , P. Koppenburg³⁶ , M. Korolev⁴² , I. Kostiuk³⁶ , O. Kot⁵¹, S. Kotriakhova , A. Kozachuk⁴² , P. Kravchenko⁴² , L. Kravchuk⁴² , M. Kreps⁵⁵ , P. Krokovny⁴² , W. Krupa⁶⁷ , W. Krzemien⁴⁰ , O.K. Kshyvanskyi⁵¹, J. Kubat²⁰, S. Kubis⁷⁷ , M. Kucharczyk³⁹ , V. Kudryavtsev⁴² , E. Kulikova⁴² , A. Kupsc⁷⁹ , B. K.  Kutsenko¹² , D. Lacarrere⁴⁷ , A. Lai³⁰ , A. Lampis³⁰ , D. Lancierini⁵⁴ , C. Landesa Gomez⁴⁵ , J.J. Lane¹ , R. Lane⁵³ , C. Langenbruch²⁰ , J. Langer¹⁸ , O. Lantwin⁴² , T. Latham⁵⁵ , F. Lazzari^33,t , C. Lazzeroni⁵² , R. Le Gac¹² , R. Lefèvre¹¹ , A. Leflat⁴² , S. Legotin⁴² , M. Lehuraux⁵⁵ , E. Lemos Cid⁴⁷ , O. Leroy¹² , T. Lesiak³⁹ , B. Leverington²⁰ , A. Li⁴ , H. Li⁷⁰ , K. Li⁸ , L. Li⁶¹ , P. Li⁴⁷ , P.-R. Li⁷¹ , Q.  Li^5,7 , S. Li⁸ , T. Li^5,d , T. Li⁷⁰ , Y. Li⁸, Y. Li⁵ , Z. Lian⁴ , X. Liang⁶⁷ , S. Libralon⁴⁶ , C. Lin⁷ , T. Lin⁵⁶ , R. Lindner⁴⁷ , V. Lisovskyi⁴⁸ , R. Litvinov³⁰ , F. L.  Liu¹ , G. Liu⁷⁰ , K. Liu⁷¹ , S. Liu^5,7 , Y. Liu⁵⁷ , Y. Liu⁷¹, Y. L.  Liu⁶⁰ , A. Lobo Salvia⁴⁴ , A. Loi³⁰ , J. Lomba Castro⁴⁵ , T. Long⁵⁴ , J.H. Lopes³ , A. Lopez Huertas⁴⁴ , S. López Soliño⁴⁵ , C. Lucarelli^25,m , D. Lucchesi^31,q , M. Lucio Martinez⁷⁶ , V. Lukashenko^36,51 , Y. Luo⁶ , A. Lupato³¹ , E. Luppi^24,l , K. Lynch²¹ , X.-R. Lyu⁷ , G. M.  Ma⁴ , R. Ma⁷ , S. Maccolini¹⁸ , F. Machefert¹³ , F. Maciuc⁴¹ , B.  Mack⁶⁷ , I. Mackay⁶² , L. M.  Mackey⁶⁷ , L.R. Madhan Mohan⁵⁴ , M. M.  Madurai⁵² , A. Maevskiy⁴² , D. Magdalinski³⁶ , D. Maisuzenko⁴² , M.W. Majewski³⁸, J.J. Malczewski³⁹ , S. Malde⁶² , L. Malentacca⁴⁷, A. Malinin⁴² , T. Maltsev⁴² , G. Manca^30,k , G. Mancinelli¹² , C. Mancuso^28,13,o , R. Manera Escalero⁴⁴, D. Manuzzi²³ , D. Marangotto^28,o , J.F. Marchand¹⁰ , R. Marchevski⁴⁸ , U. Marconi²³ , S. Mariani⁴⁷ , C. Marin Benito⁴⁴ , J. Marks²⁰ , A.M. Marshall⁵³ , G. Martelli^32,r , G. Martellotti³⁴ , L. Martinazzoli⁴⁷ , M. Martinelli^29,p , D. Martinez Santos⁴⁵ , F. Martinez Vidal⁴⁶ , A. Massafferri² , R. Matev⁴⁷ , A. Mathad⁴⁷ , V. Matiunin⁴² , C. Matteuzzi⁶⁷ , K.R. Mattioli¹⁴ , A. Mauri⁶⁰ , E. Maurice¹⁴ , J. Mauricio⁴⁴ , P. Mayencourt⁴⁸ , M. Mazurek⁴⁰ , M. McCann⁶⁰ , L. Mcconnell²¹ , T.H. McGrath⁶¹ , N.T. McHugh⁵⁸ , A. McNab⁶¹ , R. McNulty²¹ , B. Meadows⁶⁴ , G. Meier¹⁸ , D. Melnychuk⁴⁰ , F. M.  Meng⁴ , M. Merk^36,76 , A. Merli⁴⁸ , L. Meyer Garcia⁶⁵ , D. Miao^5,7 , H. Miao⁷ , M. Mikhasenko^17,f , D.A. Milanes⁷³ , A. Minotti^29,p , E. Minucci⁶⁷ , T. Miralles¹¹ , B. Mitreska¹⁸ , D.S. Mitzel¹⁸ , A. Modak⁵⁶ , A. Mödden ¹⁸ , R.A. Mohammed⁶² , R.D. Moise¹⁶ , S. Mokhnenko⁴² , T. Mombächer⁴⁷ , M. Monk^55,1 , S. Monteil¹¹ , A. Morcillo Gomez⁴⁵ , G. Morello²⁶ , M.J. Morello^33,s , M.P. Morgenthaler²⁰ , A.B. Morris⁴⁷ , A.G. Morris¹² , R. Mountain⁶⁷ , H. Mu⁴ , Z. M.  Mu⁶ , E. Muhammad⁵⁵ , F. Muheim⁵⁷ , M. Mulder⁷⁵ , K. Müller⁴⁹ , F. Muñoz-Rojas⁹ , R. Murta⁶⁰ , P. Naik⁵⁹ , T. Nakada⁴⁸ , R. Nandakumar⁵⁶ , T. Nanut⁴⁷ , I. Nasteva³ , M. Needham⁵⁷ , N. Neri^28,o , S. Neubert¹⁷ , N. Neufeld⁴⁷ , P. Neustroev⁴², J. Nicolini^18,13 , D. Nicotra⁷⁶ , E.M. Niel⁴⁸ , N. Nikitin⁴² , P. Nogarolli³ , P. Nogga¹⁷, N.S. Nolte⁶³ , C. Normand⁵³ , J. Novoa Fernandez⁴⁵ , G. Nowak⁶⁴ , C. Nunez⁸⁰ , H. N.  Nur⁵⁸ , A. Oblakowska-Mucha³⁸ , V. Obraztsov⁴² , T. Oeser¹⁶ , S. Okamura^24,l,47 , A. Okhotnikov⁴², O. Okhrimenko⁵¹ , R. Oldeman^30,k , F. Oliva⁵⁷ , M. Olocco¹⁸ , C.J.G. Onderwater⁷⁶ , R.H. O’Neil⁵⁷ , J.M. Otalora Goicochea³ , P. Owen⁴⁹ , A. Oyanguren⁴⁶ , O. Ozcelik⁵⁷ , K.O. Padeken¹⁷ , B. Pagare⁵⁵ , P.R. Pais²⁰ , T. Pajero⁴⁷ , A. Palano²² , M. Palutan²⁶ , G. Panshin⁴² , L. Paolucci⁵⁵ , A. Papanestis⁵⁶ , M. Pappagallo^22,h , L.L. Pappalardo^24,l , C. Pappenheimer⁶⁴ , C. Parkes⁶¹ , B. Passalacqua²⁴ , G. Passaleva²⁵ , D. Passaro^33,s , A. Pastore²² , M. Patel⁶⁰ , J. Patoc⁶² , C. Patrignani^23,j , A.  Paul⁶⁷ , C.J. Pawley⁷⁶ , A. Pellegrino³⁶ , J.  Peng^5,7 , M. Pepe Altarelli²⁶ , S. Perazzini²³ , D. Pereima⁴² , H.  Pereira Da Costa⁶⁶ , A. Pereiro Castro⁴⁵ , P. Perret¹¹ , A. Perro⁴⁷ , K. Petridis⁵³ , A. Petrolini^27,n , J. P.  Pfaller⁶⁴ , H. Pham⁶⁷ , L. Pica³³ , M. Piccini³² , B. Pietrzyk¹⁰ , G. Pietrzyk¹³ , D. Pinci³⁴ , F. Pisani⁴⁷ , M. Pizzichemi^29,p , V. Placinta⁴¹ , M. Plo Casasus⁴⁵ , F. Polci^15,47 , M. Poli Lener²⁶ , A. Poluektov¹² , N. Polukhina⁴² , I. Polyakov⁴⁷ , E. Polycarpo³ , S. Ponce⁴⁷ , D. Popov⁷ , S. Poslavskii⁴² , K. Prasanth⁵⁷ , C. Prouve⁴⁵ , V. Pugatch⁵¹ , G. Punzi^33,t , S.  Qasim⁴⁹ , W. Qian⁷ , N. Qin⁴ , S. Qu⁴ , R. Quagliani⁴⁷ , R.I. Rabadan Trejo⁵⁵ , J.H. Rademacker⁵³ , M. Rama³³ , M.  Ramírez García⁸⁰ , V. Ramos De Oliveira⁶⁸ , M. Ramos Pernas⁵⁵ , M.S. Rangel³ , F. Ratnikov⁴² , G. Raven³⁷ , M. Rebollo De Miguel⁴⁶ , F. Redi^28,i , J. Reich⁵³ , F. Reiss⁶¹ , Z. Ren⁷ , P.K. Resmi⁶² , R. Ribatti^33,s , G. R.  Ricart^14,81 , D. Riccardi^33,s , S. Ricciardi⁵⁶ , K. Richardson⁶³ , M. Richardson-Slipper⁵⁷ , K. Rinnert⁵⁹ , P. Robbe¹³ , G. Robertson⁵⁸ , E. Rodrigues⁵⁹ , E. Rodriguez Fernandez⁴⁵ , J.A. Rodriguez Lopez⁷³ , E. Rodriguez Rodriguez⁴⁵ , A. Rogovskiy⁵⁶ , D.L. Rolf⁴⁷ , P. Roloff⁴⁷ , V. Romanovskiy⁴² , M. Romero Lamas⁴⁵ , A. Romero Vidal⁴⁵ , G. Romolini²⁴ , F. Ronchetti⁴⁸ , M. Rotondo²⁶ , S. R.  Roy²⁰ , M.S. Rudolph⁶⁷ , T. Ruf⁴⁷ , M. Ruiz Diaz²⁰ , R.A. Ruiz Fernandez⁴⁵ , J. Ruiz Vidal^79,aa , A. Ryzhikov⁴² , J. Ryzka³⁸ , J. J. Saavedra-Arias⁹ , J.J. Saborido Silva⁴⁵ , R. Sadek¹⁴ , N. Sagidova⁴² , D. Sahoo⁷⁴ , N. Sahoo⁵² , B. Saitta^30,k , M. Salomoni^29,p,47 , C. Sanchez Gras³⁶ , I. Sanderswood⁴⁶ , R. Santacesaria³⁴ , C. Santamarina Rios⁴⁵ , M. Santimaria^26,47 , L. Santoro ² , E. Santovetti³⁵ , A. Saputi^24,47 , D. Saranin⁴² , G. Sarpis⁵⁷ , M. Sarpis⁶¹ , C. Satriano^34,u , A. Satta³⁵ , M. Saur⁶ , D. Savrina⁴² , H. Sazak¹⁶ , L.G. Scantlebury Smead⁶² , A. Scarabotto¹⁸ , S. Schael¹⁶ , S. Scherl⁵⁹ , M. Schiller⁵⁸ , H. Schindler⁴⁷ , M. Schmelling¹⁹ , B. Schmidt⁴⁷ , S. Schmitt¹⁶ , H. Schmitz¹⁷, O. Schneider⁴⁸ , A. Schopper⁴⁷ , N. Schulte¹⁸ , S. Schulte⁴⁸ , M.H. Schune¹³ , R. Schwemmer⁴⁷ , G. Schwering¹⁶ , B. Sciascia²⁶ , A. Sciuccati⁴⁷ , S. Sellam⁴⁵ , A. Semennikov⁴² , T. Senger⁴⁹ , M. Senghi Soares³⁷ , A. Sergi²⁷ , N. Serra⁴⁹ , L. Sestini³¹ , A. Seuthe¹⁸ , Y. Shang⁶ , D.M. Shangase⁸⁰ , M. Shapkin⁴² , R. S.  Sharma⁶⁷ , I. Shchemerov⁴² , L. Shchutska⁴⁸ , T. Shears⁵⁹ , L. Shekhtman⁴² , Z. Shen⁶ , S. Sheng^5,7 , V. Shevchenko⁴² , B. Shi⁷ , Q. Shi⁷ , E.B. Shields^29,p , Y. Shimizu¹³ , E. Shmanin⁴² , R. Shorkin⁴² , J.D. Shupperd⁶⁷ , R. Silva Coutinho⁶⁷ , G. Simi^31,q , S. Simone^22,h , N. Skidmore⁵⁵ , T. Skwarnicki⁶⁷ , M.W. Slater⁵² , J.C. Smallwood⁶² , E. Smith⁶³ , K. Smith⁶⁶ , M. Smith⁶⁰ , A. Snoch³⁶ , L. Soares Lavra⁵⁷ , M.D. Sokoloff⁶⁴ , F.J.P. Soler⁵⁸ , A. Solomin^42,53 , A. Solovev⁴² , I. Solovyev⁴² , R. Song¹ , Y. Song⁴⁸ , Y. Song⁴ , Y. S.  Song⁶ , F.L. Souza De Almeida⁶⁷ , B. Souza De Paula³ , E. Spadaro Norella^28,o , E. Spedicato²³ , J.G. Speer¹⁸ , E. Spiridenkov⁴², P. Spradlin⁵⁸ , V. Sriskaran⁴⁷ , F. Stagni⁴⁷ , M. Stahl⁴⁷ , S. Stahl⁴⁷ , S. Stanislaus⁶² , E.N. Stein⁴⁷ , O. Steinkamp⁴⁹ , O. Stenyakin⁴², H. Stevens¹⁸ , D. Strekalina⁴² , Y. Su⁷ , F. Suljik⁶² , J. Sun³⁰ , L. Sun⁷² , Y. Sun⁶⁵ , D. S.  Sundfeld Lima², W. Sutcliffe⁴⁹, P.N. Swallow⁵² , F. Swystun⁵⁴ , A. Szabelski⁴⁰ , T. Szumlak³⁸ , Y. Tan⁴ , M.D. Tat⁶² , A. Terentev⁴² , F. Terzuoli^33,w , F. Teubert⁴⁷ , E. Thomas⁴⁷ , D.J.D. Thompson⁵² , H. Tilquin⁶⁰ , V. Tisserand¹¹ , S. T’Jampens¹⁰ , M. Tobin⁵ , L. Tomassetti^24,l , G. Tonani^28,o,47 , X. Tong⁶ , D. Torres Machado² , L. Toscano¹⁸ , D.Y. Tou⁴ , C. Trippl⁴³ , G. Tuci²⁰ , N. Tuning³⁶ , L.H. Uecker²⁰ , A. Ukleja³⁸ , D.J. Unverzagt²⁰ , E. Ursov⁴² , A. Usachov³⁷ , A. Ustyuzhanin⁴² , U. Uwer²⁰ , V. Vagnoni²³ , A. Valassi⁴⁷ , G. Valenti²³ , N. Valls Canudas⁴⁷ , H. Van Hecke⁶⁶ , E. van Herwijnen⁶⁰ , C.B. Van Hulse^45,y , R. Van Laak⁴⁸ , M. van Veghel³⁶ , G. Vasquez⁴⁹ , R. Vazquez Gomez⁴⁴ , P. Vazquez Regueiro⁴⁵ , C. Vázquez Sierra⁴⁵ , S. Vecchi²⁴ , J.J. Velthuis⁵³ , M. Veltri^25,x , A. Venkateswaran⁴⁸ , M. Vesterinen⁵⁵ , M. Vieites Diaz⁴⁷ , X. Vilasis-Cardona⁴³ , E. Vilella Figueras⁵⁹ , A. Villa²³ , P. Vincent¹⁵ , F.C. Volle⁵² , D. vom Bruch¹² , N. Voropaev⁴² , K. Vos⁷⁶ , G. Vouters^10,47 , C. Vrahas⁵⁷ , J. Wagner¹⁸ , J. Walsh³³ , E.J. Walton^1,55 , G. Wan⁶ , C. Wang²⁰ , G. Wang⁸ , J. Wang⁶ , J. Wang⁵ , J. Wang⁴ , J. Wang⁷² , M. Wang²⁸ , N. W.  Wang⁷ , R. Wang⁵³ , X. Wang⁸, X. Wang⁷⁰ , X. W.  Wang⁶⁰ , Z. Wang¹³ , Z. Wang⁴ , Z. Wang²⁸ , J.A. Ward^55,1 , M. Waterlaat⁴⁷, N.K. Watson⁵² , D. Websdale⁶⁰ , Y. Wei⁶ , J. Wendel⁷⁸ , B.D.C. Westhenry⁵³ , D.J. White⁶¹ , M. Whitehead⁵⁸ , A.R. Wiederhold⁵⁵ , D. Wiedner¹⁸ , G. Wilkinson⁶² , M.K. Wilkinson⁶⁴ , M. Williams⁶³ , M.R.J. Williams⁵⁷ , R. Williams⁵⁴ , F.F. Wilson⁵⁶ , W. Wislicki⁴⁰ , M. Witek³⁹ , L. Witola²⁰ , C.P. Wong⁶⁶ , G. Wormser¹³ , S.A. Wotton⁵⁴ , H. Wu⁶⁷ , J. Wu⁸ , Y. Wu⁶ , K. Wyllie⁴⁷ , S. Xian⁷⁰, Z. Xiang⁵ , Y. Xie⁸ , A. Xu³³ , J. Xu⁷ , L. Xu⁴ , L. Xu⁴ , M. Xu⁵⁵ , Z. Xu¹¹ , Z. Xu⁷ , Z. Xu⁵ , D. Yang , K.  Yang⁶⁰ , S. Yang⁷ , X. Yang⁶ , Y. Yang^27,n , Z. Yang⁶ , Z. Yang⁶⁵ , V. Yeroshenko¹³ , H. Yeung⁶¹ , H. Yin⁸ , C. Y.  Yu⁶ , J. Yu⁶⁹ , X. Yuan⁵ , E. Zaffaroni⁴⁸ , M. Zavertyaev¹⁹ , M. Zdybal³⁹ , C.  Zeng^5,7 , M. Zeng⁴ , C. Zhang⁶ , D. Zhang⁸ , J. Zhang⁷ , L. Zhang⁴ , S. Zhang⁶⁹ , S. Zhang⁶ , Y. Zhang⁶ , Y. Z.  Zhang⁴ , Y. Zhao²⁰ , A. Zharkova⁴² , A. Zhelezov²⁰ , X. Z.  Zheng⁴ , Y. Zheng⁷ , T. Zhou⁶ , X. Zhou⁸ , Y. Zhou⁷ , V. Zhovkovska⁵⁵ , L. Z.  Zhu⁷ , X. Zhu⁴ , X. Zhu⁸ , V. Zhukov¹⁶ , J. Zhuo⁴⁶ , Q. Zou^5,7 , D. Zuliani^31,q , G. Zunica⁴⁸ .

¹School of Physics and Astronomy, Monash University, Melbourne, Australia
²Centro Brasileiro de Pesquisas Físicas (CBPF), Rio de Janeiro, Brazil
³Universidade Federal do Rio de Janeiro (UFRJ), Rio de Janeiro, Brazil
⁴Center for High Energy Physics, Tsinghua University, Beijing, China
⁵Institute Of High Energy Physics (IHEP), Beijing, China
⁶School of Physics State Key Laboratory of Nuclear Physics and Technology, Peking University, Beijing, China
⁷University of Chinese Academy of Sciences, Beijing, China
⁸Institute of Particle Physics, Central China Normal University, Wuhan, Hubei, China
⁹Consejo Nacional de Rectores (CONARE), San Jose, Costa Rica
¹⁰Université Savoie Mont Blanc, CNRS, IN2P3-LAPP, Annecy, France
¹¹Université Clermont Auvergne, CNRS/IN2P3, LPC, Clermont-Ferrand, France
¹²Aix Marseille Univ, CNRS/IN2P3, CPPM, Marseille, France
¹³Université Paris-Saclay, CNRS/IN2P3, IJCLab, Orsay, France
¹⁴Laboratoire Leprince-Ringuet, CNRS/IN2P3, Ecole Polytechnique, Institut Polytechnique de Paris, Palaiseau, France
¹⁵LPNHE, Sorbonne Université, Paris Diderot Sorbonne Paris Cité, CNRS/IN2P3, Paris, France
¹⁶I. Physikalisches Institut, RWTH Aachen University, Aachen, Germany
¹⁷Universität Bonn - Helmholtz-Institut für Strahlen und Kernphysik, Bonn, Germany
¹⁸Fakultät Physik, Technische Universität Dortmund, Dortmund, Germany
¹⁹Max-Planck-Institut für Kernphysik (MPIK), Heidelberg, Germany
²⁰Physikalisches Institut, Ruprecht-Karls-Universität Heidelberg, Heidelberg, Germany
²¹School of Physics, University College Dublin, Dublin, Ireland
²²INFN Sezione di Bari, Bari, Italy
²³INFN Sezione di Bologna, Bologna, Italy
²⁴INFN Sezione di Ferrara, Ferrara, Italy
²⁵INFN Sezione di Firenze, Firenze, Italy
²⁶INFN Laboratori Nazionali di Frascati, Frascati, Italy
²⁷INFN Sezione di Genova, Genova, Italy
²⁸INFN Sezione di Milano, Milano, Italy
²⁹INFN Sezione di Milano-Bicocca, Milano, Italy
³⁰INFN Sezione di Cagliari, Monserrato, Italy
³¹INFN Sezione di Padova, Padova, Italy
³²INFN Sezione di Perugia, Perugia, Italy
³³INFN Sezione di Pisa, Pisa, Italy
³⁴INFN Sezione di Roma La Sapienza, Roma, Italy
³⁵INFN Sezione di Roma Tor Vergata, Roma, Italy
³⁶Nikhef National Institute for Subatomic Physics, Amsterdam, Netherlands
³⁷Nikhef National Institute for Subatomic Physics and VU University Amsterdam, Amsterdam, Netherlands
³⁸AGH - University of Krakow, Faculty of Physics and Applied Computer Science, Kraków, Poland
³⁹Henryk Niewodniczanski Institute of Nuclear Physics Polish Academy of Sciences, Kraków, Poland
⁴⁰National Center for Nuclear Research (NCBJ), Warsaw, Poland
⁴¹Horia Hulubei National Institute of Physics and Nuclear Engineering, Bucharest-Magurele, Romania
⁴²Affiliated with an institute covered by a cooperation agreement with CERN
⁴³DS4DS, La Salle, Universitat Ramon Llull, Barcelona, Spain
⁴⁴ICCUB, Universitat de Barcelona, Barcelona, Spain
⁴⁵Instituto Galego de Física de Altas Enerxías (IGFAE), Universidade de Santiago de Compostela, Santiago de Compostela, Spain
⁴⁶Instituto de Fisica Corpuscular, Centro Mixto Universidad de Valencia - CSIC, Valencia, Spain
⁴⁷European Organization for Nuclear Research (CERN), Geneva, Switzerland
⁴⁸Institute of Physics, Ecole Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland
⁴⁹Physik-Institut, Universität Zürich, Zürich, Switzerland
⁵⁰NSC Kharkiv Institute of Physics and Technology (NSC KIPT), Kharkiv, Ukraine
⁵¹Institute for Nuclear Research of the National Academy of Sciences (KINR), Kyiv, Ukraine
⁵²University of Birmingham, Birmingham, United Kingdom
⁵³H.H. Wills Physics Laboratory, University of Bristol, Bristol, United Kingdom
⁵⁴Cavendish Laboratory, University of Cambridge, Cambridge, United Kingdom
⁵⁵Department of Physics, University of Warwick, Coventry, United Kingdom
⁵⁶STFC Rutherford Appleton Laboratory, Didcot, United Kingdom
⁵⁷School of Physics and Astronomy, University of Edinburgh, Edinburgh, United Kingdom
⁵⁸School of Physics and Astronomy, University of Glasgow, Glasgow, United Kingdom
⁵⁹Oliver Lodge Laboratory, University of Liverpool, Liverpool, United Kingdom
⁶⁰Imperial College London, London, United Kingdom
⁶¹Department of Physics and Astronomy, University of Manchester, Manchester, United Kingdom
⁶²Department of Physics, University of Oxford, Oxford, United Kingdom
⁶³Massachusetts Institute of Technology, Cambridge, MA, United States
⁶⁴University of Cincinnati, Cincinnati, OH, United States
⁶⁵University of Maryland, College Park, MD, United States
⁶⁶Los Alamos National Laboratory (LANL), Los Alamos, NM, United States
⁶⁷Syracuse University, Syracuse, NY, United States
⁶⁸Pontifícia Universidade Católica do Rio de Janeiro (PUC-Rio), Rio de Janeiro, Brazil, associated to ^3
69School of Physics and Electronics, Hunan University, Changsha City, China, associated to ^8
70Guangdong Provincial Key Laboratory of Nuclear Science, Guangdong-Hong Kong Joint Laboratory of Quantum Matter, Institute of Quantum Matter, South China Normal University, Guangzhou, China, associated to ^4
71Lanzhou University, Lanzhou, China, associated to ^5
72School of Physics and Technology, Wuhan University, Wuhan, China, associated to ^4
73Departamento de Fisica , Universidad Nacional de Colombia, Bogota, Colombia, associated to ^15
74Eotvos Lorand University, Budapest, Hungary, associated to ^47
75Van Swinderen Institute, University of Groningen, Groningen, Netherlands, associated to ^36
76Universiteit Maastricht, Maastricht, Netherlands, associated to ^36
77Tadeusz Kosciuszko Cracow University of Technology, Cracow, Poland, associated to ^39
78Universidade da Coruña, A Coruna, Spain, associated to ^43
79Department of Physics and Astronomy, Uppsala University, Uppsala, Sweden, associated to ^58
80University of Michigan, Ann Arbor, MI, United States, associated to ^67
81Departement de Physique Nucleaire (SPhN), Gif-Sur-Yvette, France

^aUniversidade de Brasília, Brasília, Brazil
^bCentro Federal de Educacão Tecnológica Celso Suckow da Fonseca, Rio De Janeiro, Brazil
^cHangzhou Institute for Advanced Study, UCAS, Hangzhou, China
^dSchool of Physics and Electronics, Henan University , Kaifeng, China
^eLIP6, Sorbonne Université, Paris, France
^fExcellence Cluster ORIGINS, Munich, Germany
^gUniversidad Nacional Autónoma de Honduras, Tegucigalpa, Honduras
^hUniversità di Bari, Bari, Italy
ⁱUniversita degli studi di Bergamo, Bergamo, Italy
^jUniversità di Bologna, Bologna, Italy
^kUniversità di Cagliari, Cagliari, Italy
^lUniversità di Ferrara, Ferrara, Italy
^mUniversità di Firenze, Firenze, Italy
ⁿUniversità di Genova, Genova, Italy
^oUniversità degli Studi di Milano, Milano, Italy
^pUniversità degli Studi di Milano-Bicocca, Milano, Italy
^qUniversità di Padova, Padova, Italy
^rUniversità di Perugia, Perugia, Italy
^sScuola Normale Superiore, Pisa, Italy
^tUniversità di Pisa, Pisa, Italy
^uUniversità della Basilicata, Potenza, Italy
^vUniversità di Roma Tor Vergata, Roma, Italy
^wUniversità di Siena, Siena, Italy
^xUniversità di Urbino, Urbino, Italy
^yUniversidad de Alcalá, Alcalá de Henares , Spain
^zFacultad de Ciencias Fisicas, Madrid, Spain
^aaDepartment of Physics/Division of Particle Physics, Lund, Sweden
^†Deceased