Measurements of the \Pp\Pp\to\Z\Z production cross section and the \Z\to 4\ell branching fraction, and constraints on anomalous triple gauge couplings at \sqrt{s}=13\TeV

Measurements of the production cross section and the branching fraction, and constraints on anomalous triple gauge couplings at

CMS Collaboration
September 25, 2019
Abstract

Four-lepton production in proton-proton collisions, , where or , is studied at a center-of-mass energy of 13\TeVwith the CMS detector at the LHC. The data sample corresponds to an integrated luminosity of 35.9\fbinv. The ZZ production cross section, , measured using events with two opposite-sign, same-flavor lepton pairs produced in the mass region , is consistent with standard model predictions. Differential cross sections are measured and are well described by the theoretical predictions. The Z boson branching fraction to four leptons is measured to be for events with a four-lepton invariant mass in the range and a dilepton mass for all opposite-sign, same-flavor lepton pairs. The results agree with standard model predictions. The invariant mass distribution of the four-lepton system is used to set limits on anomalous ZZZ and ZZ couplings at 95% confidence level: , , , .

\cmsNoteHeader

SMP-16-017

\RCS

\RCS \RCS

\cmsNoteHeader

SMP-16-017

0.1 Introduction

Measurements of diboson production at the CERN LHC allow precision tests of the standard model (SM). In the SM, production proceeds mainly through quark-antiquark - and -channel scattering diagrams. In calculations at higher orders in quantum chromodynamics (QCD), gluon-gluon fusion also contributes via box diagrams with quark loops. There are no tree-level contributions to production from triple gauge boson vertices in the SM. Anomalous triple gauge couplings (aTGC) could be induced by new physics models such as supersymmetry [1]. Nonzero aTGCs may be parametrized using an effective Lagrangian as in Ref. [2]. In this formalism, two and two couplings are allowed by electromagnetic gauge invariance and Lorentz invariance for on-shell bosons. These are described by two CP-violating () and two CP-conserving () parameters, where or .

Previous measurements of the ZZ production cross section by the CMS Collaboration were performed for pairs of on-shell bosons, produced in the dilepton mass range 60–120\GeV [3, 4, 5, 6]. These measurements were made with data sets corresponding to integrated luminosities of 5.1\fbinvat and 19.6\fbinvat in the and decay channels, where or and , or , and with an integrated luminosity of 2.6\fbinvat in the decay channel, where or . All of them agree with SM predictions. The ATLAS Collaboration produced similar results at , 8, and 13\TeV [7, 8, 9, 10], which also agree with the SM. These measurements are important for testing predictions that were recently made available at next-to-next-to-leading order (NNLO) in QCD [11]. Comparing these predictions with data at a range of center-of-mass energies provides information about the electroweak gauge sector of the SM. Because the uncertainty of the CMS measurement at  [6] was dominated by the statistical uncertainty of the observed data, repeating and extending the measurement with a larger sample of proton-proton collision data at improves the precision of the results.

The most stringent previous limits on and aTGCs from CMS were set using the 7 and data samples: , , , and at 95% confidence level (CL) [5, 4]. Similar limits were obtained by the ATLAS Collaboration [12], who also recently produced limits using data [10].

Extending the dilepton mass range to lower values allows measurements of production, where indicates an on-shell boson or an off-shell boson. The resulting sample includes Higgs boson events in the channel, and rare decays of a single boson to four leptons. The decay was studied in detail at LEP [13] and was observed in pp collisions by CMS [14, 6] and ATLAS [15]. Although the branching fraction for this decay is orders of magnitude smaller than that for the decay, the precisely known mass of the boson makes the four-lepton mode useful for calibrating mass measurements of the nearby Higgs boson resonance.

This paper reports a study of four-lepton production (, where and indicate opposite-sign pairs of electrons or muons) at with a data set corresponding to an integrated luminosity of \fbinvrecorded in 2016. Cross sections are measured for nonresonant production of pairs of bosons, , where both bosons are produced on-shell, defined as the mass range 60–120\GeV, and resonant production. Detailed discussion of resonant Higgs boson production decaying to , is beyond the scope of this paper and may be found in Ref. [16].

0.2 The CMS detector

A detailed description of the CMS detector, together with a definition of the coordinate system used and the relevant kinematic variables, can be found in Ref. [17].

The central feature of the CMS apparatus is a superconducting solenoid of 6\unitm internal diameter, providing a magnetic field of 3.8\unitT. 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, which provide coverage in pseudorapidity in a cylindrical barrel and in two endcap regions. Forward calorimeters extend the coverage provided by the barrel and endcap detectors to . Muons are measured in gas-ionization detectors embedded in the steel flux-return yoke outside the solenoid in the range , with detection planes made using three technologies: drift tubes, cathode strip chambers, and resistive plate chambers.

Electron momenta are estimated by combining energy measurements in the ECAL with momentum measurements in the tracker. The momentum resolution for electrons with transverse momentum from decays ranges from 1.7% for nonshowering electrons in the barrel region to 4.5% for showering electrons in the endcaps [18]. Matching muons to tracks identified in the silicon tracker results in a resolution for muons with of 1.3–2.0% in the barrel and better than 6% in the endcaps. The \ptresolution in the barrel is better than 10% for muons with \ptup to 1\TeV [19].

0.3 Signal and background simulation

Signal events are generated with \POWHEG 2.0 [20, 21, 22, 23, 24] at next-to-leading order (NLO) in QCD for quark-antiquark processes and leading order (LO) for quark-gluon processes. This includes , , , and production with a constraint of applied to all pairs of oppositely charged leptons at the generator level to avoid infrared divergences. The process is simulated at LO with \MCFM v7.0 [25]. These samples are scaled to correspond to cross sections calculated at NNLO in QCD for  [11] (a scaling  factor of 1.1) and at NLO in QCD for  [26] ( factor of 1.7). The process is calculated to , where is the strong coupling constant, while the other contributing processes are calculated to ; this higher-order correction is included because the effect is known to be large [26]. Electroweak  production in association with two jets is generated with Phantom v1.2.8 [27].

A sample of Higgs boson events is produced in the gluon-gluon fusion process at NLO with \POWHEG. The Higgs boson decay is modeled with jhugen 3.1.8 [28, 29, 30]. Its cross section is scaled to the NNLO prediction with a  factor of 1.7 [26].

Samples for background processes containing four prompt leptons in the final state, and production, are produced with \MGvATNLOv2.3.3 [31]. The process is generated with \POWHEG.

Samples with aTGC contributions included are generated at LO with sherpa v2.1.1 [32]. Distributions from the sherpa samples are normalized such that the total yield of the SM sample is the same as that of the \POWHEGsample.

The \PYTHIA v8.175 [33, 34, 23] package is used for parton showering, hadronization, and the underlying event simulation, with parameters set by the CUETP8M1 tune [35], for all samples except the samples generated with sherpa, which performs these functions itself. The NNPDF 3.0 [36] set is used as the default set of parton distribution functions (PDFs). For all simulated event samples, the PDFs are calculated to the same order in QCD as the process in the sample.

The detector response is simulated using a detailed description of the CMS detector implemented with the \GEANTfourpackage [37]. The event reconstruction is performed with the same algorithms used for data. The simulated samples include additional interactions per bunch crossing, referred to as pileup. The simulated events are weighted so that the pileup distribution matches the data, with an average of about 27 interactions per bunch crossing.

0.4 Event reconstruction

