Measurement of associated production of a W boson and a charm quark in proton-proton collisions at $\sqrt{s} =$ 13 TeV

Measurements are presented of associated production of a W boson and a charm quark (W+c) in proton-proton collisions at a center-of-mass energy of 13 TeV. The data correspond to an integrated luminosity of 35.7 fb$^{-1}$ collected by the CMS experiment at the CERN LHC. The W bosons are identified by their decay into a muon and a neutrino. The charm quarks are tagged via the full reconstruction of D$^*$(2010)$^\pm$ mesons that decay via D$^*$(2010)$^\pm \to$ D$^0$ + $\pi^\pm \to$ K$^{\mp}$ + $\pi^\pm$ + $\pi^\pm$. A cross section is measured in the fiducial region defined by the muon transverse momentum $p_{T}^{\mu}>$ 26 GeV, muon pseudorapidity $|\eta^{\mu}|<$ 2.4, and charm quark transverse momentum $p_{T}^{c}>$ 5 GeV. The inclusive cross section for this kinematic range is $\sigma$(W+c) = 1026 $\pm$ 31 (stat) $\substack{+76\\-72}$ (syst) pb. The cross section is also measured differentially as a function of the pseudorapidity of the muon from the W boson decay. These measurements are compared with theoretical predictions and are used to probe the strange quark content of the proton.


Introduction
Precise knowledge of the structure of the proton, expressed in terms of parton distribution functions (PDFs), is important for interpreting results obtained in proton-proton (pp) collisions at the CERN LHC. The PDFs are determined by comparing theoretical predictions obtained at a particular order in perturbative quantum chromodynamics (pQCD) to experimental measurements. The precision of the PDFs, which affects the accuracy of the theoretical predictions for cross sections at the LHC, is determined by the uncertainties of the experimental measurements used, and by the limitations of the available theoretical calculations. The flavor composition of the light quark sea in the proton and, in particular, the understanding of the strange quark distribution is important for the measurement of the W boson mass at the LHC [1]. Therefore, it is of great interest to determine the strange quark distribution with improved precision.
Before the start of LHC data taking, information on the strange quark content of the nucleon was obtained primarily from charm production in (anti)neutrino-iron deep inelastic scattering (DIS) by the NuTeV [2], CCFR [3], and NOMAD [4] experiments. In addition, a direct measurement of inclusive charm production in nuclear emulsions was performed by the CHORUS experiment [5]. At the LHC, the production of W or Z bosons, inclusive or associated with charm quarks, provides an important input for tests of the earlier determinations of the strange quark distribution. The measurements of inclusive W or Z boson production at the LHC, which are indirectly sensitive to the strange quark distribution, were used in a QCD analysis by the AT-LAS experiment, and an enhancement of the strange quark distribution with respect to other measurements was observed [6].
The associated production of W bosons and charm quarks in pp collisions at the LHC probes the strange quark content of the proton directly through the leading order (LO) processes g + s → W + +c and g + s → W − +c, as shown in Fig. 1. The contribution of the Cabibbosuppressed processes g + d → W + +c and g + d → W − +c amounts to only a few percent of the total cross section. Therefore, measurements of associated W+c production in pp collisions provide valuable insights into the strange quark distribution of the proton. Furthermore, these measurements allow important cross-checks of the results obtained in the global PDF fits using the DIS data and measurements of inclusive W and Z boson production at the LHC.
Production of W+c in hadron collisions was first investigated at the Tevatron [7][8][9]. The first measurement of the cross section of W+c production in pp collisions at the LHC was performed by the CMS Collaboration at a center-of-mass energy of √ s = 7 TeV with an integrated luminosity of 5 fb −1 [10]. This measurement was used for the first direct determination of the strange quark distribution in the proton at a hadron collider [11]. The extracted strangeness suppression with respect to u and d quark densities was found to be in agreement with measurements in neutrino scattering experiments. The cross section for W+c production was also measured by the ATLAS experiment at √ s = 7 TeV [12] and used in a QCD analysis, which supported the enhanced strange quark content in the proton suggested by the earlier ATLAS analysis [6]. A subsequent joint QCD analysis [13] of all available data that were sensitive to the strange quark distribution demonstrated consistency between the W+c measurements by the ATLAS and CMS Collaborations. In Ref. [13], possible reasons for the observed strangeness enhancement were discussed. Recent results of an ATLAS QCD analysis [14], including measurements of inclusive W and Z boson production at √ s = 7 TeV, indicated an even stronger strangeness enhancement in disagreement with all global PDFs. In Ref. [15], possible reasons for this observation were attributed to the limitations of the parameterization used in this AT-LAS analysis [14]. The associated production of a W boson with a jet originating from a charm quark is also studied in the forward region by the LHCb experiment [16].
In this paper, the cross section for W+c production is measured in pp collisions at the LHC at √ s = 13 TeV using data collected by the CMS experiment in 2016 corresponding to an integrated luminosity of 35.7 fb −1 . The W bosons are selected via their decay into a muon and a neutrino. The charm quarks are tagged by a full reconstruction of the charmed hadrons in the process c → D * (2010) ± → D 0 + π ± slow → K ∓ + π ± + π ± slow , which has a clear experimental signature. The pion originating from the D * (2010) ± decay receives very little energy because of the small mass difference between D * (2010) ± and D 0 (1865) and is therefore denoted a "slow" pion π ± slow . Cross sections for W+D * (2010) ± production are measured within a selected fiducial phase space. The W+c cross sections are compared with theoretical predictions at next-to-leading order (NLO) QCD, which are obtained with MCFM 6.8 [17][18][19], and are used to extract the strange quark content of the proton. This paper is organized as follows. The CMS detector is briefly described in Section 2. The data and the simulated samples are described in Section 3. The event selection is presented in Section 4. The measurement of the cross sections and the evaluation of systematic uncertainties are discussed in Section 5. The details of the QCD analysis are described in Section 6. Section 7 summarizes the results.

