Charm production in charged current deep inelastic scattering at HERA

Charm production in charged current deep inelastic scattering has been measured for the first time in e±p collisions, using data collected with the ZEUS detector at HERA, corresponding to an integrated luminosity of 358 pb−1. Results are presented separately for e+p and e−p scattering at a centre-of-mass energy of s\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \sqrt{s} $$\end{document} = 318 GeV within a kinematic phase-space region of 200 GeV2 < Q2 < 60000 GeV2 and y < 0.9, where Q2 is the squared four-momentum transfer and y is the inelasticity. The measured cross sections of electroweak charm production are consistent with expectations from the Standard Model within the large statistical uncertainties.


Introduction
Measurements of heavy-flavour production serve as a good testing ground to investigate the predictive power of perturbative quantum chromodynamics (pQCD) as the large mass provides a natural hard scale. While charm production in neutral current deep inelastic scattering (NC DIS) and in photoproduction has been extensively studied at HERA, it has not been measured in charged current deep inelastic scattering (CC DIS) owing to its small cross section.
In CC DIS, single charm quarks in the final state already occur at the level of the Quark Parton Model (QPM) when either an incoming s or d quark is converted to a charm quark, or an incoming charm quark is converted to an s or d quark, as illustrated in figure 1 (i, ii). In the latter case, the single charm in the event arises from the associated charm quark in the proton remnant. In addition, single charm can arise from boson-gluon fusion (BGF) producing a cs (cd) quark pair. In this case, the incoming virtual W boson fuses with a gluon from the proton. The gluon splits into a ss (dd) or cc pair in the initial state, as shown in figure 1 (iii, iv). All these e + p processes lead to the same final state, e + p →ν e cs(d) X; this is also true for e − p, e − p → ν ec s(d) X. The characteristics of JHEP05(2019)201 (iii) (iv) Figure 1. Feynman diagrams of charm-production subprocesses in e + p collisions. The QPM process illustrated in (i) describes s(d) → c transitions. In the QPM process (ii)c →s(d), the charm in the final state arises from the associated charm quark in the proton remnant X. In the BGF processes, the incoming W boson couples to (iii) an ss(dd) or (iv) a cc pair from the gluon in the proton, producing a cs pair in the final state.
the events associated with these subprocesses and their association to particular kinematic configurations in the final state depend on the QCD scheme chosen, as detailed in the next section. The subprocess depicted in figure 1 (i) is directly sensitive to the strange-quark content of the proton and can be used to constrain it. However, the extraction of the relevant part of the cross section is model dependent.
In the SU(3) flavour model, a perfect symmetry is assumed between the three light flavours, which results in equal quark densities for the sea quark components in nucleons. This symmetry is broken if the strange-quark density is suppressed by the mass of the strange quark, as happens in the well established strange-quark suppression in fragmentation [1]. This symmetry breaking can also occur in the initial state, depending on x,

JHEP05(2019)201
the fraction of the proton momentum carried by the interacting parton. For larger values of x, some support for this has been found experimentally, such as in dimuon production in charged current by the CCFR [2] and NuTeV [3], as well as the NOMAD [4] and CHORUS [5] neutrino scattering experiments. However, the interpretation of these measurements depends on nuclear corrections and charm fragmentation and no consensus has emerged on the exact level of suppression as a function of x. Additionally, the recent highprecision measurements of inclusive W and Z production by the ATLAS collaboration [6] report an unsuppressed strange sea in the low-x regime. A similar result was obtained in a combined global QCD analysis of inclusive W and Z data from both the ATLAS and CMS experiments [7]. This observation was also supported by the analysis of the ATLAS W + c data [8]. However, the CMS W + c data [9,10] favour strangeness suppression also at low x. A re-evaluation of the LHC inclusive and W + c measurements and the neutrino scattering measurements by NOMAD [4] and CHORUS [5] has been performed [11,12], partly in an attempt to reconcile the factor-of-two discrepancy in the measured strange-quark densities. The resulting strange-quark parton distribution function (PDF) was reported to be inconsistent with the ATLAS fit [6].
This paper presents measurements of charm production in CC DIS in e ± p collisions using data from the HERA II data-taking period. The electroweak contribution to charmproduction cross sections is compared with several QCD schemes that are detailed in the following section.