All long-lived particles—electrons, muons, photons, and charged and neutral hadrons—in each collision event are identified and reconstructed with the CMS particle-flow (PF) algorithm [38] from a combination of the signals from all subdetectors. Reconstructed electrons [18] and muons [19] are considered candidates for inclusion in four-lepton final states if they have and or and .

Lepton candidates are also required to originate from the event vertex, defined as the reconstructed proton-proton interaction vertex with the largest value of summed physics object . The physics objects used in the event vertex definition are the objects returned by a jet finding algorithm [39, 40] applied to all charged tracks associated with the vertex, plus the corresponding associated missing transverse momentum [41]. The distance of closest approach between each lepton track and the event vertex is required to be less than 0.5\unitcm in the plane transverse to the beam axis, and less than 1\unitcm in the direction along the beam axis. Furthermore, the significance of the three-dimensional impact parameter relative to the event vertex, , is required to satisfy for each lepton, where is the distance of closest approach of each lepton track to the event vertex and is its associated uncertainty.

Lepton candidates are required to be isolated from other particles in the event. The relative isolation is defined as

(1)

where the sums run over the charged and neutral hadrons and photons identified by the PF algorithm, in a cone defined by around the lepton trajectory. Here is the azimuthal angle in radians. To minimize the contribution of charged particles from pileup to the isolation calculation, charged hadrons are included only if they originate from the event vertex. The contribution of neutral particles from pileup is . For electrons, is evaluated with the “jet area” method described in Ref. [42]; for muons, it is taken to be half the sum of the of all charged particles in the cone originating from pileup vertices. The factor one-half accounts for the expected ratio of charged to neutral particle energy in hadronic interactions. A lepton is considered isolated if .

The lepton reconstruction, identification, and isolation efficiencies are measured with a “tag-and-probe” technique [43] applied to a sample of data events. The measurements are performed in several bins of and . The electron reconstruction and selection efficiency in the ECAL barrel (endcaps) varies from about 85% (77%) at to about 95% (89%) for , while in the barrel-endcap transition region this efficiency is about 85% averaged over all electrons with . The muons are reconstructed and identified with efficiencies above within .

0.5 Event selection

The primary triggers for this analysis require the presence of a pair of loosely isolated leptons of the same or different flavors [44]. The highest \ptlepton must have , and the subleading lepton must have if it is an electron or if it is a muon. The tracks of the triggering leptons are required to originate within 2 mm of each other in the plane transverse to the beam axis. Triggers requiring a triplet of lower-\ptleptons with no isolation criterion, or a single high-\ptelectron or muon, are also used. An event is used if it passes any trigger regardless of the decay channel. The total trigger efficiency for events within the acceptance of this analysis is greater than 98%.

The four-lepton candidate selections are based on those used in Ref. [45]. A signal event must contain at least two candidates, each formed from an oppositely charged pair of isolated electron candidates or muon candidates. Among the four leptons, the highest \ptlepton must have , and the second-highest \ptlepton must have if it is an electron or if it is a muon. All leptons are required to be separated from each other by , and electrons are required to be separated from muons by .

Within each event, all permutations of leptons giving a valid pair of candidates are considered separately. Within each candidate, the dilepton candidate with an invariant mass closest to 91.2\GeV, taken as the nominal boson mass [46], is denoted and is required to have a mass greater than 40\GeV. The other dilepton candidate is denoted . Both and are required to be less than 120\GeV. All pairs of oppositely charged leptons in the candidate are required to have regardless of their flavor.

If multiple candidates within an event pass all selections, the one with closest to the nominal boson mass is chosen. In the rare case of further ambiguity, which may arise in less than 0.5% of events when five or more passing lepton candidates are found, the candidate that maximizes the scalar sum of the four leptons is chosen.

Additional requirements are applied to select events for measurements of specific processes. The cross section is measured using events where both and are greater than 60\GeV. The branching fraction is measured using events with , a range chosen to retain most of the decays in the resonance while removing most other processes with four-lepton final states. Decays of the bosons to leptons with subsequent decays to electrons and muons are heavily suppressed by requirements on lepton , and the contribution of such events is less than 0.5% of the total yield. If these events pass the selection requirements of the analysis, they are considered signal, while they are not considered at generator level in the cross section unfolding procedure. Thus, the correction for possible decays is included in the efficiency calculation.

0.6 Background estimate

The major background contributions arise from boson and diboson production in association with jets and from \ttbarproduction. In all these cases, particles from jet fragmentation satisfy both lepton identification and isolation criteria, and are thus misidentified as signal leptons.

The probability for such objects to be selected is measured from a sample of events, where denotes a pair of oppositely charged, same-flavor leptons that pass all analysis requirements and satisfy , where is the nominal boson mass. Each event in this sample must have exactly one additional object that passes relaxed identification requirements with no isolation requirements applied. The misidentification probability for each lepton flavor, measured in bins of lepton candidate and , is defined as the ratio of the number of candidates that pass the final isolation and identification requirements to the total number in the sample. The number of events is corrected for the contamination from production and production in which one lepton is not reconstructed. These events have a third genuine, isolated lepton that must be excluded from the misidentification probability calculation. The WZ contamination is suppressed by requiring the missing transverse momentum to be below 25\GeV. The is defined as the magnitude of the missing transverse momentum vector , the projection onto the plane transverse to the beams of the negative vector sum of the momenta of all reconstructed PF candidates in the event, corrected for the jet energy scale. Additionally, the transverse mass calculated with and the of , , is required to be less than 30\GeV. The residual contribution of and events, which may be up to a few percent of the events with passing all selection criteria, is estimated from simulation and subtracted.

To account for all sources of background events, two control samples are used to estimate the number of background events in the signal regions. Both are defined to contain events with a dilepton candidate satisfying all requirements () and two additional lepton candidates . In one control sample, enriched in events, one candidate is required to satisfy the full identification and isolation criteria and the other must fail the full criteria and instead satisfy only the relaxed ones; in the other, enriched in +jets events, both candidates must satisfy the relaxed criteria, but fail the full criteria. The additional leptons must have opposite charge and the same flavor (). From this set of events, the expected number of background events in the signal region, denoted “” in the figures, is obtained by scaling the number of observed events by the misidentification probability for each lepton failing the selection. It is found to be approximately 4% of the total expected yield. The procedure is described in more detail in Ref. [45].

In addition to these nonprompt backgrounds, and processes contribute a smaller number of events with four prompt leptons, which is estimated from simulated samples to be around 1% of the expected yield. In the selection, the contribution from these backgrounds is negligible. The total background contributions to the and signal regions are summarized in Section 0.8.

0.7 Systematic uncertainties

The major sources of systematic uncertainty and their effect on the measured cross sections are summarized in Table 0.7. In both data and simulated event samples, trigger efficiencies are evaluated with a tag-and-probe technique. The ratio of data to simulation is applied to simulated events, and the size of the resulting change in expected yield is taken as the uncertainty in the determination of the trigger efficiency. This uncertainty is around 2% of the final estimated yield. For events, the uncertainty increases to 4%.

\topcaption

The contributions of each source of systematic uncertainty in the cross section measurements. The integrated luminosity uncertainty, and the PDF and scale uncertainties, are considered separately. All other uncertainties are added in quadrature into a single systematic uncertainty. Uncertainties that vary by decay channel are listed as a range. Uncertainty Lepton efficiency 6–10% 2–6% Trigger efficiency 2–4% 2% Statistical (simulation) 1–2% 0.5% Background 0.6–1.3% 0.5–1% Pileup 1–2% 1% PDF 1% 1% , 1% 1% Integrated luminosity 2.5% 2.5%

