![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0386v1-figure-95a1d45f-10fb-4bad-94e8-b250f164a3ef.jpg)

Prog. Theor. Exp. Phys. 2014 , 00000 (14 pages) DOI: 10.1093 / ptep/0000000000

## Study of hadron interactions in a lead-emulsion target

Hirokazu Ishida , Tsutomu Fukuda , Takafumi Kajiwara , Koichi Kodama , 1 1 1 2 Masahiro Komatsu , Tomokazu Matsuo , Shoji Mikado , Mitsuhiro Nakamura , 3 1 4 3 Satoru Ogawa , Andrey Sheshukov , Hiroshi Shibuya , Jun Sudou , Taira Suzuki , 1 5 1 1 1 Yusuke Tsuchida 1

- 1 Department of Physics, Toho University, Funabashi 274-8510, Japan
- 2 Aichi University of Education, Kariya 448-8542, Japan
- 3 Nagoya University, Nagoya 464-8602, Japan
- 4 Nihon University, Narashino 275-8576, Japan
- 5 JINR - Joint Institute for Nuclear Research, Dubna 141980, Russia
- ∗ E-mail: shibuya@ph.sci.toho-u.ac.jp

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

Topological and kinematical characteristics of hadron interactions have been studied using a lead-emulsion target exposed to 2, 4 and 10 GeV/ c hadron beams. A total length of 60 m π -tracks was followed using a high speed automated emulsion scanning system. A total of 318 hadron interaction vertices and their secondary charged particle tracks were reconstructed. Measurement results of interaction lengths, charged particle multiplicity, emission angles and momenta of secondary charged particles are compared with a Monte Carlo simulation and appear to be consistent. Nuclear fragments emitted from interaction vertices were also detected by a newly developed emulsion scanning system with wide-angle acceptance. Their emission angle distributions are in good agreement with the simulated distributions. Probabilities of an event being associated with at least one fragment track are found to be greater than 50% for beam momentum P &gt; 4 GeV/ c and are well reproduced by the simulation. These experimental results validate estimation of the background due to hadron interactions in the sample of τ decay candidates in the OPERA ν µ → ν τ oscillation experiment.

Subject Index

C30, C32, H14, H16

## 1. Introduction and overview

The neutrino oscillation channel ν µ → ν τ is considered to be the dominant underlying process in the atmospheric neutrino sector. To prove this, the detection, in an initially pure ν µ beam, of the τ lepton produced in a ν τ charged current interaction is essential and provides a clear signature of the oscillation. The nuclear emulsion is the only particle detector which can detect τ decays by resolving their short flight path [1, 2]. The ongoing ν µ → ν τ oscillation experiment, OPERA [3-5], employs the emulsion cloud chamber (ECC) technique to provide both a large target mass and an excellent spatial resolution. Furthermore, ν τ appearance can be used to probe new physics, such as non-standard interactions and light sterile neutrino scenarios [6, 7]. An ECC brick is composed of layers of emulsion films interleaved with heavy material plates. In the ECC brick, a secondary hadron produced in neutrino interactions could interact with the heavy materials and mimic a hadronic decay of the τ lepton, as shown in Figure 1. Hence such a secondary hadron interaction could be a source of background

for τ decays. Expected number of hadron interactions as background in the sample of τ decay candidates is evaluated with a Monte Carlo simulation, which must be validated by experimental data 1 .

A high speed automated microscope system, S-UTS [8], has been developed to analyze particle tracks in emulsion films. This system allows to follow beam particle tracks for a long distance and measure their interactions with large statistics. In addition, a new type of automated emulsion scanning system [9] has recently been developed to detect nuclear fragments emitted with large angles.

We have studied hadron interactions in an ECC brick exposed to 2, 4 and 10 GeV/ c hadron beams and have compared experimental results with those of the Monte Carlo simulation.

Fig. 1 Left: a hadronic decay of the τ lepton produced in a ν τ charged current interaction. Right: an interaction of a secondary hadron ( h ) produced in a ν µ neutral current interaction.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0386v1-figure-470bb3a3-85b8-40f3-a956-c2ffa5b5bd2e.jpg)

## 2. Detectors and beam exposure

We exposed an ECC brick to 2, 4, 10 GeV/ c secondary hadron beams at CERN PS-T7 beam line in May 2001. The brick was composed of 29 emulsion films [10] (44 µ m thick emulsion layers on both sides of a 205 µ m thick plastic base) interleaved with 28 lead plates, 1 mm thick. The brick was 12.8 cm wide, 10.2 cm high and 3.7 cm thick. It was put on a turntable and tilted in the horizontal plane by an angle with respect to the beam of ± 50 mrad as shown in Figure 2. The relative beam momentum spread, ∆ P/P , is 1% [11]. The negatively (positively) charged beams are mainly composed of pions (pions and protons) with small contamination from electrons (positrons). The exposure conditions; composition, momentum, tilt angle, spot size, mean density integrated over the spot size are summarized in Table 1.