Charm production in CC DIS at HERA
The kinematics of lepton-proton scattering can be described in terms of the Lorentzinvariant variables x Bj , y and Q 2 . The variable Q 2 is the negative squared four-momentum of the exchange boson −q 2 = −(k − k ) 2 , where k and k are the four-momenta of the incoming and outgoing lepton, respectively. The Bjorken-x scaling variable, x Bj , is defined as x Bj = Q 2 /(2p · q), where p is the four-momentum of the incoming proton. The variable y is the inelasticity defined as y = Q 2 /(sx Bj ), where s is the squared centre-of-mass energy of the collision.
The differential cross section of charm production in CC DIS at HERA, mediated by a W boson, can be expressed in terms of the proton structure functions F 2 , xF 3 and F L as follows [13] where G F is the Fermi coupling constant, M W is the mass of the W boson and Y ± = 1 ± (1 − y) 2 . The contribution from the longitudinal structure function, F L , vanishes except at values of y ≈ 1. The basic electroweak single-charm production mechanisms have been outlined in section 1. In the leading-order plus parton-shower Monte Carlo (MC) simulation, the core electroweak matrix elements are based on the QPM graphs in figure 1 (i, ii) and BGF-like configurations in figure 1 (iii, iv) through initial-state parton JHEP05(2019)201 showering. In addition, other tree-level higher-order processes are also added through leading-log (LL) parton showering. The electroweak matrix elements involving only light quarks are complemented by occasional final-state gluon splitting into cc pairs in the parton shower, as depicted in figure 2, with a cutoff mimicking charm-mass effects. At the singleevent level, if only one of the two charm quarks (or its resulting hadron) is detected and its charge is not measured (such as in the measurement technique used in this paper), then the contribution of this final-state QCD radiation is experimentally indistinguishable from electroweak production. The experimental measurement thus refers to a sum of all these processes, which make differing contributions to different regions of phase space, but cannot be disentangled with the presently available statistics.
In fixed-order QCD calculations, the final-state gluon-splitting contribution in figure 2 is formally of next-to-next-to-leading order (NNLO, O(α 2 s )) and thus not included in the next-to-leading-order (NLO, O(α s )) QCD predictions considered in this work, even though its contribution can be substantial. Contributions from QPM-like (figure 1 (i, ii)) and BGF-like (figure 1 (iii, iv)) processes are separated by the virtuality of the quark entering the electroweak process in relation to the chosen factorisation scale. The NLO corrections to figure 1 (i, ii) arise in the form of initial-or final-state gluon radiation, or a vertex correction.
In the zero-mass variable-flavour-number scheme (ZM-VFNS) [14,15], the charm part of the structure functions F c 2 and xF c 3 can be expressed in terms of different PDFs as follows 3)

JHEP05(2019)201
in e + p collisions, and (2.5) in e − p collisions. Here C i,j is the coefficient function for parton j in structure-function F i and d, s, c and g are respectively the down, strange, charm and gluon PDFs with the argument (x Bj , Q 2 ) omitted. The parameters |V ij | are the Cabbibo-Kobayashi-Maskawa matrix elements. Part of the effects beyond NLO are resummed at next-to-leading log in the zero-mass approximation in this scheme.
In the NLO fixed-flavour-number (FFN) scheme [16,17], charm-mass effects are treated explicitly up to O(α s ) in the matrix elements. In this scheme, there is no charm-quark content in the proton, thus the charm QPM graph in figure 1 (ii) and its associated higherorder corrections do not occur. This is compensated by a correspondingly larger gluon content in the proton, such that all initial-state charm contributions irrespective of scale are treated explicitly in the BGF matrix element (figure 1 (iv)). No resummation is performed.
In the FONLL-B scheme [18,19], a general-mass variable-flavour-number scheme, charm-mass effects are accounted for by interpolating between the ZM-VFNS and FFN predictions, such that all mass effects are correctly included up to O(α s ).
The xFitter framework [20] was used to interface the theoretical predictions. Predictions in the FFN scheme were obtained from OPENQCDRAD [21] using the ABMP 16.3 NLO PDF sets [22,23]. Predictions in the FONLL-B scheme were obtained from APFEL [24] with NNPDF3.1 [25]. The total uncertainties of the FFN and FONLL-B schemes were obtained by adding in quadrature the PDF, scale and charm-mass uncertainties.
In order to study the effects of strangeness suppression, the ZM-VFNS predictions were obtained from QCDNUM [26] with HERAPDF2.0 [27]. The strange-quark fraction, f s =s/(d+s), was chosen to vary in the range between a suppressed strange sea [28,29] and an unsuppressed strange sea [6,30]. In addition, two more variations of the assumptions about the strange sea were made. Instead of assuming that the strange contribution is a fixed fraction of the d-type sea, an x-dependent shape, xs = 0.5f s tanh(−20(x − 0.07)) xD, where xD = xd + xs, was used in which high-x strangeness is highly suppressed. This shape was suggested by HERMES measurements [31,32]. The value of f s was also varied between f s = 0.3 and f s = 0.5. The ZM-VFNS prediction was also evaluated with the ATLAS-epW Z16 PDF sets [6].

