alphaFeedback

WorksPhysics and Astronomy

Search for heavy neutral leptons in top quark decays with same-sign, same-flavor dilepton final states in proton-proton collisions at s\sqrt{s} = 13 TeV

AI contributions
Not recorded for this paper
View details

Abstract

A search for a Majorana-type heavy neutral lepton (HNL) produced in top quark decays is reported. The search is based on proton-proton collision data at s\sqrt{s} = 13 TeV, corresponding to an integrated luminosity of 138 fb−1^{-1}, collected during 2016−-2018 with the CMS detector at the CERN LHC. In the type-I seesaw model, a Majorana HNL can mediate lepton-number-violating top quark decays. The search targets top quark pair production events in which one top quark decays into a bottom quark and a W boson, which subsequently decays into a charged lepton (electron or muon) and an HNL. The Majorana HNL then decays into a W boson and a charged lepton with the same sign and flavor as the one produced with the HNL. The experimental signature features a same-sign, same-flavor lepton pair and jets. The HNL masses between 20 and 100 GeV are probed. The data agree with the standard model background prediction. Observed (expected) upper limits at 95% confidence level are set on the branching fraction of the top quark decay into an HNL, a charged lepton, and a bottom quark, varying between [1.1 (0.71)] ×\times10−5^{-5} and [4.8 (3.0)] ×\times10−5^{-5} in the dielectron channel, and between [3.7 (3.4)] ×\times10−6^{-6} and [2.4 (1.3)] ×\times10−5^{-5} in the dimuon channel.

Introduction

The observation of neutrino oscillations has confirmed the nonzero masses for at least two mass eigenstates of the standard model (SM) neutrinos (νℓ\nu_\ell) [1]. The mass-squared differences measured from neutrino oscillation experiments have set lower limits on the mass values of the two heavier neutrinos at 8.7 and 50 meV [1]. Cosmological measurements of the cosmic microwave background and baryon acoustic oscillations also place an upper limit on the sum of the neutrino masses at 64 meV [2]. An important question is the origin of the smallness of the neutrino masses. Since neutrinos are assumed to be massless in the SM, several beyond-the-SM (BSM) scenarios have been suggested to resolve this issue. For example, in the type-I seesaw mechanism [3, 4, 5, 6], Majorana mass (mMm_{\mathrm{M}}) terms are introduced for right-handed chiral neutrino fields, in addition to Dirac masses (mDm_{\mathrm{D}}) originating from the Yukawa couplings of neutrino fields to the Higgs field. This results in two physical states for each neutrino flavor: a light SM-like neutrino νℓ\nu_\ell and a heavy neutral lepton (HNL) NN. The mixing between νℓ\nu_\ell and NN can be quantified by a mixing parameter VℓNV_{\ell N}, which—along with the physical masses—is determined by mDm_{\mathrm{D}} and mMm_{\mathrm{M}}. This mixing parameter determines the HNL component of the neutrino flavor eigenstates and hence determines the interaction between the HNL and the W boson. In the simplest seesaw realization, the light-neutrino masses scale approximately as mD2/mMm_{\mathrm{D}}^2/m_{\mathrm{M}}, where mMm_{\mathrm{M}} is effectively the HNL mass (mNm_N). More general realizations can allow larger HNL-SM mixing for a given mNm_N [7].

Searches for HNL span a broad mNm_N range. At low masses of O(100)\mathcal{O}(100) eV, beta-decay spectra have been examined for spectral kinks that would indicate the presence of an HNL [8, 9]. At collider-accessible mNm_N values up to O(1)\mathcal{O}(1) TeV, the CERN LHC experiments searched for the decay of the on- and off-shell W bosons, directly produced in proton-proton (pp) collisions, into an HNL (pp→W±→Nℓ±\mathrm{pp} \to\mathrm{W}^{\pm} \to N\ell^{\pm}), where the hadronic (N→ℓ±qqˉ′N \to\ell^{\pm}q\bar{q}') and leptonic (N→ℓ±ℓ′∓νN \to\ell^{\pm}\ell'^{\mp}\nu) decay chains of an HNL into quarks (qq), charged leptons (ℓ\ell) including the tau leptons (τ\tau), and νℓ\nu_\ell lead to dilepton and trilepton final states, respectively [10, 11, 12, 13]. The LHCb experiment searched for HNLs produced in B meson decays and decaying into a muon (μ\mu) and a pion [14]. Throughout the paper, ℓ\ell indicates electrons (e) and muons, excluding τ\tau leptons unless mentioned otherwise. These direct searches for HNLs have focused on analyzing their decay products or the energy spectrum of the charged lepton produced with the HNL.

The hypothetical production of HNLs via top quark (t) decays has recently received attention in the theory community and provides a complementary probe of HNLs below the W boson mass (mWm_{\mathrm{W}}) [15, 16, 17, 18, 19]. In some BSM scenarios, such as supersymmetric models with RR-parity violation or models with leptoquarks, destructive interference can reduce the sensitivity of direct-search channels, increasing the complementarity of HNL production via top quark decays. In addition, the presence of both W+\mathrm{W}^{+} and W−\mathrm{W}^{-} bosons in top quark pair (ttˉ\mathrm{t}\bar{\mathrm{t}}) events makes this production mode potentially sensitive to charge-parity violating effects in the HNL sector [19]. The ATLAS Collaboration recently searched for Majorana HNLs produced in top quark decays in the same-sign (SS), same-flavor (SF) dilepton final state, using 140 fb−1^{-1} of data at s=13\sqrt{s} = 13 TeV, and set limits for HNL masses between 15 and 75 GeV [20]. In this article, we present a search for HNLs in top quark decays using SS, SF dilepton events in the pp collision data at s=13\sqrt{s} = 13 TeV, collected using the CMS detector between 2016 and 2018, corresponding to an integrated luminosity of 138 fb−1^{-1}.2

The search targets lepton number violating (LNV) processes in ttˉ\mathrm{t}\bar{\mathrm{t}} production, where the top quark decays into an HNL, which further decays into SM particles, t→bW+→bℓ+N→bℓ+ℓ+W−→bℓ+ℓ+ffˉ′\mathrm{t} \to\mathrm{bW}^{+} \to\mathrm{b}\ell^{+}N \to\mathrm{b}\ell^{+}\ell^{+}\mathrm{W}^{-} \to\mathrm{b}\ell^{+}\ell^{+}\mathrm{f}\bar{\mathrm{f}}', while the top antiquark decays via the SM process tˉ→bWˉ−→bˉffˉ′\bar{\mathrm{t}} \to\bar{\mathrm{bW}}^{-} \to\bar{\mathrm{b}}\mathrm{f}\bar{\mathrm{f}}', where f denotes an SM fermion. Charge-conjugated processes are implied throughout the text.

If mN<mWm_{\mathrm{N}} < m_{\mathrm{W}}, the W boson produced in the decay of the HNL is off-shell; if mN>mWm_{\mathrm{N}} > m_{\mathrm{W}}, the W boson from which the HNL originates is off-shell. Representative Feynman diagrams of these processes are depicted in Figure 1. In this search, only mixings of an HNL with the electron and muon neutrinos are considered, with strengths ∣VeN∣2|V_{\mathrm{eN}}|^2 and ∣VμN∣2|V_{\mu\mathrm{N}}|^2, respectively. Each mixing parameter is considered separately, with one mixing parameter nonzero and the other set to zero.

Representative Feynman diagrams of the SM top quark decays and BSM decays involving an HNL

Figure 1. Representative Feynman diagrams of the SM top quark decays (left) and BSM decays involving an HNL (right), where ℓ\ell represents e or μ\mu, and f, f′ represent SM fermions.

Searching for LNV decay modes is of paramount interest, as they provide a characteristic signature of Majorana HNLs. The discovery of a particle with a Majorana mass would provide evidence for the existence of different mass generation mechanisms in nature. In the context of an LNV search, the hadronic decay mode N→ℓqq‾′\mathrm{N} \to\ell\mathrm{q}\overline{\mathrm{q}}' is the most powerful one, as there is no ambiguity between the lepton number conserving and violating decays. In the case of the leptonic decay mode, N→ℓℓ′ν\mathrm{N} \to\ell\ell'\nu, it is impossible to conclusively determine whether LNV occurred, since the lepton number of the undetected neutrino through a decay involving a Z boson, N→Zν→ℓℓν\mathrm{N} \to\mathrm{Z}\nu\to\ell\ell\nu, cannot be known. The possibility of an HNL mixing with more than one generation of SM neutrinos also complicates the interpretation of events in the fully leptonic decay modes, even for the W boson mediated cases.

In this search, mNm_{\mathrm{N}} values between 20 and 100 GeV are considered. For masses below mWm_{\mathrm{W}}, the HNL decay width scales approximately as mN5m_{\mathrm{N}}^5 for a fixed mixing parameter, so a smaller mNm_{\mathrm{N}} results in a longer lifetime of the HNL [16]. For mNm_{\mathrm{N}} values below 20 GeV, the decay length of an HNL cannot be neglected within the allowed mixing parameter region and must be addressed using specialized search techniques optimized for long-lived particles. For masses above mWm_{\mathrm{W}}, this HNL production channel cannot be mediated by an on-shell W boson from a top quark decay. The production rate decreases rapidly with increasing mNm_{\mathrm{N}} values, and this analysis therefore restricts the search to mN<100m_{\mathrm{N}} < 100 GeV. A machine-learning technique is used to classify events as background or signal, based on their distinguishing characteristics. The backgrounds are composed of nonprompt-lepton, conversion, charge mismeasurement, and prompt-lepton components, as discussed in § 5.

This paper is structured as follows. The CMS detector and object reconstruction are described in § 2, and the data and simulation samples are discussed in § 3. Event reconstruction and selection are presented in § 4. The background estimation techniques are summarized in § 5. The systematic uncertainties considered in this analysis are described in § 6. Lastly, the results and interpretations are presented in § 8, with an overall summary in § 9. Tabulated results are available in the HEPData record [21] for this analysis.

The CMS detector and event reconstruction

The CMS apparatus [22, 23] is a multipurpose, nearly hermetic detector, designed to trigger on [24, 25, 26] and identify e, μ\mu, photons (γ\gamma), and (charged and neutral) hadrons [27, 28, 29]. Its central feature 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. Forward calorimeters extend the pseudorapidity (η\eta) coverage provided by the barrel and endcap detectors. Muons are reconstructed using gas-ionization detectors interleaved with the layers of the steel flux-return yoke outside the solenoid. More detailed descriptions of the CMS detector, together with a definition of the coordinate system used and the relevant kinematic variables, can be found in Refs. [22, 23].

Events of interest are selected using a two-tiered trigger system. The first level (L1), composed of custom hardware processors, uses information from the calorimeters and muon detectors to select events at a rate of around 100 kHz within a fixed latency of 4 μs4\,\mu\mathrm{s} [24]. The second level, known as the high-level trigger, consists of a farm of processors running a version of the full event reconstruction software optimized for fast processing, and reduces the event rate to a few kHz before data storage [25, 26].

The global event description with the particle-flow (PF) algorithm [30] aims to reconstruct and identify each individual particle in an event, with an optimized combination of all subdetector information. In this process, the identification of the particle type (photon, electron, muon, charged hadron, neutral hadron) plays an important role in the determination of the particle direction and energy. Photons are identified as ECAL energy clusters not linked to the extrapolation of any charged-particle trajectory to the ECAL. Electrons are identified as a primary charged-particle track and potentially many ECAL energy clusters corresponding to this track extrapolation to the ECAL and to possible bremsstrahlung photons emitted along the way through the tracker material. Muons are identified as tracks in the central tracker consistent with either a track or several hits in the muon system, and associated with calorimeter deposits compatible with the muon hypothesis. Charged hadrons are identified as charged-particle tracks neither identified as electrons, nor as muons. Finally, neutral hadrons are identified as HCAL energy clusters not linked to any charged-hadron trajectory, or as a combined ECAL and HCAL energy excess with respect to the expected charged-hadron energy deposit. The primary vertex (PV) of the interaction is taken to be the vertex corresponding to the hardest scattering in the event, evaluated using tracking information alone, as described in Section 9.4.1 of Ref. [31].

Electrons reconstructed within the tracker coverage (∣η∣<2.5|\eta| < 2.5), except for the calorimeter gap region (1.442<∣η∣<1.5661.442 < |\eta| < 1.566), are used in the analysis. We consider electrons with transverse momentum pT>15 GeVp_{\mathrm{T}} > 15\,\mathrm{GeV} and require them to be isolated from other particles in the event to reject nonprompt leptons from decays in jets. Two electron identification (ID) working points are used in the analysis, referred to as loose and tight ID in the following. Electron ID is based on track quality criteria, matching between the electron trajectory and the associated ECAL cluster, and shower shape observables [27]. The relative isolation (IrelI_{\mathrm{rel}}) is defined as the scalar pTp_{\mathrm{T}} sum of all PF particles reconstructed within a cone around the lepton direction divided by the lepton pTp_{\mathrm{T}}. The cone radius is defined by ΔR=(Δη)2+(Δϕ)2\Delta R = \sqrt{(\Delta\eta)^2 + (\Delta\phi)^2}, where ϕ\phi is the azimuthal angle in radians. We use a pTp_{\mathrm{T}}-dependent cone radius of 0.2 for leptons with pT<50 GeVp_{\mathrm{T}} < 50\,\mathrm{GeV}, of 10 GeV/pT10\,\mathrm{GeV}/p_{\mathrm{T}} for 50<pT<200 GeV50 < p_{\mathrm{T}} < 200\,\mathrm{GeV}, and of 0.05 for pT>200 GeVp_{\mathrm{T}} > 200\,\mathrm{GeV}, which improves the efficiency for high-pTp_{\mathrm{T}} leptons by removing the accidental overlap with other particles [32]. The loose and tight electron IDs require Irel<0.4I_{\mathrm{rel}} < 0.4 and 0.1, respectively. Furthermore, we require that there be at most one tracker layer that contributes no hit on the track, to reduce contributions from photon conversions [27]. To pass the tight ID selection, an electron must also be compatible with originating from the PV, requiring the ratio of the three-dimensional impact parameter of the associated track to its uncertainty (SIP3D_{3D}) to be less than 4. The loose electron ID is used to veto additional leptons and to estimate the nonprompt-lepton background. The tight electron ID is used for the signal region (SR) and background validation region (VR) definitions. The tight electron ID efficiency for the signal event selection varies between 40 and 90%, depending on the data-taking era, and on the electron pTp_\mathrm{T} and η\eta. Because the electron charge is crucial and has a higher mismeasurement rate than that of muons, we use a tight charge ID criterion requiring agreement in the charge measurements from three different methods. These methods are based on the Gaussian sum filter (GSF) track curvature [33], the closest Kalman filter track curvature [30], and the Δϕ\Delta\phi between the ECAL cluster and the GSF track extrapolated to the vertex. An event is discarded if it contains an electron for which the three charge measurements do not agree.

Muons are reconstructed in the range ∣η∣<2.4|\eta| < 2.4, and we consider muons with pT>15p_\mathrm{T} > 15 GeV. Two muon ID working points are defined in the analysis, referred to as loose and tight ID in the following, serving the same purposes as the electron ID working points. Muon ID is based on the number of reconstructed energy deposits in the tracker and in the muon system, and on the fit quality of the muon track [28]. Additionally, the loose (tight) muon ID requires Irel<0.6I_\mathrm{rel} < 0.6 (0.1) and SIP3D<5\mathrm{SIP}_{3\mathrm{D}} < 5 (3). The muon tight ID efficiency for the signal event selection varies between 70 and 98% depending on the data-taking era, and on the muon pTp_\mathrm{T} and η\eta.

For each event, hadronic jets are clustered from the PF particles using the infrared- and collinear-safe anti-kTk_\mathrm{T} algorithm [34, 35] with a distance parameter of 0.4. Jet momentum is determined as the vectorial sum of all particle momenta in the jet, and is found from simulation to be, on average, within 5–10% of the true momentum over the entire pTp_\mathrm{T} spectrum and detector acceptance. The jet energy resolution amounts typically to 15–20% at 30 GeV, 10% at 100 GeV, and 5% at 1 TeV [36]. Additional pp interactions within the same or nearby bunch crossings (pileup) can contribute extra tracks and calorimetric energy depositions, increasing the apparent jet momentum. To mitigate this effect, tracks identified to be originating from pileup vertices are discarded, and an offset correction is applied to correct for the remaining pileup contributions [37]. Jet energy corrections are derived from simulation studies so that the average measured energy of jets becomes identical to that of particle-level jets. In situ measurements of the momentum balance in dijet, γ\gamma + jet, Z + jet, and multijet events are used to determine any residual differences between the jet energy scale in data and simulation, and appropriate corrections are made [36]. Jets reconstructed with pTp_\mathrm{T} larger than 25 GeV and ∣η∣|\eta| smaller than 2.4 are considered in the analysis. Additional selection criteria are applied to each jet to remove jets potentially dominated by instrumental effects or reconstruction failures [37]. Lastly, for jets with pT<50p_\mathrm{T} < 50 GeV, an additional identification algorithm based on a boosted decision tree technique, referred to as the pileup jet ID, is used to remove jets likely originating from pileup interactions [38].

To identify jets originating from the hadronization of b quarks, the DEEPJET [39, 40] b tagging algorithm is applied to each reconstructed jet in an event. The threshold on the discriminant score of the DEEPJET algorithm for a jet to be considered a b-tagged jet is chosen such that the tagging efficiency for b quark initiated jets is approximately 80% with a misidentification probability of 1% for jets initiated by light-flavor quarks or gluons. The tagging efficiencies measured in simulation are corrected as functions of the pTp_\mathrm{T}, η\eta, and flavor of the jet to match those measured in data [41].

The missing transverse momentum vector p⃗T miss\vec{p}_\mathrm{T}^{\,\mathrm{miss}} is computed as the negative vector pTp_\mathrm{T} sum of all the PF candidates in an event, and its magnitude is denoted as pTmissp_\mathrm{T}^{\mathrm{miss}} [42]. The p⃗T miss\vec{p}_\mathrm{T}^{\,\mathrm{miss}} is modified to account for corrections to the energy scale of the reconstructed jets, muons, and electrons in the event.

Data and simulated samples

The analysis is performed using a sample of pp collision data at s=13\sqrt{s}=13 TeV, collected with the CMS detector in 2016–2018 and corresponding to a total integrated luminosity of 138 fb−1^{-1} [43, 44, 45]. Events used in this analysis are selected by SF dilepton triggers. Each trigger has pTp_{\mathrm{T}} threshold requirements for the leading (pT(ℓ1)p_{\mathrm{T}}(\ell_{1})) and subleading (pT(ℓ2)p_{\mathrm{T}}(\ell_{2})) leptons. For each lepton in the trigger, a track-based IrelI_{\mathrm{rel}} threshold requirement is implemented with a fixed cone radius of 0.2 for electrons and 0.3 for muons. Depending on the data-taking era, the trigger also requires the two leptons to originate from the PV (via a requirement on the distance of closest approach, dzd_{z}, along the beam axis) and/or to pass a minimum invariant mass (m(ℓℓ)m(\ell\ell)) threshold. Additionally, two single-lepton triggers for each lepton generation are used for the nonprompt-lepton background estimation. The thresholds are summarized in Table 1.

Trigger typepT(ℓ1)p_{\mathrm{T}}(\ell_{1}) (GeV)pT(ℓ2)p_{\mathrm{T}}(\ell_{2}) (GeV)Other
Dilepton triggers
Dielectron>23>23>12>12Irel<0.20I_{\mathrm{rel}}<0.20; dz<0.2d_{z}<0.2 cm (2016)
Dimuon>17>17>8>8Irel<0.40I_{\mathrm{rel}}<0.40; m(μμ)>3.8m(\mu\mu)>3.8 GeV (2017, 2018); dz<0.2d_{z}<0.2 cm
Single-lepton triggers
Single electron>12>12, pT(e)<35p_{\mathrm{T}}(e)<35 GeV; >23>23, pT(e)>35p_{\mathrm{T}}(e)>35 GeV—Irel<0.20I_{\mathrm{rel}}<0.20, pT(j1)>30p_{\mathrm{T}}(j_{1})>30 GeV
Single muon>8>8, pT(μ)<30p_{\mathrm{T}}(\mu)<30 GeV; >17>17, pT(μ)>30p_{\mathrm{T}}(\mu)>30 GeV—Irel<0.40I_{\mathrm{rel}}<0.40

Table 1. Summary of trigger pTp_{\mathrm{T}} thresholds and additional selection requirements. The dilepton triggers are used to define the SR and VR. For the nonprompt-lepton background estimate, two prescaled single-lepton triggers are used per flavor, depending on the lepton pTp_{\mathrm{T}}.

The signal events are simulated with MADGRAPH5_aMC@NLO [46] version 2.6.5 using the Majorana HNL model [47, 48, 49] at leading order (LO) in the four-flavor scheme. The tt‾\mathrm{t}\overline{\mathrm{t}} cross section is calculated with the TOP++ version 2.0 program to next-to-next-to-leading order (NNLO) in perturbative quantum chromodynamics (QCD), including soft-gluon resummation to next-to-next-to-leading log (NNLL) order (as discussed in Ref. [50] and references therein), and assuming a top quark mass of 172.5 GeV. The signal events from tt‾\mathrm{t}\overline{\mathrm{t}} production with subsequent top quark decays into an HNL are generated with up to one additional parton in the matrix element calculations. Subsequently, jets from parton shower are merged with those from the matrix element calculations using the MLM algorithm [51]. Decays of the top quarks, W bosons, and HNLs are modeled with the same generator, considering the spin correlations among the decay products. The signal samples are simulated for mNm_{N} values of 20, 30, 40, 50, 60, 70, 75, 85, 95, and 100 GeV. The decay width of the HNL in the considered mass range is less than 0.1 MeV. The HNL is assumed to decay promptly. The nonzero HNL mixing parameter ∣VℓN∣|V_{\ell N}| is set to 0.01. In addition to the decay mode introduced in § 1, the generated signal samples include all HNL decay modes leading to final states with more than two charged leptons. We exclude τ\tau leptons from the primary HNL vertices (W→τNW\to\tau N and N→τWN\to\tau W), while allowing for τ\tau lepton production via W and Z boson decays. The full decay modes considered are:

