Design and test of an extremely high resolution Timing Counter for the MEG II experiment: preliminary results

Similar documents
An extreme high resolution Timing Counter for the MEG experiment Upgrade

An extreme high resolution Timing Counter for the MEG Upgrade

Pixelated Positron Timing Counter with SiPM-readout Scintillator for MEG II experiment

THE TIMING COUNTER OF THE MEG EXPERIMENT: DESIGN AND COMMISSIONING (OR HOW TO BUILD YOUR OWN HIGH TIMING RESOLUTION DETECTOR )

Liquid Xenon Scintillation Detector with UV-SiPM Readout for MEG Upgrade

Scintillation Tile Hodoscope for the PANDA Barrel Time-Of-Flight Detector

Production and Development status of MPPC

A prototype of fine granularity lead-scintillating fiber calorimeter with imaging read-out

Solid State Photon-Counters

Silicon PhotoMultiplier Kits

CAEN Tools for Discovery

Studies of large dynamic range silicon photomultipliers for the CMS HCAL upgrade

Time Resolution Improvement of an Electromagnetic Calorimeter Based on Lead Tungstate Crystals

Systematic study of innovative hygroscopic and non-hygroscopic crystals with SiPM array readout

Beam test of the QMB6 calibration board and HBU0 prototype

The Scintillating Fibre Tracker for the LHCb Upgrade. DESY Joint Instrumentation Seminar

Realization and Test of the Engineering Prototype of the CALICE Tile Hadron Calorimeter

Tests of Timing Properties of Silicon Photomultipliers

Single Photoelectron timing resolution of SiPM

Sensors for precision timing HEP

arxiv: v1 [physics.ins-det] 1 Nov 2015

Large photocathode 20-inch PMT testing methods for the JUNO experiment

TORCH a large-area detector for high resolution time-of-flight

PoS(PhotoDet 2012)018

THE WaveDAQ SYSTEM FOR THE MEG II UPGRADE

Report from the 2015 AHCAL beam test at the SPS. Katja Krüger CALICE Collaboration Meeting MPP Munich 10 September 2015

A flexible FPGA based QDC and TDC for the HADES and the CBM calorimeters TWEPP 2016, Karlsruhe HADES CBM

A very brief review of recent SiPM developments

arxiv: v1 [physics.ins-det] 2 Dec 2008

SciFi A Large Scintillating Fibre Tracker for LHCb

li, o p a f th ed lv o v ti, N sca reb g s In tio, F, Z stitu e tests o e O v o d a eters sin u i P r th e d est sezio tefa ectro lity stem l su

The hybrid photon detectors for the LHCb-RICH counters

HAPD and Electronics Updates

The Silicon Pixel Detector (SPD) for the ALICE Experiment

Development of Ultra-High-Density (UHD) Silicon Photomultipliers with improved Detection Efficiency

Quick Report on Silicon G-APDs (a.k.a. Si-PM) studies. XIV SuperB General Meeting LNF - Frascati

Status of the CUORE Electronics and the LHCb RICH Upgrade photodetector chain

TIMING COUNTER: status report. Flavio Gatti - Lecce, Sept 23, INFN Genova, Pavia, Roma1

A new Scintillating Fibre Tracker for LHCb experiment

This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and

The field cage for a large TPC prototype

Progress Update FDC Prototype Test Stand Development Upcoming Work

TitleLarge strip RPCs for the LEPS2 TOF. Author(s) Chu, M.-L.; Chang, W.-C.; Chen, J.- Equipment (2014), 766:

Status of the Timing Detector Plastic+SiPM Readout Option

Study of Timing and Efficiency Properties of Multi-Anode Photomultipliers

Updates on the Central TOF System for the CLAS12 detector

Performance and aging of OPERA bakelite RPCs. A. Bertolin, R. Brugnera, F. Dal Corso, S. Dusini, A. Garfagnini, L. Stanco

Photo Multipliers Tubes characterization for WA105 experiment. Chiara Lastoria TAE Benasque 07/09/2016

CMS Note Mailing address: CMS CERN, CH-1211 GENEVA 23, Switzerland

Spectroscopy on Thick HgI 2 Detectors: A Comparison Between Planar and Pixelated Electrodes

PHOTOTUBE SCANNING SETUP AT THE UNIVERSITY OF MARYLAND. Doug Roberts U of Maryland, College Park

