Loading [MathJax]/jax/output/CommonHTML/jax.js
CMS logoCMS event Hgg
Compact Muon Solenoid
LHC, CERN

CMS-SMP-22-017 ; CERN-EP-2025-013
Measurements of the inclusive W and Z boson production cross sections and their ratios in proton-proton collisions at s= 13.6 TeV
Submitted to J. High Energy Phys.
Abstract: Measurements are presented of the W and Z boson production cross sections in proton-proton collisions at a center-of-mass energy of 13.6 TeV. Data collected in 2022 and corresponding to an integrated luminosity of 5.01 fb1 with one or two identified muons in the final state are analyzed. The results for the products of total inclusive cross sections and branching fractions for muonic decays of W and Z bosons are 11.93 ± 0.08 (syst) ± 0.17 (lumi) +0.070.07 (acceptance) nb for W+ boson production, 8.86 ± 0.06 (syst) ± 0.12 (lumi) +0.050.06 (acceptance) nb for W boson production, and 2.021 ± 0.009 (syst) ± 0.028 (lumi) +0.0110.013 (acceptance) nb for the Z boson production in the dimuon mass range of 60--120 GeV, all with negligible statistical uncertainties. Furthermore, the corresponding fiducial cross sections, as well as cross section ratios for both fiducial and total phase space, are provided. The ratios include charge-separated results for W boson production (W+ and W) and the sum of the two contributions (W±), each relative to the measured Z boson production cross section. Additionally, the ratio of the measured cross sections for W+ and W boson production is reported. All measurements are in agreement with theoretical predictions, calculated at next-to-next-to-leading order accuracy in quantum chromodynamics.
Figures & Tables Summary References CMS Publications
Figures

png pdf
Figure 1:
Corrected distributions of the pT of the leading (upper left) and trailing (upper right) muon, as well as pmissT (lower) in the Z boson signal region. The distributions of both data and combined simulation in the upper panels are normalized because the scale of the signal sample is extracted later in the fit. The lower panels show the ratio of the normalized data to the normalized prediction. Overflow entries are included in the rightmost bin, respectively. Differences in the shape between data and prediction are fully covered by the systematic uncertainty band. The systematic uncertainties are described in detail in Section 9. The vertical uncertainty bars on the data represent the statistical uncertainty.

png pdf
Figure 1-a:
Corrected distributions of the pT of the leading (upper left) and trailing (upper right) muon, as well as pmissT (lower) in the Z boson signal region. The distributions of both data and combined simulation in the upper panels are normalized because the scale of the signal sample is extracted later in the fit. The lower panels show the ratio of the normalized data to the normalized prediction. Overflow entries are included in the rightmost bin, respectively. Differences in the shape between data and prediction are fully covered by the systematic uncertainty band. The systematic uncertainties are described in detail in Section 9. The vertical uncertainty bars on the data represent the statistical uncertainty.

png pdf
Figure 1-b:
Corrected distributions of the pT of the leading (upper left) and trailing (upper right) muon, as well as pmissT (lower) in the Z boson signal region. The distributions of both data and combined simulation in the upper panels are normalized because the scale of the signal sample is extracted later in the fit. The lower panels show the ratio of the normalized data to the normalized prediction. Overflow entries are included in the rightmost bin, respectively. Differences in the shape between data and prediction are fully covered by the systematic uncertainty band. The systematic uncertainties are described in detail in Section 9. The vertical uncertainty bars on the data represent the statistical uncertainty.

png pdf
Figure 1-c:
Corrected distributions of the pT of the leading (upper left) and trailing (upper right) muon, as well as pmissT (lower) in the Z boson signal region. The distributions of both data and combined simulation in the upper panels are normalized because the scale of the signal sample is extracted later in the fit. The lower panels show the ratio of the normalized data to the normalized prediction. Overflow entries are included in the rightmost bin, respectively. Differences in the shape between data and prediction are fully covered by the systematic uncertainty band. The systematic uncertainties are described in detail in Section 9. The vertical uncertainty bars on the data represent the statistical uncertainty.