pp→tt‾→(t→bℓ+N→bℓ+ℓ+ff‾′, t‾→b‾W−→b‾ff‾′),\mathrm{pp}\to\mathrm{t}\overline{\mathrm{t}}\to\left(\mathrm{t}\to\mathrm{b}\ell^{+}\mathrm{N}\to\mathrm{b}\ell^{+}\ell^{+}f\overline{f}',\ \overline{\mathrm{t}}\to\overline{\mathrm{b}}\mathrm{W}^{-}\to\overline{\mathrm{b}}f\overline{f}'\right),

where ff and f′f' can be any fermion and NN decay modes include both charged- and neutral-current. The inclusive signal cross section is set to 10 fb, and the fitted signal strength limit directly scales this cross section. The branching fraction limit of a top quark into an HNL is then obtained by dividing the fitted signal cross section limit by 2σ(ttˉ)2\sigma(\mathrm{t\bar{t}}). This branching fraction limit is then compared to the branching fraction calculated in simulation with ∣VℓN∣=0.01|V_{\ell N}| = 0.01 and converted into the mixing parameter limit using the relation B(t→bℓN)∝∣VℓN∣2\mathcal{B}(\mathrm{t} \to\mathrm{b}\ell\mathrm{N}) \propto|V_{\ell N}|^2.

Background processes involving leptons produced promptly and from photon conversions are estimated using simulation. Additionally, simulation is used to model the background from lepton charge mismeasurement, causing opposite-sign (OS) lepton pairs to be reconstructed as SS lepton pairs. These processes are generated at next-to-leading order (NLO), and their cross sections are calculated at the same order, unless a specific other order is mentioned. The following processes are simulated using MADGRAPH5_aMC@NLO: Drell–Yan (DY) production, WZ diboson production with decays into three charged leptons and a neutrino, triboson production with all possible combinations of vector bosons (VVV, V = W or Z), W boson production in association with a photon (Wγ\gamma), Higgs boson production in association with a vector boson (VH), vector boson production in association with a ttˉ\mathrm{t\bar{t}} pair (ttˉ\mathrm{t\bar{t}}V), single top quark production via the ss channel, and single top quark production in association with a Z boson (tZ). The DY processes with m(ℓℓ)m(\ell\ell) below 50 GeV and Wγ\gamma production are generated at LO. The decays in VH and ttˉ\mathrm{t\bar{t}}W processes are simulated using MADSPIN [52]. The following processes are generated using POWHEG v2.0 [53, 54, 55]: production of a pair of Z bosons (ZZ) decaying into four leptons [56], Higgs boson production in association with a ttˉ\mathrm{t\bar{t}} pair (ttˉ\mathrm{t\bar{t}}H) [57], Higgs boson production via gluon fusion (ggH) [58] and via vector boson fusion (qqH) [59], single top quark production in the tt channel and in association with a W boson [60], and ttˉ\mathrm{t\bar{t}} production [61]. For the normalization of the ZZ process, the production cross section is calculated at NNLO [62]. The decays of the Higgs boson in ggH and qqH production are simulated with JHUGEN [63], while MADSPIN is used for the top quark decay in the single top quark tt-channel production. Lastly, the diboson processes, except for the Z boson pair production with decays into four leptons and WZ diboson production with decays into three charged leptons and a neutrino, are generated at LO using PYTHIA 8.240 [64].

All simulated samples are interfaced with PYTHIA to perform parton shower and hadronization. For all simulated samples, the description of the underlying event is provided by the CP5 tune [65]. The NNPDF 3.1 NNLO [66] parton distribution function (PDF) set is used for all simulated samples. Pileup effects are included via simulated minimum bias pp collisions superimposed on the primary simulated interaction, with the multiplicity matching that in data. Simulation of the CMS detector response is performed with GEANT4 [67].

Event selection and analysis strategy

The offline selection requires an electron pair of pT>(25,15)p_{\mathrm{T}} > (25, 15) GeV or a muon pair of pT>(20,15)p_{\mathrm{T}} > (20, 15) GeV, where the pTp_{\mathrm{T}} thresholds are chosen such that the selected events lie at the trigger efficiency plateau. The baseline selection of the SR requires an event to have exactly two SF SS leptons to target the leptons from the production and decay of the HNL and vetoes the event if there is a third lepton that passes the loose lepton ID. To remove events involving the Z boson resonance with a charge-mismeasured electron, the reconstructed dielectron invariant mass (m(ee)m(ee)) must be outside a 10 GeV window around the Z boson mass (mZm_{\mathrm{Z}}) of 91.2 GeV [1]. However, this selection is not required for the dimuon channel because of the significantly lower charge mismeasurement rate for muons. The dimuon mass must be greater than 4 GeV, so that it is above the m(ℓℓ)m(\ell\ell) threshold of the trigger. The signal events are expected to have four jets from the W boson decays and two b-tagged jets from the top quark decays. However, the events are required to have at least three jets, which must include at least one b-tagged jet. Because signal jets tend to be more energetic and likely to pass the acceptance for mN>mWm_{\mathrm{N}} > m_{\mathrm{W}}, we define two baseline selections. When testing a signal hypothesis with mN<mWm_{\mathrm{N}} < m_{\mathrm{W}}, we require at least one non-b-tagged jet; when testing a hypothesis with mN>mWm_{\mathrm{N}} > m_{\mathrm{W}}, we require at least two. The selected data and background event sets for the two mass regimes can overlap. After the baseline selection, the signal acceptance times efficiency ranges from 0.987 (2.71) to 8.58 (7.54)%, depending on mNm_{\mathrm{N}}, in the electron (muon) channel. The dominant source of background events passing the SR selection comes from events with nonprompt leptons. This background is estimated by measuring the misidentification rate of a nonprompt lepton in the background measurement region (BMR) and applying it to the data in the application region (AR), as is further explained in § 5.

Because of the high multiplicity of partons in the signal process, only a small fraction of signal events has all the final-state partons within the pTp_{\mathrm{T}} and η\eta acceptances. For the simulated signal events passing the baseline selection, the efficiency ranges from 0.02 to 4% for matching all six partons at the generator level to the reconstructed jets within a ΔR\Delta R of 0.4, depending on mNm_{\mathrm{N}}. Therefore, instead of reconstructing the full signal decay products, we search for the signal by inspecting partial decay products using a multivariate analysis based on the gradient boosted decision tree (BDT) technique. Since the ttˉt\bar{t} background is the main contribution to the nonprompt-lepton background, only signal samples and dilepton and lepton+jets ttˉt\bar{t} background simulation samples are used for the training. As discussed in § 8, a training based only on ttˉt\bar{t} events still separates other backgrounds from the signal. Adding ttˉVt\bar{t}V and ttˉHt\bar{t}H backgrounds to the BDT training did not improve the expected sensitivity and led to instabilities in the nuisance parameter fit because of the small number of events in the SR. Despite the imbalance between the number of signal and background events, we used all ttˉt\bar{t} background and signal events, splitting each of them evenly: half for training and the other half for testing. Since the number of signal events passing the baseline selection varies for each mNm_{\mathrm{N}} and lepton flavor scenario, the total number of events used in the training process differs among the ten classifier models (for mN=20m_{\mathrm{N}} = 20, 50, 75, 85, and 100 GeV, each for dielectron and dimuon final states). Each model is applied to the signal sample with mNm_{\mathrm{N}} closest to the training signal sample mNm_{\mathrm{N}}.

The BDT classifier utilizes a set of input variables designed to exploit the kinematic differences between signal and background, specifically focusing on the energy scales, angular separations, and reconstructed masses associated with the final-state leptons and jets, as summarized in Table 2. For the reconstructed mass variables, jet candidates from the decays of a W boson (jWj_{\mathrm{W}}) and an HNL (jNj_{\mathrm{N}}) are determined using the combination of jets and leptons that best reconstructs the W boson and top quark masses. Applying the classifier model yields a score in the range [−1,1-1,1], where signal-like (background-like) events have scores closer to 1 (−1-1). Although the exact order of importance of the input variables differs across different mNm_{\mathrm{N}} and final states, the m(ℓℓ)m(\ell\ell) and angular separation of the two leptons are among the most important variables in all training models. These variables, along with the pTp_{\mathrm{T}} of the leading lepton and pTmissp_{\mathrm{T}}^{\mathrm{miss}} are shown in Figure 2–Figure 5. The backgrounds shown in these figures are explained in § 5. The performance of the training models is slightly better in the dimuon case, because of the larger number of events from higher muon ID efficiency and lower pTp_{\mathrm{T}} threshold. For both dielectron and dimuon final states, the performance slightly worsens as mNm_{\mathrm{N}} approaches mWm_{\mathrm{W}}, resulting in the area under the receiver operating characteristic curve decreasing from 0.90–0.93 to 0.82–0.88. As mNm_{\mathrm{N}} approaches mWm_{\mathrm{W}}, one of the two leptons becomes soft and is likely to fail the selection. As a result, the lepton that passes the selection would pair with a lepton from the SM top quark decay or with a nonprompt lepton, leading to background-like BDT outputs.

VariablemN<mWm_N < m_WmN>mWm_N > m_WDefinition
N(j)N(\mathrm{j})✓✓jet multiplicity
N(b)N(\mathrm{b})✓✓b-tagged jet multiplicity
pT(ℓa)p_T(\ell_a)✓✓atha^{\mathrm{th}} lepton pTp_T (a=1,2a = 1,2)
pTmissp_T^{\mathrm{miss}}✓✓missing transverse momentum
HTH_T✓✓scalar sum of jet pTp_T
ΔR(ℓ1,ℓ2)\Delta R(\ell_1,\ell_2)✓✓ΔR\Delta R between two leptons
ΔR(ℓa,jb)\Delta R(\ell_a,j_b)✓✓ΔR\Delta R between atha^{\mathrm{th}} lepton and bthb^{\mathrm{th}} jet (a=1,2a = 1,2, b=1,2,3b = 1,2,3)
m(ℓ1ℓ2)m(\ell_1\ell_2)✓✓dilepton invariant mass
m(ℓajN,1)m(\ell_a j_{N,1})✓✓invariant mass of atha^{\mathrm{th}} lepton and leading HNL jet candidate (a=1,2a = 1,2)
m(ℓ1ℓ2jN,1)m(\ell_1\ell_2j_{N,1})✓invariant mass of two leptons and leading HNL jet candidate
m(b1ℓ1ℓ2jN,1)m(b_1\ell_1\ell_2j_{N,1})✓✓invariant mass of two leptons, leading b-tagged jet, and leading HNL jet candidate
m(ℓajN,1jN,2)m(\ell_a j_{N,1}j_{N,2})✓invariant mass of atha^{\mathrm{th}} lepton and two HNL jet candidates (a=1,2a = 1,2)
m(b1ℓ1ℓ2jN,1jN,2)m(b_1\ell_1\ell_2j_{N,1}j_{N,2})✓invariant mass of two leptons, leading b-tagged jet, and two HNL jet candidates
mT(ℓ1p⃗Tmiss)m_T(\ell_1\vec{p}_T^{\mathrm{miss}})✓✓transverse mass of leading lepton and p⃗Tmiss\vec{p}_T^{\mathrm{miss}}
mT(b1ℓ1p⃗Tmiss)m_T(b_1\ell_1\vec{p}_T^{\mathrm{miss}})✓✓transverse mass of leading lepton leading b-tagged jet, and p⃗Tmiss\vec{p}_T^{\mathrm{miss}}
m(ℓ2jL,1)m(\ell_2j_{L,1})✓✓invariant mass of subleading lepton and leading non-b-tagged jet
m(ℓ2jW,1jW,2)m(\ell_2j_{W,1}j_{W,2})✓invariant mass of subleading lepton and two W boson jet candidates

Table 2. Input variables used for the training of classifiers. The usage in the low-mass (mN<mWm_N < m_W) and high-mass (mN>mWm_N > m_W) scenarios is shown in the second and third columns. Subscripts aa and bb denote leptons and jets, respectively, and jLj_L denotes a non-b-tagged jet. The check marks indicate the corresponding input variable is used for the training.

Pre-fit distributions of the $\Delta R$ between the two leptons

Figure 2. Pre-fit distributions of the ΔR\Delta R between the two leptons are shown for dielectron (left) and dimuon (right) events as examples of the input variables for the training. The upper row and the lower row correspond to the mN<mWm_{\mathrm{N}} < m_{\mathrm{W}} and mN>mWm_{\mathrm{N}} > m_{\mathrm{W}} selections, respectively, which differ in the required number of non-b-tagged jets. The cyan line is the ttˉ\mathrm{t\bar{t}} simulation used for the BDT training. Signal distributions are normalized to the total background. The hatched bands show the total uncertainty, including the systematic component. The last bins include the overflow entries. The lower panels show the ratio of the observed data over the background estimate.

Pre-fit distributions of the missing transverse momentum for dielectron and dimuon events

Figure 5. Pre-fit distributions of the missing transverse momentum are shown for dielectron (left) and dimuon (right) events as examples of the input variables for the training. The upper row and the lower row correspond to the mN<mWm_\mathrm{N} < m_\mathrm{W} and mN>mWm_\mathrm{N} > m_\mathrm{W} selections, respectively, which differ in the required number of non-b-tagged jets. The cyan line is the ttˉt\bar{t} simulation used for the BDT training. Signal distributions are normalized to the total background. The hatched bands show the total uncertainty, including the systematic component. The last bins include the overflow entries. The lower panels show the ratio of the observed data over the background estimate.

Pre-fit distributions of the invariant mass of the two leptons for dielectron and dimuon events

Figure 3. Pre-fit distributions of the invariant mass of the two leptons are shown for dielectron (left) and dimuon (right) events as examples of the input variables for the training. The upper row and the lower row correspond to the mN<mWm_{\mathrm{N}} < m_{\mathrm{W}} and mN>mWm_{\mathrm{N}} > m_{\mathrm{W}} selections, respectively, which differ in the required number of non-b-tagged jets. The cyan line is the ttˉ\mathrm{t\bar{t}} simulation used for the BDT training. Signal distributions are normalized to the total background. The region between two magenta lines depicts what is removed by the ∣m(ee)−91.2 GeV∣>10 GeV\lvert m(\mathrm{ee}) - 91.2\ \mathrm{GeV}\rvert> 10\ \mathrm{GeV} requirement in the dielectron channel. The blue line is the ttˉ\mathrm{t\bar{t}} simulation and DY simulation combined, without the mZm_{\mathrm{Z}} window veto requirement. In order to match the mZm_{\mathrm{Z}} window veto requirement to the bin edges, the plots on the left have the bin size of 10 GeV, except for the [70 81.2] and [101.2, 110] GeV bins. The hatched bands show the total uncertainty, including the systematic component. The last bins include the overflow entries. The lower panels show the ratio of the observed data over the background estimate.

Pre-fit distributions of the leading lepton transverse momentum for dielectron and dimuon events

Figure 4. Pre-fit distributions of the leading lepton transverse momentum are shown for dielectron (left) and dimuon (right) events as examples of the input variables for the training. The upper row and the lower row correspond to the mN<mWm_{\mathrm{N}} < m_{\mathrm{W}} and mN>mWm_{\mathrm{N}} > m_{\mathrm{W}} selections, respectively, which differ in the required number of non-b-tagged jets. The cyan line is the tt‾\mathrm{t}\overline{\mathrm{t}} simulation used for the BDT training. Signal distributions are normalized to the total background. The hatched bands show the total uncertainty, including the systematic component. The last bins include the overflow entries. The lower panels show the ratio of the observed data over the background estimate.

Background estimation

Backgrounds are classified into four categories: nonprompt leptons, lepton charge mismeasurement, conversion, and prompt SS dileptons, as listed below. Each background is discussed in detail after the list. The SR and the VRs used to validate the background estimates are summarized in Table 3.

SelectionSRNonprompt VRConversion VRChgMis VR
N(ℓ)N(\ell)2332
Lepton chargeSS——SS
Lepton pTp_{\mathrm{T}} (GeV)e1,2>25,15e_{1,2}>25,15; μ1,2>20,15\mu_{1,2}>20,15e1,2>25,15e_{1,2}>25,15; μ1,2>20,15\mu_{1,2}>20,15e1,2>25,15e_{1,2}>25,15; μ1,2>20,15\mu_{1,2}>20,15e1,2>25,15e_{1,2}>25,15
Lepton flavoree; μμ\mu\mueeμee\mu or eee; μμe\mu\mu e or μμμ\mu\mu\mueeμee\mu or eee; μμe\mu\mu e or μμμ\mu\mu\muee
N(j)N(j)≥3\geq3———
N(b)N(b)≥1\geq100—
N(JL)N(J_{\mathrm{L}})≥1\geq1 if mN<mWm_{\mathrm{N}}<m_{\mathrm{W}}; ≥2\geq2 if mN>mWm_{\mathrm{N}}>m_{\mathrm{W}}———
Mass window (GeV)∣m(ee)−mZ∣>10\lvert m(\mathrm{ee})-m_{\mathrm{Z}}\rvert>10; m(μμ)>4m(\mu\mu)>4OS SF pair; mass ∈[81.2,101.2]\in[81.2,101.2]m(ℓℓℓ)∈[81.2,101.2]m(\ell\ell\ell)\in[81.2,101.2]m(ee)∈[76.2,106.2]m(\mathrm{ee})\in[76.2,106.2]

Table 3. Summary of event selection criteria for the SR and VRs.

  1. Nonprompt-lepton backgrounds: the contribution of nonprompt leptons originating from jets is estimated using data in the AR and the matrix method [68], which uses the lepton misidentification rate as a function of lepton pTp_{\mathrm{T}} and η\eta, as measured in the BMR.

  1. Lepton charge mismeasurement backgrounds (ChgMis): a lepton charge can be mismeasured, causing an OS lepton pair to appear as an SS pair. The lepton charge mismeasurement rate (CMR) is measured using simulation and corrected to match that of the data as a function of the lepton pTp_{\mathrm{T}} and η\eta.

  1. Conversion backgrounds: leptons from photon conversion can pass the ID requirements and pair with a prompt lepton of the same charge. The contribution from this background source is estimated using simulations.

  1. Prompt-lepton backgrounds: backgrounds that contain only the prompt leptons from W or Z bosons decays are estimated using simulation.

Nonprompt-lepton backgrounds

Nonprompt-lepton backgrounds refer to processes with fewer than two prompt leptons that still pass the baseline selection by including at least one jet-induced nonprompt lepton. This background constitutes the dominant contribution, accounting for approximately 74 (87)% of the total background in the electron (muon) final state. The nonprompt-lepton backgrounds are estimated with the matrix method. This method extrapolates the number of events in the SR from AR, defined by the same baseline selection as the SR but using the loose lepton ID requirements instead of the tight ID ones. The AR thus serves as the sideband region in the lepton ID variables. The extrapolation from the AR to the SR is done through an extrapolation factor, defined as the fraction of leptons that pass the tight ID requirements relative to those that pass the loose ID ones. The extrapolation factor is measured in an independent BMR enriched with dijet events.

The extrapolation factor for nonprompt leptons is called the “misidentification rate”, while that for prompt leptons is called the “prompt rate”. The composition of the observed events (N⃗obs\vec{N}_{\mathrm{obs}}) is categorized by how many leptons pass the tight ID requirement (T) or fail the tight ID requirements but pass the loose ID requirements (L). The composition of the underlying events (N⃗origin\vec{N}_{\mathrm{origin}}) is categorized by the origin of each lepton as prompt (p) or nonprompt (n). The composition of the observed events and that of the underlying events are connected by a connection matrix (M) constructed from prompt rates (ε(p)\varepsilon(\mathrm{p})) and misidentification rates (ε(n)\varepsilon(\mathrm{n})). Each element of the connection matrix represents the probability that a specific lepton origin configuration is observed in a given ID configuration:

N⃗obs=(NTT  NTL  NLT  NLL)T,\vec{N}_{\mathrm{obs}} = (N_{\mathrm{TT}}\ \ N_{\mathrm{TL}}\ \ N_{\mathrm{LT}}\ \ N_{\mathrm{LL}})^{\mathrm{T}},
N⃗origin=(Npp  Npn  Nnp  Nnn)T,\vec{N}_{\mathrm{origin}} = (N_{\mathrm{pp}}\ \ N_{\mathrm{pn}}\ \ N_{\mathrm{np}}\ \ N_{\mathrm{nn}})^{\mathrm{T}},
N⃗obs=MN⃗origin.\vec{N}_{\mathrm{obs}} = \mathrm{M}\vec{N}_{\mathrm{origin}}.

The matrix elements of M, which connect the observed ID composition NID1ID2N_{\mathrm{ID}_1\mathrm{ID}_2} (with IDk=T\mathrm{ID}_k = \mathrm{T} or L\mathrm{L}) and the underlying origin composition No1o2N_{o_1o_2} (with ok=po_k = \mathrm{p} or n\mathrm{n}), are given by

