## QCD analysis of CMS W + charm measurements at LHC with √ s = 7 TeV and implications for strange PDF

Nijat Yalkun 1, 2, 3 and Sayipjamal Dulat 1, ∗

1 School of Physics Science and Technology, Xinjiang University, Urumqi, Xinjiang 830046 China CAS Key Laboratory of Theoretical Physics, Institute of Theoretical Physics, Chinese Academy of Sciences, Beijing 100190,China

2

3 University of Chinese Academy of Sciences (UCAS), Beijing 100049, China

## Abstract

We calculate cross-sections and cross-section ratios of a charm quark production in association with a W gauge boson at next-to-leading order QCD using MadGraph and CT10NNLO, CT14NNLO, and MSTW2008NNLO PDFs. We compare the results with measurements from the CMS detector at the LHC at a center-of-mass energy of 7TeV. Moreover, we calculate absolute and normalized differential cross-sections as well as differential cross-section ratios as a function of the lepton pseudorapidity from the W boson decay. The correlation between the CT14NNLO PDFs and predictions for W + charm data are studied as well. Furthermore, by employing the error PDF updating method proposed by the CTEQ-TEA group, we update CT14NNLO PDFs, and analyze the impact of CMS 7TeV W + charm production data to the original CT14NNLO PDFs. By comparison of the g x, Q ( ), s x, Q ( ), u x, Q ( ), d x, Q ( ), ¯( u x, Q ), and ¯ ( d x, Q ) PDFs at Q = 1 3 GeV . and Q = 100 GeV for the CT14NNLO and CT14NNLO+Wc, we see that the error band of the s x, Q ( ) PDF is reduced in the region x &lt; 0 4, and the error band of . g x, Q ( ) PDF is also slightly reduced in the region 0 01 . &lt; x &lt; 0 1. .

PACS numbers: 12.15.Ji, 12.38 Cy, 13.85.Qk

Keywords: parton distribution functions; W + charm production;

∗ Electronic address: sdulat@hotmail.com

## Contents

| I Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .   | . 1   |
|----------------------------------------------------------------------------------------------------------|-------|
| II Results . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . | . 2   |
| A Total cross-section. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .           | . 3   |
| B Absolute and normalized differential cross-section . . . . . . . . . . . . . . . . .                   | . 5   |
| C Charged cross-section ratio . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .              | . 9   |
| D Correlation between the CT14NNLO and predictions for W + charm data .                                  | . 12  |
| III Using EPUMP to study the impact of W + charm data on CT14NNLO PDFs .                                 | . 15  |
| IV Conclusions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .     | . 20  |
| Acknowledgments . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .        | . 21  |
| References . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .   | . 21  |

## I. INTRODUCTION

In the standard model (SM), the associated W + charm production in hadron collisions is described at leading order (LO) in perturbative quantum chromodynamics (QCD) by g + q → W -+ , ( c q = d, s, b ) and g + ¯ q → W + + ¯, c (¯ = q ¯ ¯ ¯ ). d, s, b Although the dquark parton distribution function (PDF) is large in the proton, the processes g + d → W -+ c and g + ¯ d → W + + ¯ contribute c only about 10% [1] to the total W + charm production rate, because it is suppressed by the small quark-mixing Cabibbo-KobayashiMaskawa (CKM) matrix element[2] | V cd | and | V ¯ ¯ cd | . The major contribution to the total W + charm production rate is due to strange quark-gluon fusion g + s → W -+ , and c g +¯ s → W + +¯. The contribution from c g + b → W -+ c and g + ¯ b → W + +¯is also heavily c suppressed by the quark mixing matrix elements ( | V cb | , | V ¯ ¯ cb | ) and the b-quark PDF. The W + charm production cross-section is therefore particularly sensitive to the proton g x, Q ( ) and s x, Q ( ) PDFs [3] and to the magnitude of the CKM matrix element V cs , where x is the momentum fraction of the proton carried by the s-quark, and Q is the hard scale. The study cited in Ref.[4] calculated the W + charm production at LO and next-to-leading order (NLO) in QCD, and found that the factorization and renormalization scale uncertainty in the NLO calculation is about 20%. Ref.[5] explored the strangeness degrees of freedom in the parton structure of the nucleon within the global analysis framework, and showed that the precise determination of the s x, Q ( ) PDF affects the W + charm cross-section. The s x, Q ( ) PDF has been determined by neutrino-nucleon deep inelastic scattering experiments at momentum transfer squared Q 2 = 10 GeV, and momentum fraction x ∼ 0 1 . [6, 7]. The Tevatron CDF[8] and D0[9] experiments have measured the cross-section for charm quark produced in association with W bosons, using muon tagging of the charm-quark jet. The ATLAS collaboration [10] measured the total cross-section, differential cross-section as a function of the pseudorapidity of the lepton from the W boson decay, and the cross-section ratio of the production of a W boson in association with a single charm quark at √ s = 7 TeV. The CMS experiment measured [11] total cross-sections ( σ W ( -+ ) c , σ ( W + + ¯)), absolute and c normalized differential cross-sections as a function of the absolute value of the pseudorapidity of the lepton from the W boson decay, and the cross-section ratio R c = σ W ( + +¯) c /σ W ( -+ ) c

at a center of mass energy 7TeV for the fiducial region defined, namely

$$p _ { T } ^ { j } & > 2 5 \, \text{GeV}, \ | \eta ^ { j } | < 2. 5, \ | \eta ^ { l } | < 2. 1, \\ p _ { T } ^ { l } & > 2 5 \, \text{GeV}, \ \text{for} \ W \to \mu \nu _ { \mu }, \\ p _ { T } ^ { l } & > 3 5 \, \text{GeV}, \ \text{for} \ W \to \mu \nu _ { \mu } \ \text{and} \ W \to e \nu _ { e }.$$

There are two different transverse momentum cuts for the charged lepton in the final state. When p l T &gt; 25 GeV, we only consider the muon decay channel( W → µν µ ) for W boson; both muon( W → µν µ ) and electron( W → eν e ) decay channels for W boson are considered, when p l T &gt; 35 GeV. This study is organized as follows: in Section II, we present our results for various latest PDF sets and compare these with the CMS measurements of the total cross-section, absolute and normalized differential cross-sections and ratios, as well as the correlation between the CT14NNLO PDFs and predictions for W + charm data. In Section III, we discuss the impact of the CMS W + charm production 7TeV data to the CT14NNLO PDFs. In Section IV, we draw our conclusions.

## II. RESULTS

