Study of W boson production in pPb collisions at sqrt(s[NN]) = 5.02 TeV

The first study of W boson production in pPb collisions is presented, for bosons decaying to a muon or electron, and a neutrino. The measurements are based on a data sample corresponding to an integrated luminosity of 34.6 inverse nanobarns at a nucleon-nucleon centre-of-mass energy of sqrt(s[NN]) = 5.02 TeV, collected by the CMS experiment. The W boson differential cross sections, lepton charge asymmetry, and forward-backward asymmetries are measured for leptons of transverse momentum exceeding 25 GeV, and as a function of the lepton pseudorapidity in the abs(eta[lab])<2.4 range. Deviations from the expectations based on currently available parton distribution functions are observed, showing the need for including W boson data in nuclear parton distribution global fits.


Introduction
Electroweak boson production in proton-nucleus and nucleus-nucleus collisions at the CERN LHC offers a unique opportunity to probe nuclear parton distribution functions (nPDFs) [1][2][3][4].Leptonic decays of electroweak bosons are of particular interest since leptons do not interact strongly with the medium produced in these collisions [5,6].As compared to those in a proton, the nPDFs are expected to be depleted (shadowing) for partons carrying small momentum fractions x 10 −2 , and enhanced (anti-shadowing) in the 5 × 10 −2 x 10 −1 range [7].However, because of the lack of available data, parton densities are less precisely known for nuclei than for nucleons.As a consequence, precise calculations describing hard processes in high-energy heavy ion collisions are limited by uncertainties in the nPDFs.For W boson production, the dominant processes at LHC energies are ud → W + and du → W − , principally reflecting interactions that take place between valence quarks and sea antiquarks.According to Ref. [4], PDF nuclear modifications could affect the yield of W bosons in pPb collisions at the LHC by as much as 15% in certain kinematic regions.Therefore, precise measurements of W boson production in heavy ion collisions might lead to an improved determination of the nPDFs.Moreover, asymmetries in the individual yields of W + and W − should permit the flavour decomposition of u and d quark distributions in nuclei.
The ATLAS [8,9] and CMS [10,11] collaborations reported the observations of Z bosons in heavy ion interactions, at a centre-of-mass energy of 2.76 TeV per nucleon pair.These data showed that the Z boson yields per nucleon-nucleon (NN) collision are essentially unmodified by the medium produced in the collisions.Although W bosons decaying to a lepton and a neutrino are more difficult to detect, their rate is about ten times larger than that of Z bosons decaying to leptonic final states.The production of W bosons in PbPb collisions was reported by CMS [12] and ATLAS [13], using data corresponding to an integrated luminosity of 7.3 µb −1 and 150 µb −1 , collected in 2010 and 2011, respectively.The W boson yield per NN collision was shown to be compatible with the one measured in pp collisions, when taking into account isospin effects arising from the mixture of protons and neutrons in the colliding nuclei.However, the presence of 10-20% nPDF effects on Z and W boson production could not be excluded due to the relatively large experimental and theoretical uncertainties of these results.
The 2013 pPb LHC run provides the best currently available data sample to look for initial-state effects (such as PDF modifications) using electroweak bosons.The NN-equivalent luminosity is of the same order of magnitude as for the 2011 PbPb run, and the production cross sections are approximately a factor of two greater owing to the increased energy, 5.02 TeV per nucleon pair.Furthermore, the asymmetry of the pPb collision system allows for the measurement of other observables such as forward-backward pseudorapidity asymmetries.This Letter reports a study of W boson production in a sample of pPb collisions corresponding to an integrated luminosity of (34.6 ± 1.2) nb −1 [14], collected by the CMS experiment.