M(ID1,ID2,o1,o2)=∏k=12[ε(ok)]δIDk,T[1−ε(ok)]δIDk,L.(1)\mathrm{M}(\mathrm{ID}_1,\mathrm{ID}_2,o_1,o_2) = \prod_{k=1}^{2} \left[\varepsilon(o_k)\right]^{\delta_{\mathrm{ID}_k,\mathrm{T}}} \left[1-\varepsilon(o_k)\right]^{\delta_{\mathrm{ID}_k,\mathrm{L}}}. \tag*{(1)}

The underlying composition can be obtained by applying the inverse of the matrix (M−1\mathrm{M}^{-1}) to the observed ID composition of events. Because the contribution from the prompt background is small, ε(p)\varepsilon(\mathrm{p}) is assumed to be unity, as its effect on the nonprompt background is negligible.

Lepton misidentification rates in data are measured in single-lepton events passing single-lepton triggers. Two different single-lepton triggers described in Table 1 are used depending on the lepton pTp_{\mathrm{T}}. Events containing only jets produced through the strong interaction, referred to as QCD multijet events, are selected for the measurement, requiring exactly one loose ID lepton and at least one jet with pT>40p_{\mathrm{T}} > 40 GeV and ΔR(ℓ,j)>1\Delta R(\ell,j) > 1. To reduce the contamination from processes with prompt leptons, such as W boson production in association with jets, pTmiss<25p_{\mathrm{T}}^{\mathrm{miss}} < 25 GeV and mT(p⃗Tℓp⃗Tmiss)<25m_{\mathrm{T}}(\vec{p}_{\mathrm{T}}^{\ell}\vec{p}_{\mathrm{T}}^{\mathrm{miss}}) < 25 GeV are also required. Residual prompt-lepton contributions after the selection are estimated using simulation and are subtracted in the calculation of the misidentification rates. The misidentification rate ranges from 0.18 (0.10) to 0.46 (0.25) for an electron (muon), depending on the lepton pTp_{\mathrm{T}} and η\eta. A closure test is performed on the ttˉ\mathrm{t\bar{t}} simulation. For this closure test, the nonprompt-lepton events estimated using misidentification rates obtained from simulation and the true nonprompt-lepton events identified using generator-level information are compared in the SR. The result of the closure test is presented in Figure 6. The closure test shows agreement within 10 (15)% for dielectron (dimuon) events. The nonprompt-lepton background estimation is validated by comparing it to the data in a separate VR with three leptons, where an OS SF lepton pair from the Z boson decay is selected by requiring ∣m(ℓℓ)−mZ∣<10|m(\ell\ell)-m_{\mathrm{Z}}| < 10 GeV, and the third lepton is assumed to be the nonprompt lepton. In the case of three SF leptons, among two possible combinations of OS lepton pairs, the lepton pair from Z boson is defined as the pair with m(ℓℓ)m(\ell\ell) closer to mZm_{\mathrm{Z}}. The pTp_{\mathrm{T}} distribution of the nonprompt lepton in the VR is shown in Figure 7, with good agreement between data and the background prediction.

Four-panel comparison of BDT classifier distributions for dielectron and dimuon events

Figure 6. Comparison of the BDT classifier distributions between the ttˉt\bar{t} simulated events with nonprompt leptons from the generator-level information passing the tight ID working point and the expected distribution of the ttˉt\bar{t} simulation for dielectron (left) and dimuon (right) events in the SR for mN=50m_\mathrm{N} = 50 GeV (upper) and for mN=100m_\mathrm{N} = 100 GeV (lower). The hatched uncertainty bands are ±10\pm10 (±15\pm15)% for the dielectron (dimuon) channel. The lower panels show the ratio of the ttˉt\bar{t} simulated events with nonprompt leptons from the generator-level information over the expected distribution.

Distributions of the transverse momentum of the electron and muon not originating from the Z boson decay

Figure 7. Distributions of pTp_{\mathrm{T}} of the electron (left) or muon (right) not originating from the Z boson decay in the trilepton VR are shown for the nonprompt-lepton background validation. The hatched bands show the total uncertainty, including the systematic component. The last bins include the overflow entries. The lower panels show the ratio of the observed data over the background estimate.

Lepton charge mismeasurement backgrounds

Lepton charge mismeasurement backgrounds are processes with two OS leptons in which the reconstructed charge of one of the leptons is mismeasured. From the earlier measurement in CMS, the CMR for muons is known to be less than 10−510^{-5} and does not produce a measurable amount of background, so it is ignored [69]. The SM processes that contribute to this background include DY, ttˉ\mathrm{t\bar{t}}, diboson, and single top quark production processes. A direct estimate of the events with two SS leptons from the respective simulated samples suffers from low event counts and results in large statistical uncertainties. Consequently, we estimate this background by measuring the CMR in data and applying it to events in the simulated samples that pass the SR selection, except for the SF OS leptons requirement. The electron CMR is determined from simulated Z →\to ee events as a function of the electron pTp_{\mathrm{T}} and η\eta, and is subsequently corrected to match the observation in data. The electron CMR in simulation is measured by comparing the number of charge-mismeasured electrons and the number of electrons with the correct charge using the generator-level information. The electron CMR in data is measured using the Z →\to ee resonance. Events must have exactly two electrons, with pTp_{\mathrm{T}} larger than 25 and 15 GeV, respectively, and pass the tight ID with the m(ee)m(\mathrm{ee}) between 60 and 120 GeV. The number of SS and OS Z →\to ee events is extracted by fitting the m(ee)m(\mathrm{ee}) distributions to the template shape from simulation.

An iterative approach is used to correct the CMR in simulation to match that in data. At each step of the iteration, the CMR ratio between data and simulation is measured and used to scale the number of simulated events with a charge-mismeasured electron. Because both electrons in each event contribute to their corresponding CMR bin depending on whether their charge is mismeasured or not, the new CMR for the simulation does not perfectly agree with the data CMR. The iteration of scaling the CMR in simulation stops when the magnitude of the corrections for the CMR falls below 1.5% for all pTp_{\mathrm{T}}-η\eta bins. After the iteration, the product of the correction factors at each step becomes the total correction factor. The total correction factor ranges from 0.74 to 1.61, depending on the electron pTp_{\mathrm{T}}, η\eta, and the data-taking era. The electron charge mismeasurement background is validated in data using charge-mismeasured DY events with tighter selection on the m(ee)m(\mathrm{ee}) than the one used in the data CMR measurement region: the two electrons from the Z boson decay are measured to have the same charge, and their invariant mass is within 15 GeV of mZm_{\mathrm{Z}}. The invariant mass distribution shows good agreement between data and the background prediction within the uncertainty, as shown in Figure 8.

Invariant mass distribution for same-sign electron pairs

Figure 8. The invariant mass distribution for the SS electron pairs in the lepton charge mismeasurement VR. The hatched bands show the total uncertainty, including the systematic component. The lower panels show the ratio of the observed data over the background estimate.

Conversion backgrounds

The conversion backgrounds include internal and external photon conversions, depending on whether the conversion occurs at the collision point or within the detector, respectively. Since the rate of external conversion is proportional to the inverse of the produced lepton mass squared [70], it is only significant for electrons. For the validation of this background, we select events with three leptons: two from a Z boson decay and an additional lepton from a photon conversion. The invariant mass of the three leptons is required to be within a 10 GeV window around mZm_{\mathrm{Z}}. In this three-lepton VR, we observe reasonable agreement between data and the background prediction within the uncertainty, as shown in Figure 9. Therefore, the conversion backgrounds are taken from the simulation.

Invariant mass distributions for conversion electron and muon validation regions

Figure 9. The invariant mass distribution of three leptons is shown for the Z boson plus conversion electron (left) or muon (right) events for the validation of the conversion background. The hatched bands show the total uncertainty, including the systematic component. The lower panels show the ratio of the observed data over the background estimate.

Prompt-lepton backgrounds

Prompt-lepton backgrounds originate from processes involving at least two SS leptons produced in the decays of V bosons, including leptons from decays of τ\tau leptons originating from these bosons. These contributions are estimated using simulated samples. Major contributions come from ttˉ\mathrm{t\bar{t}}V, ttˉ\mathrm{t\bar{t}}H, WZ, and ZZ production. Processes with small expected contributions are combined as “rare” background, dominated by the production of tZ, VVV, ggH, and qqH.

Systematic uncertainties

Various sources of systematic uncertainties are considered in the analysis, categorized as background, experimental, and theoretical uncertainties. Experimental and theoretical uncertainties are applied to both signal and simulated backgrounds. The systematic uncertainties are implemented as nuisance parameters in the analysis. Depending on their nature, the uncertainties are treated as correlated or uncorrelated between the data-taking eras. The summary of the uncertainties is shown in Table 4.

Table 4: Summary of systematic uncertainties. Each uncertainty is implemented as either a rate or shape uncertainty. Uncertainties affecting only the normalization are marked as “Rate” in the table. Those affecting BDT template shapes are indicated as “Shape.” For shape uncertainties, “—” indicates that its effect is propagated to the BDT score distribution through the corresponding up and down templates and may depend on the process, channel, and HNL mass hypothesis.

SourceAssigned variationType
Background uncertainties
Nonprompt leptons20%Rate
CMR fitting model—Shape
Charge-mismeasurement closure test6%Rate
Electron energy shifts—Shape
Conversions20%Rate
Experimental uncertainties
Lepton ID efficiency—Shape
Electron reconstruction—Shape
Electron energy correction—Shape
Muon momentum scale—Shape
Jet energy scale—Shape
Jet energy resolution—Shape
Pileup jet identification—Shape
b tagging efficiency corrections—Shape
Integrated luminosity1.6%Rate
Inelastic pp cross section4.6%Rate
Trigger efficiencies1%Rate
L1-prefiring inefficiency—Shape
Theoretical uncertainties
Signal PDF and αS\alpha_{\mathrm{S}} variations—Shape
Signal μR\mu_{\mathrm{R}} and μF\mu_{\mathrm{F}} variations—Shape
Signal parton shower modeling—Shape
ttˉ\mathrm{t}\bar{\mathrm{t}} cross section−4.4%+3.5%^{+3.5\%}_{-4.4\%}Rate
ttˉZ\mathrm{t}\bar{\mathrm{t}}\mathrm{Z} cross section11.9%Rate
ttˉW\mathrm{t}\bar{\mathrm{t}}\mathrm{W} cross section13.3%Rate
ttˉH\mathrm{t}\bar{\mathrm{t}}\mathrm{H} cross section10.0%Rate
WZ cross section12%Rate
ZZ cross section6.4%Rate
tZ cross section5.2%Rate
VVV cross section50%Rate
ggH cross section50%Rate
qqH cross section50%Rate
VH cross section50%Rate

Table 4. Summary of systematic uncertainties. Each uncertainty is implemented as either a rate or shape uncertainty. Uncertainties affecting only the normalization are marked as “Rate” in the table. Those affecting BDT template shapes are indicated as “Shape.” For shape uncertainties, “—” indicates that its effect is propagated to the BDT score distribution through the corresponding up and down templates and may depend on the process, channel, and HNL mass hypothesis.

Background systematic uncertainties

The misidentification rates in nonprompt-lepton backgrounds contribute the largest uncertainty in the total background estimation. We use two methods for assessing the uncertainties in the nonprompt-lepton background. First, from the closure test described in § 5, an uncertainty of 10 (15)% is set in the dielectron (dimuon) channel. Secondly, the uncertainty is assessed by analyzing three sources of uncertainty in the misidentification rate itself. The normalization of the prompt-lepton background is varied up and down by the theoretical cross section uncertainty when subtracting it from data. The minimum jet pTp_{\mathrm{T}} requirement for the BMR selection is varied up and down by 10 or 5 GeV, depending on the single-lepton trigger threshold. The dependence of the measured misidentification rate on the jet flavor is quantified by requiring at least one b-tagged jet among the jets. The relative contributions of these three sources of uncertainty vary with the lepton pTp_{\mathrm{T}} and η\eta. In the SR, the combined uncertainty in the nonprompt-lepton background after propagating the effects of the three uncertainty sources is within a 20% band around the nominal BDT score distribution. To conservatively account for the uncertainties assessed from the closure test and the misidentification rate, a flat 20% uncertainty is assigned to the nonprompt-lepton background. This uncertainty is treated as uncorrelated between the data-taking eras.

Three sources of uncertainty are considered in the estimation of the charge-mismeasured background. ground. The first uncertainty arises from the measurement of the correction factors for the CMR. This uncertainty is quantified by using different fitting models for data CMR extraction, and it is found to vary between 2 and 7% across the pTp_{\mathrm{T}}-η\eta bins. Since this uncertainty is dominated by the statistical component, it is treated as uncorrelated between the data-taking eras. The second uncertainty accounts for an imperfect parameterization of the CMR due to the limited sample size, and it is estimated to be 6% by performing a closure test on ttˉt\bar{t} simulated events. The test consists of estimating the number of SS dielectron events by applying the measured CMR to the OS dielectron events and comparing the result to the observed number of SS dielectron events. This uncertainty is treated as correlated between the data-taking eras. The last source of uncertainty originates from the electron energy measurement. Highly energetic bremsstrahlung can deflect the electron path, leading to a mismeasurement of the electron charge. Additionally, the emitted radiation may become spatially separated from the electron and not be taken into account in the electron energy reconstruction. To account for this, we emulate the energy underestimation in the simulated events with OS dielectron pairs by shifting the electron pTp_{\mathrm{T}} in the VR. The nominal shift is optimized by minimizing the χ2\chi^2 difference between the OS and SS m(ee)m(ee) distributions over a ±3%\pm3\% energy shift scan. We evaluate the systematic uncertainty by repeating this optimization in two kinematic regions, using events with both electrons pT>40p_{\mathrm{T}} > 40 GeV or both pT<40p_{\mathrm{T}} < 40 GeV. The maximum absolute deviation from the nominal shift is assigned as a symmetric systematic uncertainty. The uncertainty in the electron pTp_{\mathrm{T}} from this correction is between 0.7 and 1.5%, depending on the electron η\eta and data-taking era. This uncertainty is treated as correlated between the data-taking eras.

For the conversion background, a 20% uncertainty is assigned to the rate based on the agreement in the three-lepton VR. The simulation underestimates the conversion event rate in the three-lepton VR by approximately 20% in all data-taking eras for dimuon events, while the data-to-simulation difference for electron conversions varies between 0 and 20%, depending on the data-taking eras. Consequently, the muon conversion background uncertainty is treated as correlated, while the electron one is treated as uncorrelated, between the data-taking eras.

Experimental systematic uncertainties

The uncertainty in lepton ID is included, ranging from 1 to 4% depending on the lepton pTp_{\mathrm{T}} and η\eta, and is treated as correlated between the data-taking eras. Lepton ID correction factors are measured using the “tag-and-probe” method [71], where the efficiency is measured as the fraction of the reconstructed leptons from DY events that pass the ID requirement. The uncertainty is then determined by changing the minimum pTp_{\mathrm{T}} requirement on the leptons, the fit function for the m(ℓℓ)m(\ell\ell) distribution, and the Monte Carlo generator.

For electrons, there are uncertainties associated with the efficiency of electron reconstruction and energy correction [27]. These uncertainties are treated as correlated between different data-taking eras. In the case of muons, the reconstruction efficiency is close to unity, and no statistically significant evidence for the difference between data and simulation is observed. The muon energy scale uncertainty is considered, and it is treated as correlated between the data-taking eras [28].

Uncertainties in the jet energy scale and resolution are assessed by reevaluating observables in simulated events with the jet energy shifted by its uncertainty. The systematic uncertainties associated with the jet energy scale and resolution corrections are considered uncorrelated between the data-taking eras. These uncertainties are propagated to the pTmissp_{\mathrm{T}}^{\mathrm{miss}} uncertainty. Additionally, the uncertainty in the energy scale of PF particles not associated with jets, leptons, and photons, known as unclustered energy, is propagated to the pTmissp_{\mathrm{T}}^{\mathrm{miss}} uncertainty [42]. The uncertainties in the b tagging efficiency correction are evaluated by varying the correction factors for light-quark and gluon, c quark, and b quark jets in simulation [40]. Each b tagging uncertainty has both correlated and uncorrelated sources between the data-taking eras. The uncertainty in the efficiency correction factor corresponding to the pileup jet identification is also applied and treated as uncorrelated between the data-taking eras [38].

Other experimental uncertainties have small impacts. The uncertainty in the integrated luminosity [43, 44, 45] has both correlated and uncorrelated components between the data-taking eras. The overall integrated luminosity uncertainty for the 2016–2018 data set is 1.6%. The uncertainty in the number of pileup interactions is estimated by applying a ±4.6%\pm4.6\% variation to the total inelastic pp cross section used to reweight simulated events to match the pileup distribution observed in data [72] and is treated as correlated between the data-taking eras. The uncertainties in the trigger efficiencies are measured using the tag-and-probe method described earlier and are treated as correlated between the data-taking eras. During the 2016 and 2017 data-taking periods, a gradual shift in the timing of the inputs of the ECAL L1 trigger caused a trigger inefficiency, referred to as the “L1-prefiring”. An uncertainty related to the inefficiency from the L1-prefiring in 2016 and 2017 is applied and treated as uncorrelated between the data-taking eras [73].

The statistical uncertainty of the data used for the nonprompt-lepton background estimations, as well as that of the simulated samples used in all background estimations, is considered.

Theoretical systematic uncertainties

An uncertainty of −4.4+3.5%^{+3.5}_{-4.4}\% is assigned to the signal production cross section [50]. Theoretical uncertainties in the signal acceptance from the uncertainties in the PDFs and in the renormalization (μR\mu_\mathrm{R}) and factorization (μF\mu_\mathrm{F}) scales are considered. For the uncertainty arising from PDFs, the changes in the acceptance for 100 replicas of the PDF set are added in quadrature. The uncertainties in the signal acceptance arising from μR\mu_\mathrm{R} and μF\mu_\mathrm{F} are studied by varying the scales by factors of 0.5 and 2 individually or simultaneously, excluding the two variations in the opposite directions, and taking the envelope of all variations. The uncertainty in the parton shower modeling is estimated by varying initial- and final-state radiation scales separately by factors of 0.5 and 2 and is defined by the envelope obtained from these variations. Among the signal theoretical uncertainties, the ones from initial- and final-state radiation, and from the ttˉ\mathrm{t\bar{t}} cross section are dominant.

Three theoretical uncertainties in the cross sections of the background processes are taken into account: μR\mu_\mathrm{R} and μF\mu_\mathrm{F}, PDFs, and strong coupling constant (αS\alpha_\mathrm{S}), added in quadrature. Because the uncertainties resulting from higher and lower scale choices are comparable in magnitude, the larger deviation is assigned as a symmetric uncertainty. For WZ production, the radiation zero effect is considered, whose magnitude decreases in the diagrams with additional partons in higher-order QCD calculations [74]. We assign an 11% uncertainty due to the radiation zero effect, resulting in a total uncertainty of 12%. For ZZ production, we apply a cross section uncertainty of 6.4% [62]. For ttˉ\mathrm{t\bar{t}} + V/H we apply the uncertainty from Ref. [75]. For ttˉ\mathrm{t\bar{t}}H processes, the uncertainty in the Higgs boson branching fraction is also included. The total uncertainty values obtained are 11.9% for ttˉ\mathrm{t\bar{t}}Z, 13.3% for ttˉ\mathrm{t\bar{t}}W, and 10.0% for ttˉ\mathrm{t\bar{t}}H. For the tZ process, the total uncertainty is 5.2% [76]. For VVV, ggH, qqH, and V H production, a conservative 50% uncertainty is assigned. All uncertainty sources with theoretical origins are considered as fully correlated between the data-taking eras.

Theoretical systematic uncertainties

Results

To extract the final results, a simultaneous binned maximum likelihood fit with profiling of nuisance parameters is performed over the 2016–2018 data set, separately for each mNm_{\mathrm{N}} hypothesis in the SR, using the BDT classifier score distribution. The following results have been determined using the CMS statistical analysis tool COMBINE [77], which is based on the ROOFIT [78] and ROOSTATS [79] frameworks. Nuisance parameters are categorized into two groups [80]: rate uncertainties are implemented with log-normal nuisance parameters, while shape uncertainties are incorporated through variations of the BDT templates with Gaussian-constrained nuisance parameters. A simplified Barlow–Beeston approach [81], with a single nuisance parameter per bin of the distribution, is used to treat the statistical uncertainties due to the limited size of the simulated samples. The uncertainties with the greatest impact for all mass hypotheses are the statistical uncertainties of the background templates in the individual bins. With its conservative approach, the nonprompt-lepton background uncertainty gets reduced by a factor of 2–4, depending on the data-taking era, lepton flavor, and mNm_{\mathrm{N}} hypothesis. To assess the compatibility between the observed data and the background-only hypothesis, a goodness-of-fit test is performed using the χ2\chi^{2} test statistic. The expected distribution of the test statistic is derived from 2000 background-only pseudo-experiments generated from the post-fit nuisance parameter configurations. The observed data yield a pp-value [82] of 0.3 or above for all signal models, indicating that the data are compatible with the background model. The fit is performed for each mass hypothesis to assess its dependence on the HNL mass. The observed and expected BDT score distributions after the fit in the SR are shown in Figure 10–Figure 12 for various mNm_{\mathrm{N}} hypotheses. The score distributions in data are consistent with those in the background. The figures show that the BDT classifier trained on the tt‾\mathrm{t}\overline{\mathrm{t}} process also yields low-score values for the majority of other minor backgrounds. Therefore, a classifier trained on dilepton and lepton+jets tt‾\mathrm{t}\overline{\mathrm{t}} events as the only background sample is sufficient to distinguish other backgrounds from the signal.

