pith. sign in

arxiv: 2507.14637 · v2 · submitted 2025-07-19 · ✦ hep-ph · hep-ex

Determination of Fragmentation Functions from Charge Asymmetries in Hadron Production

Pith reviewed 2026-05-19 04:02 UTC · model grok-4.3

classification ✦ hep-ph hep-ex
keywords fragmentation functionsnon-singletcharge asymmetryQCDNNLOpionskaonsSIDIS
0
0 comments X p. Extension

The pith

Charge asymmetries in SIA and SIDIS isolate non-singlet fragmentation functions for pions and kaons at NNLO.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper proposes a method to extract non-singlet fragmentation functions of light charged hadrons directly from charge asymmetries observed in single-inclusive electron-positron annihilation and semi-inclusive deep-inelastic scattering. The extraction is performed at next-to-next-to-leading order in quantum chromodynamics with a full uncertainty analysis. The resulting functions display a scaling index near 0.7 at large momentum fractions and low scales, a strangeness suppression factor near 0.5, and universal behavior for light meson fragmentation. These functions supply benchmarks for non-perturbative QCD models and Monte Carlo generators while serving as input for future electron-ion collider experiments.

Core claim

We propose a novel method for extracting non-singlet (NS) fragmentation functions (FFs) of light charged hadrons from charge asymmetries measured in hadron fragmentation, using data from both single-inclusive electron-positron annihilation and semi-inclusive deep-inelastic scattering processes. We determine the NS FFs for pions and kaons at next-to-next-to-leading order in Quantum Chromodynamics, including a comprehensive uncertainty analysis. The extracted FFs reveal a scaling index of about 0.7 at large momentum fractions and low energy scales, a strangeness suppression factor of about 0.5, and universality in fragmentation of light mesons. Our findings provide a valuable benchmark for non

What carries the argument

Charge asymmetry measured simultaneously in SIA and SIDIS, which cancels singlet contributions to isolate the non-singlet fragmentation functions.

If this is right

  • The extracted non-singlet FFs provide a benchmark for testing non-perturbative QCD models and Monte Carlo event generators.
  • They serve as crucial input for future electron-ion collider analyses.
  • The scaling index of about 0.7 guides model building for fragmentation at large momentum fractions.
  • The strangeness suppression of about 0.5 and observed universality support simplified assumptions in light-meson hadronization studies.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If the isolation holds, certain global fits of parton distributions and fragmentation functions could be simplified by treating non-singlet channels separately.
  • Repeating the charge-asymmetry extraction with higher-precision data from future facilities would tighten the uncertainty bands on the scaling index.
  • The universality result suggests the same method could be applied to additional light mesons to check consistency of the suppression factor.

Load-bearing premise

Charge asymmetries in the measured SIA and SIDIS data directly isolate the non-singlet fragmentation functions with negligible contamination from singlet contributions, higher-twist effects, or process-specific corrections beyond the NNLO framework.

What would settle it

A high-precision measurement of charge asymmetries at large momentum fractions that deviates from the predictions of these extracted non-singlet FFs outside the reported uncertainty bands would indicate significant contamination or incorrect scaling.

Figures

Figures reproduced from arXiv: 2507.14637 by Bin Zhou, ChongYang Liu, Jun Gao.