Current status of Hamamatsu Si detectors mainly for High Energy Physics Experiments

arxiv:hep-ex/ v1 27 Nov 2003

Precise Digital Integration of Fast Analogue Signals using a 12-bit Oscilloscope

The CALICE test beam programme

MPPC (multi-pixel photon counter)

Development of an Abort Gap Monitor for High-Energy Proton Rings *

Spatial Response of Photon Detectors used in the Focusing DIRC prototype

Performance of a double-metal n-on-n and a Czochralski silicon strip detector read out at LHC speeds

This work was supported by FINEP (Research and Projects Financing) under contract

Status of GEM-based Digital Hadron Calorimetry

Photodetector Testing Facilities at Nevis Labs & Barnard College. Reshmi Mukherjee Barnard College, Columbia University

PMT Gain & Resolution Measurements in High Magnetic Fields

Commissioning and Initial Performance of the Belle II itop PID Subdetector

SPE analysis of high efficiency PMTs for the DEAP-3600 dark matter detector

Imaging TOP (itop), Cosmic Ray Test Stand & PID Readout Update

Hamamatsu R1584 PMT Modifications

A fast and precise COME & KISS* QDC and TDC for diamond detectors and further applications

Report from the Tracking and Vertexing Group:

Application of Hamamatsu MPPC to T2K near neutrino detectors

The TORCH PMT: A close packing, multi-anode, long life MCP-PMT for Cherenkov applications

Design of a Gaussian Filter for the J-PARC E-14 Collaboration

Concept and operation of the high resolution gaseous micro-pixel detector Gossip

3 EXPERIMENTAL INVESTIGATIONS Caroline Robson. 3.1 Aims and Objectives. 3.2 Experimental Method Set Up of the Test Stand

The 20 inch MCP-PMT R&D in China

Atlas Pixel Replacement/Upgrade. Measurements on 3D sensors

Calibration of photomultiplier tubes for the large-angle beamstrahlung detector at CESR

Lifetime of MCP-PMTs

ARDESIA: an X-ray Spectroscopy detection system for synchrotron experiments based on arrays of Silicon Drift Detectors.

DAQ Systems in Hall A

Investigation of time-of-flight PET detectors with depth encoding

Software Tools for the Analysis of the Photocathode Response of Photomultiplier Vacuum Tubes

A pixel chip for tracking in ALICE and particle identification in LHCb

Sensors for the CMS High Granularity Calorimeter

Laser Beam Analyser Laser Diagnos c System. If you can measure it, you can control it!

First LHC Beams in ATLAS. Peter Krieger University of Toronto On behalf of the ATLAS Collaboration

Institute of Electrical and Electronics Engineers (IEEE)

MCP Upgrade: Transmission Line and Pore Importance

RX40_V1_0 Measurement Report F.Faccio

... A COMPUTER SYSTEM FOR MULTIPARAMETER PULSE HEIGHT ANALYSIS AND CONTROL*

Comparison Between DRS4 Chip-Based Boards and ADCs for a Flexible PET Electronics

Development at Jefferson Lab

Review of the CMS muon detector system

Advances in multi-pixel Geiger mode APDs (Silicon Photomultipliers).

Screen investigations for low energetic electron beams at PITZ

CMS Upgrade Activities

Comparison of SONY ILX511B CCD and Hamamatsu S10420 BT-CCD for VIS Spectroscopy

Wavelength selective electro-optic flip-flop

Performance of a DC GaAs photocathode gun for the Jefferson lab FEL

Drift Tubes as Muon Detectors for ILC

Transcription:

