Search for lepton flavour violating decays of heavy resonances and quantum black holes to an \mathrm{e}\mu pair in proton-proton collisions at \sqrt{s}=8\,\text{TeV}

A search for narrow resonances decaying to an electron and a muon is presented. The mass spectrum is also investigated for non-resonant contributions from the production of quantum black holes (QBHs). The analysis is performed using data corresponding to an integrated luminosity of 19.7 collected in proton-proton collisions at a centre-of-mass energy of 8 with the CMS detector at the LHC. With no evidence for physics beyond the standard model in the invariant mass spectrum of selected pairs, upper limits are set at 95 confidence level on the product of cross section and branching fraction for signals arising in theories with charged lepton flavour violation. In the search for narrow resonances, the resonant production of a sneutrino in R-parity violating supersymmetry is considered. The sneutrino is excluded for masses below 1.28 for couplings , and below 2.30 for and . These are the most stringent limits to date from direct searches at high-energy colliders. In addition, the resonance searches are interpreted in terms of a model with heavy partners of the boson and the photon. In a framework of TeV-scale quantum gravity based on a renormalization of Newton’s constant, the search for non-resonant contributions to the mass spectrum excludes QBH production below a threshold mass of 1.99. In models that invoke extra dimensions, the bounds range from 2.36 for one extra dimension to 3.63 for six extra dimensions. This is the first search for QBHs decaying into the final state.


CERN-EP/2016-068 2019/\two@digits8/\two@digits4


Search for lepton flavour violating decays of heavy resonances and quantum black holes to an pair in proton-proton collisions at

The CMS Collaboration111See Appendix A for the list of collaboration members


Please replace the default abstract using the abstract command.

Published in the European Physical Journal C as doi:10.1140/epjc/s10052-016-4149-y.

© 2019 CERN for the benefit of the CMS Collaboration. CC-BY-3.0 license

1 Introduction

Several extensions of the standard model (SM) predict the existence of heavy, short-lived states that decay to the final state, and motivate the search for lepton flavour violating (LFV) signatures in interactions involving charged leptons. This paper reports a search for phenomena beyond the SM in the invariant mass spectrum of pairs. The analysis is based on data with an integrated luminosity of 19.7 collected in proton-proton () collisions at with the CMS detector at the CERN LHC [1]. The results are interpreted in terms of three theoretically predicted objects: a sneutrino () lightest supersymmetric particle (LSP) in R-parity violating (RPV) supersymmetry (SUSY) [2], interfering LFV and bosons [3], and quantum black holes (QBHs) [4, 5, 6].

In RPV SUSY, lepton number can be violated at tree level in interactions between fermions and sfermions, and the may be the LSP [7]. For the resonant signal, the following trilinear RPV part of the superpotential is considered: (), where , , and are generation indices, and are the doublet superfields of the leptons and quarks, and and are the singlet superfields of the charged leptons and down-like quarks. We assume that all RPV couplings vanish, except for , , and , and consider a SUSY mass hierarchy with a LSP. In this model, the can be produced resonantly in collisions via the coupling and it can decay either into an pair via the and couplings, or into a pair via the coupling. In this analysis we consider only the final state and, for simplicity, we assume .

The LFV signal is based on a model with two extra dimensions [8, 3], where the three generations of the SM arise from a single generation in higher-dimensional space-time. Flavour changing processes are introduced through the Kaluza-Klein modes of gauge fields that are not localised on a brane. In four-dimensional space-time, an effective Lagrangian can be obtained that contains two complex vector fields and . These vector fields generate transitions between the families in which the generation number changes by unity, such as the process and its charge conjugate. The structure of the terms in the Lagrangian for the production and decay of the and bosons is analogous to that describing the interactions of the boson and the photon with quarks and charged leptons, respectively. The coupling strengths and are related to their SM counterparts through a multiplicative coupling modifier . For simplicity, the masses and are assumed to be equal, and the model is referred to as the LFV model. It is characterized by the two independent parameters and .

Theories that have a fundamental Planck scale of the order of a TeV [9, 10, 11, 12, 13] offer the possibility of producing microscopic black holes [14, 15, 16] at the LHC. In contrast to semiclassical, thermal black holes, which would decay to high-multiplicity final states, QBHs are non-thermal objects expected to decay predominantly to pairs of particles. We consider the production of a , colourless, neutral QBH in a model with lepton flavour violation, in which the cross section for QBH production is extrapolated from semiclassical black holes and depends on the threshold mass for QBH production and the number of extra dimensions . For , it corresponds to a 3+1-dimensional model with low-scale quantum gravity, where a renormalization of Newton’s constant leads to a Planck scale at the TeV scale [13, 17, 18]; corresponds to the Randall–Sundrum (RS) brane world model [9, 10]; and to the Arkani-Hamed–Dimopoulos–Dvali (ADD) model [11, 12]. We consider flat-space black holes (black holes that are spherical both in the brane and in the bulk dimensions) and, in the case of RS-type black holes (), consider only the regime in which almost flat five-dimensional space is an applicable metric. This is the case for , where is the Schwarzschild radius, denotes the Anti-de Sitter curvature, and is the size of the extra dimension. The threshold is assumed to be at the Planck scale in the definition of the Particle Data Group [19] for and , whereas for both the PDG and RS definitions [4] are adopted. In this model, the branching fraction of QBH decays to the final state is 1.1%, which is twice that of the dimuon or dielectron decay modes, making the signature the most promising leptonic decay channel. While the resonant and LFV signals result in a narrow peak in the invariant mass spectrum of the pair, the mass distribution of the QBH signal is characterized by an edge at the threshold for QBH production, and a monotonically decreasing tail.

Direct searches for resonances in the invariant mass spectrum with interpretations in terms of production have been carried out by the CDF [20] and D0 [21] collaborations at the Fermilab Tevatron and most recently by the ATLAS collaboration [22] using collision data at a centre-of-mass energy of at the LHC. For couplings and , the most stringent of these limits stems from the search performed by the ATLAS collaboration, excluding at 95% confidence level (CL) a below a mass of . Low-energy muon conversion experiments [23] yield strong limits as a function of the sneutrino mass on the product of the two RPV couplings of at 90% CL [24]. In the case of the signal, searches for decays constrain the coupling modifier . For the choice , a bound of is obtained at 90% CL [25, 3]. There have been searches for QBHs decaying hadronically, by the CMS [26, 27, 28] and ATLAS [29, 30] collaborations, and in the photon plus jet, lepton plus jet, dimuon, and dielectron final states, by the ATLAS collaboration [31, 32, 33, 34]. This is the first search for QBH decays into the final state.

The search for the phenomena beyond the SM described above is carried out for invariant masses of the pair of , which is the relevant region in light of existing constraints from other direct searches. Using the same event selection, the  invariant mass spectrum is searched for two different signal shapes: the shape associated with a narrow resonance that may be interpreted in terms of any model involving a resonance decaying promptly into an electron and a muon, and the more model-specific QBH signal shape. With a relative invariant mass resolution ranging from 1.6% at to 6% at , the CMS detector is a powerful tool for searches for new physics in the  invariant mass spectrum.

2 The CMS detector

