1 Introduction

We report a measurement of the effective leptonic weak mixing angle (\(\sin ^2\theta ^{\ell }_{\text {eff}}\)) using the forward–backward asymmetry (\(A_\text {FB}\)) in Drell–Yan \(\mathrm{q} \overline{\mathrm{q}} \rightarrow \ell ^+\ell ^-\) events, where \(\ell \) stands for muon (\(\mu \)) or electron (\(\mathrm {e}\)). The analysis is based on data from the CMS experiment at the CERN LHC. At leading order (LO), lepton pairs are produced through the annihilation of a quark with its antiquark into a \(\mathrm{Z}\) boson or a virtual photon: \(\mathrm{q} \overline{\mathrm{q}} \rightarrow \mathrm{Z}/\gamma \rightarrow \ell ^+\ell ^-\). For a given dilepton invariant mass \(m_{\ell \ell }\), the differential cross section at LO can be expressed at the parton level as

$$\begin{aligned} \frac{\mathrm{d}\sigma }{\mathrm{d}(\cos \theta ^{*})} \propto 1 + \cos ^{2} \theta ^{*} + A_{4} \cos \theta ^{*}, \end{aligned}$$
(1)

where the \((1 + \cos ^{2} \theta ^{*})\) term arises from the spin-1 of the exchanged boson, and the \(\cos \theta ^{*}\) term originates from interference between vector and axial-vector contributions. The definition of \(A_\text {FB}\) is based on the angle \(\theta ^*\) of the negative lepton (\(\ell ^-\)) in the Collins–Soper [1] frame of the dilepton system:

$$\begin{aligned} A_\text {FB} =\frac{3}{8}A_4 =\frac{\sigma _\mathrm {F}-\sigma _\mathrm {B}}{\sigma _\mathrm {F} +\sigma _\mathrm {B}}, \end{aligned}$$
(2)

where \(\sigma _\mathrm {F} \) and \(\sigma _\mathrm {B} \) are, respectively, the cross sections in the forward (\(\cos \theta ^{*} >0\)) and backward (\(\cos \theta ^{*} <0\)) hemispheres. In this frame, \(\theta ^*\) is the angle of the \(\ell ^-\) relative to the axis that bisects the angle between the direction of the quark and the reversed direction of the antiquark. In proton–proton (\(\mathrm {p}\mathrm {p}\)) collisions, the direction of the quark is more likely to be in the direction of the Lorentz boost of the dilepton. Therefore, \(\cos \theta ^{*}\) can be calculated using the following variables in the laboratory frame:

$$\begin{aligned} \cos \theta ^{*} =\frac{2(P_1^+P_2^- - P_1^-P_2^+)}{\sqrt{m_{\ell \ell } ^2(m_{\ell \ell } ^2+p_{\mathrm {T},\ell \ell } ^2)}}\,\frac{p_{z,\ell \ell }}{|p_{z,\ell \ell } |}, \end{aligned}$$
(3)

where \(m_{\ell \ell }\), \(p_{\mathrm {T},\ell \ell }\), and \(p_{z,\ell \ell }\) are the mass, transverse momentum, and longitudinal momentum, respectively, of the dilepton system, and the \(P_i^\pm \) are defined in terms of the energies (\(E_i\)) and longitudinal momenta (\(p_{z, i}\)), of the negatively and positively charged leptons as \(P_i^\pm =(E_i\pm p_{z,i})/\sqrt{2}\) [1].

Fig. 1
figure 1

The dependence of \(A_\text {FB} \) on \(m_{\ell \ell } \) in dimuon events generated using pythia  8.212 [16] and the LO NNPDF3.0 [17] PDFs for dimuon rapidities of \(|y_{\ell \ell } | <2.4\). The distributions for the total production (\(\mathrm{q} \overline{\mathrm{q}} \)) and the different channels are given on the left, overlaid with results based on Eq. (6), using the definition of \(A_\text {FB} ^\text {true}(m_{\ell \ell })\) for the known quark direction. The middle panel gives the diluted \(A_\text {FB}\) using instead the direction of the dilepton boost, and the right panel shows the diluted \(A_\text {FB}\) in \(|y_{\ell \ell } | \) bins of 0.4 for all channels

A non-zero \(A_\text {FB}\) value in dilepton events arises from the vector and axial-vector couplings of electroweak bosons to fermions. At LO, these respective couplings of \(\mathrm{Z}\) bosons to fermions (f) can be expressed as:

$$\begin{aligned} v_\mathrm {f}&= T_3^\mathrm {f}-2Q_\mathrm {f} \sin ^2\theta _\mathrm {W}, \end{aligned}$$
(4)
$$\begin{aligned} a_\mathrm {f}&= T_3^\mathrm {f}, \end{aligned}$$
(5)

where \(Q_\mathrm {f} \) and \(T_3^\mathrm {f} \) are the charge and the third component of the weak isospin of the fermion, respectively, and \(\sin ^2\theta _\mathrm {W}\) refers to the weak mixing angle, which is related to the masses of the \(\mathrm {W}\) and \(\mathrm{Z}\) bosons through the relation \(\sin ^2\theta _\mathrm{W} =1-m_\mathrm {W}^2/m_\mathrm{Z} ^2\). Electroweak (EW) radiative corrections affect these LO relations. In the improved Born approximation [2, 3], some of the higher-order corrections are absorbed into an effective mixing angle. The effective weak mixing angle is based on the relation \(v_\mathrm {f}/a_\mathrm {f} =1-4|Q_\mathrm {f} |\sin ^2\theta _{\text {eff}}^\mathrm {f} \), with \(\sin ^2\theta _\text {eff} ^\mathrm {f} =\kappa _\mathrm {f} \sin ^2\theta _\mathrm {W}\), where the flavor-dependent \(\kappa _\mathrm {f} \) is determined through EW corrections. The \(A_\text {FB}\) for dilepton events is sensitive primarily to \(\sin ^2\theta ^{\ell }_{\text {eff}}\).

We measure \(\sin ^2\theta ^{\ell }_{\text {eff}}\) by fitting the mass and rapidity \((y_{\ell \ell })\) dependence of the observed \(A_\text {FB}\) in dilepton events to standard model (SM) predictions as a function of \(\sin ^2\theta ^{\ell }_{\text {eff}}\). The most precise previous measurements of \(\sin ^2\theta ^{\ell }_{\text {eff}}\) were performed by the combined LEP and SLD experiments [4]. There is, however, a known discrepancy of about 3 standard deviations between the two most precise values. Other measurements of \(\sin ^2\theta ^{\ell }_{\text {eff}}\) have also been reported by the Tevatron and LHC experiments [5,6,7,8,9,10,11,12,13,14,15].

Using the LO expressions for the \(\mathrm{Z}\) boson, virtual photon exchange, and their interference, the “true” \(A_\text {FB} \) (i.e., using the quark direction in the definition of \(\cos \theta ^{*}\)) can be evaluated as

