Differential cross section measurements for the production of a W boson in association with jets in proton-proton collisions at sqrt(s) = 7 TeV

Measurements are reported of differential cross sections for the production of a W boson, which decays into a muon and a neutrino, in association with jets, as a function of several variables, including the transverse momenta (pt) and pseudorapidities of the four leading jets, the scalar sum of jet transverse momenta (HT), and the difference in azimuthal angle between the directions of each jet and the muon. The data sample of pp collisions at a centre-of-mass energy of 7 TeV was collected with the CMS detector at the LHC and corresponds to an integrated luminosity of 5.0 inverse femtobarns. The measured cross sections are compared to predictions from Monte Carlo generators, MADGRAPH + PYTHIA and SHERPA, and to next-to-leading-order calculations from BLACKHAT + SHERPA. The differential cross sections are found to be in agreement with the predictions, apart from the pt distributions of the leading jets at high pt values, the distributions of the HT at high-HT and low jet multiplicity, and the distribution of the difference in azimuthal angle between the leading jet and the muon at low values.


Introduction
This paper reports measurements of the total and differential cross sections for W-boson production in association with jets at the LHC.Measurements of the production of vector bosons in association with jets are fundamental tests of perturbative quantum chromodynamics (pQCD).The W+jets processes also provide the main background to other, much rarer, standard model (SM) processes, such as tt [1] and single top-quark production [2], and to Higgs boson production and a variety of physics processes beyond the SM.Searches for phenomena beyond the SM are often limited by the uncertainty in the theoretical cross sections for W (and Z) + jets processes at high momentum scales and large jet multiplicities.It is crucial to perform precision measurements of W+jets production at the LHC.Leptonic decay modes of the vector boson are often used in the measurement of SM processes and in searches for new physics because they provide clean signatures with relatively low background.This paper focuses on the production of a W boson decaying into a muon and a neutrino, as part of a final-state topology characterised by one high-transverse-momentum (p T ) isolated muon, significant missing transverse energy (E miss T ), and one or more jets.The cross sections are measured as a function of the inclusive and exclusive jet multiplicities for up to six jets.Differential cross sections are measured for different inclusive jet multiplicities as a function of the transverse momentum and the pseudorapidity (η) of the jets, where η = − ln[tan(θ/2)], and θ is the polar angle measured with respect to the anticlockwise beam direction.The cross sections are also measured as a function of the difference in azimuthal angle between the direction of each jet and that of the muon, and of H T , which is defined as the scalar sum of the p T of all jets with p T > 30 GeV and |η| < 2.4.It is important to study the distribution of the jet p T and the observable H T because they are sensitive to higher order corrections, and are often used to discriminate against background in searches for signatures of physics beyond the SM.Additionally, H T is commonly used as a component of the factorisation and renormalisation scale in event generators simulating W+jets production at the LHC.Finally, the η distributions of jets and the azimuthal separations between the jets and the muon are also important, because they are sensitive to the modelling of parton emission.
These measurements use proton-proton (pp) collision data at a centre-of-mass energy of √ s = 7 TeV recorded with the CMS detector at the LHC in 2011 and corresponding to an integrated luminosity of 5.0 ± 0.1 fb −1 [3].Previous studies of leptonic decay modes at the LHC have measured the cross sections and cross section ratios for W boson production in association with jets in pp collisions with an integrated luminosity of 36 pb −1 at √ s = 7 TeV with the AT-LAS [4] and CMS [5] detectors.Measurements have also been made with pp collisions with the D0 detector [6, 7] at the Tevatron collider for integrated luminosities up to 4.2 fb −1 , as well as with the CDF detector [8] for an integrated luminosity of 320 pb −1 .In order to perform a differential measurement of the W+jets cross section, a high-purity sample of W → µν events is selected and the kinematic distributions are corrected to the particle level by means of regularised unfolding [9].This procedure corrects a measured observable for the effects of detector response, finite experimental resolutions, acceptance, and efficiencies, and therefore allows for direct comparison with theoretical predictions.The measured differential cross sections are compared to the predictions of generators such as MADGRAPH 5.1.1 [10] interfaced with PYTHIA 6.426 [11], SHERPA 1.4.0 [12][13][14][15], and BLACKHAT [16,17], interfaced to SHERPA.The BLACKHAT+SHERPA samples [18] provide parton-level predictions of W + n (n = 1-4) jets at next-to-leading order (NLO), while the MADGRAPH+PYTHIA and SHERPA samples provide final states after hadronisation.
The paper proceeds as follows: Section 2 presents the CMS detector.Section 3 describes the Monte Carlo (MC) event generators, as well as the data samples used for the analysis.The identification criteria for the final-state objects (leptons and jets) and the selection of the W → µν+ jets events are presented in Section 4. Section 5 describes the modelling of instrumental backgrounds and irreducible physics backgrounds.The procedure used for unfolding is detailed in Section 6, and Section 7 describes the systematic uncertainties.Finally, the unfolded distributions are presented in Section 8 and compared to theoretical predictions, and Section 9 summarises the results.