JHEP05(2019)201
A detailed description of the ZEUS detector can be found elsewhere [33]. A brief outline of the components that are most relevant for this analysis is given below.
Charged particles were tracked in the central tracking detector (CTD) [34][35][36], the microvertex detector (MVD) [37] and the straw-tube tracker (STT) [38]. The CTD and the MVD operated in a magnetic field of 1.43 T provided by a thin superconducting solenoid. The CTD drift chamber covered the polar-angle 1 region 15 • < θ < 164 • . The MVD silicon tracker consisted of a barrel (BMVD) and a forward (FMVD) section. The BMVD provided polar angle coverage for tracks with three measurements from 30 • to 150 • . The FMVD extended the polar-angle coverage in the forward region to 7 • . The STT covered the polar-angle region 5 • < θ < 25 • .
The high-resolution uranium-scintillator calorimeter (CAL) [39][40][41][42] consisted of three parts: the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part was subdivided transversely into towers and longitudinally into one electromagnetic section (EMC) and either one hadronic section in RCAL (RHAC) or two in BCAL and FCAL (B/FHAC1 and B/FHAC2). The smallest subdivision of the calorimeter was called a cell. The CAL energy resolutions, as measured under test-beam conditions, were σ(E)/E = 0.18/ √ E for electrons and σ(E)/E = 0.35/ √ E for hadrons, with E in GeV. The iron yoke surrounding the CAL was instrumented with proportional drift chambers to form the backing calorimeter (BAC) [43]. The BAC consisted of 5142 aluminium chambers inserted into the gaps between 7.3 cm thick iron plates (10, 9 and 7 layers in forecap, barrel and rearcap, respectively) serving as calorimeter absorber. The chambers were typically 5 m long and had a wire spacing of 1.5 cm. The anode wires were covered by 50 cm long cathode pads. The BAC was equipped with energy readout and position sensitive readout for muon tracking. The former was based on 1692 pad towers (50 × 50 cm 2 ), providing an energy resolution of ∼ 100%/ √ E, with E in GeV. The position information from the wires allowed the reconstruction of muon trajectories in two dimensions (XY in barrel and Y Z in endcaps) with spatial accuracy of a few mm.
The luminosity was measured using the Bethe-Heitler reaction ep → eγp by a luminosity detector which consisted of independent lead-scintillator calorimeter [44][45][46] and magnetic spectrometer [47] systems. The fractional systematic uncertainty on the measured luminosity was 2%.

Monte Carlo simulation
Inclusive CC DIS MC samples were generated to simulate the charm signal and lightflavour (LF) background. Neutral current DIS and photoproduction samples were used to simulate non-CC DIS backgrounds, which were found to be negligible after the CC selection defined below. The charged current events were generated with DJANGOH 1.6 [48], using the CTEQ5D PDF sets [49] including QED and QCD radiative effects at the par-1 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the nominal proton beam direction, referred to as the "forward direction", and the X axis pointing towards the centre of HERA. The coordinate origin is at the centre of the CTD. The pseudorapidity is defined as η = −ln tan θ 2 , where the polar angle, θ, is measured with respect to the Z axis.

JHEP05(2019)201
ton level. The ARIADNE 4.12 colour-dipole model [50] was used for parton showering. The Lund string model was used for hadronisation, as implemented in JETSET 7.4.1 [51]. The NC DIS events and photoproduction events were simulated by using DJANGOH and HERWIG 5.9 [52], respectively.

