\PHnumber2017–xxx \PHdateAugust 27, 2019

\CollaborationThe COMPASS Collaboration \ShortAuthorThe COMPASS Collaboration

A semi-inclusive measurement of charged hadron multiplicities in deep inelastic muon scattering off an isoscalar target was performed using data collected by the COMPASS Collaboration at CERN. The following kinematic domain is covered by the data: photon virtuality (GeV/), invariant mass of the hadronic system GeV/, Bjorken scaling variable in the range , fraction of the virtual photon energy carried by the hadron in the range , square of the hadron transverse momentum with respect to the virtual photon direction in the range 0.02 (GeV/ (GeV/). The multiplicities are presented as a function of in three-dimensional bins of , , and compared to previous semi-inclusive measurements. We explore the small- region, i.e. (GeV/), where hadron transverse momenta are expected to arise from non-perturbative effects, and also the domain of larger , where contributions from higher-order perturbative QCD are expected to dominate. The multiplicities are fitted using a single-exponential function at small to study the dependence of the average transverse momentum on , and . The power-law behaviour of the multiplicities at large is investigated using various functional forms. The fits describe the data reasonably well over the full measured range.

Keywords: TMDs, Deep inelastic scattering, hadron multiplicities, unpolarised SIDIS

\Submitted(to be submitted to Phys. Rev. D)

## The COMPASS Collaboration

M. Aghasyan\Irefntriest_i, M.G. Alexeev\Irefnturin_u, G.D. Alexeev\Irefndubna, A. Amoroso\Irefnnturin_uturin_i, V. Andrieux\Irefnnillinoissaclay, N.V. Anfimov\Irefndubna, V. Anosov\Irefndubna, A. Antoshkin\Irefndubna, K. Augsten\Irefnndubnapraguectu, W. Augustyniak\Irefnwarsaw, A. Austregesilo\Irefnmunichtu, C.D.R. Azevedo\Irefnaveiro, B. Badełek\Irefnwarsawu, F. Balestra\Irefnnturin_uturin_i, M. Ball\Irefnbonniskp, J. Barth\Irefnbonnpi, R. Beck\Irefnbonniskp, Y. Bedfer\Irefnsaclay, J. Bernhard\Irefnnmainzcern, K. Bicker\Irefnnmunichtucern, E. R. Bielert\Irefncern, R. Birsa\Irefntriest_i, M. Bodlak\Irefnpraguecu, P. Bordalo\Irefnlisbon\Arefa, F. Bradamante\Irefnntriest_utriest_i, A. Bressan\Irefnntriest_utriest_i, M. Büchele\Irefnfreiburg, V.E. Burtsev\Irefntomsk, L. Capozza\Irefnsaclay, W.-C. Chang\Irefntaipei, C. Chatterjee\Irefncalcutta, M. Chiosso\Irefnnturin_uturin_i, I. Choi\Irefnillinois, A.G. Chumakov\Irefntomsk, S.-U. Chung\Irefnmunichtu\Arefb, A. Cicuttin\Irefntriest_i\Arefictp, M.L. Crespo\Irefntriest_i\Arefictp, Q. Curiel\Irefnsaclay, S. Dalla Torre\Irefntriest_i, S.S. Dasgupta\Irefncalcutta, S. Dasgupta\Irefnntriest_utriest_i, O.Yu. Denisov\Irefnturin_i, L. Dhara\Irefncalcutta, S.V. Donskov\Irefnprotvino, N. Doshita\Irefnyamagata, Ch. Dreisbach\Irefnmunichtu, W. Dünnweber\Arefsr, R.R. Dusaev\Irefntomsk, M. Dziewiecki\Irefnwarsawtu, A. Efremov\Irefndubna\Arefo, P.D. Eversheim\Irefnbonniskp, M. Faessler\Arefsr, A. Ferrero\Irefnsaclay, M. Finger\Irefnpraguecu, M. Finger jr.\Irefnpraguecu, H. Fischer\Irefnfreiburg, C. Franco\Irefnlisbon, N. du Fresne von Hohenesche\Irefnnmainzcern, J.M. Friedrich\Irefnmunichtu, V. Frolov\Irefnndubnacern, E. Fuchey\Irefnsaclay\Arefp2i, F. Gautheron\Irefnbochum, O.P. Gavrichtchouk\Irefndubna, S. Gerassimov\Irefnnmoscowlpimunichtu, J. Giarra\Irefnmainz, F. Giordano\Irefnillinois, I. Gnesi\Irefnnturin_uturin_i, M. Gorzellik\Irefnfreiburg\Arefc, A. Grasso\Irefnnturin_uturin_i, A. Gridin\Irefndubna, M. Grosse Perdekamp\Irefnillinois, B. Grube\Irefnmunichtu, T. Grussenmeyer\Irefnfreiburg, A. Guskov\Irefndubna, D. Hahne\Irefnbonnpi, G. Hamar\Irefntriest_i, D. von Harrach\Irefnmainz, F.H. Heinsius\Irefnfreiburg, R. Heitz\Irefnillinois, F. Herrmann\Irefnfreiburg, N. Horikawa\Irefnnagoya\Arefd, N. d’Hose\Irefnsaclay, C.-Y. Hsieh\Irefntaipei\Arefx, S. Huber\Irefnmunichtu, S. Ishimoto\Irefnyamagata\Arefe, A. Ivanov\Irefnnturin_uturin_i, T. Iwata\Irefnyamagata, V. Jary\Irefnpraguectu, R. Joosten\Irefnbonniskp, P. Jörg\Irefnfreiburg, E. Kabuß\Irefnmainz, A. Kerbizi\Irefnntriest_utriest_i, B. Ketzer\Irefnbonniskp,G.V. Khaustov\Irefnprotvino, Yu.A. Khokhlov\Irefnprotvino\Arefg, Yu. Kisselev\Irefndubna, F. Klein\Irefnbonnpi, J.H. Koivuniemi\Irefnnbochumillinois, V.N. Kolosov\Irefnprotvino, K. Kondo\Irefnyamagata, K. Königsmann\Irefnfreiburg, I. Konorov\Irefnnmoscowlpimunichtu, V.F. Konstantinov\Irefnprotvino, A.M. Kotzinian\Irefnturin_i\Arefyerevan, O.M. Kouznetsov\Irefndubna, Z. Kral\Irefnpraguectu, M. Krämer\Irefnmunichtu, P. Kremser\Irefnfreiburg, F. Krinner\Irefnmunichtu, Z.V. Kroumchtein\Irefndubna\Deceased, Y. Kulinich\Irefnillinois, F. Kunne\Irefnsaclay, K. Kurek\Irefnwarsaw, R.P. Kurjata\Irefnwarsawtu, I.I. Kuznetsov\Irefntomsk, A. Kveton\Irefnpraguectu, A.A. Lednev\Irefnprotvino\Deceased, E.A. Levchenko\Irefntomsk, M. Levillain\Irefnsaclay, S. Levorato\Irefntriest_i, Y.-S. Lian\Irefntaipei\Arefy, J. Lichtenstadt\Irefntelaviv, R. Longo\Irefnnturin_uturin_i, V.E. Lyubovitskij\Irefntomsk, A. Maggiora\Irefnturin_i, A. Magnon\Irefnillinois, N. Makins\Irefnillinois, N. Makke\Irefntriest_i\Arefictp\CorAuth, G.K. Mallot\Irefncern, S.A. Mamon\Irefntomsk, B. Marianski\Irefnwarsaw, A. Martin\Irefnntriest_utriest_i, J. Marzec\Irefnwarsawtu, J. Matoušek\Irefnnntriest_utriest_ipraguecu, H. Matsuda\Irefnyamagata, T. Matsuda\Irefnmiyazaki, G.V. Meshcheryakov\Irefndubna, M. Meyer\Irefnnillinoissaclay, W. Meyer\Irefnbochum, Yu.V. Mikhailov\Irefnprotvino, M. Mikhasenko\Irefnbonniskp, E. Mitrofanov\Irefndubna, N. Mitrofanov\Irefndubna, Y. Miyachi\Irefnyamagata, A. Moretti\Irefnntriest_utriest_i, A. Nagaytsev\Irefndubna, F. Nerling\Irefnmainz, D. Neyret\Irefnsaclay, J. Nový\Irefnnpraguectucern, W.-D. Nowak\Irefnmainz, G. Nukazuka\Irefnyamagata, A.S. Nunes\Irefnlisbon, A.G. Olshevsky\Irefndubna, I. Orlov\Irefndubna, M. Ostrick\Irefnmainz, D. Panzieri\Irefnturin_i\Arefturin_p, B. Parsamyan\Irefnnturin_uturin_i, S. Paul\Irefnmunichtu, J.-C. Peng\Irefnillinois, F. Pereira\Irefnaveiro, M. Pešek\Irefnpraguecu, M. Pešková\Irefnpraguecu, D.V. Peshekhonov\Irefndubna, N. Pierre\Irefnnmainzsaclay, S. Platchkov\Irefnsaclay, J. Pochodzalla\Irefnmainz, V.A. Polyakov\Irefnprotvino, J. Pretz\Irefnbonnpi\Arefh, M. Quaresma\Irefnlisbon, C. Quintans\Irefnlisbon, S. Ramos\Irefnlisbon\Arefa, C. Regali\Irefnfreiburg, G. Reicherz\Irefnbochum, C. Riedl\Irefnillinois, N.S. Rogacheva\Irefndubna, D.I. Ryabchikov\Irefnnprotvinomunichtu, A. Rybnikov\Irefndubna, A. Rychter\Irefnwarsawtu, R. Salac\Irefnpraguectu, V.D. Samoylenko\Irefnprotvino, A. Sandacz\Irefnwarsaw, C. Santos\Irefntriest_i, S. Sarkar\Irefncalcutta, I.A. Savin\Irefndubna\Arefo, T. Sawada\Irefntaipei, G. Sbrizzai\Irefnntriest_utriest_i, P. Schiavon\Irefnntriest_utriest_i, K. Schmidt\Irefnfreiburg\Arefc, H. Schmieden\Irefnbonnpi, K. Schönning\Irefncern\Arefi, E. Seder\Irefnsaclay, A. Selyunin\Irefndubna, L. Silva\Irefnlisbon, L. Sinha\Irefncalcutta, S. Sirtl\Irefnfreiburg, M. Slunecka\Irefndubna, J. Smolik\Irefndubna, A. Srnka\Irefnbrno, D. Steffen\Irefnncernmunichtu, M. Stolarski\Irefnlisbon, O. Subrt\Irefnncernpraguectu, M. Sulc\Irefnliberec, H. Suzuki\Irefnyamagata\Arefd, A. Szabelski\Irefnnntriest_utriest_iwarsaw T. Szameitat\Irefnfreiburg\Arefc, P. Sznajder\Irefnwarsaw, M. Tasevsky\Irefndubna, S. Tessaro\Irefntriest_i, F. Tessarotto\Irefntriest_i, A. Thiel\Irefnbonniskp, J. Tomsa\Irefnpraguecu, F. Tosello\Irefnturin_i, V. Tskhay\Irefnmoscowlpi, S. Uhl\Irefnmunichtu, B.I. Vasilishin\Irefntomsk, A. Vauth\Irefncern, J. Veloso\Irefnaveiro, A. Vidon\Irefnsaclay, M. Virius\Irefnpraguectu, S. Wallner\Irefnmunichtu, T. Weisrock\Irefnmainz, M. Wilfert\Irefnmainz, J. ter Wolbeek\Irefnfreiburg\Arefc, K. Zaremba\Irefnwarsawtu, P. Zavada\Irefndubna, M. Zavertyaev\Irefnmoscowlpi, E. Zemlyanichkina\Irefndubna\Arefo, M. Ziembicki\Irefnwarsawtu

aveiroUniversity of Aveiro, Dept. of Physics, 3810-193 Aveiro, Portugal

bochumUniversität Bochum, Institut für Experimentalphysik, 44780 Bochum, Germany\Arefsl\Arefs

bonniskpUniversität Bonn, Helmholtz-Institut für Strahlen- und Kernphysik, 53115 Bonn, Germany\Arefsl

bonnpiUniversität Bonn, Physikalisches Institut, 53115 Bonn, Germany\Arefsl

brnoInstitute of Scientific Instruments, AS CR, 61264 Brno, Czech Republic\Arefsm

calcuttaMatrivani Institute of Experimental Research & Education, Calcutta-700 030, India\Arefsn

dubnaJoint Institute for Nuclear Research, 141980 Dubna, Moscow region, Russia\Arefso

freiburgUniversität Freiburg, Physikalisches Institut, 79104 Freiburg, Germany\Arefsl\Arefs

cernCERN, 1211 Geneva 23, Switzerland

liberecTechnical University in Liberec, 46117 Liberec, Czech Republic\Arefsm

lisbonLIP, 1000-149 Lisbon, Portugal\Arefsp

mainzUniversität Mainz, Institut für Kernphysik, 55099 Mainz, Germany\Arefsl

miyazakiUniversity of Miyazaki, Miyazaki 889-2192, Japan\Arefsq

moscowlpiLebedev Physical Institute, 119991 Moscow, Russia

munichtuTechnische Universität München, Physik Dept., 85748 Garching, Germany\Arefsl\Arefr

nagoyaNagoya University, 464 Nagoya, Japan\Arefsq

praguecuCharles University in Prague, Faculty of Mathematics and Physics, 18000 Prague, Czech Republic\Arefsm

praguectuCzech Technical University in Prague, 16636 Prague, Czech Republic\Arefsm

protvinoState Scientific Center Institute for High Energy Physics of National Research Center ‘Kurchatov Institute’, 142281 Protvino, Russia

saclayIRFU, CEA, Université Paris-Saclay, 91191 Gif-sur-Yvette, France\Arefss

taipeiAcademia Sinica, Institute of Physics, Taipei 11529, Taiwan\Arefstw

telavivTel Aviv University, School of Physics and Astronomy, 69978 Tel Aviv, Israel\Arefst

triest_uUniversity of Trieste, Dept. of Physics, 34127 Trieste, Italy

triest_iTrieste Section of INFN, 34127 Trieste, Italy

turin_uUniversity of Turin, Dept. of Physics, 10125 Turin, Italy

turin_iTorino Section of INFN, 10125 Turin, Italy

tomskTomsk Polytechnic University,634050 Tomsk, Russia\Arefsnauka

illinoisUniversity of Illinois at Urbana-Champaign, Dept. of Physics, Urbana, IL 61801-3080, USA\Arefsnsf

warsawNational Centre for Nuclear Research, 00-681 Warsaw, Poland\Arefsu

warsawuUniversity of Warsaw, Faculty of Physics, 02-093 Warsaw, Poland\Arefsu

warsawtuWarsaw University of Technology, Institute of Radioelectronics, 00-665 Warsaw, Poland\Arefsu

yamagataYamagata University, Yamagata 992-8510, Japan\Arefsq {Authlist}

Corresponding authors

Deceased

aAlso at Instituto Superior Técnico, Universidade de Lisboa, Lisbon, Portugal

bAlso at Dept. of Physics, Pusan National University, Busan 609-735, Republic of Korea and at Physics Dept., Brookhaven National Laboratory, Upton, NY 11973, USA

ictpAlso at Abdus Salam ICTP, 34151 Trieste, Italy

rSupported by the DFG cluster of excellence ‘Origin and Structure of the Universe’ (www.universe-cluster.de) (Germany)

p2iSupported by the Laboratoire d’excellence P2IO (France)

dAlso at Chubu University, Kasugai, Aichi 487-8501, Japan\Arefsq

xAlso at Dept. of Physics, National Central University, 300 Jhongda Road, Jhongli 32001, Taiwan

eAlso at KEK, 1-1 Oho, Tsukuba, Ibaraki 305-0801, Japan

gAlso at Moscow Institute of Physics and Technology, Moscow Region, 141700, Russia

hPresent address: RWTH Aachen University, III. Physikalisches Institut, 52056 Aachen, Germany

yerevanAlso at Yerevan Physics Institute, Alikhanian Br. Street, Yerevan, Armenia, 0036

yAlso at Dept. of Physics, National Kaohsiung Normal University, Kaohsiung County 824, Taiwan

turin_pAlso at University of Eastern Piedmont, 15100 Alessandria, Italy

iPresent address: Uppsala University, Box 516, 75120 Uppsala, Sweden

c Supported by the DFG Research Training Group Programmes 1102 and 2044 (Germany)

l Supported by BMBF - Bundesministerium für Bildung und Forschung (Germany)

s Supported by FP7, HadronPhysics3, Grant 283286 (European Union)

m Supported by MEYS, Grant LG13031 (Czech Republic)

n Supported by SAIL (CSR) and B.Sen fund (India)

o Supported by CERN-RFBR Grant 12-02-91500

pSupported by FCT - Fundação para a Ciência e Tecnologia, COMPETE and QREN, Grants CERN/FP 116376/2010, 123600/2011 and CERN/FIS-NUC/0017/2015 (Portugal)

q Supported by MEXT and JSPS, Grants 18002006, 20540299, 18540281 and 26247032, the Daiko and Yamada Foundations (Japan)

tw Supported by the Ministry of Science and Technology (Taiwan)

t Supported by the Israel Academy of Sciences and Humanities (Israel)

naukaSupported by the Russian Federation program “Nauka” (Contract No. 0.1764.GZB.2017) (Russia)

nsf Supported by the National Science Foundation, Grant no. PHY-1506416 (USA)

u Supported by NCN, Grant 2015/18/M/ST2/00550 (Poland)

## 1 Introduction

A complete understanding of the three-dimensional parton structure of a fast moving nucleon requires the knowledge of the intrinsic motion of quarks in the plane transverse to the direction of motion, both in momentum and coordinate space. While the spatial distributions of quarks in the transverse plane are described by generalised parton distributions (GPDs), the momentum distributions of quarks in the transverse plane are described by transverse-momentum-dependent (TMD) parton distribution functions (PDFs). A precise knowledge of TMD-PDFs is found [1] to be crucial in the explanation of many single-spin effects observed in hard scattering reactions [2, 3, 4], in addition to the important role they play in spin-independent processes. In a similar way, transverse-momentum-dependent fragmentation functions (TMD-FFs) are crucial for the description of hard scattering reactions involving hadron production. Both PDFs and FFs are non-perturbative quantities that are assumed to be process-independent. The simplest examples are the spin-averaged TMD-PDF ) and the spin-averaged TMD-FF , where is the Bjorken scaling variable, is the quark intrinsic transverse momentum, is the fractional energy of the final-state hadron, and is the transverse momentum of the final-state hadron relative to the direction of the fragmenting quark. After integration over and , the TMD-PDFs and TMD-FFs reduce to the standard spin-averaged collinear PDFs and FFs, where collinear means along the direction of the virtual photon. While the knowledge on the collinear PDFs and FFs is quite advanced, very little is presently known about the dependence of TMD-PDFs and TMD-FFs on and , as only sparse experimental data are available to date.