In this section, we present a detailed numerical study of the pp → W + c + X process at the LHC at a center-of-mass energy of 7TeV at NLO order QCD using the Monte-Carlo numerical calculation program MadGraph[12] with CT10NNLO[13], CT14NNLO[14], and MSTW2008NNLO[15] PDFs. PDF uncertainties on the theoretical predictions are given at 68% confidence level (C.L.). We calculate the total cross-section, differential (absolute and normalized) cross-sections, and the cross-section ratio R c = σ W ( + + ¯) c /σ W ( -+ ) with c the W → lν decay (where l = µ or e ). In our study, we use the same kinematical cuts as the CMS detector at the LHC at a center-of-mass energy of 7TeV [11], that are given in section I; both the factorization and renormalization scales are set to the value of the W boson mass µ R = µ F = m w ; charm quark mass is considered and is set to 1.550 GeV; strong interaction coupling α s is set to 0.118, and for electro-weak parameters, the W boson mass is set to 80.385 GeV; Fermi coupling is set to as 1 166 . × 10 -5 GeV -2 ; related CKM matrix elements are set to as V cd = 0 225 and . V cs = 0 974; the mass of charged light leptons is . considered and set to as m e = 0 511 MeV and . m µ = 105 658 MeV. At LO, the Feynman . diagrams for the hard scattering processes of the W + charm production pp → W + c + X

are shown in Fig.1. The main contribution for the cross-sections of W + charm production comes from strange quark and gluon scattering, the down-quark contribution is strongly Cabibbo suppressed, and the contribution from the bottom quark and gluon scattering is negligible.

FIG. 1: Possible tree-level diagrams at partonic level for W + charm production

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00026v2-figure-291743dd-7de5-493b-8b4b-3d1abaa30e9a.jpg)

## A. Total cross-section

The total cross-sections σ W ( + + ¯) and c σ W ( -+ ) of the production of a c W boson in association with a charm quark in pp collisions at √ s = 7 TeV at NLO QCD are summarized in Table I. PDF uncertainties are at 68% confidence level (C.L.), that are obtained from the error sets of the CT10NNLO, CT14NNLO, and MSTW2008NNLO. The experimental measurements from CMS collaboration at the LHC at a center-of-mass energy of 7TeV [11] are also included in this table. The comparison between theory predictions based on various PDF sets and experimental measurements are illustrated in Fig 2. The predictions obtained with the CT10NNLO PDFs are in agreement with CMS measurements. The predictions obtained with CT14NNLO agree with CMS measurements within the uncertainty range. The prediction obtained with MSTW2008NNLO is less favored. Those differences in the size of the PDF uncertainties depend on a different methodology and the parametrization of the s x, Q ( ) PDF used by different PDF sets.

TABLE I: Total cross-section of σ pp ( → W + ) c × B W ( → lν )

|              | σ ( pp → W + c ) × B ( W → lν )[pb]   | σ ( pp → W + c ) × B ( W → lν )[pb]   |
|--------------|---------------------------------------|---------------------------------------|
| PDF sets     | p l T > 25GeV                         | p l T > 35GeV                         |
| CT10NNLO     | 108.1 +6 . 6% - 5 . 5%                | 86.4 +6 . 7% - 5 . 5%                 |
| CT14NNLO     | 100.4 +7 . 1% - 10 . 0%               | 80.1 +7 . 2% - 10 . 1%                |
| MSTW2008NNLO | 98.5 +2 . 1% - 2 . 6%                 | 78.7 +2 . 1% - 2 . 6%                 |
| CMS          | 107.7 ± 3.1% (stat.) ± 6.4% (syst.)   | 84.1 ± 2.4% (stat.) ± 5.8% (syst.)    |

FIG. 2: Comparison of theoretical predictions for total cross-section σ W ( + ) computed with Madc Graph using the CT10NNLO, CT14NNLO and MSTW2008NNLO PDFs with CMS measurements. The left figure contains the prediction for the lepton from W boson decay with p l T &gt; 25GeV, while the right figure is for the case p l T &gt; 35GeV. The solid vertical line shows the central value of the measurement, the inner error band corresponds to the statistical uncertainty, and the outer error band to the sum in quadrature of the statistical and systematic uncertainties.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00026v2-figure-75fc1760-3eea-4102-a62e-1a94856f44df.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00026v2-figure-1e3c905a-e104-4675-8a81-d2856c61aad8.jpg)

Our MadGraph calculations of the d , ¯ , d s , and ¯ quarks contributions (in pb) to the LO s W +charm total cross-sections with NNLO PDFs for the leptonic decay channel W → eν are shown in Table II. The strange quark contributes the most to this W + charm production. With regard to the parametrization of the strange-quark content of the proton, CTEQ-TEA and MSTW2008 PDF groups make different assumptions in their global fits. In CT10NNLO

/negationslash and CT14NNLO, the strange is parameterized symmetrically s = ¯, and in MSTW2008, it is s parameterized asymmetrically, s -¯ = 0. s Hence corresponding theoretical predictions differ accordingly. For MSTW2008, the production of W + + c is slightly larger than the W + +¯, c as expected because of the s -¯ asymmetry. Because of the dominance of the s d quark over the ¯ -quark in the proton, the production of d W -+ c is larger than W + +¯. The numbers c in bracket correspond to p l T &gt; 35 GeV.

TABLE II: LO contributions of d , ¯ , d s and ¯ quarks to s σ pp ( → W + ) c × B W ( → lν ) within the kinematic region p jet T &gt; 25 GeV, | η jet | &lt; 2 5, . lepton pseudorapidity range | η l | &lt; 2 1 . and p l T &gt; 25(35) GeV.

| Subprocess         | CT10NNLO     | CT14NNLO     | MSTW2008NNLO   |
|--------------------|--------------|--------------|----------------|
| ¯+ s g → W + +¯ c  | 35.82(28.72) | 32.85(26.28) | 31.59(25.49)   |
| ¯ d + g → W + +¯ c | 2.33(1.89)   | 2.37(1.92)   | 2.43(1.96)     |
| s + g → W - + c    | 35.85(28.78) | 32.89(26.32) | 32.49(26.15)   |
| d + g → W - + c    | 4.50(3.73)   | 4.58(3.78)   | 4.66(3.86)     |

## B. Absolute and normalized differential cross-section