The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diameter, providing a magnetic field of 3.8 T. Within the solenoid volume are a silicon pixel and strip tracker, a lead tungstate crystal electromagnetic calorimeter (ECAL), and a brass and scintillator hadron calorimeter (HCAL), each composed of a barrel and two endcap sections. Extensive forward calorimetry complements the coverage provided by the barrel and endcap detectors. Muons are measured in gas-ionization detectors embedded in the steel flux-return yoke outside the solenoid. The silicon tracker consists of 1440 silicon pixel and 15 148 silicon strip detector modules and measures charged particles within the pseudorapidity range . The ECAL consists of 75 848 lead tungstate crystals and provides coverage for in a barrel region and in two endcap regions. Muons are measured in the range , with detection planes using three technologies: drift tubes, cathode strip chambers, and resistive plate chambers. A two-level trigger system is used by the CMS experiment. The first level is composed of custom hardware processors and uses information from the calorimeters and muon detectors to select interesting events and to reduce the event rate from the initial bunch crossing frequency of 20 MHz to a maximum of 100 kHz. The high-level trigger processor farm further decreases the event rate to 400 Hz before data storage. A detailed description of the CMS detector, together with a definition of the coordinate system used and the relevant kinematic variables, can be found in Ref. [35].

3 Event selection

The search is designed in a model-independent way by requiring only one prompt, isolated muon and one prompt, isolated electron in the event selection. This minimal selection allows for a reinterpretation of the results in terms of models with more complex event topologies than the single pair present in the signals considered in this paper.

The data sample is selected using a single-muon trigger with a minimum transverse momentum () requirement of . In order to allow the trigger to remain unprescaled, the pseudorapidity of the muons is constrained to values . Offline, each event is required to have a reconstructed collision vertex with at least four associated tracks, located less than 2 from the centre of the detector in the plane transverse to the beam and less than 24 from it in the direction along the beam. The primary vertex is defined as the vertex with the largest sum of squared transverse momenta of its associated tracks.

The reconstruction and identification of electrons and muons is carried out using standard CMS algorithms, described in more detail in Refs. [36, 37, 38, 39, 40]. Reconstruction of the muon track starts from two tracks, one built in the silicon tracker and one built in the muon system. Hits used to reconstruct the tracks in the two systems are then used to reconstruct a track spanning over the entire detector [36]. Muon candidates are required to have a transverse momentum of with a measured uncertainty of and must fall into the acceptance of the trigger of . The candidate’s track must have transverse and longitudinal impact parameters with respect to the primary vertex position of less than 0.2 and 0.5, respectively. At least one hit in the pixel detector, six or more hits in silicon-strip tracker layers, and matched segments in at least two muon detector planes are required to be associated with the reconstructed track. In order to suppress backgrounds from muons within jets, the scalar sum of all other tracks within a cone of size 0.3 in (where is the azimuthal angle in radians) around the muon candidate’s track is required to be less than 10% of the candidate’s .

In the electron reconstruction, ECAL clusters are matched to silicon pixel detector hits, which are then used as seeds for the reconstruction of tracks in the tracker. Electron candidates are built from clusters with associated tracks and must lie within the barrel or endcap acceptance regions, with pseudorapidities of 1.442 and 1.56 2.5, respectively, with a transverse energy . The transverse energy is defined as the magnitude of the projection on the plane perpendicular to the beam of the electron momentum vector normalized to the electron energy measured in the ECAL. Misidentification of jets as electrons is suppressed by requiring that the scalar sum of the of all other tracks in a cone of size 0.3 in around the electron candidate’s track is less than 5. In addition, the sum of the of calorimeter energy deposits in the same cone that are not associated with the electron candidate must be less than 3% of the candidate’s (plus a small -dependent offset). To minimise the impact of additional interactions in the same bunch crossing (pileup) on the selection efficiency, the calorimeter isolation is corrected for the average energy density in the event  [41]. Further reduction of electron misidentification is achieved by requiring the transverse profile of the energy deposition in the ECAL to be consistent with the expected electron profile, and the sum of HCAL energy deposits in a cone of size 0.15 in to be less than 5% of the electron’s ECAL energy. The transverse impact parameter of the electron candidate’s track with respect to the primary vertex must not exceed 0.02 cm and 0.05 cm, for barrel and endcap candidates, respectively, and the track must not have more than one missing hit in the layers of the pixel detector it crossed.

The trigger efficiency has been measured using the “tag-and-probe” technique in dimuon events from Z decays described in  [38, 36, 39]. The trigger efficiency for muons that pass the selection requirements is within , within , and within . The muon identification efficiency, including the isolation requirement, is measured with the tag-and-probe technique applied to muons from boson decays using tracks in the inner silicon tracker as probes. The same efficiency of (syst) is obtained in the three pseudorapidity regions , , and , with corresponding efficiency ratios between data and the simulation of (syst), (syst), and (syst). A range up to has been probed with the tag-and-probe method and the muon identification efficiencies remain constant within the statistical precision, as do the corresponding efficiency ratios between data and simulation. The evolution of the muon reconstruction and identification efficiencies and the muon trigger efficiency for muon is based on simulation. Using dielectron events from boson decays [37], the total efficiency to reconstruct and select electrons with is found to be % (syst) in the barrel region and  (syst) in the endcaps. According to Monte Carlo (MC) simulation, the variation of these efficiencies with electron is less than 1% in the barrel and 2% in the endcaps. The corresponding efficiency ratios for between data and simulation are (syst) in the barrel and (syst) in the endcaps. These efficiencies and efficiency ratios have been measured up to an electron of 1 in the barrel and 500 in the endcap regions.

In the event selection, at least one isolated muon and one isolated electron that both pass the identification criteria described above are required. After the application of all efficiency scale factors that correct the simulation to the efficiencies measured in data, the combined dilepton reconstruction and identification efficiency for RPV signal events within the detector acceptance is expected to be % at and the full selection efficiency including the trigger requirement is 71.2%. The MC simulation predicts that this efficiency is constant within 3% for masses between 200 and 3 . The electron and the muon are not required to have opposite charge, in order to avoid a loss in signal efficiency due to possible electron charge misidentification at high electron . Since highly energetic muons can produce bremsstrahlung resulting in an associated supercluster in the calorimeter in the direction of the muon’s inner track, they can be misidentified as electrons. Therefore, an electron candidate is rejected if there is a muon with greater than within of the candidate. Only one pair per event is considered. For about of the events passing the event selection there is more than one pair in the event, in which case the pair with the highest invariant mass is selected.

4 Signal simulation