The CMS detector
The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diameter, providing a magnetic field of 3.8 T. Within the solenoid volume are a silicon pixel and strip tracker, a lead tungstate crystal electromagnetic calorimeter (ECAL), and a brass and scintillator hadron calorimeter (HCAL), each composed of a barrel and two endcap sections. Forward calorimeters extend the pseudorapidity coverage provided by the barrel and endcap detectors. Muons are detected in gas-ionization chambers embedded in the steel flux-return yoke outside the solenoid.
The silicon tracker measures charged particles within the pseudorapidity range |η| < 2.5. It consists of 1440 silicon pixel and 15 148 silicon strip detector modules. For nonisolated particles of 1 < p T < 10 GeV and |η| < 1.4, the track resolutions are typically 1.5% in p T and 25-90  µm in the transverse (longitudinal) impact parameter [20]. The reconstructed vertex with the largest value of summed physics-object p 2 T is taken to be the primary pp interaction vertex. The physics objects are the jets, clustered using the jet finding algorithm [21,22] with the tracks assigned to the vertex as inputs, and the associated missing transverse momentum, taken as the negative vector sum of the p T of those jets. Muons are measured in the pseudorapidity range |η| < 2.4, with detection planes made using three technologies: drift tubes, cathode strip chambers, and resistive-plate chambers. The single muon trigger efficiency exceeds 90% over the full η range, and the efficiency to reconstruct and identify muons is greater than 96%. Matching muons to tracks measured in the silicon tracker results in a relative trans-verse momentum resolution, for muons with p T up to 100 GeV, of 1% in the barrel and 3% in the endcaps. The p T resolution in the barrel is better than 7% for muons with p T up to 1 TeV [23]. A more detailed description of the CMS detector, together with a definition of the coordinate system used and the relevant kinematic variables, can be found in Ref. [24].