The absolute and normalized differential cross-sections are obtained by MadGraph using the same setup as the CMS collaboration at the LHC and a center-of-mass energy of 7TeV. In Fig.3 and Fig.4, we compare the absolute and normalized differential cross-sections in bins of lepton pseudo-rapidity with CMS measurements. The absolute and normalized differential cross-sections with PDF uncertainty at 68% C.L. are summarized in Table III and IV, and the CMS 7TeV measurement with statistical and systematic uncertainty is given in the last column. There is good agreement between theoretical predictions and measured distributions. The comparisons among predictions from various PDFs may lead to different conclusions. For instance, the predictions based on CT14NNLO and MSTW2008NNLO PDFs are smaller than the predictions based on CT10NNLO PDFs, and PDF uncertainties of CT14NNLO PDFs are much larger then the PDF uncertainties of CT10NNLO and MSTW2008NNLO.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00026v2-figure-bf0994c0-a6c2-4f05-a799-47794075643a.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00026v2-figure-880ec9ed-1643-4d4d-a0fb-ade309cd33e7.jpg)

l

FIG. 3: Comparison of theoretical predictions for differential cross-sections, dσ W ( + ) c /d η | | , as a function of absolute value of pseudorapidity of lepton from W-boson decay, with CMS measurements. Theoretical predictions at NLO are calculated using MadGraph with CT10NNLO, CT14NNLO and MSTW2008NNLO PDFs. The left figure shows the predictions for the lepton from W boson decay with p l T &gt; 25GeV, and the right with p l T &gt; 35GeV. Error bars on the theoretical predictions show 68% C.L..

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00026v2-figure-8f95c822-c244-4b7e-9c7f-e9aa7148be37.jpg)

l

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00026v2-figure-d73bc3a8-089c-4ca5-be6e-55fff7e0ba6d.jpg)

l

FIG. 4: Comparison of theoretical predictions for normalized differential cross-sections, dσ W ( + c /σd η ) | | , as a function of absolute value of pseudorapidity of lepton from W-boson decay, with CMS measurements. Theoretical predictions at NLO are calculated using MadGraph with CT10NNLO, CT14NNLO and MSTW2008NNLO PDFs. The left figure shows the predictions for the lepton from the W boson decay with p l T &gt; 25GeV, and the right with p l T &gt; 35GeV.

TABLE III: Theory predictions of differential cross-sections dσ pp ( → W + ) c × B W ( → lν ) /d η | l | with PDF uncertainty at 68 % C.L., last three column for CMS 7TeV measurement with statistical and systematic uncertainty.

| p l T > 25GeV    | p l T > 25GeV           | p l T > 25GeV            | p l T > 25GeV           | p l T > 25GeV            |
|------------------|-------------------------|--------------------------|-------------------------|--------------------------|
| | η l |          | CT10NNLO                | CT14NNLO                 | MSTW2008NNLO            | CMS measurement          |
| [0 , 0 . 35]     | 65 . 7 +6 . 1% - 5 . 7% | 60 . 5 +6 . 8% - 9 . 4%  | 59 . 2 +2 . 3% - 2 . 7% | 68 . 7 ± 3 . 9% ± 6 . 7% |
| [0 . 35 , 0 . 7] | 62 . 8 +6 . 2% - 5 . 6% | 58 . 1 +6 . 9% - 9 . 5%  | 57 . 0 +2 . 3% - 2 . 7% | 59 . 9 ± 4 . 2% ± 6 . 7% |
| [0 . 7 , 1 . 1]  | 56 . 9 +6 . 4% - 5 . 5% | 52 . 6 +7 . 0% - 9 . 8%  | 51 . 6 +2 . 2% - 2 . 6% | 56 . 7 ± 4 . 2% ± 6 . 7% |
| [1 . 1 , 1 . 6]  | 46 . 8 +7 . 0% - 5 . 7% | 43 . 8 +7 . 4% - 10 . 4% | 42 . 9 +2 . 2% - 2 . 7% | 44 . 8 ± 4 . 2% ± 7 . 1% |
| [1 . 6 , 2 . 1]  | 33 . 9 +8 . 5% - 6 . 2% | 32 . 0 +8 . 7% - 11 . 3% | 31 . 5 +2 . 4% - 2 . 9% | 35 . 1 ± 4 . 8% ± 6 . 8% |

p l T &gt; 35GeV

| | η l |          | CT10NNLO CT14NNLO       | CT10NNLO CT14NNLO        | MSTW2008NNLO            | CMS measurement          |
|------------------|-------------------------|--------------------------|-------------------------|--------------------------|
| [0 , 0 . 35]     | 53 . 4 +6 . 2% - 5 . 8% | 49 . 2 +6 . 9% - 9 . 5%  | 48 . 1 +2 . 3% - 2 . 7% | 52 . 3 ± 3 . 3% ± 6 . 1% |
| [0 . 35 , 0 . 7] | 50 . 9 +6 . 3% - 5 . 7% | 46 . 9 +6 . 9% - 9 . 7%  | 46 . 1 +2 . 3% - 2 . 7% | 49 . 2 ± 3 . 3% ± 6 . 1% |
| [0 . 7 , 1 . 1]  | 45 . 7 +6 . 6% - 5 . 6% | 42 . 2 +7 . 0% - 10 . 0% | 41 . 6 +2 . 2% - 2 . 6% | 45 . 5 ± 3 . 3% ± 5 . 9% |
| [1 . 1 , 1 . 6]  | 37 . 1 +7 . 2% - 5 . 7% | 34 . 5 +7 . 5% - 10 . 6% | 34 . 0 +2 . 2% - 2 . 7% | 34 . 2 ± 3 . 5% ± 6 . 1% |
| [1 . 6 , 2 . 1]  | 26 . 1 +8 . 9% - 6 . 3% | 24 . 6 +9 . 0% - 11 . 4% | 24 . 3 +2 . 4% - 3 . 0% | 26 . 6 ± 3 . 8% ± 6 . 4% |

TABLE IV: Theory predictions of normalized differential cross-sections (1 /σ W ( + )) c dσ W ( + ) c /d η | | with PDF uncertainty at 68 % C.L., three column depicts CMS 7TeV measurement with statistical and systematic uncertainty.