The lepton identification, isolation, and track reconstruction efficiencies in simulation are corrected with scaling factors derived with a tag-and-probe method and applied as a function of lepton and . To estimate the uncertainties associated with the tag-and-probe technique, the total yield is recomputed with the scaling factors varied up and down by the tag-and-probe fit uncertainties. The uncertainties associated with lepton efficiency in the () signal regions are found to be 6(10)% in the , 3(6)% in the , and 2(7)% in the final states. These uncertainties are higher for events because the leptons generally have lower \pt, and the samples used in the tag-and-probe method have fewer events and more contamination from nonprompt leptons in this low-\ptregion.

Uncertainties due to the effect of factorization () and renormalization () scale choices on the acceptance are evaluated with \POWHEGand \MCFMby varying the scales up and down by a factor of two with respect to the default values . All combinations are considered except those in which and differ by a factor of four. Parametric uncertainties (PDF) are evaluated according to the pdf4lhc prescription [47] in the acceptance calculation, and with NNPDF3.0 in the cross section calculations. An additional theoretical uncertainty arises from scaling the \POWHEG simulated sample from its NLO cross section to the NNLO prediction, and the \MCFM samples from their LO cross sections to the NLO predictions. The change in the acceptance corresponding to this scaling procedure is found to be 1.1%. All these theoretical uncertainties are added in quadrature.

The largest uncertainty in the estimated background yield arises from differences in sample composition between the control sample used to calculate the lepton misidentification probability and the control sample. A further uncertainty arises from the limited number of events in the sample. A systematic uncertainty of 40% is applied to the lepton misidentification probability to cover both effects. The size of this uncertainty varies by channel, but is less than 1% of the total expected yield.

The uncertainty in the integrated luminosity of the data sample is 2.5% [48].

0.8 Cross section measurements

The distributions of the four-lepton mass and the masses of the and candidates are shown in Fig. 1. The expected distributions describe the data well within uncertainties. The SM predictions include nonresonant predictions, production of the SM Higgs boson with mass 125\GeV [49], and resonant production. The backgrounds estimated from data and simulation are also shown. The reconstructed invariant mass of the candidates, and a scatter plot showing the correlation between and in data events, are shown in Fig. 2. In the scatter plot, clusters of events corresponding to , , and production can be seen.

Figure 1: Distributions of (left) the four-lepton invariant mass and (right) the dilepton invariant mass of all bosons in selected four-lepton events. Both selected dilepton candidates are included in each event. In the distribution, bin contents are normalized to a bin width of 25\GeV; horizontal bars on the data points show the range of the corresponding bin. Points represent the data, while filled histograms represent the SM prediction and background estimate. Vertical bars on the data points show their statistical uncertainty. Shaded grey regions around the predicted yield represent combined statistical, systematic, theoretical, and integrated luminosity uncertainties.
Figure 2: (Left): the distribution of the reconstructed mass of , the dilepton candidate closer to the nominal \cPZ boson mass. Points represent the data, while filled histograms represent the SM prediction and background estimate. Vertical bars on the data points show their statistical uncertainty. Shaded grey regions around the predicted yield represent combined statistical, systematic, theoretical, and integrated luminosity uncertainties. (Right): the reconstructed plotted against the reconstructed in data events, with distinctive markers for each final state. For readability, only every fourth event is plotted.

The four-lepton invariant mass distribution below 100\GeVis shown in Fig. 3 (left). Figure 3 (right) shows plotted against for events with between 80 and 100\GeV, and the observed and expected event yields in this mass region are given in Table 0.8. The yield of events in the final state is significantly lower than in the final state because minimum thresholds are higher for electrons than for muons, and inefficiencies in the detection of low- leptons affect electrons more strongly than they affect muons.

Figure 3: (Left): the distribution of the reconstructed four-lepton mass for events selected with . Points represent the data, while filled histograms represent the SM prediction and background estimate. Vertical bars on the data points show their statistical uncertainty. Shaded grey regions around the predicted yield represent combined statistical, systematic, theoretical, and integrated luminosity uncertainties. (Right): the reconstructed plotted against the reconstructed for all data events selected with between 80 and 100\GeV, with distinctive markers for each final state.
\topcaption

The observed and expected yields of four-lepton events in the mass region and estimated yields of background events, shown for each final state and summed in the total expected yield. The first uncertainty is statistical, the second one is systematic. The systematic uncertainties do not include the uncertainty in the integrated luminosity. Final Expected Background Total Observed state expected Total

The reconstructed four-lepton invariant mass is shown in Fig. 4 (left) for events with two on-shell bosons. Figure 4 (right) shows the invariant mass distribution for all  boson candidates in these events. The corresponding observed and expected yields are given in Table 0.8.

Figure 4: Distributions of (left) the four-lepton invariant mass and (right) dilepton candidate mass for four-lepton events selected with both bosons on-shell. Points represent the data, while filled histograms represent the SM prediction and background estimate. Vertical bars on the data points show their statistical uncertainty. Shaded grey regions around the predicted yield represent combined statistical, systematic, theoretical, and integrated luminosity uncertainties. In the distribution, bin contents are normalized to the bin widths, using a unit bin size of 50\GeV; horizontal bars on the data points show the range of the corresponding bin.
\topcaption

The observed and expected yields of events, and estimated yields of background events, shown for each final state and summed in the total expected yield. The first uncertainty is statistical, the second one is systematic. The systematic uncertainties do not include the uncertainty in the integrated luminosity. Decay Expected Background Total Observed channel expected Total

The observed yields are used to evaluate the and production cross sections from a combined fit to the number of observed events in all the final states. The likelihood is a combination of individual channel likelihoods for the signal and background hypotheses with the statistical and systematic uncertainties in the form of scaling nuisance parameters. The fiducial cross section is measured by scaling the cross section in the simulation by the ratio of the measured and predicted event yields given by the fit.

The definitions for the fiducial phase spaces for the and cross section measurements are given in Table 0.8. In the case, the bosons used in the fiducial definition are built by pairing final-state leptons using the same algorithm as is used to build boson candidates from reconstructed leptons. The generator-level leptons used for the fiducial cross section calculation are “dressed” by adding the momenta of generator-level photons within to their momenta.

\topcaption

Fiducial definitions for the reported cross sections. The common requirements are applied for both measurements. Cross section measurement Fiducial requirements Common requirements , , , , (any opposite-sign same-flavor pair)

The measured cross sections are

(2)

The fiducial cross section can be compared to calculated at NLO in QCD with \POWHEGusing the same settings as used for the simulated sample described in Section 0.3, with dynamic scales . The uncertainties correspond to scale and PDF variations, respectively. The fiducial cross section can be compared to calculated with \POWHEGand \MCFMusing the same settings as the simulated samples, or to computed with matrix at NNLO. The \POWHEGand matrix calculations used dynamic scales , while the contribution from \MCFMwas computed with dynamic scales .

The fiducial cross section is scaled to using the acceptance correction factor , estimated with \POWHEG. This factor corrects the fiducial cross section to the phase space with only the 80–100\GeVmass window and requirements, and also includes a correction, , for the contribution of nonresonant four-lepton production to the signal region. The uncertainty takes into account the interference between doubly- and singly-resonant diagrams. The measured cross section is

(3)

The branching fraction for the decay, , is measured by comparing the cross section given by Eq. (3) with the cross section, and is computed as

(4)

where is the cross section times branching fraction calculated at NNLO with fewz v2.0 [50] in the mass range 60–120\GeV. Its uncertainty includes PDF uncertainties and uncertainties in , the charm and bottom quark masses, and the effect of neglected higher-order corrections to the calculation. The factor corrects for the difference in boson mass windows and is estimated using \POWHEG. Its uncertainty includes scale and PDF variations. The nominal to dilepton branching fraction is 0.03366 [46]. The measured value is