The CMS detector
The CMS detector, presented in detail elsewhere [19], can be described with a cylindrical coordinate system with the +z axis directed along the anticlockwise beam axis.The detector consists of an inner tracking system and calorimeters (electromagnetic, ECAL, and hadron, HCAL) surrounded by a 3.8 T solenoid.The inner tracking system consists of a silicon pixel and strip tracker, providing the required granularity and precision for the reconstruction of vertices of charged particles in the range 0 ≤ φ < 2π in azimuth and |η| < 2.5.The crystal ECAL and the brass/scintillator sampling HCAL are used to measure the energies of photons, electrons, and hadrons within |η| < 3.0.The HCAL, when combined with the ECAL, measures jets with a resolution ∆E/E ≈ 100%/ E [GeV] ⊕ 5% [20].The three muon systems surrounding the solenoid cover a region |η| < 2.4 and are composed of drift tubes in the barrel region (|η| < 1.2), cathode strip chambers in the endcaps (0.9 < |η| < 2.4), and resistive-plate chambers in both the barrel region and the endcaps (|η| < 1.6).Events are recorded based on a trigger decision using information from the CMS detector subsystems.The first level (L1) of the CMS trigger system, composed of custom hardware processors, uses information from the calorimeters and muon detectors to select the most interesting events in a fixed time interval of less than 4 µs.The high-level trigger (HLT) processor further decreases the event rate from 100 kHz at L1 to roughly 300 Hz.

Data and simulation samples
Events are retained if they pass a trigger requiring one isolated muon with p T > 24 GeV and |η| < 2.1.Signal and background simulated samples are produced and fully reconstructed using a simulation of the CMS detector based on GEANT4 [21], and simulated events are required to pass an emulation of the trigger requirements applied to the data.These simulations include multiple collisions in a single bunch crossing (pileup) by matching the number of pileup events delivered by the LHC during the data-taking periods.
A W → ν+jets signal sample is generated with MADGRAPH 5.1.1 and is used to determine the detector response in the unfolding procedure described in Section 6. Parton showering and hadronisation of the MADGRAPH samples is performed with PYTHIA 6.424 using the Z2 tune [22].The detector response is also determined using a different W+jets event sample generated with SHERPA 1.3.0[12][13][14][15], and is used in the evaluation of systematic uncertainties due to the unfolding of the data.
The main sources of background are the production of tt, single top-quark, Z/γ * +jets, dibosons (ZZ/WZ/WW) + jets, and multijet production.With the exception of multijet production, all backgrounds are estimated from simulation.The simulated samples of tt and Z/γ * +jets are generated with MADGRAPH 5.1.1;single top-quark samples (s-, t-, and tW-channels) are generated with POWHEG [23][24][25][26]; VV samples, where V represents either a W boson or a Z boson, are generated with PYTHIA version 6.424 using the Z2 tune [22].The simulations with MAD-GRAPH and PYTHIA use the CTEQ6L1 parton distribution functions (PDF) [27].The simulation with SHERPA uses the CTEQ6.6mPDF, and the simulations with POWHEG use the CTEQ6m PDF.

Object identification and event selection
Muon candidates are reconstructed as tracks in the muon system that are matched to tracks reconstructed in the inner tracking system [34].Muon candidates are required to have p T > 25 GeV, and to be reconstructed within the fiducial volume used for the high-level trigger muon selection, i.e. within |η| < 2.1.This ensures that the offline event selection requirements are as stringent as the trigger.In addition, an isolation requirement is applied to the muon candidates by demanding that the relative isolation is less than 0.15, where the relative isolation is defined as the sum of the transverse energy deposited in the calorimeters (ECAL and HCAL) and of the p T of charged particles measured with the tracker in a cone of ∆R = √ (∆φ) 2 + (∆η) 2 = 0.3 around the muon candidate track (excluding this track), divided by the muon candidate p T .To ensure a precise measurement of the transverse impact parameter of the muon track relative to the interaction point, only muon candidates with tracks containing more than 10 hits in the silicon tracker and at least one hit in the pixel detector are considered.To reject muons from cosmic rays, the transverse impact parameter of the muon candidate with respect to the primary vertex is required to be less than 2 mm.
Jets and E miss T are reconstructed using the CMS particle-flow algorithm [35,36].Jets are reconstructed using the anti-k T [37,38] algorithm with a distance parameter of 0.5.The jet energy is calibrated using the p T balance of dijet and γ+jet events [39] to account for the following effects: nonuniformity and nonlinearity of the ECAL and HCAL energy response to neutral hadrons, the presence of extra particles from pileup interactions, the thresholds used in jet constituent selection, reconstruction inefficiencies, and possible biases introduced by the clustering algorithm.Only jets with p T > 30 GeV, |η| < 2.4, and a spatial separation of ∆R > 0.5 from the muon are considered.To reduce the contamination from pileup jets, jets are required to be associated to the same primary vertex as the muon.The vertex associated to each jet is the one that has the largest number of p T -weighted tracks in common with the jet.The contamination from pileup jets is estimated with the signal simulation, with pileup events simulated with PYTHIA, and found to be less than 1%.
The missing momentum vector, p miss T , is defined as the negative of the vectorial sum of the transverse momenta of the particles reconstructed with the particle-flow algorithm, and the E miss T is defined as the magnitude of the p miss T vector.The measurement of the E miss T in simulation is sensitive to the modelling of the calorimeter response and resolution and to the description of the underlying event.To account for these effects, the E miss T in W+jets simulation is corrected for the differences in the detector response between data and simulation, using a method detailed in Ref. [40].A recoil energy correction is applied to the W+jets simulation on an event-by-event basis, using a sample of Z → µµ events in data and simulation.The transverse recoil vector, defined as the negative vector sum of the missing transverse energy and the transverse momenta of the lepton(s), is divided into components parallel and perpendicular to the boson direction.The mean and the width of the transverse recoil vector components are parameterised as a function of the Z-boson p T in data and simulation.The ratio of the data and simulation parameterisations is used to adjust the transverse recoil vector components in each simulated event, and a new E miss T is computed using the corrected recoil components.
Events are required to contain exactly one muon satisfying the conditions described above and one or more jets with p T > 30 GeV.Events are required to have M T > 50 GeV, where M T , the transverse mass of the muon and missing transverse energy, is defined as where p µ T is the muon p T and ∆φ is the difference in azimuthal angle between the muon momentum direction and the p miss T vector.