| p l T > 25GeV    | p l T > 25GeV            | p l T > 25GeV            | p l T > 25GeV            | p l T > 25GeV             |
|------------------|--------------------------|--------------------------|--------------------------|---------------------------|
| | η l |          | CT10NNLO                 | CT14NNLO                 | MSTW2008NNLO             | CMS measurement           |
| [0 , 0 . 35]     | 0 . 607 +1 . 2% - 2 . 2% | 0 . 602 +1 . 4% - 1 . 5% | 0 . 601 +0 . 7% - 0 . 6% | 0 . 638 ± 2 . 5% ± 1 . 9% |
| [0 . 35 , 0 . 7] | 0 . 581 +0 . 9% - 1 . 7% | 0 . 578 +1 . 1% - 1 . 2% | 0 . 578 +0 . 5% - 0 . 4% | 0 . 556 ± 2 . 9% ± 2 . 2% |
| [0 . 7 , 1 . 1]  | 0 . 527 +0 . 3% - 0 . 7% | 0 . 524 +0 . 4% - 0 . 5% | 0 . 524 +0 . 2% - 0 . 2% | 0 . 527 ± 2 . 8% ± 2 . 1% |
| [1 . 1 , 1 . 6]  | 0 . 433 +1 . 3% - 0 . 8% | 0 . 436 +0 . 8% - 0 . 9% | 0 . 435 +0 . 3% - 0 . 4% | 0 . 416 ± 2 . 9% ± 2 . 2% |
| [1 . 6 , 2 . 1]  | 0 . 314 +4 . 2% - 2 . 3% | 0 . 318 +3 . 0% - 2 . 5% | 0 . 320 +1 . 1% - 1 . 3% | 0 . 326 ± 3 . 7% ± 2 . 8% |

p l T &gt; 35GeV

| | η l |          | CT10NNLO CT14NNLO MSTW2008NNLO   | CT10NNLO CT14NNLO MSTW2008NNLO   | CT10NNLO CT14NNLO MSTW2008NNLO   | CMS measurement           |
|------------------|----------------------------------|----------------------------------|----------------------------------|---------------------------|
| [0 , 0 . 35]     | 0 . 618 +1 . 3% - 2 . 3%         | 0 . 615 +1 . 5% - 1 . 6%         | 0 . 611 +0 . 7% - 0 . 6%         | 0 . 622 ± 2 . 1% ± 1 . 6% |
| [0 . 35 , 0 . 7] | 0 . 589 +1 . 0% - 1 . 8%         | 0 . 586 +1 . 1% - 1 . 3%         | 0 . 585 +0 . 6% - 0 . 5%         | 0 . 585 ± 2 . 4% ± 1 . 7% |
| [0 . 7 , 1 . 1]  | 0 . 529 +0 . 4% - 0 . 7%         | 0 . 527 +0 . 4% - 0 . 5%         | 0 . 528 +0 . 3% - 0 . 2%         | 0 . 541 ± 2 . 2% ± 1 . 7% |
| [1 . 1 , 1 . 6]  | 0 . 429 +1 . 4% - 0 . 8%         | 0 . 431 +0 . 9% - 1 . 0%         | 0 . 432 +0 . 4% - 0 . 4%         | 0 . 407 ± 2 . 5% ± 2 . 0% |
| [1 . 6 , 2 . 1]  | 0 . 302 +4 . 7% - 2 . 5%         | 0 . 307 +3 . 4% - 2 . 7%         | 0 . 309 +1 . 3% - 1 . 5%         | 0 . 316 ± 3 . 2% ± 2 . 2% |

## C. Charged cross-section ratio

We calculated total ( σ W ( -+ ) c , σ ( W + + ¯)) and differential (absolute and normalized) c cross-sections independently under the same conditions in Subsections II A and II B. The CMS[11] collaboration introduced the charged cross-section ratio,

$$R ^ { c } = \frac { \sigma ( W ^ { + } + \bar { c } ) } { \sigma ( W ^ { - } + c ) }.$$

The advantage of using this ratio is that many of the theoretical and experimental uncertainties can cancel. The comparison of the total cross-section ratio and differential cross-section ratio with PDF uncertainty at 68% C.L. with CMS data are shown in Fig.5 and Fig.6, the left column corresponds to p l T &gt; 25 GeV, and right one is for p l T &gt; 35 GeV. The total cross-section ratio and differential cross-section ratio are also summarized in Table V and Table VI. From tables V and VI, we see that the total cross-section ratio, differential

cross-section ratio, and the associated PDF uncertainties are different for the CT10, CT14, and MSTW2008 PDF sets. These differences arise from the parametrization assumptions in each global analysis. For example, the CT10 and CT14 PDF sets assume s x, Q ( ) = ¯( s x, Q ), cross-section ratios almost exclusively are determined by the d -¯ asymmetry and with d a very small PDF uncertainty. In contrast, the MSTW08 PDF set assumes asymmetric strangeness s x, Q ( ) -¯( s x, Q ) = 0, that yields a larger PDF uncertainty in the prediction.

/negationslash

FIG. 5: Comparison of theoretical prediction of total cross-section ratio σ W ( + +¯) c /σ W ( -+ ) for c three different PDF sets with CMS measurements.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00026v2-figure-238606cb-51f0-4582-9f77-d6186c6fd748.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00026v2-figure-adc5dcc3-6b77-4d5c-96b2-c72cd1e80d17.jpg)

-

-

l

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00026v2-figure-c6f85605-eb2b-4a26-8ff9-da289f860cbc.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00026v2-figure-9a6d292b-66b5-481a-94d3-ac6c75b380f6.jpg)

l

FIG. 6: Comparison of theoretical prediction of differential cross-section ratio for three different PDF sets with CMS data.

TABLE V: Theoretical prediction of total cross-section ratio of σ W ( + +¯) c /σ W ( -+ ) c

|              | σ ( W + +¯) c /σ ( W - + c )    | σ ( W + +¯) c /σ ( W - + c )    |
|--------------|---------------------------------|---------------------------------|
| PDF sets     | p l T > 25GeV                   | p l T > 35GeV                   |
| CT10NNLO     | 0.944 +0 . 3% - 0 . 3%          | 0.942 +0 . 4% - 0 . 3%          |
| CT14NNLO     | 0.946 +0 . 4% - 0 . 7%          | 0.940 +0 . 4% - 0 . 7%          |
| MSTW2008NNLO | 0.920 +2 . 2% - 2 . 7%          | 0.916 +2 . 3% - 2 . 7%          |
| CMS          | 0.954 2.5% (stat.) 0.4% (syst.) | 0.938 2.0% (stat.) 0.6% (syst.) |

±

±

±

±

TABLE VI: Theoretical prediction of differential cross-section ratio with PDF uncertainty at 68 % C.L., last column depicts CMS 7TeV measurement with statistical and systematic uncertainty.