$$\begin{aligned} A_\text {FB} ^{\text {true}}(m_{\ell \ell })= & {} a_\ell a_\mathrm{q}(8v_\ell v_\mathrm{q}- Q_\mathrm{q}KD_m)\nonumber \\&\times [16(v_\ell ^2+a_\ell ^2)(v_\mathrm{q}^2+a_\mathrm{q}^2)-8 v_\ell v_\mathrm{q}Q_\mathrm{q}KD_m\nonumber \\&+ Q_\mathrm{q}^2K^2(D_m^2+\Gamma ^2_\mathrm{Z}/m_\mathrm{Z} ^2)]^{-1}, \end{aligned}$$
(6)

where the subscript \(\mathrm{q}\) refers to the participating quark, \(K=8\sqrt{2}\pi \alpha /G_\mathrm {F} m_\mathrm{Z} ^2\), \(D_m=1-m_\mathrm{Z} ^2/m_{\ell \ell } ^2\), \(\alpha \) is the electromagnetic coupling, \(G_\mathrm {F} \) is the Fermi constant, and \(\Gamma _\mathrm{Z} \) is the full decay width of the \(\mathrm{Z}\) boson. A strong dependence of \(A_\text {FB}\) on \(m_{\ell \ell }\) originates from axial and vector interference. The \(A_\text {FB}\) is negative at small \(m_{\ell \ell }\) and positive at large values, crossing \(A_\text {FB} =0\) slightly below the \(\mathrm{Z}\) boson peak.

In collisions of hadrons, \(A_\text {FB}\) is sensitive to parton distribution functions (PDFs) for two reasons. First, the different couplings of \(\mathrm{u}\)- and \(\mathrm{d}\)-type quarks to EW bosons generate different \(A_\text {FB}\) values in the corresponding production channels, which means that the average depends on the relative contributions of \(\mathrm{u}\)- and \(\mathrm{d}\)-type quarks to the total cross section. Second, the definition of \(A_\text {FB}\) in \(\mathrm {p}\mathrm {p}\) collisions is based on the sign of \(y_{\ell \ell }\), which relies on the fact that on average the dilepton pairs are Lorentz-boosted in the quark direction. Therefore, a non-zero average \(A_\text {FB}\) originates only from valence-quark production channels and is diluted by events where the antiquark carries a larger momentum than the quark. A dependence of the “true” and diluted \(A_\text {FB}\) on dilepton mass for different \(\mathrm{q} \overline{\mathrm{q}} \) production channels and their sum is shown in Fig. 1.

The dilution of \(A_\text {FB}\) depends strongly on \(y_{\ell \ell }\), as shown in Fig. 1. At zero rapidity, the quark and antiquark carry equal momenta, and the dilution is maximal, resulting in \(A_\text {FB} =0\). The \(A_\text {FB}\) is measured in 12 bins of dilepton mass, covering the range \(60<m_{\ell \ell } <120\,\text {GeV} \), and 6 \(|y_{\ell \ell } | \) bins of equal size for \(|y_{\ell \ell } | <2.4\). The boundaries in the dilepton mass are at: 60, 70, 78, 84, 87, 89, 91, 93, 95, 98, 104, 112, and 120\(\,\text {GeV}\). The mass bins are chosen such that near \(m_\mathrm{Z} \) the bin widths are larger than the mass resolution in any of the ranges of \(y_{\ell \ell }\). Smaller and larger mass bins are chosen such that all mass bins contain enough events to perform a meaningful independent measurement. The weak dependence of \(A_\text {FB}\) on \(p_{\mathrm {T},\ell \ell }\) is included in the SM predictions. The uncertainty originating from modeling of \(p_{\mathrm {T},\ell \ell }\) is very small and included in the theoretical estimates.

2 The CMS detector

The central feature of the CMS apparatus is a superconducting solenoid of 6\(\text {\,m}\) internal diameter, providing a magnetic field of 3.8\(\text {\,T}\). 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 reside within the solenoid volume. Forward calorimeters extend the pseudorapidity \(\eta \) coverage provided by the barrel and endcap detectors. Muons are measured in gas-ionization detectors embedded in the steel flux-return yoke outside the solenoid. A more detailed description of the CMS detector can be found in Ref. [18].

Muons are measured in the range \(|\eta | < 2.4\), using detection planes based on the drift-tube, cathode-strip chamber, or resistive-plate chamber technologies. Matching muons to tracks measured in the silicon tracker provides a relative transverse momentum resolution for muons with \(20<p_{\mathrm {T}} < 100\,\text {GeV} \) of 1.3–2.0% in the barrel, and less than 6% in the endcaps. The \(p_{\mathrm {T}}\) resolution in the barrel is smaller than 10% for muons with \(p_{\mathrm {T}}\) up to 1\(\,\text {TeV}\)  [19].

The electromagnetic calorimeter consists of 75 848 lead tungstate crystals that provide a coverage of \(|\eta | < 1.48 \) in the barrel region and \(1.48< |\eta | < 3.00\) in the two endcap regions. Preshower detectors consisting of two planes of silicon sensors, interleaved with a total of 3 radiation lengths of lead, are located in front of each endcap detector. The electron momentum is obtained by combining the energy measurement in the ECAL with that in the tracker. The momentum resolution for electrons with \(p_{\mathrm {T}} \approx 45\,\text {GeV} \) from \(\mathrm{Z} \rightarrow \mathrm {e}\mathrm {e}\) decays, ranges from 1.7% for nonshowering electrons in the barrel region, to 4.5% for showering electrons in the endcaps [20].

Events of interest are selected using a two-tiered trigger system [21]. The first level, consisting of custom hardware processors, uses information from the calorimeters and muon detectors to select events at a rate of about 100\(\text {\,kHz}\) within a time interval of less than 4\(\,\mu \text {s}\). 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, that reduces the event rate to about 1\(\text {\,kHz}\) before data storage.

3 Data and simulated events

The measurement is based on \(\mathrm {p}\mathrm {p}\) collisions at \(\sqrt{s}=8\,\text {TeV} \) recorded by the CMS Experiment in 2012, corresponding to integrated luminosities of 18.8 and \(19.6{{\,\text {fb}^{-1}}} \) for muon and electron channels, respectively.