Estimation of the backgrounds and selection efficiencies
All background sources except for the multijet production are modelled with simulation.The simulated event samples are corrected for differences between data and simulation in muon identification efficiencies and event trigger efficiency.A "tag-and-probe" method [34] is used to determine the differences between simulation and data for the efficiency of the trigger and for the muon identification and isolation criteria.This method uses Z → µµ events from both data and simulated samples where the "tag" muon is required to pass the identification and isolation criteria.The efficiency measurements use the "probe" muon, which is required to pass minimal quality criteria.Trigger efficiency corrections are determined as a function of the muon η, and are in general less than 5%.Muon isolation and identification efficiency corrections are determined as a function of the muon p T and η, and are generally less than 2%.Corrections to the simulation are applied on an event-by-event basis in the form of event weights.
The dominant background to W+jets production is tt production, which has a larger total contribution than that of the W+jets signal in events with four or more jets.In order to reduce the level of tt contamination, a veto is applied to events with one or more b-tagged jets.Heavyflavour tagging is based on a tag algorithm [41] that exploits the long lifetime of b-quark hadrons.This algorithm calculates the signed impact parameter significance of all tracks in the jet that satisfy high-quality reconstruction and purity criteria, and orders the tracks by decreasing significance.The impact parameter significance is defined as the ratio of the impact parameter to its estimated uncertainty.For jets with two or more significant tracks, a highefficiency b-quark discriminator is defined as the significance of the second-most significant track.The size of the tt background is illustrated in Fig. 1, before and after the implementation of the b-jet veto, using the event selection described in Section 4. The expected contributions for the different processes in Fig. 1 are shown as a function of the jet multiplicity, along with the observed data.Differences in the tagging and mistagging rates between data and simulation are measured as a function of the jet p T in multijet and tt events [41], and are used to correct the tagging rates of the jets in simulation.For jet multiplicities of 1 to 6, the b-jet veto eliminates 44 -84% of the predicted tt background, while eliminating 3 -26% of the predicted W+jets signal.
The shape and normalisations of the Z/γ * +jets and tt predictions are cross-checked in selected data samples.The Z+jets background is compared to data in a Z-boson-dominated data sample that requires two well-identified, isolated muons.The tt background is compared to data in a control region requiring at least two b-tagged jets.Background estimations from simulation and from data control samples agree within the uncertainties described in Section 7.
The multijet background is estimated using a data control sample with an inverted muon isolation requirement.Contributions from other sources to the multijet control region are subtracted, including the dominant contribution from W+jets.In order to improve the estimation of W+jets in the multijet control region, the W+jets contribution is first normalised to data in the M T > 50 GeV region with the muon isolation condition applied.The multijet contribution is then rescaled according to a muon misidentification rate measured in the multijet enriched sideband region of M T < 50 GeV.The multijet contribution to the total background estimate is 32.7% (1.9%) for selections requiring exactly 1 (4) jets, or 2.6% (0.3%) of the total SM prediction.
Events/bin  The W+jets contribution is modelled with MADGRAPH 5.1.1+PYTHIA6.424.The solid band indicates the total statistical and systematic uncertainty in the W+jets signal and background predictions, as detailed in Section 7.This includes uncertainties in the jet energy scale and resolution, the muon momentum scale and resolution, the pileup modelling, the b-tagging correction factors, the normalisations of the simulations, and the efficiencies of reconstruction, identification, and trigger acceptance.A substantial reduction in the expected tt background is observed in the right plot.