(5)

where the theoretical uncertainty includes the uncertainties in , , and . This can be compared with , computed with \MGvATNLO, and is consistent with the CMS and ATLAS measurements at and 13\TeV [14, 15, 6].

The total production cross section for both dileptons produced in the mass range 60–120\GeVand is found to be

(6)

The measured total cross section can be compared to the theoretical value of calculated with a combination of \POWHEGand \MCFMwith the same settings as described for . It can also be compared to \unitpb, calculated at NNLO in QCD via matrix v1.0.0_beta4 [11, 51], or \unitpb, calculated with \MCFMat NLO in QCD with additional contributions from LO diagrams. Both values are calculated with the NNPDF3.0 PDF sets, at NNLO and NLO, respectively, and fixed scales set to .

This measurement agrees with the previously published cross section measured by CMS at 13\TeV [6] based on a 2.6\fbinvdata sample collected in 2015:

(7)

The two measurements can be combined to yield the “2015+2016 cross section”

(8)

The combination was performed once considering the experimental uncertainties to be fully correlated between the 2015 and 2016 data sets, and once considering them to be fully uncorrelated. The results were averaged, and the difference was added linearly to the systematic uncertainty in the combined cross section.

The total cross section is shown in Fig. 5 as a function of the proton-proton center-of-mass energy. Results from CMS [3, 4] and ATLAS [7, 8, 10] are compared to predictions from matrix and \MCFMwith the NNPDF3.0 PDF sets and fixed scales . The matrix prediction uses PDFs calculated at NNLO, while the \MCFMprediction uses NLO PDFs. The uncertainties are statistical (inner bars) and statistical and systematic added in quadrature (outer bars). The band around the matrix predictions reflects scale uncertainties, while the band around the \MCFMpredictions reflects both scale and PDF uncertainties.

Figure 5: The total ZZ cross section as a function of the proton-proton center-of-mass energy. Results from the CMS and ATLAS experiments are compared to predictions from matrix at NNLO in QCD, and \MCFMat NLO in QCD. The \MCFMprediction also includes gluon-gluon initiated production at LO in QCD. Both predictions use NNPDF3.0 PDF sets and fixed scales . Details of the calculations and uncertainties are given in the text. The ATLAS measurements were performed with a boson mass window of 66–116\GeV, and are corrected for the resulting 1.6% difference. Measurements at the same center-of-mass energy are shifted slightly along the horizontal axis for clarity.

The measurement of the differential cross sections provides detailed information about kinematics. The observed yields are unfolded using the iterative technique described in Ref. [52]. Unfolding is performed with the RooUnfold package [53] and regularized by stopping after four iterations. Statistical uncertainties in the data distributions are propagated through the unfolding process to give the statistical uncertainties on the normalized differential cross sections.

The three decay channels, , , and , are combined after unfolding because no differences are expected in their kinematic distributions. The generator-level leptons used for the unfolding are dressed as in the fiducial cross section calculation.

The differential distributions normalized to the fiducial cross sections are presented in Figs. 68 for the combination of the 4, 4, and 22 decay channels. The fiducial cross section definition includes and selections on each lepton, and the 60–120\GeVmass requirement, as described in Table 0.8 and Section 0.4. Figure 6 shows the normalized differential cross sections as functions of the mass and of the system, Fig. 7 shows them as functions of the \pt of all \cPZ bosons and the \pt of the leading lepton in each event, and Fig. 8 shows the angular correlations between the two bosons. The data are corrected for background contributions and compared with the theoretical predictions from \POWHEGand \MCFM, \MGvATNLOand \MCFM, and matrix. The bottom part of each plot shows the ratio of the measured to the predicted values. The bin sizes are chosen according to the resolution of the relevant variables, while also keeping the statistical uncertainties at a similar level in all bins. The data are well reproduced by the simulation except in the low regions, where data tend to have a steeper slope than the prediction.

Figure 6: Differential cross sections normalized to the fiducial cross section for the combined 4, 4, and 22 decay channels as a function of mass (left) and (right) of the system. Points represent the unfolded data; the solid, dashed, and dotted histograms represent the \POWHEG+\MCFM, \MGvATNLO+\MCFM, and matrix predictions for signal, respectively, and the bands around the predictions reflect their combined statistical, scale, and PDF uncertainties. \PYTHIA v8 was used for parton showering, hadronization, and underlying event simulation in the \POWHEG, \MGvATNLO, and \MCFM samples. The lower part of each plot represents the ratio of the measured cross section to the theoretical distributions. The shaded grey areas around the points represent the sum in quadrature of the statistical and systematic uncertainties, while the crosses represent the statistical uncertainties only.
Figure 7: Normalized differential cross sections as a function of the of (left) all bosons and (right) the leading lepton in events. Other details are as described in the caption of Fig. 6.
Figure 8: Normalized differential cross sections as a function of (left) the azimuthal separation of the two bosons and (right) between the -bosons. Other details are as described in the caption of Fig. 6.

Figure 9 shows the normalized differential four-lepton cross section as a function of , subject only to the common requirements of Table 0.8. This includes contributions from the \cPZ and Higgs boson resonances and continuum \cPZ\cPZ production.

Figure 9: The normalized differential four-lepton cross section as a function of the four-lepton mass, subject only to the common requirements of Table 0.8. SM production is included, simulated with \POWHEG. Other details are as described in the caption of Fig. 6.

0.9 Limits on anomalous triple gauge couplings

The presence of aTGCs would increase the yield of events at high four-lepton masses. Figure 10 presents the distribution of the four-lepton reconstructed mass of events with both \cPZ bosons in the mass range 60–120\GeVfor the combined , 4, and channels. This distribution is used to set the limits on possible contributions from aTGCs. Two simulated samples with nonzero aTGCs are shown as examples, along with the SM distribution simulated by both \SHERPAand \POWHEG.

Figure 10: Distribution of the four-lepton reconstructed mass for the combined , 4, and channels. Points represent the data, the filled histograms represent the SM expected yield including signal and irreducible background predictions from simulation and the data-driven background estimate. Unfilled histograms represent examples of aTGC signal predictions (dashed), and the \SHERPASM prediction (solid), included to illustrate the expected shape differences between the \SHERPAand \POWHEGpredictions. Vertical bars on the data points show their statistical uncertainty. The \SHERPAdistributions are normalized such that the SM sample has the same total yield as the \POWHEGsample predicts. Bin contents are normalized to the bin widths, using a unit bin size of 50\GeV; horizontal bars on the data points show the range of the corresponding bin. The last bin includes the “overflow” contribution from events at masses above 1.2\TeV.