Experimental methods
The direction of the proton beam was initially opposite to the positive direction of the CMS longitudinal axis [15], and was reversed after 60% of the data were taken.The beam energies were 4 TeV for protons and 1.58 TeV per nucleon for lead nuclei, resulting in a centre-of-mass energy per nucleon pair of √ s NN = 5.02 TeV.As a result of the energy difference of the colliding beams, the NN centre-of-mass frame in pPb collisions was not at rest with respect to the laboratory frame.Massless particles emitted at pseudorapidity η in the NN centre-of-mass frame are detected at η lab = η − 0.465 (first proton beam orientation) and η lab = η + 0.465 (second proton beam orientation) in the CMS coordinate system, as defined in Ref. [15].The results presented hereafter are expressed in the usual convention of the proton-going side defining the positive pseudorapidity.It coincides with the CMS convention in the second period of data taking, the first one being reversed before summing yields from the two beam configurations.
A detailed description of the CMS detector can be found elsewhere [15].Its central feature is a superconducting solenoid of 6 m internal diameter, providing a magnetic field of 3.8 T.
Within the field volume are the 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.Selection criteria similar to the ones developed in Ref. [16] are applied to the pPb sample to remove events with electromagnetic, beam-gas, or multiple collisions (pileup).The W boson yields are corrected for the induced (4.0 ± 0.5)% signal loss.
The primary signature of a W boson is a high transverse momentum (p T ) lepton.The current analysis is restricted to leptons of p T greater than 25 GeV/c.The muon analysis is based on a sample triggered by requiring a single muon with p T above 12 GeV/c, while the electron analysis uses an ECAL-triggered sample with a transverse energy threshold of 15 GeV.Leptons are reconstructed with the same algorithms as in proton-proton collisions [17,18], and standard selection criteria are applied, as in Refs.[12,19].A special electron charge determination, as described in Ref. [20], is used in order to reduce the electron charge misidentification to a sub-percent level.Events are reconstructed using particle-flow (PF) techniques [21,22], which reconstruct and classify individual particles with an optimised combination of all subdetector information.
Two criteria are used to remove specific background sources.First, events with two oppositely charged leptons, with the second lepton p T greater than 15 (10) GeV/c for muons (electrons) are removed, since they correspond to well-identified processes like Drell-Yan, Z boson or high-p T quarkonium production.Second, the leptons are required to be isolated, in order to reduce the contamination coming from jet fragmentation.The energies of all PF candidates are summed within a cone centred around the lepton, with the exception of the lepton itself.The lepton is considered isolated if the total transverse energy in the cone is small compared to its transverse momentum.For muons, a cone of radius ∆R = √ (∆η)2 + (∆φ) 2 = 0.3 is used, where ∆η and ∆φ are the pseudorapidity and azimuthal distances to the lepton.The candidate is rejected if the in-cone transverse energy is greater than 10% of the muon p T .For electrons, a cone of ∆R = 0.4 is used, and only particles with p T greater than 1 GeV/c are summed, to reduce the underlying-event enhanced contribution.The electron candidate is rejected if the resulting transverse energy is greater than 11.5% (9.5%) of the electron p T , for the ECAL barrel (endcaps).
An important characteristic of events containing a W → ν decay is the missing transverse energy (E T / ) associated with the undetected neutrino.It is computed as the magnitude of the vectorial sum of transverse momenta of all the PF candidates in the event.The analysis is performed using ten lepton pseudorapidity bins, each 0.5 wide except for the most forward and backward regions (2 < |η lab | < 2.4).After having applied the lepton selection criteria, examples of the resulting E T / distributions are shown in Fig. 1 for µ + and e + , in the most central (−0.5 < η lab < 0.0) and furthest forward (2.0 < η lab < 2.4) ranges.The distributions for other bins and for the negative leptons are similar.
events within the −0.5 < η lab < 0.0 (top) and 2.0 < η lab < 2.4 (bottom) ranges.Binned fits to the data (red points) are performed with four contributions, stacked from bottom to top: multijet (QCD, blue), W + → τ + ν (brown), Z → (white) and W + → + ν (yellow).The η lab regions are defined such that the proton is moving towards positive η lab values.Error bars represent statistical uncertainties.The lower panels display the data divided by the result of the fit, with the band representing the statistical uncertainties on the sum of the fit components, for each E T / bin.
To extract the number of events with a lepton coming from a W boson, binned fits of these distributions are performed, including the signal and main background contributions, in each η lab bin.The E T / shapes assumed for the electroweak processes, namely the W ± → ± ν signal as well as background from W ± → τ ± ν and Z → + − , are determined by the simulations described hereafter, taking into account the acceptance and efficiency.Their relative normalization is given by the unmodified theoretical cross sections (as computed in Ref. [23]).A maximal 20% variation of the W/Z normalization ratio is taken into account, due to potentially different nuclear modifications of the Z and W bosons, and resulting in a 1-3% systematic uncertainty in the extracted W yields.The noticeable difference between the E T / distributions for the Z → e + e − and Z → µ + µ − processes in the forward region (bottom plots of Fig. 1) results from the greater ECAL coverage allowing missed electrons with 2.4 < |η lab | < 3.0 to be accounted for in the E T / calculation.The shape of the QCD multijet background is modelled by the functional form ).It is shown to reproduce the E T / shape of data events containing non-isolated leptons, with the E T / 0 , α, and β parameters, which are observed to depend mildly and linearly on the cone/lepton transverse energy ratio.These fitted parameters are then extrapolated to the isolated lepton signal regime and the resulting function is used as the QCD background shape.The multijet background contribution is larger in the electron channel because the misidentified lepton rate is higher, particularly due to a contribution from photon-jet events.Contributions from other sources, such as tt production and high-p T quarkonia, were found to be negligible.
A small charge misidentification correction (less than 0.2%) is applied to the electron yields; this correction is negligible for muons.All fits are of good quality, as illustrated by the bottom panels of Fig. 1 that show the ratio of the data to the fit outcome.The observed numbers of leptons coming from W boson decays over the entire pseudorapidity range are: 11660 ± 111 µ + , 9459 ± 99 µ − , 9892 ± 116 e + , and 7872 ± 101 e − , where the uncertainty is statistical, determined by the fit procedure.
In order to correct for inefficiencies in the lepton trigger, reconstruction, and selection, the electroweak processes W → ν have been simulated using the PYTHIA 6.424 generator [24] with a mixture of pp and pn interactions corresponding to pPb collisions.The detector response to each PYTHIA signal event is simulated with GEANT4 [25] and then embedded in a minimum bias pPb background event.These background events are produced with the HIJING event generator [26] and passed through GEANT4 as well.Each simulation is done twice, once for each proton beam direction, and includes a boost to reproduce the 0.465 rapidity shift.The embedding is done at the level of detector hits, and the signal and background events share the same generated vertex location.The embedded event is then processed through the trigger emulation and the full event reconstruction chain.The resulting reconstructed events are then reweighted to match the distributions observed in data of the event vertex and activity (as measured in the HF calorimeters).The obtained efficiencies vary with η lab (with higher efficiencies at mid-rapidity), from 59% to 89% for muons, and from 51% to 84% for electrons.
The various components of the single-lepton efficiency are also directly computed from pPb data, using Z → samples, and techniques described in Ref. [23].These efficiencies are then compared to the corresponding efficiencies computed from simulations.In the case of trigger and reconstruction efficiencies, they are found to be consistent.The isolation criterion rejects more leptons in data, because the local activity of the underlying event is greater than in the simulation.To account for such discrepancies, the efficiency from W → ν simulation is multiplied by correction factors, which are determined as the ratio of the single-lepton efficiencies measured in Z → data to those estimated in simulations.The so-called "tag-and-probe" method used for this estimation is described in Ref. [27].These correction factors are computed in bins of η lab and for positively and negatively charged muons separately.In the electron case, the low statistical precision motivates a correction factor estimated for electrons and positrons combined.
The total systematic uncertainty in the lepton yields is estimated by adding the different contributions in quadrature.The η lab -dependent sources of systematic uncertainty arise from the method used for the estimation of multijet background (0.1-2.0% for muons, 0.5-3.8% for elec-trons), the normalization of the electroweak background (1-3% for muons and electrons), the efficiency correction factors (2.2-7.5% for muons, 2.6-7.4% for electrons), and the energy scale of electrons (0.1-2.0%).The uncertainty in the momentum scale of muons is found to be negligible.The integrated luminosity measurement uncertainty (3.5% [14]) affects only the W boson production cross sections and cancels in the asymmetry measurements, as does the additional global uncertainty arising from the efficiency of the filter rejecting pileup events (0.5% for both channels).Though the common electron/positron correction factors cancel, a residual systematic uncertainty of 3% is assigned to the charge asymmetry, based on simulation studies and η lab -integrated efficiencies determined from Z → e + e − data.No other systematic uncertainty cancellations are assumed for the asymmetry results.