The unfolding procedure
For the measurement of cross sections, the particle level is defined by a W boson, which decays into a muon and a muon neutrino, produced in association with one or more jets.Decays of the W boson to a tau and a tau neutrino, resulting in a final state with a muon, are not included in the particle level definition.Kinematic thresholds on the particle-level muon, M T , and jets are identical to those applied to the reconstructed objects.The particle-level selection also includes the requirement of exactly one muon with p T > 25 GeV and |η| < 2.1, and M T > 50 GeV.The particle-level E miss T is defined as the negative of the vectorial sum of the transverse momenta of all visible final state particles.To account for final-state radiation, the momenta of all photons in a cone of ∆R < 0.1 around the muon are added to that of the muon.Jets are clustered using the anti-k T [37] algorithm with a distance parameter of 0.5.Clustering is performed using all particles after decay and fragmentation, excluding neutrinos and the muon from the W-boson decay.Additionally, jets are required to have p T > 30 GeV and |η| < 2.4, and to be separated from the muon by ∆R > 0.5.
The reconstructed distributions are corrected to the particle level with the method of regularised singular value decomposition (SVD) [9] unfolding, using the ROOUNFOLD toolkit [42].For each distribution, the total background, including the multijet estimate from data and all simulated processes except the W boson signal, is subtracted from the data before unfolding.A response matrix, defining the migration probability between the particle-level and reconstructed quantities, as well as the overall reconstruction efficiency, is computed using W+jets events simulated with MADGRAPH+PYTHIA.For a given particle-level quantity Q with a corresponding reconstructed quantity Q , the migration probability from an interval a < Q < b to an interval c < Q < d is defined as the fraction of events with a < The unfolding of the jet multiplicity is performed with a response defined by the number of particle-level jets versus the number of reconstructed jets.For particle-level jet multiplicities of 1 to 6, 4% to 51% of simulated events exhibit migration to different values of reconstructed jet multiplicity.The unfolding of the kinematic distributions of the n th jet is performed with a response defined by the kinematic quantity of the n th -highest-p T particle-level jet versus that of the n th -highest-p T reconstructed jet.To achieve a full migration from the selection of reconstructed events to the particle-level phase space, no matching between reconstructed and particle-level jets is applied.Since a reconstructed event may contain a final state muon from the W → τν process, the contamination from the W → τν process is estimated with simulation to be at the 1% level, and is removed during the unfolding procedure.
The b-jet veto is treated as an overall event selection condition in the same manner as the basic kinematic conditions on the M T and muon p T .Events failing this condition are treated as nonreconstructed in the unfolding response, so that the cross section obtained after unfolding is valid for W-boson decays with associated jets of any flavour.

Systematic uncertainties
The sources of systematic uncertainties considered in this analysis are described below.The entirety of the unfolding procedure is repeated for each systematic variation, and the unfolded data results with these variations are compared with the central (unvaried) results to extract the uncertainties in the unfolded data distributions.
In most distributions, the dominant sources of systematic uncertainty include the jet energy scale and resolution uncertainties, which affect the shape of all reconstructed distributions, as well as the overall event acceptance.The jet energy scale uncertainties are estimated by assigning a p T -and η-dependent uncertainty in jet energy corrections as discussed in Ref. [39], and by varying the jet p T by the magnitude of the uncertainty.The uncertainties in jet energy resolution are assessed by increasing the p T difference between the reconstructed and particlelevel jets by an η-dependent value [39].The jet energy uncertainties are determined by varying the p T of the jets in data rather than in simulation.
Muon momentum scale and resolution uncertainties also introduce uncertainties in the overall event acceptance.A muon momentum scale uncertainty of 0.2% and a muon momentum resolution uncertainty of 0.6% are assumed [34].The effects of these uncertainties are assessed by directly varying the momentum scale and randomly fluctuating the muon momentum in the simulation.
Variations for uncertainties in the energy and momentum scales and resolutions affect the size and shape of the background distribution to be subtracted from the data distribution, as well as the acceptance of W+jets simulated events, which define the response matrix used for unfolding.The variations are also propagated to the measurement of E miss T , which affects the acceptance of the M T > 50 GeV requirement.
Another important source of systematic uncertainty is the choice of the generator used in the unfolding procedure.The size of this uncertainty is assessed by repeating the unfolding procedure with a response trained on a separate simulated sample generated with SHERPA 1.3.0.The absolute value of the difference between the data unfolded with a response matrix trained on SHERPA and with a response matrix trained on MADGRAPH is treated as a symmetric uncertainty on the measurement.
Other minor sources of systematic uncertainty include the uncertainties in the background normalisation, the b-tagging efficiency, the modelling of the Wb contribution in the signal simulation, integrated luminosity, the pileup modelling, the trigger and object identification efficiencies, and the finite number of simulated events used to construct the response.Background normalisation uncertainties are determined by varying the cross sections of the backgrounds within their theoretical uncertainties.These uncertainties are 4.4% for tt, 4.3% for Z+jets, 6% for single top-quark and 4% for diboson production.The estimate of the multijet background has an uncertainty based on the limited number of events in the multijet sample and in the control regions where the multijet sample normalisation is calculated.For the b-tagging algorithm used to veto events containing b jets, uncertainties in the data/simulation ratio of the b-tagging efficiencies are applied.For jets with p T > 30 GeV, these uncertainties range from 3.1% to 10.5%.An additional uncertainty is ascribed to the normalisation of the Wb content in the simulation by examining the agreement between data and simulation as a function of jet multiplicity in a control region requiring exactly one b-tagged jet.The uncertainty in the integrated luminosity is 2.2% [3].An uncertainty in the modelling of pileup in simulation is determined by varying the number of simulated pileup interactions by 5% to account for the uncertainty in the luminosity and the uncertainty in the total inelastic cross section [43], as determined by a comparison of the number of reconstructed vertices in Z → µµ events in data and simulation.Uncertainties in the differences between data and simulation efficiencies of the trigger, muon isolation, and muon identification criteria are generally less than 1%.An additional uncertainty due to the finite number of simulated events used to construct the response matrix is calculated by randomly varying the content of the response matrix according to a Poisson uncertainty in each bin.
The effect of the systematic variations on the measured cross section as a function of the exclusive jet multiplicity is illustrated in Fig. 2. The uncertainties given in Fig. 2 are the total uncertainty for each jet multiplicity.The corresponding ranges of systematic uncertainty across bins of jet p T are given in Table 1.