Figure 1
Figure 1. Figure 1: FIG. 1: The global and individual [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: Comparison of SLD measurements of kaon charge [PITH_FULL_IMAGE:figures/full_fig_p004_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: Comparisons of experimental data and predictions from nominal fit and fit fixing [PITH_FULL_IMAGE:figures/full_fig_p008_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6: Same as Fig [PITH_FULL_IMAGE:figures/full_fig_p009_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7: Same as Fig [PITH_FULL_IMAGE:figures/full_fig_p010_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8: Same as Fig [PITH_FULL_IMAGE:figures/full_fig_p010_8.png] view at source ↗
read the original abstract

We propose a novel method for extracting non-singlet (NS) fragmentation functions (FFs) of light charged hadrons from charge asymmetries measured in hadron fragmentation, using data from both single-inclusive electron-positron annihilation and semi-inclusive deep-inelastic scattering processes. We determine the NS FFs for pions and kaons at next-to-next-to-leading order in Quantum Chromodynamics, including a comprehensive uncertainty analysis. The extracted FFs reveal a scaling index of about 0.7 at large momentum fractions and low energy scales, a strangeness suppression factor of about 0.5, and universality in fragmentation of light mesons. Our findings provide a valuable benchmark for testing non-perturbative QCD models and Monte Carlo event generators, and serve as crucial input for future electron-ion colliders.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript proposes extracting non-singlet fragmentation functions for pions and kaons from charge asymmetries in SIA and SIDIS data at NNLO in QCD. It performs a fit with comprehensive uncertainty analysis and reports a scaling index of approximately 0.7 at large z and low scales, a strangeness suppression factor of about 0.5, and evidence for universality in light-meson fragmentation. The results are positioned as benchmarks for non-perturbative QCD models and input for future electron-ion colliders.

Significance. If the central isolation of NS FFs holds, the work supplies useful NNLO constraints on fragmentation at low scales and could serve as a benchmark for Monte Carlo generators. The inclusion of NNLO coefficient functions and uncertainty quantification strengthens the potential utility for EIC phenomenology.

major comments (2)
  1. [§3 (charge-asymmetry definition and isolation)] The extraction rests on the claim that charge asymmetries isolate the non-singlet combination with negligible singlet mixing or higher-twist contamination. At the low scales and large z where the scaling index ~0.7 is reported, this isolation is least secure because SIDIS asymmetries can receive target-fragmentation or power-suppressed contributions that do not cancel symmetrically; NNLO evolution alone does not suppress them. A quantitative estimate or dedicated test of these residuals is required to support both the functional form and the universality conclusion.
  2. [§4.2 (fit results and parameter extraction)] The strangeness suppression factor of ~0.5 and the scaling index are obtained from the global fit to the chosen SIA+SIDIS data sets. The sensitivity of these parameters to data selection, kinematic cuts, and the precise form of the input parameterization at the initial scale must be shown explicitly; without this, the quoted values remain tied to the specific choices even though the two processes are formally independent.
minor comments (2)
  1. [Abstract and §4] The abstract states that a comprehensive uncertainty analysis is performed; the main text should explicitly list the sources (data, parameterization, higher-order truncation) and show how they are propagated into the final bands on the scaling index and suppression factor.
  2. [§5 (discussion)] A short comparison table placing the new NS FFs against existing global fits (e.g., DSS, NNPDF) at a common scale would help readers assess the numerical impact of the charge-asymmetry method.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the detailed and constructive report. We address the major comments below, providing clarifications and indicating where revisions will be made to strengthen the manuscript.

read point-by-point responses
  1. Referee: [§3 (charge-asymmetry definition and isolation)] The extraction rests on the claim that charge asymmetries isolate the non-singlet combination with negligible singlet mixing or higher-twist contamination. At the low scales and large z where the scaling index ~0.7 is reported, this isolation is least secure because SIDIS asymmetries can receive target-fragmentation or power-suppressed contributions that do not cancel symmetrically; NNLO evolution alone does not suppress them. A quantitative estimate or dedicated test of these residuals is required to support both the functional form and the universality conclusion.

    Authors: We acknowledge the importance of verifying the isolation of the non-singlet combination, particularly at low scales and large z. While the charge asymmetry is designed to suppress singlet contributions and symmetric higher-twist effects, we agree that a more explicit test is valuable. In the revised manuscript, we will add a dedicated subsection presenting a quantitative assessment by varying the kinematic cuts (e.g., minimum Q²) and examining the impact on the extracted scaling index and universality. This will help quantify any residual power-suppressed contributions. revision: yes

  2. Referee: [§4.2 (fit results and parameter extraction)] The strangeness suppression factor of ~0.5 and the scaling index are obtained from the global fit to the chosen SIA+SIDIS data sets. The sensitivity of these parameters to data selection, kinematic cuts, and the precise form of the input parameterization at the initial scale must be shown explicitly; without this, the quoted values remain tied to the specific choices even though the two processes are formally independent.

    Authors: We thank the referee for highlighting the need for explicit sensitivity studies. The manuscript includes an uncertainty analysis that accounts for variations in the parameterization, but we agree that additional explicit demonstrations would enhance clarity. In the revision, we will include supplementary material or figures showing the dependence of the strangeness suppression factor and scaling index on different data selections, kinematic cuts, and alternative input forms at the initial scale. revision: yes

Circularity Check

0 steps flagged

No significant circularity; extraction is data-driven and self-contained

full rationale

The paper proposes extracting non-singlet fragmentation functions directly from measured charge asymmetries in SIA and SIDIS data at NNLO, with results such as the scaling index ~0.7 and strangeness suppression ~0.5 obtained via fit to experimental datasets plus uncertainty analysis. No load-bearing step reduces by construction to the inputs: the isolation of NS combinations is an assumption justified by the processes' kinematics and coefficient functions rather than a self-definition or renamed fit; there are no self-citations invoked as uniqueness theorems, no ansatze smuggled via prior work, and no predictions presented as independent of the fitted data. The central claims remain independent of any circular reduction and rest on external experimental inputs.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The extraction rests on perturbative QCD validity at NNLO, the isolation of NS components via charge asymmetries, and the assumption of universality across processes. No new particles or forces are introduced.

free parameters (2)
  • scaling index = 0.7
    Fitted value of approximately 0.7 describing the large-z behavior at low scales
  • strangeness suppression factor = 0.5
    Fitted value of approximately 0.5 for kaon relative to pion fragmentation
axioms (2)
  • domain assumption Perturbative QCD calculations at NNLO accurately describe the hard scattering and evolution for the relevant processes
    Invoked for the determination of FFs from data
  • domain assumption Charge asymmetries in SIA and SIDIS are dominated by non-singlet fragmentation functions
    Core premise of the proposed extraction method

pith-pipeline@v0.9.0 · 5661 in / 1503 out tokens · 76123 ms · 2026-05-19T04:02:00.145770+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

75 extracted references · 75 canonical work pages · 32 internal anchors

  1. [1]

    S. M. Berman, J. D. Bjorken, and J. B. Kogut, Phys. Rev. D 4, 3388 (1971)

  2. [2]

    R. D. Field and R. P. Feynman, Nucl. Phys. B 136, 1 (1978)

  3. [3]

    R. P. Feynman, R. D. Field, and G. C. Fox, Phys. Rev. D 18, 3320 (1978)

  4. [4]

    J. C. Collins, D. E. Soper, and G. F. Sterman, Adv. Ser. Direct. High Energy Phys. 5, 1 (1989), hep-ph/0409313

  5. [5]

    Parton Fragmentation Functions

    A. Metz and A. Vossen, Prog. Part. Nucl. Phys. 91, 136 (2016), 1607.02521

  6. [6]

    Borsa, M

    I. Borsa, M. Stratmann, W. Vogelsang, D. de Florian, and R. Sassot, Phys. Rev. Lett. 133, 151901 (2024), 2407.11635

  7. [7]

    Cruz-Martinez, T

    J. Cruz-Martinez, T. Hasenack, F. Hekhorn, G. Magni, E. R. Nocera, T. R. Rabemananjara, J. Rojo, T. Sharma, and G. van Seeventer (2025), 2503.11814

  8. [8]

    Abdul Khalek, J

    R. Abdul Khalek, J. J. Ethier, J. Rojo, and G. van Weelden, JHEP 09, 183 (2020), 2006.14629

  9. [9]

    K. F. Muzakka et al., Phys. Rev. D 106, 074004 (2022), 2204.13157

  10. [10]

    K. J. Eskola, P. Paakkinen, H. Paukkunen, and C. A. Salgado, Eur. Phys. J. C 82, 413 (2022), 2112.12462

  11. [11]

    Helenius, M

    I. Helenius, M. Walt, and W. Vogelsang, Phys. Rev. D 105, 094031 (2022), 2112.11904

  12. [12]

    Electron Ion Collider: The Next QCD Frontier - Understanding the glue that binds us all

    A. Accardi et al., Eur. Phys. J. A 52, 268 (2016), 1212.1701

  13. [13]

    D. P. Anderle et al., Front. Phys. (Beijing) 16, 64701 (2021), 2102.09222

  14. [14]

    Global Analysis of Fragmentation Functions for Protons and Charged Hadrons

    D. de Florian, R. Sassot, and M. Stratmann, Phys. Rev. D 76, 074033 (2007), 0707.1506

  15. [15]

    Determination of fragmentation functions and their uncertainties

    M. Hirai, S. Kumano, T. H. Nagai, and K. Sudoh, Phys. Rev. D 75, 094009 (2007), hep-ph/0702250

  16. [16]

    AKK Update: Improvements from New Theoretical Input and Experimental Data

    S. Albino, B. A. Kniehl, and G. Kramer, Nucl. Phys. B 803, 42 (2008), 0803.2768

  17. [17]

    Charged hadron fragmentation functions from collider data

    V. Bertone, N. P. Hartland, E. R. Nocera, J. Rojo, and L. Rottoli (NNPDF), Eur. Phys. J. C 78, 651 (2018), [Erratum: Eur.Phys.J.C 84, 155 (2024)], 1807.03310

  18. [18]

    R. A. Khalek, V. Bertone, and E. R. Nocera (MAP (Multi-dimensional Analyses of Partonic distributions)), Phys. Rev. D 104, 034007 (2021), 2105.08725

  19. [19]

    Moffat, W

    E. Moffat, W. Melnitchouk, T. C. Rogers, and N. Sato (Jefferson Lab Angular Momentum (JAM)), Phys. Rev. D 104, 016015 (2021), 2101.04664

  20. [20]

    J. Gao, C. Liu, X. Shen, H. Xing, and Y. Zhao, Phys. Rev. Lett. 132, 261903 (2024), 2401.02781

  21. [21]

    J. Gao, C. Liu, X. Shen, H. Xing, and Y. Zhao, Phys. Rev. D 110, 114019 (2024), 2407.04422

  22. [22]

    J. Gao, C. Liu, M. Li, X. Shen, H. Xing, Y. Zhao, and Y. Zhou (2025), 2503.21311

  23. [23]

    A determination of the fragmentation functions of pions, kaons, and protons with faithful uncertainties

    V. Bertone, S. Carrazza, N. P. Hartland, E. R. Nocera, and J. Rojo (NNPDF), Eur. Phys. J. C 77, 516 (2017), 1706.07049

  24. [24]

    First global QCD analysis of charged hadron fragmentation functions and their uncertainties at next-to-next-to-leading order

    M. Soleymaninia, M. Goharipour, and H. Khanpour, Phys. Rev. D 98, 074002 (2018), 1805.04847

  25. [25]

    Borsa, R

    I. Borsa, R. Sassot, D. de Florian, M. Stratmann, and W. Vogelsang, Phys. Rev. Lett. 129, 012002 (2022), 2202.05060

  26. [26]

    Abdul Khalek, V

    R. Abdul Khalek, V. Bertone, A. Khoudli, and E. R. Nocera (MAP (Multi-dimensional Analyses of Partonic distributions)), Phys. Lett. B 834, 137456 (2022), 2204.10331

  27. [27]

    J. Gao, X. Shen, H. Xing, Y. Zhao, and B. Zhou (2025), 2502.17837

  28. [28]

    R. J. Holt and C. D. Roberts, Rev. Mod. Phys. 82, 2991 (2010), 1002.4666

  29. [29]

    T. Ito, W. Bentz, I. C. Cloet, A. W. Thomas, and K. Yazaki, Phys. Rev. D 80, 074008 (2009), 0906.5362

  30. [30]

    Heinzl, Lect

    T. Heinzl, Lect. Notes Phys. 572, 55 (2001), hep- th/0008096

  31. [31]

    Pion Structure Function in the Nambu and Jona-Lasinio model

    T. Shigetani, K. Suzuki, and H. Toki, Phys. Lett. B 308, 383 (1993), hep-ph/9402286

  32. [32]

    J. F. Gunion, 11th International Symposium on Multi- particle Dynamics (1980), SLAC-PUB-2607

  33. [33]

    S. J. Brodsky, M. Burkardt, and I. Schmidt, Nucl. Phys. B 441, 197 (1995), hep-ph/9401328

  34. [34]

    Yuan, Phys

    F. Yuan, Phys. Rev. D 69, 051501 (2004), hep- ph/0311288

  35. [35]

    Factorization of Large-x Quark Distributions in a Hadron

    X.-d. Ji, J.-P. Ma, and F. Yuan, Phys. Lett. B 610, 247 (2005), hep-ph/0411382

  36. [36]

    C. Shi, C. Mezrag, and H.-s. Zong, Phys. Rev. D 98, 054029 (2018), 1806.10232

  37. [37]

    Z. F. Cui, M. Ding, J. M. Morgado, K. Raya, D. Binosi, L. Chang, J. Papavassiliou, C. D. Roberts, J. Rodr´ ıguez- Quintero, and S. M. Schmidt, Eur. Phys. J. A 58, 10 (2022), 2112.09210

  38. [38]

    Goyal, S.-O

    S. Goyal, S.-O. Moch, V. Pathak, N. Rana, and V. Ravin- dran, Phys. Rev. Lett. 132, 251902 (2024), 2312.17711

  39. [39]

    Bonino, T

    L. Bonino, T. Gehrmann, and G. Stagnitto, Phys. Rev. Lett. 132, 251901 (2024), 2401.16281

  40. [40]

    Bonino, T

    L. Bonino, T. Gehrmann, M. L¨ ochner, K. Sch¨ onwald, and G. Stagnitto (2025), 2504.05376

  41. [41]

    Bonino, T

    L. Bonino, T. Gehrmann, M. L¨ ochner, K. Sch¨ onwald, and G. Stagnitto (2025), 2506.19926. 6

  42. [42]

    P. J. Rijken and W. L. van Neerven, Phys. Lett. B 386, 422 (1996), hep-ph/9604436

  43. [43]

    P. J. Rijken and W. L. van Neerven, Phys. Lett. B 392, 207 (1997), hep-ph/9609379

  44. [44]

    P. J. Rijken and W. L. van Neerven, Nucl. Phys. B 487, 233 (1997), hep-ph/9609377

  45. [45]

    QCD Corrections to Semi-Inclusive Hadron Production in Electron-Positron Annihilation at Two Loops

    A. Mitov and S.-O. Moch, Nucl. Phys. B 751, 18 (2006), hep-ph/0604160

  46. [46]

    G. Soar, S. Moch, J. A. M. Vermaseren, and A. Vogt, Nucl. Phys. B 832, 152 (2010), 0912.0369

  47. [47]

    A. A. Almasy, S. Moch, and A. Vogt, Nucl. Phys. B 854, 133 (2012), 1107.2263

  48. [48]

    Xu and H

    Z. Xu and H. X. Zhu (2024), 2411.11595

  49. [49]

    C.-Q. He, H. Xing, T.-Z. Yang, and H. X. Zhu (2025), 2503.20441

  50. [50]
  51. [51]
  52. [52]

    Multiplicities of charged kaons from deep-inelastic muon scattering off an isoscalar target

    C. Adolph et al. (COMPASS), Phys. Lett. B 767, 133 (2017), 1608.06760

  53. [53]

    G. D. Alexeev et al. (2024), 2410.12005

  54. [54]

    Allen et al

    P. Allen et al. (Aachen-Bonn-CERN-Munich-Oxford), Nucl. Phys. B 214, 369 (1983)

  55. [55]

    Production of pi+, pi-, K+, K-, p and p-bar in Light (uds), c and b Jets from Z0 Decays

    K. Abe et al. (SLD), Phys. Rev. D 69, 072003 (2004), hep-ex/0310017

  56. [56]

    Production of pi+, K+, K0, K*0, phi, p and Lambda-0 in Hadronic Z0 Decays

    K. Abe et al. (SLD), Phys. Rev. D 59, 052001 (1999), hep-ex/9805029

  57. [57]

    Zhou and J

    B. Zhou and J. Gao (2025), 2507.13185

  58. [58]

    Fully differential VBF Higgs production at NNLO

    M. Cacciari, F. A. Dreyer, A. Karlberg, G. P. Salam, and G. Zanderighi, Phys. Rev. Lett. 115, 082002 (2015), [Er- ratum: Phys.Rev.Lett. 120, 139901 (2018)], 1506.02660

  59. [59]

    Next-to-Next-to-Leading Order Evolution of Non-Singlet Fragmentation Functions

    A. Mitov, S. Moch, and A. Vogt, Phys. Lett. B 638, 61 (2006), hep-ph/0604053

  60. [60]

    On Third-Order Timelike Splitting Functions and Top-Mediated Higgs Decay into Hadrons

    S. Moch and A. Vogt, Phys. Lett. B 659, 290 (2008), 0709.3899

  61. [61]

    Chen, T.-Z

    H. Chen, T.-Z. Yang, H. X. Zhu, and Y. J. Zhu, Chin. Phys. C 45, 043101 (2021), 2006.10534

  62. [62]

    G. P. Salam and J. Rojo, Comput. Phys. Commun. 180, 120 (2009), 0804.3755

  63. [63]

    C. Liu, X. Shen, B. Zhou, and J. Gao, JHEP 09, 108 (2023), 2305.14620

  64. [64]

    Zhou and J

    B. Zhou and J. Gao, JHEP 02, 003 (2025), 2407.10059

  65. [65]
  66. [66]

    See Supplemental Material for detaild comparisions of the theory to data

  67. [67]

    Novikov et al., Phys

    I. Novikov et al., Phys. Rev. D 102, 014040 (2020), 2002.02902

  68. [68]

    L. Kotz, A. Courtoy, P. Nadolsky, F. Olness, and M. Ponce-Chavez, Phys. Rev. D 109, 074027 (2024), 2311.08447

  69. [69]

    Soft-Gluon Resummation and the Valence Parton Distribution Function of the Pion

    M. Aicher, A. Schafer, and W. Vogelsang, Phys. Rev. Lett. 105, 252003 (2010), 1009.2481

  70. [70]

    P. C. Barry, C.-R. Ji, N. Sato, and W. Melnitchouk (Jef- ferson Lab Angular Momentum (JAM)), Phys. Rev. Lett. 127, 232001 (2021), 2108.05822

  71. [71]

    Xing, W.-H

    H.-Y. Xing, W.-H. Bian, Z.-F. Cui, and C. D. Roberts (2025), 2504.08142

  72. [72]

    R. D. Ball et al. (NNPDF), Eur. Phys. J. C 82, 428 (2022), 2109.02653

  73. [73]

    Bailey, T

    S. Bailey, T. Cridge, L. A. Harland-Lang, A. D. Mar- tin, and R. S. Thorne, Eur. Phys. J. C 81, 341 (2021), 2012.04684

  74. [74]

    A comprehensive guide to the physics and usage of PYTHIA 8.3

    C. Bierlich et al., SciPost Phys. Codeb. 2022, 8 (2022), 2203.11601

  75. [75]

    Sjostrand, Comput

    T. Sjostrand, Comput. Phys. Commun. 82, 74 (1994). 7 Supplemental Material Comparison between Theory and Experimental Data This analysis incorporates charged hadron multiplicity data from several key experiments. The COMPASS collab- oration provides measurements on multiplicities of charged pions and kaons in kinematic bins of Bjorken-x from 0.004 to 0.4 ...