Results
Figure 2 shows the production cross sections for pPb → W ± + X → ± ν + X as a function of the charged lepton pseudorapidity in the laboratory frame, with the lepton having p T > 25 GeV/c.The cross sections are determined by dividing the efficiency-corrected lepton yields by the integrated luminosity.Since the cross sections measured in the electron and muon channels are found to be in good agreement with each other, they are combined using the BLUE method [28].Figure 3 compares the combined cross sections with next-to-leading-order (NLO) perturbative QCD predictions provided by the authors of Ref. [4] using CT10 [29] proton parton distribution functions (PDF) without or with EPS09 [30] nPDF corrections, termed CT10 and CT10+EPS09, respectively.Their uncertainties are estimated as prescribed in Refs.[29,30].Table 1 gives the measured cross sections for each channel separately and combined, as a function of the lepton pseudorapidity, for positive and negative leptons.The theoretical predictions and their uncertainties (coming from the PDF set and from the renormalisation and factorisation scales) are also given.The agreement between the data and both theoretical predictions is within the uncertainties, although a small excess of W − candidates appears at negative η lab , i.e. in the Pb ion beam direction.
The comparison between the CT10 and CT10+EPS09 calculations shows that the predicted modifications of the PDFs are of the same order as the theoretical uncertainties.This indicates that cross sections alone lack discriminating power, and motivates the study of various asymmetries of the + and − cross sections.The interest in such asymmetries is twofold.First, some of the experimental (e.g.integrated luminosity) and theoretical (e.g.scale dependence) uncertainties cancel in such asymmetries.Second, the various asymmetries exhibit different sensitivities to the nuclear modifications of the PDFs, as discussed below.The lepton charge asymmetry, defined as (N + − N − )/(N + + N − ) with N ± being the efficiencycorrected lepton yields, is shown in Fig. 4, as a function of η lab , and compared to the theoretical predictions.For η lab > −1, both calculations reproduce the present measurements.For η lab < −1, however, the two calculations overpredict the asymmetry values.A possible physical origin of this disagreement could be a different modification of u and d quark distributions in nuclei.In proton-(anti)proton collisions, the W-boson charge asymmetry is known to be a sensitive probe of the down-to-up quark PDF ratio in a proton, d p /u p [20,31,32].Similarly, this asymmetry in pPb collisions measured in the lead fragmentation region (i.e.η lab < 0.465) probes these quark densities in a nucleon inside the lead nucleus.Assuming the standard isospin symmetry (u p = d n , u n = d p ), one can define a similar ratio, d p/A /u p/A = d p /u p × R d /R u , where R i are the nPDF ratios, R u ≡ u p/A /u p and R d ≡ d p/A /d p .The typical quark momentum fraction probed in the Pb nucleus is given by x M W / √ s NN × exp(−η lab + 0.465) (assuming that the W boson rapidity is similar to that of the lepton), therefore x 0.02 − 0.20 in the range −2 < η lab < 0. In most global fit analyses of the nPDFs (as in the case of EPS09), it is assumed that the nuclear ratios respect the isospin symmetry, namely R u = R d , essentially to minimise the number of free parameters in the fits.However, no physical reason prevents Table 1: Production cross section for pPb → W + X → ν + X for positively (top) and negatively (bottom) charged leptons of p T larger than 25 GeV/c, in nanobarns, as a function of the lepton pseudorapidity.Values are given first for muons and electrons separately, then combined.Quoted uncertainties are first statistical, then systematic.Theoretical predictions with (CT10+EPS09) and without (CT10) PDF nuclear modifications are also given, with their uncertainties.The global normalization uncertainty of 3.5% is not included in the listed uncertainties.103.6 +6.0 nuclear modifications to be different for up and down quark PDFs.For example, it is known that the shapes of the up and down quark distributions in protons are different [33].Furthermore, the present disparity between data and theory is unlikely to come from the proton PDF assumption, given the excellent agreement of lepton charge asymmetry measured in pp collisions by CMS [32] and ATLAS [34] with NLO calculations using CT10 parton densities.
A traditional way to probe nuclear parton densities is to compare the forward and backward W yields, that are respectively sensitive to the nPDFs at small and large x.The forward-backward asymmetries N ± (+η lab )/N ± (−η lab ) are shown in Fig. 5, separately for the positively and negatively charged leptons, and compared to the same predictions as mentioned above.Given the experimental accuracy and the magnitude of the differences between the two sets of predictions, the measurements have a potential to discriminate between them.However, although the negative lepton decay channel appears to slightly favour the CT10+EPS09 prediction over the CT10 calculation, the positive lepton channel does not, thus no firm conclusion can be drawn.
Another asymmetry variable, (N + (+η lab ) − N + (−η lab ))/(N − (+η lab ) − N − (−η lab )), was proposed in Ref. [4] to reach maximum sensitivity to nuclear modifications of PDFs.However, this asymmetry probability distribution shows a very non-Gaussian behaviour, when its denominator approaches zero, and its sign can be flipped within the uncertainty.A different asymmetry is proposed here, N (+η lab )/N (−η lab ), a forward-backward asymmetry of the charge-summed W bosons, which achieves a similar sensitivity.As in the case of the charge asymmetry, this asymmetry can be related to the nuclear modifications of the PDFs within the lead nucleus.Here, forward (backward) W boson production is sensitive to the PDFs of the sea quark at x ∼ 10 −3 (valence quark at x ∼ 10 −1 ) in the lead nucleus.Therefore, the forwardbackward ratio probes the small-x modification of the lead nucleus PDF (shadowing) over the large-x modifications (anti-shadowing).This asymmetry is shown in Fig. 6, and deviates from unmodified PDFs, more clearly favouring CT10+EPS09 over CT10.In order to quantify the agreement between the data and the expectation from the CT10 and CT10+EPS09 calculations, a χ 2 test is performed for each of the above (correlated) variables.The few correlations in experimental uncertainties described above, only relevant for W ± boson cross sections but not for asymmetries, are taken into account, as well as the correlations in theoretical uncertainties.The resulting χ 2 values and probabilities are given in Table 2.The CT10+EPS09 calculations provide a better description of the data, with still a relatively low probability for the lepton charge asymmetry, because of the backward region.