One of the most powerful tools to assess TMD-PDFs and TMD-FFs is the semi-inclusive measurement of deep inelastic scattering (SIDIS), , where one hadron is detected in coincidence with the scattered lepton in the final state. According to the QCD factorisation theorem [5, 6] the deep inelastic scattering (DIS) process is considered to proceed via two independent sub-processes, i.e. the elementary QED process is followed by the hadronisation of the struck quark. The outgoing hadrons provide information about the original transverse motion of the quark in the nucleon via their transverse momentum vector . The latter is defined with respect to the virtual photon direction. The SIDIS cross section can be written as a convolution of a ‘hard’ scattering cross section, which is calculable in perturbative QCD (pQCD), with the non-perturbative TMD-PDFs and TMD-FFs. It depends on five kinematic variables. Two variables describe inclusive DIS, i.e. the negative square of the four-momentum transfer and the Bjorken scaling variable , where and denote the four-momenta of the virtual photon and the nucleon, respectively. Three more variables describe the final-state hadrons, i.e. the fraction of the virtual photon energy that is carried by a hadron, , the magnitude of the transverse momentum of a hadron and its azimuthal angle in the system of virtual photon and nucleon. Here, denotes the four-momentum of the hadron. In the present analysis, the dependence on is disregarded. When integrating it over, the differential cross section for spin-independent SIDIS reads as follows in twist two “TMD factorisation scheme” [7, 8]:

