Study of charged hadron multiplicities in charged-current neutrino-lead interactions in the OPERA detector

The OPERA experiment was designed to search for $\nu_{\mu} \rightarrow \nu_{\tau}$ oscillations in appearance mode through the direct observation of tau neutrinos in the CNGS neutrino beam. In this paper, we report a study of the multiplicity of charged particles produced in charged-current neutrino interactions in lead. We present charged hadron average multiplicities, their dispersion and investigate the KNO scaling in different kinematical regions. The results are presented in detail in the form of tables that can be used in the validation of Monte Carlo generators of neutrino-lead interactions.

Abstract The OPERA experiment was designed to search for ν µ → ν τ oscillations in appearance mode through the direct observation of tau neutrinos in the CNGS neutrino beam. In this paper, we report a study of the multiplicity of charged particles produced in charged-current neutrino interactions in lead. We present charged hadron average multiplicities, their dispersion and investigate the KNO scaling in different kinematical regions. The results are presented in detail in the form of tables that can be used in the validation of Monte Carlo generators of neutrino-lead interactions.

Introduction
The multiplicity distribution of charged hadrons is an important characteristic of the hadronic final states in hard scattering processes. Since it reflects the dynamics of the interaction, it has been extensively studied in cosmic rays, fixed target and collider experiments [1][2][3][4][5]. These data are useful to improve models of particle production which are used in Monte Carlo (MC) event generators.
In this paper, we report the result on charged hadron production initiated in charged-current ν µ interactions in the OPERA target. The basic unit that constitutes the target is the Emulsion Cloud Chamber (ECC) detector which is a stack of nuclear emulsion films acting as high precision trackers interleaved with lead plates that provide a massive target. The excellent spatial resolution of nuclear emulsion allows the determination of the event topology and the measurement of charged particle trajectories. Therefore, it is well suited for the investigation of the moments of the charged particles multiplicity distribution. However, only few studies of charged particle multiplicity in neutrino-nucleon interactions were made using the nuclear emulsion technology [5][6][7].
In the following, a short description of the experimental setup and of the procedure used to locate neutrino interactions in the target is given; the data sample and the analysis procedure are presented. Then, multiplicity moments and investigation of KNO scaling [8] in different kinematical regions are presented in a form suitable for the validation of MC generators of neutrino−lead interactions.

Experimental Procedure
The OPERA experiment has been designed for the observation of ν µ → ν τ oscillations in appearance mode in the CNGS (CERN Neutrinos to Gran Sasso) neutrino beam. The detector was located at the underground Gran Sasso Laboratory (LNGS) in Italy. It was exposed to the CNGS muon Corresponding authors. E-mail: ali.murat.guler@cern.ch, Cagin.Gunes@ankara.edu.tr neutrino beam with mean energy 17 GeV at a distance of 730 km. OPERA reported the discovery of ν τ appearance with a significance of 5.1σ [10]. The OPERA detector is a hybrid setup consisting of electronic detectors and a massive lead-emulsion target segmented into ECC units, called bricks. The detector is composed of two identical Super-Modules, each of which has a target section followed by a muon spectrometer which is composed of a dipole magnet instrumented with resistive plate chambers and drift tubes. Each target section has 31 brick walls interleaved with orthogonal pairs of scintillator strip planes that compose the Target Tracker (TT). A detailed description of the OPERA detector can be found in [11]. TT hit patterns are used to identify the bricks possibly containing the neutrino interaction vertex [12]. The most probable brick is then extracted from the target and its CS films scanned. If a signal compatible with the TT predictions is found, the brick is disassembled and its films analysed. Once the vertex has been located, a surrounding volume of about 2 cm 3 is scanned to determine the event topology. Otherwise, the procedure is repeated in the next brick in the probability ranking. Track recognition in an emulsion layer is based on 16 tomographic images taken by the sensor of an automated microscope and equally spaced through the 44-µm depth of the layer. A sequence of aligned grains in a layer forms a micro-track ( Figure 1a) and the association of two matching micro-tracks on each side of the plastic base in a film constitutes a base-track ( Figure 1b). Track positions and slopes are determined by a linear fit through base-tracks in the analysed volume. The details of the event analysis procedure are described in [13].