## 3. Measurement and analysis method

## 3.1. Beam tracks

We scanned the whole area of all the emulsion films by using S-UTS. Firstly track segments, so called micro tracks, were detected on each layer of emulsion films. Positions ( x, y ) and

1 The momentum region relevant to the hadron interaction background in the OPERA experiment is around 2 GeV/ c to 10 GeV/ . c

Table 1 Summary of beam exposures. For each beam spot, its composition, its momentum P , its angle with respect to the perpendicular of the emulsion film, its spot size, and its mean density integrated over the spot size are tabulated. For each spot, a tilt of -50 mrad corresponds to the positively charged hadron beam ( π + , p ) and +50 mrad to the negatively charged hadron beam ( π -).

|   Spot | Composition   |   P [GeV/ c ] | Angle [mrad]   | Spot size φ [mm]   | Mean density [/cm   |
|--------|---------------|---------------|----------------|--------------------|---------------------|
|      1 | π -           |            10 | + 50           | ∼ 3 . 0            | 3 . 0 × 10 3        |
|      2 | π ± , p       |             2 | ∓ 50           | ∼ 10 . 0           | 14 . 0 × 10 3       |
|      3 | π ± , p       |             4 | ∓ 50           | ∼ 6 . 0            | 10 . 0 × 10 3       |
|      4 | π ± , p       |             4 | ∓ 50           | ∼ 7 . 5            | 1 . 5 × 10 3        |
|      5 | π ± , p       |             2 | ∓ 50           | ∼ 15 . 0           | 2 . 0 × 10 3        |

slopes (tan θ x , tan θ y ) of the micro tracks were measured. Base tracks were reconstructed by connecting two corresponding micro tracks (reconstructed in emulsion layers) across the plastic base. The slope acceptance of S-UTS is tan θ &lt; 0 6, where . θ is the track angle with respect to the perpendicular of the emulsion film (the z axis). The track finding efficiency, the probability to find a base track in a film, is evaluated to be 94 8 . ± 0 2% by examining . whether a base track exists or not in the middle of 5 consecutive films.

The positions ( x, y ) and slopes (tan θ x , tan θ y ) of base tracks were measured with respect to the microscope's coordinate system. After connection of base tracks among emulsion films are performed, we adjust rotation, slant, parallel translation and distance between every pair of adjacent two emulsion films, so that differences of track positions and slopes in two films are minimized. Residuals of positions and slopes are ∼ 4 µ m and ∼ 0 004 respectively. .

Good quality beam tracks are selected by requiring they have track segments found in all 3 most upstream films and their slope in the most upstream film is compatible with the average slope within 0.004 (1 σ ) where the average is calculated from all the measured beam track slopes. For further analysis to be compared with the simulation, good negatively charged pion ( π -) tracks with these conditions were selected in the area of 1 0 . × 1 0 cm . 2 from 1st,

Fig. 2 Left: schematic view of the ECC brick structure. Right: beam spot positions on the ECC brick.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0386v1-figure-70c35444-1022-49c8-a210-048612145255.jpg)

4th and 5th spots 2 , each corresponding to 10, 4 and 2 GeV/ , as shown in Figure 2. On the c other hand, positively charged hadron beams are mainly composed of π + s and protons, which cannot be differentiated in this momentum range. Relative abundance is not known precisely. Therefore positively charged hadrons are not easy to be compared with the simulation.

## 3.2. Interactions

The beam tracks in the ECC brick are composed of successive track segments and some of them are missing due to inefficiency or hadron interactions. The reconstructed tracks are followed down until no track segments found in 3 consecutive films (Figure 3, Left), then a beam particle is supposed to interact at the lead plate just downstream the film where the last track segment is found. The last track segment and its extrapolated position in the film downstream the supposed interaction position, the vertex film, are visually inspected. If the last track segment is confirmed that it is a real track and no track is found in the vertex film, the interaction is then considered as confirmed. If a track exists near the predicted position in the vertex film, we measure its position and slopes manually and the beam track is followed further downstream. For beam tracks escaping from the downstream face of the brick, a straight line fit is applied to detect a kink topology somewhere in the brick. If a kink with θ kink &gt; 20 mrad is found, it is also considered as an interaction.

## 3.3. Secondary particles

Each confirmed interaction is subject to the secondary particle track search inside the area of 3 0 . × 3 0 cm . 2 in 6 downstream films of the interaction vertex. Secondary particle tracks are sought under the following conditions, as shown in Figure 3 right.