Comparison of observed and predicted post-fit BDT score distributions in the dielectron and dimuon channels

Figure 10. Comparison of the observed (points) and predicted (filled histograms) post-fit BDT score distributions in the dielectron (left) and dimuon (right) channels. The BDT score distributions correspond to classifiers trained with signal mass hypotheses of 20 (upper) and 50 GeV (lower). The expected signal distributions extracted with the same classifier are overlaid on the figures. They are normalized assuming B(t→bℓN)=2×10−4\mathcal{B}(t \to b\ell N) = 2 \times10^{-4}, a value chosen to ensure that the signal distributions are visible across all channels and mNm_N hypotheses. The hatched bands show the total uncertainty, including the systematic component. The lower panels show the ratio of the observed data over the background estimate.

Post-fit BDT score distributions in the dielectron and dimuon channels

Figure 12. Comparison of the observed (points) and predicted (filled histograms) post-fit BDT score distributions in the dielectron (left) and dimuon (right) channels. The BDT score distributions correspond to classifiers trained with signal mass hypothesis of 100 GeV. The expected signal distributions extracted with the same classifier are overlaid on the figures. They are normalized assuming B(t→bℓN)=2×10−4\mathcal{B}(\mathrm{t} \to\mathrm{b}\ell\mathrm{N}) = 2 \times10^{-4}, a value chosen to ensure that the signal distributions are visible across all channels and mNm_{\mathrm{N}} hypotheses. The hatched bands show the total uncertainty, including the systematic component. The lower panels show the ratio of the observed data over the background estimate.

Based on the classifier score distribution, an upper limit at 95% confidence level (CL) on the signal strength is determined using the CLs_{\mathrm{s}} criterion [83, 84] in a modified frequentist approach. The test statistic used is the ratio of profiled likelihoods defined according to the LHC-style prescription from Ref. [77]. In the conversion to the branching fraction B(t→bℓN)\mathcal{B}(\mathrm{t}\to\mathrm{b}\ell\mathrm{N}), the tt‾\mathrm{t}\overline{\mathrm{t}} production cross section at NNLO+NNLL QCD of 833.9 pb is used, and the signal decay rates are calculated at LO. The observed and expected limits on B(t→bℓN)\mathcal{B}(\mathrm{t}\to\mathrm{b}\ell\mathrm{N}) for the mNm_{\mathrm{N}} values considered in the search are shown in Figure 13. The expected upper limits on B(t→bℓN)\mathcal{B}(\mathrm{t}\to\mathrm{b}\ell\mathrm{N}) vary from [7.1 (3.4)]×10−6[7.1\ (3.4)]\times10^{-6} to [3.0 (1.3)]×10−5[3.0\ (1.3)]\times10^{-5} in the dielectron (dimuon) channel. The observed upper limits vary from [11 (3.7)]×10−6[11\ (3.7)]\times10^{-6} to [4.8 (2.4)]×10−5[4.8\ (2.4)]\times10^{-5} in the dielectron (dimuon) channel, and the results are consistent with the background-only hypothesis. The upper limits increase rapidly when the mNm_{\mathrm{N}} value approaches the mWm_{\mathrm{W}} value because of the vanishing lepton signal acceptance and the challenging signal classification. For this reason, we do not present an interpolation of the measured limits in the mNm_{\mathrm{N}} range [75] [85] GeV. The decreasing lepton acceptance is also responsible for the increase in the expected limits shown for mN<50m_{\mathrm{N}}<50 GeV. The upper limits on the branching fraction for the electron neutrino mixing scenario are larger than those for the muon neutrino mixing scenario by approximately a factor of two. This originates from the lower trigger acceptance and ID efficiency of electrons compared to those of muons.

Expected and observed upper limits on the branching fraction as functions of the heavy-neutrino mass

Figure 13. The expected and observed upper limits at 95% CL on the branching fraction, B(t→bℓN)\mathcal{B}(\mathrm{t} \to\mathrm{b}\ell\mathrm{N}), as functions of mNm_{\mathrm{N}}. The left and right figures correspond to the electron and muon neutrino mixing scenarios, respectively.

For a straightforward comparison with earlier type-I seesaw model search results, the upper limits on the B(t→bℓN)\mathcal{B}(\mathrm{t}\to\mathrm{b}\ell\mathrm{N}) are converted into the upper limits on the model parameters, ∣VeN∣2\lvert V_{\mathrm{eN}}\rvert^{2} and ∣VμN∣2\lvert V_{\mu\mathrm{N}}\rvert^{2}, as shown in Figure 14. The expected upper limits on ∣VeN∣2\lvert V_{\mathrm{eN}}\rvert^{2} (∣VμN∣2\lvert V_{\mu\mathrm{N}}\rvert^{2}) vary between [12 (5.6)]×10−5[12\ (5.6)]\times10^{-5} and [9.9 (4.3)]×10−2[9.9\ (4.3)]\times10^{-2}. The observed upper limits on ∣VeN∣2\lvert V_{\mathrm{eN}}\rvert^{2} (∣VμN∣2\lvert V_{\mu\mathrm{N}}\rvert^{2}) vary between

Expected and observed upper limits on electron and muon neutrino mixing

Figure 14. The expected and observed upper limits at 95% CL on the model parameters, ∣VeN∣2|V_{eN}|^2 and ∣VμN∣2|V_{\mu N}|^2, as functions of mNm_N. The left and right figures correspond to the electron and muon neutrino mixing scenarios, respectively. For comparison, expected limits from the previous CMS searches in charged-current Drell–Yan production [10, 11] and the ATLAS ttˉt\bar{t} production search [20] are overlaid. The orange dash-dotted line is the SS dilepton channel search in CMS, the magenta dashed line is the trilepton channel search in CMS, and the blue dotted line is the ATLAS search.

Comparison of observed and predicted post-fit BDT score distributions in the dielectron and dimuon channels

Figure 11. Comparison of the observed (points) and predicted (filled histograms) post-fit BDT score distributions in the dielectron (left) and dimuon (right) channels. The BDT score distributions correspond to classifiers trained with signal mass hypotheses of 75 (upper) and 85 GeV (lower). The expected signal distributions extracted with the same classifier are overlaid on the figures. They are normalized assuming B(t→bℓN)=2×10−4B(t \to b\ell N) = 2 \times10^{-4}, a value chosen to ensure that the signal distributions are visible across all channels and mNm_N hypotheses. The hatched bands show the total uncertainty, including the systematic component. The lower panels show the ratio of the observed data over the background estimate.

between [16 (4.2)]×10−5[16\ (4.2)] \times10^{-5} and [11 (3.9)]×10−2[11\ (3.9)] \times10^{-2}. The existing results from the direct search channels in CMS and from the ATLAS search for HNLs produced in top quark decays are compared with those from this search. The expected results in this search are comparable to the existing results from the dilepton channel and weaker by approximately a factor of ten than those from the trilepton channel for mN<mWm_N < m_W. At higher mass values, this search becomes less sensitive compared to results from the direct search channels. This is because the signal decay process, t→bℓNt \to b\ell N, is no longer mediated by an on-shell WW boson, i.e., t→bWt \to bW and W→ℓNW \to\ell N, meaning the corresponding partial decay width, Γ(t→bℓN)\Gamma(t \to b\ell N), is significantly reduced. For 20≤mN≤75 GeV20 \leq m_N \leq75\ \mathrm{GeV}, the expected limits on ∣VeN∣2|V_{eN}|^2 and ∣VμN∣2|V_{\mu N}|^2 from this search are stronger by approximately a factor of two than those from the ATLAS search for HNLs produced in top quark decays [20]. The size of the difference varies with mNm_N.

Summary

A search for heavy neutral leptons (HNLs) produced in top quark decays in proton-proton collisions at s=13 TeV\sqrt{s} = 13\ \mathrm{TeV}, corresponding to an integrated luminosity of 138 fb−1138\ \mathrm{fb}^{-1} collected with the CMS detector at the LHC, has been presented. The analysis targets HNL masses between 20 and 100 GeV in the type-I seesaw framework. Events are selected by requiring jets and a pair of same-sign, same-flavor leptons. A multivariate classifier is used to enhance the sensitivity of the analysis. No significant excess over the background prediction is observed. Observed (expected) upper limits at 95% confidence level are set on the branching fraction of the top quark decay into an HNL, a charged lepton, and a bottom quark, varying between [1.1 (0.71)]×10−5[1.1\ (0.71)] \times10^{-5} and [4.8 (3.0)]×10−5[4.8\ (3.0)] \times10^{-5} in the dielectron channel, and between [3.7 (3.4)]×10−6[3.7\ (3.4)] \times10^{-6} and [2.4 (1.3)]×10−5[2.4\ (1.3)] \times10^{-5} in the dimuon channel. The results are also interpreted as upper limits on the squared mixing amplitudes between the HNL and the standard model electron and muon neutrinos, denoted ∣VeN∣2|V_{eN}|^2 and ∣VμN∣2|V_{\mu N}|^2, respectively. The observed upper limits range from 16×10−516 \times10^{-5} to 11×10−211 \times10^{-2} for ∣VeN∣2|V_{eN}|^2 and from 4.2×10−54.2 \times10^{-5} to 3.9×10−23.9 \times10^{-2} for ∣VμN∣2|V_{\mu N}|^2.

Acknowledgments

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 and other centres 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, the CMS detector, and the supporting computing infrastructure provided by the following funding agencies: SC (Armenia), BMFWF and FWF (Austria); FNRS and FWO (Belgium); CNPq, CAPES, FAPERJ, FAPERGS, and FAPESP (Brazil); MES and BNSF (Bulgaria); CERN; CAS, MoST, and NSFC (China); MINCIENCIAS (Colombia); MSES and CSF (Croatia); RIF (Cyprus); SENESCYT (Ecuador); ERC PRG and PSG, TARISTU24-TK10 and MoER TK202 (Estonia); Academy of Finland, MEC, and HIP (Finland); CEA and CNRS/IN2P3 (France); SRNSF (Georgia); BMFTR, DFG, and HGF (Germany); GSRI (Greece); MATE and NKFIH (Hungary); DAE and DST (India); IPM (Iran); SFI (Ireland); INFN (Italy); MSIT and NRF (Republic of Korea); MES (Latvia); LMTLT (Lithuania); MOE and UM (Malaysia); BUAP, CINVESTAV, CONACYT, LNS, SEP, and UASLP-FAI (Mexico); MOS (Montenegro); MBIE (New Zealand); PAEC (Pakistan); MSHE, NSC, and NAWA (Poland); FCT (Portugal); MESTD (Serbia); MICIU/AEI and PCTI (Spain); MOSTR (Sri Lanka); Swiss Funding Agencies (Switzerland); NSTC (Taiwan); MHESI (Thailand); TUBITAK and TENMAK (Türkiye); NASU (Ukraine); STFC (United Kingdom); DOE and NSF (USA).

Individuals have received support from the Marie-Curie programme and the European Research Council and Horizon 2020 Grant, contract Nos. 675440, 724704, 752730, 758316, 765710, 824093, 101115353, 101002207, 101001205, and COST Action CA16108 (European Union); the Leventis Foundation; the Alfred P. Sloan Foundation; the Alexander von Humboldt Foundation; the Science Committee, project no. 22rl-037 (Armenia); the Fonds pour la Formation à la Recherche dans l’Industrie et dans l’Agriculture (FRIA) and Fonds voor Wetenschappelijk Onderzoek contract No. 1228724N (Belgium); the Beijing Municipal Science & Technology Commission, No. Z191100007219010, the Fundamental Research Funds for the Central Universities, the Ministry of Science and Technology of China under Grant No. 2023YFA1605804, the Natural Science Foundation of China under Grant No. 12535004, and USTC Research Funds of the Double First-Class Initiative No. YD2030002017 (China); the Ministry of Education, Youth and Sports (MEYS) of the Czech Republic; the Shota Rustaveli National Science Foundation (Georgia); the Deutsche Forschungsgemeinschaft (DFG), among others, under Germany’s Excellence Strategy – EXC 2121 “Quantum Universe” – 390833306, and under project number 400140256 - GRK2497; the Hellenic Foundation for Research and Innovation (HFRI), Project Number 2288 (Greece); the Hungarian Academy of Sciences, the New National Excellence Program - ÚNKP, the NKFIH research grants K 131991, K 138136, K 143460, K 143477, K 147557, K 146913, K 146914, K 147048, TKP2021-NKTA-64, 2026-4.1.1-MISSZIÓK-2026-00005 and 2025-1.1.5-NEMZ_KI-2025-00004, and MATE KKP and KKPcs Research Excellence and Flagship Research Groups grants (Hungary); the Council of Science and Industrial Research, India; ICSC – National Research Centre for High Performance Computing, Big Data and Quantum Computing, FAIR – Future Artificial Intelligence Research, and CUP I53D23001070006 (Mission 4 Component 1), funded by the NextGenerationEU program, the Italian Ministry of University and Research (MUR) under Bando PRIN 2022 – CUP I53C24002390006, PRIN PRIMULA 2022RBYK7T (Italy); the Latvian Council of Science; the Ministry of Science and Higher Education, project no. 2022/WK/14, and the National Science Centre, contracts Opus 2021/41/B/ST2/01369, 2021/43/B/ST2/01552, 2023/49/B/ST2/03273, and the NAWA contract BPN/PPO/2021/1/00011 (Poland); the Fundaçao para a Ciência e a Tecnologia (Portugal); the National Priorities Research Program by Qatar National Research Fund; MICIU/AEI/10.13039/501100011033, ERDF/EU, “European Union NextGenerationEU/PRTR”, projects PID2022-142604OB-C21, PID2022-139519OB-C21, PID2023-147706NB-I00, PID2023-148896NB-I00, PID2023-146983NB-I00, PID2023-147115NB-I00, PID2023-148418NB-C41, PID2023-148418NB-C42, PID2023-148418NB-C43, PID2023-148418NB-C44, PID2024-158190NB-C22, RYC2021-033305-I, RYC2024-048719-I, CNS2023-144781, CNS2024-154769 and Plan de Ciencia, Tecnología e Innovación de Asturias, Spain; the Chulalongkorn Academic into Its 2nd Century Project Advancement Project, the National Science, Research and Innovation Fund program IND FF 68 369 2300 097, and the Program Management Unit for Human Resources & Institutional Development, Research and Innovation, grant B39G680009 (Thailand); the Eric & Wendy Schmidt Fund for Strategic Innovation through the CERN Next Generation Triggers project under grant agreement number SIF-2023-004; the Kavli Foundation; the Nvidia Corporation; the SuperMicro Corporation; the Welch Foundation, contract C-1845; and the Weston Havens Foundation (USA).

Data availability

Release and preservation of data used by the CMS Collaboration as the basis for publications is guided by the CMS data preservation, re-use and open access policy.

The CMS Collaboration

Yerevan Physics Institute, Yerevan, Armenia

A. Gevorgyan, A. Hayrapetyan, V. Makarenko, A. Tumasyan1

Marietta Blau Institute for Particle Physics, Vienna, Austria

W. Adam, L. Benato, T. Bergauer, M. Dragicevic, S. Gundacker, A.K. Guven, P.S. Hussain, M. Jeitler2, N. Krammer, A. Li, D. Liko, M. Matthewman, J. Schieck2, R. Schöfbeck2, M. Shooshtari, M. Sonawane, N. Van Den Bossche, W. Waltenberger, C.E. Wulz2

Universiteit Antwerpen, Antwerpen, Belgium

T. Janssen, D. Ocampo Henao, T. Van Laer, P. Van Mechelen

Vrije Universiteit Brussel, Brussel, Belgium

D. Ahmadi, J. Bierkens, N. Breugelmans, S. Dansana, A. De Moor, M. Delcourt, S. Duponcheel, C. Gupta, F. Heyen, Y. Hong, K. Kang, P. Kashko, S. Lowette, I. Makarenko, S. Nandakumar, J. Niedziela, S. Tavernier, M. Tytgat3, G.P. Van Onsem, S. Van Putte, T. Wybouw

Université Libre de Bruxelles, Bruxelles, Belgium

A. Beshr, B. Bilin, F. Caviglia Roman, B. Clerbaux, A.K. Das, I. De Bruyn, G. De Lentdecker, E. Ducarme, H. Evard, L. Favart, I. Kalaitzidou, A. Khalilzadeh, A. Malara, A. Potrebko, M.A. Shahzad, L. Thomas, M. Vanden Bemden, C. Vander Velde, P. Vanlaer, C. Yuan, F. Zhang

Ghent University, Ghent, Belgium

A. Cauwels, M. De Coen, D. Dobur4, C. Giordano, G. Gokbulut, K. Kaspar, D. Kavtaradze, D. Marckx, A. Mehta, K. Skovpen, A.M. Tomaru, J. van der Linden, J. Vandenbroeck

Université Catholique de Louvain, Louvain-la-Neuve, Belgium

H. Aarup Petersen, A. Benecke, A. Bethani, G. Bruno, A. Cappati, J. De Favereau De Jeneret, C. Delaere, F. Gameiro Casalinho, A. Giammanco, A.O. Guzel, M. Hussain, Z. Lawrence, V. Lemaitre, J. Lidrych, P. Malek, S. Turkcapar

Centro Brasileiro de Pesquisas Fisicas, Rio de Janeiro, Brazil

G. Alves, E. Coelho, M.V. Gonçalves Sales, C. Hensel, D. Matos Figueiredo, T. Menezes De Oliveira, C. Mora Herrera, P. Rebello Teles, M. Soeiro, E.J. Tonelli Manganote5, A. Vilela Pereira

Universidade do Estado do Rio de Janeiro, Rio de Janeiro, Brazil

W.L. Aldá Júnior, M. Barroso Ferreira Filho, H. Brandao Malbouisson, W. Carvalho, J. Chinellato6, G. Correia Silva, M. Costa Reis, E.M. Da Costa, D. Da Silva Dalto, G.G. Da Silveira7, D. De Jesus Damiao, S. Fonseca De Souza, R. Gomes De Souza, S. Jesus, T. Laux Kuhn7, K. Maslova Gioseffi Defante, K. Mota Amarilo, L. Mundim, H. Nogima, J.P. Pinheiro, A. Santoro, A. Sznajder, M. Thiel, F. Torres Da Silva De Araujo8, D. Torres Machado

Universidade Estadual Paulista (a), Universidade Federal do ABC (b), São Paulo, Brazil

C.A. Bernardes, L. Calligaris, J. Carvalho Leite, F. Damas, E. De Moraes Gregores, B. Lopes Da Costa, I. Maietto Silverio, P.G. Mercandante, S.F. Novaes, S. Padula, M. Pereira Coelho, V. Scheurer, T. Tomei

Institute for Nuclear Research and Nuclear Energy, Bulgarian Academy of Sciences, Sofia, Bulgaria

A. Aleksandrov, G. Antchev, P. Danev, R. Hadjiiska, P. Iaydjiev, M. Shopova, G. Sultanov

University of Sofia, Sofia, Bulgaria

A. Dimitrov, L. Litov, B. Pavlov, P. Petkov, A. Petrov

Instituto de Alta Investigación, Universidad de Tarapacá, Arica, Chile

S. Keshri, D.N. Laroze Navarrete, M. Meena, S. Thakur

Universidad Técnica Federico Santa María, Valparaiso, Chile

W. Brooks

Beihang University, Beijing, China

T. Cheng, L. Tan, L. Wang, L. Yuan

Department of Physics, Tsinghua University, Beijing, China

J. Gu, Z. Hu, Z. Liang, J. Liu, Y. Wang, H. Yang, S. Zhang, Y. Zhao

Institute of High Energy Physics, Beijing, China

N. Bi{sup}9, G.M. Chen{sup}9, H.S. Chen{sup}9, M. Chen{sup}9, Y. Chen, B. Hou{sup}9, Q. Hou, F. Iemmi, C.H. Jiang, P.z. Lai, H. Liao, G. Liu, Z. Liu{sup}10, S. Song{sup}9, J. Tao, C. Wang{sup}9, J. Wang, A. Zada{sup}9`, H. Zhang, J. Zhao

State Key Laboratory of Nuclear Physics and Technology, Peking University, Beijing, China

Y. Ban, A. Carvalho Antunes De Oliveira, S. Deng, X. Geng, B. Guo, Q. Guo, Z. He, C. Jiang, A. Levin, C. Li, L. Li, Q. Li, Y. Mao, S. Qian, S.J. Qian, X. Qin, C. Quaranta, X. Sun, D. Wang, J. Wang, T. Yang, M. Zhang, M. Zhang, Y. Zhao, C. Zhou

State Key Laboratory of Nuclear Physics and Technology, Institute of Quantum Matter, South China Normal University, Guangzhou, China, Guangzhou, China

X. Hua, S. Yang

Sun Yat-Sen University, Guangzhou, China

Z. You

University of Science and Technology of China, Hefei, China

N. Lu

Nanjing Normal University, Nanjing, China

G. Bauer{sup}11,12, L. Chen, Z. Cui{sup}12, B. Li{sup}13, H. Wang, X. Wang, K. Yi{sup}14, J. Zhang, F. Zhu

Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao, China

C. Li

Institute of Modern Physics and Key Laboratory of Nuclear Physics and Ion-beam Application (MOE) - Fudan University, Shanghai, China

Y. Li, Z. Wang, Y. Zhou{sup}15`