Summary
The first measurement of W boson production in pPb collisions has been reported, using the electron and muon decay modes for leptons of p T above 25 GeV/c and |η lab | < 2.4.The differential cross sections as a function of the lepton pseudorapidity agree with theoretical predictions assuming both unmodified (CT10) and modified (CT10+EPS09) nPDFs, except in the

Figure 1 :
Figure1: Missing transverse energy distribution for W + → µ + ν (left) and W + → e + ν (right) events within the −0.5 < η lab < 0.0 (top) and 2.0 < η lab < 2.4 (bottom) ranges.Binned fits to the data (red points) are performed with four contributions, stacked from bottom to top: multijet (QCD, blue), W + → τ + ν (brown), Z → (white) and W + → + ν (yellow).The η lab regions are defined such that the proton is moving towards positive η lab values.Error bars represent statistical uncertainties.The lower panels display the data divided by the result of the fit, with the band representing the statistical uncertainties on the sum of the fit components, for each E T / bin.

Figure 2 :
Figure 2: Production cross sections for W + → + ν (left) and W − → − ν (right), as a function of the lepton pseudorapidity.Error bars represent the statistical uncertainties, while brackets show statistical and systematic uncertainties summed in quadrature.The global luminosity uncertainty of ±3.5% is not included.To improve visibility, the muon (electron) measurements, in red circles (blue squares), have been shifted by −0.05 (+0.05) in pseudorapidity.