| p l T > 25GeV    | p l T > 25GeV            | p l T > 25GeV            | p l T > 25GeV            | p l T > 25GeV             |
|------------------|--------------------------|--------------------------|--------------------------|---------------------------|
| | η l |          | CT10NNLO                 | CT14NNLO                 | MSTW2008NNLO             | CMS measurement           |
| [0 , 0 . 35]     | 0 . 966 +0 . 2% - 0 . 2% | 0 . 977 +0 . 2% - 0 . 4% | 0 . 958 +1 . 2% - 2 . 8% | 1 . 013 ± 5 . 1% ± 0 . 5% |
| [0 . 35 , 0 . 7] | 0 . 968 +0 . 2% - 0 . 2% | 0 . 962 +0 . 2% - 0 . 4% | 0 . 953 +1 . 4% - 2 . 7% | 0 . 960 ± 5 . 5% ± 0 . 5% |
| [0 . 7 , 1 . 1]  | 0 . 959 +0 . 3% - 0 . 3% | 0 . 947 +0 . 3% - 0 . 6% | 0 . 928 +1 . 9% - 2 . 6% | 0 . 897 ± 5 . 7% ± 0 . 9% |
| [1 . 1 , 1 . 6]  | 0 . 935 +0 . 5% - 0 . 4% | 0 . 935 +0 . 5% - 0 . 9% | 0 . 898 +2 . 8% - 3 . 1% | 1 . 062 ± 5 . 7% ± 1 . 3% |
| [1 . 6 , 2 . 1]  | 0 . 881 +1 . 0% - 0 . 8% | 0 . 898 +1 . 0% - 1 . 5% | 0 . 852 +4 . 1% - 4 . 6% | 0 . 776 ± 7 . 5% ± 2 . 1% |

p l T &gt; 35GeV

| | η l |          | CT10NNLO CT14NNLO MSTW2008NNLO   | CT10NNLO CT14NNLO MSTW2008NNLO   | CT10NNLO CT14NNLO MSTW2008NNLO   | CMS measurement           |
|------------------|----------------------------------|----------------------------------|----------------------------------|---------------------------|
| [0 , 0 . 35]     | 0 . 966 +0 . 2% - 0 . 2%         | 0 . 970 +0 . 2% - 0 . 4%         | 0 . 955 +1 . 2% - 2 . 9%         | 0 . 993 ± 4 . 1% ± 0 . 7% |
| [0 . 35 , 0 . 7] | 0 . 968 +0 . 2% - 0 . 2%         | 0 . 957 +0 . 3% - 0 . 5%         | 0 . 949 +1 . 4% - 2 . 7%         | 0 . 977 ± 4 . 0% ± 0 . 7% |
| [0 . 7 , 1 . 1]  | 0 . 953 +0 . 3% - 0 . 3%         | 0 . 947 +0 . 3% - 0 . 6%         | 0 . 931 +2 . 0% - 2 . 6%         | 0 . 927 ± 4 . 3% ± 0 . 9% |
| [1 . 1 , 1 . 6]  | 0 . 931 +0 . 5% - 0 . 5%         | 0 . 930 +0 . 6% - 1 . 0%         | 0 . 892 +3 . 0% - 3 . 3%         | 0 . 948 ± 4 . 9% ± 1 . 1% |
| [1 . 6 , 2 . 1]  | 0 . 877 +1 . 1% - 0 . 8%         | 0 . 882 +1 . 2% - 1 . 7%         | 0 . 836 +4 . 5% - 5 . 0%         | 0 . 784 ± 6 . 4% ± 1 . 4% |

## D. Correlation between the CT14NNLO and predictions for W + charm data

One way to determine the sensitivity of a specific data point to some PDF f i ( x, Q ) at a given x and Q is to compute a correlation cosine between the theoretical prediction for this point and the PDFs of various flavors [16-18]. Therefore, we will study the correlations between CT14NNLO PDFs of various flavors at specific x and each data point of CMS 7TeV W + charm production with transverse momentum of the charged lepton from W boson decay at the p l T &gt; 35 GeV region. However first we briefly provide the definition of the correlation cosine. If there are two variables X /vector ( a j ) and Y ( /vector a j ) in the parameter space, where a j are the PDF parameters, then the correlation cosine can be expressed as:

$$\L o s \phi = \frac { \vec { \nabla } X \cdot \vec { \nabla } Y } { | \vec { \nabla } X | | \vec { \nabla } Y | },$$

where /vector ∇ X and /vector ∇ Y are gradient of the variables X a ( j ) and Y ( a j ). For X a ( i ), i th component of gradient vector is

$$( \vec { \nabla } X ) _ { j } = \frac { \partial X } { \partial a _ { j } } = \frac { 1 } { 2 } ( X ^ { + } ( a _ { j } ) - X ^ { - } ( a _ { j } ) ),$$

where X + ( a j ) and X -( a j ) are computed from the two sets of PDFs along the positive and negative direction of the i -th eigenvector. The quantity cosφ characterizes whether the variables X and Y are correlated ( cosφ ∼ 1), anti-correlated ( cosφ ∼ -1) or not correlated ( cosφ ∼ 0).

Fig.7 shows the correlation cosine between each data point and CT14NNLO PDF flavors at Q = 1 3 GeV and . Q = 100 GeV. CMS 7TeV W + charm data contains 17 data points, and thus there are 17 lines for certain flavors in Fig.7. Correlations of s x, Q ( ) PDF and g x, Q ( ) PDF with each data point are given at the first row in Fig.7. Correlations of u x, Q ( ) PDF and ¯( u x, Q ) PDF with each data point are given at the second row in Fig.7. Finally correlations of d x, Q ( ) PDF and ¯ ( d x, Q ) PDF with each data point are given at the third row in Fig.7. In each subfigure, the correlation between one of the PDF flavors with each data point is distinguished by different type of line. Solid, long-dashed-dotted, dotted, shortdashed, and short-dashed-dotted lines correspond to correlation of differential cross-section, differential cross-section ratio, normalized differential cross-section, total cross-section, and total cross-section ratio data respectively. As we discussed in Section II, differential crosssection, differential cross-section ratio and normalized differential cross-section data has included five data points that are measured by five rapidity bin ranges. The lines with darker color correspond to higher rapidity bin range.

FIG. 7: Correlation cos φ between CT14NNO PDFs at the specific x value on the horizontal axis and CT14NNLO predictions at NLO for CMS 7 TeV W + charm production at Q = 1 3 GeV(left . panel) and Q = 100 GeV(right panel).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00026v2-figure-3facdcae-9ec1-43f5-89b3-4634f340a363.jpg)

In the case of the total cross-section, differential cross-section and ratio, s x, Q ( ) PDF correlations are most significant (cos φ ∼ 1) at x from few times 10 -2 to few times 10 -1 , when Q = 1 3 GeV and . Q = 100 GeV respectively. However at other x range, the s x, Q ( ) PDFs correlation is not very strong. There are no clear relations between the rapidity bin range and correlation cosine, however it can be seen from two the subfigures at the first row