Results
The cross sections for exclusive and inclusive jet multiplicities are given in Fig. 3.In Figs.4-7 the differential cross sections are presented.The measured W+jets cross sections are compared to the predictions from several generators with an overall normalisation determined by the NNLO inclusive cross sections calculated with FEWZ [28].We consider W+jets signal processes generated with MADGRAPH 5.1.1 using the CTEQ6L1 PDF set, with SHERPA 1.4.0 using the CT10 [44,45] PDF set, and with BLACKHAT+SHERPA [16] using the CT10 PDF set.The SHERPA sample is a separate sample from that used for the evaluation of uncertainties in Section 7. The MADGRAPH and SHERPA predictions provide leading-order (LO) matrix element (ME) calcu- The dominant systematic uncertainties in the measurement of the W+jets cross section as a function of the exclusive jet multiplicity.The systematic uncertainties displayed include the jet energy scale and resolution (JES, JER), the choice of generator used in the unfolding procedure (Generator), the statistical uncertainty propagated through the unfolding procedure (Statistical), the uncertainty due to a finite number of simulated events used to construct the response (MC stat.), and all other systematic uncertainties (All other) detailed in Section 7. The uncertainties presented here correspond to the weighted average of the values shown in Table 1.
Table 1: Ranges of uncertainties for the measurement of dσ/dp T of the n th jet in events with n or more jets.The uncertainties displayed include the statistical uncertainty propagated through the unfolding procedure (Statistical), the jet energy scale and resolution (JES, JER), the choice of generator used in the unfolding procedure (Generator), the uncertainty due to a finite number of simulated events used to construct the response (MC Stat.), and all other systematic uncertainties (Other) detailed in Section 7. lations at each jet multiplicity, which are then combined into inclusive samples by matching the ME partons to particle jets.Parton showering (PS) and hadronisation of the MADGRAPH sample is performed with PYTHIA 6.426 using the Z2 tune.The MADGRAPH +PYTHIA calculation includes the production of up to four partons.The matching between jets and partons is performed with the kT-MLM matching scheme [46] with a 20 GeV matching threshold.The factorisation scale for each event is chosen to be the transverse mass computed after k T -clustering of the event down to a 2→2 topology.The renormalisation scale for the event is the k T computed at each vertex splitting.The predictions from SHERPA include the production of up to four partons.The matching between jets and partons is performed with the CKKW matching scheme [46], and the default factorisation and renormalisation scales are used.
The predictions from MADGRAPH+PYTHIA and SHERPA are shown with statistical uncertainties only.These MADGRAPH+PYTHIA and SHERPA samples are processed through the RIVET toolkit [47] in order to create particle level distributions, which can be compared with the unfolded data.The BLACKHAT+SHERPA samples represent fixed-order predictions at the level of ME partons of W + n jets at NLO accuracy, for n = 1, 2, 3, and 4 jets.Each measured distribution for a given inclusive jet multiplicity is compared with the corresponding fixed-order prediction from BLACKHAT+SHERPA.The choice of renormalisation and factorisation scales for BLACK-HAT+SHERPA is Ĥ T /2, where Ĥ T ≡ ∑ m p m T + E W T , and m represents the final state partons.Before comparing to data, a nonperturbative correction is applied to the BLACKHAT+SHERPA distributions to account for the effects of multiple-parton interactions and hadronisation.The nonperturbative correction is determined using MADGRAPH 5.1.1 interfaced to PYTHIA 6.426 and turning on and off the hadronisation and multiple-parton interactions.The magnitude of the nonperturbative correction is typically 1-5%, and is calculated for each bin of each measured distribution.The model dependence of the nonperturbative correction is negligible [48].The BLACKHAT+SHERPA prediction also includes uncertainties due to the PDF and variations of the factorisation and renormalisation scales.The nominal prediction is given by the central value of the CT10 PDF set, and the PDF uncertainty considers the envelope of the error sets of CT10, MSTW2008nlo68cl [49], and NNPDF2.1 [50] according to the PDF4LHC prescription [51,52].The factorisation and renormalisation scale uncertainty is determined by varying the scales simultaneously by a factor 0.5 or 2.0.
The unfolded exclusive and inclusive jet multiplicity distributions, shown in Fig. 3, are found to be in agreement, within uncertainties, with the predictions of the generators and with the NLO calculation of BLACKHAT+SHERPA.Table 2 details the measured cross sections as a function of the inclusive and exclusive jet multiplicity.
The jet p T unfolded distributions for inclusive jet multiplicities from 1 to 4 are shown in Fig. 4. The predictions of BLACKHAT+SHERPA are in agreement with the measured distributions within the systematic uncertainties, while MADGRAPH+PYTHIA is observed to overestimate the yields up to 54% (45%) for the first (second) leading jet p T distributions at high-p T values.The predictions from SHERPA are found to agree well for the second-, third-, and fourth-leading jet p T distributions, while an excess of slightly more than one standard deviation can be seen at high-p T values for the leading jet p T distribution.Similar observations hold for MADGRAPH+PYTHIA and SHERPA predictions in the H T distributions for inclusive jet multiplicities of 1-4, as shown in Fig. 5. Since the BLACKHAT+SHERPA NLO prediction for H T (≥1 jet) is a fixed-order prediction with up to two real partons, contributions from higher jet multiplicities are missing, which results in an underestimation in the tail of the distribution [53].Similar observations have been made with W+jets measurements at D0 [7] and ATLAS [4].In general, SHERPA models the H T distributions better than other generators.
The distributions of the jet η and of the difference in azimuthal angle between each jet and the muon are shown in Figs. 6 and 7, respectively.The measurements of the jet η agree with predictions from all generators, with MADGRAPH+PYTHIA and BLACKHAT+SHERPA performing best.The measurements of the ∆φ between the leading jet and the muon are underestimated by as much as 38% by BLACKHAT+SHERPA, with similar, but smaller, underestimations in predictions from MADGRAPH+PYTHIA and SHERPA.
Examples of the variation in the BLACKHAT+SHERPA prediction due to the choice of PDF are given in Fig. 8, in which the predictions with the MSTW2008nlo68cl and NNPDF2.1 PDF sets are compared to the prediction from the CT10 PDF set, and to the measurements from data.The distributions determined with the different PDF sets are consistent with one another.The cross section measurement for the exclusive and inclusive jet multiplicities, compared to the predictions of MADGRAPH 5.1.1 + PYTHIA 6.426, SHERPA 1.4.0, and BLACK-HAT+SHERPA (corrected for hadronisation and multiple-parton interactions).Black circular markers with the grey hatched band represent the unfolded data measurement and its uncertainty.Overlaid are the predictions together with their statistical uncertainties (Theory stat.).The BLACKHAT+SHERPA uncertainty also contains theoretical systematic uncertainties (Theory syst.)described in Section 8.The lower plots show the ratio of each prediction to the unfolded data.