Zhejiang University - Department of Physics, Zhejiang, China

Z. Lin, C. Lu, M. Xiao{sup}16`

Universidad de Los Andes, Bogota, Colombia

C. Avila, A. Cabrera, C. Florez, J.A. Reyes Vega

Universidad de Antioquia, Medellin, Colombia C. Rendón, M. Rodriguez, A.A. Ruales Barbosa, J.D. Ruiz Alvarez

University of Split, Faculty of Electrical Engineering, Mechanical Engineering and Naval Architecture, Split, Croatia

N. Godinovic, D. Lelas, I. Puljak, A. Sculac

University of Split, Faculty of Science, Split, Croatia

M. Kovac, A. Petkovic, T. Sculac

Institute Rudjer Boskovic, Zagreb, Croatia

P. Bargassa, V. Brigljevic, D. Ferencek, K. Jakovcic, A. Starodumov, T. Susa

University of Cyprus, Nicosia, Cyprus

A. Attikis, S. Konstantinou, C. Leonidou, L. Paizanos, F. Ptochos, P.A. Razis, H. Saka, A. Stepennov

Charles University, Prague, Czech Republic

M. Finger Jr., A. Kveton

Escuela Politecnica Nacional, Quito, Ecuador

E. Acurio

Universidad San Francisco de Quito, Quito, Ecuador

E. Carrera Jarrin

Academy of Scientific Research and Technology of the Arab Republic of Egypt, Egyptian Network of High Energy Physics, Cairo, Egypt

H. Abdalla{sup}17, Y. Assran{sup}18,19, B. El-mahdy{sup}17`

Center for High Energy Physics (CHEP-FU), Fayoum University, El-Fayoum, Egypt

A. Lotfy, M. Mahmoud, H. Mohammed, Y. Mohammed

National Institute of Chemical Physics and Biophysics, Tallinn, Estonia

K. Jaffel, M. Kadastik, T. Lange, C. Nielsen, J. Pata, M. Raidal, N. Seeba, L. Tani

Department of Physics, University of Helsinki, Helsinki, Finland

E. Brü̈cken, A. Milieva, K. Osterberg, M. Voutilainen

Helsinki Institute of Physics, Helsinki, Finland

F.I. Garcia Fuentes, T. Hilden, P. Inkaew, K.T.S. Kallonen, R. Kumar Verma, T. Lampén, K. Lassila-Perini, B. Lehtela, S. Lehti, T. Lindén, N.R. Mancilla Xinto, M. Myllymäki, M.m. Rantanen, S. Saariokari, N.T. Toikka, J. Tuominiemi, E. Veikkola

Lappeenranta-Lahti University of Technology, Lappeenranta, Finland

N. Bin Norjoharuddeen, H. Kirschenmann, P.R. Luukka, H. Petrow

IRFU, CEA, Université Paris-Saclay, Gif-sur-Yvette, France

M. Besancon, F. Couderc, M. Dejardin, D. Denegri, P. Devouge, J.L. Faure, F. Ferri, P. Gaigne, S. Ganjour, P. Gras, F. Guilloux, G. Hamel de Monchenault, M. Kumar, V. Lohezic, Y. Maidannyk, J. Malcles, F. Orlandi, L. Portales, S. Ronchi, M.Ö. Sahin, P. Simkina, M. Titov

Laboratoire Leprince-Ringuet, CNRS/IN2P3, Ecole Polytechnique, Institut Polytechnique de Paris, Palaiseau, France

R. Amella Ranz, F. Beaudette, K. Biriukov, P. Busson, F. Cetorelli, C. Charlot, M. Chiusi, T.D. Cuisset, O. Davignon, A. De Wit, T. Debnath, I.T. Ehle, S. Ghosh, A. Gilbert, R. Granier de Cassagnac, M. Manoni, M. Nguyen, S. Obraztsov, C. Ochando, L.m. Rabour, R. Salerno, J.B. Sauvan, Y. Sirois, G. Sokmen, Y. Song, L. Urda Gómez, B. Voirin, A. Zabi, A. Zghiche

Institut Pluridisciplinaire Hubert Curien (IPHC), Université de Strasbourg, CNRS/IN2P3, Strasbourg, France

J.L. Agram20, J. Andrea, D. Bloch, E.C. Chabert, C. Collard, G. Coulon, C. Eschenlauer, S. Falke, U. Goerlach, A.C. Le Bihan, G. Saha, A. Savoy-Navarro21, P. Vaucelle

Institut de Physique des 2 Infinis de Lyon (IP2I), Villeurbanne, France

D. Amram, S. Beauceron, B. Blancon, G. Boudoul, N. Chanon, D. Contardo, J. Daniel, P. Depasse, H. El Mamouni, J. Fay, E. Fillaudeau, S. Gascon, M. Gouzevitch, C. Greenberg, B. Ille, E. Jourd’Huy, M. Lethuillier, K. Long, B. Massoteau, L. Mirabito, A. Purohit, M. Vander Donckt, C. Verollet

Centre de Calcul de l’Institut National de Physique Nucléaire et de Physique des Particules, CNRS/IN2P3, Villeurbanne, France

A. Di Florio, G. Mauceri, B. Orzari

Georgian Technical University, Tbilisi, Georgia

A. Khvedelidze22, I. Lomidze, Z. Tsamalaidze22

RWTH Aachen University, I. Physikalisches Institut, Aachen, Germany

K.F. Adamowicz, V. Botta, S. Consuegra Rodríguez, L. Feld, K. Klein, M. Lipinski, P. Nattland, V. Oppenländer, A. Pauls, D. Pérez Adán

RWTH Aachen University, III. Physikalisches Institut A, Aachen, Germany

C. Daumann, S. Diekmann, E. Ehlert, N. Eich, D. Eliseev, F. Engelke, J. Erdmann, M. Erdmann, M.Z. Farkas, B. Fischer, T. Hebbeker, K. Hoepfner, A. Jung, N. Kumar, F. Mausolf, M. Merschmeyer, A. Meyer, A. Pozdnyakov, H. Reithler, U. Sarkar, V. Sarkisovi, A. Schmidt, J.G. Schulz, C. Seth, A. Sharma, J.L. Spah, V. Vaulin, U. Willemsen, S. Zaleski, F.P. Zinn

RWTH Aachen University, III. Physikalisches Institut B, Aachen, Germany

M.R. Beckers, G. Flügge, N. Hoeflich, T. Kress, A. Nowack, O. Pooth, A. Stahl

University of Hamburg, Hamburg, Germany

A.R. Alves Andrade, M. Antonello, S. Bollweg, M. Bonanomi, L. Ebeling, K. El Morabit, Y. Fischer, M. Frahm, P.P. Gadow, E. Garutti, A. Grohsjean, A.A. Guvenli, J. Haller, D. Hundhausen, M. Jalalvandi, G. Kasieczka, P. Keicher, R. Klanner, W. Korcari, T. Kramer, C.c. Kuo, J. Lange, M.y. Lee, A. Lobanov, J. Matthiesen, L. Moureaux, K. Nikolopoulos, K.J. Pena Rodriguez, N. Prouvost, B. Raciti, M. Rieger, D. Savoiu, P. Schleper, M. Schröder, J. Schwandt, M. Sommerhalder, H. Stadie, G. Steinbrück, J. Sun, T. von Schwartz, R. Ward, B. Wiederspan, M. Wolf, C. Yede

Deutsches Elektronen-Synchrotron, Hamburg, Germany

A. Abel, A. Akhkil, M. Aldaya Martin, J. Alimena, Y. An, I. Andreyev, J. Bach, S. Baxter, H. Becerril Gonzalez, O. Behnke, F. Blekman23, K. Borras24, L. Braga Da Rosa, A. Campbell, C. Cazzaniga, S. Chatterjee, L.X. Coll Saravia, G. Eckerlin, D. Eckstein, E. Gallo23, A. Geiser, M. Guthoff, A. Hinzmann, U. Husemann, M. Kasemann, C. Kleinwort, R. Kogler, M. Komm, D. Krücker, F. Labe, W. Lange, D. Leyva Pernia, J.h. Li, K.y. Lin, K. Lipka25, W. Lohmann26, J. Malvaso, R. Mankel, I.A. Melzer-Pellmann, M. Mendizabal Morentin, A.B. Meyer, G. Milella, M.N.J. Momed, K. Moral Figueroa, A. Mussgiller, L.P. Nair, A. Nürnberg, J. Park, F. Preau, E. Ranken, A. Raspereza, D. Rastorguev, L. Rygaard, M. Scham27,24, C. Schwanenberger23, D. Schwarz, P. Schütze, D. Selivanova, K. Sharko, M. Shchedrolosiev, A. Sritharan, D. Stafford, M. Torkian, S. Vashishtha, R. Walsh, D. Wang, Q. Wang, K. Wichmann, C. Wissing, S. Zakharov, A. Zimmermane Castro Santos

Institut für Experimentelle Teilchenphysik, Karlsruhe, Germany

J. Ahäuser, A. Brusamolino, E. Butz, Y.M. Chen, T. Chwalek, A. Dierlamm, G.G. Dincer, U. Elicabuk, N. Faltermann, M. Giffels, A. Gottmann, F. Hartmann28, F. Hummer, J. Kieseler, M. Klute, H.A. Krause, R. Kunnilan Muhammed Rafeek, O. Lavoryk, J.M. Lawhorn, S. Maier, N. Meenamthuruthil Radhakrishnan, T. Mehner, M. Molch, A.A. Monsch, M. Mormile, T. Müller, M. Presilla, G. Quast, K. Rabbertz, B. Regnery, R. Schmieder, T. Selezneva, N. Shadskiy, L. Sowa, L. Stockmeier, M. Toms, B. Topko, N. Trevisani, C. Verstege, T. Voigtländer, R.F. Von Cube, J. Von Den Driesch, J.H. Voss, C. Winter, R. Wolf, W.D. Zeuner, X. Zuo

Institute of Nuclear and Particle Physics (INPP), NCSR Demokritos, Aghia Paraskevi, Greece

G. Anagnostou, G. Daskalakis, A. Kyriakis

National and Kapodistrian University of Athens, Athens, Greece

P. Iosifidou, P. Katris, M. Kotsarini, G. Melachroinos, Z. Painesıs, N. Plastiras, N. Saoulidou, K. Theofilatos, E. Tzovara, K. Vellidis, I. Zisopoulos

National Technical University of Athens, Athens, Greece

T. Chatzistavrou, G. Karapostoli, K. Kousouris, K. Paschos, L.P. Rouseliotaki, E. Siamarkou, A. Taxeidi, G. Tsipolitis

University of Ioánnina, Ioánnina, Greece

I. Evangelou, C. Foudas, P. Katsoulis, P. Kokkas, P.G. Kosmoglou Kioseoglou, N. Manthos, I. Papadopoulos, J. Strologas

Department of Physics, School of Sciences Democritus, University of Thrace, Kavala, Greece

E. Tziaferi

HUN-REN Wigner Research Centre for Physics, Budapest, Hungary

C. Hajdu, D. Horvath29,30, Á. Kadlecsik, C. Lee, K. Márton, A.J. Rád31, F. Sikler, V. Veszpremi

MTA-ELTE Lendület CMS Particle and Nuclear Physics Group, Eötvös Loránd University, Budapest, Hungary

G. Balint, M. Csanad, K. Farkas, A. Fehérkuti32, M.M.A. Gadallah33, M. León Coello, G. Pasztor, G.I. Veres

Faculty of Informatics, University of Debrecen, Debrecen, Hungary, Debrecen, Hungary

B. Ujvari, G. Zilizi

HUN-REN ATOMKI - Institute of Nuclear Research, Debrecen, Hungary

G. Bencze, S. Czellar, J. Molnar, Z. Szillasi

Karoly Robert Campus, MATE Institute of Technology, Gyongyos, Hungary

T.F. Csorgo32, F. Nemes32, T. Novak, I. Szanyi34 Indian Institute of Science (IISC), Bangalore, India J.R. Komaragiri

Indian Institute of Technology Bhubaneswar, Bhubaneswar, India S. Bahinipati, R. Raturi

Panjab University, Chandigarh, India S. Bansal, V. Bhatnagar, B. Chauhan, S. Chauhan, N. Dhingra35, A. Kaur, H. Kaur, S. Kumar, T. Sheokand, A. Singla, K. Verma

University of Delhi, Delhi, India A. Bhardwaj, A. Chhetri, B.C. Choudhary, A. Kumar, A. Kumar, M. Naimuddin, S. Phor, C. Prakash, K. Ranjan, M.K. Saini

Indian Institute of Technology Mandi (IIT-Mandi), Himachal Pradesh, India M. Kumari, N. Neeraj, P. Palni, S. Rana, A. Rathore, A. Sarkar

University of Hyderabad, Hyderabad, India S. Acharya36, B. Gomber, S.K. Satapathy

Indian Institute of Technology Kanpur, Kanpur, India S. Mukherjee

Saha Institute of Nuclear Physics, HBNI, Kolkata, India S. Bhattacharya, S. Das Gupta, S. Dutta, S. Dutta, S. Sarkar

Indian Institute of Technology Madras, Madras, India M.M. Ameen, P.K. Behera, S. Chatterjee, G. Dash, A. Dattamunsi, P. Jana, P. Kalbhor, S. Kamble, P.R. Pujahari, A.K. Sikdar, R.K. Singh, A. Swain, P. Verma, S. Verma, A. Vijay

Indian Institute of Science Education and Research Mohali, Mohali, India A. Chauhan, S. Nayak, H. Rajpoot, B.K. Sirasva

Tata Institute of Fundamental Research-A, Mumbai, India L. Bhatt, S. Dugad, T. Mishra, G.B. Mohanty, M. Shelake, P. Suryadevara

Tata Institute of Fundamental Research-B, Mumbai, India A. Bala, S. Banerjee, S. Barman37, R.M. Chatterjee, J. Chhikara, M. Guchait, S. Jain, A. Jaiswal, S. Kumar, M. Maity37, G. Majumder, K. Mazumdar, L. Panwar38, R. Pramanik, R. Saxena, P. Sharma, A. Thachayath

National Institute of Science Education and Research, Jatni, Khorda, Odisha 752050, India Homi Bhabha National Institute, Training School Complex, Anushakti Nagar, Mumbai 400094, India, Odisha, India R. Kumar Agrawal, D. Maity39, P. Mal, K. Naskar39, A. Nayak39, K. Pal, P. Sadangi, S. Shuchi, S.K. Swain, S. Varghese39, D. Vats39

Indian Institute of Science Education and Research (IISER), Pune, India S. Dube, P. Hazarika, A. Laha, R. Sharma, S. Sharma, K.Y. Vaish

Indian Institute of Technology Hyderabad, Telangana, India C. Agrawal, B. Babu, S. Ghosh

Isfahan University of Technology, Isfahan, Iran H. Bakhshiansohi40, A. Jafari41, V. Sedighzadeh Dalavi, M. Zeinali42 Institute for Research in Fundamental Sciences (IPM), Tehran, Iran S. Bashiri, S. Chenarani43, S.M. Etesami, Y. Hosseini, M. Khakzad, E. Khazaie, M. Mohammadi Najafabadi, M. Nourbakhsh, S. Tizchang44

University College Dublin, Dublin, Ireland M. Felcini, M. Grunewald

INFN Sezione di Baria^{a}, Università di Barib^{b}, Politecnico di Baric^{c}, Bari, Italy M. Abbresciaa,b^{a,b}, M. Buonsantea,b^{a,b}, A. Colaleoa,b^{a,b}, D. Creanzaa,c^{a,c}, N. De Filippisa,c^{a,c}, M. De Palmaa,b^{a,b}, W. Elmetenaweea,b^{a,b}45, N. Ferraraa,c^{a,c}, L. Fiorea^{a}, L. Generosoa,b^{a,b}, L. Longoa^{a}, M. Loukaa,b^{a,b}, G. Maggia,c^{a,c}, M. Maggia^{a}, S. Mya,b^{a,b}, F. Nennaa,b^{a,b}, S. Nuzzoa,b^{a,b}, A. Pellecchiaa,b^{a,b}, A. Pompilia,b^{a,b}, F.M. Procaccia,b^{a,b}, G. Pugliesea,c^{a,c}, R. Radognaa,b^{a,b}, D. Ramosa^{a}, A. Ranieria^{a}, L. Silvestrisa^{a}, F.M. Simonea,b^{a,b}, A. Stamerraa,b^{a,b}, Ü. Sözbilira^{a}46, F. Tenchinia,b^{a,b}, D. Troianoa,b^{a,b}, R. Vendittia,b^{a,b}, P. Verwilligena^{a}, A. Zazaa,b^{a,b}

INFN Sezione di Bolognaa^{a}, Università di Bolognab^{b}, Bologna, Italy G. Abbiendia^{a}, S. Balduccia,b^{a,b}, C. Battilanaa,b^{a,b}, D. Bonacorsia,b^{a,b}, P. Capiluppia,b^{a,b}, M. Cruciania,b^{a,b}, M. Cuffiania,b^{a,b}, G.M. Dallavallea^{a}, T. Diotalevia,b^{a,b}, F. Fabbria^{a}, A. Fanfania,b^{a,b}, R. Farinellia^{a}, D. Fasanellaa^{a}, L. Ferraginaa,b^{a,b}, P. Giacomellia^{a}, C. Grandia^{a}, L. Guiduccia,b^{a,b}, M. Lorussoa,b^{a,b}, L. Lunertia^{a}, S. Marcellinia^{a}, G. Masettia^{a}, F. Navarriaa,b^{a,b}, G. Paggia,b^{a,b}, A. Perrottaa^{a}, A. Rossia,b^{a,b}, S. Rossi Tisbenia,b^{a,b}, T. Rovellia,b^{a,b}

INFN Sezione di Cataniaa^{a}, Università di Cataniab^{b}, Catania, Italy S. Costaa,b^{a,b}47, A. Di Mattiaa^{a}, A. Lapertosaa^{a}, R. Potenzaa,b^{a,b}, A. Tricomia,b47

INFN Sezione di Firenzea^{a}, Università di Firenzeb^{b}, Firenze, Italy J. Altorka,b^{a,b}, G. Barbaglia^{a}, A. Calandraa,b^{a,b}, B. Camaiania,b^{a,b}, A. Cassesea^{a}, R. Ceccarellia^{a}, V. Ciullia,b^{a,b}, C. Civinina^{a}, R. D’Alessandroa,b^{a,b}, L. Damentia,b^{a,b}, E. Focardia,b^{a,b}, T. Kelloa^{a}, G. Latinoa,b^{a,b}, P. Lenzia,b^{a,b}, M. Lizzoa^{a}, M. Meschinia^{a}, S. Paolettia^{a}, A. Papanastassioua,b^{a,b}, G. Sguazzonia^{a}, L. Viliania^{a}

INFN Laboratori Nazionali di Frascati, Frascati, Italy L. Benussi, S. Colafranceschi48, S. Meola49, D. Piccolo

INFN Sezione di Genovaa^{a}, Università di Genovab^{b}, Genova, Italy M. Alves Gallo Pereiraa^{a}, F. Ferroa^{a}, E. Robuttia^{a}, S. Tosia,b^{a,b}

INFN Sezione di Milano-Bicoccaa^{a}, Università di Milano-Bicocca, Milanob^{b}, Milano-Bicocca, Italy A. Benagliaa^{a}, F. Brivioa^{a}, V. Camagnia,b^{a,b}, F. De Guioa,b^{a,b}, M.E. Dinardoa,b^{a,b}, P. Dinia^{a}, S. Gennaiaa^{a}, R. Gerosaa,b^{a,b}, A. Ghezziaa,b^{a,b}, P. Govonia,b^{a,b}, L. Guzzia^{a}, G. Lavizzaria,b, M.T. Lucchinia,b^{a,b}, M. Malbertia, S. Malvezzia, A. Massironia, L. Moronia, M. Paganonia,b, S. Palluottoa,b, D. Pedrinia, A. Peregoa,b, T. Tabarelli de Fatisa,b^{a,b}

INFN Sezione di Napolia^{a}, Università di Napoli “Federico II”b^{b}, Università della Basilicata (Potenza)c^{c}, Scuola Superiore Meridionale (SSM)d^{d}, Napoli, Italy S. Buontempoa^{a}, F. Confortinia,b^{a,b}, C. Di Fraiaa,b^{a,b}, F. Fabozzia,c^{a,c}, A.O.M. Iorioa,b^{a,b}, L. Listaa,b^{a,b}50, P. Paoluccia^{a}28, B. Rossia^{a}

INFN Sezione di Padovaa^{a}, Università di Padovab^{b}, Università degli Studi di Cagliaric^{c}, Padova, Italy P. Azzia^{a}51, N. Bacchettaa^{a}, D. Biselloa,b^{a,b}, L. Borellaa^{a}, P. Bortignona,c^{a,c} G. Bortolatoa,b, A.C.M. Bullaa,c, R. Carlina,b, P. Checchiaa, T. Dorigoa,52, F. Gasparinia,b, U. Gasparinia,b, P. Gruttaa, N. Laia, E. Lusiania, M. Margonia,b, A.T. Meneguzzoa,b, M. Missirolia, J. Pazzinia,b, F. Primaveraa,b, P. Ronchesea,b, R. Rossina,b, F. Simonettoa,b, M. Toffanoa, M. Tosia,b, A. Triossia,b, S. Venturaa, M. Zanettia,b, P. Zottoa,b, A. Zucchettaa,b, G. Zumerlea,b