The angular acceptance of secondary particle track search is tan θ &lt; 0 6, which is limited . by the ordinary S-UTS performance. For a secondary particle candidate found in at least 1 of the first 3 downstream films including the vertex film, we require the minimum distance MD of the candidate track with respect to the parent beam track to be less than (10 + 0 01 . × ∆ ) z µ m taking into account both the measurement error and the error due to the multiple Coulomb scattering (MCS) of the tracks, where the ∆ z ( µ m) is the z distance of the MD position from the nearest emulsion layer. The scanning data have a lot of fake tracks due to low energy Compton electrons and random noises. In order to reject such fake tracks, we demand the condition that base tracks exist in ≥ 3 out of the 6 downstream films and they can be reconstructed as a single secondary track. If base tracks matching with the secondary track candidate exist in ≥ 2 films out of the 3 upstream films, it is considered as a passing-through track which is not relevant to the interaction vertex.

Once a secondary track candidate is found, it is checked by visual inspection in 3 films: 1 upstream and 2 downstream films of the interaction vertex, in order to eliminate the possibility of a passing-through track or an electron track originated from γ conversion in a downstream lead plate. The efficiency of the secondary particle track search was evaluated to be 99.2% by applying the S-UTS track finding efficiency to the conditions in the whole procedure described above.

2 The 2nd and 3rd beam spots were not used because their mean densities were too high for the analysis.

Fig. 3 Left: detection method of beam interactions. Beam tracks are followed down until they are not found in three consecutive films. Extrapolated area (a), where a matching track is sought, becomes larger in the downstream films. Right: detection method of secondary particles. (b) A base track which satisfies the condition, MD&lt; (10 + 0 01 . × ∆ ) z µ m, must exist in at least 1 of the first 3 downstream films. (c) The base track must exist in at least 3 out of the 6 downstream films. (d) If the track exists in ≥ 2 out of 3 upstream films, it is considered as a passing-through track.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0386v1-figure-886b426a-1da6-4212-9940-dfd1c5996b97.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0386v1-figure-4c767d8a-63e4-4865-81d6-8910c447cbfc.jpg)

Momenta of secondary tracks are estimated by measuring their multiple Coulomb scattering in the brick. There are two methods in the multiple scattering measurements, the coordinate method [12] and the angular method [13]. We use the coordinate method because it is more accurate when small multiple scattering signal is expected although it requires precise alignment among the emulsion films. We define the second difference δ i in one projection by:

$$\delta _ { i } = x _ { i + 2 } - x _ { i + 1 } - \frac { x _ { i + 1 } - x _ { i } } { z _ { i + 1 } - z _ { i } } \cdot ( z _ { i + 2 } - z _ { i + 1 } ) \quad \ \ ( 1 )$$

where the base track position at i th film ( i ∈ 1 , ..., n ) is ( x , z i i ), as shown in Figure 4. A similar expression can be written in the other projection, ( y , z i i ), and data in the two projections are used as independent measurements of the second difference. The multiple scattering signal ∆ sig is related to the momentum p by [14, 15]:

$$\Delta _ { \text{sig} } = \frac { t } { 2 \sqrt { 3 } } \frac { 0. 0 1 3 6 \, \text{GeV} } { p \, \beta \, c } \cdot \sqrt { \frac { t } { X _ { 0 } } } \left \{ 1 + 0. 0 3 8 \ln \left ( \frac { t } { X _ { 0 } } \right ) \right \}$$

where βc , t , and X 0 are the particle velocity, distance traversed, and radiation length in the material respectively. On the other hand, it is measured from second differences δ i and their errors /epsilon1 i :

$$< \delta _ { i } ^ { 2 } > = \Delta _ { s i g } ^ { 2 } + < \epsilon _ { i } ^ { 2 } >$$

where &lt;&gt; means the average over the available data. Since the multiple scattering measurement is statistical estimation of the particle momentum, it needs a sufficient number of data. To obtain reliable estimation, we apply the coordinate method only when at least 14 emulsion films are available for the position measurements.

Fig. 4 Schematic view of the multiple Coulomb scattering (MCS) measurement. The second differences of cell length = 1 and 2 are defined in top and bottom figures, respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0386v1-figure-a44a1477-3960-4e89-9e43-af874796b931.jpg)

## 3.4. Nuclear fragments

If a secondary particle has a value of β &lt; 0 7, the particle is observed as a heavily ionizing . track or a nuclear fragment, as shown in Figure 5. Nuclear fragments emitted from hadron interactions have also been sought by a newly developed automatic emulsion scanning system [9]. The acceptance to the track slopes are much improved up to | tan θ | &lt; . 3 0. Since nuclear fragments are emitted almost isotropically, the new scanning system with a wider angular acceptance is suitable in order to detect the nuclear fragments with a good efficiency. The system has a large field of view, 352 × 282 µ m , and has a track finding efficiency of 2 practically 100% for nuclear fragments [9].