Summary
Measurements of the cross sections and differential cross sections for a W boson produced in association with jets in pp collisions at a centre-of-mass energy of 7 TeV have been presented.The data were collected with the CMS detector during the 2011 pp run of the LHC, and correspond to an integrated luminosity of 5.0 fb −1 .Cross sections have been determined using the muon decay mode of the W boson and were presented as functions of the jet multiplicity, the transverse momenta and pseudorapidities of the four leading jets, and the difference in azimuthal angle between each jet and the muon.The results, corrected for all detector effects by means of regularised unfolding, have been compared with particle-level simulated predictions from pQCD.
Predictions from generators, MADGRAPH+PYTHIA and SHERPA, and NLO calculations from BLACKHAT+SHERPA, describe the jet multiplicity within the uncertainties.The cross section as a function of the p T of the leading jet is overestimated by MADGRAPH+PYTHIA and SHERPA, especially at high-p T .Some overestimation from MADGRAPH+PYTHIA can also be observed in the second-and third-leading jet p T distributions.The cross sections as a function of p T predicted by BLACKHAT+SHERPA agree with the measurements within uncertainties.Similar levels of agreement have been observed in the distributions of H T for N jet ≥ 1, 2, 3 and 4. The predictions from BLACKHAT+SHERPA underestimate the measurement of the cross section as  The differential cross section measurement for the leading four jets transverse momentum, compared to the predictions of MADGRAPH 5.1.1 + PYTHIA 6.426, SHERPA 1.4.0, and BLACKHAT+SHERPA (corrected for hadronisation and multiple-parton interactions).Black circular markers with the grey hatched band represent the unfolded data measurement and its uncertainty.Overlaid are the predictions together with their statistical uncertainties (Theory stat.).The BLACKHAT+SHERPA uncertainty also contains theoretical systematic uncertainties (Theory syst.)described in Section 8.The lower plots show the ratio of each prediction to the unfolded data.The differential cross section measurement for H T for inclusive jet multiplicities 1-4, compared to the predictions of MADGRAPH 5.1.1 + PYTHIA 6.426, SHERPA 1.4.0, and BLACK-HAT+SHERPA (corrected for hadronisation and multiple-parton interactions).Black circular markers with the grey hatched band represent the unfolded data measurement and its uncertainty.Overlaid are the predictions together with their statistical uncertainties (Theory stat.).The BLACKHAT+SHERPA uncertainty also contains theoretical systematic uncertainties (Theory syst.)described in Section 8.The lower plots show the ratio of each prediction to the unfolded data.The differential cross section measurement for the pseudorapidity of the four leading jets, compared to the predictions of MADGRAPH 5.1.1 + PYTHIA 6.426, SHERPA 1.4.0, and BLACKHAT+SHERPA (corrected for hadronisation and multiple-parton interactions).Black circular markers with the grey hatched band represent the unfolded data measurement and its uncertainty.Overlaid are the predictions together with their statistical uncertainties (Theory stat.).The BLACKHAT+SHERPA uncertainty also contains theoretical systematic uncertainties (Theory syst.)described in Section 8.The lower plots show the ratio of each prediction to the unfolded data.The differential cross section measurement in ∆φ(jet n , µ), for n = 1 -4, compared to the predictions of MADGRAPH 5.1.1 + PYTHIA 6.426, SHERPA 1.4.0, and BLACKHAT+SHERPA (corrected for hadronisation and multiple-parton interactions).Black circular markers with the grey hatched band represent the unfolded data measurement and its uncertainty.Overlaid are the predictions together with their statistical uncertainties (Theory stat.).The BLACK-HAT+SHERPA uncertainty also contains theoretical systematic uncertainties (Theory syst.)described in Section 8.The lower plots show the ratio of each prediction to the unfolded data.a function of H T for N jet ≥ 1, since the contribution from W+≥3 jets is missing from an NLO prediction of W+≥1 jet.The distributions of ∆φ between the leading jet and the muon are underestimated by all predictions for ∆φ values near zero, with the largest disagreement visible in BLACKHAT+SHERPA.The distributions of ∆φ between the second-, third-, and fourth-leading jets and the muon agree with all predictions within uncertainties.No significant disagreement was found in the distributions of η of the four leading jets.