INFN Sezione di Paviaa, Università di Paviab, Pavia, Italy

S. A. AbuZeida,53, C. Aimèa, A. Braghieria, M. Brunoldia,b, P. Montagnaa,b, M. Pelliccionia,b, V. Rea, C. Riccardia,b, P. Salvinia, I. Vaiaa,b, P. Vituloa,b

INFN Sezione di Perugiaa, Università di Perugiab, Perugia, Italy

S. Ajmala,b, M.E. Asciotia,b, G.M. Bileia, W.D. Buitrago Ceballosa,b, C. Carrivalea,b, D. Ciangottinia,b, L. Della Pennaa,b, L. Fanòa,b, V. Mariania,b, M. Menichellia, F. Moscatellia,54, F. Napolitanoa, A. Rossia,b, A. Santocchiaa,b, D. Spigaa, T. Tedeschia,b

INFN Sezione di Pisaa, Università di Pisab, Scuola Normale Superiore di Pisac, Università di Sienad, Pisa, Italy

C.A. Alexea,c, P. Asenovaa,b, P. Azzurria, G. Bagliesia, L. Bianchinia,b, T. Boccalaia, E. Bossinia, D. Bruschinia,c, R. Castaldia, F. Cattafestaa,c, M.A. Ciocciad, M. Cipriania,b, R. Dell’Orsoa, S. Dhaniaa,d, S. Donatoa,b, A. Feliziania,d, R. Fortia,b, A. Giassia, F. Ligabuea,c, A.C. Marinia,b, A. Messineoa,b, S. Mishraa, V.K. Muraleedharan Nair Bindhua,b, S. Nandanaa, F. Pallaa, M. Riggirelloa,c, A. Rizzia,b, G. Rolandia,c, A. Scribanoa, P. Solankia,b, P. Spagnoloa, R. Tenchinia, G. Tonellia,b, N. Turinia,d, F. Vasellia,c, A. Venturia, P.G. Verdinia

INFN Sezione di Romaa, Sapienza Università di Romab, Roma, Italy

P. Akrapa,b, S.C. Beheraa, F. Cavallaria, L. Cunqueiro Mendeza,b, F. De Riggia,b, D. Del Rea,b, M. Del Vecchioa,b, E. Di Marcoa, M. Diemoza, F. Erricoa, L. Frosinaa,b, R. Gargiuloa,b, B. Harikrishnana,b, F. Lombardia,b, L. Martikainena,b, G. Organtinia,b, N. Palmeria,b, R. Paramattia,b, T. Paulettoa,b, S. Rahatloua,b, C. Rovellia, F. Santanastasioa,b, L. Soffia, V. Vladimirova,b

INFN Sezione di Torinoa, Università di Torinob, Università del Piemonte Orientale (Novara)c, Torino, Italy

N. Amapanea,b, R. Arcidiaconoa,c, S. Argiroa,b, M. Arneodoa,c, N. Bartosika,c, F. Bashira,b, R. Bellana,b, A. Belloraa,b, C. Biinoa, C. Borcaa,b, L. Bulajaa,b, N. Cartigliaa, M. Costaa,b, R. Covarellia,b, N. Demariaa, E. Ferrandoa,b, L. Fincoa, M. Grippoa,b, B. Kiania,b, L. Lanteria,b, F. Luongoa,b, C. Mariottia,55, S. Masellia, A. Meccaa,b, L. Menzioa,b, P. Meridiania, E. Migliorea,b, M. Montenoa, M.M. Obertinoa,b, G. Ortonaa, L. Pachera,b, N. Pastronea, M. Ruspaa,c, F. Sivieroa,b, V. Solaa,b, A. Solanoa,b, A. Staianoa, C. Tarriconea,b, M. Tornagoa,b, D. Trocinoa, G. Umoreta,b, E. Vlasovb, R. Whitea,b

INFN Sezione di Triestea, Università di Triesteb, Trieste, Italy

J. Babbara,b,56, S. Belfortea, V. Candelisea,b, M. Casarsaa, F. Cossuttia, K. De Leoa, G. Della Riccaa,b, R. Delli Gattia,b, C. Giraldina,b

Kyungpook National University, Daegu, Korea

S. Dogra, J. Hong, J. Kim, J. Kim, T. Kim, D. Lee, H. Lee, J. Lee, S.W. Lee, C.S. Moon, Y.D. Oh, S. Sekmen, B. Tae, Y.C. Yang

Department of Mathematics and Physics - Gangneung-Wonju National University, Gangneung, Korea

M.S. Kim

Chonnam National University, Institute for Universe and Elementary Particles, Kwangju, Korea

G. Bak, P. Gwak, H. Kim, H. Lee, S. Lee, D.H. Moon, J. Seo

Department of Physics, Chung-Ang University, Seoul, Korea

K. Lee, Y. Lee

Hanyang University, Seoul, Korea

E. Asilar, F. Carnevali, J. Choi57, T.J. Kim, Y. Ryou, J. Song, T. Yang

Korea University, Seoul, Korea

S. Ha, B.S. Hong, J. Kim, K. Lee, K. Lee, S. Lee, J. Padmanaban, B.A.N. Putra, J. Yoo

Kyung Hee University, Department of Physics, Seoul, Korea

J. Goh, J. Shin, S. Yang

Sejong University, Seoul, Korea

L. Kalipoliti, Y. Kang, H. Kim, Y. Kim, B. Ko, S. Lee

Seoul National University, Seoul, Korea

J. Choi, J. Choi, W. Jun, H. Kim, J. Kim, J. Kim, T. Kim, Y. Kim, Y.W. Kim, S. Ko, H. Lee, J. Lee, J. Lee, B.H. Oh, J. Shin, U. Yang, I. Yoon

University of Seoul, Seoul, Korea

W. Heo, W. Jang, D. Kim, S. Kim, Y. Roh, I. J. Watson

Yonsei University, Department of Physics, Seoul, Korea

S. Calzaferri, G. Cho, Y. Eo, K. Hwang, H. Jang, B. Kim, D. Kim, S. Kim, J.S.H. Lee, G. Mocellin, H.D. Yoo

Sungkyunkwan University, Suwon, Korea

Y. Lee, I. Yu

College of Engineering and Technology, American University of the Middle East (AUM), Dasman, Kuwait

T. Beyrouthy, Y. Gharbia

Kuwait University - College of Science - Department of Physics, Safat, Kuwait

F. Alazemi

Riga Technical University, Riga, Latvia

K. Dreimanis, O.M. Eberlins, A. Gaile, J.K. Heikkilä, M. Klevs, C. Munoz Diaz, D. Osite, G. Pikurs, R. Plese, M. Seidel, D. Sidiropoulos Kontos

University of Latvia (LU), Riga, Latvia

N.R. Strautnieks

Vilnius University, Vilnius, Lithuania

M. Ambrozas, A. Juodagalvis, S. Nargelas, S. Nayak, G. Tamulaitis

National Centre for Particle Physics, Universiti Malaya, Kuala Lumpur, Malaysia

I. Yusuff58, Z. Zolkapli

University of Sonora (UNISON), Hermosillo, Mexico J.P. Barajas Ibarria, J.F. Benitez, A. Castaneda Hernandez, A. Cota Rodriguez, L.E. Cuevas Picos, H.A. Encinas Acosta, L.G. Gallegos Mariñez, J.A. Murillo Quijada, L. Valencia Palomo

Centro de Investigacion y de Estudios Avanzados del IPN, Mexico City, Mexico

H. Castilla-Valdez, H. Crotte Ledesma, R. Lopez-Fernandez, J. Mejia Guisao, R. Reyes-Almanza, A. Sánchez Hernández

Universidad Iberoamericana, Mexico City, Mexico

C. Oropeza Barrera, D.L. Ramirez Guadarrama, M. Ramírez García

Benemerita Universidad Autonoma de Puebla, Puebla, Mexico

I. Bautista, F.E. Neri Huerta, I. Pedraza, H.A. Salazar Ibarguen, C. Uribe Estrada

University of Montenegro, Podgorica, Montenegro

I. Bubanja, J. Mijuskovic, N. Raicevic

National Centre for Physics, Quaid-I-Azam University, Islamabad, Pakistan

A. Ahmad, M.I. Asghar, A. Awais, M.I.M. Awan, W.A. Khan, I. Sohail

AGH University of Krakow, Krakow, Poland

Z. Abdy, V. Avati, L. Forthomme, L. Grzanka, M. Malawski, K. Piotrzkowski

National Centre for Nuclear Research, Swierk, Poland

H. Awedikian, M. Bluj, M. Ghimiray, M. Górski, M. Kazana, M. Szleper, P. Zalewski

Institute of Experimental Physics, Faculty of Physics, University of Warsaw, Warsaw, Poland

K. Bunkowski, K. Doroba, A. Kalinowski, M. Konecki, J. Krolikowski, W. Matyszkiewicz, A. Muhammad, S. Slawinski

Warsaw University of Technology, Warsaw, Poland

P. Fokow, K. Pozniak, W. Zabolotny

Laboratório de Instrumentação e Física Experimental de Partículas, Lisboa, Portugal

M. Araujo, C. Beirão Da Cruz E Silva, A. Boletti, M. Bozzo, T. Camporesi55,59{}^{55,59}

G. Da Molin, M. Gallinaro, R. Guitton, J. Hollar, H. Legoinha, N. Leonardo60{}^{60}, G.B. Marozzo, A. Petrilli, M. Pisano, J. Seixas, J. Varela, J.W. Wulff

Faculty of Physics, University of Belgrade, Belgrade, Serbia

P. Adzic, L. Markovic, P. Milenovic, V. Milosevic

Vinca Institute of Nuclear Science, Belgrade, Serbia

D. Devetak, M. Dordevic, J. Milosevic, L. Nadderd, V. Rekovic, M. Stojanovic

Centro de Investigaciones Energéticas Medioambientales y Tecnológicas (CIEMAT), Madrid, Spain

M. Alcalde Martinez, J. Alcaraz Maestre, J.A. Brochero Cifuentes, M. Cepeda, M. Cerrada, N. Colino, B. De La Cruz, A. Escalante Del Valle, C. Fernandez Bedoya, D. Fernández Del Val, J.P. Fernández Ramos, J. Flix, M.C. Fouz, C. Garcia Sanchez, M. Gonzalez Hernandez, O. Gonzalez Lopez, S. Goy Lopez, J.M. Hernandez, M.I. Josa, J. Llorente Merino, O. Manzanilla, C. Martin Perez, E. Martin Viscasillas, D. Moran, C.M. Morcillo Perez, Á. Navarro Tobar, J. Puerta Pelayo, A.M. Pérez-Calero Yzquierdo, I. Redondo, D.D. Redondo Ferrero, E. Sanchez Berenguer, J. Vazquez Escobar

Universidad Autónoma de Madrid, Madrid, Spain

J.F. de Troconiz

Universidad de Oviedo, Instituto Universitario de Ciencias y Tecnologías Espaciales de Asturias (ICTEA), Oviedo, Spain

E. Aller Gutierrez, B. Alvarez Gonzalez, J. Ayllon Torresano, A. Cardini, J. Cuevas, J. Del Riego Badas, D. Estrada Acevedo, J. Fernandez Menendez, S. Folgueras, L. Garcia Diaz, I. Gonzalez Caballero, P. Leguina, M. Obeso Menendez, E. Palencia Cortezon, J. Prado Pico, S. Sanchez Cruz, A. Soto Rodríguez, P. Vischia

Instituto de Física de Cantabria (IFCA), CSIC-Universidad de Cantabria, Santander, Spain

S. Blanco Fernández, I.J. Cabrillo, A. Calderon, M. Caserta, J. Duarte Campderros, M. Fernandez, G. Gomez, A. Gomez Carrera, C. Lasaosa García, R. Lopez Ruiz, C. Martinez Rivero, P. Martinez Ruiz del Arbol, F. Matorras, P. Matorras Cuevas, E. Navarrete Ramos, J. Piedra Gomez, C. Quintana San Emeterio, V. Rodriguez, L. Scodellaro, I. Vila, R. Vilar Cortabitarte, J.M. Vizan Garcia

University of Colombo, Colombo, Sri Lanka

B. Kailasapathy61

University of Ruhuna, Department of Physics, Matara, Sri Lanka

W.G. Dharmaratna62, N. Perera

CERN, European Organization for Nuclear Research, Geneva, Switzerland

D. Abbaneo, C. Amendola, R. Ardino, E. Auffray, J. Baechler, G. Bardelli, D. Barney, J. Bendavid, I. Bestintzanos, M. Bianco, A. Bocci, G. Boldrini, L. Borgonovi, C. Botta, A. Bragagnolo, C.E. Brown, C. Caillol, G. Cerminara, P. Connor, K. Cormier, D. D’Enterria, A. Dabrowski, P. Das, A. David Tinoco Mendes, M.M. Defranchis, M. Deile, M. Dobson, L. Favilla, P.J. Fernández Manteca, E. Fialova, B.A. Fontana Santos Alves, E. Fontanesi, W. Funk, A. Gaddi, S. Giani, D. Gigi, K. Gill, S. Giorgetti, F. Glege, M. Glowacki, A. Gruber, J. Hegeman, R. Hofsaess, B. Huber, T. James, P. Janot, L. Jeppe, O. Kaluzinska, O. Karacheban26, G. Karathanasis, S. Laurila, P. Lecoq, E. Leutgeb, J. Leon Holgado, C. Lourenco, A.m. Lyon, M. Magherini, L. Malgeri, E. Manca, M. Mannelli, F. Meijers, J.A. Merlin, S. Mersi, E. Meschi, M. Migliorini, F. Monti, F. Moortgat, M.C. Muehlnikel, M. Mulders, M. Musich, I. Neutelings, S. Orfanelli, F. Pantaleo, M. Pari, F. Pereira Carneiro, G. Petrucciani, A. Pfeiffer, M. Pierini, M. Pitt, H. Qu, W. Redjeb, A. Reimers, B. Ribeiro Lopes, F. Riti, P. Rosado, M. Rovere, H. Sakulin, R. Salvatico, S. Scarfi, S.F. Schaefer, M. Selvaggi, P. Silva, P. Sphicas63, A.G. Stahl Leiton, A. Steen, S. Summers, G. Terragni, D. Treille, P. Tropea, E. Vernazza, M. Vojinovic, J. Wanczyk64, S. Wuchterl, M. Zarucki, P. Zehetner, P. Zejdl, G. Zevi Della Porta

PSI Center for Neutron and Muon Sciences, Villigen, Switzerland

L. Caminada65, W. Erdmann, R. Horisberger, Q. Ingram, H.C. Kaestli, D. Kotlinski, C. Lange, U. Langenegger, A. Nigamova, L. Noehte65, L. Redard-Jacot65, T. Rohe, A. Samalan

ETH Zurich - Institute for Particle Physics and Astrophysics (IPA), Zurich, Switzerland

T.K. Aarrestad, M. Backhaus, A. Belvedere, T. Bevilacqua65, G. Bonnemelli, K. Datta, P. De Bryas Dexmiers D’Archiacchiae64, A. De Cosa, G. Dissertori, M. Dittmar, M. Donegà, F. Glessgen, C. Grab, T.G. Harte, N. Härringer, B. Kaynak, M. Koppel, W. Lustermann, M. Malucchi, R.A. Manzoni, L. Marchese, F. Nessi-Tedaldi, F. Pauss, A.A. Petre, J. Prendi, B. Ristic, S. Rohletter, P.M. Sander, R. Seidita, A. Tarabini, C.Z. Tee, D. Valsecchi, P.H. Wagner, R. Wallny

Universität Zürich, Zurich, Switzerland

C. Amsler66, F. Bilandzija, P. Bärtschi, M.F. Canelli, G. Celotto, Z. Ghafoor, T.A. Goldschmidt, V. Guglielmi, A. Jofrehei, B. Kilminster, T.H. Kwok, S. Leontsinis, V. Lukashenko, A. Macchiolo, F. Meng, J. Motta, P. Robmann, E. Shokr, F. Stäger, R. Tramontano, P. Viscone

Çukurova University, Adana, Türkiye

D. Agyel, F. Dolek, I. Dumanoglu67, Y. Guler68, E. Gurpinar Guler68, A. Kayis Topaksu, G. Onengut, K. Ozdemir69, B. Tali70, U.G. Tok, E. Uslan

Hacettepe University, Ankara, Türkiye

S. Sen

Bogazici University, Istanbul, Türkiye

B. Akgun, I.O. Atakisi71, E. Gulmez, M. Kaya72, O. Kaya73, M.A. Sarkisla74, S. Tekten75

Istanbul Technical University, Istanbul, Türkiye

D. Boncukcu, A. Cakir, K. Cankocak67,76, M. Gumustekin, A.D. Gungordu

Istanbul University, Istanbul, Türkiye

B. Hacisahinoglu, I. Hos77, S. Ozkorucuklu, O. Potok, H. Sert, C. Simsek, C. Zorbilmez

Yildiz Technical University, Istanbul, Türkiye

S. Cerci, C. Dozen78, E. Iren79, B. Isildak, E. Simsek, D. Sunar Cerci, T. Yetkin78

National Central University, Chung-Li, Taiwan

D. Bhowmik, Y.h. Chou, C.M. Kuo, P.K. Rout, S. Taj, P.C. Tiwari38

National Taiwan University (NTU), Taipei, Taiwan

L. Ceard, K.F. Chen, Z.g. Chen, A. De Iorio, G.W.S. Hou, H.w. Hsia, T.h. Hsu, S. Karmakar, F. Khuzaimah, G. Kole, Y.y. Li, R.S. Lu, E. Paganis, X.f. Su, J. Thomas-Wilsker, L.s. Tsai, D. Tsionou, H.y. Wu, E. Yazgan

High Energy Physics Research Unit, Department of Physics, Faculty of Science, Chulalongkorn University, Bangkok, Thailand

C. Asawatangtrakuldee, N. Srimanobhas

Tunis El Manar University, Tunis, Tunisia

Y. Maghrbi

Institute for Scintillation Materials of National Academy of Science of Ukraine, Kharkiv, Ukraine

O. Dadazhanova, B. Grynyov

National Science Centre, Kharkiv Institute of Physics and Technology, Kharkiv, Ukraine

K. Klimenko, O. Kurov, L. Levchuk, S. Lukyanenko, A. Pristavka, D. Soroka

University of Bristol, Bristol, United Kingdom

J.J. Brooke, A. Bundock, F.J.J. Bury, E. Clement, D. Cussans, D. Dharmender, H. Flacher, J. Goldstein, H.F. Heath, M.I. Holmberg, A. Karakoulaki, L. Kreczko, S. Paramesvaran, L. Robertshaw, M.S. Sanjrani40, J. Segal, V.J. Smith Rutherford Appleton Laboratory, Didcot, United Kingdom

A. Ball, K.W. Bell, A. Belyaev80, C. Brew, R.M. Brown, D.J. Cockerill, A. Elliot, K.V. Ellis, J. Gajownik, K. Harder, S. Harper, J. Linacre, K. Manolopoulos, M. Moallemi, D.M. Newbold, E. Olaiya, D. Petyt, T. Reis, A.R. Sahasransu, T. Schuh, C. Shepherd-Themistocleous, I.R. Tomalin, K.C. Whalen, T. Williams

Imperial College, London, United Kingdom

I. Andreou, S. Awan, R. Bainbridge, P. Bloch, O. Buchmuller, C.A. Carrillo Montoya, D. Colling, A. Cox, I. Das, P. Dauncey, G. Davies, A. De Roeck, M. Della Negra, S. Fayer, G. Fedi, G. Hall, H.R. Hoorani, A. Howard, G. Iles, C.R. Knight, P. Krueper, J. Langford, K.H. Law, L. Lyons, A.M. Magnan, B. Maier, S. Mallios, A. Mastronikolis, J. Nash81, M. Pesaresi, P.B. Pradeep, E.V. Protopapa, B.C. Radburn-Smith, A. Richards, A. Rose, T.B. Runting, L. Russell, K. Savva, R. Schmitz, C. Seez, R. Shukla, A. Tapper, T. Travis, K. Uchida, G.P. Uttley, T. Virdee28, N. Wardle, D. Winterbottom, J. Xiao

Brunel University, Uxbridge, United Kingdom

J. Cole, L. Juckett, A. Khan, P. Kyberd, I. Reid

The University of Alabama, Tuscaloosa, Alabama, USA

B. Bam, A. Buchot Perraguin, S. Campbell, R. Chudasama, S. Cooper, C. Crovella, G. Fidalgo, S.V. Gleyzer, R. Kaur, A. Khukhunaishvili, K. Matchev, E. Pearson, P. Rumerio82, E. Usai

University of California, Davis, Davis, California, USA

S. Abbott, S. Baradia, B. Barton, R. Breedon, H. Cai, M. Calderon De La Barca Sanchez, E. Cannaert, M. Chertok, M. Citron, J. Conway, P.T. Cox, F. Eble, R. Erbacher, C. Fairchild, O. Kukral, S. Ostrom, I. Salazar Segovia, J.H. Steenis, J.S. Tafoya Vargas, W. Wei, S. Yoo

University of California, San Diego, La Jolla, California, USA

A. Aportela, A. Arora, J.G. Branson, S. Cittolin, B. D’Anzi, D. Diaz, J. Duarte, L. Giannini, Y. Gu, J. Guiang, V. Krutelyov, R. Lee, J. Letts, H. Li, R. Marroquin Solares, M. Masciovecchio, F. Mokhtar, S. Morovic, S. Mukherjee, M. Pieri, D. Primosch, M. Quinnan, V. Sharma, M. Tadel, E. Vourliotis, F. Würthwein, A. Yagil, Z. Zhao