Reconstruction of kinematic variables
Charged current DIS at HERA produces a neutrino in the final state. The neutrino then escapes the ZEUS detector, resulting in a lack of information on the leptonic final state. Thus, the Lorentz-invariant kinematic variables must be defined with the hadronic final state. In the present analysis, this is done with the Jacquet-Blondel method, which assumes the four-momentum of the exchange-boson q to be equal not only to the difference in leptonic four-momentum k − k but also to that in hadronic four-momentum p − p . Then, the invariant variables described in section 2 can be reconstructed as is the hadronic E− P z variable with the sum extending over the energies, E i , and the longitudinal components of the momentum, p z,i of the reconstructed hadronic final-state particles, i. The quantity p T,h = | i p T,i | is the total transverse momentum of the hadronic final state with p T,i being the transverse-momentum vector of the particle i. The mean value of the difference between the true and reconstructed kinematic variables was found to be within ≈ 1% in the MC simulation study.

CC DIS selection
The ZEUS online three-level trigger system loosely selected CC DIS candidates based on calorimeter and tracking information [53,54]. The triggered events were then required to pass the following offline selection criteria to reject non-CC DIS events: • a kinematic selection cut was implemented at 200 GeV 2 < Q 2 JB < 60000 GeV 2 and y JB < 0.9 to confine the sample into a region with good resolution of the kinematic quantities and small background; • a characteristic of CC DIS events is the large missing transverse momentum, p T,miss , in the calorimeter due to the undetected final-state neutrino. Events were required to have p T,miss > 12 GeV and p T,miss > 10 GeV, where p T,miss is the missing transverse momentum, excluding measurements taken from the CAL cells adjacent to the forward beam hole;

JHEP05(2019)201
• further background rejection is discussed in detail in a dedicated study of CC DIS at ZEUS in the e + p scattering periods [55]. In addition, the remaining cosmic muons were removed by requiring the number of fired calorimeter cells N cell > 40 and comparing fractions of energy deposited in the EMC and HAC. Events with energy deposited in the RCAL, E RCAL > 2 GeV, were rejected if E RHAC /E RCAL > 0.5. Events with energy in the BCAL, A total of 4093 events in e + p data and 8895 events in e − p data passed these selection criteria. Comparisons of data and MC at the event-level selection stage are shown in figures 3 and 4 for e + p and e − p, respectively. The MC distribution is consistent with the data in both the e + p and e − p periods. From MC studies, the charm contribution to the CC events is expected to be about 25% in the e + p periods and 12% in the e − p periods and similar for both periods in terms of numbers of events.

Charm selection and signal extraction
Charm quarks in CC DIS events were tagged by using an inclusive lifetime method [56,57]. In CC DIS at HERA, LF production has the highest production rate and is the major source of background. The lifetime method uses the measurement of the decay length of the heavyflavour (HF) particle to discriminate between signal and background contributions. The underlying principle of this method [56] is that ground-state HF particles travel on average a measurable distance before they decay at a secondary vertex.
Jets were reconstructed from energy-flow objects [58,59], which combine the information from calorimetry and tracking, corrected for energy loss in the detector material. The k T clustering algorithm [60] was used with a radius parameter R = 1 in the longitudinally invariant mode [61,62]. The E-recombination scheme, which produces massive jets whose four-momenta are the sum of the four-momenta of the clustered objects, was used. Events were selected if they contained at least one jet with transverse energy, E jet T , greater than 5 GeV and within the jet pseudorapidity range −2.5 < η jet < 2.0 (1.5). 2 These selection criteria constrained the kinematic phase-space region of this analysis, along with the kinematic selection criteria at the event-level selection stage.
Tracks from the selected jets were required to have a transverse momentum, p trk T > 0.5 GeV, and the total number of hits in the MVD, N trk MVD ≥ 4 to reduce the effect of multiple scattering and ensure a good spatial resolution. If more than two such tracks were associated with the jet, a secondary-vertex candidate was fitted from the selected tracks using a deterministic annealing filter [63][64][65]. This fit provided the vertex position and its error matrix as well as the hadronic invariant mass, M secvtx , of the charged tracks associated with the reconstructed vertex. The charged-pion mass was assumed for all tracks when JHEP05(2019)201 calculating the vertex mass. The secondary-vertex candidates were required to satisfy the following criteria: • |z secvtx | < 30 cm, • M secvtx < 6 GeV, • ∆x 2 + ∆y 2 < 1 cm, where N trk secvtx is the number of tracks used to reconstruct the vertex, χ 2 /N dof is the goodness of the vertex fitting, z secvtx is the Z-coordinate of the secondary vertex and ∆x, ∆y are the X-and Y -displacement of the secondary vertex from the primary interaction vertex. These selection criteria ensure a good fit quality and high acceptance of the CTD and MVD for tracks used to reconstruct the vertices. The requirement on the track multiplicity was implemented in order to reduce the number of background vertices. Figures 5 and 6 show the distributions of the chosen jets and secondary-vertex candidates for the e + p and e − p periods, respectively. The transverse decay length of the selected secondary vertices was projected onto the jet axis. Due to the finite resolution of the MVD and the prompt production of LF particles, the distributions of the 2D decay length (L xy ) and the significance of the decay length (S = L xy /δ Lxy ) for LF jets were symmetric. In contrast, the distributions for HF jets, in this case containing charmed particles, were asymmetric, as illustrated in figures 7 and 8  (a, b). A very small contribution from beauty is also shown; this is treated as background. This enabled the LF background to be suppressed by subtracting the negative decay-length distribution from the positive decay-length distribution. The region around |L xy | = 0 or |S| = 0 is dominated by LF production, resulting in a large statistical uncertainty of the distribution due to subtraction of two large numbers. To optimise the precision of the extracted signal, vertex candidates were required to satisfy a significance threshold, |S| > 2. Figures 7 and 8 (c, d) illustrate the shape of the variable distributions after the background subtraction. The surviving events after the decay-length subtraction were used to extract charm cross sections in two bins of Q 2 .