Candidates for the dimuon channel are collected using an isolated single-muon trigger with a \(p_{\mathrm {T}}\) threshold of 24\(\,\text {GeV}\) and \(|\eta | <2.4\). At the beginning of data taking, the muon trigger was restricted to \(|\eta | <2.1\). We do not use these events, and the integrated luminosity in the dimuon analysis is therefore somewhat smaller than for dielectrons. Background contamination is reduced by applying identification and isolation criteria to the reconstructed muons. First, muon tracks are required to be reconstructed independently in the inner tracker and in the outer muon detectors. A global fit to the momentum, including both tracker and muon detector hits, must have a fitted \(\chi ^2/\text {dof}<10\), where dof stands for the degrees of freedom. Muon tracks are required to pass within a transverse distance of \(0.2\,\text {cm} \) from the primary vertex, defined as the \(\mathrm {p}\mathrm {p}\) vertex with the largest \(\sum p_{\mathrm {T}} ^2\) of its associated tracks. Muon candidates are rejected if the scalar-\(p_{\mathrm {T}}\) sum of all tracks within a cone of \(\Delta R=\sqrt{\smash [b]{(\Delta \eta )^2+(\Delta \phi )^2}}=0.3\) around the muon is larger than 10% of the \(p_{\mathrm {T}}\) of the muon (this is referred to as track isolation, with \(\phi \) being the azimuth in radians). The track isolation requirement is insensitive to contributions from additional soft \(\mathrm {p}\mathrm {p}\) interactions (pileup). An event is selected when there are at least two isolated muons, with the leading muon (i.e., the one with largest \(p_{\mathrm {T}}\)) having \(p_{\mathrm {T}} >25\,\text {GeV} \), and the next-to-leading muon having \(p_{\mathrm {T}} >15\,\text {GeV} \). At least one muon with \(p_{\mathrm {T}} >25\,\text {GeV} \) is required to trigger the event. For the Drell–Yan signal, the two leptons are required to have opposite sign (OS).

Dielectron candidates are collected using a single-electron trigger with a \(p_{\mathrm {T}}\) threshold of 27\(\,\text {GeV}\) and \(|\eta | <2.5\). Variables pertaining to the energy distribution in electromagnetic showers and to impact parameters of inner tracks are used to separate prompt electrons from electrons originating from photon conversions in detector material. The jet background from SM events produced through quantum chromodynamics (QCD) is referred to as multijet production. A particle-flow (PF) event reconstruction algorithm is used to identify different particle types (photons, electrons, muons, and charged and neutral hadrons [22]). The scalar-\(p_{\mathrm {T}}\) sum of all PF particles in a cone of \(\Delta R<0.3\) around the electron direction is required to be less than 15% of the electron \(p_{\mathrm {T}}\), which reduces the background from hadrons in multijet events that are reconstructed incorrectly as electrons. This sum is corrected for contributions from pileup [20]. The electron momentum is evaluated by combining the energy in the ECAL with the momentum in the tracker. To ensure good reconstruction, the coverage is restricted to \(|\eta | <2.4\), excluding the transition region of \(1.44<|\eta | <1.57\) between the ECAL barrel and endcap detectors, as electron reconstruction in this region is not optimal. Dielectron candidates are selected when at least two OS electrons pass all quality requirements. The leading and next-to-leading electrons must have respectively \(p_{\mathrm {T}} >30\) and \(>20\,\text {GeV} \), with the triggering electron always required to have \(p_{\mathrm {T}} >30\,\text {GeV} \).

A total of about 8.2 million dimuon and 4.9 million dielectron candidate events are selected for further analysis. The number of dielectron events is smaller because of the higher \(p_{\mathrm {T}}\) thresholds and more stringent selection criteria implemented in electron selections. The \(\mathrm{Z}/\gamma \rightarrow \mathrm {\mu ^+}\mathrm {\mu ^-} \) and \(\mathrm{Z}/\gamma \rightarrow \mathrm {e}^+\mathrm {e}^- \) data include small (\({<}1\%\)) background contaminations that originate from \(\mathrm{Z}/\gamma \rightarrow \mathrm {\tau }^{+}\mathrm {\tau }^{-} \), \({\mathrm{t}\overline{\mathrm{t}}}\), single top quark, and diboson (\(\mathrm {W}\) \(\mathrm {W}\), \(\mathrm {W}\) \(\mathrm{Z}\), and \(\mathrm{Z}\) \(\mathrm{Z}\)) events, as well as multijet and \(\mathrm {W}\) \(+\)jets events. Contributions from these backgrounds are subtracted from data as described below. Contamination from photon-induced background near the \(\mathrm{Z}\) boson peak is negligible [23].

Monte Carlo (MC) simulation is used to model signal and background processes. The signal as well as the single-boson and top quark backgrounds are based on next-to-leading order (NLO) matrix elements implemented in the powheg  v1 event generator [24,25,26,27] using the CT10 [28] PDFs. The generator is interfaced to pythia  6.426 [29] using the Z2* [30, 31] underlying event tune, which generates the parton showering, the hadronization, and the electromagnetic final-state radiation (FSR). The background events from \(\tau \) lepton decays are simulated with tauola 2.7 [32]. Diboson and multijet background events are generated with pythia 6 using the CTEQ6L1 PDFs [33]. Simulated minimum-bias events are superimposed on the hard-interaction events to model the effects from pileup. The detector response to all particles is simulated through Geant4  [34], and all final-state objects are reconstructed using the same algorithms used for data.

4 Corrections and backgrounds

The MC simulations are corrected to improve the modeling of the data. First, weight factors are applied to all simulated events to match the pileup distribution in data, which consists of roughly 20 interactions per crossing. These weights are based on the measured instantaneous luminosity and the total inelastic cross section that provides a good description of the average number of reconstructed vertices.

The total lepton-selection efficiency is factorized into the product of reconstruction, identification, isolation, and trigger efficiencies, with each component measured in samples of \(\mathrm{Z}/\gamma \rightarrow \ell ^+\ell ^-\) events through a “tag-and-probe” method [19, 20], in bins of lepton \(p_{\mathrm {T}}\) and \(\eta \). A charge-dependent efficiency in the muon triggering and reconstruction was observed in previous CMS measurements [35]. In the muon channel, all efficiencies are therefore determined separately for positively and negatively charged muons. The same procedures are used for data as for the simulated events, and scale factors are extracted to match the simulated event-selection efficiencies to those in the data.

The lepton momentum is calibrated using \(\mathrm{Z}/\gamma \rightarrow \ell ^+\ell ^-\) events [36]. The dominant sources of the mismeasurement of muon momentum originate from the mismodeling of tracker alignment and of the magnetic field. The correction parameters are obtained in bins of muon \(\eta \) and \(\phi \). First, the average \(1/p_{\mathrm {T}} \) values of the reconstructed muon curvature in data and simulation are corrected to the corresponding values calculated for MC generated muons. Then, using MC simulation, the resolution in the reconstructed muon momentum is parametrized as a function of the muon \(p_{\mathrm {T}}\) in bins of muon \(|\eta | \) and the number of tracker hits used in the reconstruction. Next, the correction parameters of the muon momentum scale are fine-tuned by matching the average dimuon mass in each bin of muon charge, \(\eta \), and \(\phi \) to their reference values. At this point, the “reference” distributions, which are based on the generated muons, are smeared by the reconstruction resolution derived in the previous step. Finally, the scale factors for the muon momentum resolution, in bins of muon \(|\eta | \), are determined by fitting the “reference” dimuon mass distribution to data.