png pdf
Figure 2:
Example figures to illustrate the extrapolation procedure to estimate the QCD multijet background from the QCD-enriched control region. The yield in each bin is obtained by subtracting the simulated non-QCD contributions from the measured data. At high values of the transverse mass, the statistical uncertainty of the QCD yield increases, represented by the vertical uncertainty bars, leaving more freedom for the polynomial fit.

png pdf
Figure 2-a:
Example figures to illustrate the extrapolation procedure to estimate the QCD multijet background from the QCD-enriched control region. The yield in each bin is obtained by subtracting the simulated non-QCD contributions from the measured data. At high values of the transverse mass, the statistical uncertainty of the QCD yield increases, represented by the vertical uncertainty bars, leaving more freedom for the polynomial fit.

png pdf
Figure 2-b:
Example figures to illustrate the extrapolation procedure to estimate the QCD multijet background from the QCD-enriched control region. The yield in each bin is obtained by subtracting the simulated non-QCD contributions from the measured data. At high values of the transverse mass, the statistical uncertainty of the QCD yield increases, represented by the vertical uncertainty bars, leaving more freedom for the polynomial fit.

png pdf
Figure 2-c:
Example figures to illustrate the extrapolation procedure to estimate the QCD multijet background from the QCD-enriched control region. The yield in each bin is obtained by subtracting the simulated non-QCD contributions from the measured data. At high values of the transverse mass, the statistical uncertainty of the QCD yield increases, represented by the vertical uncertainty bars, leaving more freedom for the polynomial fit.

png pdf
Figure 2-d:
Example figures to illustrate the extrapolation procedure to estimate the QCD multijet background from the QCD-enriched control region. The yield in each bin is obtained by subtracting the simulated non-QCD contributions from the measured data. At high values of the transverse mass, the statistical uncertainty of the QCD yield increases, represented by the vertical uncertainty bars, leaving more freedom for the polynomial fit.

png pdf
Figure 3:
Result of the extrapolated mT distribution for the QCD multijet background in the signal region from the QCD-enriched control region. The template normalization is left freely floating in the fit, as the final normalization is determined in the fit. The vertical uncertainty bars correspond to the extrapolation uncertainty from the fit.

png pdf
Figure 3-a:
Result of the extrapolated mT distribution for the QCD multijet background in the signal region from the QCD-enriched control region. The template normalization is left freely floating in the fit, as the final normalization is determined in the fit. The vertical uncertainty bars correspond to the extrapolation uncertainty from the fit.

png pdf
Figure 3-b:
Result of the extrapolated mT distribution for the QCD multijet background in the signal region from the QCD-enriched control region. The template normalization is left freely floating in the fit, as the final normalization is determined in the fit. The vertical uncertainty bars correspond to the extrapolation uncertainty from the fit.

png pdf
Figure 4:
The post-fit distributions of (upper left) μ+νμ, (upper right) μ¯νμ, and (lower) μ+μ signal regions. The lower panel in each plot shows the ratio of the number of events observed in data to that of the signal and background predictions. Overflow entries are included in the rightmost bin for the two upper plots. The vertical uncertainty bars on the data represent the statistical uncertainty.

png pdf
Figure 4-a:
The post-fit distributions of (upper left) μ+νμ, (upper right) μ¯νμ, and (lower) μ+μ signal regions. The lower panel in each plot shows the ratio of the number of events observed in data to that of the signal and background predictions. Overflow entries are included in the rightmost bin for the two upper plots. The vertical uncertainty bars on the data represent the statistical uncertainty.

png pdf
Figure 4-b:
The post-fit distributions of (upper left) μ+νμ, (upper right) μ¯νμ, and (lower) μ+μ signal regions. The lower panel in each plot shows the ratio of the number of events observed in data to that of the signal and background predictions. Overflow entries are included in the rightmost bin for the two upper plots. The vertical uncertainty bars on the data represent the statistical uncertainty.