Preprint typeset in JINST style - HYPER VERSION Design and test of an extremely high resolution Timing Counter for the MEG II experiment: preliminary results arxiv:32.087v [physics.ins-det] 3 Dec 203 M. De Gerone a, F. Gatti a,b, W. Ootani c, Y. Uchiyama c, M. Nishimura c, S. Shirabe d, P.W. Cattaneo e, M. Rossella e a Istituto Nazionale di Fisica Nucleare, Sezione di Genova, Via Dodecaneso 33, 646, Genova (GE), Italy b Universitá degli Studi di Genova, Via Dodecaneso 33, 646, Genova, Italy, c International Center for Elementary Particle Physics, University of Tokyo 7-3- Hongo, Bunkyo-ku, Tokyo 3-0033, Japan d Department of physics, Kyushu University 6-0- Hakozaki, Higashi-ku, Fukuoka 82-858, Japan e Istituto Nazionale di Fisica Nucleare, Sezione di Pavia, Via Agostino Bassi, 6, 2700, Pavia (PV), Italy E-mail: degerone@ge.infn.it ABSTRACT: The design and tests of Timing Counter elements for the upgrade of the MEG experiment, MEG II, is presented. The detector is based on several small plates of scintillator with a Silicon PhotoMultipliers dual-side readout. The optimisation of the single counter elements (SiPMs, scintillators, geometry) is described. Moreover, the results obtained with a first prototype tested at the Beam Test Facility (BTF) of the INFN Laboratori Nazionali di Frascati (LNF) are presented. KEYWORDS: Photon detectors for UV, visible and IR photons (solid-state) (PIN diodes, APDs, Si-PMTs, CCDs, EBCCDs etc), Scintillators, scintillation and light emission processes (solid, gas and liquid scintillators), Timing detectors. Corresponding author.

Contents. Introduction: the MEG experiment 2. The Timing Counter upgrade 3. Single counter optimisation 3 3. SiPM comparison 4 3.2 Scintillator comparison 6 4. Beam test 7 4. Setup 7 4.2 Data analysis 7 4.3 Results 8 5. Conclusions 8. Introduction: the MEG experiment The MEG experiment has been running since 2008 at the Paul Scherrer Institut (Villigen, CH), looking for the µ! eg decay. The MEG collaboration recently published the results based on the analysis of data collected in the years 2009-20: BR(µ! eg)apple 5.7 0 3 @90% C.L. [2]. While the analysis of the 202-203 data is still ongoing, an upgrade program of the MEG experiment (MEG II) has unfolded since 202 [3], aiming to improve the experiment sensitivity by an order of magnitude, down to 5 0 4. In order to reach such a sensitivity, most of the current detectors have to be re-designed or modified. In this paper, the development of an extremely high resolution detector for the measurement of the positron timing is described in details. 2. The Timing Counter upgrade The MEG detector [4] is designed to measure with the highest possible resolution the kinematic variables that define the signature of the decay µ! eg. Photons are detected by a Liquid Xenon detector placed outside the magnetic spectrometer where positrons are reconstructed (see figure ). The spectrometer is made of a superconductive magnet, a set of Drift CHambers (DCH) and the Timing Counter (TC). The DCH system, together with the specially designed field provided by the COBRA magnet, measures positron energy and emission angle, while the purpose of the TC is to measure the positron time of impact. The current Timing Counter [5] is made of two identical arrays (placed inside the magnet upand down-stream the target position) of 5 scintillating bars (Bicron BC404), with 80 4 4cm 3

Figure. Schematic view of the MEG detector: side and front views. Figure 2. Picture of the current Timing Counter. PMT and scintillator bars lodged in a black plastic socket are visible. size arranged in a barrel-like shape (see figure 2). Each bar is read-out on both sides by a fine mesh PhotoMultiplier Tube (PMT, Hamamatsu R5924). Signals from PMTs are processed to be fed into the trigger and DAQ system. The Timing Counter has been running since 2008, showing good and stable time resolution of 65 ps [6]. Some issues suggest that the design of the detector has to be changed to increase the resolution: the PMT operation in high magnetic field and helium environment deteriorates the PMT transit time spread and gain, also using fine mesh PMTs; large size scintillator bars generate uncertainties on impact point reconstruction and spread 2

of the trajectories of the optical photons inside the scintillator itself; the large amount of material crossed by the positron in the TC bar prevents the use of hits beyond the first one. These problems originate from the usage of PMTs and large size scintillator bars. Thus, the natural choice is to increase the detector granularity and upgrading the read-out system, exploiting the recent development of fast high gain solid state detector like Silicon PhotoMultipliers (SiPMs). A detector consisting of many scintillator plates (from now on: pixel) with SiPM read-out allows to overcome the limitations of the current TC: magnetic field has no influence on SiPM operation; higher granularity results in smaller uncertainties from impact position measurement; thanks to the smaller amount of material along the positron trajectory, it is possible to take advantage of the information coming from all the pixels crossed by the particle. Figure 3. Schematic view of the new timing counter design. On the left: overview of the detector. On the right: detail of the single counter configuration. The last point is quite remarkable, because the time resolution is expected to improve as / p N hit, where N hit is the number of pixels crossed by the positron. Moreover, the small size of the single element results in a more flexible configuration of the detector, allowing the possibility to tailor the position and the density of the pixels along the detector. Also, very short rise time scintillator (like Bicron BC422, see section 3.2) can be used even in presence of a short attenuation length. Single counter good performances have already been proved [7, 0]. In the following, the research and development work on several prototypes in order to choose the best material and device is presented. 3. Single counter optimisation The optimisation of the single counter configuration started from a systematic study among the SiPMs and the scintillators available to compare the properties relevant for our application. 3

