NLO + parton-shower generator for Wc production in the POWHEG BOX RES

We present the implementation of a next-to-leading-order plus parton-shower event generator for the hadronic production of a heavy charm quark accompanied by a leptonically-decaying $W$ boson in the POWHEG BOX RES framework. We consider both signatures, i.e. $pp\to c\, \ell^- \bar{\nu}_{\ell}$ and $pp\to \bar{c} \,\ell^+ {\nu}_{\ell}$, and we include exactly off-shell and spin-correlation effects, as well as off-diagonal Cabibbo-Kobayashi-Maskawa contributions. We present particle-level results, obtained interfacing our code with the Herwig7.2 and Pythia8.3 shower Monte Carlo event generators, including hadronization and underlying-event effects, and compare them against the data collected by the CMS Collaboration at $\sqrt{s}=13$ TeV.


Introduction
The production of a W boson in association with a charm quark (c), tagged by the full reconstruction of the charmed D hadron from the clear experimental signature of its decay products, is one of the main probe of the strange-quark content of the colliding particles at an hadron collider, such as the LHC.In fact, the leading contribution to W c production (W − c and W + c) comes from the scattering of a strange quark (s or s) and a gluon, from the two incoming hadrons, since the corresponding Cabibbo-Kobayashi-Maskawa (CKM) matrix element is the dominant one for this channel.The two leading-order (LO) Feynman diagrams for W c production are depicted in Fig. 1.
As first illustrated in Ref. [1], this process can be used to put constraints on the strangequark content of the colliding hadrons, i.e. protons at the Large Hadron Collider (LHC).Several measurements of W c production have already been performed at the LHC as it went through the increasing energy upgrade from 7 to 13 TeV, by the ATLAS [2,3], the CMS [4][5][6] and the LHCb Collaborations [7].
The next-to-leading order (NLO) QCD cross section for W c production has been known for a while, and studied both at the Tevatron [8] and at the LHC [9].More recently, the complete set of next-to-next-to-leading order (NNLO) QCD corrections to the dominant CKM-diagonal contribution have been computed in Ref. [10] and further extended with full CKM dependence, including the dominant NLO electro-weak corrections, in Ref. [11], in the massless-charm limit.
The first dedicated NLO QCD + parton shower (PS) generator for W c production was done in Ref. [12], within the PowHel event generator, based on the POWHEG method [13], with the W boson decaying leptonically, including charm-quark mass effects in the hardscattering matrix elements and a non-diagonal CKM matrix. 1  In this paper we describe the implementation of a NLO + parton-shower generator for the hadro-production of a massive charm quark accompanied by a leptonically-decaying W boson in the POWHEG BOX RES framework [17][18][19].We include exactly off-shell and spin-correlation effects, as well as off-diagonal CKM contributions.
The POWHEG BOX RES framework is able to deal with radiation off resonances.Although for the case at hand of QCD corrections to a leptonically decaying W boson, the machinery of radiation off resonances is not necessary, since no QCD radiation can be emitted from the leptons, we have implemented W c production in this framework for several reasons: 1. Better handling of final-state radiation from heavy quarks, as radiation collinear to heavy partons is dealt with an appropriate importance sampling. 2  2. The future inclusion of NLO electro-weak corrections is facilitated, since the process is already in the right framework to deal with photon radiation, not only from quarks, but also from the charged lepton in the W decay, where the virtuality of the resonance has to be preserved during the generation of the hardest radiation performed by the POWHEG BOX RES.
3. All the processes for which the NNLO+PS accuracy has been reached [21][22][23][24] using the MiNNLO PS formalism [25,26] are implemented in the POWHEG BOX RES framework.It would then be easier to implement also the NNLO QCD corrections to W c production in this framework too.
The paper is organized as follows.In Sect. 2 we give some technical details about the process we are studying at LO and NLO, we describe the change of scheme of our calculation and the default choice of the renormalization and factorization scales that we have used.In addition we give the numerical value of the input parameters we have used in our simulation, and we present the validation of the fixed-order NLO results.In Sect. 3 we give details of the matching of the NLO amplitudes with the POWHEG method, as implemented in the POWHEG BOX RES framework.We discuss the effects of the extra degrees of freedom offered by this matching and their numerical impact.In Sect. 4 we present a few interesting kinematic distributions computed with events obtained by completing the first radiation emission performed by the POWHEG BOX RES code, with three showering models, 1 Predictions for W c hadro-production, at NLO QCD accuracy, matched to PS, according to the MC@NLO matching framework [14], can also be obtained using the MadGraph5 aMC@NLO framework [15].In addition, one can also use the authomatic interface between the POWHEG BOX and MadGraph5 aMC@NLO [16] to build the code, according to the POWHEG matching framework. 2The first treatment of radiation from heavy parton in the POWHEG BOX appeared in Ref. [20].
-2 -implemented in Herwig7.2 and Pythia8.3, and we discuss the differences.We also present results with full hadronization in place and in the presence of underlying events, and we compare them with the experimental data collected by the CMS Collaboration at 13 TeV.
The POWHEG BOX RES framework, together with the Wc generator, can be downloaded at http://powhegbox.mib.infn.it.