University of California, Los Angeles, California, USA

K. Adamidis, H. Ancelin, M. Bachtis, D. Campos, R. Cousins, S. Crossley, G. Flores Avila, J. Hauser, M. Ignatenko, M.A. Iqbal, T. Lam, Y.f. Lo, A. Nunez Del Prado, D. Saltzberg, V. Valuev

California Institute of Technology, Pasadena, California, USA

A. Albert, S. Bhattacharya, A. Bornheim, O. Cerri, Z. Hao, R. Kansal, L. Mori, H.B. Newman, G. Reales Gutiérrez, T. Sievert, P. Simmerling, E. Sledge, M. Spiropulu, C. Sun, J.R. Vlimant, R.A. Wynne, S. Xie, R.Y. Zhu

University of California, Riverside, Riverside, California, USA

R. Clare, J.W. Gary, G. Hanson

University of California, Santa Barbara - Department of Physics, Santa Barbara, California, USA

A. Barzdukas, L. Brennan, C. Campagnari, S. Carron Montero83, K. Downham, C. Grieco, J.S. Guo, M.M. Hussain, D. Imani, J. Incandela, A. Krishna, M.W.K. Lai, P. Masterson, J.J.H. Ockenfuss, J. Richman, S.N. Santpur, D. Stuart, T.Á. Vámi, X. Yan, D. Zhang

University of Colorado Boulder, Boulder, Colorado, USA

J.P. Cumalat, W.T. Ford, J. Fraticelli, A. Hart, M. Herrmann, S. Kwan, J. Pearkes, N. Schonbeck, K. Stenson, K. Ulmer, S.R. Wagner, N. Zipper, D. Zuolo

The Catholic University of America, Washington, DC, USA

R. Bartek, A. Dominguez, S. Raj, B. Sahu, A.E. Simsek, B. Singhal, S.S. Yu

University of Florida, Gainesville, Florida, USA

C. Aruta, P. Avery, C. Basile, D. Bourilkov, P. Chang, V. Cherepanov, M. Dittrich, R.D. Field, C. Huh, E. Koenig, M. Kolosova, J. Konigsberg, A. Korytov, G. Mitselmakher, K. Mohrman, A. Muthirakalayil Madhu, N. Rawal, S. Rosenzweig, Y. Takahashi, J. Wang

Florida Institute of Technology, Melbourne, Florida, USA

B. Alsufyani, S. Das, S. Demarest, L. Hasa, M. Hohlmann, M. Lavinsky, E. Yanes

Florida State University, Tallahassee, Florida, USA

T. Adams, A. Al Kadhim, D. Alam, A. Askew, S. Bower, R. Goff, R. Hashmi, A. Hassani, T. Kolberg, G. Martinez, M. Mazza, H. Prosper, P.R. Prova, R. Yohay

Fermi National Accelerator Laboratory, Batavia, Illinois, USA

M. Albrow, M. Alyari, O. Amram, G. Apollinari, A. Apresyan, L.A. Bauerdick, D. Berry, J. Berryhill, P.C. Bhat, K. Burkett, J.N. Butler, A. Canepa, G.B. Cerati, H. Cheung, F. Chlebana, C. Cosby, G. Cummings, I. Dutta, V.D. Elvira, J. Freeman, A. Gandrakota, Z. Gecse, L. Gray, D. Green, A. Grummer, S. Grunendahl, D. Guerrero, O. Gutsche, R.M. Harris, J. Hirschauer, V. Innocente, B. Jayatilaka, S. Jindariani, M. Johnson, R.S. Kim, S. Lammel, D. Lincoln, R. Lipton, T. Liu, K. Maeshima, D. Mason, P. McBride, P. Merkel, S. Mrenna, S. Nahn, J. Ngadiuba, D. Noonan, S. Norberg, V. Papadimitriou, N. Pastika, K. Pedro, C. Pena⁸⁴, C.E. Perez Lara, V. Perovic, F. Ravera, A. Reinsvold Hall⁸⁵, L. Ristori, M. Safdari, E. Sexton-Kennedy, E. Smith, N. Smith, A. Soha, L. Spiegel, S. Stoynev, J. Strait, L. Taylor, S. Tkaczyk, N.V. Tran, L. Uplegger, E.W. Vaandering, C. Wang, I. Zoi

University of Illinois Chicago, Chicago, Illinois, USA

M.R. Adams, N. Barnett, A. Baty, C. Bennett, N. Brandman-hughes, R. Cavanaugh, S.J. Das, R. Escobar Franco, O. Evdokimov, C.E. Gerber, H. Gupta, M. Hawksworth, A. Hingrajiya, D.J. Hofman, Z. Huang, J.h. Lee, C. Mills, S. Nanda, G. Nigmatkulov, B. Ozek, V. Pant, T. Phan, D. Pilipovic, R. Pradhan, E. Prifti, T. Roy, D. Shekar, N. Singh, F. Strug, A. Thielen, M. Tonjes, N. Varelas, M.A. Wadud, A. Wang, J. Yoo

Northwestern University, Evanston, Illinois, USA

S. Dittmer, K.A. Hahn, S. King, D. Li, M. Mcginnis, Y. Miao, D.G. Monk, M.H. Schmitt, A. Taliercio, M. Velasco, J. Wang, D. Wilbern

Purdue University Northwest, Hammond, Indiana, USA

N. Parashar, A. Pathak, E. Shumka

University of Notre Dame, Notre Dame, Indiana, USA

G. Agarwal, R. Band, S. Castells, A. Das, A. Datta, A. Ehnis, R. Goldouzian, M. Hildreth, T. Ivanov, C. Jessop, K. Lannon, J. Lawrence, D. Lutton, J. Mariano, N. Marinelli, P. Mastrapasqua, A. Masud, T. McCauley, C. Mcgrady, C. Moore, Y. Musienko22, H. Nelson, M. Osherson, A. Piccinelli, R. Ruchti, A. Townsend, Y. Wan, M. Wayne, H. Yockey

Purdue University, West Lafayette, Indiana, USA

S. Chandra, A. Gu, L. Gutay, L. He, M. Huwiler, M. Jones, A.W. Jung, I.G. Karslioglu, D. Kondratyev, J. Li, M. Liu, M. Macedo, G. Negro, N. Neumeister, G. Paspalaki, S. Piperov, N.R. Saha, J.F. Schulte, R. Sharma, F. Wang, A.L. Wesolek, A. Wildridge, W. Xie, Y. Yao, Y. Zhong

The University of Iowa, Iowa City, Iowa, USA

M. Alhusseini, D. Blend, K. Dilsiz86, O.K. Köseyan, A. Mestvirishvili59, O. Neogi, H. Ogul87, Y. Onel, A. Penzo, C. Snyder

The University of Kansas, Lawrence, Kansas, USA

A. Abreu, L.F. Alcerro Alcerro, J. Anguiano, S. Arteaga Escatel, P. Baringer, A. Bean, R. Bhattacharya, M. Chukwuka, Z. Flowers, D. Grove, J. King, G. Krintiras, M. Lazarovits, C. Le Mahieu, J. Marquez, M. Murray, M. Nickel, E. Reynolds, C. Rogan, C. Royon, S. Rudrabhatla, S. Sanders, J.A. Velazquez Corral, G. Wilson

Kansas State University, Manhattan, Kansas, USA

A. Ahmad, B. Allmond, N. Islam, A. Ivanov, K. Kaadze, Y. Maravin, J. Natoli, G.G. Reddy, D. Roy, G. Sorrentino

Johns Hopkins University, Baltimore, Maryland, USA

B. Blumenfeld, J. Davis, A. Gritsan, Z. Huang, L. Kang, P. Maksimovic, N. Pinto, M. Roguljic, S. Sekhar, M.V. Srivastav, M. Swartz

University of Maryland, College Park, Maryland, USA

Z. Alton, D. Baden, A. Belloni, J. Bistany-riebman, S.C. Eno, N.J. Hadley, S. Jabeen, R.G. Kellogg, T. Koeth, B. Kronheim, J. Lee, P. Major, A. Mignerey, C. Palmer, C. Papageorgakis, M.M. Paranjpe, E. Popova88, A. Shevelev, M. Wrotny, L. Zhang

Boston University, Boston, Massachusetts, USA

S. Cholak, Z. Demiragli, C. Erice, C. Fangmeier, C. Fernandez Madrazo, J. Fulcher, J. Garcia De Castro, F. Golf, S. Jeon, G. Linney, J. O’Cain, I. Reed, J. Rohlf, D. Sperka, I. Suarez, A. Tsatsos, E. Wurtz, A.G. Zecchinelli

Northeastern University, Boston, Massachusetts, USA

A. Aarif, G. Alverson, E. Barberis, S. Bein, J. Bonilla, B. Bylsma, M. Campana, R. Clark, Y. Han, I. Israr, M. Lu, N. Manganelli, R. Mccarthy, D.M. Morse, T. Orimoto, L. Skinnari, C.S. Thoreson, E. Tsai, D. Wood

Massachusetts Institute of Technology, Cambridge, Massachusetts, USA

C. Baldenegro Barrera, H. Bossi, S. Bright-Thorne, I.A. Cali, Y.C. Chen, P.C. Chou, M. D’Alfonso, K. Devereaux, J. Eysmans, C. Freer, G. Gomez Ceballos, M. Goncharov, G. Grosso, P. Harris, D. Hoang, A. Holterman, G.M. Innocenti, K. Ivanov, G. Kopp, D. Kovalskyi, J. Lang, L. Lavezzo, Y.J. Lee, P. Lugato, C. Mcginn, E. Moreno, A. Novak, M.I. Park, C. Paus, C. Reissel, C. Roland, G. Roland, S. Rothman, T.a. Sheng, G. Stephans, D. Walter, J. Wang, Z. Wang, B. Wyslouch, T. Yang, K. Yoon Wayne State University, Detroit, Michigan, USA

P.E. Karchin

University of Minnesota, Minneapolis, Minnesota, USA

A. Alpana, B. Crossman, W.J. Jackson, C. Kapsiak, D. Mahon, J. Mans, B. Marzocchi, R. Rusack, O. Sancar, R. Saradhy, N. Strobbe

Bethel University, St. Paul, Minnesota, USA

J.M. Hogan

University of Nebraska-Lincoln, Lincoln, Nebraska, USA

K. Bloom, D.R. Claes, S.V. Dixit, G. Haza, J. Hossain, C. Joo, I. Kravchenko, K.H.M. Kwok, Y. Mehra, J. Morris, A. Rohilla, J.E. Siado, A. Vagnerini, A. Wightman

Rutgers, The State University of New Jersey, Piscataway, New Jersey, USA

B. Chiarito, J.P. Chou, S. Donnelly, D. Gadkari, Y. Gershtein, E. Halkiadakis, C. Houghton, D. Jaroslawski, A. Kaur, A. Kobert, A. Lath, J. Martins, P. Meltzer, M. Perez Prada, K. Ramdin, B. Rand, J. Reichert, P. Saha, S. Salur, S. Somalwar, R. Stone, S.A. Thayil, S. Thomas, J. Vora

Princeton University, Princeton, New Jersey, USA

H. Bouchamaoui, G. Dezoort, P. Elmer, A. Frankenthal, M. Galli, B. Greenberg, K. Kennedy, Y. Lai, D. Lange, A. Loeliger, D. Marlow, I. Ojalvo, J. Olsen, F. Simpson, D. Stickland, C. Tully, S. Yoon

State University of New York at Buffalo, Buffalo, New York, USA

H. Bandyopadhyay, I. Iashvili, A. Kalogeropoulos, A. Kharchilava, A. Mandal, C. McLean, D. Nguyen, O. Poncet, S. Rappoccio, H. Rejeb Sfar, W. Terrill, D. Yu

Cornell University, Ithaca, New York, USA

J. Alexander, X. Chen, G. De Castro, J. Dickinson, A. Duquette, J. Fan, X. Fan, J. Grassi, P. Kotamnives, K. Krzyzanska, J. Monroy, G. Niendorf, M. Oshiro, J.R. Patterson, A. Ryd, J. Thom, H.A. Weber, B. Weiss, P. Wittich, Y. Wu, R. Zou, L. Zygala

University of Rochester, Rochester, New York, USA

A. Bodek, R. Demina, A. Garcia-Bellido, H.S. Hare, O. Hindrichs, Y.w. Kao, N. Parmar, P. Parygin{sup}88`, H. Seo, R. Taus, Y.h. Yu

The Ohio State University, Columbus, Ohio, USA

M. Carrigan, R. De Los Santos, L.S. Durkin, C. Hill, M. Joyce, L. Nestor, D.A. Wenzl, B.L. Winer, B. Yates

Carnegie Mellon University, Pittsburgh, Pennsylvania, USA

J. Alison, S. An, M. Cremonesi, V. Dutta, E.Y. Ertorer, T. Ferguson, T.A. Gómez Espinosa, A. Harilal, A. Kallil Tharayil, M. Kanemura, A. Khanal, C. Liu, M. Marchegiani, P. Meiring, S. Murthy, P. Palit, K. Park, M. Paulini, A. Roberts, Y. Zhou

University of Puerto Rico, Mayaguez, Puerto Rico, USA

S. Malik, R. Sharma

Brown University, Providence, Rhode Island, USA

G. Barone, G. Benelli, D. Cutts, S. Ellis, S. Gottlieb, L. Gouskos, M. Hadley, L. Hay, U. Heintz, K.W. Ho, R. Jain, T. Kwon, L. Lambrecht, G. Landsberg, K.T. Lau, M. LeBlanc, J. Luo, S. Mondal, J. Roloff, T. Russell, S. Sagir89, X. Shen, M. Stamenkovic, S. Sunnarborg, J. Tang, N. Venkatasubramanian

University of Tennessee, Knoxville, Tennessee, USA

A. Abdelhamid, D. Ally, A.G. Delannoy, J. Dervan, S. Fiorendi, J. Harris, T. Holmes, A.R. Kanuganti, N. Karunarathna, J. Lawless, L. Lee, E. Nibigira, B. Skipworth, S. Spanier, C. Thompson, A. Vendrasco

Vanderbilt University, Nashville, Tennessee, USA

U. Acharya, E. Appelt, Y. Chen, S. Greene, A. Gurrola, W. Johns, R. Kunnawalkam Elayavalli, A. Melo, D. Rathjens, F. Romeo, I. Shvetsov, S. Tuo, J. Velkovska, J. Zhang

Texas A&M University, College Station, Texas, USA

D. Aebi, M. Ahmad, T. Akhter, K. Androsov, A. Basnet, A. Bolshov, O. Bouhali90, A. Cagnotta, S. Cooperstein, V. D’Amante, R. Eusebi, P. Flanagan, J. Gilmore, Y. Guo, T. Kamon, R. Mueller, G. Pizzati, A. Safonov

Rice University, Houston, Texas, USA

D. Acosta, A. Agrawal, C. Arbour, T. Carnahan, K.M. Ecklund, F.J. Geurts, I. Krommydas, N. Lewis, W. Li, J. Lin, X. Liu, C. Loizides, O. Miguel Colin, B.P. Padley, R. Redjimi, J. Rotter, C. Vico Villalba, M. Wulansatiti, E. Yigitbasi, Y. Zhang

Texas Tech University, Lubbock, Texas, USA

N. Akchurin, J. Damgov, Y. Feng, N. Gogate, W. Jin, S.W. Lee, C. Madrid, S. Magedov, A. Mankel, T. Peltola, I. Volobouev

Baylor University, Waco, Texas, USA

S. Abdullin, A. Brinkerhoff, E. Collins, M.R. Darwish, J. Dittmann, K. Hatakeyama, V. Hegde, J. Hiltbrand, B. McMaster, J. Samudio, S. Sawant, C. Sutantawibul, J. Wilson

University of Virginia, Charlottesville, Virginia, USA

B. Cardwell, H. Chung, B. Cox, J. Hakala, G. Hamilton Ilha Machado, R. Hirosky, M. Jose, A. Ledovskoy, C. Mantilla, R. Menon Raghunandanan, C. Neu, C. Ramón Álvarez, Z. Wu

University of Wisconsin - Madison, Madison, Wisconsin, USA

A. Aravind, S. Banerjee, K. Black, T. Bose, E. Chavez, R. Cruz, S. Dasu, P. Everaerts, C. Galloni, H. He, M. Herndon, A. Herve, C.K. Koraka, S. Lomte, R. Loveless, J. Marquez, A. Mohammadi, S. Mondal, T. Nelson, G. Parida, D. Pinna, A. Savin, V. Sharma, R. Simeon, W.H. Smith, D. Teague, M. Thakore, A. Thete, A. Warden

Authors affiliated with an international laboratory covered by a cooperation agreement with CERN

S. Afanasiev, V. Alexakhin, Y. Andreev, D. Budkouski, R. Chistov, M. Danilov, T. Dimova, I. Gorbunov, A. Kamenev, V. Karjavine, O. Kodolova91, V. Korenkov, I. Korsakov, A. Kozyrev, A. Lanev, A. Malakhov, V. Matveev, A. Nikitenko92,91, V. Palichik, V. Perelygin, S. Polikarpov, O. Radchenko, M. Savina, V. Shalaev, S. Shmatov, S. Shulha, Y. Skovpen, K. Slizhevskiy, V. Smirnov, O. Teryaev, A. Toropin, N. Voytishin, A. Zarubin, I. Zhizhin Authors affiliated with an institute formerly covered by a cooperation agreement with CERN

L. Dudko, V. Kim{sup}22`, V. Murzin, V. Oreshkin, D. Sosnov ⁴⁶ Also at Kocaeli University, Kocaeli, Türkiye

⁴⁷ Also at Centro Siciliano di Fisica Nucleare e di Struttura della Materia, Catania, Italy

⁴⁸ Also at James Madison University, Harrisonburg, Virginia, USA

⁴⁹ Also at Università degli Studi Guglielmo Marconi, Roma, Italy

⁵⁰ Also at Scuola Superiore Meridionale, Università di Napoli ’Federico II’, Napoli, Italy

⁵¹ Also at Fermi National Accelerator Laboratory, Batavia, Illinois, USA

⁵² Also at Lulea University of Technology, Lulea, Sweden

⁵³ Also at Ain Shams University, Cairo, Egypt

⁵⁴ Also at Consiglio Nazionale delle Ricerche - Istituto Officina dei Materiali, Perugia, Italy

⁵⁵ Also at Boston University, Boston, Massachusetts, USA

⁵⁶ Also at UPES - University of Petroleum and Energy Studies, Dehradun, India

⁵⁷ Also at Institut de Physique des 2 Infinis de Lyon (IP2I), Villeurbanne, France

⁵⁸ Also at Department of Applied Physics, Faculty of Science and Technology, Universiti Kebangsaan Malaysia, Bangi, Malaysia

⁵⁹ Also at Georgian Technical University, Tbilisi, Georgia

⁶⁰ Also at Departamento de Física Instituto Superior Técnico Universidade de Lisboa, Lisbon, Portugal

⁶¹ Also at Trincomalee Campus, Eastern University, Sri Lanka, Nilaveli, Sri Lanka

⁶² Also at Saegis Campus, Nugegoda, Sri Lanka

⁶³ Also at National and Kapodistrian University of Athens, Athens, Greece

⁶⁴ Also at Ecole Polytechnique Fédérale Lausanne, Lausanne, Switzerland

⁶⁵ Also at Universität Zürich, Zurich, Switzerland

⁶⁶ Also at Stefan Meyer Institute for Subatomic Physics (SMI), Vienna, Austria

⁶⁷ Also at Near East University, Research Center of Experimental Health Science, Mersin, Türkiye

⁶⁸ Also at Konya Technical University, Konya, Türkiye

⁶⁹ Also at Izmir Bakircay University Faculty of Engineering and Architecture, Izmir, Türkiye

⁷⁰ Also at Adiyaman University, Adiyaman, Türkiye

⁷¹ Also at Istanbul Sabahattin Zaim University, Istanbul, Türkiye

⁷² Also at Marmara University, Istanbul, Türkiye

⁷³ Also at Milli Savunma University, Naval Academy, Istanbul, Türkiye

⁷⁴ Also at The Science and Technological research Council of Türkiye, Informatics and Information Security Research Center, Gebze/Kocaeli, Türkiye

⁷⁵ Also at Kafkas University, Kars, Türkiye

⁷⁶ Now at Istanbul Okan University, Istanbul, Türkiye

⁷⁷ Also at Istanbul University - Cerrahpasa, Faculty of Engineering, Istanbul, Türkiye

⁷⁸ Also at Istinye University, Istanbul, Türkiye

⁷⁹ Also at Mimar Sinan University, Istanbul, Istanbul, Türkiye

⁸⁰ Also at School of Physics and Astronomy, University of Southampton, Southampton, United Kingdom

⁸¹ Also at Monash University, Faculty of Science, Clayton, Australia

⁸² Also at Università di Torino, Torino, Italy

⁸³ Also at California Lutheran University, Thousand Oaks, California, USA

⁸⁴ Also at California Institute of Technology, Pasadena, California, USA

⁸⁵ Also at United States Naval Academy - Physics Department, Annapolis, Maryland, USA

⁸⁶ Also at Bingol University, Bingol, Türkiye

⁸⁷ Also at Sinop University, Sinop, Türkiye

⁸⁸ Now at another institute formerly covered by a cooperation agreement with CERN