Analysis
During the physics runs between 2008 and 2012, OPERA collected data corresponding to 1.8x10 20 protons on target. The electronic detectors recorded 19 505 neutrino interactions in the target fiducial volume. The search of the neutrino vertex in the first and second most probable bricks plus some additional selections (see [13] for more details) resulted in a sample of 5603 located events out of which 4406 have an identified muon. For the present measurement an unbiased sub-sample of 818 events occurring in the lead with a negatively charged muon identified by the muon spectrometer was selected in order to measure the track and vertex parameters in the target including a detailed check of the nuclear break-up and evaporation processes.
It is further required that W 2 , the square of the invariant mass of the hadronic system measured with the electronic detector, is larger than 1 GeV 2 /c 4 in order to eliminate the quasi-elastic contribution. The final data sample contains 795 events with an identified muon. The contamination from ν µ neutral-current interactions in the final data sample is estimated by MC simulations to be less than 1%.
Selected ν µ CC events are inspected carefully and tracks are classified as being left by a minimum ionisation particle (mip), grey and black depending on their ionisation features. The mip tracks are left by the muon and by highly relativistic charged hadrons resulting from the cascade of interactions generated inside the target nucleus by the primary hadrons emitted with the muon at the neutrino-nucleon interaction. Black tracks are produced by low energy fragments (protons, deuterons, alpha particles and heavier fragments) emitted from the excited target nucleus. Black tracks are classified as backward or forward based on the emission direction. Grey tracks are left by slow particles which are interpreted as being recoil nucleons emitted during the nuclear cascade [14].
The black tracks are easy to recognise visually since they are heavily ionising, they have short path lengths and stop within two lead plates. The separation between mip and grey tracks is based on the Pulse Height Volume (PHV) [15,16] which is defined as the sum of the number of pixels associated with each track in all sixteen layers of images [17]. PHV indicates the track width and is therefore a measure of the grain density of a track that reflects the energy deposition of a particle in the emulsion film. The PHV distribution of muon tracks is shown in Figure 2. All muon tracks have a PHV below 85 and mip tracks are defined by requiring a PHV smaller than 85, the tracks with a higher PHV being classified as grey. The mip, grey and black track multiplicities are shown in Figures 3, 4 and 5, respectively. The av-erage numbers of mip and grey tracks in ν µ CC events are n mip = 2.94 ± 0.05 and n g = 0.22 ± 0.01, respectively. The average number of backward and forward black tracks are measured to be n b B = 0.15 ± 0.01 and n b F = 0.38 ± 0.02, respectively. The charged hadrons multiplicity n ch is defined as the number of mip tracks excluding the muon track, its average being n ch = n mip − 1 = 1.94 ± 0.05. The distribution of mip tracks as a function of their emission angle with respect to the beam axis is given separately for muons and charged hadrons in Table 1.

Efficiency Estimation
Reconstruction and location efficiencies are computed using the standard OPERA simulation framework. The neutrino fluxes and spectra are based on the FLUKA simulation [18] of the CNGS beam-line. The neutrino interactions in the detector are generated using the NEGN generator [19]. MC-generated ν µ CC events are processed through the full OPERA simulation chain, from the event classification and brick finding provided by the electronics detectors to the CS analysis and event location and analysis in the brick. The location efficiency, shown in Table 2, of the ν µ CC events is estimated as a function of W 2 and of the charged hadron multiplicity. Since the event location is done using mip tracks, the location efficiency does not depend on the black and grey track multiplicities at the neutrino interaction vertex. Figure 6 shows the good agreement between the charged hadrons multiplicities obtained for observed and MC simulated data. The location efficiency correction is hereafter applied to the measured data distributions.  In ν µ CC interactions, neglecting the Fermi motion, the full kinematics of the event can be reconstructed from the measurement of the muon momentum p µ , its angle θ µ with respect to the beam axis, and E had , the energy transfer to the hadronic system, obtained by measuring the energy deposition in the TT (detailed information is available in [11]) where E ν and E µ are the energy of the incoming neutrino and the muon, Q 2 ν the squared four-momentum transfer, m N and m µ are the mass of the nucleon and muon respectively. The average neutrino energy of the data sample is E ν = (19.6 ± 0.3) GeV, the mean-square momentum transferred to the hadronic system is Q 2 ν = (5.7 ± 0.3) GeV 2 /c 2 , and the mean-square of the invariant mass of the hadronic system is W 2 = (16.9 ± 0.6) GeV 2 /c 4 .