Contributing processes and technical details
In this section we discuss the implementation of the processes pp → − ν c and pp → + ν c at NLO QCD accuracy.We consider the c quark to be massive, with all the other lighter quarks and anti-quarks treated as massless.At LO, the partonic subprocesses that contribute to W c production are where the channels involving a down (anti-)quark are Cabibbo-suppressed with respect to those containing the strange (anti-)quark.We are neglecting contributions arising from a bottom quark in the initial state: these contributions are suppressed both by the bottom parton distribution function (PDF), and also by the smallness of the CKM matrix element |V cb |.
The tree-level Feynman diagrams contributing to the process g s(d) → − ν c are represented in Fig. 1: this process can proceed either via an s-channel exchange of a s(d) quark (left panel), or via a t-channel exchange of a c quark (right panel).The mass of the charm quark ensures that the latter contribution is finite also when its transverse momentum is vanishing.
-3 -In the following we discuss only W − production, since a similar discussion is also valid for W + production, after charge conjugation of the involved subprocesses.Together with the QCD virtual corrections to the subprocesses in Eq. (2.1), not depicted here, also the real corrections contribute at the next-to-leading order.Examples of Feynman diagrams contributing to the real corrections for the process pp → − ν c j are shown in Fig. 2. We can separate the real-correction contributions into two classes: 1.The corrections of the type pp → − ν c j, with j = g, u(ū), d( d), s(s), that are singular in the limit of the unresolved jet 2) where q = u, d, s.Examples of these contributions are depicted in Fig. 2, diagrams (a), (b) and (c).We will refer to these terms with R in the following.
2. The corrections of the type q q → − ν c q , where q = u, d, s and q = s, d, which do not have any singularities associated with the extra light-parton emission, and are called regular, in the POWHEG jargon and will be indicated with R reg .Again we do not consider contributions with a final-state b quark.An example of this type of contributions is depicted in Fig. 2 (d).
Notice that the square of the diagrams contributing to the regular processes (Fig. 2 (d)), as well as similar diagrams with two gluons in the initial state (that interfere with the one in (c)), suffer from another potential source of divergence, due to the presence of the charm anti-quark in the s channel.In fact, when the invariant mass of the system comprising the s( d) quark and the two leptons is equal to the charm mass m c , these diagrams are divergent, due to the fact that the c-quark propagator goes on shell.In order to avoid this singularity, we can impose a technical cut on the invariant mass of the off-shell W boson in the theoretical simulation.In order to make a closer contact to experimentally accessible quantities, a more realistic cut is in general imposed on the transverse mass of the W system, m T W , defined as where ∆ϕ is the azimuthal separation between the transverse momenta of the charged lepton and the missing energy of the neutrino.We neglect contributions with two charmed quark in the final state, such as ūd → − ν cc, coming from the splitting g → cc, since these contributions are subtracted from experimental data in studies aimed at extracting information on the strange-quark parton distribution function, as illustrated in Sect.4.2.
The Born and virtual matrix elements, as well as the Born colour-and spin-correlated amplitudes, necessary to the POWHEG BOX RES to automatically implement the POWHEG formalism, have been computed analytically. 3The one-loop amplitudes have been validated against the Gosam code [27,28].The matrix elements for the real contributions were instead generated using Madgraph4 [29] and its interface with the POWHEG BOX [30].