A similar procedure is followed for electrons to reduce the small residual difference between the data and MC simulation. Unlike for muons, the measured electron energy is dominated by the calorimeter, and the corrections are extracted identically for electrons and positrons. The electron energy-scale parameters are fine-tuned by correcting the average dielectron mass in each bin of electron \(\eta \) and \(\phi \) to the corresponding “reference” values. Here, the “reference” distributions are based on the generated electrons (post FSR), combined with the FSR photons in a cone, and smeared by the reconstructed energy resolution.

The EW and top quark backgrounds are estimated using MC simulations based on the cross sections calculated at next-to-the-next-to-leading order in QCD [37, 38] and normalized to the integrated luminosity. We use cross sections calculated at NLO for the diboson backgrounds. The multijet background in dimuon events, dominated by muons from heavy-flavor hadron decays, is evaluated using same-sign (SS) dimuon events. A small EW and top quark contamination is evaluated in an MC simulation and subtracted from the SS sample. The distributions are then scaled by roughly a factor of 2, estimated from simulated events, to obtain the multijet contamination in the signal OS dimuon sample. The multijet background in the dielectron analysis is evaluated using the SS sample in combination with the \(\mathrm {e}\mu \) events to subtract the contribution from the OS events caused by the misidentification of charge. The distributions used to estimate the background from jets misidentified as leptons (that include the multijet and \(\mathrm {W}\)+jet events) are obtained from the SS \(\mathrm {e}\mu \) sample. These distributions are used to fit the dielectron mass distribution in the SS events in each \(y_{\ell \ell } \) bin to extract the normalization of this background.

Fig. 2
figure 2

Dimuon (left) and dielectron (right) mass distributions in three representative bins in rapidity: \(|y_{\ell \ell } | <0.4\) (upper), \(0.8<|y_{\ell \ell } | <1.2\) (middle), and \(1.6<|y_{\ell \ell } | <2.0\) (lower)

Fig. 3
figure 3

The muon (left) and electron (right) \(\cos \theta ^{*}\) distributions in three representative bins in rapidity: \(|y_{\ell \ell } | <0.4\) (upper), \(0.8<|y_{\ell \ell } | <1.2\) (middle), and \(1.6<|y_{\ell \ell } | <2.0\) (lower). The small contributions from backgounds are included in the predictions

The dilepton mass and \(\cos \theta ^{*}\) distributions in three of the six rapidity bins are shown in Figs. 2 and 3, respectively. The figures include lepton momentum and efficiency corrections, background samples normalized as described above, and the signal normalized to the total expected number of events in the data.

5 Weighted \(A_\text {FB}\) measurement

As introduced in Sect. 1, the LO angular distribution of dilepton events has a \((1+\cos ^2\theta ^*)\) term that arises from the spin-1 of the exchanged boson and a \(\cos \theta ^{*}\) term that originates from the interference between vector and axial-vector contributions. However, there is also a \((1-3\cos ^2\theta ^*)\) NLO term that originates from the \(p_{\mathrm {T}}\) of the interacting partons [39]. Each \((m_{\ell \ell },y_{\ell \ell })\) bin of the dilepton pair at NLO therefore has an angular distribution in \(\cos \theta ^{*}\) that follows the form [39]:

$$\begin{aligned} \frac{1}{\sigma }\frac{\mathrm{d}\sigma }{\mathrm{d}\cos \theta ^{*}} = \frac{3}{8}\left[ 1+\cos ^2\theta ^*+\frac{A_0}{2}(1-3\cos ^2\theta ^*) + A_4\cos \theta ^{*} \right] . \end{aligned}$$
(7)

The \(A_\text {FB} \) value in each \((m_{\ell \ell },y_{\ell \ell })\) bin is calculated using the “angular event weighting” method, described in Ref. [40], in which each event with a \(\cos \theta ^{*}\) value (denoted as “c”), is reflected in the denominator (D) and numerator (N) weights through:

$$\begin{aligned} w_\mathrm {D}= & {} \frac{1}{2}\frac{c^2}{(1+c^2+h)^3}, \end{aligned}$$
(8)
$$\begin{aligned} w_\mathrm {N}= & {} \frac{1}{2}\frac{|c |}{(1+c^2+h)^2}, \end{aligned}$$
(9)

where \(h=0.5A_0(1-3c^2)\). Here, as a baseline we use the \(p_{\mathrm {T},\ell \ell }\)-averaged \(A_0\) value of about 0.1 in each measurement \((m_{\ell \ell },y_{\ell \ell })\) bin, as predicted by the signal MC simulation. Using the weighted sums N and D for forward (\(\cos \theta ^{*} >0\)) and backward (\(\cos \theta ^{*} <0\)) events, we obtain

$$\begin{aligned}&D_\mathrm {F} =\sum _{c>0}w_\mathrm {D}, \quad D_\mathrm {B} =\sum _{c<0}w_\mathrm {D}, \end{aligned}$$
(10)
$$\begin{aligned}&N_\mathrm {F} =\sum _{c>0}w_\mathrm {N}, \quad N_\mathrm {B} =\sum _{c<0}w_\mathrm {N}, \end{aligned}$$
(11)

from which the weighted \(A_\text {FB}\) of Eq. (2) can be written as:

$$\begin{aligned} A_\text {FB} =\frac{3}{8}\frac{N_\mathrm {F}-N_\mathrm {B}}{D_\mathrm {F} +D_\mathrm {B}} . \end{aligned}$$
(12)

The statistical uncertainty in this weighted \(A_\text {FB}\) value takes into account correlations among the numerator and denominator sums. For data, the background contribution in the event-weighted sums are subtracted before calculating \(A_\text {FB}\). In the full phase space, the values of the weighted and the nominal \(A_\text {FB}\), calculated as an asymmetry between the total event counts in the forward and backward hemispheres, are the same. Since the acceptances of the forward and backward events are equal for same values of \(|\cos \theta ^{*} |\), the fiducial values of the event-weighted \(A_\text {FB}\) are also the same as in the full phase space, while the nominal \(A_\text {FB}\) values are smaller because of the limited acceptance at large \(\cos \theta ^{*}\). This feature makes an event-weighted \(A_\text {FB}\) less sensitive than the nominal \(A_\text {FB}\) to the specific modeling of the acceptance. In addition, because the event-weighted \(A_\text {FB}\) exploits the full distribution in \(\cos \theta ^{*}\), as opposed to only its sign in the nominal \(A_\text {FB}\), it therefore provides a smaller statistical uncertainty.

6 Extraction of \(\sin ^2\theta ^{\ell }_{\text {eff}}\)