(1) |

with

(2) |

Here, is the lepton energy fraction that is carried by the virtual photon and is the centre-of-mass energy, which are related to and through . The hadron transverse momentum is related to and by [7]. An important consequence of the factorisation theorem is that the fragmentation function is independent of , and the parton distribution function is independent of , while both depend on .

In addition to azimuthal asymmetries in spin-independent SIDIS [9], the most relevant experimental observable to investigate spin-averaged TMD-PDFs and TMD-FFs is the differential hadron multiplicity as a function of , which is defined in Eq. 3 below. ‘Soft’ non-perturbative processes are expected to generate relatively small values of with an approximately Gaussian distribution in [10]. Hard QCD processes are expected to generate large non-Gaussian tails for 1 (GeV/). They are expected to play an important role in the interpretation of the results reported here, which reach values of up to 3 (GeV/).

Transverse-momentum-dependent distributions of charged hadrons in DIS were first measured by the EMC collaboration [11] at CERN, followed by measurements by ZEUS [12] and H1 [13, 14] at HERA. These measurements only provided data in a limited dimensional space. Only new-generation experiments provided higher statistics, thereby opening the way to analyse and present the results in several dimensions simultaneously. Recent results were obtained by several fixed-target experiments using various targets and complementary energy regimes, i.e. HERMES [15] at DESY and COMPASS [16] at CERN.