The RPV and QBH signal samples are generated with the CalcHEP (v. 3.4.1) event generator [42]. A cross section calculation at next-to-leading order (NLO) in perturbative QCD is used for the RPV signal [43], in which the factorization and renormalization scales are set to and the CTEQ6M [44] set of parton distribution functions (PDF) is used. The invariant mass distributions of reconstructed pairs from simulated QBH signal samples are presented in Fig. 1 for different signal masses and numbers of extra dimensions. A more detailed description of the implemented QBH model including the dependence of the spectrum from QBH decays on the model parameters is presented in Ref. [45]. The LFV signal events are produced with the MadGraph (v. generator [46]. The effects of the interference resulting from the mass degeneracy on the cross section and signal acceptance are taken into account, and the coupling parameters of the model are taken to be the same as in Ref. [3]. All signal samples use the CTEQ6L1 [44] PDF, pythia (v. 6.426) [47] for hadronization with the underlying event tune Z2*, and are processed through a simulation of the full CMS detector based on Geant4 (v. 9.4) [48]. The pythia Z2* tune is derived from the Z1 tune [49], which uses the CTEQ5L PDF set, whereas Z2* adopts CTEQ6L.

The total acceptance times efficiency for each of the three signal models considered in this analysis is determined using MC simulation with selection efficiencies corrected to the values measured in data. The signal acceptance, as defined by the selection on the lepton and applied to the generated leptons in the signal simulation, and the product of acceptance and selection efficiency, are shown in Tables 2 and 4, evaluated for selected signal masses. The acceptance of the RPV model is that of a generic spin-0 resonance. In the case of the LFV model, the acceptance is more model-specific due to the interference between the and the . This interference shapes the distributions of the leptons in the final state, which leads to a smaller acceptance compared to a generic spin-1 resonance. Table 6 lists the parameterizations of the acceptance times efficiency as a function of signal mass for the RPV and LFV resonance signals, resulting from fits in the mass range from to . These parameterizations are used later in the statistical interpretation of the resonance search.

Figure 1: Invariant mass distributions of reconstructed pairs from simulated QBH signal events that pass the event selection, normalized to unit area. The steps at the threshold masses are smeared out by the detector resolution.
(TeV) (TeV)
Table 2: Signal acceptance () and the product of acceptance and efficiency () for different signal masses, for the RPV and LFV models. The acceptance is defined as the fraction of signal events in the simulation passing the selection on lepton and applied to the generated leptons.
(TeV) (TeV)
Table 4: Signal acceptance () and the product of acceptance and efficiency () for different threshold masses , for the QBH models with and extra dimensions. The acceptance is defined as the fraction of signal events in the simulation passing the selection on lepton and applied to the generated leptons.
Model Functional form of
Table 6: Parametrization of the product of signal acceptance and efficiency () as a function of signal mass , for the RPV and LFV models. The value of is expressed in units of .

5 Background estimation

The SM backgrounds contributing to the final state can be divided into two classes of events. The first class comprises events with at least two prompt, isolated leptons. The second class consists of events with either jets or photons that are misidentified as isolated leptons, and events with jets containing non-prompt leptons. This second class of background is referred to as ”non-prompt background” in this paper. The expected SM background from processes with two prompt leptons is obtained from MC simulations. It consists mostly of events from production and production; the former process is dominant at lower masses and the latter becomes equally important above . Other background processes estimated from MC simulation are the additional diboson processes and , single top production, and Drell–Yan (DY) events with subsequent decay of the pair into an electron and a muon. The , , and simulated samples are generated using powheg (v. 1.0) [50, 51, 52] with the CT10 PDF [53], and the DY, , and background samples are generated using the MadGraph (v. event generator with the CTEQ6L1 PDF. All background samples use pythia (v. 6.426) for hadronization with the underlying event tune . The generated events are processed through a full simulation of the CMS detector based on Geant4 (v. 9.4). Pileup interactions are included in the simulation and event-dependent weights are applied in order to reproduce the number of interactions expected for the measured instantaneous luminosity. After this procedure, the distribution of the number of vertices per event observed in data is well described by the simulation. The simulated samples are normalized to the integrated luminosity of the data sample, . The cross sections are calculated to next-to-next-to-leading order (NNLO) accuracy in perturbative QCD for [54] and DY [55] and to NLO accuracy for the  [56], , , and  [57] processes.

The main sources of non-prompt background in the selection arise from +jet and production with a jet or photon that are misidentified as an electron. The +jet, QCD multijet, and processes yield subleading contributions to the background with non-prompt leptons. The background is estimated from simulation based on the MadGraph (v. event generator. A background estimation based on control samples in data, using the jet-to-electron misidentification rate (MR) method explained below, is used to determine the distributions from +jet and QCD multijet production. The measurement of the jet-to-electron misidentification rate has been carried out in the context of Ref. [40]. It starts from a sample collected using a prescaled single electromagnetic cluster trigger, in which the presence of an electron candidate with relaxed electron identification criteria is required. The events of the sample must have no more than one reconstructed electron with , in order to suppress the contribution from decays. The misidentification measurement can be biased by selecting genuine electrons from +jet events or converted photons from +jet events. Processes that can give a single electron, such as , , , , , and   where, if a second electron is produced, it fails to be reconstructed, give another less significant source of contamination. Simulated samples are used to correct for this contamination and its effect on the MR. After these corrections, the electron MR, measured in bins of and , is the number of electrons passing the full selection over the number of electron candidates in the sample.

Using the measured electron MR, the +jet and QCD multijet contributions can be estimated from a sample with a muon passing the single-muon trigger and the full muon selection, and an electron candidate satisfying the relaxed selection requirements but failing the full electron selection. Each event in the sample is weighted by the factor to determine the overall contribution of the jet backgrounds. Contributions from processes other than +jet and QCD multijet are subtracted from the sample to which the MR is applied, to avoid double counting. This subtraction is based on MC simulated background samples. A systematic uncertainty of 30% is applied to the jet background estimate, based on cross-checks and closure tests. An uncertainty of 50% is assigned to the background estimate for the process, which is taken from simulation at leading order (LO) in perturbative QCD.

6 Results

After the event selection, 28 925 events are observed in data. The invariant mass distribution is shown in Fig. 2, together with the corresponding cumulative distribution. A comparison of the observed and expected event yields is given in Table 8. The dominant background process is , which contributes 69% of the total background yield after selection, followed by production, contributing 11%. The two selected leptons carry opposite measured electric charge in 26 840 events and carry the same charge in 2085 events. According to the background estimation, events with same-charge pairs are expected, most of which stem from the +jet process, followed by and diboson production /.

Figure 2: The invariant mass distribution of selected pairs (left), and the corresponding cumulative distribution, where all events above the mass value on the -axis are summed (right). The points with error bars represent the data and the stacked histograms represent the expectations from SM processes. The label ’Jets’ refers to the estimate of the +jet and QCD multijet backgrounds from data. The ratio of the data to the background for each bin is shown at the bottom. The horizontal lines on the data points indicate the bin width.
Total Invariant mass ranges in units of
Total bkg
Table 8: The number of observed events compared to the background expectation in five invariant mass ranges and in the full invariant mass range. The yields obtained from simulations are normalized according to their expected cross sections. The background label ’Jets’ refers to the estimate of the +jet and QCD multijet backgrounds from data.

The systematic uncertainties assigned to backgrounds obtained from simulation include the integrated luminosity (2.6%) [58] and the acceptance times efficiency (5%). The latter is based on the uncertainties in the various efficiency scale factors that correct the simulation to the efficiencies measured in data. According to simulation, the evolution of the lepton selection efficiencies from the pole, where they are measured, to high lepton is covered within this uncertainty. The uncertainty in the muon momentum scale is 5% per . Electron energy scale uncertainties are 0.6% in the barrel and 1.5% in the endcap. These momentum and energy scale uncertainties cumulatively lead to an uncertainty in the total background yield of 2% at and 3.5% at . Uncertainties in the electron and muon resolutions have a negligible impact on the total background yield. The uncertainty associated with the choice of PDF in the background simulation is evaluated according to the PDF4LHC prescription [59, 60] and translates into an uncertainty in the background yield ranging from 5% at to 9% at . Among the uncertainties in the cross sections used for the normalization of the various simulated background samples, the 5% uncertainty in the NNLO QCD cross section of the dominant background [54] is the most relevant. Further uncertainties associated with the modelling of the shape of the invariant mass distribution are taken into account for the two leading backgrounds: (higher-order corrections on the top- description discussed in [61]) and (scale uncertainties studied with the powheg generator). These lead to an uncertainty in the total background yield of up to 13% at . A further systematic uncertainty arises from the limited sizes of the simulated background samples at high invariant mass, where the background expectation is small. Taking all systematic uncertainties into account, the resulting uncertainty in the background yield ranges from 9% at to 18% at .

As shown in the cumulative invariant mass distribution in Fig. 2, we observe a deficit in data compared to the background expectation for . In this invariant mass region, 17 events are observed and the background estimate yields  (syst) events. Combining the systematic and statistical uncertainties, the local significance of this discrepancy is below 2.

No significant excess with respect to the expectation is found in the measured invariant mass distribution, and we set limits on the product of signal cross section and branching fraction for signal mass hypotheses above . Two types of signal shapes are considered for the limit setting: a narrow resonance and the broader invariant mass spectrum from QBH decays. The RPV and signals both result in a narrow resonance. For coupling values not excluded by existing searches, the intrinsic widths of these signals are small compared to the detector resolution. Therefore, Gaussian functions are used to model the signal shapes. For each probed resonance signal mass, the two parameters, acceptance times efficiency (Table 6) and invariant mass resolution, define the signal shape used for limit setting. The invariant mass resolution is derived from fits of Gaussian distributions to the  invariant mass spectra from MC simulated signal samples and ranges from 1.6% at a resonance mass of to 6% at . For high values of pair invariant mass, it is dominated by the resolution on the measurement of the muon , which ranges from about 2% at to 6% at and 10% at . These values are obtained from MC simulations and agree within the uncertainties with measurements using cosmic ray muons. This model of the narrow resonance allows for a scan of the invariant mass spectrum with a fine spacing of the signal mass hypothesis that corresponds to the invariant mass resolution.

Unlike the and signals, the QBH signal exhibits a broader shape with a sharp edge at the threshold mass and a tail towards higher masses (Fig. 1). The QBH signal shapes are obtained directly from simulated samples.

The systematic uncertainties in the signal entering the limit calculation are the 2.6% uncertainty in the integrated luminosity, the 5% uncertainty in the product of acceptance and efficiency, and the relative uncertainty in the mass resolution, which ranges from 2% at to 40% at . The uncertainty in the signal acceptance times efficiency is dominated by the uncertainty in the trigger, lepton reconstruction, and identification efficiencies, and includes the subleading PDF uncertainty in the signal acceptance.

Upper limits at CL on the product of cross section and branching fraction are determined using a binned likelihood Bayesian approach with a positive, uniform prior for the signal cross section [62]. The signal and background shapes enter the likelihood with a binning of 1, well below the invariant mass resolution for masses above 200 GeV. For the resonant signals and , search regions in the invariant mass spectrum are defined as times the invariant mass resolution evaluated at the hypothetical resonance mass. Only events in these search regions enter the binned likelihood in the limit calculation. The impact of a further broadening of the signal window size on the median expected limit has been found to be negligible within the uncertainties. For mass hypotheses above 800, the upper bound of the search region is dropped. In the case of the QBH signal, the search region is defined by a lower bound at , where is the invariant mass resolution, and there is no upper bound. The nuisance parameters associated with the systematic uncertainties are modelled with log-normal distributions, and a Markov Chain MC method is used for integration. For each mass hypothesis considered, the posterior probability density function is derived as a function of the signal cross section times branching fraction and yields the 95% CL upper limit on this parameter of interest.

The CL limits on the signal cross section times branching fraction for the RPV resonance signal are shown in Fig. 3 (left). The signal cross section shown is calculated at NLO in perturbative QCD with the RPV couplings set to and . For these couplings, a lower mass limit of 1.28 is obtained. At this mass, the observed limit on the cross section times branching fraction is 0.25 fb. For a comparison with earlier searches at hadron colliders [20, 22], the two coupling benchmarks , and , are considered. For RPV couplings and , we set a mass limit of 2.30, and improve the lower bound of 2.0 previously set [22]. The lower bound on the signal mass for and is 2.16. In the narrow width approximation, the cross section times branching fraction scales with the RPV couplings as:

Using this relation and the observed upper cross section bounds, we derive the limit contour in the parameter plane as a function of a fixed value of .

Figure 3: Left: The CL upper limit on the product of signal cross section and branching fraction for the RPV signal as a function of the mass of the resonance . Right: The 95% CL limit contours for the RPV signal in the () parameter plane. The values of the parameter are fixed to (red dashed and dotted), (green small-dashed), (blue dashed), and (black solid). The regions above the curves are excluded.
Figure 4: Left: The CL exclusion limit on the product of signal cross section and branching fraction for the signal as a function of the mass . Right: The 95% CL limit contour for the signal in the () parameter plane.
Figure 5: The CL exclusion limit on the product of signal cross section and branching fraction for the QBH signal as a function of the threshold mass . The limits have been calculated using the signal shape of the QBH model without extra dimensions (). For signal masses , the change in the QBH signal shape for different numbers of extra dimensions has a negligible impact on the limit.

For the results presented in Fig. 3 (right), values of the couplings and up to 0.2 and 0.07 are considered, respectively. The ratio of decay width to mass of the sneutrino is less than 0.5% for these coupling values and finite-width effects are small. Searches for resonant dijet production [27, 29] that cover the sneutrino decay to a pair via the coupling do not exclude this region of parameter space. In the model considered here with resonant production of the , we do not reach the sensitivity of muon conversion experiments, which lead to a bound on the coupling product of at 90% CL, assuming . For comparison, with a signal mass of and the assumption , we obtain a limit of at 90% CL. We present results in terms of the product of the production cross section and branching fraction of the that do not depend on a specific production mechanism of the sneutrino.

The CL limits on the signal cross section times branching fraction for the signal, which exhibits a different acceptance from the spin-0 resonance in the RPV model, are presented in Fig. 4 (left). For the coupling modifier , a lower bound on the signal mass of 1.29 is obtained. Figure 4 (right) shows the corresponding limit contour in the parameter plane. Since this resonance is produced dominantly in the initial state, the bound from searches for muon conversion is not as strong as for the RPV signal, but searches for decays yield a stringent exclusion limit of at 90% CL. This can be compared to our bound of at 90% CL for .

In the QBH search, we set limits on the mass threshold for QBH production, , in models with to extra dimensions. The CL limits on the signal cross section times branching fraction for the QBH signal are shown in Fig. 5. For in a model with a Planck scale at the TeV scale from a renormalization of the gravitational constant, we exclude QBH production below a threshold mass of 1.99. For , two signal cross sections are considered with the Schwarzschild radius evaluated in the RS and PDG conventions. The resulting limits on are 2.36 and 2.81, respectively. For ADD-type black holes with , we obtain lower bounds on ranging from 3.15 for to 3.63 for . A summary of the 95% CL lower mass limits set for all signal models is presented in Table 10.

Signal model Lower limit signal mass (TeV)
observed expected
RPV () 1.28 1.24
RPV () 2.16 2.16
RPV () 2.30 2.30
LFV () 1.29 1.25
QBH   1.99 1.99
        QBH   (RS) 2.36 2.36
           QBH   (PDG) 2.81 2.81
QBH   3.15 3.15
QBH   3.34 3.34
QBH   3.46 3.46
QBH   3.55 3.55
QBH   3.63 3.63
Table 10: The CL observed and expected lower bounds on the signal masses of sneutrinos in RPV SUSY, resonances in the LFV model, and QBHs, each with subsequent decay into an pair. For the QBH signal with , two signal cross sections are considered with the Schwarzschild radius evaluated in either the Randall-Sundrum (RS) or the Particle Data Group (PDG) convention.

7 Summary

A search has been reported for heavy states decaying promptly into an electron and a muon using 19.7 of proton-proton collision data recorded with the CMS detector at the LHC at a centre-of-mass energy of . Agreement is observed between the data and the standard model expectation with new limits set on resonant production of sneutrinos in R-parity violating supersymmetry with subsequent decay into pairs. For couplings and , sneutrino lightest supersymmetric particles for masses below are excluded at 95% CL. For couplings and , masses below are excluded. These are the most stringent limits from direct searches at high-energy colliders. For the signal model, a lower mass limit of is set at 95% CL for the coupling modifier . This direct search for resonant production of an pair at the TeV scale does not reach the sensitivity of dedicated low-energy experiments, but complements such indirect searches and can readily be interpreted in terms of different signals of new physics involving a heavy state that decays promptly into an electron and a muon. Lower bounds are set on the mass threshold for the production of quantum black holes with subsequent decay into an pair in models with zero to six extra dimensions, assuming the threshold mass to be at the Planck scale, ranging from to . These are the first limits on quantum black holes decaying into final states.


We congratulate our colleagues in the CERN accelerator departments for the excellent performance of the LHC and thank the technical and administrative staffs at CERN and at other CMS institutes for their contributions to the success of the CMS effort. In addition, we gratefully acknowledge the computing centres and personnel of the Worldwide LHC Computing Grid for delivering so effectively the computing infrastructure essential to our analyses. Finally, we acknowledge the enduring support for the construction and operation of the LHC and the CMS detector provided by the following funding agencies: BMWFW and FWF (Austria); FNRS and FWO (Belgium); CNPq, CAPES, FAPERJ, and FAPESP (Brazil); MES (Bulgaria); CERN; CAS, MoST, and NSFC (China); COLCIENCIAS (Colombia); MSES and CSF (Croatia); RPF (Cyprus); MoER, ERC IUT and ERDF (Estonia); Academy of Finland, MEC, and HIP (Finland); CEA and CNRS/IN2P3 (France); BMBF, DFG, and HGF (Germany); GSRT (Greece); OTKA and NIH (Hungary); DAE and DST (India); IPM (Iran); SFI (Ireland); INFN (Italy); MSIP and NRF (Republic of Korea); LAS (Lithuania); MOE and UM (Malaysia); CINVESTAV, CONACYT, SEP, and UASLP-FAI (Mexico); MBIE (New Zealand); PAEC (Pakistan); MSHE and NSC (Poland); FCT (Portugal); JINR (Dubna); MON, RosAtom, RAS and RFBR (Russia); MESTD (Serbia); SEIDI and CPAN (Spain); Swiss Funding Agencies (Switzerland); MST (Taipei); ThEPCenter, IPST, STAR and NSTDA (Thailand); TUBITAK and TAEK (Turkey); NASU and SFFR (Ukraine); STFC (United Kingdom); DOE and NSF (USA).

Individuals have received support from the Marie-Curie programme and the European Research Council and EPLANET (European Union); the Leventis Foundation; the A. P. Sloan Foundation; the Alexander von Humboldt Foundation; the Belgian Federal Science Policy Office; the Fonds pour la Formation à la Recherche dans l’Industrie et dans l’Agriculture (FRIA-Belgium); the Agentschap voor Innovatie door Wetenschap en Technologie (IWT-Belgium); the Ministry of Education, Youth and Sports (MEYS) of the Czech Republic; the Council of Science and Industrial Research, India; the HOMING PLUS programme of the Foundation for Polish Science, cofinanced from European Union, Regional Development Fund; the Mobility Plus programme of the Ministry of Science and Higher Education (Poland); the OPUS programme of the National Science Center (Poland); MIUR project 20108T4XTM (Italy); the Thalis and Aristeia programmes cofinanced by EU-ESF and the Greek NSRF; the National Priorities Research Program by Qatar National Research Fund; the Rachadapisek Sompot Fund for Postdoctoral Fellowship, Chulalongkorn University (Thailand); the Chulalongkorn Academic into Its 2nd Century Project Advancement Project (Thailand); and the Welch Foundation, contract C-1845.


Appendix A The CMS Collaboration

Yerevan Physics Institute, Yerevan, Armenia
V. Khachatryan, A.M. Sirunyan, A. Tumasyan Institut für Hochenergiephysik der OeAW, Wien, Austria
W. Adam, E. Asilar, T. Bergauer, J. Brandstetter, E. Brondolin, M. Dragicevic, J. Erö, M. Flechl, M. Friedl, R. Frühwirth\@textsuperscript1, V.M. Ghete, C. Hartl, N. Hörmann, J. Hrubec, M. Jeitler\@textsuperscript1, V. Knünz, A. König, M. Krammer\@textsuperscript1, I. Krätschmer, D. Liko, T. Matsushita, I. Mikulec, D. Rabady\@textsuperscript2, N. Rad, B. Rahbaran, H. Rohringer, J. Schieck\@textsuperscript1, R. Schöfbeck, J. Strauss, W. Treberer-Treberspurg, W. Waltenberger, C.-E. Wulz\@textsuperscript1 National Centre for Particle and High Energy Physics, Minsk, Belarus
V. Mossolov, N. Shumeiko, J. Suarez Gonzalez Universiteit Antwerpen, Antwerpen, Belgium
S. Alderweireldt, T. Cornelis, E.A. De Wolf, X. Janssen, A. Knutsson, J. Lauwers, S. Luyckx, M. Van De Klundert, H. Van Haevermaet, P. Van Mechelen, N. Van Remortel, A. Van Spilbeeck Vrije Universiteit Brussel, Brussel, Belgium
S. Abu Zeid, F. Blekman, J. D’Hondt, N. Daci, I. De Bruyn, K. Deroover, N. Heracleous, J. Keaveney, S. Lowette, L. Moreels, A. Olbrechts, Q. Python, D. Strom, S. Tavernier, W. Van Doninck, P. Van Mulders, G.P. Van Onsem, I. Van Parijs Université Libre de Bruxelles, Bruxelles, Belgium
P. Barria, H. Brun, C. Caillol, B. Clerbaux, G. De Lentdecker, W. Fang, G. Fasanella, L. Favart, R. Goldouzian, A. Grebenyuk, G. Karapostoli, T. Lenzi, A. Léonard, T. Maerschalk, A. Marinov, L. Perniè, A. Randle-conde, T. Seva, C. Vander Velde, P. Vanlaer, R. Yonamine, F. Zenoni, F. Zhang\@textsuperscript3 Ghent University, Ghent, Belgium
K. Beernaert, L. Benucci, A. Cimmino, S. Crucy, D. Dobur, A. Fagot, G. Garcia, M. Gul, J. Mccartin, A.A. Ocampo Rios, D. Poyraz, D. Ryckbosch, S. Salva, M. Sigamani, M. Tytgat, W. Van Driessche, E. Yazgan, N. Zaganidis Université Catholique de Louvain, Louvain-la-Neuve, Belgium
S. Basegmez, C. Beluffi\@textsuperscript4, O. Bondu, S. Brochet, G. Bruno, A. Caudron, L. Ceard, C. Delaere, D. Favart, L. Forthomme, A. Giammanco, A. Jafari, P. Jez, M. Komm, V. Lemaitre, A. Mertens, M. Musich, C. Nuttens, L. Perrini, K. Piotrzkowski, A. Popov\@textsuperscript5, L. Quertenmont, M. Selvaggi, M. Vidal Marono Université de Mons, Mons, Belgium
N. Beliy, G.H. Hammad Centro Brasileiro de Pesquisas Fisicas, Rio de Janeiro, Brazil
W.L. Aldá Júnior, F.L. Alves, G.A. Alves, L. Brito, M. Correa Martins Junior, M. Hamer, C. Hensel, A. Moraes, M.E. Pol, P. Rebello Teles Universidade do Estado do Rio de Janeiro, Rio de Janeiro, Brazil
E. Belchior Batista Das Chagas, W. Carvalho, J. Chinellato\@textsuperscript6, A. Custódio, E.M. Da Costa, D. De Jesus Damiao, C. De Oliveira Martins, S. Fonseca De Souza, L.M. Huertas Guativa, H. Malbouisson, D. Matos Figueiredo, C. Mora Herrera, L. Mundim, H. Nogima, W.L. Prado Da Silva, A. Santoro, A. Sznajder, E.J. Tonelli Manganote\@textsuperscript6, A. Vilela Pereira Universidade Estadual Paulista ,  Universidade Federal do ABC ,  São Paulo, Brazil
S. Ahuja, C.A. Bernardes, A. De Souza Santos, S. Dogra, T.R. Fernandez Perez Tomei, E.M. Gregores, P.G. Mercadante, C.S. Moon\@textsuperscript7, S.F. Novaes, Sandra S. Padula, D. Romero Abad, J.C. Ruiz Vargas Institute for Nuclear Research and Nuclear Energy, Sofia, Bulgaria
A. Aleksandrov, R. Hadjiiska, P. Iaydjiev, M. Rodozov, S. Stoykova, G. Sultanov, M. Vutova University of Sofia, Sofia, Bulgaria
A. Dimitrov, I. Glushkov, L. Litov, B. Pavlov, P. Petkov Institute of High Energy Physics, Beijing, China
M. Ahmad, J.G. Bian, G.M. Chen, H.S. Chen, M. Chen, T. Cheng, R. Du, C.H. Jiang, D. Leggat, R. Plestina\@textsuperscript8, F. Romeo, S.M. Shaheen, A. Spiezia, J. Tao, C. Wang, Z. Wang, H. Zhang State Key Laboratory of Nuclear Physics and Technology, Peking University, Beijing, China
C. Asawatangtrakuldee, Y. Ban, Q. Li, S. Liu, Y. Mao, S.J. Qian, D. Wang, Z. Xu Universidad de Los Andes, Bogota, Colombia
C. Avila, A. Cabrera, L.F. Chaparro Sierra, C. Florez, J.P. Gomez, B. Gomez Moreno, J.C. Sanabria University of Split, Faculty of Electrical Engineering, Mechanical Engineering and Naval Architecture, Split, Croatia
N. Godinovic, D. Lelas, I. Puljak, P.M. Ribeiro Cipriano University of Split, Faculty of Science, Split, Croatia
Z. Antunovic, M. Kovac Institute Rudjer Boskovic, Zagreb, Croatia
V. Brigljevic, K. Kadija, J. Luetic, S. Micanovic, L. Sudic University of Cyprus, Nicosia, Cyprus
A. Attikis, G. Mavromanolakis, J. Mousa, C. Nicolaou, F. Ptochos, P.A. Razis, H. Rykaczewski Charles University, Prague, Czech Republic
M. Bodlak, M. Finger\@textsuperscript9, M. Finger Jr.\@textsuperscript9 Academy of Scientific Research and Technology of the Arab Republic of Egypt, Egyptian Network of High Energy Physics, Cairo, Egypt
Y. Assran\@textsuperscript10\@textsuperscript11, S. Elgammal\@textsuperscript10, A. Ellithi Kamel\@textsuperscript12\@textsuperscript12, M.A. Mahmoud\@textsuperscript13\@textsuperscript10 National Institute of Chemical Physics and Biophysics, Tallinn, Estonia
B. Calpas, M. Kadastik, M. Murumaa, M. Raidal, A. Tiko, C. Veelken Department of Physics, University of Helsinki, Helsinki, Finland
P. Eerola, J. Pekkanen, M. Voutilainen Helsinki Institute of Physics, Helsinki, Finland
J. Härkönen, V. Karimäki, R. Kinnunen, T. Lampén, K. Lassila-Perini, S. Lehti, T. Lindén, P. Luukka, T. Peltola, J. Tuominiemi, E. Tuovinen, L. Wendland Lappeenranta University of Technology, Lappeenranta, Finland
J. Talvitie, T. Tuuva DSM/IRFU, CEA/Saclay, Gif-sur-Yvette, France
M. Besancon, F. Couderc, M. Dejardin, D. Denegri, B. Fabbro, J.L. Faure, C. Favaro, F. Ferri, S. Ganjour, A. Givernaud, P. Gras, G. Hamel de Monchenault, P. Jarry, E. Locci, M. Machet, J. Malcles, J. Rander, A. Rosowsky, M. Titov, A. Zghiche Laboratoire Leprince-Ringuet, Ecole Polytechnique, IN2P3-CNRS, Palaiseau, France
A. Abdulsalam, I. Antropov, S. Baffioni, F. Beaudette, P. Busson, L. Cadamuro, E. Chapon, C. Charlot, O. Davignon, N. Filipovic, R. Granier de Cassagnac, M. Jo, S. Lisniak, L. Mastrolorenzo, P. Miné, I.N. Naranjo, M. Nguyen, C. Ochando, G. Ortona, P. Paganini, P. Pigard, S. Regnard, R. Salerno, J.B. Sauvan, Y. Sirois, T. Strebler, Y. Yilmaz, A. Zabi Institut Pluridisciplinaire Hubert Curien, Université de Strasbourg, Université de Haute Alsace Mulhouse, CNRS/IN2P3, Strasbourg, France
J.-L. Agram\@textsuperscript14, J. Andrea, A. Aubin, D. Bloch, J.-M. Brom, M. Buttignol, E.C. Chabert, N. Chanon, C. Collard, E. Conte\@textsuperscript14, X. Coubez, J.-C. Fontaine\@textsuperscript14, D. Gelé, U. Goerlach, C. Goetzmann, A.-C. Le Bihan, J.A. Merlin\@textsuperscript2, K. Skovpen, P. Van Hove Centre de Calcul de l’Institut National de Physique Nucleaire et de Physique des Particules, CNRS/IN2P3, Villeurbanne, France
S. Gadrat Université de Lyon, Université Claude Bernard Lyon 1,  CNRS-IN2P3, Institut de Physique Nucléaire de Lyon, Villeurbanne, France
S. Beauceron, C. Bernet, G. Boudoul, E. Bouvier, C.A. Carrillo Montoya, R. Chierici, D. Contardo, B. Courbon, P. Depasse, H. El Mamouni, J. Fan, J. Fay, S. Gascon, M. Gouzevitch, B. Ille, F. Lagarde, I.B. Laktineh, M. Lethuillier, L. Mirabito, A.L. Pequegnot, S. Perries, J.D. Ruiz Alvarez, D. Sabes, L. Sgandurra, V. Sordini, M. Vander Donckt, P. Verdier, S. Viret Georgian Technical University, Tbilisi, Georgia
T. Toriashvili\@textsuperscript15 Tbilisi State University, Tbilisi, Georgia
L. Rurua RWTH Aachen University, I. Physikalisches Institut, Aachen, Germany
C. Autermann, S. Beranek, L. Feld, A. Heister, M.K. Kiesel, K. Klein, M. Lipinski, A. Ostapchuk, M. Preuten, F. Raupach, S. Schael, J.F. Schulte, T. Verlage, H. Weber, V. Zhukov\@textsuperscript5 RWTH Aachen University, III. Physikalisches Institut A,  Aachen, Germany
M. Ata, M. Brodski, E. Dietz-Laursonn, D. Duchardt, M. Endres, M. Erdmann, S. Erdweg, T. Esch, R. Fischer, A. Güth, T. Hebbeker, C. Heidemann, K. Hoepfner, S. Knutzen, P. Kreuzer, M. Merschmeyer, A. Meyer, P. Millet, S. Mukherjee, M. Olschewski, K. Padeken, P. Papacz, T. Pook, M. Radziej, H. Reithler, M. Rieger, F. Scheuch, L. Sonnenschein, D. Teyssier, S. Thüer RWTH Aachen University, III. Physikalisches Institut B,  Aachen, Germany
V. Cherepanov, Y. Erdogan, G. Flügge, H. Geenen, M. Geisler, F. Hoehle, B. Kargoll, T. Kress, A. Künsken, J. Lingemann, A. Nehrkorn, A. Nowack, I.M. Nugent, C. Pistone, O. Pooth, A. Stahl Deutsches Elektronen-Synchrotron, Hamburg, Germany
M. Aldaya Martin, I. Asin, N. Bartosik, O. Behnke, U. Behrens, K. Borras\@textsuperscript16, A. Burgmeier, A. Campbell, C. Contreras-Campana, F. Costanza, C. Diez Pardos, G. Dolinska, S. Dooling, T. Dorland, G. Eckerlin, D. Eckstein, T. Eichhorn, G. Flucke, E. Gallo\@textsuperscript17, J. Garay Garcia, A. Geiser, A. Gizhko, P. Gunnellini, J. Hauk, M. Hempel\@textsuperscript18, H. Jung, A. Kalogeropoulos, O. Karacheban\@textsuperscript18, M. Kasemann, P. Katsas, J. Kieseler, C. Kleinwort, I. Korol, W. Lange, J. Leonard, K. Lipka, A. Lobanov, W. Lohmann\@textsuperscript18, R. Mankel, I.-A. Melzer-Pellmann, A.B. Meyer, G. Mittag, J. Mnich, A. Mussgiller, S. Naumann-Emme, A. Nayak, E. Ntomari, H. Perrey, D. Pitzl, R. Placakyte, A. Raspereza, B. Roland, M.Ö. Sahin, P. Saxena, T. Schoerner-Sadenius, C. Seitz, S. Spannagel, N. Stefaniuk, K.D. Trippkewitz, R. Walsh, C. Wissing University of Hamburg, Hamburg, Germany
V. Blobel, M. Centis Vignali, A.R. Draeger, J. Erfle, E. Garutti, K. Goebel, D. Gonzalez, M. Görner, J. Haller, M. Hoffmann, R.S. Höing, A. Junkes, R. Klanner, R. Kogler, N. Kovalchuk, T. Lapsien, T. Lenz, I. Marchesini, D. Marconi, M. Meyer, D. Nowatschin, J. Ott, F. Pantaleo\@textsuperscript2, T. Peiffer, A. Perieanu, N. Pietsch, J. Poehlsen, D. Rathjens, C. Sander, C. Scharf, P. Schleper, E. Schlieckau, A. Schmidt, S. Schumann, J. Schwandt, V. Sola, H. Stadie, G. Steinbrück, F.M. Stober, H. Tholen, D. Troendle, E. Usai, L. Vanelderen, A. Vanhoefer, B. Vormwald Institut für Experimentelle Kernphysik, Karlsruhe, Germany
C. Barth, C. Baus, J. Berger, C. Böser, E. Butz, T. Chwalek, F. Colombo, W. De Boer, A. Descroix, A. Dierlamm, S. Fink, F. Frensch, R. Friese, M. Giffels, A. Gilbert, D. Haitz, F. Hartmann\@textsuperscript2, S.M. Heindl, U. Husemann, I. Katkov\@textsuperscript5, A. Kornmayer\@textsuperscript2, P. Lobelle Pardo, B. Maier, H. Mildner, M.U. Mozer, T. Müller, Th. Müller, M. Plagge, G. Quast, K. Rabbertz, S. Röcker, F. Roscher, M. Schröder, G. Sieber, H.J. Simonis, R. Ulrich, J. Wagner-Kuhr, S. Wayand, M. Weber, T. Weiler, S. Williamson, C. Wöhrmann, R. Wolf Institute of Nuclear and Particle Physics (INPP),  NCSR Demokritos, Aghia Paraskevi, Greece
G. Anagnostou, G. Daskalakis, T. Geralis, V.A. Giakoumopoulou, A. Kyriakis, D. Loukas, A. Psallidas, I. Topsis-Giotis National and Kapodistrian University of Athens, Athens, Greece
A. Agapitos, S. Kesisoglou, A. Panagiotou, N. Saoulidou, E. Tziaferi University of Ioánnina, Ioánnina, Greece
I. Evangelou, G. Flouris, C. Foudas, P. Kokkas, N. Loukas, N. Manthos, I. Papadopoulos, E. Paradas, J. Strologas Wigner Research Centre for Physics, Budapest, Hungary
G. Bencze, C. Hajdu, A. Hazi, P. Hidas, D. Horvath\@textsuperscript19, F. Sikler, V. Veszpremi, G. Vesztergombi\@textsuperscript20, A.J. Zsigmond Institute of Nuclear Research ATOMKI, Debrecen, Hungary
N. Beni, S. Czellar, J. Karancsi\@textsuperscript21, J. Molnar, Z. Szillasi\@textsuperscript2 University of Debrecen, Debrecen, Hungary
M. Bartók\@textsuperscript22, A. Makovec, P. Raics, Z.L. Trocsanyi, B. Ujvari National Institute of Science Education and Research, Bhubaneswar, India
S. Choudhury\@textsuperscript23, P. Mal, K. Mandal, D.K. Sahoo, N. Sahoo, S.K. Swain Panjab University, Chandigarh, India
S. Bansal, S.B. Beri, V. Bhatnagar, R. Chawla, R. Gupta, U.Bhawandeep, A.K. Kalsi, A. Kaur, M. Kaur, R. Kumar, A. Mehta, M. Mittal, J.B. Singh, G. Walia University of Delhi, Delhi, India
Ashok Kumar, A. Bhardwaj, B.C. Choudhary, R.B. Garg, S. Malhotra, M. Naimuddin, N. Nishu, K. Ranjan, R. Sharma, V. Sharma Saha Institute of Nuclear Physics, Kolkata, India
S. Bhattacharya, K. Chatterjee, S. Dey, S. Dutta, N. Majumdar, A. Modak, K. Mondal, S. Mukhopadhyay, A. Roy, D. Roy, S. Roy Chowdhury, S. Sarkar, M. Sharan Bhabha Atomic Research Centre, Mumbai, India
R. Chudasama, D. Dutta, V. Jha, V. Kumar, A.K. Mohanty\@textsuperscript2, L.M. Pant, P. Shukla, A. Topkar Tata Institute of Fundamental Research, Mumbai, India
T. Aziz, S. Banerjee, S. Bhowmik\@textsuperscript24, R.M. Chatterjee, R.K. Dewanjee, S. Dugad, S. Ganguly, S. Ghosh, M. Guchait, A. Gurtu\@textsuperscript25, Sa. Jain, G. Kole, S. Kumar, B. Mahakud, M. Maity\@textsuperscript24, G. Majumder, K. Mazumdar, S. Mitra, G.B. Mohanty, B. Parida, T. Sarkar\@textsuperscript24, N. Sur, B. Sutar, N. Wickramage\@textsuperscript26 Indian Institute of Science Education and Research (IISER),  Pune, India
S. Chauhan, S. Dube, A. Kapoor, K. Kothekar, S. Sharma Institute for Research in Fundamental Sciences (IPM),  Tehran, Iran
H. Bakhshiansohi, H. Behnamian, S.M. Etesami\@textsuperscript27, A. Fahim\@textsuperscript28, M. Khakzad, M. Mohammadi Najafabadi, M. Naseri, S. Paktinat Mehdiabadi, F. Rezaei Hosseinabadi, B. Safarzadeh\@textsuperscript29, M. Zeinali University College Dublin, Dublin, Ireland
M. Felcini, M. Grunewald INFN Sezione di Bari , Università di Bari , Politecnico di Bari ,  Bari, Italy
M. Abbrescia, C. Calabria, C. Caputo, A. Colaleo, D. Creanza, L. Cristella, N. De Filippis, M. De Palma, L. Fiore, G. Iaselli, G. Maggi, M. Maggi, G. Miniello, S. My, S. Nuzzo, A. Pompili, G. Pugliese, R. Radogna, A. Ranieri, G. Selvaggi, L. Silvestris\@textsuperscript2, R. Venditti INFN Sezione di Bologna , Università di Bologna ,  Bologna, Italy
G. Abbiendi, C. Battilana\@textsuperscript2, D. Bonacorsi, S. Braibant-Giacomelli, L. Brigliadori, R. Campanini, P. Capiluppi, A. Castro, F.R. Cavallo, S.S. Chhibra, G. Codispoti, M. Cuffiani, G.M. Dallavalle, F. Fabbri, A. Fanfani, D. Fasanella, P. Giacomelli, C. Grandi, L. Guiducci, S. Marcellini, G. Masetti, A. Montanari, F.L. Navarria, A. Perrotta, A.M. Rossi, T. Rovelli, G.P. Siroli, N. Tosi\@textsuperscript2 INFN Sezione di Catania , Università di Catania ,  Catania, Italy
G. Cappello, M. Chiorboli, S. Costa, A. Di Mattia, F. Giordano, R. Potenza, A. Tricomi, C. Tuve INFN Sezione di Firenze , Università di Firenze ,  Firenze, Italy
G. Barbagli, V. Ciulli, C. Civinini, R. D’Alessandro, E. Focardi, V. Gori, P. Lenzi, M. Meschini, S. Paoletti, G. Sguazzoni, L. Viliani\@textsuperscript2 INFN Laboratori Nazionali di Frascati, Frascati, Italy
L. Benussi, S. Bianco, F. Fabbri, D. Piccolo, F. Primavera\@textsuperscript2 INFN Sezione di Genova , Università di Genova ,  Genova, Italy
V. Calvelli, F. Ferro, M. Lo Vetere, M.R. Monge, E. Robutti, S. Tosi INFN Sezione di Milano-Bicocca , Università di Milano-Bicocca ,  Milano, Italy
L. Brianza, M.E. Dinardo, S. Fiorendi, S. Gennai, R. Gerosa, A. Ghezzi, P. Govoni, S. Malvezzi, R.A. Manzoni\@textsuperscript2, B. Marzocchi, D. Menasce, L. Moroni, M. Paganoni, D. Pedrini, S. Ragazzi, N. Redaelli, T. Tabarelli de Fatis INFN Sezione di Napoli , Università di Napoli ’Federico II’ , Napoli, Italy, Università della Basilicata , Potenza, Italy, Università G. Marconi , Roma, Italy
S. Buontempo, N. Cavallo, S. Di Guida\@textsuperscript2, M. Esposito, F. Fabozzi, A.O.M. Iorio, G. Lanza, L. Lista, S. Meola\@textsuperscript2, M. Merola, P. Paolucci\@textsuperscript2, C. Sciacca, F. Thyssen INFN Sezione di Padova , Università di Padova , Padova, Italy, Università di Trento , Trento, Italy
P. Azzi\@textsuperscript2, N. Bacchetta, L. Benato, D. Bisello, A. Boletti, A. Branca, R. Carlin, P. Checchia, M. Dall’Osso\@textsuperscript2, T. Dorigo, U. Dosselli, F. Gasparini, U. Gasparini, A. Gozzelino, K. Kanishchev, S. Lacaprara, M. Margoni, A.T. Meneguzzo, F. Montecassiano, J. Pazzini\@textsuperscript2, N. Pozzobon, P. Ronchese, F. Simonetto, E. Torassa, M. Tosi, M. Zanetti, P. Zotto, A. Zucchetta\@textsuperscript2, G. Zumerle INFN Sezione di Pavia , Università di Pavia ,  Pavia, Italy
A. Braghieri, A. Magnani, P. Montagna, S.P. Ratti, V. Re, C. Riccardi, P. Salvini, I. Vai, P. Vitulo INFN Sezione di Perugia , Università di Perugia ,  Perugia, Italy
L. Alunni Solestizi, G.M. Bilei, D. Ciangottini\@textsuperscript2, L. Fanò, P. Lariccia, G. Mantovani, M. Menichelli, A. Saha, A. Santocchia INFN Sezione di Pisa , Università di Pisa , Scuola Normale Superiore di Pisa ,  Pisa, Italy
K. Androsov\@textsuperscript30, P. Azzurri\@textsuperscript2, G. Bagliesi, J. Bernardini, T. Boccali, R. Castaldi, M.A. Ciocci