Searched area is defined as 3 5 . × 2 5 mm . 2 of both upstream and downstream films of each interaction vertex and its angular acceptance is | tan θ | &lt; . 3 0. Since the nuclear fragments have low energy (few tens of MeV), they are expected to suffer large multiple scattering. Therefore the loose condition, IP &lt; (100 + 0 01 . × ∆ ) z µ m is imposed. Here the 'Impact Parameter' IP is defined as the perpendicular distance between the path of a track and the vertex which is reconstructed by the beam and secondary tracks. The ∆ z ( µ m) is the distance between the vertex and the nearest emulsion layer where the nuclear fragment track is measured. A narrower condition, IP &lt; (50 + 0 01 . × ∆ ) z µ m, is imposed in the case of 10 GeV/ c interactions because the beam density is high. To remove fake tracks, each found track candidate is required to be a base track and is confirmed by visual inspection.

## 4. FLUKA simulation

A FLUKA 3 [16, 17] based Monte Carlo simulation was employed for comparison. 30,000 negatively charged beam pions each were generated for momentum of P =2, 4 and 10 GeV/ . c The generated beam pions, of which directions were spread using observed distributions in this experiment, entered the brick. Some of them interacted with the brick material according to their cross sections and produced secondary particles including nuclear fragments by

3 The FLUKA version 2011.1 was used in this study.

Fig. 5 Photograph of a hadron interaction in an emulsion layer. Nuclear fragments are observed as black or gray tracks. Association of such highly ionizing particles is evidence of a hadron interaction.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0386v1-figure-80376d2c-aa67-4932-89ec-0259bba0f2d9.jpg)

using the PEANUT model [18] of hadron interactions. Exposed density, position and slope measurement errors, detection efficiency of base tracks were adjusted to be the same as real experimental conditions.

Then the simulation analysis has been performed by the same methods and tools that are used for the real data. Secondary particles with β ≥ 0 7 ( . β &lt; 0 7) are defined as relativistic . particles (non-relativistic particles, i.e. nuclear fragments).

## 5. Results and discussion

## 5.1. Interaction length

In total, 318 interactions (77 for 2 GeV/ c , 68 for 4 GeV/ , 173 for 10 GeV/ ) have been found c c from the interaction measurements. Out of the found interactions of 2 GeV/ c (10 GeV/ ), c one (three) took place in the base of a film and one (one) in an emulsion layer of a film. All the others occurred in the lead plates. The statistics of the interaction measurements are summarized in Table 2 and evaluated interaction lengths are also presented. The interaction length, λ , is calculated as follows:

$$\lambda & = - \frac { L } { \ln ( 1 - \frac { N } { N _ { 0 } } ) } \\. \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde{ \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilicon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, \ - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \, \tilde { \epsilon } \, - \,$$

where L is the thickness of a lead plate and an emulsion film ( L = 1293 µ m), N is the number of found interactions and N 0 is the sum of the numbers of followed tracks in all analyzed films. Figure 6 shows momentum dependence of the interaction length, which is in good agreement with the Monte Carlo simulation within the statistical errors.

## 5.2. Topological characteristics

The multiplicity and kink angle distributions of relativistic charged secondary particles are shown in Figure 7. Topological characteristics are summarized in Table 3 and they are compared with those of simulated data. As the beam momentum increases, the average number of secondary particles &lt; n &gt; becomes larger. In the case of inelastic interactions, the higher is the momentum, the smaller is the emission angle. This tendency can be seen for 4 GeV/ c and 10 GeV/ c interactions, while it is not so for 2 GeV/ c interactions and

Table 2 Results of the interaction measurements and of the simulated data (MC) analysis. Number of tracks followed, total track length followed in the ECC brick, number of interactions in the ECC brick, evaluated interaction length for each momentum beam are presented.

| P [GeV/ c ]   | 2                        | 2 (MC)                 | 4                        | 4 (MC)                 | 10                       | 10 (MC)                |
|---------------|--------------------------|------------------------|--------------------------|------------------------|--------------------------|------------------------|
| Tracks        | 584                      | 11301                  | 913                      | 9260                   | 2205                     | 13746                  |
| Total L [mm]  | 8506                     | 191794                 | 12620                    | 162009                 | 38534                    | 240922                 |
| Interactions  | 77                       | 1544                   | 68                       | 773                    | 173                      | 1040                   |
| λ [mm]        | 109 . 8 +14 . 1 - 11 . 4 | 123 . 6 +3 . 3 - 3 . 1 | 184 . 9 +24 . 2 - 20 . 1 | 208 . 9 +7 . 8 - 7 . 3 | 222 . 5 +18 . 4 - 15 . 8 | 231 . 0 +7 . 4 - 7 . 0 |

others. This might be explained by the fact that fraction of elastic interactions is larger in 2 GeV/ c interactions [14]. For the average multiplicity &lt; n &gt; , the agreement between the experimental data and the simulated data is not so good in 2 GeV/ c interactions. The other experimental data agree well with the simulation data.

## 5.3. Nuclear fragment association

We present measurement results of nuclear fragments for 1-prong and 3-prong events, which are relevant to topologies of the τ hadronic decay, and corresponding results for simulated data in Table 4. Here numbers of events associated with at least one heavily ionizing particle, pointing to the vertex, found in either the upstream film or the downstream film of the vertex are listed. The probability is the fraction of fragment associated events in the 1-prong or 3-prong events. Figure 8 shows the association probability as a function of the beam momentum, where data of 1-prong and 3-prong events are merged, together with simulation results. At least one nuclear fragment particle is associated with about a half of events for

Fig. 6 Interaction length as a function of beam momentum. Black dots with error bars (outlined circles) show experimental data (simulated data).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0386v1-figure-7dc8a78d-25b5-44b4-a16a-8740226daa81.jpg)