The present paper reports on a new COMPASS measurement of transverse-momentum-dependent multiplicities of charged hadrons and extends the results of our earlier publication on transverse-momentum-dependent distributions of charged hadrons [16]. The present measurement enlarges the kinematic coverage in up to 0.4 instead of 0.12, in up to 81 (GeV/) instead of 10 (GeV/) and in up to 3 (GeV/) instead of about 1 (GeV/) with significantly reduced systematic uncertainties on the normalisation of the -integrated multiplicities [16]. The data reported here represent the most precise results on differential charged hadron multiplicities available at this energy scale. This measurement is unique as its high statistics allows us to analyse the -dependence of charged-hadron multiplicities in four variables simultaneously.

The paper is organised as follows. Section 2 briefly describes the experimental apparatus. Details about the data analysis are given in Sec. 3. The measured charged-hadron multiplicities are presented and compared to previous measurements in Sec. 4. In Sec. 5 fits to the results are presented and discussed. The results are summarised in Sec. 6.

## 2 Experimental setup

The set-up of the COMPASS experiment is shortly described in this section. A more detailed description can be found in Ref. [17]. It is a fixed-target experiment, which uses the CERN Super Proton Synchrotron M2 beamline that is able to deliver high-energy hadron and muon beams. The data were collected in 2006 using a naturally polarised beam of 160 GeV/ with a momentum spread of 5%. The intensity was s with a spill length of 4.8 s and a cycle time of 16.8 s. The momentum of each incoming muon was measured before the COMPASS experiment with a precision of 0.3%. The trajectory of each incoming muon was measured before the target in a set of silicon and scintillating fibre detectors. The muons were impinging on a longitudinally polarised solid-state target located inside a superconducting magnet. The target consisted of three cells that were located along the beam one after the other. It was filled with LiD beads immersed in a liquid He/He mixture. The admixtures of H, He and Li in the target lead to an effective excess of neutrons of about . To first approximation, it can be regarded as an isoscalar deuteron target and will be referred to as such in the following. The polarisation of the middle cell (60 cm length) was opposite to that of the two outer cells (30 cm long each), and the polarisation was reversed once per day. In order to obtain spin-independent results, the target polarisation was averaged over by combining the data from all three target cells. Since the data taking in the two target polarisation states was well balanced and remaining polarisation-dependent effects are very small, this procedure ensures that for the data analysis the target can be considered as unpolarised.

The COMPASS two-stage spectrometer was designed to reconstruct scattered muons and produced hadrons in a wide range of momentum and polar angle, where the latter reaches up to 180 mrad. Particle tracking is performed by a variety of tracking detectors that are located before and after the two spectrometer magnets. The direction of the reconstructed tracks at the interaction point is determined with a precision of 0.2 mrad. The momentum resolution is 1.2% in the first spectrometer stage and 0.5% in the second one. The trigger is made by hodoscope systems supplemented by hadron calorimeters. Muons are identified downstream of hadron absorbers.

## 3 Multiplicity and data analysis

### 3.1 Multiplicity extraction

The differential multiplicity for charged hadrons, where h denotes a long-lived charged hadron (, , K, K, p or ), is defined as the ratio between the differential semi-inclusive cross section and the differential inclusive cross section

(3) |

Hadron multiplicities are measured in the four-dimensional (, , , ) space. The bin limits in the four variables are presented in Table 1.

The data used in the present analysis were collected during six weeks in 2006. The data analysis comprises event and hadron selection, the correction for radiative effects, the determination of and the correction for the kinematic and geometric acceptance of the experimental set-up as well as for detector inefficiencies, detector resolutions and bin migration, and the correction for diffractive vector-meson production. Differential hadron multiplicities are evaluated as the ratio of hadron yields in every interval of () and the number of DIS events in every interval of () corrected as described above

(4) |

Here, and denote the correction factors accounting for radiative effects in the inclusive and in the semi-inclusive case, respectively, accounts for acceptance effects, and denotes the correction factor accounting for the diffractive vector-meson contribution in the case of an inclusive (semi-inclusive) measurement. The dependence is omitted for simplicity as it enters all terms. All corrections are evaluated in the four-dimensional bins except , and , which are evaluated only in bins of and . Further kinematic dependences of upon and will be discussed in Sec. 3.2.

bin limits | |||||||||

0.003 | 0.008 | 0.013 | 0.02 | 0.032 | 0.055 | 0.1 | 0.21 | 0.4 | |

(GeV/) | 1.0 | 1.7 | 3.0 | 7.0 | 16 | 81 | |||

0.2 | 0.3 | 0.4 | 0.6 | 0.8 | |||||

(GeV/) | 0.02 | 0.04 | 0.06 | 0.08 | 0.10 | 0.12 | 0.14 | 0.17 | 0.196 |

0.23 | 0.27 | 0.30 | 0.35 | 0.40 | 0.46 | 0.52 | 0.60 | 0.68 | |

0.76 | 0.87 | 1.00 | 1.12 | 1.24 | 1.38 | 1.52 | 1.68 | 1.85 | |

2.05 | 2.35 | 2.65 | 3.00 |

### 3.2 Event and hadron selection

The present analysis uses events taken with ’inclusive triggers’, i.e. the trigger decision is based on scattered muons only. The selected events are required to have a reconstructed interaction vertex associated with an incident and a scattered muon track. This vertex has to lie inside a fiducial target volume. The incident muon energy is constrained to the range from 140 GeV to 180 GeV. In addition to the kinematic constraints given by the spectrometer acceptance, the selected events are required to have 1 (GeV/c) and GeV/. These requirements select the DIS regime and exclude the nucleon resonance region. The relative virtual-photon energy is constrained to the range to exclude kinematic regions where the momentum resolution degrades and radiative effects are most pronounced [18]. In the range , the total number of inclusive DIS events is , which corresponds to an integrated luminosity of 0.54 fb. The (, ) distribution of this selected ‘DIS sample’ is shown in Fig. 1, where a strong correlation between and is observed as expected in fixed-target experiments.

For a selected DIS event, all reconstructed tracks associated with the primary interaction vertex are considered. Hadron tracks must be detected in detectors located before and after the magnet in the first stage of the spectrometer. The fraction of the virtual-photon energy transferred to a final-state hadron is constrained to . The lower limit excludes the target fragmentation region, while the upper one removes muons wrongly identified as hadrons and excludes the region with larger contributions from diffractive production. This selection yields the ‘hadron sample’ with a total of and positively and negatively charged hadrons, respectively.