The decoupling and MS schemes
When performing a fixed-order calculation with massive quarks, one can define al least two consistent renormalization schemes that describe the same physics: the usual MS scheme, where all flavours are treated on equal footing, and a mixed scheme [31], called decoupling scheme [32], in which the n lf light flavours are subtracted in the MS scheme, while the heavy-flavour loop is subtracted at zero momentum.In this scheme, the heavy flavour decouples at low energies.This is the scheme we have used in the renormalization of our analytic calculation, with n lf = 3.
If the decoupling scheme is chosen, then the strong coupling constant α S runs with three light flavours, and the parton distribution functions should not include the charm quark in the evolution.
To make contact with other results expressed in terms of the MS strong coupling constant, running with five light flavours, and with PDFs with five flavours, we prefer to change our renormalization scheme and to switch to the MS one.The procedure for such a switch is well known, and was discussed in Ref. [33].This procedure, applied to the b quark, is also illustrated in Appendix A of Ref. [34].
For the case under study, since the LO process contains only one power of α S , if we want to express our calculation in terms of a coupling constant defined in the MS scheme with five active flavours, we need to add the following contribution to the cross section where B is the Born cross section, m b the bottom mass and µ R is the renormalization scale.
A corresponding modification has to be applied if we want to employ PDFs with five active flavours.The gluon PDF, which appears in the LO cross section, induces a correction where µ F is the factorization scale.The quark PDF receives only corrections that starts at order α 2 S , and can thus be neglected.By adding Eqs.(2.6) and (2.7), we get the conversion factor for the two schemes that turns out to be independent from the bottom and charm masses, and different from zero only if the renormalization and factorization scales are different.

Renormalization and factorization scale settings
In the POWHEG BOX there is the option to use a different renormalization and factorization scale when computing the real contribution (and/or the subtraction terms, i.e. their soft and collinear limits) or the Born and the virtual ones.The only requirement is that the scale used in the evaluation of the real contributions (and possibly of the subtraction terms) must tend to the scale used in the Born amplitude, when the emitted parton gets unresolved.
The default central scale used in our simulation in the evaluation of the Born, virtual, real and subtraction terms is evaluated with the kinematics of the underlying-Born configuration.The regular real term R reg has no underlying Born configuration, and the scale H T is computed using the kinematics of the real emission contribution, adding also the transverse momentum of the massless emitted parton.
Other options contemplate the use of the scale (2.11) also in the real term and/or in the subtraction terms, 4 or a fixed scale, 5 such as µ = m W .The renormalization and factorization scales, µ R and µ F , are defined multiplying µ by the factors ξ R and ξ F respectively.The uncertainty due to missing higher order corrections is estimated considering the 7-point scale variations (1,2), (1, 0.5), (2, 1), (0.5, 1), (2, 2), (0.5, 0.5). (2.12) We employ the NNPDF3.1 nlo pdfas [35,36] PDF set, which is also used to determine the running of the strong coupling. 6