We extract \(\sin ^2\theta ^{\ell }_{\text {eff}}\) by fitting the \(A_\text {FB}\) \((m_{\ell \ell },y_{\ell \ell })\) distribution in data with the theoretical predictions. The default signal distributions are based on the powheg  v2 event generator using the NNPDF3.0 PDFs [17]. The powheg generator is interfaced with pythia  8 [16] and the CUETP8M1 [31] underlying event tune to provide parton showering and hadronization, including electromagnetic FSR. The dependence on \(\sin ^2\theta ^{\ell }_{\text {eff}}\), on the renormalization and factorization scales, and on the PDFs is modeled through the powheg MC generator that provides matrix-element-based, event-by-event weights for each change in these parameters. The distributions are modified to different values of \(\sin ^2\theta ^{\ell }_{\text {eff}}\) by weighting each event in the full simulation by the ratio of \(\cos \theta ^{*}\) distributions obtained with the modified and default configurations in each \((m_{\ell \ell },y_{\ell \ell })\) bin. The uncertainties in the simulation of the detector have a small effect because \(A_\text {FB}\) is extracted through the angular event-weighting technique that is insensitive to efficiency and acceptance.

Table 1 Summary of statistical uncertainties in \(\sin ^2\theta ^{\ell }_{\text {eff}}\). The statistical uncertainties in the lepton-selection efficiency and in the calibration coefficients in data are included in the estimates
Fig. 4
figure 4

Comparison between data and best-fit \(A_\text {FB}\) distributions in the dimuon (upper) and dielectron (lower) channels. The best-fit \(A_\text {FB}\) value in each bin is obtained via linear interpolation between two neighboring templates. Here, the templates are based on the central prediction of the NLO NNPDF3.0 PDFs. The error bars represent the statistical uncertainties in the data

Table 1 summarizes the statistical uncertainty in the extracted \(\sin ^2\theta ^{\ell }_{\text {eff}}\) in the muon and electron channels and in their combination. Comparisons between the data and best-fit distributions are shown in Fig. 4. The statistical uncertainties are evaluated through the bootstrapping technique [41], and take account of correlations among the measured \(A_\text {FB}\), lepton selection efficiencies, and calibration coefficients introduced through the repeated use of the same dilepton events. We generate 400 pseudo-experiments that provide an accurate estimate of the statistical uncertainties and correlations. In each pseudo-experiment, every event in the data is replicated n times, where n is a random number sampled from a Poisson distribution with a mean of unity. All steps of the analysis, including extraction of muon selection efficiencies, calibration coefficients, and a measurement of \(A_\text {FB}\), are performed for each pseudo-experiment. The statistical uncertainties in electron-selection efficiencies and calibration coefficients, which have no charge dependence, are small and are evaluated separately.

7 Experimental systematic uncertainties

The experimental sources of systematic uncertainty reflect the statistical uncertainties in the simulated events, corrections to lepton-selection efficiency, and to the lepton-momentum scale and resolution, background subtraction, and modeling of pileup. For electrons, the selection efficiencies, which have no dependence on charge, cancel to first order, since we are using the angular event-weighting technique.

7.1 Statistical uncertainties in MC simulated events

To reduce the statistical uncertainties associated with the limited number of events in the signal MC samples, which include simulation of detector response and lepton reconstruction, the generated \(\cos \theta ^{*}\) distributions in each \((m_{\ell \ell },y_{\ell \ell })\) bin within the acceptance of the detector is reweighted to much larger MC samples, generated without simulating detector response or lepton reconstruction. This makes the fluctuations in the generated \(\cos \theta ^{*}\) distributions negligible, and therefore the statistical uncertainties in the reconstructed \(A_\text {FB}\) values become dominated by fluctuations in the simulated detector response and lepton reconstruction. These uncertainties are evaluated using the bootstrapping [41] method in both dimuon and dielectron channels, described in Sect. 6, by reweighting the generated \(\cos \theta ^{*}\) distributions in each of the bootstrap samples. The total statistical uncertainties in the simulated events also include contributions from uncertainties in the measured lepton-selection efficiencies and calibration coefficients.

7.2 Lepton selection efficiencies

Several sources of uncertainty are considered in measuring of efficiencies. The statistical uncertainties in the lepton-selection efficiencies, evaluated through studies of pseudo-experiments, are included in the combined statistical uncertainty of the measured \(\sin ^2\theta ^{\ell }_{\text {eff}}\).

Combined scale factors for muon reconstruction, identification, and isolation efficiencies are changed by 0.5%, and trigger-selection efficiency scale factors by 0.2%, coherently for all bins for both positive and negative lepton charges. These take into account uncertainties associated with the tag-and-probe method, and are evaluated by changing signal and background models for dimuon mass distributions, levels of backgrounds, the dimuon mass range, and binning used in the fits. These uncertainties are considered fully correlated between the two charges, and therefore have a negligible impact on the measurement of \(\sin ^2\theta ^{\ell }_{\text {eff}}\). In addition, we assign the difference between the offline efficiencies obtained by fitting the dimuon mass distributions to extract the signal yields, and those found using simple counting method, as additional systematic uncertainties. The total systematic uncertainty in \(\sin ^2\theta ^{\ell }_{\text {eff}}\) originating from the muon selection efficiency is ±0.00005.

In a similar way as for muons, the scale factors for electron reconstruction, identification, and trigger-selection efficiencies are changed coherently within their uncertainties in all \((p_{\mathrm {T}},\eta )\) bins, and the corresponding changes in the resulting \(\sin ^2\theta ^{\ell }_{\text {eff}}\) are assigned as systematic uncertainties. The total uncertainty in \(\sin ^2\theta ^{\ell }_{\text {eff}}\) originating from all electron efficiency-related systematic sources is \(\pm 0.00004\).

7.3 Lepton momentum calibration

The statistical uncertainties in the parameters used to calibrate lepton momentum, described in Sect. 4, are included in the combined statistical uncertainty. The theoretical uncertainties, discussed in Sect. 8, are also propagated to the reference distributions used to extract the coefficients in the lepton momentum calibration.