Data and Monte Carlo samples and signal definition
Candidate events for the muon decay channel of the W boson are selected by a muon trigger [25] that requires a reconstructed muon with p µ T > 24 GeV. The presence of a high-p T neutrino is implied by the missing transverse momentum, p miss T , which is defined as the negative vector sum of the transverse momenta of the reconstructed particles.
Muon candidates and p miss T are reconstructed using the particle-flow (PF) algorithm [26], which reconstructs and identifies each individual particle with an optimized combination of information from the various elements of the CMS detector. The energy of photons is obtained directly from the ECAL measurement. The energy of electrons is determined from a combination of the electron momentum at the primary interaction vertex determined by the tracking detector, the energy of the corresponding ECAL cluster, and the energy sum of all bremsstrahlung photons spatially compatible with originating from the electron track. The muon momentum is obtained from the track curvature in both the tracker and the muon system, and identified by hits in multiple stations of the flux-return yoke. The energy of charged hadrons is determined from a combination of their momentum measured in the tracker and the matching ECAL and HCAL energy deposits, corrected for both zero-suppression effects and the response function of the calorimeters to hadronic showers. Finally, the energy of neutral hadrons is obtained from the corresponding corrected ECAL and HCAL energy.
The D * (2010) ± meson candidates are reconstructed from tracks formed by combining the measurements in the silicon pixel and strip detectors through the CMS combinatorial track finder [20].
The signal and background processes are simulated using Monte Carlo (MC) generators to estimate the acceptance and efficiency of the CMS detector. The corresponding MC events are passed through a detailed GEANT4 [27] simulation of the CMS detector and reconstructed using the same software as the real data. The presence of multiple pp interactions in the same or adjacent bunch crossing (pileup) is incorporated by simulating additional interactions (both in-time and out-of-time with respect to the hard interaction) with a vertex multiplicity that matches the distribution observed in data. The simulated samples are normalized to the integrated luminosity of the data using the generated cross sections. To simulate the signal, inclusive W+jets events are generated with MADGRAPH5 aMC@NLO (v2.2.2) [28] using the NLO matrix elements, interfaced with PYTHIA8 (8.2.12) [29] for parton showering and hadronization. A matching scale of 10 GeV is chosen, and the FXFX technique [30] is applied for matching and merging. The factorization and renormalization scales, µ 2 r and µ 2 f , are set to µ 2 r = µ 2 f = m 2 W + p 2 T,W . The proton structure is described by the NNPDF3.0nlo [31] PDF set. To enrich the sample with simulated W+c events, an event filter that requires at least one muon with p µ T > 20 GeV and |η µ | < 2.4, as well as at least one D * (2010) ± meson, is applied at the generator level.
Several background contributions are considered, which are described in the following. An inclusive W+jets event sample is generated using the same settings as the signal events, but without the event filter, to simulate background contributions from W events that do not contain D * (2010) ± mesons. Events originating from Drell-Yan (DY) with associated jets are simu-lated with MADGRAPH5 aMC@NLO (v2.2.2) with µ 2 r and µ 2 f set to m 2 Z + p 2 T,Z . Events originating from top quark-antiquark pair (tt) production are simulated using POWHEG (v2.0) [32], whereas single top quark events are simulated using POWHEG (v2.0) [33,34] or POWHEG (v1.0) [35], depending on the production channel. Inclusive production of WW, WZ, and ZZ bosons and contributions from the inclusive QCD events are generated using PYTHIA8. The CUETP8M1 [36] underlying event tune is used in PYTHIA8 for all, except for the tt sample, where the CUETP8M2T4 [37] tune is applied.
The dominant background originates from processes like u + d → W + + g * → W + + cc or d + u → W − + g * → W − + cc, with c quarks produced in gluon splitting. In the W+c signal events the charges of the W boson and the charm quark have opposite signs. In gluon splitting, an additional c quark is produced with the same charge as the W boson. At the generator level, an event is considered as a W+c event if it contains at least one charm quark in the final state. In the case of an odd number of c quarks, the c quark with the highest p T and a charge opposite to that of the W boson is considered as originating from a W+c process, whereas the other c quarks in the event are labeled as originating from gluon splitting. In the case of an even number of c quarks, all are considered to come from gluon splitting. Events containing both c and b quarks are considered to be W+c events, since c quarks are of higher priority in this analysis, regardless of their momentum or production mechanism. Events containing no c quark and at least one b quark are classified as W + b. Otherwise, an event is assigned to the W + udsg category.
The contribution from gluon splitting can be significantly reduced using data. Events with the same charge sign for both the W boson and charm quark, which correlates to the charge sign of the D * (2010) ± meson, are background, which is due to gluon splitting. Since the gluon splitting background for opposite charge pairs is identical, it can be removed by subtracting the same-sign distribution from the signal. The measurement is performed in the central kinematic range and is not sensitive to the contributions of processes c + g → W + s with a spectator charm quark.
For validation and tuning of MC event generators using a Rivet plugin [38], the W+D * (2010) ± measurement is performed. This requires a particle-level definition without constraints on the origin of D * (2010) ± mesons. Therefore, any contributions from B meson decays and other hadrons, though only a few pb, are included as signal for this part of the measurement.

Event selection
The associated production of W bosons and charm quarks is investigated using events, where W → µ + ν µ and the c quarks hadronize into a D * (2010) ± meson. The reconstruction of the muons from the W boson decays and of the D * (2010) ± candidates is described in detail in the following.