that each data point in various rapidity bins has a strong correlation with s x, Q ( ) PDF at the x -region mentioned above. There is other information for the normalized differential cross-section, which includes five data points that are partially correlated and partially anticorrelated, and represented with each flavor, as illustrated in each subfigure of Fig.7 with five dashed-lines (red or blue) that are inconsistent with other types of lines. Correlated data points prefer PDFs to become larger and anti-correlated data points prefer to PDFs to become smaller. Hence the impact of normalized cross-section data on PDFs is neutralized. Not considering the correlations with normalized cross-section data, the d (anti) quark and u (anti) PDFs are anti-correlated at x &lt; 0 4 region, . and small at all x -regions when both Q = 1 3 GeV and . Q = 100 GeV, which can be seen in subfigures in the second and third rows of Fig.7. Thus those data points without normalized cross-section data prefer d (anti) quark and u (anti) PDFs to become smaller in the x &lt; 0 4 region. . Gluon's correlation is likewise miniscule in the x -region which can be seen in the subfigures in the first row of Fig.7. One can conclude that CMS 7TeV W + charm data have larger impact on s x, Q ( ) PDF in CT14NNLO than other flavors.

## III. USING EPUMP TO STUDY THE IMPACT OF W + CHARM DATA ON CT14NNLO PDFS

The study cited in Ref.[19] presented a software package, ePump (error PDF updating method package), that can be used to update or optimize a set of PDFs, including the best-fit PDF set and error PDFs, and to update any other set of observables. Furthermore, Ref. [20] and Ref.[21] cite interesting further studies using ePump. In this section, we use ePump to analyze the impact of CMS 7TeV W + charm production measurements on the CT14NNLO PDFs. To update CT14NNLO PDFs, we use the CMS 7TeV total cross-section(one data point), differential cross-section(five data points), total cross-section ratio(one data point) and differential cross-section ratio(five data points), as well as combined data sets and their NLO QCD predictions from MadGraph as ePump inputs. CT14NNLO+sig, CT14NNLO+dsig, CT14NNLO+R, CT14NNLO+dR, and CT14NNLO+Wc in Figs.8-9 are the ePump-updated PDFs by total cross-section data, differential cross-section data, total cross-section ratio data, differential cross-section ratio data, and combined CMS 7TeV W + charm data. The weight factor for each data is three in our ePump studies. A weight larger

than one is equivalent to having more data points with the same experimental uncertainties or, alternatively, to reducing the experiment uncertainties by a factor of the square root of the weight. In the combined data, we excluded the normalized differential cross-section data to avoid double counting. After updating, the relative changes in CT14NNLO ensembles are best visualized by comparing their PDF error band and PDF ratio, in which ratio plot is obtained by dividing the error set and best fit of updated PDFs by the best fit of original CT14NNLO PDFs. In Figs.8-9, we show the impact of W + c data and combined data on CT14NNLO PDFs, namely, Comparisons of CT14NNLO PDF (light blue) and ePumpupdated PDFs CT14NNLO+sig (orange), CT14NNLO+dsig (magenta), CT14NNLO+R (purple), CT14NNLO+dR (green), and CT14NNLO+Wc(blue) at Q = 1 3 GeV (left col-. umn) and Q = 100 GeV (right column). Flavors g x, Q ( ), s x, Q ( ), u x, Q ( ), ¯( u x, Q ), d x, Q ( ), and ¯ ( d x, Q ) are shown. The PDF uncertainty bands are 90% C.L..

The change in g x, Q ( ) PDF mostly comes from the differential cross-section data. The changes in the best fit value and uncertainty of the ePump-updated g x, Q ( ) PDF are visualized in the first row of Fig.8, compared to CT14NNLO PDFs. The central value of the updated g x, Q ( ) PDF in the range 10 -1 &lt; x &lt; 0 5 remains almost unchanged, and it is . increased slightly in the range 10 -2 &lt; x &lt; 10 -1 , compared to CT14NNLO for Q = 1 3GeV, . whereas it is decreased by large factors for x &lt; 10 -3 and x &gt; 0 5. . For Q = 100GeV, the central value of the updated g x, Q ( ) PDF in the range x &lt; 0 5 remains almost unchanged, . and some variations of the best fit g x, Q ( ) in the region x &gt; 0 5 are observable. . The g x, Q ( ) PDF is not efficiently determined in very small and very large x -regions, however all remain within the error bands of the PDFs. The error band of the updated g x, Q ( ) PDF is slightly reduced in the range 10 -2 &lt; x &lt; 10 -1 for Q = 1 3GeV. In other regions, the updated . g x, Q ( ) PDFs uncertainty bands are comparable to that of CT14NNLO;

The s x, Q ( ) PDF is most sensitive to CMS 7TeV W + charm data. Total and differential cross-section data are responsible for most of the changes in the central values and uncertainties of the updated s x, Q ( ) PDF. The changes in central value and uncertainty of the ePump-updated s x, Q ( ) PDF is visualized in the second row of Fig.8, compared to CT14NNLO PDFs. After updating the PDF by combined data, the best fit s x, Q ( ) PDF is increased slightly in the region x &lt; 0 4 for . Q = 1 3GeV. There is a large change of the . s x, Q ( ) PDF in the region x &gt; 0 5 for both . Q = 1 3 GeV and . Q = 100 GeV, however well within the error bands of PDFs.

FIG. 8: Comparison of 90% C.L. g x, Q ( ) PDF (first row) and s x, Q ( ) PDF (second row) uncertainties from CT14NNLO and CT14NNLO+sig, CT14NNLO+dsig, CT14NNLO+R, CT14NNLO+dR, CT14NNLO+Wc. Shaded area represents error bands of CT14NNLO PDFs. The area between solid line represents error bands of ePump-updated PDFs, which are distinguished by different colors, and dotted line depicts the best fit PDFs. Both best fit PDFs and error bands are normalized to the corresponding central CT14NNLO PDFs.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00026v2-figure-06fa1345-8847-466c-88f3-8a915bf26cc8.jpg)

The LO contributions of s , ¯, s d , ¯ quarks to the d W + charm production cross-section are shown in Table II. d and ¯ quarks contributions are significantly small relative to d s and ¯ s quarks contributions, because dg → W -+ c and ¯ dg → W + +¯processes are suppressed by c the CKM matrix element. However we can not neglect d and ¯ quark contributions to the d W + charm cross-section, d quark contribution is about 11 % of the W -+ c productions cross-section and the ¯ quark contribution is about 6 % of the d W -+ c productions crosssection. Therefor, CMS 7TeV W + charm data can have an impact on both d x, Q ( ) and ¯ ( d x, Q ) PDFs. Fig. 9 shows the changes in central values and uncertainties of the updated d x, Q ( ) and ¯ ( d x, Q ) PDFs for both Q = 1 3 GeV and . Q = 100 GeV. Most of the changes in