3. SiPM comparison Silicon Photomultipliers are good candidate for the new TC, thanks to their characteristics: good time resolution, quite high gain, compactness. We tested different models from Hamamatsu Photonics, Advansid, Ketek and SensL. All these devices share some features: they have a size 3 3 mm 3, in order to be easily coupled to few mm thick scintillator pixels, and a good sensitivity in the near ultraviolet range, in order to match common plastic scintillators emission spectra. The SiPM models under test are summarised in table. Table. Summary of SiPMs model tested. Manufacturer Model Type Note S0362-33-050C Conventional (Old) MPPC Ceramic package S093-050P surface mount Hamamatus Photonics S2572-050C(X) New (standard type) MPPC Metal quench resistor S2572-020C(X) 25 µm pitch S2652-050C(X) Trench-type MPPC Metal quench resistor 3X3MM50UMLCT-B Improved fill factor Advansid NUV type Ketek PM3350 prototype-a Trench Type SensL MicroFB-30050-SMT B-Type Fast output For each model, the noise level (dark count rate and cross-talk) together with the PDE has been evaluated. Moreover, also the breakdown dependence on temperature has been evaluated. Finally, the timing resolution has been measured on a prototype pixel with fixed sizes. Setup SiPMs are put in a thermal chamber, which allows to keep the device at fixed temperature (23 C in the following measurements). Signals are transmitted on a coaxial cable to a voltage amplifier (developed at PSI, based on MAR-6SM amplifier [7]), then they are sampled at 5 Gs/s by a waveform digitiser (DRS4 evaluation board [8, 9] also developed at PSI). Dark noise and cross-talk The noise level of the device is evaluated by looking at the waveforms acquired by a random trigger. The dark count rate is calculated from the probability of observing zero photo-electron P(Np.e.= 0) in a fixed time window. The result is shown in figure 4(a) as a function of the applied over-voltage. The cross-talk probability is calculated from the P(Np.e. 2)/P(Np.e. ) ratio including a correction for the accidental coincidence of dark pulses. The result is shown in figure 4(b). The cross-talk probability almost linearly increases with the over-voltage. The standard-type SiPMs, namely SiPMs without a trench structure, turned out to have worse performance with respect to the trench type whose improved structure strongly reduces the noise level. Anyway, the typical energy release in a pixel should guarantee an adequate signal-to-noise ratio also for those SiPMs with higher dark count and cross-talk rates. PDE The PDE for Near UltraViolet (NUV) light is measured with a LED whose wavelength (370 40 nm) approximately matches the scintillator emission peak. The LED intensity is adjusted in such a way that the average number of observed photo-electrons ranges between 0.5 to.0. The relative PDE is then calculated from P(Np.e.= 0) in accordance with Poisson statistics, 4