Table 3 Topological characteristics of experimental data and simulated data (MC). Number of events, average charged particle multiplicity &lt; n &gt; , number of 1-prong events, average kink angle for 1-prong events &lt; θ kink &gt; , number of 3-prong events, average of kink angle average for 3-prong events &lt; θ kink &gt; for 2, 4, 10 GeV/ c π -interactions are summarized.

| P [GeV/ c ]      | 2    |   2 (MC) |     4 |   4 (MC) |     10 |   10 (MC) |
|------------------|------|----------|-------|----------|--------|-----------|
| Events           | 77   |  1544    | 68    |   773    | 173    |   1040    |
| < n >            | 0.48 |     0.65 |  0.93 |     1.03 |   2.45 |      2.41 |
| 1-prong events   | 33   |   915    | 29    |   372    |  26    |    216    |
| < θ kink > [rad] | 0.13 |     0.13 |  0.23 |     0.21 |   0.2  |      0.15 |
| 3-prong events   | 0    |     3    |  2    |    32    |  44    |    246    |
| < θ kink > [rad] | -    |     0.36 |  0.32 |     0.3  |   0.27 |      0.24 |

Fig. 7 Multiplicity distributions of relativistic particle tracks in the forward hemisphere (Top) and kink angle distributions of 1-prong events (Bottom) for experimental data (dots with error bars) and simulated data (histogram). The simulated distributions are normalized to the real data.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0386v1-figure-4d82bb6a-3865-43f0-b99c-daae48aa2aac.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0386v1-figure-5ccf8cde-9e47-4487-8473-7db60b6a9339.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0386v1-figure-46e5bde5-304e-4236-821f-387825bd0cef.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0386v1-figure-b4a5eed3-9be7-436d-b2fb-0f72283e6741.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0386v1-figure-57f7fd58-fd29-467c-b464-009209428864.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0386v1-figure-9a44b0a9-5c2f-475a-9749-3c767f69143f.jpg)

beam momentum greater than 4 GeV/ . The experimental data of the fragment associac tion probability are well reproduced by the simulation with differences less than 10%. The number of background nuclear fragment tracks was estimated to be 0.035 tracks/event by scanning around dummy vertices with the same conditions. Figure 9 shows good agreement in multiplicity and polar angle distributions of nuclear fragments between experimental data (dots with error bars) and simulated data (histogram).

Table 4 Results of nuclear fragment search (experimental data and simulated data). Lower limit of the probability shows the value at 90% confidence level.

| P [GeV/ c ]              | 2                     | 2   | 4                     | 4                     | 10                     | 10                    |
|--------------------------|-----------------------|-----|-----------------------|-----------------------|------------------------|-----------------------|
| Prong                    | 1                     | 3   | 1                     | 3                     | 1                      | 3                     |
| Events                   | 32                    | 0   | 29                    | 2                     | 25                     | 41                    |
| Fragment associated      | 10                    | 0   | 16                    | 2                     | 15                     | 27                    |
| Probability [%]          | 31 . 3 +9 . 1 - 6 . 9 | -   | 55 . 2 +8 . 6 - 9 . 3 | > 46 . 5              | 60 . 0 +8 . 9 - 10 . 2 | 65 . 9 +6 . 5 - 8 . 0 |
| Events (MC)              | 908                   | 0   | 372                   | 32                    | 214                    | 246                   |
| Fragment associated (MC) | 213                   | 0   | 219                   | 17                    | 127                    | 169                   |
| Probability (MC) [%]     | 23 . 5 +1 . 5 - 1 . 3 | -   | 58 . 9 +2 . 5 - 2 . 6 | 53 . 1 +8 . 4 - 8 . 8 | 59 . 3 +3 . 3 - 3 . 5  | 68 . 7 +2 . 8 - 3 . 1 |

Fig. 8 Association probability of nuclear fragments as a function of beam momentum. Black dots with error bars (outlined circles) show experimental data (simulated data).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0386v1-figure-ab4efd42-cb16-44b0-ab80-f93e767de992.jpg)

## 5.4. Kinematical characteristics