The corrections for QED higher-order effects are applied on an event-by-event basis taking into account the target composition. They are computed as a function of and according to the scheme described in Ref. [19]. For the hadron yields, the correction is calculated by excluding the elastic and quasi-elastic tails. The correction factors and are evaluated in bins of and . They are found to be smaller than 12 % for and are smaller than 5% elsewhere. An attempt to evaluate the smearing due to radiative effects as a function of and was performed using a Monte Carlo (MC) simulation, where radiative effects were simulated using the RADGEN generator [20]. A possible impact on the (, ) dependence of the results due to radiative effects is accounted for in the systematic uncertainties of the -dependence of the multiplicities.

### 3.3 Acceptance correction

The hadron multiplicities must be corrected for geometric and kinematic acceptances of the experimental set-up as well as for detector inefficiencies and resolutions, and for bin migration. The correction for a possible misidentification of electrons as hadrons is included in the acceptance correction. The full correction factor is evaluated using a MC simulation of the muon-deuteron deep inelastic scattering processes. Events are generated using the LEPTO generator [21], where the parton hadronisation mechanism is simulated using the JETSET package [22] with the tuning from Ref. [23]. Secondary hadron interactions are simulated using the FLUKA package [24]. The experimental set-up is simulated using the GEANT3 toolkits [25] and the MC data are reconstructed using the same software that was used for the experimental data [17]. The kinematic distributions of the experimental data are fairly well reproduced by the MC simulation.

In order to minimise a possible dependence on the physics generator used in the simulation and to exclude kinematic regions with large acceptance corrections, a four-dimensional evaluation of the acceptance correction factor is performed in narrow kinematic bins. In each (, ) kinematic bin, where r denotes the reconstructed values of the variables, the acceptance correction is calculated as the ratio of reconstructed () and generated () hadron yields, where both are evaluated using the simulated DIS sample after reconstruction

(5) |

An advantage of this definition is that the correction for muon acceptance cancels as it enters both numerator and denominator. The generated values of kinematic variables are used for the generated particles and the reconstructed values of kinematic variables are used for the reconstructed particles. All reconstructed MC events and particles are subject to the same kinematic and geometric selection criteria as the data, while the generated ones are subject to kinematic requirements only. The acceptance correction factor exhibits an almost flat behaviour as a function of and in most (, ) bins, except at high for (GeV/), where it remains larger than 0.4. Elsewhere, its average value is above and close to 0.6 for (GeV/) and is less than or equal to 0.6 for (GeV/). As an example, Fig. 2 shows the acceptance as a function of for positively charged hadrons. The two panels show the two bins between 0.4 and 0.8, with two bins in (, ) in each case. The acceptance correction factors for positively and negatively charged hadrons are found to be very similar, with differences on the level of 0.02-0.04.

### 3.4 Diffractive vector meson contribution

The final-state hadron(s) selected as described above may also originate from diffractive production of vector mesons (, , ) that decay into lighter hadrons (, K, p) [4, 26, 27]. This process, which can be described by the fluctuation of the virtual photon into a vector meson that subsequently interacts diffractively with the nucleon through multiple gluon exchange, is different from the interaction of the virtual photon with a single quark in the DIS process. The fraction of selected final-state hadrons originating from diffractive vector-meson decays and their contribution to the SIDIS yields are estimated in each kinematic bin using two Monte Carlo simulations. The first one uses the LEPTO generator to simulate SIDIS events, and the other one uses the HEPGEN generator [28] to simulate diffractively produced and events. Further channels, which are characterised by smaller cross sections, are not taken into account. Events with diffractive dissociation of the target nucleon represent about 25% of those with the nucleon staying intact and are also simulated. The simulation of these events includes nuclear effects, i.e. coherent production and nuclear absorption as described in Ref. [28]. The contribution of pions originating from decay to the hadron sample increases with , and reaches up to % for close to 1. For kaons, the contribution from decay is concentrated in the range , where it reaches up to 15%. The correction factors are separately evaluated for the DIS sample and the hadron sample

(6) |

(7) |

Here, denotes the fraction of diffractively produced vector-mesons present in the DIS sample, while and denote the fraction of and decay products in the hadron sample, respectively. The fraction of pions, kaons, and protons in the latter sample, which is denoted by , amounts to about 75%, 20% and 5%. The fractions () and ( and , ) are evaluated as functions of , , and .

In the following, the general behaviour of some of the above discussed correction factors is illustrated. The correction factor to account for diffractive production in the DIS yield is shown in Fig. 3 as a function of in the five bins. It reaches a maximum value of about 4% in the lowest bin. The correction factor for the contribution of diffractively produced mesons to the pion sample, (), is shown in Fig. 4 as a function of in the four bins for the lowest bin, where it has the largest value. It reaches a maximum value of about 25% for (GeV/) in the highest bins, i.e. , and decreases to few percent at small . The correction factor for the contribution of diffractively produced mesons to the kaon sample, (), is shown in Fig. 4. In this case, the maximum correction of about 35% is reached at very small in the middle bin, i.e. .

### 3.5 Systematic uncertainties

The dominant contributions to the systematic uncertainties originate from the uncertainties on the determination of the acceptance correction factor and those of the diffractive vector-meson contribution. The uncertainty on the acceptance calculation is evaluated by varying in the MC simulation both the PDF set and the JETSET parameters describing the hadronisation mechanism. The acceptance correction is estimated for each MC sample and the largest deviation with respect to the values obtained using the MC simulation described in section 3.3 is quoted as a systematic uncertainty. The validity of the correction for the electron contamination is confirmed by comparing the simulated and measured electron distributions for momenta below 8 GeV/, where electrons are identified using the RICH detector. In order to check a possible dependence on the target cell, in which the event vertex is located, the multiplicities are independently measured from the three target cells. Results from upstream and downstream target cells agree within 2-3%, while the agreement is better than 1% with the middle target cell. These differences are well covered in the acceptance correction uncertainty. A total uncertainty of 5% is estimated for the multiplicities.

The cross section for exclusive production of calculated in HEPGEN is normalised to the phenomenological model of Ref. [29]. The theoretical uncertainty on the predicted cross section in a kinematic region close to COMPASS kinematics amounts to about 30%. This results in an uncertainty on the diffractive vector-meson correction factor, which amounts up to mainly at small values of , and , and large values of .

Nuclear effects may be caused by the presence of He/He and Li in the target. The EMC Collaboration has studied in detail such nuclear effects in a similar kinematic range using carbon, copper and tin targets [11]. A -dependent decrease of 5% was observed for the multiplicities obtained using copper compared to the ones obtained using deuterium. While the effect was larger for tin, no such effect was found for carbon, so that possible nuclear effects in the present experiment are expected to be very small and are hence neglected. When comparing the results obtained from the data taken in six different weeks, no difference is observed.

All contributions to the systematic uncertainties are added in quadrature and yield a total systematic uncertainty of , except at large and at large ( (GeV/)) where it reaches about 10%. The total systematic uncertainties are shown as bands in Figs. 5 to 8. Systematic uncertainties in other figures are not shown.

## 4 Measured hadron multiplicities and comparison with other experiments

### 4.1 Results