The invariant mass distributions are interpolated from the \SHERPAsimulations for different values of the anomalous couplings in the range between 0 and 0.015. For each distribution, only one or two couplings are varied while all others are set to zero. The measured signal is obtained from a comparison of the data to a grid of aTGC models in the and parameter planes. Expected signal values are interpolated between the 2D grid points using a second-degree polynomial, since the cross section for the signal depends quadratically on the coupling parameters. A binned profile likelihood method, Wald Gaussian approximation, and Wilk’s theorem are used to derive one-dimensional limits at a 95% confidence level (CL) on each of the four aTGC parameters, and two-dimensional limits at a 95% CL on the pairs (, ) and (, [46, 54, 55]. When the limits are calculated for each parameter or pair, all other parameters are set to their SM values. The systematic uncertainties described in Section 0.7 are treated as nuisance parameters with log-normal distributions. No form factor is used when deriving the limits so that the results do not depend on any assumed energy scale characterizing new physics. The constraints on anomalous couplings are displayed in Fig. 11. The curves indicate 68 and 95% confidence levels, and the solid dot shows the coordinates where the likelihood reaches its maximum. Coupling values outside the contours are excluded at the corresponding confidence levels. The limits are dominated by statistical uncertainties.

Figure 11: Two-dimensional observed 95% CL limits (solid contour) and expected 68 and 95% CL limits (dashed contour) on the and aTGCs. The left (right) plot shows the exclusion contour in the parameter planes. The values of couplings outside of contours are excluded at the corresponding confidence level. The solid dot is the point at which the likelihood is at its maximum. The solid lines at the center show the observed one-dimensional 95% CL limits for (horizontal) and (vertical). No form factor is used.

The observed one-dimensional 95% CL limits for the and anomalous coupling parameters are:

(9)

These are the most stringent limits to date on anomalous and trilinear gauge boson couplings, improving on the previous strictest results from CMS [5] by factors of two or more and constraining the coupling parameters more than the corresponding ATLAS results [10].

One way to impose unitarity on the aTGC models is to restrict the range of four-lepton invariant mass used in the limit calculation. The limits will then depend on the “cutoff” value used. The computation of the one-dimensional limits is repeated for different maximum allowed values of , and the results are presented in Fig. 12 as a function of this cutoff.

Figure 12: Expected and observed one-dimensional limits on the four aTGC parameters, as a function of an upper cutoff on the invariant mass of the four-lepton system. No form factor is used.

0.10 Summary

A series of measurements of four-lepton final states in proton-proton collisions at have been performed with the CMS detector at the CERN LHC. The measured cross section is for boson masses in the range . The measured branching fraction for boson decays to four leptons is for four-lepton mass in the range and dilepton mass for all oppositely charged same-flavor lepton pairs. Normalized differential cross sections were also measured. All results agree well with the SM predictions. Improved limits on anomalous and triple gauge couplings were established, the most stringent to date.

Acknowledgements.
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 centers and personnel of the Worldwide LHC Computing Grid for delivering so effectively the computing infrastructure essential to our analyses. Finally, we acknowledge the enduring support for the construction and operation of the LHC and the CMS detector provided by the following funding agencies: BMWFW and FWF (Austria); FNRS and FWO (Belgium); CNPq, CAPES, FAPERJ, and FAPESP (Brazil); MES (Bulgaria); CERN; CAS, MoST, and NSFC (China); COLCIENCIAS (Colombia); MSES and CSF (Croatia); RPF (Cyprus); SENESCYT (Ecuador); MoER, ERC IUT, and ERDF (Estonia); Academy of Finland, MEC, and HIP (Finland); CEA and CNRS/IN2P3 (France); BMBF, DFG, and HGF (Germany); GSRT (Greece); OTKA and NIH (Hungary); DAE and DST (India); IPM (Iran); SFI (Ireland); INFN (Italy); MSIP and NRF (Republic of Korea); LAS (Lithuania); MOE and UM (Malaysia); BUAP, CINVESTAV, CONACYT, LNS, SEP, and UASLP-FAI (Mexico); MBIE (New Zealand); PAEC (Pakistan); MSHE and NSC (Poland); FCT (Portugal); JINR (Dubna); MON, RosAtom, RAS, RFBR and RAEP (Russia); MESTD (Serbia); SEIDI, CPAN, PCTI and FEDER (Spain); Swiss Funding Agencies (Switzerland); MST (Taipei); ThEPCenter, IPST, STAR, and NSTDA (Thailand); TUBITAK and TAEK (Turkey); NASU and SFFR (Ukraine); STFC (United Kingdom); DOE and NSF (USA). Individuals have received support from the Marie-Curie program and the European Research Council and Horizon 2020 Grant, contract No. 675440 (European Union); the Leventis Foundation; the A. P. Sloan Foundation; the Alexander von Humboldt Foundation; the Belgian Federal Science Policy Office; the Fonds pour la Formation à la Recherche dans l’Industrie et dans l’Agriculture (FRIA-Belgium); the Agentschap voor Innovatie door Wetenschap en Technologie (IWT-Belgium); the Ministry of Education, Youth and Sports (MEYS) of the Czech Republic; the Council of Science and Industrial Research, India; the HOMING PLUS program of the Foundation for Polish Science, cofinanced from European Union, Regional Development Fund, the Mobility Plus program of the Ministry of Science and Higher Education, the National Science Center (Poland), contracts Harmonia 2014/14/M/ST2/00428, Opus 2014/13/B/ST2/02543, 2014/15/B/ST2/03998, and 2015/19/B/ST2/02861, Sonata-bis 2012/07/E/ST2/01406; the National Priorities Research Program by Qatar National Research Fund; the Programa Severo Ochoa del Principado de Asturias; the Thalis and Aristeia programs cofinanced by EU-ESF and the Greek NSRF; the Rachadapisek Sompot Fund for Postdoctoral Fellowship, Chulalongkorn University and the Chulalongkorn Academic into Its 2nd Century Project Advancement Project (Thailand); the Welch Foundation, contract C-1845; and the Weston Havens Foundation (USA).

References

.11 The CMS Collaboration

Yerevan Physics Institute, Yerevan, Armenia
A.M. Sirunyan, A. Tumasyan \cmsinstskipInstitut für Hochenergiephysik, Wien, Austria
W. Adam, F. Ambrogi, E. Asilar, T. Bergauer, J. Brandstetter, E. Brondolin, M. Dragicevic, J. Erö, M. Flechl, M. Friedl, R. Frühwirth\cmsAuthorMark1, V.M. Ghete, J. Grossmann, J. Hrubec, M. Jeitler\cmsAuthorMark1, A. König, N. Krammer, I. Krätschmer, D. Liko, T. Madlener, I. Mikulec, E. Pree, D. Rabady, N. Rad, H. Rohringer, J. Schieck\cmsAuthorMark1, R. Schöfbeck, M. Spanring, D. Spitzbart, W. Waltenberger, J. Wittmann, C.-E. Wulz\cmsAuthorMark1, M. Zarucki \cmsinstskipInstitute for Nuclear Problems, Minsk, Belarus
V. Chekhovsky, V. Mossolov, J. Suarez Gonzalez \cmsinstskipUniversiteit Antwerpen, Antwerpen, Belgium
E.A. De Wolf, D. Di Croce, X. Janssen, J. Lauwers, M. Van De Klundert, H. Van Haevermaet, P. Van Mechelen, N. Van Remortel \cmsinstskipVrije Universiteit Brussel, Brussel, Belgium
S. Abu Zeid, F. Blekman, J. D’Hondt, I. De Bruyn, J. De Clercq, K. Deroover, G. Flouris, D. Lontkovskyi, S. Lowette, S. Moortgat, L. Moreels, Q. Python, K. Skovpen, S. Tavernier, W. Van Doninck, P. Van Mulders, I. Van Parijs \cmsinstskipUniversité Libre de Bruxelles, Bruxelles, Belgium
H. Brun, B. Clerbaux, G. De Lentdecker, H. Delannoy, G. Fasanella, L. Favart, R. Goldouzian, A. Grebenyuk, G. Karapostoli, T. Lenzi, J. Luetic, T. Maerschalk, A. Marinov, A. Randle-conde, T. Seva, C. Vander Velde, P. Vanlaer, D. Vannerom, R. Yonamine, F. Zenoni, F. Zhang\cmsAuthorMark2 \cmsinstskipGhent University, Ghent, Belgium
A. Cimmino, T. Cornelis, D. Dobur, A. Fagot, M. Gul, I. Khvastunov, D. Poyraz, C. Roskas, S. Salva, M. Tytgat, W. Verbeke, N. Zaganidis \cmsinstskipUniversité Catholique de Louvain, Louvain-la-Neuve, Belgium
H. Bakhshiansohi, O. Bondu, S. Brochet, G. Bruno, C. Caputo, A. Caudron, S. De Visscher, C. Delaere, M. Delcourt, B. Francois, A. Giammanco, A. Jafari, M. Komm, G. Krintiras, V. Lemaitre, A. Magitteri, A. Mertens, M. Musich, K. Piotrzkowski, L. Quertenmont, M. Vidal Marono, S. Wertz \cmsinstskipUniversité de Mons, Mons, Belgium
N. Beliy \cmsinstskipCentro Brasileiro de Pesquisas Fisicas, Rio de Janeiro, Brazil
W.L. Aldá Júnior, F.L. Alves, G.A. Alves, L. Brito, M. Correa Martins Junior, C. Hensel, A. Moraes, M.E. Pol, P. Rebello Teles \cmsinstskipUniversidade do Estado do Rio de Janeiro, Rio de Janeiro, Brazil
E. Belchior Batista Das Chagas, W. Carvalho, J. Chinellato\cmsAuthorMark3, A. Custódio, E.M. Da Costa, G.G. Da Silveira\cmsAuthorMark4, D. De Jesus Damiao, S. Fonseca De Souza, L.M. Huertas Guativa, H. Malbouisson, M. Melo De Almeida, C. Mora Herrera, L. Mundim, H. Nogima, A. Santoro, A. Sznajder, E.J. Tonelli Manganote\cmsAuthorMark3, F. Torres Da Silva De Araujo, A. Vilela Pereira \cmsinstskipUniversidade Estadual Paulista ,  Universidade Federal do ABC ,  São Paulo, Brazil
S. Ahuja, C.A. Bernardes, T.R. Fernandez Perez Tomei, E.M. Gregores, P.G. Mercadante, S.F. Novaes, Sandra S. Padula, D. Romero Abad, J.C. Ruiz Vargas \cmsinstskipInstitute for Nuclear Research and Nuclear Energy of Bulgaria Academy of Sciences
A. Aleksandrov, R. Hadjiiska, P. Iaydjiev, M. Misheva, M. Rodozov, M. Shopova, S. Stoykova, G. Sultanov \cmsinstskipUniversity of Sofia, Sofia, Bulgaria
A. Dimitrov, I. Glushkov, L. Litov, B. Pavlov, P. Petkov \cmsinstskipBeihang University, Beijing, China
W. Fang\cmsAuthorMark5, X. Gao\cmsAuthorMark5 \cmsinstskipInstitute of High Energy Physics, Beijing, China
M. Ahmad, J.G. Bian, G.M. Chen, H.S. Chen, M. Chen, Y. Chen, C.H. Jiang, D. Leggat, H. Liao, Z. Liu, F. Romeo, S.M. Shaheen, A. Spiezia, J. Tao, C. Wang, Z. Wang, E. Yazgan, H. Zhang, J. Zhao \cmsinstskipState Key Laboratory of Nuclear Physics and Technology, Peking University, Beijing, China
Y. Ban, G. Chen, Q. Li, S. Liu, Y. Mao, S.J. Qian, D. Wang, Z. Xu \cmsinstskipUniversidad de Los Andes, Bogota, Colombia
C. Avila, A. Cabrera, L.F. Chaparro Sierra, C. Florez, C.F. González Hernández, J.D. Ruiz Alvarez \cmsinstskipUniversity of Split, Faculty of Electrical Engineering, Mechanical Engineering and Naval Architecture, Split, Croatia
B. Courbon, N. Godinovic, D. Lelas, I. Puljak, P.M. Ribeiro Cipriano, T. Sculac \cmsinstskipUniversity of Split, Faculty of Science, Split, Croatia
Z. Antunovic, M. Kovac \cmsinstskipInstitute Rudjer Boskovic, Zagreb, Croatia
V. Brigljevic, D. Ferencek, K. Kadija, B. Mesic, A. Starodumov\cmsAuthorMark6, T. Susa \cmsinstskipUniversity of Cyprus, Nicosia, Cyprus
M.W. Ather, A. Attikis, G. Mavromanolakis, J. Mousa, C. Nicolaou, F. Ptochos, P.A. Razis, H. Rykaczewski \cmsinstskipCharles University, Prague, Czech Republic
M. Finger\cmsAuthorMark7, M. Finger Jr.\cmsAuthorMark7 \cmsinstskipUniversidad San Francisco de Quito, Quito, Ecuador
E. Carrera Jarrin \cmsinstskipAcademy of Scientific Research and Technology of the Arab Republic of Egypt, Egyptian Network of High Energy Physics, Cairo, Egypt
Y. Assran\cmsAuthorMark8\cmsAuthorMark9, M.A. Mahmoud\cmsAuthorMark10\cmsAuthorMark9, A. Mahrous\cmsAuthorMark11 \cmsinstskipNational Institute of Chemical Physics and Biophysics, Tallinn, Estonia
R.K. Dewanjee, M. Kadastik, L. Perrini, M. Raidal, A. Tiko, C. Veelken \cmsinstskipDepartment of Physics, University of Helsinki, Helsinki, Finland
P. Eerola, J. Pekkanen, M. Voutilainen \cmsinstskipHelsinki Institute of Physics, Helsinki, Finland
J. Härkönen, T. Järvinen, V. Karimäki, R. Kinnunen, T. Lampén, K. Lassila-Perini, S. Lehti, T. Lindén, P. Luukka, E. Tuominen, J. Tuominiemi, E. Tuovinen \cmsinstskipLappeenranta University of Technology, Lappeenranta, Finland
J. Talvitie, T. Tuuva \cmsinstskipIRFU, CEA, Université Paris-Saclay, Gif-sur-Yvette, France
M. Besancon, F. Couderc, M. Dejardin, D. Denegri, J.L. Faure, F. Ferri, S. Ganjour, S. Ghosh, A. Givernaud, P. Gras, G. Hamel de Monchenault, P. Jarry, I. Kucher, E. Locci, M. Machet, J. Malcles, G. Negro, J. Rander, A. Rosowsky, M.Ö. Sahin, M. Titov \cmsinstskipLaboratoire Leprince-Ringuet, Ecole polytechnique, CNRS/IN2P3, Université Paris-Saclay, Palaiseau, France
A. Abdulsalam, I. Antropov, S. Baffioni, F. Beaudette, P. Busson, L. Cadamuro, C. Charlot, R. Granier de Cassagnac, M. Jo, S. Lisniak, A. Lobanov, J. Martin Blanco, M. Nguyen, C. Ochando, G. Ortona, P. Paganini, P. Pigard, S. Regnard, R. Salerno, J.B. Sauvan, Y. Sirois, A.G. Stahl Leiton, T. Strebler, Y. Yilmaz, A. Zabi, A. Zghiche \cmsinstskipUniversité de Strasbourg, CNRS, IPHC UMR 7178, F-67000 Strasbourg, France
J.-L. Agram\cmsAuthorMark12, J. Andrea, D. Bloch, J.-M. Brom, M. Buttignol, E.C. Chabert, N. Chanon, C. Collard, E. Conte\cmsAuthorMark12, X. Coubez, J.-C. Fontaine\cmsAuthorMark12, D. Gelé, U. Goerlach, M. Jansová, A.-C. Le Bihan, N. Tonon, P. Van Hove \cmsinstskipCentre de Calcul de l’Institut National de Physique Nucleaire et de Physique des Particules, CNRS/IN2P3, Villeurbanne, France
S. Gadrat \cmsinstskipUniversité de Lyon, Université Claude Bernard Lyon 1,  CNRS-IN2P3, Institut de Physique Nucléaire de Lyon, Villeurbanne, France
S. Beauceron, C. Bernet, G. Boudoul, R. Chierici, D. Contardo, P. Depasse, H. El Mamouni, J. Fay, L. Finco, S. Gascon, M. Gouzevitch, G. Grenier, B. Ille, F. Lagarde, I.B. Laktineh, M. Lethuillier, L. Mirabito, A.L. Pequegnot, S. Perries, A. Popov\cmsAuthorMark13, V. Sordini, M. Vander Donckt, S. Viret \cmsinstskipGeorgian Technical University, Tbilisi, Georgia
A. Khvedelidze\cmsAuthorMark7 \cmsinstskipTbilisi State University, Tbilisi, Georgia
Z. Tsamalaidze\cmsAuthorMark7 \cmsinstskipRWTH Aachen University, I. Physikalisches Institut, Aachen, Germany
C. Autermann, S. Beranek, L. Feld, M.K. Kiesel, K. Klein, M. Lipinski, M. Preuten, C. Schomakers, J. Schulz, T. Verlage \cmsinstskipRWTH Aachen University, III. Physikalisches Institut A,  Aachen, Germany
A. Albert, E. Dietz-Laursonn, D. Duchardt, M. Endres, M. Erdmann, S. Erdweg, T. Esch, R. Fischer, A. Güth, M. Hamer, T. Hebbeker, C. Heidemann, K. Hoepfner, S. Knutzen, M. Merschmeyer, A. Meyer, P. Millet, S. Mukherjee, M. Olschewski, K. Padeken, T. Pook, M. Radziej, H. Reithler, M. Rieger, F. Scheuch, D. Teyssier, S. Thüer \cmsinstskipRWTH Aachen University, III. Physikalisches Institut B,  Aachen, Germany
G. Flügge, B. Kargoll, T. Kress, A. Künsken, J. Lingemann, T. Müller, A. Nehrkorn, A. Nowack, C. Pistone, O. Pooth, A. Stahl\cmsAuthorMark14 \cmsinstskipDeutsches Elektronen-Synchrotron, Hamburg, Germany
M. Aldaya Martin, T. Arndt, C. Asawatangtrakuldee, K. Beernaert, O. Behnke, U. Behrens, A. Bermúdez Martínez, A.A. Bin Anuar, K. Borras\cmsAuthorMark15, V. Botta, A. Campbell, P. Connor, C. Contreras-Campana, F. Costanza, C. Diez Pardos, G. Eckerlin, D. Eckstein, T. Eichhorn, E. Eren, E. Gallo\cmsAuthorMark16, J. Garay Garcia, A. Geiser, A. Gizhko, J.M. Grados Luyando, A. Grohsjean, P. Gunnellini, M. Guthoff, A. Harb, J. Hauk, M. Hempel\cmsAuthorMark17, H. Jung, A. Kalogeropoulos, M. Kasemann, J. Keaveney, C. Kleinwort, I. Korol, D. Krücker, W. Lange, A. Lelek, T. Lenz, J. Leonard, K. Lipka, W. Lohmann\cmsAuthorMark17, R. Mankel, I.-A. Melzer-Pellmann, A.B. Meyer, G. Mittag, J. Mnich, A. Mussgiller, E. Ntomari, D. Pitzl, A. Raspereza, B. Roland, M. Savitskyi, P. Saxena, R. Shevchenko, S. Spannagel, N. Stefaniuk, G.P. Van Onsem, R. Walsh, Y. Wen, K. Wichmann, C. Wissing, O. Zenaiev \cmsinstskipUniversity of Hamburg, Hamburg, Germany
S. Bein, V. Blobel, M. Centis Vignali, T. Dreyer, E. Garutti, D. Gonzalez, J. Haller, A. Hinzmann, M. Hoffmann, A. Karavdina, R. Klanner, R. Kogler, N. Kovalchuk, S. Kurz, T. Lapsien, I. Marchesini, D. Marconi, M. Meyer, M. Niedziela, D. Nowatschin, F. Pantaleo\cmsAuthorMark14, T. Peiffer, A. Perieanu, C. Scharf, P. Schleper, A. Schmidt, S. Schumann, J. Schwandt, J. Sonneveld, H. Stadie, G. Steinbrück, F.M. Stober, M. Stöver, H. Tholen, D. Troendle, E. Usai, L. Vanelderen, A. Vanhoefer, B. Vormwald \cmsinstskipInstitut für Experimentelle Kernphysik, Karlsruhe, Germany
M. Akbiyik, C. Barth, S. Baur, E. Butz, R. Caspart, T. Chwalek, F. Colombo, W. De Boer, A. Dierlamm, B. Freund, R. Friese, M. Giffels, A. Gilbert, D. Haitz, F. Hartmann\cmsAuthorMark14, S.M. Heindl, U. Husemann, F. Kassel\cmsAuthorMark14, S. Kudella, H. Mildner, M.U. Mozer, Th. Müller, M. Plagge, G. Quast, K. Rabbertz, M. Schröder, I. Shvetsov, G. Sieber, H.J. Simonis, R. Ulrich, S. Wayand, M. Weber, T. Weiler, S. Williamson, C. Wöhrmann, R. Wolf \cmsinstskipInstitute of Nuclear and Particle Physics (INPP),  NCSR Demokritos, Aghia Paraskevi, Greece
G. Anagnostou, G. Daskalakis, T. Geralis, V.A. Giakoumopoulou, A. Kyriakis, D. Loukas, I. Topsis-Giotis \cmsinstskipNational and Kapodistrian University of Athens, Athens, Greece
G. Karathanasis, S. Kesisoglou, A. Panagiotou, N. Saoulidou \cmsinstskipNational Technical University of Athens, Athens, Greece
K. Kousouris \cmsinstskipUniversity of Ioánnina, Ioánnina, Greece
I. Evangelou, C. Foudas, P. Kokkas, S. Mallios, N. Manthos, I. Papadopoulos, E. Paradas, J. Strologas, F.A. Triantis \cmsinstskipMTA-ELTE Lendület CMS Particle and Nuclear Physics Group, Eötvös Loránd University, Budapest, Hungary
M. Csanad, N. Filipovic, G. Pasztor, G.I. Veres\cmsAuthorMark18 \cmsinstskipWigner Research Centre for Physics, Budapest, Hungary
G. Bencze, C. Hajdu, D. Horvath\cmsAuthorMark19, Á. Hunyadi, F. Sikler, V. Veszpremi, G. Vesztergombi\cmsAuthorMark18, A.J. Zsigmond \cmsinstskipInstitute of Nuclear Research ATOMKI, Debrecen, Hungary
N. Beni, S. Czellar, J. Karancsi\cmsAuthorMark20, A. Makovec, J. Molnar, Z. Szillasi \cmsinstskipInstitute of Physics, University of Debrecen, Debrecen, Hungary
M. Bartók\cmsAuthorMark18, P. Raics, Z.L. Trocsanyi, B. Ujvari \cmsinstskipIndian Institute of Science (IISc),  Bangalore, India
S. Choudhury, J.R. Komaragiri \cmsinstskipNational Institute of Science Education and Research, Bhubaneswar, India
S. Bahinipati\cmsAuthorMark21, S. Bhowmik, P. Mal, K. Mandal, A. Nayak\cmsAuthorMark22, D.K. Sahoo\cmsAuthorMark21, N. Sahoo, S.K. Swain \cmsinstskipPanjab University, Chandigarh, India
S. Bansal, S.B. Beri, V. Bhatnagar, R. Chawla, N. Dhingra, A.K. Kalsi, A. Kaur, M. Kaur, R. Kumar, P. Kumari, A. Mehta, J.B. Singh, G. Walia \cmsinstskipUniversity of Delhi, Delhi, India
Ashok Kumar, Aashaq Shah, A. Bhardwaj, S. Chauhan, B.C. Choudhary, R.B. Garg, S. Keshri, A. Kumar, S. Malhotra, M. Naimuddin, K. Ranjan, R. Sharma \cmsinstskipSaha Institute of Nuclear Physics, HBNI, Kolkata, India
R. Bhardwaj, R. Bhattacharya, S. Bhattacharya, U. Bhawandeep, S. Dey, S. Dutt, S. Dutta, S. Ghosh, N. Majumdar, A. Modak, K. Mondal, S. Mukhopadhyay, S. Nandan, A. Purohit, A. Roy, D. Roy, S. Roy Chowdhury, S. Sarkar, M. Sharan, S. Thakur \cmsinstskipIndian Institute of Technology Madras, Madras, India
P.K. Behera \cmsinstskipBhabha Atomic Research Centre, Mumbai, India
R. Chudasama, D. Dutta, V. Jha, V. Kumar, A.K. Mohanty\cmsAuthorMark14, P.K. Netrakanti, L.M. Pant, P. Shukla, A. Topkar \cmsinstskipTata Institute of Fundamental Research-A, Mumbai, India
T. Aziz, S. Dugad, B. Mahakud, S. Mitra, G.B. Mohanty, N. Sur, B. Sutar \cmsinstskipTata Institute of Fundamental Research-B, Mumbai, India
S. Banerjee, S. Bhattacharya, S. Chatterjee, P. Das, M. Guchait, Sa. Jain, S. Kumar, M. Maity\cmsAuthorMark23, G. Majumder, K. Mazumdar, T. Sarkar\cmsAuthorMark23, N. Wickramage\cmsAuthorMark24 \cmsinstskipIndian Institute of Science Education and Research (IISER),  Pune, India
S. Chauhan, S. Dube, V. Hegde, A. Kapoor, K. Kothekar, S. Pandey, A. Rane, S. Sharma \cmsinstskipInstitute for Research in Fundamental Sciences (IPM),  Tehran, Iran
S. Chenarani\cmsAuthorMark25, E. Eskandari Tadavani, S.M. Etesami\cmsAuthorMark25, M. Khakzad, M. Mohammadi Najafabadi, M. Naseri, S. Paktinat Mehdiabadi\cmsAuthorMark26, F. Rezaei Hosseinabadi, B. Safarzadeh\cmsAuthorMark27, M. Zeinali \cmsinstskipUniversity College Dublin, Dublin, Ireland
M. Felcini, M. Grunewald \cmsinstskipINFN Sezione di Bari , Università di Bari , Politecnico di Bari ,  Bari, Italy
M. Abbrescia, C. Calabria, A. Colaleo, D. Creanza, L. Cristella, N. De Filippis, M. De Palma, F. Errico, L. Fiore, G. Iaselli, S. Lezki, G. Maggi, M. Maggi, G. Miniello, S. My, S. Nuzzo, A. Pompili, G. Pugliese, R. Radogna, A. Ranieri, G. Selvaggi, A. Sharma, L. Silvestris\cmsAuthorMark14, R. Venditti, P. Verwilligen \cmsinstskipINFN Sezione di Bologna , Università di Bologna ,  Bologna, Italy
G. Abbiendi, C. Battilana, D. Bonacorsi, S. Braibant-Giacomelli, R. Campanini, P. Capiluppi, A. Castro, F.R. Cavallo, S.S. Chhibra, G. Codispoti, M. Cuffiani, G.M. Dallavalle, F. Fabbri, A. Fanfani, D. Fasanella, P. Giacomelli, C. Grandi, L. Guiducci, S. Marcellini, G. Masetti, A. Montanari, F.L. Navarria, A. Perrotta, A.M. Rossi, T. Rovelli, G.P. Siroli, N. Tosi \cmsinstskipINFN Sezione di Catania , Università di Catania ,  Catania, Italy
S. Albergo, S. Costa, A. Di Mattia, F. Giordano, R. Potenza, A. Tricomi, C. Tuve \cmsinstskipINFN Sezione di Firenze , Università di Firenze ,  Firenze, Italy
G. Barbagli, K. Chatterjee, V. Ciulli, C. Civinini, R. D’Alessandro, E. Focardi, P. Lenzi, M. Meschini, S. Paoletti, L. Russo\cmsAuthorMark28, G. Sguazzoni, D. Strom, L. Viliani\cmsAuthorMark14 \cmsinstskipINFN Laboratori Nazionali di Frascati, Frascati, Italy
L. Benussi, S. Bianco, F. Fabbri, D. Piccolo, F. Primavera\cmsAuthorMark14 \cmsinstskipINFN Sezione di Genova , Università di Genova ,  Genova, Italy
V. Calvelli, F. Ferro, E. Robutti, S. Tosi \cmsinstskipINFN Sezione di Milano-Bicocca , Università di Milano-Bicocca ,  Milano, Italy
A. Benaglia, L. Brianza, F. Brivio, V. Ciriolo, M.E. Dinardo, S. Fiorendi, S. Gennai, A. Ghezzi, P. Govoni, M. Malberti, S. Malvezzi, R.A. Manzoni, D. Menasce, L. Moroni, M. Paganoni, K. Pauwels, D. Pedrini, S. Pigazzini\cmsAuthorMark29, S. Ragazzi, T. Tabarelli de Fatis \cmsinstskipINFN Sezione di Napoli , Università di Napoli ’Federico II’ , Napoli, Italy, Università della Basilicata , Potenza, Italy, Università G. Marconi , Roma, Italy
S. Buontempo, N. Cavallo, S. Di Guida\cmsAuthorMark14, F. Fabozzi, F. Fienga, A.O.M. Iorio, W.A. Khan, L. Lista, S. Meola\cmsAuthorMark14, P. Paolucci\cmsAuthorMark14, C. Sciacca, F. Thyssen \cmsinstskipINFN Sezione di Padova , Università di Padova , Padova, Italy, Università di Trento , Trento, Italy
P. Azzi\cmsAuthorMark14, N. Bacchetta, L. Benato, D. Bisello, A. Boletti, R. Carlin, A. Carvalho Antunes De Oliveira, P. Checchia, M. Dall’Osso, P. De Castro Manzano, T. Dorigo, U. Dosselli, U. Gasparini, A. Gozzelino, S. Lacaprara, P. Lujan, M. Margoni, A.T. Meneguzzo, N. Pozzobon, P. Ronchese, R. Rossin, F. Simonetto, E. Torassa, S. Ventura, M. Zanetti, P. Zotto \cmsinstskipINFN Sezione di Pavia , Università di Pavia ,  Pavia, Italy
A. Braghieri, A. Magnani, P. Montagna, S.P. Ratti, V. Re, M. Ressegotti, C. Riccardi, P. Salvini, I. Vai, P. Vitulo \cmsinstskipINFN Sezione di Perugia , Università di Perugia ,  Perugia, Italy
L. Alunni Solestizi, M. Biasini