png pdf
Figure 4-c:
The post-fit distributions of (upper left) μ+νμ, (upper right) μ¯νμ, and (lower) μ+μ signal regions. The lower panel in each plot shows the ratio of the number of events observed in data to that of the signal and background predictions. Overflow entries are included in the rightmost bin for the two upper plots. The vertical uncertainty bars on the data represent the statistical uncertainty.

png pdf
Figure 5:
Comparison of measured product of fiducial (left) and total (right) cross sections and branching fractions with different theoretical predictions at NNLO+NNLL QCD accuracy. The vertical uncertainty bars on the markers represent the total uncertainty of the theoretical prediction.

png pdf
Figure 5-a:
Comparison of measured product of fiducial (left) and total (right) cross sections and branching fractions with different theoretical predictions at NNLO+NNLL QCD accuracy. The vertical uncertainty bars on the markers represent the total uncertainty of the theoretical prediction.

png pdf
Figure 5-b:
Comparison of measured product of fiducial (left) and total (right) cross sections and branching fractions with different theoretical predictions at NNLO+NNLL QCD accuracy. The vertical uncertainty bars on the markers represent the total uncertainty of the theoretical prediction.

png pdf
Figure 6:
Comparison of measured fiducial (left) and total (right) ratios with different theoretical predictions at NNLO+NNLL QCD accuracy. The vertical uncertainty bars on the markers represent the total uncertainty of the theoretical prediction.

png pdf
Figure 6-a:
Comparison of measured fiducial (left) and total (right) ratios with different theoretical predictions at NNLO+NNLL QCD accuracy. The vertical uncertainty bars on the markers represent the total uncertainty of the theoretical prediction.

png pdf
Figure 6-b:
Comparison of measured fiducial (left) and total (right) ratios with different theoretical predictions at NNLO+NNLL QCD accuracy. The vertical uncertainty bars on the markers represent the total uncertainty of the theoretical prediction.

png pdf
Figure 7:
Comparison of measured products of total cross section and branching fractions for W and Z boson production at different center-of-mass energies with the corresponding theoretical prediction at NNLO+NNLL QCD accuracy obtained from DYTURBO. The uncertainties in the theoretical prediction include variations of the renormalization and factorization scales, as well as the PDF uncertainty evaluated with the NNPDF 3.1 set. The vertical uncertainty bars on the markers represent the total uncertainty of the measurement.
Tables

png pdf
Table 1:
Predictions for the product of the fiducial and total inclusive cross sections and branching fractions. The first uncertainty is the PDF uncertainty, the second is the QCD scale uncertainty of the calculation, and the third is the integration uncertainty.

png pdf
Table 2:
Predictions for the fiducial and total inclusive cross section ratios. The first uncertainty is the PDF uncertainty, the second is the QCD scale uncertainty of the calculation, and the third is the integration uncertainty.

png pdf
Table 3:
Post-fit uncertainties in percent for the fiducial cross section measurement. For completeness, also the integrated luminosity and statistical uncertainty are given.

png pdf
Table 4:
Post-fit uncertainties in percent for the fiducial cross section ratio measurement. For completeness, also the statistical uncertainty is given.

png pdf
Table 5:
Pre-fit event yields in the fiducial region. A dash indicates that the corresponding contribution is found to be negligible in this signal region.

png pdf
Table 6:
Post-fit event yields in the fiducial region. The post-fit uncertainties include only statistical and systematic uncertainties, but not the uncertainty in the luminosity. The individual uncertainties in the event yields for a given process are derived by taking the full covariance matrix into account. A dash indicates that the corresponding contribution is found to be negligible in this signal region.

png pdf
Table 7:
Results for the product of the fiducial and total inclusive cross sections and branching fractions measurements. For the measured values the quoted uncertainty represents the systematic uncertainty, while the statistical uncertainty is negligible. For the acceptance predictions, as explained in Section 5, the first uncertainty is the PDF uncertainty, the second is the scale uncertainty, and the third is the integration uncertainty of the calculation.