d x, Q ( ) and ¯ ( d x, Q ) PDFs come from total and differential cross-section data of CMS 7TeV W + charm production. After updating the CT14NNLO PDFs by combined data, the error band of d x, Q ( ) PDF is slightly reduced in the region 10 -2 &lt; x &lt; 10 -1 for both Q = 1 3GeV . and Q = 100GeV. The best fit d x, Q ( ) PDF decreased a little bit in this region. For both Q = 1 3GeV and . Q = 100GeV, the central value and uncertainty of the ePump-updated ¯ ( d x, Q ) PDF are close to that of CT14NNLO in the region x &lt; 0 5; . The best fit ¯ ( d x, Q ) PDF is increased significantly in the region x &gt; 0 5, however well within the error bands of . PDFs.

FIG. 9: Comparison of 90% C.L. d x, Q ( ) PDF(first row) and ¯ ( d x, Q ) PDF(second row) uncertainties from CT14NNLO and CT14NNLO+sig, CT14NNLO+dsig, CT14NNLO+R, CT14NNLO+dR, CT14NNLO+Wc. Shaded area stands for error bands of CT14NNLO PDFs. The area between solid line stands for error bands of the epump updated PDFs, which are distinguished by different colors, and dotted line stands for best fit PDFs.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00026v2-figure-fa2c4a91-c7ad-4a5f-aa64-7988eaeaaf98.jpg)

At the LO, u quark does not contribute to the W + charm production cross-section, however it does so beyond the LO. In our ePump study, we employ the theoretical prediction

of the W +charm production cross-section at NLO QCD. Therefore, we compare the ePumpupdated PDFs via CMS 7TeV W + charm data and CT14NNLO PDFs to see the impact on u x, Q ( ) and ¯( u x, Q ) PDFs in CT14NNLO for both Q = 1 3 GeV and . Q = 100 GeV. We found that the central value and uncertainties of the ePump updated u x, Q ( ) and ¯( u x, Q ) are almost unchanged.

In Fig.10, we compared the s x, Q ( ) PDF from CT14NNLO, and ePump-updated s x, Q ( ) PDF from combined CMS 7TeV W + charm data with weights three amd ten. Fig. 10 shows that the s x, Q ( ) PDF error band greatly decrease at x &lt; 0 4 for . Q = 1 3 GeV and . Q = 100 GeV when the weight factor is increased from three to ten. The best fit s x, Q ( ) PDF increases in the region x &lt; 0 1. . At large x , the best fit s x, Q ( ) PDF decreases significantly for both Q = 1 3 GeV and . Q = 100 GeV, however remains well within the error bands of PDFs.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00026v2-figure-1045d3ed-e8dd-470c-9371-7b075f7fe8b1.jpg)

x

FIG. 10: Comparison of 90% C.L. s x, Q ( ) PDF uncertainties at Q = 1 3 GeV and . Q = 100 GeV from CT14NNLO and CT14NNLO+Wc3 and CT14NNLO+Wc10 corresponding to s x, Q ( ) PDF updated by combined CMS 7TeV W + charm data with weight factors three and ten.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00026v2-figure-f688577d-3cd0-441b-8598-f2b84899e190.jpg)

Using ePump [19], we updated CT14NNLO PDFs. One might want to know how the inclusion of the CMS 7TeV W + charm data in the global PDF fits would modify the prediction and uncertainties for any other set of observables including the original observables that were used for updating the CT14NNLO PDFs. ePump can also directly update predictions and uncertainties for any observables after including new data. In Fig.11, we compare the predictions from CT14NNLO PDFs and CT14NNLO+Wc PDFs, obtained by updating the CT14NNLO with CMS 7TeV W + charm data using ePump with the CMS 7TeV

W + charm data. Fig.11 shows that the uncertainties decreased after the update, and the predicted central value is also closer to the data.

FIG. 11: Comparison of data and theoretical predictions obtained from updated and original CT14NNLO PDFs, PDF uncertainty at 68% C.L.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00026v2-figure-1d01313f-7b88-4696-8ae9-70fb29e6cef8.jpg)

## IV. CONCLUSIONS

In this study we calculated the total and differential cross-sections and cross-section ratios using the MadGraph up to O α ( 2 s ) with a massive charm quark m c = 1 55 GeV for . three NNLO PDF sets: MSTW2008, CT10, and CT14NNLO. Subsequently, we compared the experimental measurements of W + charm production at √ s = 7 TeV at LHC. In our calculation, we use the same kinematic cuts as experimental measurements: p jet T &gt; 25 GeV, | η jet | &lt; 2 5, and . | η l | &lt; 2 1 GeV, and two different transverse momentum cuts . p l T &gt; 25GeV for the W → µν channel and p l T &gt; 35 GeV for the W → µν and W → eν channels. In our calculation, both the factorization and the renormalization scales are set to the value of the

W boson mass, and α s ( M Z ) is set to the central value provided by the respective PDF groups. Our results are summarized in Tables I - VI and in Figures 2 - 6, where the central value of the prediction and the PDF uncertainty are given. The theoretical predictions from various PDFs agree well with experimental measurements. However, there are some differences depending on the PDFs used in the calculations. For example, unlike the assumption in MSTW20018 NNLO PDFs, the CT10 and CT14 assume s = ¯ in the proton, yielding to a s total and differential cross-sections ratio dominated by the d -¯ asymmetry. The total and d differential cross-sections are larger for the W -+ c production than for W + + , because c the former process involves a d , whereas the latter involves ¯ (sea) antiquark. d Hence, both total and differential cross-section ratios are smaller than 1 0. .

Fig. 7 shows that the observable from the CMS 7TeV W + charm production has a strong correlation with the strange(anti) quark PDFs, therefore these measurements also provide a direct constraint on the strange(anti) quark content of the proton.

Furthermore, using the ePump updating method, and CMS 7TeV W + charm production data at lepton transverse momentum p l T &gt; 35 GeV, we find that these data sets mainly reduce the s x, Q ( ) PDF error band and increase magnitude of its best fit in the x &lt; 0 4 . region for both Q = 1 3 GeV and . Q = 100 GeV. In Fig.11, we also compare the predictions from CT14NNLO PDFs and CT14NNLO+Wc PDFs.

## Acknowledgments

We thank Tie-Jiun Hou many helpful discussions. This work is supported by the National Natural Science Foundation of China under the Grant No. 11965020.