and thus the measured PDE value does not include the effect of cross-talk nor after-pulsing. The result is shown in figure 4(c). The highest PDE is obtained with Hamamatsu S2572 model, with 50 µm pitch cell. A more detailed description of noise and PDE studies can be found in [2]. ) 3 Dark count rate (Mcps/mm.2 0.8 Old MPPC Old MPPC SMD New MPPC (50µm) New MPPC (25µm) New MPPC (Trench) KETEK Prototype-A AdvanSiD NUV SensL B-series 0.6 0.4 0.2 0 0 2 4 6 8 0 Over voltage (V) (a) Cross-talk probability 0.9 0.8 0.7 0.6 0.5 Old MPPC Old MPPC SMD New MPPC (50µm) New MPPC (25µm) New MPPC (Trench) KETEK Prototype-A AdvanSiD NUV SensL B-series 0.4 0.3 0.2 0. 0 0 2 4 6 8 0 Over voltage (V) (b) Relative PDE.4.2 Resolution (ps) 00 90 80 70 Old MPPC New MPPC (50µ m) New MPPC (25µm) KETEK Prototype-A SensL B-series New MPPC (Trench) AdvanSiD NUV 0.8 0.6 0.4 0.2 Old MPPC Old MPPC SMD New MPPC (50µm) New MPPC (25µ m) New MPPC (Trench) KETEK Prototype-A 60 50 40 AdvanSiD NUV SensL B-series 0 0 2 4 6 8 0 Over voltage (V) 30 2 4 6 8 Over voltage/sipm (V) (c) (d) Figure 4. Summary of the SiPMs comparison. a) Dark Count. b) Cross Talk Probability c) PDE d) Time Resolution. All results are given as a function of the applied over-voltage. Breakdown voltage versus temperature dependence The BreakDown voltage (BD) versus temperature dependence has been measured by plotting the I-V characteristic of each SiPM at different temperatures (see figure 5) in the range 20 45, showing a linear dependence on temperature, Time resolution The basic setup for the timing resolution measurement is the same described above. A scintillator pixel with size 60 30 5 mm 3 is read-out on each side by an array of 3 SiPMs connected in series. SiPMs are coupled to the pixel with optical grease. A 35 ns coaxial cable (7.5 m) transports signals to amplifiers, simulating the final experimental conditions. Counters are excited by using a 90 Sr b-source, providing electrons with 2.2 MeV endpoint energy. An 5

I-V curve VS Temperature: Advansid Current [µa] 0 ADV 5 ADV 20 ADV 25 ADV 30 ADV 35 ADV 40 ADV 45-0 -2 0 24 26 28 30 Voltage [V] Figure 5. I-V curves for different temperatures acquired with Advansid NUV SiPM. external reference counter (RC) made of a small piece of scintillator (BC422, size: 5 5 5 mm 3 ) coupled to a Hamamatsu S0362-33-050C SiPM is used for triggering purposes. The timing is extracted by applying a software constant fraction discrimination on the recorded waveform with discriminating fraction in the range 5 0% depending on SiPM model. The time resolution of the system is evaluated as the width of the distribution DT = T re f (T 0 + T )/2, being T re f and T i the time measured by the reference counter and each SiPM array respectively. The summary of the results is showed in figure 4(d) as a function of the applied over-voltage. 3.2 Scintillator comparison Three types of ultra fast plastic scintillator from Saint-Gobain Crystals, BC48, BC420 and BC422, were tested. The main characteristics of each scintillator are summarised in table 2, where also the characteristics of the BC404 are listed. The test was performed using 60 30 5 mm 3 pixels. The best resolution is obtained with BC422, the one with the fastest rise time. Table 2. Summary of the properties of plastic scintillators tested. Measured time resolutions on a 60 30 5 mm 3 sample are also listed. Properties BC404 BC48 BC420 BC422 Light Yield (% Anthracene) 68 67 64 55 Rise Time (ns) 0.7 0.5 0.5 0.35 Decay time (ns).8.4.5.6 Wavelength peak (nm) 408 39 39 370 Attenuation length (cm) 40 00 0 8 Measured resolution (ps) - 48 5 43 6

4. Beam test In order to test the detector in experimental conditions similar to the final one and check the multiple hit scheme, a small prototype was built and tested at the Beam Test Facility (BTF) at the INFN Laboratori Nazionali di Frascati []. The BTF beam can be tuned in such a way to provide electrons with energy similar to the MEG signal (48 MeV in our test) with average bunch multiplicity lower than. We decide to test both Hamamatsu and Advansid counters, the ones with the best trade off between time resolution and temperature dependence. 4. Setup We prepared two sets of pixel prototypes with BC48 scintillator, with 90 40 5 mm 3 sizes, equipped with Hamamatsu S2572-050C(X) (8 counters) and Advansid NUV (6 counters) SiPMs. Pixels are wrapped with 3M Radiant Mirror Film. Pixels are mounted on a moving stage that controls the movement in the plane perpendicular to the beam. The whole system is mounted on an optical bench enclosed in a shielded black box. The same reference counter described in section 3. is placed along the beam trajectory in front of the pixels. A lead glass calorimeter is placed behind the pixels for beam monitoring. The whole system is aligned to the beam line by using a laser tracker. Signals from SiPMs are fed into six DRS4 evaluation boards and sampled at 2.5 Gs/s. 4.2 Data analysis Charge analysis Events are selected by cutting on the charge distribution of the first two pixels. An example of distribution is showed in figure 6, where the bunch multiplicity is clearly visible. Moreover, we applied also a cut on the reference counter charge spectrum, by selecting the events around the Landau peak of the charge distribution. Figure 6. Charge distribution of the first couple of pixels. The selected events (single electron bunch) are marked in red. 7