When evaluating the average dimuon masses to extract the \((\eta ,\phi )\) dependent corrections, the dimuon mass window is restricted to \(86<m_{\mu \mu }<96\,\text {GeV} \). This range of \(\pm 5\,\text {GeV} \) centered at 91\(\,\text {GeV}\) is changed from \(\pm 2.5\) to \(\pm 10\,\text {GeV} \) in steps of 0.5\(\,\text {GeV}\), and the full calibration sequence is repeated each time. Similarly, a dimuon mass window of \(\pm 10\) (i.e., 81–101)\(\,\text {GeV}\), used in the dimuon fits to obtain the resolution-correction factors, is changed from \(\pm 5\) to \(\pm 25\,\text {GeV} \) in steps of 1\(\,\text {GeV}\). For each of these modifications, the maximum deviation in the extracted \(\sin ^2\theta ^{\ell }_{\text {eff}}\) relative to the nominal configuration is taken as a systematic uncertainty. The total experimental systematic uncertainty in \(\sin ^2\theta ^{\ell }_{\text {eff}}\) originating from the muon-momentum calibration, evaluated by adding individual uncertainties in quadrature, is ±0.00008. The effects due to PDF uncertainties in the calibration coefficients were found to be negligible. In studies of the impact of the value of \(\sin ^2\theta ^{\ell }_{\text {eff}}\) used to generate the reference distributions for muon-momentum calibration over the range of \(\Delta \sin ^2\theta ^{\ell }_{\text {eff}} =0.02000\), the extracted result changes at most by ±0.00008 due to the changes made in the muon-calibration parameters. Since the uncertainty in \(\sin ^2\theta ^{\ell }_{\text {eff}}\) is much smaller than \(\pm 0.02000\), we conclude that this effect is negligible.

Similarly, the windows in the dielectron invariant mass used to extract the electron momentum-correction factors are changed to estimate the corresponding systematic uncertainty. And consider additional independent sources of systematic uncertainty from the modeling of pileup, background estimation, and bias in the dielecton mass-fitting procedure. The size of the EW corrections in the extracted electron energy-calibration coefficients is estimated by modifying reference dielectron mass distributions through the weight factors obtained with zgrad [42]. All these systematic uncertainties are found to be rather small. The dominant uncertainty originates from the full corrections to the electron energy resolution, which improve the agreement between data and simulated dielectron mass distributions. The total systematic uncertainty in the extracted value of \(\sin ^2\theta ^{\ell }_{\text {eff}}\) due to both the electron energy scale and resolution is \(\pm 0.00019\).

7.4 Background

The systematic uncertainties in the estimated background are evaluated as follows. The normalizations of the top quark and \(\mathrm{Z}/\gamma \rightarrow \mathrm {\tau }^{+}\mathrm {\tau }^{-} \) backgrounds are changed respectively by 10 and 20%, covering the maximum deviations between the data and simulation observed in the \(\mathrm {e}\mu \) control region. The uncertainty in the multijet and \(\mathrm {W}\)+jets background is estimated by changing them by ±100%. Changing the diboson background prediction by 100% provides a negligible change in the result (\({<}0.00001\)). Changing all EW and top quark backgrounds by the uncertainty in the integrated luminosity of 2.6% [43] also produces a negligible change in the result (\({<}0.00001\)). The total systematic uncertainty in the measured \(\sin ^2\theta ^{\ell }_{\text {eff}}\) from the uncertainty in the background estimation is \(\pm 0.00003\) and \(\pm 0.00005\) in the dimuon and dielectron channels, respectively.

7.5 Pileup

To take into account the uncertainty originating from differences in pileup between data and simulation, we change the total inelastic cross section by \(\pm 5\%\), and recompute the expected pileup distribution in data. The analysis is repeated and the difference relative to the central value is taken as the systematic uncertainty. These uncertainties are respectively \(\pm 0.00003\) and \(\pm 0.00002\) in the dimuon and dielectron channels.

All the above systematic uncertainties are summarized in Table 2.

8 Theoretical systematic uncertainties

We investigate sources of systematic uncertainty in modeling the MC templates. For each change in the model, we rederive the reference distributions described in Sect. 4 to adjust the lepton momentum calibration coefficients. As a baseline, the signal MC events are weighted to match the \(p_{\mathrm {T},\ell \ell }\) distribution in each \(|y_{\ell \ell } | \) bin in the data. The difference relative to the result obtained without applying the weight factors, which is 0.00003 in both channels, is assigned as a systematic uncertainty associated with the modeling of \(p_{\mathrm {T},\ell \ell }\).

Table 2 Summary of experimental systematic uncertainties in \(\sin ^2\theta ^{\ell }_{\text {eff}}\)

The renormalization and factorization scales, \(\mu _\mathrm {R} \) and \(\mu _\mathrm {F} \), are each changed independently by a factor of 2, up and down, such that their ratio is within \(0.5<\mu _\mathrm {R}/\mu _\mathrm {F} <2.0\). The maximum deviation among these six variants relative to the nominal choice (excluding the two opposite changes) is assigned as a systematic uncertainty associated with the missing higher-order QCD correction terms.

In addition, we use a multi-scale improved NLO (MiNLO [44]) calculation for the \(\mathrm{Z}\)+1 jet partonic final state (henceforth referred to as “\(\mathrm{Z}\)+j”), interfaced with pythia  8 for parton showering, FSR, and hadronization, to assess the uncertainty from the missing higher-order QCD terms and modeling of the angular coefficients. The MiNLO \(\mathrm{Z}\)+j process has NLO accuracy for both \(\mathrm{Z}\)+0 and \(\mathrm{Z}\)+1 jet events, which provides a better description of the dependence of the angular coefficients on \(p_{\mathrm {T},\ell \ell }\).

Systematic uncertainties in modeling electromagnetic FSR are estimated by comparing results obtained with distributions based on pythia  8 and photos  2.15 [45,46,47] for the modeling of FSR. Electroweak effects from the difference between the \(\mathrm{u}\) and \(\mathrm{d}\) quarks and leptonic effective mixing angles, are estimated by changing \(\sin ^2\theta _\text {eff} ^\mathrm{u}\) and \(\sin ^2\theta _\text {eff} ^\mathrm{d}\) by 0.0001 and 0.0002 [42], respectively, relative to \(\sin ^2\theta ^{\ell }_{\text {eff}}\). The \(\sin ^2\theta ^{\ell }_{\text {eff}}\) extracted using the corresponding distributions is shifted by 0.00001.

The underlying event tune parameters [31] are changed by their uncertainties, and \(\sin ^2\theta ^{\ell }_{\text {eff}}\) is extracted also using the corresponding distributions. The maximum difference from the default tune is taken as the corresponding uncertainty. The systematic uncertainties from these and all the above sources, are summarized in Table 3.