Figure 1 :
Figure 1: The jet multiplicity in data and simulation before (left) and after (right) the b-jet veto.The W+jets contribution is modelled with MADGRAPH 5.1.1+PYTHIA6.424.The solid band indicates the total statistical and systematic uncertainty in the W+jets signal and background predictions, as detailed in Section 7.This includes uncertainties in the jet energy scale and resolution, the muon momentum scale and resolution, the pileup modelling, the b-tagging correction factors, the normalisations of the simulations, and the efficiencies of reconstruction, identification, and trigger acceptance.A substantial reduction in the expected tt background is observed in the right plot.

Figure 2 :
Figure2: The dominant systematic uncertainties in the measurement of the W+jets cross section as a function of the exclusive jet multiplicity.The systematic uncertainties displayed include the jet energy scale and resolution (JES, JER), the choice of generator used in the unfolding procedure (Generator), the statistical uncertainty propagated through the unfolding procedure (Statistical), the uncertainty due to a finite number of simulated events used to construct the response (MC stat.), and all other systematic uncertainties (All other) detailed in Section 7. The uncertainties presented here correspond to the weighted average of the values shown in Table1.

Figure 3 :
Figure3: The cross section measurement for the exclusive and inclusive jet multiplicities, compared to the predictions of MADGRAPH 5.1.1 + PYTHIA 6.426, SHERPA 1.4.0, and BLACK-HAT+SHERPA (corrected for hadronisation and multiple-parton interactions).Black circular markers with the grey hatched band represent the unfolded data measurement and its uncertainty.Overlaid are the predictions together with their statistical uncertainties (Theory stat.).The BLACKHAT+SHERPA uncertainty also contains theoretical systematic uncertainties (Theory syst.)described in Section 8.The lower plots show the ratio of each prediction to the unfolded data.