The measured multiplicities of charged hadrons are presented in the four bins ranging from to in Figs. 58 as a function of the hadron transverse momentum in bins of and . Error bars showing the statistical uncertainties on the points are too small to be visible. The systematic uncertainties are given as bands at the bottom. All multiplicities presented in the following figures are corrected for diffractive vector-meson production. The results amount to a total of 4918 experimental data points. Their numerical values are available on HepData [30] with and without correction for diffractive vector-meson production. It should be noted that a few (, ) kinematic bins are discarded in the lowest (Fig. 5) and the highest (Fig. 8) bins of because of low statistical precision as well as large acceptance correction factors (Sec. 3.3). The average values of and in the various kinematic bins are evaluated using the DIS sample and are given in Table 2. The results obtained by integrating the multiplicities presented here over are in very good agreement with those of Ref. [26], where the multiplicities of charged pions are measured as a function of in a restricted momentum range based on an independent analysis of the same data.

(GeV/) bins | 1-1.7 | 1.7-3 | 3-7 | 7-16 | 16-81 | 1-1.7 | 1.7-3 | 3-7 | 7-16 | 16-81 |

bins | ||||||||||

0.003 - 0.008 |
0.0062 | 0.0074 | 1.25 | 1.82 | ||||||

0.008 - 0.013 | 0.010 | 0.011 | 1.32 | 2.12 | ||||||

0.013 - 0.02 | 0.016 | 0.016 | 0.017 | 1.30 | 2.28 | 3.51 | ||||

0.02 - 0.032 | 0.025 | 0.025 | 0.026 | 1.29 | 2.29 | 4.10 | ||||

0.032-0.055 | 0.039 | 0.042 | 0.043 | 0.045 | 1.39 | 2.29 | 4.52 | 8.33 | ||

0.055-0.1 | 0.068 | 0.075 | 0.077 | 2.47 | 4.65 | 9.30 | ||||

0.1-0.21 | 0.133 | 0.149 | 0.157 | 5.29 | 9.78 | 19.9 | ||||

0.21-0.4 | 0.254 | 0.291 | 11.04 | 22.1 |

Multiplicities are larger for positively than for negatively charged hadrons. This difference significantly increases as increases and shows a weak variation with . It is observed to also depend on and it increases in the range of large , i.e. , which confirms the observations made in Ref. [26]. Besides their magnitude, the dependence of the multiplicities shows a significant variation with at fixed (as well as with at fixed ) for any interval of . These observations are separately illustrated in Figs. 9 and 10 and discussed in detail in the following.

The comparison between the multiplicities of positively and negatively charged hadron is illustrated as a function of and in Fig. 9 for . On the top row, and multiplicities are presented at (GeV/) in the smallest and the largest bins with average values and , respectively. In the right column, and multiplicities are similarly presented at in the smallest and the largest bins with average values (GeV/) and (GeV/), respectively. At fixed , the ratio of to multiplicities ranges from about 1 in the first bin to about 1.3 in the last bin. This increase as a function of confirms the expectation from valence u-quark dominance, i.e. the dominance of scattering off -quarks. At fixed , the ratio of to multiplicities decreases from 1.3 in the first bin to about 1.2 in the last bin. While no significant difference is observed in the -dependence of and multiplicities, the -dependence of the multiplicities is observed to flatten at large values of , where contributions from higher-order QCD processes like QCD Compton and photon-gluon fusion (PGF) are expected to dominate. The data suggest that flattening occurs both as increases (at fixed ) and when decreases (at fixed ).

In Figure 10, the comparison between and multiplicities is illustrated as a function of . The multiplicities are presented as a function of in the four intervals in a given (, ) bin with average values and (GeV/). A high bin is chosen, where the difference in the magnitude of the multiplicities is most recognisable. The ratio of to ranges from about 1.1 in the first bin to about 2 in the last bin, reflecting the fact that part of the negative hadrons (K and ) can not be produced by the favoured fragmentation of a nucleon valence quark, which enhances the expected flavour dependence of TMD-FFs. Another feature of the data is the variation of the -dependence with increasing for both small and large values of . In particular, the data show a tendency to flatten at large as decreases, which emphases a significant -dependence of the hadron transverse momentum with respect to the transverse momentum of the fragmenting quark, .

Another intriguing effect is observed in the kinematic domain 1 (GeV/) (GeV/) and , in the range of small . Charged hadron multiplicities do not exhibit an exponential form in in this kinematic region and show an unexpected flat dependence at very small values of . This effect is also present in the earlier published [16] distributions of charged hadrons as a function of . It is illustrated in Fig. 11, which shows the multiplicity of positive hadrons as a function of up to 0.8 (GeV/) at (GeV/) and (left-hand side) and at (GeV/) and (right-hand side). It should be noted that this particular kinematic region suffers from the highest contribution of the decay products to the charged hadron sample (Fig. 4, blue curve) evaluated using the MC simulation. This effect is further discussed in Sec. 5.

The multiplicities shown in Figs 5-11 agree with the previous measurement of hadron distributions performed by COMPASS [16]. However, as mentioned in Sec. 1, this measurement considerably extends the kinematic range and reduces the statistical and systematic uncertainties, in particular the uncertainties on the normalisation of the -integrated multiplicities.

### 4.2 Comparison with other measurements

The multiplicities presented above are compared in Figs. 12-14 to results from previous semi-inclusive measurements in similar kinematic regions. The experiments are compared in Table 3.

EMC [11] | HERMES [15] | JLAB [31] | COMPASS [16] | COMPASS (This paper) | |

Target | p/d | p/d | d | d | d |

Beam energy (GeV) | 100-280 | 27.6 | 5.479 | 160 | 160 |

Hadron type | ,K | ||||

Observable | |||||

(GeV/) | 2/3/4/5 | 1 | 2 | 1 | 1 |

(GeV/) | - | 10 | 4 | 25 | 25 |

range | [0.2,0.8] | [0.1,0.85] | [0.1,0.9] | [0.1,0.9] | [0.1,0.9] |

range | [0.01,1] | [0.023,0.6] | [0.2,0.6] | [0.004,0.12] | [0.003,0.4] |

range (GeV/) | [0.081, 15.8] | [0.0047,0.9] | [0.004,0.196] | [0.02,0.72] | [0.02,3] |

In order to compare the present COMPASS results on TMD hadron multiplicities with the corresponding ones by EMC [11], our data sample is reanalysed in bins of and according to the binning given in Ref. [11]. The EMC measurements are performed in slightly different kinematic ranges in and , as shown in Tab. 3. While for the measurement described in this paper a deuteron target was used, EMC used proton and deuteron targets and also four different beam energies, which led to four different kinematic ranges. The comparison shown in Fig. 12, where the sum of and multiplicities is presented as a function of in four bins in the range , demonstrates good agreement between COMPASS and EMC results. According to the study in Ref. [10], the -dependence of the EMC data could be explained in the simple collinear parton model up to 8 (GeV/) in .