png pdf
Table 8:
Ratios of the measured product of the fiducial and total inclusive cross sections and branching fractions along with the corresponding acceptance predictions. Since some contributions of the systematic uncertainty, most prominently the luminosity uncertainty, cancel out in the ratios, the statistical component becomes relevant. For the acceptance predictions, as explained in Section 5, the first uncertainty is the PDF uncertainty, the second is the scale uncertainty, and the third is the integration uncertainty of the calculation.
Summary
Measurements of fiducial and total inclusive W and Z boson production cross sections multiplied by the Wμνμ and Zμ+μ branching fractions, respectively, in proton-proton collisions at 13.6 TeV are presented. Muon final states are studied in data samples collected with the CMS detector corresponding to an integrated luminosity of 5.01 ± 0.07 fb1. The measured Z boson cross section is defined as the inclusive Drell--Yan Z/γ production of muon pairs in the invariant mass range of 60--120 GeV, where the production of a Z boson is the dominant contribution.
References
1 ATLAS Collaboration Measurement of W±-boson and Z-boson production cross-sections in pp collisions at s= 2.76 TeV with the ATLAS detector EPJC 79 (2019) 901 1907.03567
2 CMS Collaboration Study of Z production in PbPb and pp collisions at {\sqrt{\smash[b]{s_{_{\mathrm{NN}}}}}= 2.76 TeV in the dimuon and dielectron decay channels JHEP 03 (2015) 022 CMS-HIN-13-004
1410.4825
3 ATLAS Collaboration Measurements of W and Z boson production in pp collisions at s= 5.02 TeV with the ATLAS detector EPJC 79 (2019) 128 1810.08424
4 CMS Collaboration Measurement of the inclusive cross sections for W and Z boson production in proton-proton collisions at s= 5.02 and 13 TeV Submitted to JHEP, 2024 CMS-SMP-20-004
2408.03744
5 ATLAS Collaboration Precision measurement and interpretation of inclusive W+, W and Z/γ production cross sections with the ATLAS detector EPJC 77 (2017) 367 1612.03016
6 CMS Collaboration Measurements of inclusive W and Z cross sections in pp collisions at s= 7 TeV JHEP 01 (2011) 080 CMS-EWK-10-002
1012.2466
7 CMS Collaboration Measurement of the inclusive W and Z production cross sections in pp collisions at s= 7 TeV with the CMS experiment JHEP 10 (2011) 132 CMS-EWK-10-005
1107.4789
8 LHCb Collaboration Measurement of the forward W boson cross-section in pp collisions at s= 7 TeV JHEP 12 (2014) 079 1408.4354
9 LHCb Collaboration Measurement of the cross-section for Ze+e production in pp collisions at s= 7 TeV JHEP 02 (2013) 106 1212.4620
10 LHCb Collaboration Measurement of the forward Z boson production cross-section in pp collisions at s= 7 TeV JHEP 08 (2015) 039 1505.07024
11 LHCb Collaboration Measurement of forward Weν production in pp collisions at s= 8 TeV JHEP 10 (2016) 030 1608.01484
12 LHCb Collaboration Measurement of the forward Z boson production cross-section in pp collisions at s= 13 TeV JHEP 09 (2016) 136 1607.06495
13 LHCb Collaboration Measurement of forward W and Z boson production in pp collisions at s= 8 TeV JHEP 01 (2016) 155 1511.08039
14 CMS Collaboration Measurement of inclusive W and Z boson production cross sections in pp collisions at s= 8 TeV PRL 112 (2014) 191802 CMS-SMP-12-011
1402.0923
15 ATLAS Collaboration Measurement of W± and Z-boson production cross sections in pp collisions at s= 13 TeV with the ATLAS detector PLB 759 (2016) 601 1603.09222
16 ATLAS Collaboration Measurement of vector boson production cross sections and their ratios using pp collisions at s= 13.6 TeV with the ATLAS detector PLB 854 (2024) 138725 2403.12902
17 CMS Collaboration HEPData record for this analysis link
18 CMS Collaboration The CMS experiment at the CERN LHC JINST 3 (2008) S08004
19 CMS Collaboration Development of the CMS detector for the CERN LHC \mboxRun 3 JINST 19 (2024) P05064 CMS-PRF-21-001
2309.05466
20 CMS Collaboration Performance of the CMS Level-1 trigger in proton-proton collisions at s= 13 TeV JINST 15 (2020) P10017 CMS-TRG-17-001
2006.10165
21 CMS Collaboration The CMS trigger system JINST 12 (2017) P01020 CMS-TRG-12-001
1609.02366
22 CMS Collaboration Performance of the CMS high-level trigger during LHC \mboxRun 2 JINST 19 (2024) P11021 CMS-TRG-19-001
2410.17038
23 CMS Collaboration Particle-flow reconstruction and global event description with the CMS detector JINST 12 (2017) P10003 CMS-PRF-14-001
1706.04965
24 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
link
25 CMS Collaboration Performance of missing transverse momentum reconstruction in proton-proton collisions at s= 13 TeV using the CMS detector JINST 14 (2019) P07004 CMS-JME-17-001
1903.06078
26 J. Alwall et al. The automated computation of tree-level and next-to-leading order differential cross sections, and their matching to parton shower simulations JHEP 07 (2014) 079 1405.0301
27 R. Frederix and S. Frixione Merging meets matching in MC@NLO JHEP 12 (2012) 061 1209.6215
28 P. Nason A new method for combining NLO QCD with shower Monte Carlo algorithms JHEP 11 (2004) 040 hep-ph/0409146
29 S. Frixione, P. Nason, and C. Oleari Matching NLO QCD computations with parton shower simulations: the POWHEG method JHEP 11 (2007) 070 0709.2092
30 S. Alioli, P. Nason, C. Oleari, and E. Re A general framework for implementing NLO calculations in shower Monte Carlo programs: the POWHEG \textscbox JHEP 06 (2010) 043 1002.2581
31 R. Frederix, E. Re, and P. Torrielli Single-top t-channel hadroproduction in the four-flavors scheme with POWHEG and aMC@NLO JHEP 09 (2012) 130 1207.5391
32 E. Re Single-top Wt-channel production matched with parton showers using the POWHEG method EPJC 71 (2011) 1547 1009.2450
33 P. Nason and G. Zanderighi W+W, WZ and ZZ production in the POWHEG -\textscbox-v2 EPJC 74 (2014) 2702 1311.1365
34 P. Artoisenet, R. Frederix, O. Mattelaer, and R. Rietkerk Automatic spin-entangled decays of heavy resonances in Monte Carlo simulations JHEP 03 (2013) 015 1212.3460
35 C. Bierlich et al. A comprehensive guide to the physics and usage of PYTHIA8.3 SciPost Phys. Codeb. 8 (2022) 2203.11601
36 CMS Collaboration Extraction and validation of a new set of CMS PYTHIA8 tunes from underlying-event measurements EPJC 80 (2020) 4 CMS-GEN-17-001
1903.12179
37 NNPDF Collaboration Parton distributions from high-precision collider data EPJC 77 (2017) 663 1706.00428
38 GEANT4 Collaboration GEANT 4---a simulation toolkit NIM A 506 (2003) 250
39 CMS Collaboration Performance of the CMS muon detector and muon reconstruction with proton-proton collisions at s= 13 TeV JINST 13 (2018) P06015 CMS-MUO-16-001
1804.04528
40 CMS Collaboration Electron and photon reconstruction and identification with the CMS experiment at the CERN LHC JINST 16 (2021) P05014 CMS-EGM-17-001
2012.06888
41 CMS Collaboration ECAL 2016 refined calibration and \mboxRun 2 summary plots CMS Detector Performance Note CMS-DP-2020-021, 2020
CDS
42 CMS Collaboration High-precision measurement of the W boson mass with the CMS experiment at the LHC Submitted to Nature, 2024 CMS-SMP-23-002
2412.13872
43 S. Camarda et al. \textscDYTurbo: Fast predictions for Drell--Yan processes EPJC 80 (2020) 251 1910.07049
44 S. Camarda, L. Cieri, and G. Ferrera Drell--Yan lepton-pair production: qT resummation at N\textsuperscript3LL accuracy and fiducial cross sections at N\textsuperscript3LO PRD 104 (2021) L111503 2103.04974
45 S. Camarda, L. Cieri, and G. Ferrera Fiducial perturbative power corrections within the qT subtraction formalism EPJC 82 (2022) 575 2111.14509
46 S. Camarda, L. Cieri, and G. Ferrera Drell--Yan lepton-pair production: qT resummation at N\textsuperscript4LL accuracy PLB 845 (2023) 138125 2303.12781
47 CMS Collaboration Performance of the CMS muon trigger system in proton-proton collisions at s= 13 TeV JINST 16 (2021) P07001 CMS-MUO-19-001
2102.04790
48 M. J. Oreglia A study of the reactions ψγγψ PhD thesis, Stanford University, . SLAC-R-236, 1980
link
49 J. E. Gaiser Charmonium spectroscopy from radiative decays of the J/ψ and ψ PhD thesis, Stanford University, . SLAC-R-255, 1982
link
50 CMS Collaboration The CMS statistical analysis and combination tool: \textsccombine Comput. Softw. Big Sci. 8 (2024) 19 CMS-CAT-23-001
2404.06614
51 W. Verkerke and D. Kirkby The \textscRooFit toolkit for data modeling in th International Conference on Computing in High Energy and Nuclear Physics (CHEP ): La Jolla CA, United States, March 24--28, . . . [eConf C0303241 MOLT007], 2003
Proc. 1 (2003) 3
physics/0306116
52 L. Moneta et al. The \textscRooStats project in th International Workshop on Advanced Computing and Analysis Techniques in Physics Research (ACAT ): Jaipur, India, February 22--27, . . . [PoS (ACAT) 057], 2010
Proc. 1 (2010) 3
1009.1003
53 CMS Collaboration First measurement of the top quark pair production cross section in proton-proton collisions at s= 13.6 TeV JHEP 08 (2023) 204 CMS-TOP-22-012
2303.10680
54 J. Campbell, T. Neumann, and Z. Sullivan Single-top-quark production in the t-channel at NNLO JHEP 02 (2021) 040 2012.01574
55 J. S. Conway Incorporating nuisance parameters in likelihoods for multisource spectra in orkshop on Statistical Issues Related to Discovery Claims in Search Experiments and Unfolding (PHYSTAT ): Geneva, Switzerland, January 17--20,, 2011
Proc. 2011 (2011) W
1103.0354
56 R. Barlow and C. Beeston Fitting using finite Monte Carlo samples Comput. Phys. Commun. 77 (1993) 219
57 CMS Collaboration Precision luminosity measurement in proton-proton collisions at s= 13 TeV in 2015 and 2016 at CMS EPJC 81 (2021) 800 CMS-LUM-17-003
2104.01927
58 CMS Collaboration Luminosity measurement in proton-proton collisions at 13.6 TeV in 2022 at CMS CMS Physics Analysis Summary, 2024
CMS-PAS-LUM-22-001
CMS-PAS-LUM-22-001
59 T.-J. Hou et al. New CTEQ global analysis of quantum chromodynamics with high-precision data from the LHC PRD 103 (2021) 014013 1912.10053
60 S. Bailey et al. Parton distributions from LHC, HERA, Tevatron and fixed target data: MSHT20 PDFs EPJC 81 (2021) 341 2012.04684
Compact Muon Solenoid
LHC, CERN