⁸⁹ Also at Karamanoğlu Mehmetbey University, Karaman, Türkiye ^{90}Also at Hamad Bin Khalifa University (HBKU), Doha, Qatar

^{91}Also at Yerevan Physics Institute, Yerevan, Armenia

^{92}Also at Imperial College, London, United Kingdom

References

  1. [1]Particle Data Group, “Review of particle physics”, Int. J. Mod. Phys. A 41 (2026) 2630011, doi:10.1142/S0217751X26300115.DOI
  2. [2]W. Elbers et al., “Constraints on neutrino physics from DESI DR2 BAO and DR1 full shape”, Phys. Rev. D 112 (2025) 083513, doi:10.1103/w9pk-xsk7, arXiv:2503.14744.DOI
  3. [3]P. Minkowski, “μ → eγ at a rate of one out of 10⁹ muon decays?”, Phys. Lett. B 67 (1977) 421, doi:10.1016/0370-2693(77)90435-X.
  4. [4]M. Gell-Mann, P. Ramond, and R. Slansky, “Complex spinors and unified theories”, in Supergravity: Proceedings of the Supergravity Workshop at Stony Brook, p. 315. North-Holland, 1979. arXiv:1306.4669.
  5. [5]T. Yanagida, “Horizontal symmetry and masses of neutrinos”, Prog. Theor. Phys. 64 (1980) 1103, doi:10.1143/PTP.64.1103.DOI
  6. [6]R. N. Mohapatra and G. Senjanović, “Neutrino mass and spontaneous parity nonconservation”, Phys. Rev. Lett. 44 (1980) 912, doi:10.1103/PhysRevLett.44.912.DOI
  7. [7]A. M. Abdullahi et al., “The present and future status of heavy neutral leptons”, J. Phys. G 50 (2023) 020501, doi:10.1088/1361-6471/ac98f9, arXiv:2203.08039.DOI
  8. [8]M. Galeazzi, F. Fontanelli, F. Gatti, and S. Vitale, “Limits on the existence of heavy neutrinos in the range 50–1000 eV from the study of the ¹⁸⁷Re beta decay”, Phys. Rev. Lett. 86 (2001) 1978, doi:10.1103/PhysRevLett.86.1978.DOI
  9. [9]K. H. Hiddemann, H. Daniel, and O. Schwentker, “Limits on neutrino masses from the tritium beta spectrum”, J. Phys. G 21 (1995) 639, doi:10.1088/0954-3899/21/5/008.DOI
  10. [10]CMS Collaboration, “Search for heavy Majorana neutrinos in same-sign dilepton channels in proton-proton collisions at √s = 13 TeV”, JHEP 01 (2019) 122, doi:10.1007/JHEP01(2019)122, arXiv:1806.10905.DOI
  11. [11]CMS Collaboration, “Search for heavy neutral leptons in final states with electrons, muons, and hadronically decaying tau leptons in proton-proton collisions at √s = 13 TeV”, JHEP 06 (2024) 123, doi:10.1007/JHEP06(2024)123, arXiv:2403.00100.DOI
  12. [12]ATLAS Collaboration, “Search for heavy neutral leptons in decays of W bosons produced in 13 TeV pp collisions using prompt and displaced signatures with the ATLAS detector”, JHEP 10 (2019) 265, doi:10.1007/JHEP10(2019)265, arXiv:1905.09787.DOI
  13. [13]CMS Collaboration, “Review of searches for vector-like quarks, vector-like leptons, and heavy neutral leptons in proton-proton collisions at √s = 13 TeV at the CMS experiment”, Phys. Rept. 1115 (2025) 570, doi:10.1016/j.physrep.2024.09.012, arXiv:2405.17605.DOI
  14. [14]LHCb Collaboration, “Search for heavy neutral leptons in B-meson decays”, JHEP 03 (2026) 178, doi:10.1007/JHEP03(2026)178, arXiv:2512.14551.
  15. [15]S. Bar-Shalom et al., “Majorana neutrinos and lepton-number-violating signals in top-quark and W-boson rare decays”, Phys. Lett. B 643 (2006) 342, doi:10.1016/j.physletb.2006.10.060, arXiv:hep-ph/0608309.
  16. [16]Z. Si and K. Wang, “GeV Majorana neutrinos in top-quark decay at the LHC”, Phys. Rev. D 79 (2009) 014034, doi:10.1103/PhysRevD.79.014034, arXiv:0810.5266.DOI
  17. [17]N. Quintero, G. Lopez Castro, and D. Delépine, “Lepton number violation in top quark and neutral B meson decays”, Phys. Rev. D 84 (2011) 096011, doi:10.1103/PhysRevD.84.096011, arXiv:1108.6009. [Erratum: doi:10.1103/PhysRevD.86.079905].DOI
  18. [18]N. Liu et al., “Top quark as a probe of heavy Majorana neutrino at the LHC and future colliders”, Phys. Rev. D 101 (2020) 071701, doi:10.1103/PhysRevD.101.071701, arXiv:1910.00749.DOI
  19. [19]P.-C. Lu, Z.-G. Si, Z. Wang, and X.-H. Yang, “CP violation in top quark decay via heavy Majorana neutrinos at the LHC”, Phys. Rev. D 104 (2021) 115003, doi:10.1103/PhysRevD.104.115003, arXiv:2110.10463.DOI
  20. [20]ATLAS Collaboration, “Search for heavy right-handed Majorana neutrinos in the decay of top quarks produced in proton-proton collisions at √s = 13 TeV with the ATLAS detector”, Phys. Rev. D 110 (2024) 112004, doi:10.1103/PhysRevD.110.112004, arXiv:2408.05000.DOI
  21. [21]HEPData record for this analysis, 2026. doi:10.17182/hepdata.170072.
  22. [22]CMS Collaboration, “The CMS experiment at the CERN LHC”, JINST 3 (2008) S08004, doi:10.1088/1748-0221/3/08/S08004.
  23. [23]CMS Collaboration, “Development of the CMS detector for the CERN LHC Run 3”, JINST 19 (2024) P05064, doi:10.1088/1748-0221/19/05/P05064, arXiv:2309.05466.DOI
  24. [24]CMS Collaboration, “Performance of the CMS level-1 trigger in proton-proton collisions at √s = 13 TeV”, JINST 15 (2020) P10017, doi:10.1088/1748-0221/15/10/P10017, arXiv:2006.10165.DOI
  25. [25]CMS Collaboration, “The CMS trigger system”, JINST 12 (2017) P01020, doi:10.1088/1748-0221/12/01/P01020, arXiv:1609.02366.DOI
  26. [26]CMS Collaboration, “Performance of the CMS high-level trigger during LHC Run 2”, JINST 19 (2024) P11021, doi:10.1088/1748-0221/19/11/P11021, arXiv:2410.17038.DOI
  27. [27]CMS Collaboration, “Electron and photon reconstruction and identification with the CMS experiment at the CERN LHC”, JINST 16 (2021) P05014, doi:10.1088/1748-0221/16/05/P05014, arXiv:2012.06888.DOI
  28. [28]CMS Collaboration, “Performance of the CMS muon detector and muon reconstruction with proton-proton collisions at √s = 13 TeV”, JINST 13 (2018) P06015, doi:10.1088/1748-0221/13/06/P06015, arXiv:1804.04528.DOI
  29. [29]CMS Collaboration, “Description and performance of track and primary-vertex reconstruction with the CMS tracker”, JINST 9 (2014) P10009, doi:10.1088/1748-0221/9/10/P10009, arXiv:1405.6569.DOI
  30. [30]CMS Collaboration, “Particle-flow reconstruction and global event description with the CMS detector”, JINST 12 (2017) P10003, doi:10.1088/1748-0221/12/10/P10003, arXiv:1706.04965.DOI
  31. [31]CMS Collaboration, “Technical proposal for the Phase-II upgrade of the Compact Muon Solenoid”, CMS Technical Proposal CERN-LHCC-2015-010, CMS-TDR-15-02, 2015. https://cds.cern.ch/record/2020886.
  32. [32]K. Rehermann and B. Tweedie, “Efficient identification of boosted semileptonic top quarks at the LHC”, JHEP 03 (2011) 059, doi:10.1007/JHEP03(2011)059, arXiv:1007.2221.DOI
  33. [33]W. Adam, R. Frühwirth, A. Strandlie, and T. Todorov, “Reconstruction of electrons with the Gaussian-sum filter in the CMS tracker at the LHC”, J. Phys. G 31 (2005) N9, doi:10.1088/0954-3899/31/9/N01.DOI
  34. [34]M. Cacciari, G. P. Salam, and G. Soyez, “The anti-kₜ jet clustering algorithm”, JHEP 04 (2008) 063, doi:10.1088/1126-6708/2008/04/063, arXiv:0802.1189.DOI
  35. [35]M. Cacciari, G. P. Salam, and G. Soyez, “FastJet user manual”, Eur. Phys. J. C 72 (2012) 1896, doi:10.1140/epjc/s10052-012-1896-2, arXiv:1111.6097.DOI
  36. [36]CMS Collaboration, “Jet energy scale and resolution in the CMS experiment in pp collisions at 8 TeV”, JINST 12 (2017) P02014, doi:10.1088/1748-0221/12/02/P02014, arXiv:1607.03663.DOI
  37. [37]CMS Collaboration, “Pileup mitigation at CMS in 13 TeV data”, JINST 15 (2020) P09018, doi:10.1088/1748-0221/15/09/P09018, arXiv:2003.00503.DOI
  38. [38]CMS Collaboration, “Performance of the pile up jet identification in CMS for Run 2”, CMS Detector Performance Note CMS-DP-2020-020, 2020. https://cds.cern.ch/record/2715906.
  39. [39]E. Bols et al., “Jet flavour classification using DEEPJET”, JINST 15 (2020) P12012, doi:10.1088/1748-0221/15/12/P12012, arXiv:2008.10519.DOI
  40. [40]CMS Collaboration, “Performance of the DEEPJET b tagging algorithm using 41.9/fb of data from proton-proton collisions at 13 TeV with Phase 1 CMS detector”, CMS Detector Performance Note CMS-DP-2018-058, 2018. https://cds.cern.ch/record/2646773.
  41. [41]CMS Collaboration, “Identification of heavy-flavour jets with the CMS detector in pp collisions at 13 TeV”, JINST 13 (2018) P05011, doi:10.1088/1748-0221/13/05/P05011, arXiv:1712.07158.DOI
  42. [42]CMS Collaboration, “Performance of missing transverse momentum reconstruction in proton-proton collisions at √s = 13 TeV using the CMS detector”, JINST 14 (2019) P07004, doi:10.1088/1748-0221/14/07/P07004, arXiv:1903.06078.DOI
  43. [43]CMS Collaboration, “Precision luminosity measurement in proton-proton collisions at √s = 13 TeV in 2015 and 2016 at CMS”, Eur. Phys. J. C 81 (2021) 800, doi:10.1140/epjc/s10052-021-09538-2, arXiv:2104.01927.DOI
  44. [44]CMS Collaboration, “CMS luminosity measurement for the 2017 data-taking period at √s = 13 TeV”, CMS Physics Analysis Summary CMS-PAS-LUM-17-004, 2018. https://cds.cern.ch/record/2621960/.
  45. [45]CMS Collaboration, “CMS luminosity measurement for the 2018 data-taking period at √s = 13 TeV”, CMS Physics Analysis Summary CMS-PAS-LUM-18-002, 2019. https://cds.cern.ch/record/2676164/.
  46. [46]J. Alwall et al., “The automated computation of tree-level and next-to-leading order differential cross sections, and their matching to parton shower simulations”, JHEP 07 (2014) 079, doi:10.1007/JHEP07(2014)079, arXiv:1405.0301.DOI
  47. [47]D. Alva, T. Han, and R. Ruiz, “Heavy Majorana neutrinos from Wγ fusion at hadron colliders”, JHEP 02 (2015) 072, doi:10.1007/JHEP02(2015)072, arXiv:1411.7305.DOI
  48. [48]C. Degrande, O. Mattelaer, R. Ruiz, and J. Turner, “Fully-automated precision predictions for heavy neutrino production mechanisms at hadron colliders”, Phys. Rev. D 94 (2016) 053002, doi:10.1103/PhysRevD.94.053002, arXiv:1602.06957.DOI
  49. [49]A. Alloul et al., “FEYN-RULES 2.0 — a complete toolbox for tree-level phenomenology”, Comput. Phys. Commun. 185 (2014) 2250, doi:10.1016/j.cpc.2014.04.012, arXiv:1310.1921.DOI
  50. [50]M. Czakon and A. Mitov, “TOP++: a program for the calculation of the top-pair cross-section at hadron colliders”, Comput. Phys. Commun. 185 (2014) 2930, doi:10.1016/j.cpc.2014.06.021, arXiv:1112.5675.DOI
  51. [51]S. Hoeche et al., “Matching parton showers and matrix elements”, in Proc. HERA and the LHC: A Workshop on the Implications of HERA for LHC Physics, p. 288. 2005. arXiv:hep-ph/0602031. doi:10.5170/CERN-2005-014.288.
  52. [52]P. Artoisenet, R. Frederix, O. Mattelaer, and R. Rietkerk, “Automatic spin-entangled decays of heavy resonances in Monte Carlo simulations”, JHEP 03 (2013) 015, doi:10.1007/JHEP03(2013)015, arXiv:1212.3460.DOI
  53. [53]P. Nason, “A new method for combining NLO QCD with shower Monte Carlo algorithms”, JHEP 11 (2004) 040, doi:10.1088/1126-6708/2004/11/040, arXiv:hep-ph/0409146.DOI
  54. [54]S. Frixione, P. Nason, and C. Oleari, “Matching NLO QCD computations with parton shower simulations: the POWHEG method”, JHEP 11 (2007) 070, doi:10.1088/1126-6708/2007/11/070, arXiv:0709.2092.DOI
  55. [55]S. Alioli, P. Nason, C. Oleari, and E. Re, “A general framework for implementing NLO calculations in shower Monte Carlo programs: the POWHEG BOX”, JHEP 06 (2010) 043, doi:10.1007/JHEP06(2010)043, arXiv:1002.2581.DOI
  56. [56]T. Melia, P. Nason, R. Rontsch, and G. Zanderighi, “W⁺W⁻, WZ and ZZ production in the POWHEG BOX”, JHEP 11 (2011) 078, doi:10.1007/JHEP11(2011)078, arXiv:1107.5051.DOI
  57. [57]H. B. Hartanto, B. Jager, L. Reina, and D. Wackeroth, “Higgs boson production in association with top quarks in the POWHEG BOX”, Phys. Rev. D 91 (2015) 094003, doi:10.1103/PhysRevD.91.094003, arXiv:1501.04498.DOI
  58. [58]E. Bagnaschi, G. Degrassi, P. Slavich, and A. Vicini, “Higgs production via gluon fusion in the POWHEG approach in the SM and in the MSSM”, JHEP 02 (2012) 088, doi:10.1007/JHEP02(2012)088, arXiv:1111.2854.DOI
  59. [59]P. Nason and C. Oleari, “NLO Higgs boson production via vector-boson fusion matched with shower in POWHEG”, JHEP 02 (2010) 037, doi:10.1007/JHEP02(2010)037, arXiv:0911.5299.DOI
  60. [60]E. Re, “Single-top Wt-channel production matched with parton showers using the POWHEG method”, Eur. Phys. J. C 71 (2011) 1547, doi:10.1140/epjc/s10052-011-1547-z, arXiv:1009.2450.DOI
  61. [61]S. Frixione, P. Nason, and G. Ridolfi, “A positive-weight next-to-leading-order Monte Carlo for heavy flavour hadroproduction”, JHEP 09 (2007) 126, doi:10.1088/1126-6708/2007/09/126, arXiv:0707.3088.DOI
  62. [62]F. Cascioli et al., “ZZ production at hadron colliders in NNLO QCD”, Phys. Lett. B 735 (2014) 311, doi:10.1016/j.physletb.2014.06.056, arXiv:1405.2219.DOI
  63. [63]Y. Gao et al., “Spin determination of single-produced resonances at hadron colliders”, Phys. Rev. D 81 (2010) 075022, doi:10.1103/PhysRevD.81.075022, arXiv:1001.3396.DOI
  64. [64]T. Sjöstrand et al., “An introduction to PYTHIA 8.2”, Comput. Phys. Commun. 191 (2015) 159, doi:10.1016/j.cpc.2015.01.024, arXiv:1410.3012.
  65. [65]CMS Collaboration, “Extraction and validation of a new set of CMS PYTHIA8 tunes from underlying-event measurements”, Eur. Phys. J. C 80 (2020) 4, doi:10.1140/epjc/s10052-019-7499-4, arXiv:1903.12179.DOI
  66. [66]NNPDF Collaboration, “Parton distributions from high-precision collider data”, Eur. Phys. J. C 77 (2017) 663, doi:10.1140/epjc/s10052-017-5199-5, arXiv:1706.00428.
  67. [67]GEANT4 Collaboration, “GEANT4—a simulation toolkit”, Nucl. Instrum. Meth. A 506 (2003) 250, doi:10.1016/S0168-9002(03)01368-8.
  68. [68]CMS Collaboration, “Search for new physics with same-sign isolated dilepton events with jets and missing transverse energy at the LHC”, JHEP 06 (2011) 077, doi:10.1007/JHEP06(2011)077, arXiv:1104.3168.DOI
  69. [69]CMS Collaboration, “Measurement of the lepton charge asymmetry in inclusive W production in pp collisions at √s = 7 TeV”, JHEP 04 (2011) 050, doi:10.1007/JHEP04(2011)050, arXiv:1103.3470.DOI
  70. [70]H. Bethe and W. Heitler, “On the stopping of fast particles and on the creation of positive electrons”, Proc. Roy. Soc. Lond. A 146 (1934) 83, doi:10.1098/rspa.1934.0140.DOI
  71. [71]CMS Collaboration, “Measurement of tracking efficiency”, CMS Physics Analysis Summary CMS-PAS-TRK-10-002, 2010. https://cds.cern.ch/record/1279139.
  72. [72]CMS Collaboration, “Measurement of the inelastic proton-proton cross section at √s = 13 TeV”, JHEP 07 (2018) 161, doi:10.1007/JHEP07(2018)161, arXiv:1802.02613.DOI
  73. [73]CMS Collaboration, “Performance of the CMS electromagnetic calorimeter in pp collisions at √s = 13 TeV”, JINST 19 (2024) P09004, doi:10.1088/1748-0221/19/09/P09004, arXiv:2403.15518.DOI
  74. [74]M. Grazzini, S. Kallweit, D. Rathlev, and M. Wiesemann, “W±Z production at hadron colliders in NNLO QCD”, Phys. Lett. B 761 (2016) 179, doi:10.1016/j.physletb.2016.08.017, arXiv:1604.08576.DOI
  75. [75]LHC Higgs Cross Section Working Group, “Handbook of LHC Higgs cross sections: 4. Deciphering the nature of the Higgs sector”, CERN Yellow Rep. Monogr. 2 (2017) 1, doi:10.23731/CYRM-2017-002, arXiv:1610.07922.DOI
  76. [76]C. Degrande et al., “Single-top associated production with a Z or H boson at the LHC: the SMEFT interpretation”, JHEP 10 (2018) 005, doi:10.1007/JHEP10(2018)005, arXiv:1804.07773.DOI
  77. [77]CMS Collaboration, “The CMS statistical analysis and combination tool: Combine”, Comput. Softw. Big Sci. 8 (2024) 19, doi:10.1007/s41781-024-00121-4, arXiv:2404.06614.DOI
  78. [78]W. Verkerke and D. P. Kirkby, “The RooFit toolkit for data modeling”, in Proc. 13th International Conference for Computing in High-Energy and Nuclear Physics (CHEP 2003). 2003. arXiv:physics/0306116.
  79. [79]L. Moneta, K. Cranmer, G. Schott, and W. Verkerke, “The RooStats project”, PoS ACAT2010 (2011) 057, doi:10.22323/1.093.0057.DOI
  80. [80]ATLAS and CMS Collaborations, and LHC Higgs Combination Group, “Procedure for the LHC Higgs boson search combination in summer 2011”, CMS Note CMS-NOTE-2011-005, ATL-PHYS-PUB-2011-11, 2011. https://cds.cern.ch/record/1379837.
  81. [81]R. J. Barlow and C. Beeston, “Fitting using finite Monte Carlo samples”, Comput. Phys. Commun. 77 (1993) 219, doi:10.1016/0010-4655(93) 90005-W.DOI
  82. [82]L. Demortier, “P values and nuisance parameters”, in Proc. PHYSTAT LHC Workshop on Statistical Issues for LHC Physics, p. 23. Geneva, Switzerland, 2008. doi:10.5170/CERN-2008-001.
  83. [83]T. Junk, “Confidence level computation for combining searches with small statistics”, Nucl. Instrum. Meth. A 434 (1999) 435, doi:10.1016/S0168-9002(99) 00498-2, arXiv:hep-ex/9902006.DOI
  84. [84]A. L. Read, “Presentation of search results: the CL_s technique”, J. Phys. G 28 (2002) 2693, doi:10.1088/0954-3899/28/10/313.DOI

Paper details

vlnc:W_w2wq2ugz
v1 · latestHistory (1)

published · Fresh — added to the Hub within the last 30 days

added ·

Usage

views · 0

downloads · 0

A view counts once per browser per day, a download each time the PDF is requested.

CC BY 4.0

Original source

  • Physics and Astronomy

Not recorded for this paper · View statement →

Looking for related papers…
Contents