Charm cross section
The lifetime method used in this analysis tags charm quarks regardless of their origin. Thus, the selected reactions include charm production from final-state gluon splitting, such as shown in figure 2, which is here denoted by QCD charm, in addition to the electroweak (EW) charm production discussed in section 2. In the present analysis, charm production was measured inclusively for 200 GeV 2 < Q 2 < 60000 GeV 2 and y < 0.9. Additionally, to reflect the detector acceptance, a visible phase-space region was defined as: 200 GeV 2 < Q 2 < 60000 GeV 2 , y < 0.9, E jet T > 5 GeV and −2. and absence of a charm-charge determination prevented an experimental separation of the different theoretical contributions. The visible charm-jet cross section, σ c,vis , was initially measured as follows: 1) where N data is the reconstructed number of charm-jet candidates in the data after the S + − S − subtraction, N MC bg is the background contribution and N MC c is the charm/anticharm contribution estimated from the MC. Here σ MC c,vis is the cross section of jets that σ MC c EW + σ(g → cc) 2 31 57 10 Table 1. MC contributions (%) of charm subprocesses to σ MC c,vis and σ MC c EW as predicted by ARI-ADNE. The first two columns (d → c and s → c for e + p collisions, for example) reflect the contributions from the QPM processes described in figure 1 (i) and a higher-order correction described in figure 1 (iii). The contribution of the final-state gluon splitting described in figure 2 enters the fourth column (g → cc).
are generated in the MC within the visible kinematic region and associated to a generated charm or anti-charm quark when ∆φ 2 + ∆η 2 < 1, where ∆φ and ∆η are, respectively, the azimuthal angle and pseudorapidity difference between the jet and the charm quark. Each charm quark was associated to the jet with the highest E jet T satisfying the above criteria and each such jet entered the visible cross section. The different processes contributing to σ MC c,vis as predicted by MC are given in table 1. The EW contribution in the charm-quark signal, σ EW c,vis , should be evaluated by subtracting the QCD contribution from gluon splitting (figure 2). However, the prediction from ARIADNE 4.12, like any prediction from gluon splitting in the massless mode with cutoff, cannot be considered to be reliable. Since the contribution predicted by ARIADNE (see table 1) is both small and imprecise, it was not subtracted but rather included in the systematic uncertainties. The visible jet cross section was extrapolated and converted to the total EW cross section via a factor C ext , calculated from the ratio of the number of charm events generated in the full kinematic range, N EW gen , to the number of charm jets of EW origin within the visible kinematic region, N EW vis : The resulting total EW charm cross section, σ c EW , is then given by This is predicted by the ARIADNE MC to be approximately 9 pb.

JHEP05(2019)201 7 Systematic uncertainties
Although the statistical power of the current data is limited, it is important for future studies to understand the limitations of the current method by careful evaluation of the systematic uncertainties. The sources of uncertainty and their estimated effects on the total EW charm cross sections provided in parentheses (δσ e + p , δσ e − p ) are: • Secondary vertex rescaling The MC samples used in this analysis produced a higher fraction of events with secondary vertices than the data. For the nominal result, N MC c and N MC bg in eq. (6.1) were reduced proportionally. For the systematic uncertainty, only N MC bg was rescaled (−1.2 pb, +0.9 pb).