Selection of W boson candidates
Events containing a W boson decay are identified by the presence of a high-p T isolated muon and p miss T . The muon candidates are reconstructed by combining the tracking information from the muon system and from the inner tracking system [23], using the CMS particle-flow algorithm. Muon candidates are required to have p µ T > 26 GeV, |η µ | < 2.4, and must fulfill the CMS "tight identification" criteria [23]. To suppress contamination from muons contained in jets, an isolation requirement is imposed: where the p T sum of PF candidates for charged hadrons (CH), neutral hadrons (NH), photons (EM) and charged particles from pileup (PU) inside a cone of radius ∆R ≤ 0.4 is used, and the factor 0.5 corresponds to the typical ratio of neutral to charged particles, as measured in jet production [26].
Events in which more than one muon candidate fulfills all the above criteria are rejected in order to suppress background from DY processes. Corrections are applied to the simulated samples to adjust the trigger, isolation, identification, and tracking efficiencies to the observed data. These correction factors are determined through dedicated tag-and-probe studies.
The presence of a neutrino in an event is assured by imposing a requirement on the transverse mass, which is defined as the combination of p µ T and p miss T : In this analysis, m T > 50 GeV is required, which results in a significant reduction of background.
The D 0 candidates are constructed by combining two oppositely charged tracks with transverse momenta p track T > 1 GeV, assuming the K ∓ and π ± masses. The D 0 candidates are further combined with a track of opposite charge to the kaon candidate, assuming the π ± mass, following the well-established procedure of Ref. [40,41]. The invariant mass of the K ∓ π ± combination is required to be |m( [39]. The candidate K ∓ and π ± tracks must originate at a fitted secondary vertex [42] that is displaced by not more than 0.1 cm in both the xy-plane and z-coordinate from the third track, which is presumed to be the π ± slow candidate. The latter is required to have p track T > 0.35 GeV and to be in a cone of ∆R ≤ 0.15 around the direction of the D 0 candidate momentum. The combinatorial background is reduced by requiring the D * (2010) ± transverse momentum p D * T > 5 GeV and by applying an isolation criterion p D * T / ∑ p T > 0.2. Here ∑ p T is the sum of transverse momenta of tracks in a cone of ∆R ≤ 0.4 around the direction of the D * (2010) ± momentum. The contribution of D * (2010) ± mesons produced in pileup events is suppressed by rejecting candidates with a z-distance > 0.2 cm between the muon and the π ± slow . After applying all selection criteria, the contribution of D 0 decays other than K ∓ π ± is negligible compared to the uncertainties.
The D * (2010) ± meson candidates are identified using the mass difference method [41] via a peak in the ∆m(D * , D 0 ) distribution. Wrong-charge combinations with K ± π ± pairs in the accepted D 0 mass range mimic the background originating from light-flavor hadrons. By subtracting the wrong-charge combinations, the combinatorial background in the ∆m(D * , D 0 ) distribution is mostly removed. The presence of nonresonant charm production in the rightcharge K ∓ π ± π ± combinations introduces a small normalization difference of ∆m(D * , D 0 ) distributions for right-and wrong-charge combinations, which is corrected utilizing fits to the ratio of both distributions.

Selection of W+c candidates
An event is selected as a W+c signal if it contains a W boson and a D * (2010) ± candidate fulfilling all selection criteria. The candidate events are split into two categories: with W ± + D * (2010) ± combinations falling into the same sign (SS) category, and W ∓ + D * (2010) ± combinations falling into the opposite sign (OS) category. The signal events consist of only OS combinations, whereas the W + cc and W + bb background processes produce the same number of OS and SS candidates. Therefore, subtracting the SS events from the OS events removes the background contributions from gluon splitting. The contributions from other background sources, such as tt and single top quark production, are negligible.
The number of W+c events corresponds to the number of D * (2010) ± mesons after the subtraction of light-flavor and gluon splitting backgrounds. The invariant mass of K ∓ π ± candidates, which are selected in a ∆m(D * , D 0 ) window of ±1 MeV, is shown in Fig. 2, along with the observed reconstructed mass difference ∆m(D * , D 0 ). A clear D 0 peak at the expected mass and a clear ∆m(D * , D 0 ) peak around the expected value of 145.4257 ± 0.0017 MeV [39] are observed. The remaining background is negligible, and the number of D * (2010) ± mesons is determined by counting the number of candidates in a window of 144 < ∆m(D * , D 0 ) < 147 MeV. Alternately, two different functions are fit to the distributions, and their integral over the same mass window is used to estimate the systematic uncertainties associated with the method chosen.

Measurement of the fiducial W+c cross section
The fiducial cross section is measured in a kinematic region defined by requirements on the transverse momentum and the pseudorapidity of the muon and the transverse momentum of the charm quark. The simulated signal is used to extrapolate from the fiducial region of the D * (2010) ± meson to the fiducial region of the charm quark. Since the D * (2010) ± kinematics is integrated over at the generator level, the only kinematic constraint on the corresponding charm quark arises from the requirement on the transverse momentum of D * (2010) ± meson. The correlation of the kinematics of charm quarks and D * (2010) ± mesons is investigated using simulation, and the requirement of p D * T > 5 GeV translates into p c T > 5 GeV. The distributions of |η µ | and p c T in the simulation are shown to reproduce very well the fixed order prediction at NLO obtained, using MCFM 6.8 [17][18][19] calculation. The kinematic range of the measured fiducial cross section corresponds to p µ T > 26 GeV, |η µ | < 2.4, and p c T > 5 GeV. The fiducial W+c cross section is determined as: where N sel is the number of selected OS − SS events in the ∆m(D * , D 0 ) distribution and S is the signal fraction. The latter is defined as the ratio of the number of reconstructed W+D * (2010) ± candidates originating from W+c to the number of all reconstructed D * (2010) ± . It is determined from the MC simulation, includes the background contributions, and varies between 0.95 and 0.99. The integrated luminosity is denoted by L int . The combined branching fraction B for the channels under study is a product of B(c → D * (2010) ± ) = 0.2429 ± 0.0049 [43] and B(D * (2010) ± → K ∓ + π ± + π ± slow ) = 0.0266 ± 0.0003 [39]. The correction factor C accounts for the acceptance and efficiency of the detector. The latter is determined using the MC simulation and is defined as the ratio of the number of reconstructed W+D * (2010) ± candidates to the number of generated W+D * (2010) ± originating from W+c events that fulfill the fiducial requirements. In the measurement of the W + +c (W + +D * (2010) − ) and W − +c (W − +D * (2010) + ) cross sections, the factor C is determined separately for different charge combinations.
The measurement of the W+c cross section relies to a large extent on the MC simulation and requires extrapolation to unmeasured phase space. To reduce the extrapolation and the corresponding uncertainty, the cross section for W+D * (2010) ± production is also determined in the fiducial phase space of the detector-level measurement, p µ T > 26 GeV, |η µ | < 2.4, |η D * | < 2.4 and p D * T > 5 GeV, in a similar way by modifying Eq. (3) as follows: only the branching fraction B = B(D * (2010) ± → K ∓ + π ± + π ± slow ) is considered and the factor C is defined as the ratio between the numbers of reconstructed and of generated W+D * (2010) ± candidates in the fiducial phase space after OS − SS subtraction.
The cross sections are determined inclusively and also in five bins of the absolute pseudorapidity |η µ | of the muon originating from the W boson decay. The number of signal (OS − SS) events in each range of |η µ | is shown in Fig. 3. Good agreement between the data and MC simulation within the statistical uncertainties is observed.

Systematic uncertainties
The efficiencies and the assumptions relevant for the measurement are varied within their uncertainties to estimate the systematic uncertainty in the cross section measurement. The resulting shift of the cross section with respect to the central result is taken as the corresponding uncertainty contribution. The various sources of the systematic uncertainties in the W+c production cross section are listed in Table 1 for both the inclusive and the differential measurements.
• Uncertainties associated with the integrated luminosity measurement are estimated as 2.5% [44].
• The uncertainty in the tracking efficiency is 2.3% for the 2016 data. It is determined using the same method described in Ref. [45], which exploits the ratio of branching fractions between the four-body and two-body decays of the neutral charm meson to all-charged final states.
• The muon systematic uncertainties are 1% each for for the muon identification and isolation, and 0.5% for the trigger and tracking corrections. These are added in quadrature and the resulting uncertainty for muons is 1.2%, which is referred to as the 'muon uncertainty'.
• The uncertainty in the determination of N sel is estimated from the difference in using a Gaussian or Crystal Ball fit [46]. The largest value of this uncertainty determined differentially, 1.5%, is considered for all.
• Uncertainties in the modeling of kinematic observables of the generated D * (2010) ± meson are estimated by reweighting the simulated p D * T and |η D * | distributions to the shape observed in data. The respective uncertainty in the inclusive cross section measurement is 0.5%. Due to statistical limitations, this uncertainty is determined inclusively in |η µ |. • The uncertainty in the difference of the normalization of the ∆m(D * , D 0 ) distributions for K ∓ π ± π ± and K ± π ± π ∓ combinations ('background normalization') is 0.5%.
• Uncertainties in the measured p miss T are estimated in Ref. [47] and result in an overall uncertainty of 0.9% for this analysis.
• Uncertainties due to the modeling of pileup are estimated by varying the total inelastic cross section used in the simulation of pileup events by 5%. The corresponding uncertainty in the W+c cross section is 2%.
• The uncertainty related to the requirement of a valid secondary vertex, fitted from the tracks associated with a D 0 candidate, is determined by calculating the D * (2010) ± reconstruction efficiency in data and MC simulation for events with and without applying this selection criterion. The number of reconstructed D * (2010) ± candidates after the SS event subtraction is compared for events with or without a valid secondary vertex along with the proximity requirement (∆ xy < 0.1 cm, ∆ z < 0.1 cm). The difference in efficiency between data and MC simulation is calculated and an uncertainty in the inclusive cross section of −1.1% is obtained. Since this variation is not symmetric, the uncertainty is one-sided.
• The PDF uncertainties are determined according to the prescription of the PDF group  Figure 3: Number of events after OS − SS subtraction for data (filled circles) and MC simulation (filled histograms) as a function of |η µ |. [31]. These are added in quadrature to the uncertainty related to the variation of α S (m Z ) in the PDF, resulting in an uncertainty of 1.2% in the inclusive cross section. • The uncertainty associated with the fragmentation of the c quark into a D * (2010) ± meson is determined through variations of the function describing the fragmentation parameter z = p D * T /p c T . The investigation of this uncertainty is inspired by a dedicated measurement of the c → D * (2010) ± fragmentation function in electronproton collisions [48], in which the fragmentation parameters in various phenomenological models were determined with an uncertainty of 10%. In the PYTHIA MC event generator, the fragmentation is described by the phenomenological Bowler-Lund function [49,50], in the form To estimate the uncertainty, the parameters a and b are varied within ±10% around their central values, following the precision achieved for the fragmentation parameters in [48]. An additional constraint on the upper boundary on the a parameter in PYTHIA is consistent with this 10% variation. The resulting uncertainty in the cross section is 3.9%.

Cross section results
The numbers of signal events and the inclusive fiducial cross sections with their uncertainties are listed in Table 2 together with the ratio of σ(W + +c)/σ(W − +c). For the differential measurement of the W+c cross section, the numbers of signal events are summarized in Table 3 together with the corrections C derived using MC simulations in each |η µ | bin. The results are presented for dσ(W+c)/d|η µ |, as well as for dσ(W + +c)/d|η µ | and for dσ(W − +c)/d|η µ |.
The measured inclusive and differential fiducial cross sections of W+c are compared to predictions at NLO (O(α 2 s )) that are obtained using MCFM 6.8. Similarly to the earlier analysis [11], the mass of the charm quark is chosen to be m c = 1.5 GeV, and the factorization and the renormalization scales are set to the value of the W boson mass. The calculation is performed for p µ T > 26 GeV, |η µ | < 2.4, and p c T > 5 GeV. In Fig. 4, the measurements of the inclusive W+c cross section and the charge ratio are compared to the NLO predictions calculated using the ABMP16nlo [51], ATLASepWZ16nnlo [14], CT14nlo [52], MMHT14nlo [53], NNPDF3.0nlo [31], and NNPDF3.1nlo [54] PDF sets. The values of the strong coupling constant α S (m Z ) are set to those used in the evaluation of a particular PDF. The details of the experimental data, used for constraining the strange quark content of the proton in the global PDFs, are given in Refs. [14,31,52,53,55]. In these references, the treatment of the sea quark distributions in different PDF sets is discussed, and the comparison of the PDFs is presented. The ABMP16nlo PDF includes the most recent data on charm quark production in chargedcurrent neutrino-nucleon DIS collected by the NOMAD and CHORUS experiments in order to improve the constraints on the strange quark distribution and to perform a detailed study of the isospin asymmetry of the light quarks in the proton sea [56]. Despite differences in the data used in the individual global PDF fits, the strangeness suppression distributions in ABMP16nlo, NNPDF3.1nlo, CT14nlo and MMHT14nlo are in a good agreement among each other and disagree with the ATLASepWZ16nnlo result [14]. The predicted inclusive cross sections are summarized in Table 4. The PDF uncertainties are calculated using prescriptions from each PDF group. For the ATLASepWZ16nnlo PDFs no respective NLO set is available and only Hessian uncertainties are considered in this paper. For other PDFs, the variation of α s (m Z ) is taken into account as well. The uncertainties due to missing higher-order corrections are estimated by varying µ r and µ f simultaneously by a factor of 2 up and down, and the resulting variation of the cross section is referred to as the scale uncertainty, ∆µ. Good agreement between NLO predictions and the measurements is observed, except for the prediction using ATLASepWZ16nnlo. For the cross section ratio σ(W + +c)/σ(W − +c), all theoretical predictions are in good agreement with the measured value. In Table 5, the theoretical predictions for dσ(W+c)/d|η µ | using different PDF sets are summarized. In Fig. 5, the measurements of differential W+c and W+D * (2010) ± cross sections are compared with the MCFM NLO calculations and with the signal MC prediction, respectively. Good agreement between the measured W+c cross section and NLO calculations is observed except for the prediction using the ATLASepWZ16nnlo PDF set. The signal MC prediction using NNPDF3.0nlo is presented with the PDF and α s uncertainties and accounts for simultaneous variations of µ r and µ f in the matrix element by a factor of 2. The W+D * (2010) ± cross section is described well by the simulation.

Impact on the strange quark distribution in the proton
The associated W+c production at 13 TeV probes the strange quark distribution directly in the kinematic range of x ≈ 0.007 at the scale of m 2 W . The first measurement of a fiducial W+c cross section in pp collisions was performed by the CMS experiment at a center-of-mass energy √ s = 7 TeV with a total integrated luminosity of 5 fb −1 [11]. The results were used in a

Figure 4:
Inclusive fiducial cross section σ(W+c) and the cross section ratio σ(W + +c)/σ(W − +c) at 13 TeV. The data are represented by a line with the statistical (total) uncertainty shown by a light (dark) shaded band. The measurements are compared to the NLO QCD prediction using several PDF sets, represented by symbols of different types. All used PDF sets are evaluated at NLO, except for ATLASepWZ16, which is obtained at NNLO. The error bars depict the total theoretical uncertainty, including the PDF and the scale variation uncertainty.  Figure 5: Left: differential cross sections of σ(W+c) production at 13 TeV measured as a function |η µ |. The data are presented by filled circles with the statistical (total) uncertainties shown by vertical error bars (light shaded bands). The measurements are compared to the QCD predictions calculated with MCFM at NLO using different PDF sets, presented by symbols of different style. All used PDF sets are evaluated at NLO, except for ATLASepWZ16, which is obtained at NNLO. The error bars represent theoretical uncertainties, which include PDF and scale variation uncertainty. Right: σ(W+D * (2010) ± ) production differential cross sections presented as a function of |η µ |. The data (filled circles) are shown with their total (outer error bars) and statistical (inner error bars) uncertainties and are compared to the predictions of the signal MC generated with MADGRAPH5 aMC@NLO and using NNPDF3.0nlo to describe the proton structure. PDF uncertainties and scale variations are accounted for and added in quadrature (shaded band). QCD analysis [10] together with measurements of neutral-and charged-current cross sections of DIS at HERA [57] and of the lepton charge asymmetry in W production at √ s = 7 TeV at the LHC [11].
The present measurement of the W+c production cross section at 13 TeV, determined as a function of the absolute pseudorapidity |η µ | of the muon from the W boson decay and p µ T > 26 GeV, is used in an NLO QCD analysis. This analysis also includes an updated combination of the inclusive DIS cross sections [58] and the available CMS measurements of the lepton charge asymmetry in W boson production at √ s = 7 TeV [11] and at √ s = 8 TeV [59]. These latter measurements probe the valence quark distributions in the kinematic range 10 −3 ≤ x ≤ 10 −1 and have indirect sensitivity to the strange quark distribution. The earlier CMS measurement [10] of W+c production at √ s = 7 TeV is also used to exploit the strange quark sensitive measurements at CMS in a joint QCD analysis. The correlations of the experimental uncertainties within each individual data set are taken into account, whereas the CMS data sets are treated as uncorrelated to each other. In particular, the measurements of W+c production at 7 and 13 TeV are treated as uncorrelated because of the different methods of charm tagging and the differences in reconstruction and event selection in the two data sets.
The Thorne-Roberts [70,71] general mass variable flavor number scheme at NLO is used for the treatment of heavy quark contributions with heavy quark masses m b = 4.5 GeV and m c = 1.5 GeV, which correspond to the values used in the signal MC simulation in the cross section measurements. The renormalization and factorization scales are set to Q, which denotes the four-momentum transfer for the case of the DIS data and the mass of the W boson for the case of the muon charge asymmetry and the W+c measurement. The strong coupling constant is set to α s (m Z ) = 0.118. The Q 2 range of HERA data is restricted to Q 2 ≥ Q 2 min = 3.5 GeV 2 to ensure the applicability of pQCD over the kinematic range of the fit. The procedure for the determination of the PDFs follows the approach used in the earlier CMS analyses [11,59]. In the following, a similar PDF parameterization is used as in the most recent CMS QCD analysis [59] of inclusive W boson production.
The parameterized PDFs are the gluon distribution, xg, the valence quark distributions, xu v , xd v , the u-type, xu, and xd-type anti-quark distributions, with xs (xs) denoting the strange (anti-)quark distribution. The initial scale of the QCD evolution is chosen as Q 2 0 = 1.9 GeV 2 . At this scale, the parton distributions are parameterized as: The normalization parameters A u v , A d v , A g are determined by the QCD sum rules, the B parameter is responsible for small-x behavior of the PDFs, and the parameter C describes the shape of the distribution as x → 1. The strangeness fraction f s = s/(d + s) is a free parameter in the fit.
The strange quark distribution is determined by fitting the free parameters in Eqs. (4)- (9). The constraint A u = A d ensures the same normalization for u and d densities at x → 0. It is assumed that xs = xs.
In the earlier CMS analysis [11], the assumption B u = B d was applied. An alternative assumption B u = B d led to a significant change in the result, which was included in the parameterization uncertainty. In the present analysis, the B parameters of the light sea quarks are independent from each other, B u = B d = B s , following the suggestion of Ref. [15].
For all measured data, the predicted and measured cross sections together with their corresponding uncertainties are used to build a global χ 2 , minimized to determine the initial PDF parameters [61,62]. The quality of the overall fit can be judged based on the global χ 2 divided by the number of degrees of freedom, n dof . For each data set included in the fit, a partial χ 2 divided by the number of measurements (data points), n dp , is provided. The correlated part of χ 2 quantifies the influence of the correlated systematic uncertainties in the fit. The global and partial χ 2 values for each data set are listed in Table 6, illustrating a general agreement among all the data sets. Total The PDF uncertainties are investigated according to the general approach of HERAPDF 1.0 [57]. The experimental PDF uncertainties arising from the uncertainties in the measurements are estimated by using the Hessian method [72], adopting the tolerance criterion of ∆χ 2 = 1. The experimental uncertainties correspond to 68% confidence level. Alternatively, the experimental uncertainties in the measurements are propagated to the extracted QCD fit parameters using the MC method [73,74]. In this method, 426 replicas of pseudodata are generated, with measured values for the cross sections allowed to vary within the statistical and systematic uncertainties. For each of them, the PDF fit is performed and the uncertainty is estimated as the root-mean-square around the central value. Because of possible nonGaussian tails in the PDF uncertainties, the MC method is usually considered to be more robust and to give more realistic uncertainties, in particular for PDFs not strongly constrained by the measurements, e.g., in the case of too little or not very precise data. In Fig. 6, the distributions of the strange quark content s(x, Q 2 ), and of the strangeness suppression factor r s (x, µ 2 f ) = (s + s)/(u + d) are presented. In Fig. 7 the strangeness suppression factor is shown in comparison with the ATLASepWZ16nnlo and the ABMP16nlo, similar to Fig. 1 in Ref. [15]. Whereas the CMS result for r s (x) is close to the ABMP16nlo PDF, it shows a significant difference with regard to the ATLASepWZ16nnlo PDF for x > 10 −3 . The significant excess of the strange quark content in the proton reported by ATLAS [14] is not observed in the present analysis.
To investigate the impact of model assumptions on the resulting PDFs, alternative fits are performed, in which the heavy quark masses are varied as 4.3 ≤ m b ≤ 5.0 GeV, 1.37 ≤ m c ≤ 1.55 GeV, and the value of Q 2 min imposed on the HERA data is varied in the interval 2.5 ≤ Q 2 min ≤ 5.0 GeV 2 . Also, the variations in PDF parameterization, following Ref. [59] are investigated. These variations do not alter the results for the strange quark distribution or the suppression factor significantly, compared to the PDF fit uncertainty. Since each global PDF group is using their own assumptions for the values of heavy quark masses and cutoffs on the DIS data, these model variations are not quantified further. To compare the results of the present PDF fit with the earlier determination of the strange quark content in the proton at CMS [11], the "free-s" parameterization of Ref.
[11] is used. There, a flexible form [70,71] for the gluon distribution was adopted, allowing the gluon to be negative. Furthermore, the condition B u = B d = B s was applied in the central parameterization, while B d = B s was used to estimate the parameterization uncertainty. A complete release of the condition B u = B d = B s was not possible because of limited data input, in contrast to the current analysis. The same PDF parameterization was used in the ATLASepWZ16nnlo analysis [14]. The results are presented in Fig. 8. The central value obtained of the s quark distribution is well within the experimental uncertainty of the results at √ s = 7 TeV, while the PDF uncertainty is reduced.

Summary
Associated production of W bosons with charm quarks in proton-proton collisions at √ s = 13 TeV is measured using the data collected by the CMS experiment in 2016 and corresponding to an integrated luminosity of 35.7 fb −1 . The W boson is detected via the presence of a high-p T muon and missing transverse momentum, suggesting the presence of a neutrino. The charm quark is identified via the full reconstruction of the D * (2010) ± meson decaying to D 0 + π ± slow → K ∓ + π ± + π ± slow . Since in W+c production the W boson and the c quark have opposite charge, contributions from background processes, mainly c quark production from gluon splitting, are largely removed by subtracting the events in which the charges of the W boson and of the D * (2010) ± meson have the same sign. The fiducial cross sections are measured in the kinematic range of the muon transverse momentum p µ T > 26 GeV, pseudorapidity |η µ | < 2.4, and transverse momentum of the charm quark p c T > 5 GeV. The fiducial cross section of W+D * (2010) ± production is measured in the kinematic range p µ T > 26 GeV, |η µ | < 2.4, transverse momentum of the D * (2010) ± meson p D * T > 5 GeV and |η D * | < 2.4, and compared to the Monte Carlo prediction. The measurements are performed inclusively and in five bins of |η µ |.
The measurements are in good agreement with the theoretical predictions at next-to-leading order (NLO) for different sets of parton distribution functions (PDF), except for the one using the ATLASepWZ16nnlo PDF. To illustrate the impact of these measurements in the determination of the strange quark distribution in the proton, the data is used in a QCD analysis at NLO together with inclusive DIS measurements and earlier results from CMS on W+c production and the lepton charge asymmetry in W boson production. The strange quark distribution and the strangeness suppression factor r s (x, µ 2 f ) = (s + s)/(u + d) are determined and agree with earlier results obtained in neutrino scattering experiments. The results do not support the hypothesis of an enhanced strange quark contribution in the proton quark sea reported by ATLAS [14].

Acknowledgments
We congratulate our colleagues in the CERN accelerator departments for the excellent performance of the LHC and thank the technical and administrative staffs at CERN and at other CMS institutes for their contributions to the success of the CMS effort. In addition, we gratefully acknowledge the computing centers and personnel of the Worldwide LHC Computing Grid for delivering so effectively the computing infrastructure essential to our analyses. Finally, we acknowledge the enduring support for the construction and operation of the LHC and the CMS detector provided by the following funding agencies: [14] ATLAS Collaboration, "Precision measurement and interpretation of inclusive W + , W − and Z/γ * production cross sections with the ATLAS detector", Eur. Phys. J. C 77 (2017) 367, doi:10.1140/epjc/s10052-017-4911-9, arXiv:1612.03016.