Multiplicity Distributions
The average charged hadrons multiplicity as a function of W 2 is presented in Figure 7. The data is well described by a linear function in lnW 2 : The values of the fitted parameters are a = −0.19 ± 0.18 and b = 0.76 ± 0.07. Values for different W 2 bin intervals are shown in Table 3. A comparison with other experiments is given in Table 4.  Fig. 7 The average charged hadron multiplicity distributions as a function of lnW 2 .

Dispersion
One of the characteristics of the multiplicity distribution which is of considerable theoretical interest is its dispersion. In this section we investigate its dependence on the average multiplicity. The dispersion D ch is defined as D ch = n 2 ch − n ch 2 . For independent particle production it would follow a Poisson distribution with D ch = n ch . However, it was observed that charged particles production in hadronic interactions satisfies an empirical parameterisation [20]: (2) Figure 8 shows the dependence of the dispersion on the average multiplicity n ch with a linear fit superimposed. The values of the fitted parameters are A = 0.59 ± 0.12 and B = 0.46 ± 0.06. They are shown in Table 5 together with those obtained in other experiments for comparison.   Table 3 The charged hadron multiplicity n ch distribution as a function of W 2 (errors shown are statistical only). Koba, Nielsen, and Olesen have shown that the multiplicity distribution P(n ch ) scaled by the average multiplicity n ch is asymptotically independent of the primary energy and depends only on variable z = n ch n ch : n ch . P(n ch ) KNO scaling is derived from Feynman scaling, i.e. based on the assumption that the rapidity density dn ch dy reaches its limit value at y = 0 above a certain energy which corresponds to an asymptotic scaling of the total multiplicity as n ch ∝ ln √ s. KNO scaling implies that the intercept A in Eq. 2 be compatible with 0, which is not the case at low to medium energies for all kinds of interactions. Buras et al. [21] have where the reaction dependent and energy independent parameter α is chosen in order to provide an extension of the KNO scaling to low energies. In the current analysis, this implies that α = A B = 1.28. Hence (n ch − α) P(n ch ) = Ψ (z ) A tentative explanation for a non-zero value for α has been proposed in terms of a leading particle effect in interactions of hadrons [22] and of neutrinos [23] as well as resulting from the heavy nuclear targets in neutrinos experiments using emulsion [5]. Figure 9 shows the distributions obtained for Ψ as a function of z for three different intervals of W 2 for ν µ -lead CC interactions. The fitted curve superimposed on the data results from a parameterisation of Ψ (z ) Ψ (z ) = (Az + B(z ) 3 +C(z ) 5 + D(z ) 7 )e −(Ez +F(z ) 3 ) (6) The fit parameters are given in Table 6. The data shows good agreement with approximate KNO scaling. Such observation has been made by other experiments [5], [23][24][25][26].  Fig. 9 KNO scaling distribution. The curve superposed on the data results from a fit using the parameterisation defined in Eq. 6. Table 6 The best-fitted parameters of the function given by Eq. 6.

Conclusion
In this article, we studied the characteristics of the multiplicity distribution of charged hadrons in neutrino-lead interactions in the OPERA detector with the main objective to aid in tuning the models used in MC event generators. For this purpose, the results are presented in detail in the form of tables. They can be summarized as follow: i The dependence of the average multiplicity n ch on ln W 2 is approximately linear. ii The dependence on the charged hadrons multiplicity n ch of its dispersion D ch is approximately linear. iii Approximate KNO scaling is valid for the charged hadrons multiplicity.