In Figure 13, the multiplicities of positively charged hadrons are compared in the four bins of to the multiplicities of positively charged pions measured by the HERMES Collaboration [15], where both were corrected for diffractive vector-meson contribution. The measurements by HERMES cover the kinematic range (GeV/) and . For this comparison, the COMPASS h multiplicities are integrated over in the closest possible range and also over . It should be noted that the two experiments cover different ranges in . While the highest value reached by HERMES is (GeV/), COMPASS reaches (GeV/). Despite this difference, a reasonable agreement in the magnitude of the measured multiplicities is found for and small . Most likely due to the differences in kinematic coverage, the agreement between the two sets is rather modest, and the data sets exhibit different dependences upon . In addition, a dip is observed in the HERMES data at very small transverse momenta, i.e. (GeV/). This dip, which is not observed in the shown -integrated distribution, appears to be very similar to the trend shown in Fig. 11 by the COMPASS data at low .

In Figure 14, the h multiplicities are compared to the semi-inclusive cross section measured by the E00-18 experiment [31] at Jefferson Lab. The measurement by the E00-18 was performed at and in the range 2 (GeV/) 4 (GeV/). The COMPASS results are given at similar (, ) values, i.e. , , and span the range 7 (GeV/ 16 (GeV/. Similar to the case of the comparison of COMPASS and HERMES data shown in Fig. 13, here the observed different -dependence could be due to the different values of the two measurements.

## 5 Fits of the measured hadron multiplicities

### 5.1 The range of small

The -dependence of the cross section for semi-inclusive measurements of hadron leptoproduction was empirically reasonably well described by a Gaussian parameterisation for the - and - dependence of TMD-PDFs and TMD-FFs in the range of small , i.e. (GeV/). This Gaussian parameterisation leads to a -dependence of the multiplicities of the form

(8) |

where the normalisation coefficient and the average transverse momentum , i.e. the absolute value of the inverse slope of the exponent in Eq. 8, are functions of , and .

A fairly good description of SIDIS [1] data was reached with the Gaussian parameterisation without considering neither the nor the quark flavour dependence of TMD-FFs. Recent semi-inclusive measurements of transverse-momentum-dependent hadron multiplicities [15] and distributions [16] aimed at an extraction of both and . These two observables, however, were found to be too strongly anti-correlated to be disentangled [8, 32, 33]. In order to extract them, a combined analysis of both the differential transverse-momentum-dependent hadron multiplicities and the spin-independent azimuthal asymmetries in SIDIS may be required. In the following we will discuss separately fits in the region of small and in the full range of accessible by COMPASS, i.e. 0.02 (GeV/ 3 (GeV/.

The hadron multiplicities presented in Figs. 5-8 are fitted in each (, , ) kinematic bin in the range 0.02 (GeV/) (GeV/) using the single-exponential function given in Eq. 8. Using only statistical uncertainties in the fit, reasonable values of per degree of freedom () are obtained in all (, ) bins, except for low values of and small values of , i.e. , where the values are significantly larger than 3 in most of the bins. Including the systematic uncertainties in the fit by adding them in quadrature to the statistical ones significantly improves the values of , whereas the fitted parameters remain unchanged. The -dependence of obtained from the fits is shown in Fig. 15 for in the five bins available in a given bin. A non-linear dependence of on is observed in the range of small and , in contrast to the range of large and where it becomes linear. In addition, significantly increases with at fixed and , especially at high . The h multiplicities have larger values of than the h ones at large , while no significant difference is observed at small . This conclusion confirms the one made in our previous publication [16], where a detailed study of the kinematic dependence of was presented and discussed.

As mentioned earlier in Sec. 4.2, the kinematic region of small and large , i.e. (GeV/) and , shows an intriguing effect in the range of small . As can be seen from Fig. 11, in this range and multiplicities do not exhibit an exponential form in and show an unexpected flat dependence at very small values of . Figure 16 shows the multiplicity of positively charged hadrons as a function of up to (GeV/) at (GeV/) and . While a single-exponential function reasonably describes the -dependence for , the experimental data clearly deviate from this functional form as increases, with values increasing from 1.8 in the smallest bin to 4.6 in the largest one. As an example, Fig. 16 shows h multiplicities at larger , i.e. (GeV/) and , where the single-exponential function fits the data well in all bins.

The measured charged-hadron multiplicities show that in the range of small , i.e. for (GeV/, the simple parameterisation using a single-exponential function describes the -dependence of the results quite well for not too large values of . For increasing , the -dependence of the multiplicities changes as can be seen in Fig. 9. A more complex parameterisation appears to be necessary to fit the data, as shown in Ref. [35] .

### 5.2 The full measured range

Up to now, only one study [10] was performed to describe the full range in using a Gaussian parameterisation for the - and - dependence of TMD-PDFs and TMD-FFs in the range GeV/, and calculating pQCD higher order collinear contributions in the range GeV/. A reasonable description of semi-inclusive hadron multiplicities and cross sections measured by the EMC [11] and ZEUS [12] Collaborations, respectively, was achieved. Below, we attempt to describe the -dependence of the above presented charged-hadron multiplicities over the full -range explored by COMPASS, i.e. 0.02 (GeV/ (GeV/, using the following two parameterisations

(9) |

(10) |

The first function () is defined as the sum of two single-exponential functions (Eq. 9). While and denote the normalisation coefficients, and denote the inverse slope coefficients of the first and the second exponential function, respectively. All coefficients depend on , and . Figure 17 shows in a typical (, , ) bin the multiplicities of positively charged hadrons as a function of fitted using . As described above for Ref. [10], the two exponential functions in our parameterisation can be attributed to two completely different underlying physics mechanisms that overlap in the region (GeV/). Figure 18 shows, as an example, multiplicities of positively charged hadrons as a function of , measured at (GeV/) for two bins of with average values and , in the four bins. Only statistical uncertainties are shown and used in the fit. Values of of about 1 are obtained in all (, , ) bins, except for a few (6 out of 81) bins, where where values as small as 0.52 and as large as 2.52 are obtained. The normalisation coefficients and are found to have a strong variation with and and a rather weak variation with , reflecting the (, ) dependence of collinear PDFs and the -dependence of collinear FFs. The inverse slope has an average value of about 0.23 (GeV/) for (GeV/) and about 0.28 (GeV/) for larger values of . Its dependence on is discussed below using Fig. 19. The inverse slope has an average value of about 0.6 (GeV/) and shows a rather weak variation with and .

The so-called Tsallis function [36], see Eq. 10, describes the two different kinds of power-law behaviour in the two regions of through a single function. The advantage of this function is that it provides both the inverse slope parameter that characterises the small- range and the exponent that parameterises the power-law tail at large . The charged-hadron multiplicities (Figs. 5-8) are fitted in each (, , ) bin using only statistical uncertainties. Reasonable values of are obtained in most bins except for 11 out of 81 bins where they are larger than 2 reaching up to . The exponent parameter has an average value of about 1.2. The exponent strongly depends on with a weaker dependence on and no variation with , while its -dependence is observed to increase with . The inverse slope parameter ranges between 0.15 (GeV/) and 0.4 (GeV/) and shows a significant non-linear dependence on over the full -range.

The inverse slopes , and , which were obtained using the fitting functions given in equations 8, 9 and 10 respectively, are presented and compared to each other as a function of in bins in Fig. 19. A weak non-linear dependence on is observed at small and , which becomes more pronounced at larger values of . The inverse slope T reproduces the same -dependence as that of described in Fig. 8. It is observed to be in fair agreement with except for .

A comparison between the data and the fit function is shown in Fig. 20 in a typical kinematic bin with (GeV/) and . The upper panel shows the multiplicities of positive hadrons as a function of and the corresponding fit function, and the lower panel shows the ratio between the data and the fit. A comparison between the two fitting functions and is shown in Fig. 20 for the same (, , ) bin. The -dependence of multiplicities is equally well described by the two functions and , as can be seen from the ratio in Fig. 20. The same agreement is obtained for negatively charged hadrons.

## 6 Summary

We have measured differential multiplicities of charge-separated hadrons in semi-inclusive measurements using muons of 160 GeV/ impinging on an isoscalar (deuteron) target. Using a high-statistics data set collected in 2006, the measurement covers a wide kinematic domain of (GeV/), GeV/, , and 0.02 (GeV/) 3 (GeV/). The results are presented as a function of the square of the hadron transverse momentum in three-dimensional bins of , and , which leads to a total of 4918 experimental data points. The numerical values are available on HepData [30] with and without subtraction of the estimated contribution of diffractive vector-meson production in SIDIS.

The h multiplicities are only slightly larger than the h ones in most of the bins, while for large and this difference increases. No significant difference between h and h is observed in the shape of the -dependence of the multiplicity. Both h and h multiplicities are observed to flatten at very small values of in the kinematic region of low and and large , where contributions from diffractive vector-meson production are the highest. Our results are compared to earlier measurements of hadron multiplicities and cross sections by EMC, HERMES and JLab. Good agreement was found with EMC for (GeV/), although the EMC data were collected at different beam energies and with different targets. In order to compare with HERMES, we have integrated our multiplicities over the phase space that is common to both experiments. While reasonable agreement is obtained at small and , differences are observed for large and where neither magnitudes nor -dependences agree. The semi-inclusive cross section measured by the E00-18 experiment at JLab shows fair agreement with COMPASS h multiplicities, albeit with some discrepancy in the -dependence that might be explained by the difference in the kinematic ranges of the measurements.

In the range of small-, i.e. (GeV/, the measured multiplicities were successfully fitted using a single Gaussian parameterisation. A non-linear -dependence of the average transverse momentum is observed in the range of small and , which confirms the conclusions of Ref. [16], while it is almost linear for large values of and . In order to fit the multiplicities over the full -range measured by COMPASS, a more complex functional form is required, i.e. either a sum of two Gaussian functions or the so-called Tsallis function. All fits reproduce the data well and their inverse slopes agree well with one another already when using only statistical uncertainties in the fits.

## 7 Acknowledgements

We gratefully acknowledge the support of the CERN management and staff and the skill and effort of the technicians of our collaborating institutes. This work was made possible by the financial support of our funding agencies.

## References

- [1] M. Anselmino, M. Boglione, U. D’Alesio, A. Kotzinian, F. Murgia and A. Prokudin, Phys. Rev. D 71 (2005) 074006.
- [2] D. L. Adams et al. [E581 and E704 Collaborations], Phys. Lett. B 261 (1991) 201.
- [3] J. Adams et al. [STAR Collaboration], Phys. Rev. Lett. 92 (2004) 171801.
- [4] A. Airapetian et al. [HERMES Collaboration], Phys. Rev. D 64 (2001) 097101.
- [5] J. C. Collins, L. Frankfurt and M. Strikman, Phys. Rev. D 56 (1997) 2982.
- [6] X. d. Ji, J. P. Ma and F. Yuan, Phys. Lett. B 597 (2004) 299.
- [7] A. Bacchetta, M. Diehl, K. Goeke, A. Metz, P. J. Mulders and M. Schlegel, JHEP 0702 (2007) 093.
- [8] M. Anselmino, M. Boglione, J. O. Gonzalez Hernandez, S. Melis and A. Prokudin, JHEP 1404 (2014) 005.
- [9] C. Adolph et al. [COMPASS Collaboration], Nucl. Phys. B 886 (2014) 1046.
- [10] M. Anselmino, M. Boglione, A. Prokudin and C. Turk, Eur. Phys. J. A 31 (2007) 373.
- [11] J. Ashman et al. [European Muon Collaboration], Z. Phys. C 52 (1991) 361.
- [12] M. Derrick et al. [ZEUS Collaboration], Z. Phys. C 70 (1996) 1.
- [13] C. Adloff et al. [H1 Collaboration], Nucl. Phys. B 485 (1997) 3.
- [14] F. D. Aaron et al. [H1 Collaboration], Phys. Lett. B 665 (2008) 139.
- [15] A. Airapetian et al. [HERMES Collaboration], Phys. Rev. D 87 (2013) 074029.
- [16] C. Adolph et al. [COMPASS Collaboration], Eur. Phys. J. C 73 (2013) 2531; Erratum: [Eur. Phys. J. C 75 (2015) 94]
- [17] P. Abbon et al. [COMPASS Collaboration], Nucl. Instrum. Meth. A 577 (2007) 455.
- [18] N. Makke, CERN-THESIS-2011-279.
- [19] A. A. Akhundov, D. Y. Bardin, L. Kalinovskaya and T. Riemann, Fortsch. Phys. 44 (1996) 373.
- [20] I. Akushevich, H. Bottcher and D. Ryckbosch, In *Hamburg 1998/1999, Monte Carlo generators for HERA physics* 554-565 [hep-ph/9906408].
- [21] G. Ingelman, A. Edin and J. Rathsman, Comput. Phys. Commun. 101 (1997) 108.
- [22] T. Sjostrand, hep-ph/9508391.
- [23] C. Adolph et al. [COMPASS Collaboration], Phys. Lett. B 718 (2013) 922.
- [24] A. Ferrari, P. R. Sala, A. Fasso and J. Ranft, CERN-2005-010, SLAC-R-773, INFN-TC-05-11.
- [25] R. Brun, M. Caillat, M. Maire, G. N. Patrick and L. Urban, CERN-DD/85/1.
- [26] C. Adolph et al. [COMPASS Collaboration], Phys. Lett. B 764 (2017) 1.
- [27] C. Adolph et al. [COMPASS Collaboration], Phys. Lett. B 767 (2017) 133.
- [28] A. Sandacz and P. Sznajder, arXiv:1207.0333 [hep-ph].
- [29] S. V. Goloskokov and P. Kroll, Eur. Phys. J. C 53 (2008) 367.
- [30] The Durham HepData Project, url = ”http://durpdf.dur.ac.uk/”.
- [31] R. Asaturyan et al., Phys. Rev. C 85 (2012) 015202.
- [32] A. Signori, A. Bacchetta, M. Radici and G. Schnell, JHEP 1311 (2013) 194.
- [33] M. Radici, talk at QCD evolution 2014, https://www.jlab.org/conferences/qcd2014/.
- [34] A. D. Martin, W. J. Stirling, R. S. Thorne and G. Watt, Eur. Phys. J. C 63 (2009) 189.
- [35] C. A. Aidala, B. Field, L. P. Gamberg and T. C. Rogers, Phys. Rev. D 89 (2014) 094002.
- [36] G. Wilk and Z. Wlodarczyk, Eur. Phys. J. A 40 (2009) 299.