Figure 4 :
Figure4: The differential cross section measurement for the leading four jets transverse momentum, compared to the predictions of MADGRAPH 5.1.1 + PYTHIA 6.426, SHERPA 1.4.0, and BLACKHAT+SHERPA (corrected for hadronisation and multiple-parton interactions).Black circular markers with the grey hatched band represent the unfolded data measurement and its uncertainty.Overlaid are the predictions together with their statistical uncertainties (Theory stat.).The BLACKHAT+SHERPA uncertainty also contains theoretical systematic uncertainties (Theory syst.)described in Section 8.The lower plots show the ratio of each prediction to the unfolded data.

Figure 5 :
Figure5: The differential cross section measurement for H T for inclusive jet multiplicities 1-4, compared to the predictions of MADGRAPH 5.1.1 + PYTHIA 6.426, SHERPA 1.4.0, and BLACK-HAT+SHERPA (corrected for hadronisation and multiple-parton interactions).Black circular markers with the grey hatched band represent the unfolded data measurement and its uncertainty.Overlaid are the predictions together with their statistical uncertainties (Theory stat.).The BLACKHAT+SHERPA uncertainty also contains theoretical systematic uncertainties (Theory syst.)described in Section 8.The lower plots show the ratio of each prediction to the unfolded data.

Figure 6 :
Figure6: The differential cross section measurement for the pseudorapidity of the four leading jets, compared to the predictions of MADGRAPH 5.1.1 + PYTHIA 6.426, SHERPA 1.4.0, and BLACKHAT+SHERPA (corrected for hadronisation and multiple-parton interactions).Black circular markers with the grey hatched band represent the unfolded data measurement and its uncertainty.Overlaid are the predictions together with their statistical uncertainties (Theory stat.).The BLACKHAT+SHERPA uncertainty also contains theoretical systematic uncertainties (Theory syst.)described in Section 8.The lower plots show the ratio of each prediction to the unfolded data.

Figure 7 :
Figure7: The differential cross section measurement in ∆φ(jet n , µ), for n = 1 -4, compared to the predictions of MADGRAPH 5.1.1 + PYTHIA 6.426, SHERPA 1.4.0, and BLACKHAT+SHERPA (corrected for hadronisation and multiple-parton interactions).Black circular markers with the grey hatched band represent the unfolded data measurement and its uncertainty.Overlaid are the predictions together with their statistical uncertainties (Theory stat.).The BLACK-HAT+SHERPA uncertainty also contains theoretical systematic uncertainties (Theory syst.)described in Section 8.The lower plots show the ratio of each prediction to the unfolded data.

Figure 8 :
Figure8: The ratio of the cross section measurements, and of the predictions of BLACK-HAT+SHERPA using the MSTW2008nlo68cl and NNPDF PDF sets, to the prediction of BLACK-HAT+SHERPA using the CT10 PDF set, for four different quantities.Black circular markers represent the unfolded data measurement and its uncertainty.The grey hatched band indicates the PDF uncertainty on the BLACKHAT+SHERPA predictions with CT10.Triangular (square) markers indicate the ratio of the prediction MSTW2008nlo68cl (NNPDF) to the CT10 prediction, including the PDF uncertainty for the given PDF set.

Table 2 :
Cross section measurements with statistical and systematic uncertainties for inclusive and exclusive jet multiplicities up to 6 jets.FWO (Belgium); CNPq, CAPES, FAPERJ, and FAPESP (Brazil); MES (Bulgaria); CERN; CAS, MoST, and NSFC (China); COLCIENCIAS (Colombia); MSES and CSF (Croatia); RPF (Cyprus); MoER, ERC IUT and ERDF (Estonia); Academy of Finland, MEC, and HIP (Finland); CEA and CNRS/IN2P3 (France); BMBF, DFG, and HGF (Germany); GSRT (Greece); OTKA and NIH (Hungary); DAE and DST (India); IPM (Iran); SFI (Ireland); INFN (Italy); NRF and WCU (Republic of Korea); LAS (Lithuania); MOE and UM (Malaysia); CINVESTAV, CONACYT, SEP, Individuals have received support from the Marie-Curie programme and the European Research Council and EPLANET (European Union); the Leventis Foundation; the A. P. Sloan Foundation; the Alexander von Humboldt Foundation; the Belgian Federal Science Policy Office; the Fonds pour la Formation à la Recherche dans l'Industrie et dans l'Agriculture (FRIA-Belgium); the Agentschap voor Innovatie door Wetenschap en Technologie (IWT-Belgium); the Ministry of Education, Youth and Sports (MEYS) of the Czech Republic; the Council of Science and Industrial Research, India; the HOMING PLUS programme of Foundation for Polish Science, cofinanced from European Union, Regional Development Fund; the Compagnia di San Paolo (Torino); the Thalis and Aristeia programmes cofinanced by EU-ESF and the Greek NSRF; and the National Priorities Research Program by Qatar National Research Fund.