We also separately study the modeling of the \(A_0\) angular coefficient, which is included in the definition of \(A_\text {FB}\). As a baseline, the \(p_{\mathrm {T},\ell \ell }\)-averaged \(A_0\) value in each measurement \((m_{\ell \ell },y_{\ell \ell })\) bin is used in the definition of the weighted \(A_\text {FB}\). Several other options are studied: (i) the LO expression: \(A_0=p_{\mathrm {T},\ell \ell } ^2/(p_{\mathrm {T},\ell \ell } ^2+m_{\ell \ell } ^2)\), (ii) the \(p_{\mathrm {T},\ell \ell } \)-dependent \(A_0\) in each \((m_{\ell \ell },y_{\ell \ell })\) bin as predicted in the baseline NLO powheg simulation, (iii) the \(p_{\mathrm {T},\ell \ell } \)-dependent \(A_0\) predicted in the MiNLO \(\mathrm{Z}\)+j powheg generator, and (iv) \(A_0\) set to 0. The same definition is used for data and simulation, and the extracted \(\sin ^2\theta ^{\ell }_{\text {eff}}\) is identical within \(\pm 0.00002\) of the default. In addition, we weight the \(|\cos \theta ^{*} |\) distribution from the MiNLO \(\mathrm{Z}\)+j MC sample to match the dependence of \(A_0\) on \(p_{\mathrm {T},\ell \ell }\) in each \((m_{\ell \ell },y_{\ell \ell })\) bin to the corresponding values of the baseline MC simulation. The change in the resulting \(\sin ^2\theta ^{\ell }_{\text {eff}}\) is also negligible.

Table 3 Summary of the theoretical uncertainties for the dimuon and dielectron channels, as discussed in the text

9 Uncertainties in the PDFs

The observed \(A_\text {FB}\) values depend on the size of the dilution effect, as well as on the relative contributions from \(\mathrm{u}\) and \(\mathrm{d}\) valence quarks to the total dilepton production cross section. The uncertainties in the PDFs translate into sizable changes in the observed \(A_\text {FB}\) values. However, changes in PDFs affect the \(A_\text {FB} (m_{\ell \ell },y_{\ell \ell })\) distribution in a different way than changes in \(\sin ^2\theta ^{\ell }_{\text {eff}}\).

Changes in PDFs produce large changes in \(A_\text {FB}\), when the absolute values of \(A_\text {FB}\) are large, i.e., at large and small dilepton mass values. In contrast, the effect of changes in \(\sin ^2\theta ^{\ell }_{\text {eff}}\) are largest near the \(\mathrm{Z}\) boson peak, and are significantly smaller at high and low masses. Because of this behavior, which is illustrated in Fig. 5, we apply a Bayesian \(\chi ^2\) reweighting method to constrain the PDFs [48,49,50], and thereby reduce their uncertainties in the extracted value of \(\sin ^2\theta ^{\ell }_{\text {eff}}\).

Fig. 5
figure 5

Distribution in \(A_\text {FB}\) as a function of dilepton mass, integrated over rapidity (top), and in six rapidity bins (bottom) for \(\sin ^2\theta ^{\ell }_{\text {eff}} =0.23120\) in powheg. The solid lines in the bottom panel correspond to six changes at \(\sin ^2\theta ^{\ell }_{\text {eff}} \) around the central value, corresponding to: \(\pm 0.00040\), \(\pm 0.00080\), and \(\pm 0.00120\). The dashed lines refer to the \(A_\text {FB}\) predictions for 100 NNPDF3.0 replicas. The shaded bands illustrate the standard deviation in the NNPDF3.0 replicas

Fig. 6
figure 6

The upper panel in each figure shows a scatter plot in \(\chi ^2_\text {min} \) vs. the best-fit \(\sin ^2\theta ^{\ell }_{\text {eff}}\) for 100 NNPDF replicas in the muon channel (upper left), electron channel (upper right), and their combination (below). The corresponding lower panels have the projected distributions in the best-fit \(\sin ^2\theta ^{\ell }_{\text {eff}}\) for the nominal (open circles) and weighted (solid circles) replicas

As a baseline, we use the NLO NNPDF3.0 PDFs. In the Bayesian \(\chi ^2\) reweighting method, PDF replicas that offer good descriptions of the observed \(A_\text {FB}\) distribution are assigned large weights, and those that poorly describe the \(A_\text {FB}\) are given small weights. Each weight factor is based on the best-fit \(\chi ^2_{\text {min},i}\) value obtained by fitting the \(A_\text {FB}\) (\(m_{\ell \ell }\),\(y_{\ell \ell }\)) distribution with a given PDF replica i:

$$\begin{aligned} w_i = \frac{\mathrm{e}^{-\frac{\chi ^2_{\text {min},i}}{2}}}{\frac{1}{N}\sum _{i=1}^N \mathrm{e}^{-\frac{\chi ^2_{\text {min},i}}{2}}}, \end{aligned}$$
(13)

where N is the number of replicas in a set of PDFs. The final result is then calculated as a weighted average over the replicas: \(\sin ^2\theta ^{\ell }_{\text {eff}} =\sum _{i=1}^{N} w_i s_i/N\), where \(s_i\) is the best-fit \(\sin ^2\theta ^{\ell }_{\text {eff}}\) value obtained for the ith replica.

Figure 6 shows a scatter plot of the \(\chi ^2_{\text {min}}\) vs. the best-fit \(\sin ^2\theta ^{\ell }_{\text {eff}}\) value for the 100 NNPDF3.0 replicas for the \(\mu \mu \) and \(\mathrm {e}\mathrm {e}\) samples, and for the combined dimuon and dielectron results. All sources of statistical and experimental systematic uncertainties are included in a \(72{\times }72\) covariance matrices for data and template \(A_\text {FB}\) distributions. The \(\chi ^2(s)\) is defined as:

$$\begin{aligned} \chi ^2(s)= (\varvec{D}-\varvec{T}(s))^{T}\varvec{V}^{-1}(\varvec{D}-\varvec{T}(s)), \end{aligned}$$
(14)

where \(\varvec{D}\) represents the measured \(A_\text {FB}\) values for data in 72 bins, \(\varvec{T}(s)\) denotes the theoretical predictions for \(A_\text {FB}\) as a function of s, or \(\sin ^2\theta ^{\ell }_{\text {eff}} \), and \(\varvec{V}\) represents the sum of the covariance matrices for the data and templates. As illustrated in these figures, the extreme PDF replicas from either side are disfavored by both the dimuon and dielectron data. For each of the NNPDF3.0 replicas, the muon and electron results are combined using their respective best-fit \(\chi ^2\) values, \(\sin ^2\theta ^{\ell }_{\text {eff}}\), and their fitted statistical and experimental systematic uncertainties.

Figure 7 shows the extracted \(\sin ^2\theta ^{\ell }_{\text {eff}}\) in the muon and electron decay channels and their combination, with and without constraining the uncertainties in the PDFs. The corresponding numerical values are also listed in Table 4. After Bayesian \(\chi ^2\) reweighting, the PDF uncertainties are reduced by about a factor of 2. It should be noted that the Bayesian \(\chi ^2\) reweighting technique works well when the replicas span the optimal value on both of its sides. In addition, the effective number of replicas after \(\chi ^2\) reweighting, \(n_\text {eff} =N^2/\sum _{i=1}^{N}w_i^2\), should also be large enough to give a reasonable estimate of the average value and its standard deviation. There are 39 effective replicas after the \(\chi ^2\) reweighting (\(n_\text {eff} =39\)). Including the corresponding statistical uncertainty of 0.00005, the total PDF uncertainty becomes 0.00031. As a cross-check, we perform the analysis with the corresponding set of 1000 NNPDF3.0 replicas in the dimuon channel, and find good consistency between the two results.