In OPERA, for τ candidates decaying into a single charged hadron, kinematical selection cuts; p &gt; 2 GeV/ c and p T &gt; . 0 6 GeV/ c are required [19, 20] 4 . Scatter plots of p T versus p of the secondary particles of 1-prong hadron interactions are shown together with those of simulated data in Figure 10. For comparison, scatter plots of p T versus p of the secondary particles of 10 GeV/ c 3-prong interactions are also shown in Figure 11. Statistics are limited because only secondary tracks for which more than 14 lead plates are available as scattering materials are able to be momentum-reconstructed. In this condition, error of the momentum measurement was estimated to be σ p = 45(55)% at p ≤ 4 GeV /c ( p &gt; 4 GeV /c ) by using the multiple scattering measurements of beam pions with known momentum, i.e. 2, 4, 10 GeV/ , c

4 The cut on p T is 0.3 GeV/ c if there is a converted γ attached to the secondary vertex. In this analysis, no attempt is made to find γ s because the numbers of lead plates and emulsion films are limited.

Fig. 9 Top: multiplicity distributions of nuclear fragments for experimental data (dots with error bars) and simulated data (histogram). Bottom: polar angle distributions of nuclear fragments for experimental data (dots with error bars) and simulated data (histogram).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0386v1-figure-cae04d51-e1ec-44ec-921a-721f5d086bae.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0386v1-figure-1df2c48b-7e7e-4b90-a063-ad951a62ceac.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0386v1-figure-104e5cb1-4c98-49b9-b7dc-62121c748c03.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0386v1-figure-f1661189-41ba-41f4-9804-f54ab16d21e9.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0386v1-figure-1c6656d8-cc8a-445a-81c1-98a40d29460a.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0386v1-figure-a4363a17-aae1-41e5-b726-360575d4ee64.jpg)

Table 5 Fraction of secondary particles in the selection domain of p &gt; 2 GeV/ c and p T &gt; . 0 6 GeV/ . c Numbers of secondary tracks for which the momentum is measurable and numbers of tracks in the selection domain are summarized. Upper limits are at 90% confidence level.

| P [GeV/ c ]                     | 2                    | 2        | 4                    | 4       | 10                     | 10                   |
|---------------------------------|----------------------|----------|----------------------|---------|------------------------|----------------------|
| Prong                           | 1                    | 3        | 1                    | 3       | 1                      | 3                    |
| Tracks p measurable             | 12                   | 0        | 7                    | 0       | 10                     | 54                   |
| Tracks in selection domain      | 0                    | 0        | 0                    | 0       | 1                      | 1                    |
| Fraction [%]                    | < 16 . 3             | -        | < 25 . 0             | -       | 10 . 0 +16 . 0 - 3 . 0 | 1 . 9 +3 . 7 - 0 . 4 |
| Tracks p measurable (MC)        | 393                  | 1        | 127                  | 34      | 93                     | 299                  |
| Tracks in selection domain (MC) | 2                    | 0        | 9                    | 0       | 13                     | 19                   |
| Fraction (MC) [%]               | 0 . 5 +0 . 7 - 0 . 2 | < 69 . 0 | 7 . 1 +3 . 0 - 1 . 7 | < 6 . 4 | 14 . 0 +4 . 4 - 2 . 9  | 4 . 8 +1 . 3 - 0 . 9 |

under the same condition. In Figure 10, the dark area defines the domain in which τ decay candidates are selected and the hatched area defines the domain rejected by the selection cuts ( θ kink &gt; . 0 02 rad) and S-UTS angle acceptance (tan θ &lt; 0 6). Table 5 summarizes the . fraction of secondary particles in the selection domain for the events in which secondary particle momenta are measurable. Comparisons between experimental and simulated data show generally good agreement, however the final statistics of the tracks in the selection domain are not sufficient to obtain quantitative results on how good the agreement is. We will therefore consider an alternative method to estimate the goodness of the agreement in the next subsection.

Fig. 10 Top: Scatter plot of p T vs. p of the secondary particle of 1-prong interactions (Experimental data). Bottom: Scatter plot of p T vs. p of the secondary particle of 1-prong interactions (Simulated data). On the figures the dark area defines the domain in which τ decay candidates are selected and the hatched area defines the domain rejected by the selection cuts.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0386v1-figure-ff837019-02d1-4089-82c5-f215b830a735.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0386v1-figure-a3d1bc07-2f3f-4ba3-8856-18c99fa79c0d.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0386v1-figure-f74b5ff8-61d7-4b6f-a4ae-f37ac81a57bf.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0386v1-figure-624b81c6-c055-41d5-b4f3-dd678b86bfc3.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0386v1-figure-7ffa8f13-d974-4ef4-9659-b16df0d9afc2.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0386v1-figure-5dd729ad-1f37-473f-adac-0d29741295fb.jpg)