Figure 3 :
Figure3: Production cross sections for W + → + ν (left) and W − → − ν (right), as a function of the lepton pseudorapidity.Error bars represent the statistical uncertainties, while brackets show statistical and systematic uncertainties summed in quadrature.The global luminosity uncertainty of ±3.5% is not displayed.Theoretical predictions with (CT10+EPS09, dashed green line) and without (CT10, solid red line) PDF nuclear modifications are also shown, with the uncertainty bands.The bottom panels show the ratio of the data (black points) and CT10+EPS09 (dashed green line) to the CT10 baseline.All theory uncertainty bands include scale and PDF uncertainties, except the EPS09 of the bottom panels which only includes the EPS09 PDF uncertainties.

Figure 4 :Figure 5 :
Figure 4: Lepton charge asymmetry, (N + − N − )/(N + + N − ), as a function of the lepton pseudorapidity.Error bars represent the statistical uncertainties, while brackets show statistical and systematic uncertainties summed in quadrature.Theoretical predictions with (CT10+EPS09, dashed green line) and without (CT10, solid red line) PDF nuclear modifications are also shown, with their uncertainty bands.

Figure 6 :
Figure 6: The forward-backward asymmetry of charge-summed W bosons, as a function of the lepton pseudorapidity.Error bars represent the statistical uncertainties, while brackets show statistical and systematic uncertainties summed in quadrature.Theoretical predictions with (CT10+EPS09, dashed green line) and without (CT10, solid red line) PDF nuclear modifications are also shown, with their uncertainty bands.
The silicon tracker consists of 66 M pixel and 10 M strip sensor elements, and measures charged-particle trajectories in the pseudorapidity range |η lab | < 2.5.Outside of the solenoid, muons are detected in the |η lab | < 2.4 range, with gas-ionization detector planes based on three technologies: drift tubes, cathode strip chambers, and resistive-plate chambers.Electrons are identified in the ECAL, which is made of 75 848 lead tungstate crystals and covers |η lab | < 1.48 in the barrel and 1.48 < |η lab | < 3.00 in the two endcap regions.The CMS apparatus also has extensive forward calorimetry, including two steel/quartz-fiber Cherenkov hadron forward (HF) calorimeters, which cover the 2.9 < |η lab | < 5.2 range.For online event selection, CMS uses a two-level trigger system.

Table 2 :
Values of the χ 2 test between the measurements and the theoretical predictions, with (CT10+EPS09) or without (CT10) nuclear modifications of the PDFs.The probability (Prob.) to measure a value greater to that measured in data is also given for ten degrees of freedom in the case of the first three observables and five degrees of freedom for the three others observables.Pb ion beam direction), where a hint of an enhancement is seen for the W − bosons.In the same region, the related lepton charge asymmetry deviates slightly from the predictions, something that could potentially arise from different nuclear modifications of the up and down quark PDFs.In a related observation, forward-backward asymmetries show a deviation from unmodified PDFs.Taken together, these measurements show the need for including W boson data in nuclear parton distribution global fits.and Industrial Research, India; the HOMING PLUS programme of the Foundation for Polish Science, cofinanced from European Union, Regional Development Fund; the Compagnia di San Paolo (Torino); the Consorzio per la Fisica (Trieste); MIUR project 20108T4XTM (Italy); the Thalis and Aristeia programmes cofinanced by EU-ESF and the Greek NSRF; and the National Priorities Research Program by Qatar National Research Fund.