Fig. 7
figure 7

The extracted values of \(\sin ^2\theta ^{\ell }_{\text {eff}}\) in the muon and electron channels, and their combination. The horizontal bars include statistical, experimental, and PDF uncertainties. The PDF uncertainties are obtained both without (top) and with (bottom) using the Bayesian \(\chi ^2\) weighting

Table 4 The central value and the PDF uncertainty in the measured \(\sin ^2\theta ^{\ell }_{\text {eff}}\) in the muon and electron channels, and their combination, obtained without and with constraining PDFs using Bayesian \(\chi ^2\) reweighting
Fig. 8
figure 8

Extracted values of \(\sin ^2\theta ^{\ell }_{\text {eff}}\) from the dimuon data for different sets of PDFs with the nominal (top) and \(\chi ^2\)-reweighted (bottom) replicas. The horizontal error bars include contributions from statistical, experimental, and PDF uncertainties

Fig. 9
figure 9

Comparison of the measured \(\sin ^2\theta ^{\ell }_{\text {eff}}\) in the muon and electron channels and their combination, with previous LEP, SLD, Tevatron, and LHC measurements. The shaded band corresponds to the combination of the LEP and SLD measurements

We have also studied the PDFs represented by Hessian eigenvectors using the CT10 [28], CT14 [51], and MMHT2014 [52] PDFs in an analysis performed in the dimuon channel. First, we generate the replica predictions (i) for each observable O for the Hessian eigensets (k):

$$\begin{aligned} O_i = O_0 + \frac{1}{2} \sum _{k=0}^{n} (O_{2k+1}-O_{2k+2})R_{ik}, \end{aligned}$$
(15)

where n is the number of eigenvector axes, and the \(R_{ik}\) are random numbers sampled from the normal distribution with a mean of 0 and a standard deviation of unity. Then, the same technique is applied as used in the NNPDF analysis. The results of fits for these PDFs are summarized in Fig. 8. After Bayesian \(\chi ^2\) reweighting the central predictions for all PDFs are closer to each other, and the corresponding uncertainties are significantly reduced. The result using CT14 is within about 1/3 of the PDF uncertainty of the NNPDF3.0 result in the muon channel, whereas the MMHT2014 set yields a smaller \(\sin ^2\theta ^{\ell }_{\text {eff}}\) value by about one standard deviation. Some of these differences can be reduced by adding more data (e.g. including the electron channel, which is not considered in this check). Some can be attributed to the residual differences in the valence and sea quark distributions, which are not fully constrained using the \(A_\text {FB}\) distributions alone. For example, we find that the NLO NNPDF3.0 PDF set yields a very good description for the published 8 TeV CMS muon charge asymmetry (\(\chi ^2\) of 4.6 for 11 dof). In contrast, the \(\chi ^2\) values with the CT14 and MMHT2014 PDF sets are 21.3 and 21.4, respectively. We also constructed a combined set from same number of replicas of NNPDF3.0, CT14, and MMHT2014 PDFs, and after including the data from the \(\mathrm {W}\) charge asymmetry in the PDF reweighting, we find the combined weighted average in the dimuon channel differs from the NNPDF3.0 result by only 0.00009, and the standard deviation only increases from 0.00032 to 0.00036. Consequently, for our quoted results we use only the NNPDF3.0 PDF set, which is used in both dimuon and dielectron analyses.

As an additional test, for the case of Hessian PDFs (including the Hessian NNPDF3.0 [53]) we perform a simultaneous \(\chi ^2\) fit for \(\sin ^2\theta ^{\ell }_{\text {eff}}\) and all PDF nuisance parameters representing the variations for each eigenvector. As expected for Gaussian distributions, we obtain the same central values and the total uncertainties that are extracted from Bayesian reweighting of the corresponding set of replicas.

Finally, as a cross-check, we also repeat the measurement using different mass windows for extracting \(\sin ^2\theta ^{\ell }_{\text {eff}}\), and for constraining the PDFs. Specifically, we first use the central five bins, corresponding to the dimuon mass range of \(84<m_{\mu \mu }<95\,\text {GeV} \), to extract \(\sin ^2\theta ^{\ell }_{\text {eff}}\). Then, we use predictions based on the extracted \(\sin ^2\theta ^{\ell }_{\text {eff}}\) in the lower three \((60< m_{\mu \mu } <84\,\text {GeV})\) and the higher four \((95<m_{\mu \mu }<120\,\text {GeV})\) dimuon mass bins, to constrain the PDFs. We find that the statistical uncertainty increases by only about 10%, and the PDF uncertainty increases by only about 6% relative to the uncertainties obtained when using the full mass range to extract the \(\sin ^2\theta ^{\ell }_{\text {eff}}\) and simultaneously constrain the PDFs. The test thereby confirms that the PDF uncertainties are constrained mainly by the high- and low-mass bins, and that we obtain consistent results with these two approaches.

10 Summary

The effective leptonic mixing angle, \(\sin ^2\theta ^{\ell }_{\text {eff}}\), has been extracted from measurements of the mass and rapidity dependence of the forward–backward asymmetries \(A_\text {FB}\) in Drell–Yan \(\mu \mu \) and \(\mathrm {e}\mathrm {e}\) production. As a baseline model, we use the powheg event generator for the inclusive \(\mathrm {p}\mathrm {p}\rightarrow \mathrm{Z}/\gamma \rightarrow \ell \ell \) process at leading electroweak order, where the weak mixing angle is interpreted through the improved Born approximation as the effective angle incorporating higher-order corrections. With more data and new analysis techniques, including precise lepton-momentum calibration, angular event weighting, and additional constraints on PDFs, the statistical and systematic uncertainties are significantly reduced relative to previous CMS measurements. The combined result from the dielectron and dimuon channels is:

$$\begin{aligned} \sin ^2\theta ^{\ell }_{\text {eff}}= & {} 0.23101 \pm 0.00036\,\text {(stat)} \pm 0.00018\,\text {(syst)} \nonumber \\&\pm 0.00016\,\text {(theo)} \pm 0.00031\,(\text {PDF}), \end{aligned}$$
(16)

or summing the uncertainties in quadrature,

$$\begin{aligned} \sin ^2\theta ^{\ell }_{\text {eff}} =0.23101\pm 0.00053. \end{aligned}$$
(17)

A comparison of the extracted \(\sin ^2\theta ^{\ell }_{\text {eff}}\) with previous results from LEP, SLC, Tevatron, and LHC, shown in Fig. 9, indicates consistency with the mean of the most precise LEP and SLD results, as well as with the other measurements.