## 5.5. Goodness of agreement between the experimental and the simulated data

The numbers of secondary particles in the domain where τ decay candidates are selected are too few to evaluate the goodness of the agreement between the experimental and the simulated data. So we try to widen the selection condition for comparison. We select 1-prong hadron interaction events with θ kink &gt; 0.30, 0.15 and 0.06 rad for 2 GeV/ , 4 GeV/ c c and 10 GeV/ , respectively. These angles correspond to c p T &gt; 0 6 GeV/ . . c

There are some events in the domains defined by these conditions, as shown in Table 6. We use these samples to check the agreement between the experimental and the simulated data.

Fig. 11 Left: Scatter plot of p T vs. p of the secondary particle of 10 GeV/ c 3-prong interactions (Experimental data). Right: Scatter plot of p T vs. p of the secondary particle of 10 GeV/ c 3-prong interactions (Simulated data). On the figures the dark area defines the domain in which τ decay candidates are selected and the hatched area defines the domain rejected by the selection cuts.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0386v1-figure-db523777-d5cf-4e68-8778-741264ac9cf7.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0386v1-figure-3e63e18a-543e-44a0-9734-8eb466bd3ca9.jpg)

Table 6 Number of selected events for each momentum beam. We select 1-prong events with θ kink &gt; 0.30, 0.15 and 0.06 rad for 2 GeV/ , 4 GeV/ c c and 10 GeV/ , respectively. For c the simulated data (MC), numbers of events normalized to the total track length followed in this experiment are shown. Relative differences δ sys between experimental data and simulated data are less than 30% and are within statistical errors δ stat of experimental data although statistics is not sufficient for 2 GeV/ c .

| P [GeV/ c ]      |   2 |   2 (MC) |   δ sys [%] |   4 |   4 (MC) |   δ sys [%] |   10 |   10 (MC) |   δ sys [%] |
|------------------|-----|----------|-------------|-----|----------|-------------|------|-----------|-------------|
| Events           |  77 |     68.5 |          12 |  68 |     60.2 |          13 |  173 |     166   |           4 |
| 1-prong          |  33 |     40.6 |         -19 |  29 |     29.1 |           0 |   26 |      34.6 |         -25 |
| After θ kink cut |   6 |      7.5 |         -20 |  19 |     15.3 |          24 |   17 |      23.7 |         -28 |

In Table 6, numbers of simulated data (MC) were normalized to the total track length followed for each beam momentum in this experiment. Then the statistical error δ stat and the relative difference δ sys are defined as follows,