• EW charm fraction
The MC predictions of the QCD contribution (figure 2) shown in table 1 of +6% for e + p collisions and +12% for e − p collisions were taken as systematic uncertainty (−0.6 pb, −1.1 pb).
• CC DIS selection The uncertainty due to the CC selection cuts was estimated by varying these cuts as in the previous ZEUS analysis [66] (±0.2 pb, ±0.1 pb).
• Jet energy scale The part of the transverse jet energy measured in the calorimeter in the MC was varied by its estimated uncertainty of ±3% (±0.0 pb, ±0.1 pb).
These uncertainties were added in quadrature. The uncertainty in the ZEUS luminosity measurement is ±2% and was not included in the results.
In addition, the effect of the significance cut, |S| > 2, was studied. Small changes in the value of the significance cut resulted in large changes of the extracted signal. This was found to be due to statistical fluctuations in the number of events in the region close to the |S| lower cut value. From a dedicated study, the effects on the cross sections were found to be as large as ±5 pb. As this result was still strongly affected by statistical fluctuations, which have been included in the quoted statistical uncertainty, it was not included in the systematic uncertainty.
Additionally, the uncertainty in the secondary-vertex selection method was estimated by reducing the requirement on the number of tracks, N trk secvtx , from three to two. The effects on the cross sections were found to be as large as +3 pb. This was again strongly affected by statistical fluctuations and not included in the systematic uncertainty.

Results
The charm-jet cross sections in CC DIS in e ± p collisions were measured in the visible kinematic phase space of 200 GeV 2 < Q 2 < 60000 GeV 2 , y < 0.9, E jet T > 5 GeV and −2.5 < η jet < 2.0 to be where the superscript ± denotes the charge of the incoming lepton. In addition, the cross sections were obtained for two separate Q 2 bins, 200 GeV 2 < Q 2 < 1500 GeV 2 and 1500 GeV 2 < Q 2 < 60000 GeV 2 , and are shown in figure 9. The total electroweak charm cross sections were found, following eq. (6.3), to be The QCD contribution to charm production was introduced as an additional systematic uncertainty. Theory predictions obtained at NLO QCD with the FFN and FONLL-B schemes are compared to the data in bins of Q 2 in figure 10.   strange-quark fraction are given in table 5. A further reduction of the theory uncertainty can be achieved in the future by including NNLO corrections [67].
The theory predictions in table 3 suggest that the most interesting subprocess, namely the QPM process depicted in figure 1 (i), contributes about 30 − 50% to the final EW cross section, depending on the kinematic range and QCD scheme used. In general, the data are well described by the theory predictions, however the large experimental uncertainties prevent a discrimination between the different models.   figure 1 (ii) enters in the third column (c →s(d)) with a higher-order correction from the BGF process in figure 1 (iv). For the FFN scheme, the process described in figure 1 (ii) does not participate. Thus the content of the third column is provided by the BGF process of figure 1 (iv) only.    Table 5. The NLO ZM-VFNS predictions with varied strange-quark fraction f s . Additionally, two x-dependent strange quark fractions were used as suggested by the HERMES collaboration. The ZM-VFNS predictions were also evaluated with the ATLAS-epW Z16 PDF set with an unsuppressed strange-quark content.

Summary and outlook
Measurements of charm production in charged current deep inelastic scattering in e ± p collisions have been performed based on HERA II data with an integrated luminosity of 358 pb −1 , which corresponds to e + p collisions with an integrated luminosity of 173 pb −1 and e − p collisions with an integrated luminosity of 185 pb −1 . Visible charm-jet cross sections for each lepton beam type were measured within a kinematic region 200 GeV 2 < Q 2 < 60000 GeV 2 , y < 0.9, E jet T > 5 GeV and −2.5 < η jet < 2.0. They were extrapolated to the EW cross sections given in the kinematic range 200 GeV 2 < Q 2 < 60000 GeV 2 and y < 0.9. Theoretical predictions with several assumptions about the strange-quark content of the proton and using different heavy-flavour schemes were found to be consistent with the data within the large experimental uncertainties. The analysis presented here shows the potential of DIS measurements to increase the knowledge about the strange-quark content of the proton. Future lepton-ion collider projects such as the electron-ion collider [68] or LHeC [69] will have much higher luminosity than HERA, accompanied by improved vertex detection capabilities. These projects should then be able to make an important contribution to the knowledge of the strange-quark content of the proton.

JHEP05(2019)201
Open Access. This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.