Numerical inputs
We have simulated W c production in pp collisions at √ s = 13 TeV.The input parameters we have used are The electromagnetic coupling is obtained in the G µ scheme and the values of relevant CKM matrix elements are In our calculation we have neglected the mass of the leptons.Only at the stage of event generation a leptonic mass is assigned to the charged lepton once the full kinematics of the event is generated. 7The assignment of the mass is done through a reshuffling procedure that preserves the total energy and momentum, and the virtuality of the resonances.The values of the lepton masses we have used are given by m e = 0.511 × 10 −3 GeV, m µ = 0.1057 GeV, m τ = 1.777GeV. (2.16) Furthermore, following the discussion that precedes Eq. (2.4), we always apply a technical loose cut, at the generation stage, on the W -boson virtuality m ν > 2 GeV .
(2.17  In this section we compare our fixed-order results with those obtained with the MCFM code [37,38].In order to perform the comparison we employ a fixed renormalization and factorization scale set to As we are setting µ F = µ R , the term ∆V defined in Eq. (2.8) vanishes and there is no correction factor.We define the fiducial cross section imposing the following cuts on the transverse momentum and rapidity of the charged lepton, as well on the transverse momentum of the charm quark These cuts corresponds to the ones employed by the CMS Collaboration to define the inclusive cross section in the analysis described in Ref. [6] at √ s = 13 TeV.
As discussed in Sect.2, around Eqs. (2.4) and (2.5), in order to avoid the charmpropagator singularity in diagrams like the one in Fig. 2 (d), we also implement a cut on the transverse mass of the W system m T W > 20 GeV. (2.20) In with MCFM.In all cases a sub-permille agreement is found.Good agreement is found also at the differential level, as it can be seen in Fig. 3, where we illustrated the transverse mass of the reconstructed W boson (left panel) and the transverse momentum of the charmed quark (right panel).

Details of the POWHEG matching
In this section we briefly recall the expression for the differential cross section generated with the POWHEG method [13] for the process under consideration, focusing for simplicity only on the pp → ν c signature, in order to discuss further sources of theoretical uncertainty intrinsic to the POWHEG method.
The NLO cross section reads where B, V , R and R reg are the Born, virtual, real and regular real matrix elements, multiplied by the appropriate PDF and flux factors.As stated in Sect.2, we indicate with R reg the terms arising from the flavour configurations of Eq. (2.3), which do not contain any QCD singularity, while the processes in Eq. (2.2) are designated with R. The Born and virtual cross sections are evaluated using a three-body phase space Φ 3 , while, in the real contributions, an additional light parton is present.The matrix elements depend on the µ R scale, as the LO starts with one power of the strong coupling constant α S , while the µ F dependence arises from the incoming-parton PDFs.To build the POWHEG cross section, the real term R is split into a singular R s and a finite contribution R f , and the B term is defined as where Φ 1 is the radiation phase space, defined in terms of three radiation variables.The probability of not having any resolved emissions up to a scale k T is given by where the factorization and the renormalization scales are evaluated at the scale k T , which corresponds to the transverse momentum of the light parton, if it is collinear to an initialstate parton, while, for gluon emissions from the final-state charm quark, it is given by where p g and p c are the gluon and charm four-momenta and E g and E c their energies in the partonic center-of-mass frame.
The POWHEG cross section finally reads where k min T is an infrared cutoff.There are several degrees of freedom in the implementation of Eq. (3.5).One of them is in the choice of the renormalization and factorization scales µ R and µ F in the R f and R reg terms, as well as in the term R s appearing in B in Eq. (3.2), as discussed in Sect.2.2.Another one is in the choice of how to separate R into a singular part R s and a finite one R f .Without loss of generality, one can separate the two contributions introducing what is called a "damping function", f (k T ), that satisfies underlying-Born kinematics f 2 (h = 0.1) R 0 real kinematics 1 Table 2: Setups for the implementation of the POWHEG cross section given in Eq. (3.5).
The second column denotes the kinematic configuration used to evaluate the renormalization and factorization scales in the computation of the real and subtraction terms, while the last column shows which profile function is used to evaluate Eq. (3.7).The functional forms of f 1 and f 2 are given in Eqs.(3.8) and (3.9) respectively.
For the process we are studying, we have considered two different forms for f where µ is defined in Eq. (2.9), and 0 ≤ h ≤ 1/2 is a free parameter.The value of x is set to 0 if k T < m c , so that no finite contribution is present with k T below the charm mass. 8he functional form f 2 in Eq. (3.9), with h = 0.3, is the default damping function in the Herwig7.2event generator [39], while f 1 in Eq. (3.8) corresponds to the default option of the POWHEG BOX RES.In Fig. 4 we provide a graphical representation of f 1 , f 2 (h = 0.1) and f 2 (h = 0). 9n the following we investigate the dependence of several inclusive (with respect to the QCD radiation) observables, such as the leptonic observables and the rapidity and the transverse momentum of the jet containing the charm quark.
We reconstruct jets using the Fastjet [40] implementation of the anti-k T algorithm [41] with R = 0.5.To this aim, we have considered the setups listed in Table 2 to compute the POWHEG cross section given in Eq. (3.5).The central value of the renormalization and factorization scales is set to H T /2 (see Eq. (2.9)).
In the following figures we will compare some kinematic distributions computed with six different setups: two fixed-order NLO results, one computed using the underlying-Born kinematics (NLO B ) for the renormalization and factorization scale in the real contributions, and the other one using the full real-event kinematics, indicated with NLO R .The other four distributions are computed analysing the POWHEG-generated events at the "Les Houches" (LHE) level, i.e. events with the first hardest radiation generated according to Eq. (3.5), adopting the setup scales and damping functions of Table 2.For inclusive observables, the LHE results should agree with the fixed-order NLO ones up to NNLO terms.2. The grey bands are obtained performing the 7-point scale variations in Eq. (2.12) of the NLO R results.The statistical errors from the integration procedure are also shown.
In Fig. 5 we plot the transverse mass of the W boson (left panel) and the rapidity of the charmed jet (right panel).The NLO B results (blue line) are roughly 3% higher than the NLO R ones (black line), although always contained in the grey scale-uncertainty band of NLO R , computed performing the 7-point scale variations in Eq. (2.12).Furthermore, all the distributions at the LHE level agree with the corresponding fixed order calculation, within the scale-variation band.
The transverse momentum of the W boson, p ⊥W , and of the charmed jet, p ⊥jc , are illustrated on the left and right panel of Fig. 6, respectively.Comparing the two NLO predictions computed with different scales, one can see that the use of the underlying-Born scales favours a harder spectrum, although the result always lies within the scale-variation band of the NLO R one, except for the very last bin.This is due to the fact that, for large c-jet transverse momentum, it is quite likely to produce a c → cg splitting that enhances the value of H T .As such, the underlying-Born kinematics yields a lower scale and hence a -12 - larger value for α S , that enhances the cross section.The same behaviour can be seen also at the LHE level.In the experimental analyses, such as the one in Ref. [6], an important role is played by the measurement of the transverse momentum and rapidity of the charged lepton.In Fig. 7 we then plot these quantities, after applying the cuts of Eq. (2.19), that are the same as those used by the CMS Collaboration, as well as a cut on m T W , defined in Eq. (2.20).As far as the fixed-order NLO results is concerned, the NLO B curves lie a few percent above the NLO R one, both for the transverse momentum and the rapidity.This is again due to the higher value of α S when the renormalization scale is evaluated with the underlying-Born kinematics.Differences up to 10% can also be seen among the LHE curves, but, in general, all these curves lie within the grey uncertainty band of the NLO R result, so that they are consistent among each other.
We thus conclude that no significant differences are found among the several damping factors, particularly when realistic analysis cuts are applied, and that the main theoretical uncertainties arise from scale variations. 10

NLO + parton-shower results
In this section, we present full results after the completion of the POWHEG shower performed by general-purpose Monte Carlo (MC) programs, such as Herwig7.2 [42,43] and Pythia8.3 [44,45], including hadronization and underlying-event effects.
We have showered the Les Houches events that we have produced using the POWHEG BOX RES code, employing the Herwig7.2angular-ordered shower (that we label as "QTilde" shower), and the default Pythia8.3 shower with fully local recoil [46] (that we denote as "Dipole" shower), as well as the Pythia implementation of the Vincia shower [47,48].For all the different types of shower we use the default setting of their parameters, and we only change the charm mass to agree with the value we have used to generate the W c sample (see Eq.(2.13)).In Sect.4.1 we provide additional details of the showers we have employed, while in Sect.4.2 we investigate the impact of several ingredients provided by the general-purpose MC generators on the simulation of W c events.Finally, in Sect.

NLO + PS results with hadronization and underlying-event effects
In this section we discuss the impact of the parton shower, of the hadronization and of the underlying events on the results generated using our simulation in the POWHEG BOX RES.
The results presented in this section do not include the QED shower.We focus on the µ + ν µ c signature at √ s = 13 TeV.We use the values of the parameters described in Sect.2.3, and the B 2 setting of Table 2, i.e. the central values for the renormalization and factorization scales are computed using the underlying-Born kinematics, according to Eq. (2.10), and as damping function we use the f 2 profile function of Eq. (3.9) with h = 0.1.We begin by analyzing the effect of the parton shower on the POWHEG BOX distributions (PWG).If, at the end of the shower, an event contains more than one muon (or anti-muon) passing the charged leptonic cuts, only the hardest one is considered to reconstruct the W kinematics.
Since the aim of our analysis is to enhance the sensitivity of the measured cross section to the strange quark PDF, we proceed as follows: if more than one c quark or charmed hadron is present, they are all considered, with the exception of hadrons containing a cc pair, which are counted as not charmed.In addition, when in an event the charmed-quark and the muon charges have the same sign (SS), the weight of the event is subtracted from the differential cross sections, while, when they have opposite sign (OS), the weight is added.This procedure effectively removes the background due to events where the finalstate charm comes from gluon splitting g → cc, and not from an initial-state strange quark emitting a W boson. Charmed quarks originated from the gluon splitting have the same probability to be SS or OS, so their contribution cancels, when the above procedure is applied.
Jets are reconstructed using the anti-k T algorithm with R = 0.5.The charm content of a jet is defined as N c − N c, being N c (N c) the number of c(c) quarks clustered in the jet, and the event weight that we use to fill a jet-related distribution is further weighted by the factor ±(N c − N c), where we use the + sign when the muon has negative charge, and the − sign otherwise.
In the following, we consider both inclusive results (with the technical cut on the W boson virtuality of Eq. (2.17) always in place), as well as more exclusive ones, where we apply the following acceptance cuts used by the CMS Collaboration in their analyses.At the inclusive level, the parton shower acts non-trivially on the transverse-momentum distribution of the charm quark and of the charmed jet, as portrayed in Fig. 8.In particular, one can notice that the hard tail is depleted, and the region around 5 GeV is instead increased.The effect is milder for the c-jet, as radiation inside the cone with jet radius R = 0.5 is still collected in the jet, but more evident for the hadron: in fact, in the high transverse-momentum region, the ratio of the LHE results with the showered ones reaches values around 50% for the quark/hadron, and 10% for the charmed jet.
In Fig. 11 we then compare the full simulations obtained with Herwig7.2 and Pythia8.3 against the PWG standalone results.With the exception of the very low transversemomentum region, Herwig7.2 and Pythia8.3 distributions are in very good agreement.In particular we notice that now the charmed-hadron distributions, obtained with all the three showers, are essentially indistinguishable for p ⊥c > 10 GeV, while in Fig. 8 the Herwig7.2result was quite different from the Pythia8.3 ones.The improved agreement at the particle level is not unexpected, as the results without the inclusion of the hadronization are very sensitive to the shower cutoff, which is tuned alongside the hadronization parameters.In Fig. 12 we plot the pseudo-rapidity11 (left panel) and the transverse momentum (right panel) of the muon, after applying the cuts in Eq. (4.1).Leptonic distributions are indirectly affected by the shower, that induces a modification of the momenta of the charmed quarks/hadrons.In particular, the softening of the charm quarks induced by the hadronization reduces the cross section.This reduction is clearly visible on the left panel of Fig. 12, with an almost uniform reduction of the cross section in the whole pseudo-rapidity range, of roughly 15%, with respect to the LHE result.In addition, all the three shower results are in very good agreement with each others.As far as the transverse momentum of the muon is concerned, shown on the right panel, we observe a significant shape variation as the bulk of the distribution, peaked around p ⊥µ ≈ m W /2, gets depleted.Differences among the three showers amount to few percents, and are much smaller than the scale-variation band.
been tuned in Ref. [51] using LEP data at the Z mass.Due to the fact that B mesons primarily decay into D ones, the charm hadronization is further affected by its interplay with the bottom hadronization.Thus, a dedicated tune will certainly allow for a better comparison with the data.

Conclusions
In this article we have presented the implementation of a NLO + parton-shower generator for massive-charm production in association with a leptonically-decaying W boson in hadron collisions in the POWHEG BOX RES framework.This process is particularly crucial for the determination of the strange-quark content of the proton.
All spin-correlation effects and off-shell contributions in the W decay are properly included, and also non-diagonal elements in the Cabibbo-Kobayashi-Maskawa matrix.
The leading-order and the one-loop matrix elements have been computed analytically, allowing for a very fast and numerically-stable evaluation of the amplitudes.
Using the flexibility of the POWHEG method on the possibility to separate the realradiation contribution into a finite and a divergent part, that is then exponentiated in the Sudakov form factor, we have implemented different ways to separate the real term and assessed the sensitivity of several kinematic distributions to this extra degree of freedom.
The main source of theoretical uncertainty remains the dependence on the factorization and renormalization scales, of the order of 10% for most kinematic distributions.
Events generated with our code have been interfaced with the Herwig7.2QTilde shower, as well with the Pythia8.3 default shower and the Vincia one.We have investigated the effect of the parton shower, hadronization and of the underlying-event activity on several kinematic distributions.In particular hadronization corrections turned out to have the highest impact on the differential cross sections.
Particle-level results obtained with the three shower models we have used agree very well with each other, once hadronization corrections and the underlying-event simulation are also included.
Finally we compared our results with the CMS data for the W D * signature, finding good agreement, within the scale-variation bands.

Figure 1 :
Figure 1: Feynman diagrams contributing to the partonic process g s(d) → − ν c at LO, in the s channel (a), on the left, and t channel (b), on the right.The heavy charm quark is represented as a double line.

Figure 2 :
Figure 2: Examples of Feynman diagrams contributing to the real corrections for the process pp → − ν c j, where j = g, u, d, s, ū, d, s.The last Feynman diagram (d) is an example of "regular" contribution, i.e. no singularities are associated with this process.

Figure 3 :
Figure 3: Differential distribution for the transverse mass of the reconstructed W − boson for the processes pp → e − νe c (left) and the transverse momentum of the heavy anti-charm quark for the processes pp → e + ν e c (right) at NLO accuracy.The input parameters are discussed in Sect.2.3.The ratio between predictions obtained with MCFM (blue) and the POWHEG BOX RES (red) is also shown.The statistical error bars from the Monte Carlo integrator are also shown.

Figure 4 :
Figure 4: Comparison of three damping functions used to separate the real contribution into a singular and a finite part, as detailed in Sect.3.

Figure 5 :
Figure 5: The transverse mass distribution of the reconstructed W boson (left) and the rapidity of the charmed jet (right) at NLO (NLO B and NLO R ) and at the Les Houches event level, for the choices of the scales and damping functions listed in Table2.The grey bands are obtained performing the 7-point scale variations in Eq. (2.12) of the NLO R results.The statistical errors from the integration procedure are also shown.

Figure 6 :
Figure 6: Same as Fig. 5 but for the transverse momentum of the reconstructed W boson (left) and of the charmed jet (right).

Figure 7 :
Figure 7: Same as Fig. 5 but for the transverse momentum (left) and rapidity (right) of the charged lepton, with the cuts in Eq. (2.19) and (2.20) in place.

4 . 1 3 4. 1 . 1
4.3 we compare our results against the experimental data at √ s = 13 GeV, taken by the CMS Collaboration.The POWHEG BOX RES matching with Herwig7.2 and Pythia8.QTilde shower The relevant parameters to shower the POWHEG-generated events with the QTilde shower are taken from the LHE-POWHEG.ininput card distributed in the public release of the Herwig7.2code.In particular the options set /Herwig/Shower/ShowerHandler:MaxPtIsMuF Yes set /Herwig/Shower/ShowerHandler:RestrictPhasespace Yes instruct the shower to veto all the emissions with transverse momentum larger than the scalup variable, which is the hardness of the POWHEG emission.We set the charm mass to 1.5 GeV set /Herwig/Particle/c:NominalMass 1.5*GeV set /Herwig/Particle/cbar:NominalMass 1.5*GeV

Figure 12 :
Figure 12: Same as Fig. 11 for the muon pseudo-rapidity (left) and transverse momentum (right), and using the cuts in Eq. (4.1).

Table 1 :
19ble 1 we show the integrated cross section at LO and NLO obtained with our code and σ [fb]e − νe c @ LO e + ν e c @ LO e − νe c @ NLO e + ν e c @ NLO Fiducial cross section obtained with MCFM and our generator (PWG RES) for the processes pp → e − νe c and e + ν e c at LO (left columns) and NLO (right columns) at √ s = 13 TeV within the fiducial cuts of Eqs.(2.19) and (2.20), with µ F = µ R = m W .The input parameters are discussed in Sect.2.3.In the last row we show the ratio between the MCFM and the POWHEG BOX RES predictions.