The timing resolution is then evaluated by taking the width of the DT distribution, defined in two different ways: DT = T ref N DT = N N Â i= N Â T j j= N T i, (4.) N Â i= T i, (4.2) where T ref and T i, j is the time measured by the reference counter and by the pixels respectively. In formula 4.2 the sum is made over two different subgroups of pixels. In both cases, we can evaluate the timing resolution as a function of the number of hits used in the time averaging. DRS calibrations Dedicated runs were taken to evaluate the contribution from the electronics jitter. It was found to be 8.7 ps and 6.2 ps for pixels whose arrays are read-out by the same or different boards, respectively. 4.3 Results We checked the multiple hit scheme, with the same method described in section 3. studying the time resolution versus the number of hits. As expected, the best result is obtained with the largest number of hits, with s(dt ) < 30 ps. Preliminary resolutions are summarised in figure 7, compared with the expected / p N hit behaviour, which is also shown. Figure 7. Summary of the obtained resolution as a function of the number of hit. The expected resolution and the estimated / p N hit behaviour are also showed. 5. Conclusions We presented the R&D work on the upgrade of the Timing Counter for the MEG II experiment. The basic concepts of the new design, namely the good time resolution achievable with small 8

scintillator counters read-out by SiPMs and the improvement of the overall time resolution by averaging the time measurements over multiple hits has been tested. Optimising the choice among different types of SiPM and scintillators leads to obtain extremely good time resolution with a single counter down to s(dt ) 43 ps. A beam test performed at the Beam Test Facility in Frascati proved experimentally the multiple hit scheme. Analysis is still ongoing, a prelimiary resolution s(dt ) < 30 ps with eight counters is measured. Acknowledgments The authors would like to thank the Beam Test Facility crew, the mechanical and electronics workshops at INFN Section of Genova the Paul Scherrer Institute detector group for their valuable help. References [] G. Blankenburg et al., Neutrino masses and LFV from minimal breaking of U(3) 5 and U(2) 5 flavor symmetries, [arxiv:204.0688v2] [hep-ph]. [2] J. Adam et al., [MEG Collaboration], New Constraint on the Existence of the µ +! e + g Decay, Phys.Rev.Lett. 0 (203) 2080. [3] A.M. Baldini et al., [MEG Collaboration], MEG upgrade proposal, [arxiv:30.7225] [physics.ins-det]. [4] J. Adam et al., [MEG Collaboration], The MEG detector for µ +! e + g decay search, Eur. Phys. J. C 73 (203) 2365. [5] M. De Gerone et al., Development and Commissioning of the Timing Counter for the MEG Experiment, IEEE Trans. Nucl. Sci. 59 (202) 379. [6] M. De Gerone et al., The MEG timing counter calibration and performance, Nucl. Inst. Meth. A 638 (20) 4. [7] A. Stoykov et al., A time resolution study with a plastic scintillator read out by a Geiger-mode Avalanche Photodiode, Nucl. Inst. Meth. A, 695 (202) 202. [8] http://www.psi.ch/drs/evaluation-board [9] S. Ritt et al., Application of the DRS Chip for Fast Waveform Digitizing, Nucl. Inst. Meth. A, 623 (200) 486. [0] W. Ootani, Development of Pixelated Scintillation Detector for Highly Precise Time Measurement in MEG Upgrade, Nucl. Inst. Meth. A, 732 (203) 46. [] G. Mazzitelli et al., Commissioning of the DAFNE beam test facility, Nucl. Inst. Meth. A, 55 (2003) 524. [2] Y. Uchiyama [MEG collaboration], Nuclear Science Symposium Conference Record, IEEE, Seoul, Korea, 203., in press 9