[1] W. J. Stirling and E. Vryonidou, 'Charm production in association with an electroweak gauge boson at the LHC', Phys. Rev. Lett. 109 (2012) 082002. arXiv:1203.6781 .

[2] N. Cabibbo, 'Unitary Symmetry and Leptonic Decays,' Phys. Rev. Lett. 10 (1963) 531-533.

[3] U. Baur, F. Halzen, S. Keller, M.L. Mangano, and K. Riesselmann, Phys. Lett. B 318, 544 (1993), 544-548, arXiv:hep-ph/9308370 .

[4] S. Keller, W. T. Giele, and E. Laenen, Phys. Lett. B 372, 141 (1996).

- [5] H. L. Lai, P. M. Nadolsky, J. Pumplin, D. Stump, W. K. Tung, and C. P. Yuan, 'The Strange parton distribution of the nucleon: Global analysis and applications', JHEP 04 (2007) 089, arXiv:hep-ph/0702268 .
- [6] NuTeV Collaboration, M. Goncharov et al., 'Precise measurement of dimuon production cross-sections in muon neutrino Fe and muon anti-neutrino Fe deep inelastic scattering at the Tevatron', Phys. Rev. D64 (2001) 112006, arXiv:hep-ex/0102049 .
- [7] NuTeV Collaboration, D. Mason et al., Measurement of the Nucleon StrangeAntistrange 'Asymmetry at Next-to-Leading Order in QCD from NuTeV Dimuon Data', Phys. Rev. Lett. 99 (2007) 192001
- [8] T. Aaltonen et al. , CDF Collaboration, 'First measurement of the production of a W boson in association with a single charm quark in pp ¯ collisions at √ s = 1.96-TeV', Phys. Rev. Lett. 100 (2008) 091803, arXiv:0711.2901 .
- [9] V. M. Abazov et al. , D0 Collaboration, 'Measurement of the ratio of the pp ¯ → W c + -jet cross-section to the inclusive pp ¯ → W + jets cross-section', Phys. Lett. B666 (2008) 23-30, arXiv:0803.2259 .
- [10] ATLAS Collaboration, 'Measurement of the production of a W boson in association with a charm quark in pp collisions at 7TeV with the ATLAS detector', arXiv:1402.6263 .
- [11] S. Chatrchyan et al. , CMS Collaboration, 'Measurement of associated W + charm production in pp collisions at √ s = 7TeV', JHEP 02 (2014) 013, arXiv:1310.1138 .
- [12] J. Alwall, M. Herquet, F. Maltoni, O. Mattelaer, and T. Stelzer, 'MadGraph 5 : Going Beyond', JHEP 06 (2011) 128, arXiv:1106.0522 .
- [13] H.-L. Lai, M. Guzzi, J. Huston, Z. Li, P. M. Nadolsky, J. Pumplin, and C. P. Yuan, 'New parton distributions for collider physics', Phys. Rev. D82 (2010) 074024, arXiv:1007.2241 .
- [14] S. Dulat, T.-J. Hou, J. Gao, M. Guzzi, J. Huston, P. Nadolsky, J. Pumplin, C. Schmidt, D. Stump, and C. P. Yuan, 'New parton distribution functions from a global analysis of quantum chromodynamics', Phys. Rev. D93 (2016) 033006, arXiv:1506.07443 .
- [15] A. D. Martin, W. J. Stirling, R. S. Thorne, and G. Watt, 'Parton distributions for the LHC,' Eur. Phys. J. C63 (2009) 189-285, arXiv:0901.0002 .
- [16] J. Pumplin, D. Stump, R. Brock, D. Casey, J. Huston, J. Kalk, H. L. Lai and W. K. Tung, 'Uncertainties of predictions from parton distribution functions. 2. The Hessian method', Phys. Rev. D65 (2001) 014013, arXiv:hep-ph/0101032 .

- [17] P. M. Nadolsky and Z. Sullivan, 'PDF uncertainties in WH production at Tevatron', eConf C 010630 (2001) P510, arXiv:hep-ph/0110378 .
- [18] P. M. Nadolsky, H.-L. Lai, Q.-H. Cao, J. Huston, J. Pumplin, D. Stump, W.-K. Tung, and C. P. Yuan, 'Implications of CTEQ global analysis for collider observables,' Phys. Rev. D78 (2008) 013004, arXiv:0802.0007 .
- [19] C. Schmidt, J. Pumplin, C. P. Yuan, 'Updating and optimizing error parton distribution function sets in the Hessian approach', Phys. Rev. 98 , (2018) 094005, arXiv:1806.07950 .
- [20] C. Willis, R. Brock, D. Hayden, T. J. Hou, J. Isaacson, C. Schmidt and C. P. Yuan, 'New method for reducing parton distribution function uncertainties in the high-mass Drell-Yan spectrum', Phys. Rev. 99 , (2019) 054004, arXiv:1809.09481 .
- [21] T. J. Hou, Z. Yu, S. Dulat, C. Schmidt and C.-P. Yuan, 'Updating and Optimizing Error PDFs in the Hessian Approach. Part II', arXiv:1907.12177 .
Nijat Yalkun Sayipjamal Dulat
QCD analysis of CMS W + charm measurements at LHC with \sqrt s = 7TeV and implications for the strange PDF
We calculate cross-sections and cross-section ratios of a charm quark
production in association with a $W$ gauge boson at next-to-leading order QCD
using MadGraph and CT10NNLO, CT14NNLO, and MSTW2008NNLO PDFs. We compare the
results with measurements from the CMS detector at the LHC at a center-of-mass
energy of 7TeV. Moreover, we calculate absolute and normalized differential
cross-sections as well as differential cross-section ratios as a function of
the lepton pseudorapidity from the $W$ boson decay. The correlation between the
CT14NNLO PDFs and predictions for $W+$ charm data are studied as well.
Furthermore, by employing the error PDF updating method proposed by the
CTEQ-TEA group, we update CT14NNLO PDFs, and analyze the impact of CMS 7TeV
$W+$ charm production data to the original CT14NNLO PDFs. By comparison of the
$g(x,Q)$, $s(x,Q)$, $u(x,Q)$, $d(x,Q)$, $\bar u(x,Q)$, and $\bar d(x,Q)$ PDFs
at $Q=1.3$ GeV and $Q = 100$ GeV for the CT14NNLO and CT14NNLO+Wc, we see that
the error band of the $s(x,Q)$ PDF is reduced in the region $x<0.4$, and the
error band of $g(x,Q)$ PDF is also slightly reduced in the region $0.01 < x <
0.1$.