First neutrino event detection with nuclear emulsion at J-PARC neutrino beamline

Transcription

1 Prog. Theor. Exp. Phys. 2015, (26 pages) DOI: /ptep/ First neutrino event detection with nuclear emulsion at J-PARC neutrino beamline arxiv: v2 [hep-ex] 17 May 2017 T. Fukuda 1,, S. Aoki 2, S. Cao 3, N. Chikuma 4, Y. Fukuzawa 1, M. Gonin 5, T. Hayashino 6, Y. Hayato 7, A. Hiramoto 6, F. Hosomi 4, K. Ishiguro 1, S. Iori 8, T. Inoh 8, H. Kawahara 1, H. Kim 9, N. Kitagawa 1, T. Koga 4, R. Komatani 1, M. Komatsu 1, A. Matsushita 1, S. Mikado 10, A. Minamino 11, H. Mizusawa 8, K. Morishima 1, T. Matsuo 8, T. Matsumoto 8, Y. Morimoto 8, M. Morishita 1, K. Nakamura 6, M. Nakamura 1, Y. Nakamura 1, N. Naganawa 1, T. Nakano 1, T. Nakaya 6, Y. Nakatsuka 1, A. Nishio 1, S. Ogawa 8, H. Oshima 8, B. Quilain 6, H. Rokujo 1, O. Sato 1, Y. Seiya 9, H. Shibuya 8, T. Shiraishi 1, Y. Suzuki 1, S. Tada 1, S. Takahashi 2, K. Yamada 2, M. Yoshimoto 1 and M. Yokoyama 4 1 Nagoya University, Nagoya , Japan 2 Kobe University, Kobe , Japan 3 High Energy Accelerator Research Organization (KEK), Tsukuba , Japan 4 University of Tokyo, Tokyo , Japan 5 Ecole Polytechnique, IN2P3-CNRS, Laboratoire Leprince-Ringuet, Palaiseau, France 6 Kyoto University, Kyoto , Japan 7 University of Tokyo, ICRR, Kamioka , Japan 8 Toho University, Funabashi , Japan 9 Osaka City University, Osaka , Japan 10 Nihon University, Narashino , Japan 11 Yokohama National University, Yokohama , Japan tfukuda@flab.phys.nagoya-u.ac.jp Precise neutrino nucleus interaction measurements in the sub-multi GeV region are important to reduce the systematic uncertainty in future neutrino oscillation experiments. Furthermore, the excess of ν e interactions, as a possible interpretation of the existence of a sterile neutrino has been observed in such an energy region. The nuclear emulsion technique can measure all the final state particles with low energy threshold for a variety of targets (Fe, C, H 2 O, and so on). Its sub-µm position resolution allows measurements of the ν e cross-section with good electron/gamma separation capability. We started a new experiment at J-PARC to study sub-multi GeV neutrino interactions by introducing the nuclear emulsion technique. The J-PARC T60 experiment has been implemented as a first step of such a project. Systematic neutrino event analysis with full scanning data in the nuclear emulsion detector was performed for the first time. The first neutrino event detection and its analysis is described in this paper Subject Index xxxx, xxx 1. Introduction Currently, many experimental neutrino oscillation projects are running or planned to search for the CP violation, and/or to probe the neutrino mass hierarchy [1 9]. Precise cross-section measurements of neutrino-nucleus interactions in the sub-multi GeV region are important to c The Author(s) Published by Oxford University Press on behalf of the Physical Society of Japan. This is an Open Access article distributed under the terms of the Creative Commons Attribution License ( which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

2 reduce the systematic uncertainty, specially for neutrino mixing angle θ 23, in these experiments. Furthermore, one possible explanation of the MiniBooNE anomaly [10] in this energy region is a result of the existence of so-called sterile neutrino indicated by LSND [11]. Therefore that is also important to verify the MiniBooNE anomaly. Nuclear emulsion is a three-dimensional solid tracking detector with sub-micron positional resolution. Thanks to its high spatial resolution, the nuclear emulsion technology contributed to the discovery of the pion [12], of the charmed particle in cosmic rays [13], the direct observation of ν τ [14], the discovery of ν τ appearance in a ν µ beam [15] and so on. This high spatial resolution allows measurements of all the final state particles with low-energy threshold. Nuclear emulsion has 4π track detection capability [16]. This is also very useful to investigate the event topology in low energy neutrino interactions. Furthermore, good e/γ separation capability allows measurement of ν e CC interactions with a strong suppression of the background from ν µ NC interactions with π 0 production. There is much flexibility for target material selection because the detector, a so-called Emulsion Cloud Chamber (ECC), is constructed as a sandwich structure of thin nuclear emulsion films and the target materials. The J-PARC T60 experiment [17] was proposed to study the feasibility and check the detector performance of nuclear emulsion at the J-PARC neutrino beamline. The first results of beam exposure in 2015 will be presented in this paper. Accelerator-based neutrino experiments with nuclear emulsion were established from the end of 1970s [18 20]. The experimental concept is that of a hybrid analysis with the emulsion and the electronic detectors, so-called Emulsion-Counter Hybrid Experiment. Only a small area of emulsions that is predicted by the electronic detectors has been analyzed, or only a small number of tracks that are predicted have been traced back because of the limitation of the track readout speed in the emulsions. The high-speed scanning system of the emulsion has been developed [21, 22] and the track readout of the entire area in the emulsion has recently become available. This allows one to analyze the emulsions independently of the prediction from the electronic detectors. In this paper, the first systematic neutrino event analysis with scanning all data in the nuclear emulsion will be also described. 2. Neutrino Detector 2.1. Nuclear Emulsion Films Nuclear emulsion gel was produced at the facility of Nagoya University. In this emulsion gel production, high-sensitivity emulsion gel, containing 55% AgBr crystal by volume ratio, was produced [23]. Two emulsion layers, each 50 µm thick, were formed on both faces of a 180-µm thick polystyrene base. The size of sheets was 100 mm 125 mm. The long-term variation of the emulsion sensitivity and background noise was investigated using cut pieces of emulsion sheets. The emulsion sensitivity was measured as grain density (GD), the number of grains per 100 µm along a minimum ionizing particle (MIP) trajectory. The noise of the emulsion was measured as fog density (FD), the number of randomly distributed grains per 1000 µm 3 in the emulsion. If an emulsion is left undeveloped for a long time after track recording, the silver atoms tend to oxidize into silver ions, owing to the presence of water and oxygen in the material. Some of the latent image cores are destroyed and fewer grains are produced in the development process. This effect is called fading. Variations of GD and FD over 133 days are shown in Fig.1 under four different temperature and humidity conditions. We found that GD and FD remain at a safe levels, even after 133 days at 23 C. 2/26

3 Fig. 1 Initial and long-term performance of emulsion gels. The GD and FD is kept at a safe levels ( GD 35, FD 10 ) for 133 days at less than 23 C. The fading effect is suppressed in low-temperature and low-humidity conditions Emulsion Cloud Chamber (ECC) The ECC is constituted by alternatively stacking 41 emulsion films and 40 stainless steel plates (SUS304, 500 µm thickness) as shown in Fig.2. These are inserted in a paper box and then vacuum-packed with a light-shield aluminum laminated bag at 40% R.H. condition. The total target mass is approximately 2.2 kg (1.2 kg inside the fiducial volume). Fig. 2 The Emulsion Cloud Chamber (ECC) 3/26

4 2.3. Emulsion Shifter Precise tracking detectors with timing information, e.g. the silicon microstrip detector or the scintillating fiber detector, have been applied to connect tracks between the main emulsion detector and the electronic detectors, which have muon-identification capability in the conventional Emulsion-Counter Hybrid Experiment. In the present work, an emulsion device with timing information played this role. The emulsion detector intrinsically has no timing information. Therefore, the timing information is produced from the emulsion tracks by track coincidence between the films on the a mechanically controlled moving stage [24, 25], a so-called the emulsion shifter. An emulsion shifter (Fig.3) was reused from a balloon experiment with nuclear emulsion [26]. This equipment gives a time-stamp to the emulsion track by using a clock-based multi-stage emulsion shifter technique. Three stages play the roles of the hour, minute and second hands of a watch, respectively. The emulsion tracks that are coincident on the all three stages obtained the time-stamp corresponding to the each stage position. The timing resolution of this emulsion shifter is approximately 7.9 seconds, which was investigated using tracks passing through the emulsion shifter and electronic detector located at the back of the emulsion detectors. The details for this equipment and the analysis are described in [26, 27]. Fig. 3 Emulsion shifter 2.4. Interactive Neutrino GRID (INGRID) INGRID is one of the near detectors for the T2K experiment [28]. The tracking device is a 1-cm-thickness plastic scintillator (5 cm width) read out with wave length-shifting fibers. The constitution of the detector is a sandwich structure of the tracking planes and 6.5-cmthickness iron blocks. The detector was assembled from a total of nine iron layers (58.5 cm) and 11 vertical and horizontal tracking planes. The size of an module is 1.2 m 1.2 m 0.9 m. INGRID consists of 16 identical modules arranged in horizontal and vertical arrays around the beam center. The original purpose of INGRID is to measure the neutrino beam direction and verify the number of neutrino interactions per protons on target. In the J-PARC T60 experiment, INGRID is utilized to identify the muons from neutrino nucleus interactions in the ECC. As shown in Fig.4, the emulsion shifter is placed between the ECC and INGRID. The particles that cross the emulsion shifter are supplied time-stamps. The tracks are then matched between the emulsion detector and INGRID using timing information. 4/26

5 Fig. 4 Detector concept 2.5. Installation at J-PARC The emulsion detectors and the assembling frame were prepared at Toho University. These detectors were transported to J-PARC before the installation. The emulsion detectors were installed in front of an INGRID module adjacent to the on-axis module of the Neutrino Monitor Building in J-PARC on 14th Jan. 2015, as shown in Fig.5. The operation of the emulsion shifter has been monitored by verifying the stage reproducibility, which was sent to Kobe University with a monitoring laptop PC. Fig. 5 Detector installation at J-PARC neutrino beamline 5/26

6 3. Neutrino Beam Exposure The neutrino beam exposure started in 15th Jan Unfortunately, beam exposure was stopped on 16th Jan., then re-started from 25th Feb. and finished on 1st Apr. The beam exposure was performed with the anti-neutrino beam mode after re-start of beam exposure. The accumulated number of protons on target is approximately The ratio of neutrino and anti-neutrino interactions is 1:3 because of the contamination of neutrinos in anti-neutrino beam. The total number of neutrino or anti-neutrino events is estimated as approximately 44 events (Fiducial Volume: 24 events, including two events at the region of the emulsion layer). After beam exposure, the stages of the emulsion shifter were moved to a reference position and reference cosmic-ray tracks were accumulated for one week. The emulsion detectors were uninstalled on 8th Apr. and all films were immediately developed at Toho University. The temperature at the experimental site was approximately 22 C during the exposure, as shown in Fig.6. Therefore, the tracks that are recorded in emulsion films were kept in good condition, as shown in Fig.1. The operation of the emulsion shifter was well performed [27]. Fig. 6 Temperature at the experimental site 6/26

7 4. Event Analysis 4.1. Scanning The entire areas of all the emulsion films for the experiment ( 1.2 m 2 ; for two emulsion layers of all the emulsion films, including six films for the emulsion shifter) were scanned using the high-speed automated microscope system, Hyper Track Selector (HTS) [22] (Fig.7). This system has been developed in Nagoya University to quickly analyze particle tracks in emulsion films and allows to readout tracks at a scanning speed of 5,000 cm 2 /h. HTS takes 16 tomographic images in an emulsion layer and recognizes a series of grains on a straight line as a track, by taking a coincidence among these images. In this time, the slope acceptance of HTS is set at tanθ < 1.7 ( θ 60 degree), where θ is the track angle with respect to the perpendicular of the emulsion film (the z-axis). Fig. 7 Hyper Track Selector (HTS) 4.2. Analysis in ECC Track Reconstruction. First, track segments, so-called micro tracks, were detected on each layer of the emulsion films (Fig.8). The positions (x, y), slopes (tanθ x, tanθ y ), number of hits in an emulsion layer, Pulse Height (PH), and number of pixels constructed track data, Volume Pulse Height (VPH) of the micro tracks, were measured by HTS. The number of grains belonging to the track (blackness) is well correlated with the PH and the VPH. As shown in Fig.8, the tracks connecting the positions of micro tracks on both side of the plastic base are called base tracks. Sub-MeV or few MeV electrons from the environment are recorded in the emulsion because nuclear emulsion records all tracks from immediately after production until development with no dead time. Therefore, the automated scanning system recognizes fake tracks made of low-energy electrons and/or random noise grains (fog). The chance coincidence of these fake tracks become the background of 7/26

8 base tracks. However, such background base tracks can be rejected by evaluating their blackness (VPH) and linearity (the angle difference between base track and micro tracks) [29]. Fig.9 shows the blackness vs. the linearity of base tracks for a small-angle region ( tanθ < 0.1 : left figure) and large-angle region ( tanθ > 1.4 : right figure). The VPH and angular accuracy of base tracks depend on their angle. The smaller the track angle, the better the signal-to-noise ratio (S/N) of the base tracks. So the discrimination line is set for each track angle as shown by the red line in Fig.9. The tracks under the red line are rejected as the fake track. Fig. 8 Schematic view of the track reconstruction in an emulsion film. Fig. 9 Track linearity and blackness (left: tanθ < 0.1, right: tanθ > 1.4). The sum of VPH of micro tracks, which is a reconstructed base track is shown on the vertical axis. The square-root of the sum of squares of the angle difference between the base track and micro tracks is shown on the horizontal axis. The blackness and linearity of fake tracks tend to be thin and bad, respectively. The red line shows the discrimination line and is set for each 100 mrad of track angle. The positions (x, y) and slopes (tanθ x, tanθ y ) of the base tracks were measured with respect to the microscope s coordinate system. After connecting the base tracks among the 8/26

9 emulsion films, the rotation, slant, parallel translation, and distance between every pair of adjacent emulsion films are adjusted, so that the differences in the track positions and slopes between the two films are minimized after the reconstruction procedure. These connection elements of base tracks are called linklets. All base tracks in all films are reconstructed by tracking linklets, as the chain. Finally, the track density in a film is 4,000 tracks / cm 2. The track efficiency after reconstruction was evaluated at each measured angle for all plates, as shown in Fig.10. In the small-angle region, the track efficiency is approximately 96 %. However, in the large-angle region, the efficiency loss caused is found. Fig. 10 Track efficiency after reconstruction. When the track efficiency is estimated at one plate, having base tracks at upstream and downstream plates is required to select penetrating tracks in the target plate. 9/26

10 Proton Identification. The momentum of charged particles can be estimated in the ECC by measuring their Multiple Coulomb Scattering (MCS) [30]. When a particle of charge z, momentum p, velocity βc traverses a material of depth x and radiation length X 0, the deviation of the distribution of scattering angle is given by Eq.1. θ 0 = 13.6MeV βcp x z ( ln( x )) (1) X 0 X 0 The GD is nearly proportional to its energy loss de/dx in nuclear emulsion [31]. The average number of hit layers (PH) or the average number of pixels associated to the reconstructed track (VPH) in scanned track data corresponds to the GD of tracks [32]. In Fig.11, the measured momentum and the average VPH of each track that is connected to more than 10 segments is shown. The proton-like tracks are found in the low-momentum and high-de/dx region. Fig. 11 Proton-like tracks are seen in the low pβ and high average VPH region Track Selection. The reconstructed tracks are selected to search for neutrino event vertices as shown the following steps. (1) Fiducial area selection : the acceptance of the average track angle is less than tanθ = 1.6 and the searched position area is set at x: mm and y: 6 76 mm to avoid the bad film condition area (thin thickness of emulsion layer). The plate ID of the most upstream and downstream plates in the neutrino beam direction is given as PL41 and PL01, respectively. The number of tracks and the number of tracks that just start and connect to downstream plates after fiducial selections at each plate is described in Fig.12-(a). As shown by the red line in Fig.12-(a), there are many tracks at the upstream of the detector that are identified as penetrating tracks (cosmic ray muons or muons produced from outside of the detector). 10/26

11 (2) Penetrate track veto : tracks that have at least one segment on three plates at the most upstream position of the ECC (PL41, PL40, PL39) are rejected as penetrate tracks. (3) Edge out track veto : tracks are rejected as edge out tracks in the case that the position extrapolated from the most upstream plate to its upstream plate using their position and angle is out of the fiducial area. The number of tracks and the number of starting tracks after veto selections at each plate is described in Fig.12-(b). The angle distribution of penetrate and edge out tracks is shown in Fig.13-(a). This is consistent with cosmic rays from the sky. Figure13-(b) shows the track angle distribution after (1) (3) selections. The low-energy track or the electron/positron from gamma-ray or low efficiency tail from large angle tracks still remain as background for starting tracks in the ECC. (4) Final track selection for vertex search : the tracks that have at least three segments and for which the track angle is less than tanθ = 1.5 ( θ 56 degree) at the most upstream plate are selected to avoid accidental disconnecting tracks in the large-angle region and low reliability tracks. Finally, tracks are selected as the candidate tracks from neutrino interactions. The number of tracks and the number of starting tracks after these selections at each plate is described in Fig.12-(c). Additionally, the track angle distribution after this selection is shown in Fig.13-(c). Fig. 12 (a): Number of tracks (black) and starting tracks (red) at each plate after fiducial selections. (b): Number of tracks (black) and starting tracks (red) at each plate after veto selections. (c): Number of tracks (black), number of starting tracks (blue), and number of starting tracks that have at least three segments (red) at each plate after final track selection. 11/26

12 Fig. 13 (a): Angle distribution of penetrate and edge out tracks. (b): Angle distribution of remaining tracks after penetrate and edge out track cut. (c): Angle distribution of finally selected tracks for neutrino event search Event Selection. The event selection is performed with the following procedure. The events are then categorized into the types of (a) (h), as shown in Fig.14. (1) Vertex Search: the most upstream base tracks converge at a point within 3 σ position difference that is calculated from their angular accuracy and the MCS assumed the momentum of 400 MeV/c. For example, the position allowance for the minimum distance between two tracks (tanθ = 0.5, respectively) is approximately 20 µm at the level of 500 µm of iron above the surface of emulsion film. In this search, the tracks that emit to the forward direction are selected as vertex tracks. Most of starting tracks are the accidental disconnecting tracks or the stopped cosmicray tracks from the back of the ECC. So the positions of these tracks are random and fake vertices are created with their accidental chance coincidence. Such background vertices can be emulated by shuffling only track positions randomly in each plate and making vertices. Repeating 10 times procedure, the number of background is estimated as the average value in Fig.14. (2) Track Multiplicity: as a first step, reconstructed vertices are categorized based on their track multiplicity. The three categories are three and more tracks, two tracks and single track. This three and more tracks is categorized as Category-(a). 12/26

13 In Category-(a), seven vertices composed of 25 tracks are found and the estimated number of background vertex is less than 0.1. Thus, the vertices of Category-(a) are real vertex candidates because of low background. In the two-track vertex category, there are 269 vertices composed of 537 tracks including one duplicated case, as shown in Fig.15, corresponding to background vertices. Then, tracks are categorized as single track vertex event (Category-(c)). Single track vertices are mainly considered to be disconnecting of low-energy cosmic ray muons or one of an electron/positron pair production, as shown in Fig.13-(c). (3) Track Blackness: two-track vertices are categorized according to track blackness. Black tracks are defined as tracks whose average VPH is over 120, and minimum ionizing particle (MIP) tracks are defined as tracks whose average VPH is less than 120, as shown in Fig.11. When two-track vertices are reconstructed by only black tracks, they belong to Category-(b). In Category-(b), 18 vertices composed of 36 tracks are found, corresponding to 0.5 background vertices. Thus, they are also real vertex candidates. When one black track and one MIP track made a vertex, there are 20 vertices with 40 tracks and the number of the estimated background is (4) Additional Black Search: the short black tracks made by only one or two plates are added as vertex tracks. This procedure is done by the short black tracks both the forward and backward directions. When the short black track is searched around 20 vertices that are reconstructed by one black track and one MIP track, four new vertices with multiple black tracks are found with less than 0.1 background (Category-(e)). Other vertices are stored in Category-(d). 231 two track vertices that are reconstructed by only MIPs are also attached to additional black tracks. One new vertex with two MIPs and one short black track is found. These vertices are stored in Category-(f). (5) Track Parallelity: Fig.16-(a) shows the opening angle of 230 MIP vertices with no additional black tracks. The opening angle of ten times of background data is shown in Fig.16-(b). According to the comparison between real and background data, an excess of small opening angle is found in real data. They are almost parallel tracks, so are explained as the electron/positron pair production from gamma-ray. Gamma-ray candidate events are selected by using their opening angle and minimum distance of vertex tracks as shown in Fig.16-(c). 125 vertices are selected as gamma-ray events and the estimated number of background is 0.5 as shown in Fig.16-(d). Gamma-ray events are stored in Category-(g) and the other vertices are stored in Category-(h). Through the event selections, Categories (a), (b), (e), (f), and (g) are selected as a vertex events with high reliability. The 125 events in Category-(g) are gamma-ray candidates because they have no black track and their opening angle is small. The 18 vertices in Category-(b) are reconstructed by only black tracks, so it is considered that these events are from neutron interactions or a part of neutrino neutral current interactions. The neutrino candidate events are stored in Categories (a), (e), and (f). The 12 vertices in these categories have some MIP and black tracks. On the other hand, there may be some neutrino events in Category-(d) because the number of vertex signals is 5.2 (= ). The near vertices in Category-(h) are accidental chance coincidence of misconnected cosmic-ray muons or one of electron/positron pair. It is considered that they are same group as Category-(c). 13/26

14 Fig. 14 Flow chart for event selection 14/26

15 Fig. 15 Two vertices are reconstructed by three tracks Fig. 16 Opening angle between two tracks in real data (a) and 10 times of the estimated background data (b). The opening angle vs. minimum distance between two tracks in real data (c) and 10 times of the estimated background data (d). Gamma-ray events are selected within the dotted red line in (c) and (d). 15/26

16 4.3. Gamma-ray events In this section, Category-(g), or gamma-ray candidate events, will be discussed. Only forward going gamma-ray is reconstructed in this category. The angle distribution of gamma-ray candidate events is shown in Fig.17-(a). The incident direction of tracks from gamma-ray is forward and upper side and it is consistent with the direction of cosmic-rays (Fig.17- (b)). Additionally, the number of converted gamma-ray is high in the upstream region of the detector, as shown in Fig.18. This is evidence that the gamma-rays come from outside the detector. The detected electron-positron pairs from gamma-ray in the ECC is shown in Fig.19. It is considered that cosmic gamma-rays from the backward and upper side directions make the single-track vertex in the Category-(c) and background vertices. Fig. 17 Angle distribution of electron/positron (red or blue) from gamma-ray. Fig. 18 Converted plate distribution of gamma-rays. 16/26

17 Fig. 19 The 3 dimensional (x-y, x-z, y-z) and zooming view. Neutrino beam direction is same as z axis. 17/26

18 4.4. Black track vertex events In this section, a discussion of 18 found vertices that are reconstructed by only black tracks in Category-(b) will be presented. The chance coincidence background is estimated as 0.5 events. Therefore, almost events are real vertices. These events are considered as the neutron interactions or a part of neutrino neutral current interactions with no MIP tracks. Additional short or backward black tracks attached to these vertices were also searched for. As shown in table.1, new black tracks were found in more than half of events. The event feature is shown in Fig.20. The emission angle and range distributions are shown in Fig.21-(a) and (b), respectively. Under the current selection criteria, the existence of two black tracks that are emitted to the forward direction is required. Thus, this becomes one of the biases to drop neutron events from the backward direction. Table 1 Black track vertices 2 track vertex 3 track vertex 4 track vertex Fig. 20 Three-dimensional view of black vertex events. Neutrino beam direction is parallel to the z-axis. 18/26

19 Fig. 21 Emission angle distribution (a) and the range distribution (b) of tracks in black vertex event Neutrino Candidate Events The other converged tracks, for example, low-energy tracks, tracks emitted to the backward direction, or electron/positron pair from vertices are searched for in the 12 vertices found in Categories (a), (e), and (f). The detected particle momenta are then measured with MCS and the timing information is obtained with the emulsion shifter to connect the INGRID and identify muons [27]. In this section, the event summary is shown in table 2 and then some of the detailed event features will be described. Table 2 Event summary Event number VTX Plate material Multiplicity (MIP) Multiplicity (Black) Multiplicity (e + e ) 1 8 iron iron iron iron iron iron emulsion emulsion iron iron iron iron /26

20 The position distribution of 12 vertices is shown in the Fig.22. The impact parameter (IP) and emission angle for each MIP and black tracks is also shown in the Fig.23. Fig. 22 Three-dimensional view of neutrino candidate events. Neutrino beam direction is parallel to the z-axis. Fig. 23 Impact Parameter and Emission angle for MIP and black 20/26

21 Event.7. The vertex is in emulsion region at plate.27 and its vertex point is (x, y) = ( µm, µm) in the ECC coordinate system (the origin (0.0 µm, 0.0 µm) is the corner at PL01). One electron/positron pair from a gamma-ray in Category-(g) points to this vertex point. The detailed track information is listed in the table 3. There are two MIP tracks and two black tracks in the vertex (Fig.24-(a)). Also, this event was created at the emulsion region, so we could see more detail with microscope. As shown in Fig.24-(b), the additional very large angle tracks and black tracks, shown as red arrows, were observed. Additionally, this event was given the time stamp by analyzing the emulsion shifter, then the timing matching between the emulsion detector and INGRID succeeded uniquely. As shown in Fig.25 and Fig.26, a muon track is found and the event feature is well matched. Furthermore, the measured track momenta in ECC are consistent with the measured track ranges in INGRID. The momenta of muon candidate track is also enough to pass through INGRID. PL segment id tanθ x tanθ y (µ m) PH AV. VPH AV. Table 3 Track information in Event.7 (GeV/c) PID IP pβ proton proton muon MIP e ± e Fig. 24 Event.7 21/26

22 Fig. 25 Event.7: Hybrid analysis with INGRID Fig. 26 Event.7: Hybrid analysis with INGRID 22/26

23 Event.3. The vertex is steel region at upstream of plate.14 and its vertex point is (x, y) = ( µm, µm) in the ECC coordinate system. There are three MIP tracks and one black track at the vertex. The details of the event tracks are listed in table 4. This event is also matched between the emulsion detector and INGRID with the timing information, as shown in Fig.28 and Fig.29. The measured track momenta are consistent with the measured range in INGRID. Unfortunately, since two tracks penetrate through INGRID, it is difficult to decide which is the muon track from the information of INGRID. PL segment id tanθ x tanθ y (µ m) PH AV. VPH AV. Table 4 Track information in Event.3 (GeV/c) PID IP pβ muon/pion pion/muon proton MIP e ± Fig. 27 Event.3 23/26

24 Fig. 28 Event.3: Hybrid analysis with INGRID Fig. 29 Event.3: Hybrid analysis with INGRID 24/26

25 5. Future Prospects 5.1. Improvement of the Systematic Analysis Further neutrino events are going to be searched for in Categories (c), (d), and (h). It is important to reject many fake starting tracks to find new neutrino events. To do that, the energy threshold to connect tracks between emulsion films will be lowered. Another method is to reconstruct the electron/positron pairs of gamma-rays and neutron-interacted vertices from the backward direction. Then, the search for additional much larger angle black tracks will be also effective to identify the interaction point with large-angle track scanning [16]. Furthermore, when some tracks reached the emulsion shifter, it will be possible to take the coincidence with timing information for background rejection Future Experimental Plan In 2016, anti-neutrino beam was exposed to 60 kg steel target ECC and film scanning was done. Currently the neutrino event analysis for this experiment is in progress to investigate the detector performance with higher statistics. Approximately 3,000 anti-neutrino interaction events are expected to accumulate in this ECC. Moreover, not only anti-muon neutrino events but also anti-electron neutrino events will be detected. Furthermore, the neutrino analysis will be proceeded in a few-kg water target emulsion chamber that is assembled with the vacuum-packed emulsion films and the frame-type spacers to fill the water as a target. The status and results of these analyses will be reported in the future. Then, a large-scale water target emulsion chamber will be constructed to study the neutrino water interactions in detail. 6. Conclusion We proposed a new experimental project with a nuclear emulsion detector to study low-energy neutrino nucleus interactions in detail for future neutrino oscillation physics experiments. The J-PARC T60 experiment has been implemented as a first step of this program to check the feasibility and confirm the neutrino event analysis. In this work, the systematic neutrino event analysis with full ECC scanning data has been developed and the first neutrino event detection has been successfully demonstrated. Moreover, the hybrid analysis with INGRID has been performed for the first time. We will polish the systematic neutrino analysis chain, develop the water target emulsion chamber and steadily expand the detector scale to reach the desired physics output. Acknowledgment We appreciate the support provided by the J-PARC Neutrino Group / Accelerator Group and the T2K Collaboration. We acknowledge the support from the Japan Society for the Promotion of Science (JSPS) through their grants (JSPS KAKENHI Grant Number , , , , , , , 16H00873). References [1] K. Abe et al. (T2K Collaboration), Phys. Rev. D 88, (2013) [2] K. Abe et al. (Hyper-Kamiokande Proto-Collaboration), PTEP 2015 (2015) 053C02. [3] P. Adamson et al. (NovA Collaboration), Phys. Rev. Lett. 116 (2016) [4] R. Acciarri at al. (DUNE Collaboration), arxiv: [physics.ins-det]. [5] F. An et al. (JUNO Collaboration), arxiv: [physics.ins-det]. 25/26

26 [6] S. B. Kim, arxiv: [physics.ins-det]. [7] M. Athar et al. (INO Collaboration), ino/openreports/inoreport.pdf. [8] M. G. Aartsen et al. (IceCube PINGU Collaboration), arxiv: [physics.ins-det]. [9] S. Adrian-Martinez et al. (KM3NeT Collaboration), arxiv: [astro-ph.im]. [10] A. A. Aguilar-Arevalo et al. (MiniBooNE Collaboration), Phys. Rev. Lett. 110 (2013) [11] A. Aguilar et al. (LSND Collaboration), Phys. Rev. D 64 (2001) [12] C. F. Powell et al., Nature 159 (1947) 694. [13] K. Niu et al., Prog. Theor. Phys. 46 (1971) [14] K. Kodama et al. (DONUT Collaboration), Phys. Lett. B 504 (2001) 218. [15] N. Agafonova et al. (OPERA Collaboration), Phys. Rev. Lett. 115 (2015) [16] T. Fukuda et al., 2013 JINST 8 P01023; T. Fukuda et al., 2014 JINST 9 P [17] S. Aoki et al., KEK/J-PARC-PAC ; S. Aoki et al., JPS Conf. Proc. 8, (2015). [18] A. L. Lead et al., Phys. Rev. D 19 (1979) [19] N. Ushida et al. (Fermilab E531 Collaboration), Nucl. Instrum. Meth. A 224 (1984) 50. [20] E. Eskut et al. (CHORUS Collaboration), Nucl. Instrum. Meth. A 401 (1997) 7. [21] K. Niwa, K. Hoshino and K. Niu, Proc. Int. Cosmic Ray Symp. High Energy Phenomena, Tokyo, p. 149 (1974); S. Aoki et al., Nucl. Instrum. Meth. B 51 (1990) 466; T. Nakano, Ph.D. Thesis, Nagoya University, Japan, 1997 (in Japanese); T. Nakano and K. Morishima, J. SPSTJ. 72 (2008) 229; K. Morishima and T. Nakano, 2010 JINST 5 P [22] M. Yoshimoto, T. Nakano, R. Komatani and H. Kawahara, arxiv: [physics.ins-det]. [23] N. Naganawa and K. Kuwabara, Abstract for SPIJ s annual Conference (2011) (in Japanese). [24] K. Kodama et al., Advances in Space Research 37 (2006) [25] S. Takahashi et al., Nucl. Instrum. Meth. A 620 (2010) 192. [26] H. Rokujo et al., Nucl. Instrum. Meth. A 701 (2013) 127. [27] K. Yamada et al. (J-PARC T60 Collaboration), arxiv: [physics.ins-det]. [28] K. Abe et al. (T2K Collaboration), Nucl. Instrum. Meth. A 694 (2012) 211. [29] T. Fukuda et al., 2010 JINST 5 P [30] N. Agafonova et al. (OPERA Collaboration), New J. Phys. 14 (2012) [31] T. Fukuda, OPERA Public Note, No. 179 (2015), OPERA note 179.pdf; [32] T. Toshito et al., Nucl. Instrum. Methods A 516 (2004) /26