Measurement of $W^{\pm}$-boson and $Z$-boson production cross-sections in $pp$ collisions at $\sqrt{s}=2.76$ TeV with the ATLAS detector
ATLAS Collaboration

TL;DR
This paper reports measurements of W and Z boson production cross-sections at 2.76 TeV using ATLAS data, providing key insights into electroweak processes at this energy.
Contribution
First measurement of W and Z boson production cross-sections at 2.76 TeV with detailed analysis and comparison to theoretical predictions.
Findings
Total cross-sections for W+ and W- bosons measured with uncertainties.
Z boson production cross-section determined with precision.
Ratios and asymmetries analyzed to reduce systematic uncertainties.
Abstract
The production cross-sections for and bosons are measured using ATLAS data corresponding to an integrated luminosity of 4.0 pb collected at a centre-of-mass energy TeV. The decay channels and are used, where can be an electron or a muon. The cross-sections are presented for a fiducial region defined by the detector acceptance and are also extrapolated to the full phase space for the total inclusive production cross-section. The combined (average) total inclusive cross-sections for the electron and muon channels are: \begin{eqnarray} \sigma^{\text{tot}}_{W^{+}\rightarrow \ell \nu}& = & 2312 \pm 26\ (\text{stat.})\ \pm 27\ (\text{syst.}) \pm 72\ (\text{lumi.}) \pm 30\ (\text{extr.})\text{pb} \nonumber, \\ \sigma^{\text{tot}}_{W^{-}\rightarrow \ell \nu}& = & 1399 \pm 21\ (\text{stat.})\ \pm…
| Process | Generator | Generator QCD precision |
| Signal Samples | ||
| Powheg-Box +Pythia 8 | NLO | |
| Powheg-Box +Pythia 8 | NLO | |
| Background Samples | ||
| Powheg-Box +Pythia 8 | NLO | |
| Powheg-Box +Pythia 8 | NLO | |
| Powheg-Box +Pythia 6 | NLO | |
| Herwig | LO | |
| Herwig | LO | |
| Herwig | LO | |
| Pythia 8 | LO | |
| Pythia 8 | LO | |
| -boson fiducial region | -boson fiducial region |
|---|---|
| GeV | GeV |
| | | |
| GeV | GeV |
| GeV |
| Lepton trigger | 0.14 | 0.13 | 1.07 | 1.07 | 0.03 | |
|---|---|---|---|---|---|---|
| Lepton reconstr. and ident. | 2.31 | 2.33 | 4.55 | 0.30 | 0.32 | 0.62 |
| Lepton isolation | 0.71 | 0.71 | 1.41 | 0.51 | 0.51 | 1.01 |
| Lepton scale and resolution | 0.44 | 0.43 | 0.34 | 0.05 | 0.05 | 0.04 |
| Recoil scale and resolution | 0.25 | 0.20 | – | 0.22 | 0.22 | – |
| 0.22 | 0.29 | 0.11 | 0.11 | 0.20 | 0.06 | |
| MC statistical uncertainty | 0.24 | 0.31 | 0.30 | 0.24 | 0.34 | 0.43 |
| Total | 2.5 | 2.5 | 4.8 | 1.3 | 1.3 | 1.3 |
| Source | Muon channel | Electron channel | ||||
| Muon trigger | A | A | A | – | – | – |
| Muon reconstruction/ID | A | A | A | – | – | – |
| Muon energy scale/resolution | A | A | A | – | – | – |
| Muon isolation | A | A | A | – | – | – |
| Electron trigger | – | – | – | A∗ | A∗ | A∗ |
| Electron reconstruction/ID | – | – | – | A | A | A |
| Electron energy scale/resolution | – | – | – | A | A | A |
| Electron isolation | – | – | – | A | A | A |
| Recoil related | – | A | A | – | A | A |
| EW background | A | B | B | A | B | B |
| Top-quark background | A | A | A | A | A | A |
| Multijet background | – | A | A | – | A | A |
| A | A | A | A | A | A | |
| Measurement | Observed | Background | Background | Background-subtracted |
|---|---|---|---|---|
| Channel | candidates | (EW + top) | (Multijet) | data |
| 3914 | ||||
| 2209 | ||||
| 4365 | ||||
| 2460 | ||||
| 430 | – | |||
| 646 | – |
| Value stat. syst. lumi. ( extr.) | Value stat. syst. lumi. ( extr.) | |
| [pb] | ||
| [pb] | ||
| [pb] | ||
| [pb] | ||
| [pb] | ||
| [pb] |
| Value stat. syst. lumi. ( extr.) | Value stat. syst. lumi. ( extr.) | |
|---|---|---|
| [pb] | ||
| [pb] | ) | ) |
| [pb] | ||
| [pb] | ) | |
| [pb] | ||
| [pb] | ) | |
| Quantity | Predicted cross-section [pb] |
|---|---|
| Predicted cross-section PDF uncertainty [pb] | |||||
|---|---|---|---|---|---|
| Quantity | MMHT14 | NNPDF31 | ATLASepWZ12 | HERAPDF20 | ABMP16 |
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
\AtlasTitle
Measurement of -boson and -boson production cross-sections in collisions at TeV with the ATLAS detector \AtlasAbstract The production cross-sections for and bosons are measured using ATLAS data corresponding to an integrated luminosity of 4.0 pb*-1* collected at a centre-of-mass energy TeV. The decay channels and are used, where can be an electron or a muon. The cross-sections are presented for a fiducial region defined by the detector acceptance and are also extrapolated to the full phase space for the total inclusive production cross-section. The combined (average) total inclusive cross-sections for the electron and muon channels are:
[TABLE]
Measured ratios and asymmetries constructed using these cross-sections are also presented. These observables benefit from full or partial cancellation of many systematic uncertainties that are correlated between the different measurements.
\AtlasRefCodeSTDM-2018-06 \PreprintIdNumberCERN-EP-2019-095 \LEcontactSteve Lloyd [email protected] \HepDataRecordhttps://www.hepdata.net/record/ins1742785 \AtlasJournalEur. Phys. J. C \AtlasJournalRefEur. Phys. J. C 79 (2019) 901 \AtlasDOI10.1140/epjc/s10052-019-7399-7 \AtlasCoverSupportingNoteSupport Notehttps://cds.cern.ch/record/2296490 \AtlasCoverCommentsDeadline1 July 2019 \AtlasCoverAnalysisTeam F. Balli, M. Boonekamp, C. Debenedetti, J. Ferrando, K. Gasnikova, A. Glazov, U. Klein, K. Lohwasser, E Lobodzinska, E. Rizvi, K. Schmieden,F. Sforza, G. Siragusa, A. Trofymov, M. Vincter, N. Zakharchuk \AtlasCoverEdBoardMemberP. Krieger (chair), D. Fassouliotis, C. Sawyer \AtlasCoverEgroupEditorsatlas-stdm-2018-06-contact-editors@cern.ch \AtlasCoverEgroupEdBoardatlas-stdm-2018-06-editorial-board@cern.ch
\size@chapter\sectfont
Contents
@afterheading@starttoc
toc
1 Introduction
The processes that produce and bosons111In this paper it is implicit that boson refers to bosons. in collisions via Drell–Yan annihilation are two of the simplest at hadron colliders to describe theoretically. At lowest order in quantum chromodynamics (QCD), -boson production proceeds via and -boson production via . Therefore, precision measurements of these production cross-sections yield important information about the parton distribution functions (PDFs) for quarks inside the proton. Factorisation theory allows PDFs to be treated separately from the perturbative QCD high-scale collision calculation as functions of the event energy scale, , and the momentum fraction of the parton, , for each parton flavour. Usually PDFs are defined for a particular starting scale and can be evolved to other scales via the DGLAP equations [1, 2, 3, 4]. Measurements of on-shell -boson production probe the PDFs in a range of that lies close to . The range of that is probed depends on the centre-of-mass energy, , of the protons and the rapidity coverage of the detector. Each measurement of these production cross-sections at a new value of thus provides information complementary to previous measurements. The combinations of initial partons participating in the production processes of ,, and bosons are different, so each process provides complementary information about the products of different quark PDFs.
This paper presents the first measurements of the production cross-sections for , and bosons in collisions at TeV. The data were collected by the ATLAS detector at the Large Hadron Collider (LHC) [5] in 2013 and correspond to an integrated luminosity of 4.0 pb*-1*. To provide further sensitivity to PDFs, and to reduce the systematic uncertainty in the predictions, ratios of these cross-sections and the charge asymmetry for -boson production are also presented. The measurements are performed for leptonic (electron or muon) decays of the and bosons, in a defined fiducial region, and also extrapolated to the total cross-section.
Previous measurements of the -boson and -boson production cross-sections in collisons at the LHC were performed by the ATLAS, CMS and LHCb Collaborations at TeV [6], 7 TeV [7, 8, 9, 10, 11], 8 TeV [12, 13, 14, 15, 16] and 13 TeV [17, 18, 19], and by the PHENIX and STAR Collaborations at the RHIC at GeV [20, 21] and 510 GeV [22]. This is the first measurement at 2.76 TeV. Other measurements of these processes were performed in collisons at TeV and 1.96 TeV by the CDF [23, 24, 25, 26, 27] and D0 [28] Collaborations, and at GeV and GeV by the UA1 [29] and UA2 [30] Collaborations.
2 ATLAS detector
The ATLAS detector [31] at the LHC covers nearly the entire solid angle around the collision point. It consists of an inner tracking detector surrounded by a thin superconducting solenoid, electromagnetic (EM) and hadronic calorimeters, and a muon spectrometer (MS) incorporating three large superconducting toroid magnets. The inner-detector system (ID) is immersed in a axial magnetic field and provides charged-particle tracking in the pseudorapidity range .222ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector and the -axis along the beam pipe. The -axis points from the IP to the centre of the LHC ring, and the -axis points upwards. Cylindrical coordinates are used in the transverse plane, being the azimuthal angle around the -axis. The pseudorapidity is defined in terms of the polar angle as . Angular distance is measured in units of .
The high-granularity silicon pixel detector covers the vertex region and typically provides three measurements per track. It is followed by the silicon microstrip tracker, which usually provides eight measurements from eight strip layers. These silicon detectors are complemented by the transition radiation tracker (TRT), which enables radially extended track reconstruction up to . The TRT also provides electron identification information based on the fraction of hits (typically 30 in total) above a higher energy-deposit threshold associated with the presence of transition radiation.
The calorimeter system covers the pseudorapidity range . Within the region , EM calorimetry is provided by barrel and endcap high-granularity lead/liquid-argon (LAr) sampling calorimeters, with an additional thin LAr presampler covering that is used to correct for energy loss in material upstream of the calorimeters. Hadronic calorimetry in this region is provided by the steel/scintillator-tile calorimeter, segmented into three barrel structures with , and two copper/LAr hadronic endcap calorimeters. The solid angle coverage is completed with forward copper/LAr and tungsten/LAr calorimeter modules optimised for EM and hadronic measurements, respectively.
The muon spectrometer comprises separate trigger and high-precision tracking chambers measuring the deflection of muons in a magnetic field generated by superconducting air-core toroids. The precision chamber system covers the region with three layers of monitored drift tubes, complemented by cathode strip chambers in the forward region, where the backgrounds are highest. The muon trigger system covers the range with resistive plate chambers in the barrel and thin gap chambers in the endcap regions.
The ATLAS detector selected events using a three-level trigger system [32]. The first-level trigger is implemented in hardware and used a subset of detector information to reduce the event rate to a design value of at most . This was followed by two software-based triggers that together reduced the event rate to about .
3 Data and simulation samples
The data used in this measurement were collected in February 2013 during a period when proton beams at the LHC were collided at a centre-of-mass energy of 2.76 TeV. During this running period a typical value of the instantaneous luminosity was cm*-2s-1*, significantly lower than in 7, 8 and 13 TeV data-taking conditions. The typical value of the mean number of collisions per proton bunch crossing (pile-up) was 0.3. Only data from stable collisions when the ATLAS detector was fully operational are used, yielding a data sample corresponding to an integrated luminosity of 4.0 pb*-1*.
Samples of Monte Carlo (MC) simulated events are used to estimate the signals from -boson and -boson production, and the backgrounds containing prompt leptons: electroweak-diboson production and top-quark pair () production. Background contributions arising from multijet events that do not contain prompt leptons are estimated directly from data, with simulated events used to cross-check these estimations in the muon channel.
Production of single and bosons was simulated using Powheg-Box v1 r1556 [33, 34, 35, 36]. The parton showering was performed using Pythia 8.17 [37]. The PDF set used for the simulation was CT10 [38], and the parton shower parameter values were those of the AU2 tune [39]. Additional quantum electrodynamics (QED) emissions from electroweak (EW) vertices and charged leptons were simulated using Photos++ v3.52 [40]. Additional samples of simulated -boson events generated with Sherpa 2.1 [41] are used to estimate uncertainties arising from the choice of event generator model. In these Sherpa samples, simulation of -boson production in association with up to two additional partons was performed at next-to-leading order (NLO) in QCD while production of bosons in association with three or four additional partons was performed at leading order (LO) in QCD. The sample cross-sections were normalised to next-to-next-to-leading-order (NNLO) QCD predictions for the total cross-sections described in Section 8.
Powheg-Box v1 r2330 was used to generate samples [42]. These samples had parton showering performed using Pythia 6.428 [43] with parameters corresponding to the Perugia2011C tune [44]. The CT10 PDF set was used. Additional QED final-state radiative corrections were applied using Photos++ v3.52 and -lepton decays were performed using Tauola v25feb06 [45]. Single production of top quarks is a negligible contribution to this analysis, compared with production, so no such samples were generated.
Production of two massive electroweak bosons () was simulated using Herwig 6.5 [46], with multiparton interactions modelled using Jimmy 4.13 [47]. The CTEQ6L1 PDF set [48] and AUET2 tune [49] were used for these samples.
Multijet production containing heavy-flavour final states, arising from the production of or pairs, were simulated using Pythia 8.186. The CTEQ6L1 PDF set and AU2 tune were used. Events were required to contain an electron or muon with transverse momentum GeV and .
The detector response to generated events was simulated by passing the events through a model of the ATLAS detector [50] based on Geant4 [51]. Additional minimum-bias events generated using Pythia 8.17 and the A2 set of tuned parameters, were overlaid in such a way that the distribution of for simulated events reproduced that in the real data. The resulting events were then passed through the same reconstruction software as the real data.
The simulated samples used for the baseline analysis are summarised in Table 1, which shows the generator used for each process together with the order in QCD at which they were generated.
4 Event selection
This section describes the selection of events consistent with the production of bosons or bosons. The -boson selection requires events to contain a single charged lepton and large missing transverse momentum. The -boson selection requires events to contain two charged leptons with opposite charge and the same flavour.
Events were selected by triggers that required at least one charged electron (muon) with p_{\text{T}}>15\leavevmode\nobreak\GeV (10 GeV). These thresholds yield an event sample with a uniform efficiency as a function of the and requirements used subsequently to select the final event sample. The hard-scatter vertex, defined as the vertex with highest sum of squared track transverse momenta (for tracks with MeV), is required to have at least three associated tracks.
Electrons are reconstructed from clusters of energy in the EM calorimeter that are matched to a track reconstructed in the ID. The electron is required to have GeV and (excluding the transition region between barrel and endcap calorimeters of ). Each electron must satisfy a set of identification criteria designed to suppress misidentified photons or jets. Electrons are required to satisfy the medium selection, following the definition provided in Ref. [52]. This includes requirements on the shower shape in the EM calorimeter, the leakage of the shower into the hadronic calorimeter, the number of hits measured along the track in the ID, and the quality of the cluster–track matching. A Gaussian sum filter [53] algorithm is used to re-fit the tracks and improve the estimated electron track parameters. To suppress background from misidentified objects such as jets, the electron is required to be isolated using calorimeter-based criteria. The sum of the transverse energies of clusters lying within a cone of size around the centroid of the electron cluster and excluding the core333The core of the shower is the contribution within around the cluster barycentre. must be less than 10% of the electron .
Muon candidates are reconstructed by combining tracks reconstructed in the ID with tracks reconstructed in the MS [54]. They are required to have GeV and . The muon candidates are also required to be isolated, by requiring that the scalar sum of the of additional tracks within a cone of size around the muon is less than 80% of the muon .
The missing transverse momentum vector [55] () is calculated as the negative vector sum of the transverse momenta of electrons and muons, and of the transverse momentum of the recoil. The magnitude of this vector is denoted by . The recoil vector is obtained by summing the transverse momenta of all clusters of energy measured in the calorimeter, excluding those within of the lepton candidate. The momentum vector of each cluster is determined by the magnitude and coordinates of the energy deposits; the cluster is assumed to be massless. Cluster energies are initially measured assuming that the energy deposition occurs only through EM interactions, and are then corrected for the different calorimeter responses to hadrons and electromagnetically interacting particles, for losses due to dead material, and for energy that is not captured by the clustering process [56]. The definition of the recoil does not make use of reconstructed jets, to avoid threshold effects. The procedure used to calibrate the recoil closely follows that used in the recent ATLAS measurement of the -boson mass [57], first correcting the modelling of the overall recoil in simulation and then applying corrections for residual differences in the recoil response and resolution that are derived from -boson data and transferred to the -boson sample.
The -boson selection requires events to contain exactly one lepton (electron or muon) candidate and have GeV. The lepton must match a lepton candidate that met the trigger criteria. The transverse mass, , of the -boson candidate in the event is calculated using the lepton candidate and according to . The transverse mass in -boson production events is expected to exhibit a Jacobian peak around the -boson mass. Thus, requiring that suppresses background processes. After these requirements there are 3914 events in the channel, 2209 events in the channel, 4365 events in the channel, and 2460 events in the channel.
The -boson selection requires events to contain exactly two lepton candidates with the same flavour and opposite charge. At least one lepton must match a lepton candidate that met the trigger criteria. Background processes are suppressed by requiring that the invariant mass of the lepton pair satisfies GeV. After these requirements there are 430 events in the channel, and 646 events in the channel.
5 Background estimation
The background processes that contribute to the sample of events passing the -boson and -boson selections can be separated into two categories: those estimated from MC simulation and theoretical calculations, and those estimated directly from data. The main backgrounds that contribute to the event sample passing the -boson selection are processes with a -lepton decaying into an electron or muon plus neutrinos, leptonic -boson decays where only one lepton is reconstructed, and multijet processes. The main background contribution to the event sample passing the -boson selection is production of two massive electroweak bosons.
The backgrounds arising from , , diboson production, and production are estimated from the simulated samples described in Section 3. Predictions of the backgrounds to the -boson and -boson production measurements arising from multijet production suffer from large theoretical uncertainties, and therefore the contribution to this background in the -boson measurement is estimated from data. This is achieved by constructing a shape template for the background using a discriminating variable in a control region and then performing a template fit to the same distribution in the signal region to extract the background contribution. The choice of template variable is motivated by the difference between signal and background and by the available number of events. Previous ATLAS measurements at 7 TeV [7] and 13 TeV [18] found that multijet production makes a background contribution of less than 0.1% for -boson measurements; this is therefore neglected.
Electron candidates in multijet background events are typically misidentified candidates produced when jets mimic the signature of an electron, for example when a neutral pion and a charged pion overlap in the detector. Additional candidates can arise from ’non-prompt’ electrons produced when a photon converts, and in decays of heavy-flavour hadrons. To construct a control region for the multijet template, a selection is used that differs from the -boson selection described in Section 4 in only two respects: the medium electron identification criteria are inverted (while keeping the looser identification criteria) and the requirement is removed. By construction, this control region is statistically independent of the -boson signal region. A template for the shape of the multijet background in the distribution is then obtained from that distribution in the control region after subtraction of expected contributions from the signal and other backgrounds determined using MC samples. The normalisation of the multijet background template in the signal region is extracted by performing a fit of the distribution (applying all signal criteria except the requirement on ) to a sum of the templates for the multijet background, the signal, and all other backgrounds. The normalisation of the signal is allowed to vary freely in the fit as is the multijet background; however, the other backgrounds are only allowed to vary from their expected values by up to 5%, corresponding to the largest level of variation in predicted electroweak-boson production cross-sections obtained from varying the choice of PDF. The normalisation from this fit can then be used together with the inverted selection to construct multijet background distributions in any other variable that is not correlated with the electron identification criteria.
Muon candidates in multijet background events are typically ‘non-prompt’ muons produced in the decays of hadrons. The multijet background contribution to the selection is estimated by using the same method as described for the selection. In this case the control region is defined by inverting the isolation requirement and removing the requirement on . The distribution used for the fits is .
The overall number of multijet background events is estimated from a fit to the total -boson sample. Comparisons between the fitted distributions and data for and are shown in Figure 1. Fits to the separate -boson and -boson samples are used in the evaluation of the systematic uncertainties, as described in Section 7. The final estimated multijet contributions are events for and and events for and . The relative contribution of the multijet events (1%) is lower than in 13 TeV (4%) and 7 TeV (3%) data. This is in agreement with expectations for this lower pile-up running, where the resolution in is improved compared to the higher pile-up running.
6 Correction for detector effects
The measurements in this paper are performed within specific fiducial regions and extrapolated to the total -boson or -boson phase space. The fiducial regions are defined by the kinematic and geometric selection criteria given in Table 2; in simulations these are applied at the generator level before the emission of QED final-state radiation from the decay lepton(s) (QED Born level).
The fiducial -boson/-boson production cross-section is obtained from the number of observed events meeting the selection criteria after background contributions are subtracted, , using the following formula:
[TABLE]
where is the total integrated luminosity of the data samples used for the analysis. The factor is the ratio of the number of generated events that satisfy the final selection criteria after event reconstruction to the number of generated events within the fiducial region. It includes the efficiency for triggering, reconstruction and identification of events falling within the acceptance. The different components of the efficiency are calculated using a mixture of MC simulation and measurements from data.
The total -boson and boson production cross-sections are obtained using the following formula:
[TABLE]
The factor is the per-lepton branching fraction of the vector boson. The factor is the acceptance for -boson events being studied. It is defined as the fraction of generated events that satisfy the fiducial requirements. This acceptance is determined using MC signal samples, corrected to the generator QED Born level, and is used to extrapolate the measured cross-section in the fiducial region to the full phase space. The central values of are around 0.6 for these measurements, compared with 0.5 at TeV and 0.4 at TeV, so the fiducial region is closer to the full phase space in this measurement than for those at higher centre-of-mass energies. This is due to a combination of higher thresholds for leptons in other measurements, and more-central production of vector bosons at lower . The values of are approximately 0.67 for the channels and 0.75 for the channels. The values of are 0.55 for the channel and 0.79 for . The values are a little higher than for previous measurements at TeV and TeV.
7 Systematic uncertainties
The systematic uncertainty in the electron reconstruction and identification efficiency is estimated using the tag-and-probe method in 8 TeV data [52, 58] and extrapolated to the 2.76 TeV dataset. The extrapolation procedure results in absolute increases of 2%, due to uncertainties in the effect of the differing pile-up conditions in 2.76 TeV data relative to the 8 TeV data, as well as a different setting of the noise filtering in the LAr calorimeter of the 2.76 TeV data relative to the 8 TeV data. These uncertainties were estimated using a comparison between 7 TeV and 8 TeV data and MC samples, after having established that the central values of the efficiencies are the same for different centre-of-mass energies when the same LAr filter settings are used. A similar methodology had been used for internal estimates of the electron efficiency performance at 13 TeV before the start of Run-2 data taking and was found to give a good prediction of the efficiencies in data as well as a conservative estimate of the uncertainties. Transverse-momentum-dependent isolation corrections, calculated with the tag-and-probe method in 2.76 TeV data, are very close to 1, so the systematic uncertainty in the electron isolation requirement is set to the size of the correction itself, that is % for low and % for higher . The electron energy scale has associated statistical uncertainties and systematic uncertainties arising from a possible bias in the calibration method, the choice of generator, the presampler energy scale, and imperfect knowledge of the material in front of the EM calorimeter [59]. The total energy-scale uncertainty is calculated as the sum in quadrature of these components.
Systematic uncertainties associated with the muon momentum can be divided into three major independent categories: momentum resolution of the MS track, momentum resolution of the ID track, and an overall scale uncertainty [54]. The total momentum scale/resolution uncertainty is the sum in quadrature of these components. An -independent uncertainty of approximately % in the muon trigger efficiency, determined using the tag-and-probe method [54] in 2.76 TeV data, is taken into account. Furthermore, a - and - dependent uncertainty in the identification and reconstruction efficiencies of approximately %, derived using the tag-and-probe method on 8 TeV data is applied. The uncertainty in the -dependent isolation correction in the muon channel, calculated with the tag-and-probe method in 2.76 TeV data, is about % for low and % for higher .
The luminosity uncertainty for the 2.76 TeV data is %. This is determined, following the same methodology as was used for the 7 TeV data recorded in 2011 [60], from a calibration of the luminosity scale derived from beam-separation scans performed during the 2.76 TeV operation of the LHC in 2013.
Systematic uncertainties in the arising from the smearing and bias corrections applied to obtain satisfactory modelling of the recoil [55] affect the factors in the measurement, and are taken into account.
Uncertainties arising from the choice of PDF set are evaluated using the error sets of the initial CT10 PDF set (at 90% confidence level (CL)) and from comparison with the results obtained using the central PDF sets from ABKM09 [61], NNPDF23 [62], and ATLAS-epWZ12 [63]. The effect of this uncertainty on () is estimated to be % (1.2%), and the effect on is estimated to be 1.4%. The effect on is between % and % depending on the channel.
A summary of the systematic uncertainties in the factors is shown in Table 3. The muon trigger, and electron reconstruction and identification uncertainties are dominant.
Uncertainties arising from the choice of event generator and parton shower models are estimated by comparing results obtained when using Sherpa 2.1 signal samples instead of the (nominal) Powheg-Box +Pythia 8. The effect of this uncertainty on is estimated to be 0.9%.
The systematic uncertainty in the multijet background estimation can be divided into several components: the normalisation uncertainty from the fit, the uncertainty in the modelling of electroweak processses by simulated samples in the fitted region, uncertainty from fit bias due to binning choice, and uncertainty from template shape. The scale normalisation uncertainty from the fit is approximately % for the channel. This uncertainty is neglected in the channel where the template bias is dominant. The mismodelling uncertainty is estimated by comparison of the fit results for and , and for the combined candidates. The central value used is with the uncertainties and , where is the fitted number of background events, is the fitted number of and is the fitted total number of background events. In the channel this leads to an uncertainty of % in the multijet background. In the channel the multijet template normalisation is derived from the fit in the small- region, where electroweak contributions are negligible and there are many data events, and this source of systematic error is found to be negligible. The fit-bias uncertainty arising from the choice of bin width is estimated by repeating the fit with different binnings. This component is negligible in the case and % in the case. The uncertainty due to a potential bias from template choice is estimated by employing different template selections. For the channel, different inverted-isolation criteria were investigated. The overall differences are considered negligible. For the channel, template variations were estimated from fits that use + MC samples as the multijet templates, leading to an uncertainty of ; this is the largest uncertainty in the multijet background in the channel.
Combining results and building ratios or asymmetries of results require a model for the correlations of particular systematic uncertainties between different measurements. Correlations arise mostly due to the fact that electrons, muons, and the recoil are reconstructed identically in the different measurements. Further correlations occur due to similarities in the analysis methodology such as the methods of signal and background estimation.
The systematic uncertainties from the electroweak background estimations are treated as uncorrelated between the -boson and -boson measurements, and fully correlated among different flavour decay channels of the and boson. The top-quark background is treated as fully correlated across all -boson and -boson decay channels. The multijet background and recoil-related systematic uncertainties are also treated as fully correlated between all four -boson decay channels despite there being an expected uncorrelated component, since the statistical uncertainty is dominant in this case.
The systematic uncertainties due to the choice of PDF are treated as fully correlated between all -boson and -boson channels. The uncertainties in electron and muon selection, reconstruction and efficiency are treated as fully correlated between all -boson and -boson channels.
A simplified form of the correlation model with the grouped list of the sources of systematic errors is presented in Table 4.
8 Results
The numbers of events passing the event selections described in Section 4 are presented in Table 5, together with the estimated background contributions described in Section 5. The distribution of for candidate events is shown in Figure 2, compared with the expected distribution for signal plus backgrounds, where the signal is normalised to the NNLO QCD prediction. Similarly, Figure 3 shows the distribution of for candidate events compared with the expectations for signal. In this case, the background contributions are not shown, because they would not be visible in the figure if included.
The measured fiducial () and total () cross-sections in the electron and muon channels are presented separately in Table 6. For these measurements, the dominant contribution to the systematic uncertainty arises from the luminosity determination.
The results obtained from the electron and muon final states are consistent. The fiducial measurements from electron and muon final states are combined following the procedure described in Ref. [64] and the result is extrapolated to the full phase space to obtain the total cross-section. The total -boson cross-section is calculated by summing the separate and cross-sections. The results are shown in Table 7.
Theoretical predictions of the fiducial and total cross-sections are computed for comparison with the measured cross-sections using Dynnlo 1.5 [65] which provides calculations at NNLO in the strong-coupling constant, , including the boson decays into leptons ( or ) with full spin correlations, finite width and interference effects. These calculations allow kinematic requirements to be implemented for direct comparison with experimental data. The procedure used follows that used for the previous ATLAS measurement at TeV [7].
Corrections for NLO EW effects are calculated with Fewz 3.1 [66, 67, 68, 69], for the bosons and with Sanc [70, 71] for the bosons. The calculation was done in the EW scheme [72]. The following input parameters are taken from the Particle Data Group’s Review of Particle Properties 2014 edition [73]: the Fermi constant, the masses and widths of and bosons as well as the elements of the CKM matrix. The cross-sections for vector bosons decaying into these leptonic final states are calculated such that they match the definition of the measured cross-sections in the data. Thus, from complete NLO EW corrections, the following components are included: virtual QED and weak corrections, real initial-state radiation (ISR), and interference between ISR and real final-state radiation (FSR) [74]. The calculated effect of these corrections on the cross-sections is % for , % for , and % for . Dynnlo is used for the central values of the predictions while Fewz is used for the PDF variations and all other systematic variations such as QCD scale and . The predictions are calculated using the CT14nnlo [75], NNPDF3.1 [76], MMHT14nnlo68cl [77], ABMP16 [78], HERAPDF2.0 [79], and ATLAS-epWZ12nnlo PDF sets. The dynamic scale, , and fixed scale, , are used as the nominal renormalisation, , and factorisation, , scales for and predictions, respectively.
Theoretical uncertainties in the predictions are also derived from the following sources:
PDF: these uncertainties are evaluated from the variations of the NNLO PDFs according to the recommended procedure for each PDF set. A table with all PDF uncertainties and their central values is shown in Appendix A; the PDF uncertainty from CT14nnlo was rescaled from 90% CL to 68% CL.
Scales: the scale uncertainties are defined by the envelope of the variations in which the scales are changed by factors of two subject to the constraint .
: the uncertainty due to was estimated by varying the value of used in the CT14nnlo PDF set by , corresponding to a 68% CL variation.
The statistical uncertainties in these theoretical predictions are negligible.
The numerical values of the predictions for the CT14nnlo PDF set are presented in Table 8. The predictions for the acceptance factor can differ by a few percent from those derived from simulated signal samples, this may be due to a poorer description of production of low -bosons by the fixed-order calculations. The predictions are shown in comparison with the combined -boson and -boson production measurements, and with results from and collisions at other centre-of-mass energies in Figure 4. A comparison of the measurements with predictions from various different PDF sets is presented in Figures 5 and 6. Overall there is good agreement.
Taking ratios of measurements leads to results that have significantly reduced systematic uncertainties due to full or partial cancellation of correlated systematic uncertainties, as discussed in Section 7. The ratios of the fiducial cross-sections for -boson and -boson production are presented, together with the ratio for -boson and -boson production, in Figure 7. It can be seen that the predictions from the different PDF sets are mostly in good agreement with the measurements. There is a slight (less than two standard deviations) tension between the data and the prediction using the ABMP16 PDF set. The measured values of the ratios are:
;
.
The measurement of the ratio is sensitive to the and valence quark distributions, while the ratio can place constraints on the strange quark distributions. A common alternative way of presenting this information is in terms of the charge asymmetry, , in -boson production:
[TABLE]
This observable also benefits from the cancellation of systematic uncertainties in the same way as the cross-section ratios. The measured value is:
.
The ratio of measured cross-sections in the electron and muon decay channels provides a test of lepton universality in -boson decays. The measured ratios are:
[TABLE]
[TABLE]
[TABLE]
[TABLE]
These results lie within one standard deviation of the Standard Model prediction and previous measurements by ATLAS.
9 Conclusion
This paper presents measurements of the and production cross-sections based on about 12 400 -boson and 1100 -boson candidates, after subtracting background events, reconstructed from 2.76 TeV proton–proton collision data recorded by the ATLAS detector at the LHC, corresponding to integrated luminosity of 4.0 pb*-1*. The total inclusive -boson production cross-sections for the combined electron and muon channels are
,
,
and the total inclusive -boson cross-section in the combined electron and muon channels is:
The results obtained, and the ratios and charge asymmetries constructed from them, are in agreement with theoretical calculations based on NNLO QCD.
Appendix
Appendix A Theoretical predictions
This appendix presents the theoretical predictions used for comparison with the measurements in the main body of the paper. Table 9 shows the predictions using the MMHT14nnlo68cl, NNPDF31_nnlo_as_0118, ATLASepWZ12, HERAPDF2.0, and ABMP16 PDF sets with associated PDF uncertainties.
Acknowledgements
We thank CERN for the very successful operation of the LHC, as well as the support staff from our institutions without whom ATLAS could not be operated efficiently.
We acknowledge the support of ANPCyT, Argentina; YerPhI, Armenia; ARC, Australia; BMWFW and FWF, Austria; ANAS, Azerbaijan; SSTC, Belarus; CNPq and FAPESP, Brazil; NSERC, NRC and CFI, Canada; CERN; CONICYT, Chile; CAS, MOST and NSFC, China; COLCIENCIAS, Colombia; MSMT CR, MPO CR and VSC CR, Czech Republic; DNRF and DNSRC, Denmark; IN2P3-CNRS, CEA-DRF/IRFU, France; SRNSFG, Georgia; BMBF, HGF, and MPG, Germany; GSRT, Greece; RGC, Hong Kong SAR, China; ISF and Benoziyo Center, Israel; INFN, Italy; MEXT and JSPS, Japan; CNRST, Morocco; NWO, Netherlands; RCN, Norway; MNiSW and NCN, Poland; FCT, Portugal; MNE/IFA, Romania; MES of Russia and NRC KI, Russian Federation; JINR; MESTD, Serbia; MSSR, Slovakia; ARRS and MIZŠ, Slovenia; DST/NRF, South Africa; MINECO, Spain; SRC and Wallenberg Foundation, Sweden; SERI, SNSF and Cantons of Bern and Geneva, Switzerland; MOST, Taiwan; TAEK, Turkey; STFC, United Kingdom; DOE and NSF, United States of America. In addition, individual groups and members have received support from BCKDF, CANARIE, CRC and Compute Canada, Canada; COST, ERC, ERDF, Horizon 2020, and Marie Skłodowska-Curie Actions, European Union; Investissements d’ Avenir Labex and Idex, ANR, France; DFG and AvH Foundation, Germany; Herakleitos, Thales and Aristeia programmes co-financed by EU-ESF and the Greek NSRF, Greece; BSF-NSF and GIF, Israel; CERCA Programme Generalitat de Catalunya, Spain; The Royal Society and Leverhulme Trust, United Kingdom.
The crucial computing support from all WLCG partners is acknowledged gratefully, in particular from CERN, the ATLAS Tier-1 facilities at TRIUMF (Canada), NDGF (Denmark, Norway, Sweden), CC-IN2P3 (France), KIT/GridKA (Germany), INFN-CNAF (Italy), NL-T1 (Netherlands), PIC (Spain), ASGC (Taiwan), RAL (UK) and BNL (USA), the Tier-2 facilities worldwide and large non-WLCG resource providers. Major contributors of computing resources are listed in Ref. [80].
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] V.. Gribov and L.. Lipatov “Deep inelastic e p scattering in perturbation theory” [Yad. Fiz. 15, (1972) 781] In Sov. J. Nucl. Phys. 15 , 1972, pp. 438–450
- 2[2] L.. Lipatov “The parton model and perturbation theory” [Yad. Fiz. 20, (1974) 181] In Sov. J. Nucl. Phys. 20 , 1975, pp. 94–102
- 3[3] Guido Altarelli and G. Parisi “Asymptotic freedom in parton language” In Nucl. Phys. B 126 , 1977, pp. 298–318 DOI: 10.1016/0550-3213(77)90384-4 · doi ↗
- 4[4] Yuri L. Dokshitzer “Calculation of the structure functions for deep inelastic scattering and e+ e- nnihilation by perturbation theory in quantum chromodynamics.” [Zh. Eksp. Teor. Fiz. 73, (1972) 1216] In Sov. Phys. JETP 46 , 1977, pp. 641–653
- 5[5] Lyndon Evans and Philip Bryant “LHC Machine” In JINST 3 , 2008, pp. S 08001 DOI: 10.1088/1748-0221/3/08/S 08001 · doi ↗
- 6[6] ATLAS Collaboration “Measurements of W 𝑊 W and Z 𝑍 Z boson production in p p 𝑝 𝑝 pp collisions at s = 5.02 𝑠 5.02 \sqrt{s}=5.02 Te V with the ATLAS detector” In Eur. Phys. J. C 79.2 , 2019, pp. 128 DOI: 10.1140/epjc/s 10052-019-6622-x · doi ↗
- 7[7] ATLAS Collaboration “Precision measurement and interpretation of inclusive W + superscript 𝑊 W^{+} , W − superscript 𝑊 W^{-} and Z / γ ∗ 𝑍 superscript 𝛾 Z/\gamma^{*} production cross sections with the ATLAS detector” In Eur. Phys. J. C 77.6 , 2017, pp. 367 DOI: 10.1140/epjc/s 10052-017-4911-9 · doi ↗
- 8[8] CMS Collaboration “Measurement of inclusive W 𝑊 W and Z 𝑍 Z production cross sections in p p 𝑝 𝑝 pp collisions at s = 7 𝑠 7 \sqrt{s}=7 Te V” In JHEP 10 , 2011, pp. 132 DOI: 10.1007/JHEP 10(2011)132 · doi ↗