$$\delta _ { \text{stat} } \ = \ \frac { 1 } { \sqrt { N _ { e v } } } \ \ _ { \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

$$\delta _ { s y s } \ = \ \frac { N _ { e v } - N _ { M C } } { N _ { M C } }$$

where N ev and N MC are numbers of the experimental data and the simulated data respectively. It is found that all the relative differences δ sys are less than 30% and can be understood within the statistical errors δ stat although statistics is not sufficient for 2 GeV/ c . We therefore conclude that the agreement between the experimental and the simulated data has been confirmed at the 30% level for beam momentum P ≥ 4 GeV/ . c

## 6. Conclusions

Topological and kinematical characteristics of hadron interactions have been studied by using an ECC brick exposed to 2, 4, 10 GeV/ c pion beams. High speed automated microscope system was employed to analyze the nuclear emulsion films. A total of 318 hadron interactions were found and reconstructed by following 60 m π -tracks in the brick. Secondary charged particle tracks from interaction vertices were also followed and reconstructed. Charged particle multiplicity of each event and emission angle of each secondary particle were measured and their distributions were found to be in good agreement with a FLUKA Monte Carlo simulation. Nuclear fragments were also searched for by newly developed automated microscope system with a wide view. We measured the probability for the interaction vertices to be associated with nuclear fragments and found it to be greater than 50% for beam momentum P &gt; 4 GeV/ . The experimental data of the fragment association probability c are well reproduced by the simulation with differences less than 10%. When possible, the momentum of the secondary particle was measured by using the coordinate method. Fractions of 1-prong and 3-prong hadron interactions being in the domain where τ → hadron decays are selected are measured and found to be consistent with the simulation at the 30% level for beam momentum P ≥ 4 GeV/ . We conclude that the FLUKA based simulation c reproduces well the experimental distributions of 2, 4, and 10 GeV/ c π -interactions in the

ECC brick. This result can be applied in neutrino oscillation experiments to evaluate the hadronic background contaminations on τ lepton decays. Since the background contaminations originate from all charged hadron interactions, mainly π -, π + and proton interactions, positively charged hadron interactions should also be investigated. This will be the subject of further studies.

## Acknowledgments

We wish to thank P. Vilain for his careful reading of the manuscript and for his valuable comments. We wish to express our gratitude to the colleagues of the Fundamental Particle Physics Laboratory, Nagoya University for their cooperation. For the beam exposure, we gratefully acknowledge the support of the PS staff at CERN. We warmly acknowledge the financial support from the Japan Society for the Promotion of Science (JSPS), the Promotion and Mutual Aid Corporation for Private Schools of Japan, and Japan Student Services Organization.

## References

- [1] K. Kodama et al., (DONUT Collaboration), Phys. Lett. B, 504 , 218 (2001).
- [2] K. Kodama et al., (DONUT Collaboration), Nucl. Instrum. Methods Phys. Res., Sect. A 493 , 45 (2002).
- [3] M. Guler et al., (OPERA Collaboration), An appearance experiment to search for ν µ → ν τ oscillations in the CNGS beam: experimental proposal, CERN-SPSC-2000-028, LNGS P25/2000 (2000).
- [4] M. Guler et al., (OPERA Collaboration), Status Report on the OPERA Experiment, CERN/SPSC 2001-025, LNGS-EXP 30/2001 add. 1/01 (2001).
- [5] R. Acquafredda et al., (OPERA Collaboration), J. Instrum., 4 , P04018 (2009).
- [6] O. Yasuda, arXiv:1211.7175 [hep-ph].
- [7] T. Ota, J. Phys. Conf. Ser., 408 , 012021 (2013).
- [8] K. Morishima, T. Nakano, J. Instrum., 5 , P04011 (2010).
- [9] T. Fukuda et al., J. Instrum., 8 , P01023 (2013).
- [10] T. Nakamura et al., Nucl. Instrum. Methods Phys. Res., Sect. A 556 , 80 (2006).
- [11] D.J. Simon et al., SECONDARY BEAMS FOR TESTS IN THE PS EAST EXPERIMENTAL AREA, PS/PA-EP Note 88-26 (1988), PS/PA/Note 93-21 (1993), http://ps-div.web.cern.ch/ps-div/Reports/PA9321/.
- [12] K. Kodama et al., Nucl. Instrum. Methods Phys. Res., Sect. A 574 , 192 (2007).
- [13] N. Agafonova et al., (OPERA Collaboration), New J. Phys., 14 , 013026 (2012).
- [14] J. Beringer et al. (Particle Data Group), Phys. Rev. D, 86 , 010001 (2012).
- [15] M. Kimura et al., Nucl. Instrum. Methods Phys. Res., Sect. A 711 , 1 (2013).
- [16] A. Ferrari, P.R. Sala, A. Fass`, and J. Ranft, FLUKA: a multi-particle transport code, CERN-2005-10 o (2005), INFN/TC 05/11, SLAC-R-773.
- [17] G. Battistoni, S. Muraro, P.R. Sala, F. Cerutti, A. Ferrari, S. Roesler, A. Fass`, J. Ranft, The FLUKA o code: Description and benchmarking, Proceedings of the Hadronic Shower Simulation Workshop 2006, Fermilab 6-8 September 2006, M. Albrow, R. Raja eds., AIP Conference Proceeding 896, 31-49, (2007).
- [18] F. Ballarini et al., Nuclear models in FLUKA: present capabilities, open problems and future improvements, SLAC-PUB-10813 (2004).
- [19] N. Agafonova et al., (OPERA Collaboration), Phys. Lett. B, 691 , 138 (2010).
- [20] N. Agafonova et al., (OPERA Collaboration), New J. Phys., 14 , 033017 (2012).
Hirokazu Ishida Tsutomu Fukuda Takafumi Kajiwara Koichi Kodama Masahiro Komatsu Tomokazu Matsuo Shoji Mikado Mitsuhiro Nakamura Satoru Ogawa Andrey Sheshukov Hiroshi Shibuya Jun Sudou Taira Suzuki Yusuke Tsuchida
Study of hadron interactions in a lead-emulsion target
Topological and kinematical characteristics of hadron interactions have been
studied using a lead-emulsion target exposed to 2, 4 and 10 GeV/c hadron beams.
A total length of 60 m $\pi^-$ tracks was followed using a high speed automated
emulsion scanning system. A total of 318 hadron interaction vertices and their
secondary charged particle tracks were reconstructed. Measurement results of
interaction lengths, charged particle multiplicity, emission angles and momenta
of secondary charged particles are compared with a Monte Carlo simulation and
appear to be consistent. Nuclear fragments emitted from interaction vertices
were also detected by a newly developed emulsion scanning system with
wide-angle acceptance. Their emission angle distributions are in good agreement
with the simulated distributions. Probabilities of an event being associated
with at least one fragment track are found to be greater than 50% for beam
momentum $P > 4$ GeV/c and are well reproduced by the simulation. These
experimental results validate estimation of the background due to hadron
interactions in the sample of $\tau$ decay candidates in the OPERA $\nu_{\mu}
\to \nu_{\tau}$ oscillation experiment.