arxiv: v1 [physics.ins-det] 8 Oct 2015

Transcription

1 Preprint typeset in JINST style - HYPER VERSION New High-Resolution Gadolinium-GEM Neutron Detectors for the NMX Instrument at ESS arxiv:11.236v1 [physics.ins-det] 8 Oct 21 Dorothea Pfeiffer a,b, Filippo Resnati a,b, Jens Birch c, Maddi Etxegarai a, Richard Hall-Wilton a,d, Carina Höglund a,c, Lars Hultman d, Isabel Llamas-Jansa a,e, Eraldo Oliveri b, Esko Oksanen a, Linda Robinson a, Leszek Ropelewski b, Susann Schmidt a,c, Christina Streli f, Patrik Thuiner b, f a European Spallation Source (ESS AB) P.O. Box 176, SE-221 Lund, Sweden b CERN CH-1211 Geneva 23, Switzerland c Linköping University IFM SE Linköping, Sweden d Mid-Sweden University SE-817 Sundsvall, Sweden e Institute for Energy Technology IFE NO-27 Kjeller, Norway f Vienna University of Technology 14 Vienna, Austria Dorothea.Pfeiffer@cern.ch ABSTRACT: European Spallation Source instruments like the macromolecular diffractometer, NMX, require an excellent neutron detection efficiency, high-rate capabilities, time resolution, and an unprecedented spatial resolution in the order of a few hundred micrometers over a wide angular range of the incoming neutrons. For these instruments solid converters in combination with Micro Pattern Gaseous Detectors (MPGDs) are a promising option. A GEM detector with gadolinium converter was tested on a cold neutron beam at the IFE research reactor in Norway. The µtpc analysis, proven to improve the spatial resolution in the case of 1 B converters, is extended to gadolinium based detectors. For the first time, a Gd-GEM was successfully operated to detect neutrons with an estimated efficiency of 1% at a wavelength of 2 Å and a position resolution better than 3 µm. KEYWORDS: Neutron detection; Solid converters; gadolinium; Gd; 1 B 4 C; Gas detectors; GEM; Position resolution; TPC. Corresponding author.

2 Contents 1. Introduction 1 2. Gadolinium as neutron converter 2 3. Experimental setup Beam line and detector setup Detector readout and DAQ 4. µtpc analysis. Conclusions and Outlook 1 6. Acknowledgements Introduction The European Spallation Source (ESS) [1] in Lund/Sweden is foreseen to start operations in 219 and will become the world s most powerful thermal neutron source. Its brightness will be significantly higher than the brightness at existing reactor sources like the Institut Laue-Langevin (ILL) [2] or other spallation sources like the Spallation Neutron Source (SNS) [3] and the Japan Proton Accelerator Research Complex (J-PARC) [4]. 22 neutron scattering instruments are currently planned as the baseline suite for the facility []. For each of the instruments efficient thermal neutron detectors are a crucial component [6]. Traditionally the dominant detector technology at neutron scattering instruments has been gaseous 3 He detectors [7]. In light of the 3 He crisis [8, 9] an extensive international R&D program is currently under way in order to develop efficient and cost-effective detectors based on other isotopes [1, 11, 12]. Solutions like gas proportional counters surrounded by thin films of 1 B-enriched boron carbide [13, 14, 1] have been developed for instruments with large area detectors. The ESS Macromolecular Diffractometer (NMX) [] requires three 6 x 6 cm 2 detectors with at least 3% detection efficiency and ca 2 µm spatial resolution. At reactor sources macromolecular crystallography instruments typically use neutron image plates [16] with this resolution, but spallation source instruments require time resolution that the image plates lack altogether. Scintillation based detectors [17, 18] are currently limited to ca 1 mm spatial resolution. For spallation source instruments solid neutron converters in combination with Micro Pattern Gaseous Detectors (MPGDs) [19] are a promising option to achieve the spatial resolution of the neutron image plate with time resolution, high-rate capabilities and a good neutron detection efficiency. A data analysis technique based on the principle of the Time Projection Chamber (TPC) 1

3 was developed, and has already been successfully applied to MPGDs with 1 B 4 C converters. A position resolution better than 2 µm was obtained by determining the start of the ionization track [2], but the detection efficiency was lower than %. This so-called µtpc analysis [21] is now extended to the case of Gas Electron Multipliers (GEM) detectors with more efficient gadolinium converters. This work presents the measurements carried out at the R2D2 beam line at the JEEP II research reactor at the Institute for Energy Technology (IFE) in Kjeller/Norway [22]. It consists of a description of the design, the setup, the measurement and the data analysis that was performed in order to retrieve the full information on the position of the neutron conversion. 2. Gadolinium as neutron converter Despite the very large neutron capture cross section of 1 Gd and 17 Gd, the material is not a popular converter due to the nature and the energy of the secondary particles. In fact, after the neutron capture, gadolinium releases prompt γs with an energy of up to 9 MeV and conversion electrons with energies ranging from 29 kev to 2 kev. Figure 1 shows the simulated kinetic energy distribution of the conversion electrons after the neutron capture in a 2 µm thick natural Gd converter and upon arrival in the drift volume. Natural gadolinium contains 14.8% of 1 Gd and 1.6% of 17 Gd, the remainder are Gd isotopes without significant cross section for thermal neutrons. Since the conversion electrons loose energy between the neutron capture location and the beginning of the drift, the discrete energy distribution present in the converter smears out and turns into a more continuous spectrum in the drift. Conversion Electrons [%] Conversion Electrons [%] Kinetic Energy [kev] Kinetic Energy [kev] (a) In 2 µm thick Gd converter (b) Entering drift volume (backwards configuration) Figure 1: Kinetic energy of conversion electrons (bin size 1 kev) as simulated with Geant4. Figure 4a demonstrates that the Gd conversion electrons on average deposit only 7 kev in the drift space by ionizing the Ar/CO 2. This is considerably less energy compared to the secondaries from 3 He, 6 Li, and 1 B converters. In addition to the small amount of deposited energy by the conversion electrons, gadolinium is also a high Z material. Gadolinium based detectors are therefore more sensitive to gamma background, when only the signal amplitude is used for the signal discrimination. The gamma sensitivity and advanced discrimination techniques though will not be discussed in this paper. Further, the range of the conversion electrons of between 2 mm and cm in Ar/CO 2 7/3 seems to be in contrast with the spatial resolution requirements. In the past, a CsI 2

4 layer applied to the gadolinium was used as an electron emitter to localize the electron energy loss and, therefore, improve the spatial resolution [23, 24]. The disadvantage in this case is the very small amount of primary electrons and deposited energy that requires a detector with very large gain. In this work a different approach was followed: the conversion electrons from the neutron capture ionize a macroscopic portion of the active volume, and an offline analysis takes care of the reconstruction of the neutron interaction point from the three dimensional topology of the event. The challenge is here that in contrast to the α particles and Li ions created in the 1 B 4 C converter, the gadolinium conversion electrons have not only on average a far larger range, but also do not leave a straight ionization track. Figure 2 shows the typical curved track of a 7 kev electron starting at (,,) in Ar/CO 2 7/3 as simulated with Degrad [2]. 2 z [mm] 2 z [mm] z [mm] y [mm] 1 x [mm] x [mm] y [mm] (a) 3D (b) Projection zx (c) Projection zy Figure 2: Example track of 7 kev electron in Ar/CO 2 7/3 originating in (,,) as simulated with Degrad. The Gd and detector simulations were carried out with Geant4 [26] version 1.1 using the high-precision neutron data G4NDL4..PUB, the gamma level data PhotonEvaporation3.1 and the flag G4NEUTRONHP_USE_ONLY_PHOTONEVAPORATION. The wavelength of the neutrons in the simulation was set to 2 Å, and the used geometry is depicted in figure 3b. Given the large thickness of the cathode, the converter was used in backwards configuration, i.e. the neutron beam, impinging orthogonally to the detector, crossed the readout board and the GEMs before reaching the gadolinium cathode. In this way the conversion electrons do not need to traverse the entire gadolinium thickness in order to reach the active volume, which leads to a higher neutron detection efficiency. On the other hand the scattering of the neutrons in the material of the readout decreases the position resolution and the efficiency. To consider the scattering of the cold neutrons, the hydrogen atoms in the Kapton of the GEM foils and the FR4 of the readout board were replaced with the TS_H_of_Polyethylene atom. For this atom neutron scattering data based on the evaluated nuclear data files ENDF/B-VII exists. Geant4 will use this scattering data, if the G4NeutronHP- ThermalScattering model is registered as model for the neutron elastic physics process. The simulations of the detector under test indicate that due to the scattering only 8% of the 2 Å neutrons reach the converter, where all of them are captured. In 1% of the capture reactions 3

5 one or more conversion electrons are created, of which in 1% of the cases at least one reaches the drift space1. The expected efficiency in backwards mode including scattering effects is therefore.8 times.1 or 12%. 3. Experimental setup 3.1 Beam line and detector setup For the detector tests at JEEP II, the monochromator of the R2D2 beam line was optimized to deliver mono-energetic neutrons of 2 Å, with a maximum flux of 11 khz at a beam size of.2 x 2 cm2. A multi-slit Gd2 O3 Soller collimator that limits the angular divergence of the beam to 3 is inserted into the R2D2 shutter. Figure 3a shows the beam line behind the shutter. It was equipped with a beam monitor, two sets of movable collimating slits made of borated Al blades, and a calibrated 3 He tube. HV1: 4 V HV2: 33 V Gd sheet (2 um) 1 MΩ kω 1 MΩ conversion electron 1 MΩ kω 1 mm 2 mm 2 mm 2 mm 1 MΩ 1 MΩ 4 kω neutron 1 MΩ (a) Detector under test at R2D2 beam line. (b) Schematic representation of the detector. Figure 3: Measurement setup at R2D2 beam line with collimation slits and Triple-GEM detector in reflecting configuration. The detector under test was a standard 1 x 1 cm2 Triple-GEM [29] detector flushed with Ar/CO2 7/3 mixture at 1 l/h at room temperature and atmospheric pressure. It consisted of three standard GEM foils stacked at a distance of 2 mm from each other and powered with a resistor chain. The current through the resistor-divider was set to 73 µa in order to obtain an effective gain of about. The conversion volume or drift space was 1 mm long, and as cathode served a 2 µm thick sheet of natural gadolinium. The high voltage to the cathode was provided with a second power supply, and a drift field of 7 V/cm was chosen to avoid electron attachment to electronegative impurities and the subsequent loss of primary ionization electrons, while keeping the drift velocity smaller than 2. cm/µs [3]. 1 The values given in the literature for the number of conversion electrons per captured neutron and the neutron detection efficiency of Gd based neutron detectors vary considerably. A probability of 9.9% for the emission of conversion electrons after the capture of thermal neutrons on natural Gd [27] and a measured detection efficiency for thermal neutrons of 17% in backwards configuration [28] can be found. 4

6 The efficiency measurement of the detector was carried out using a.2 x 2 cm 2 collimated beam. First the 3 He tube was installed directly behind the second set of slits, and a rate of 11 khz was measured with the 3 He tube. After the removal of the 3 He tube, a Cd sheet was placed in front of the GEM detector, and a rate of 3.7 khz measured with the GEM detector. Since each neutron captured in cadmium releases about 9 MeV of gamma energy, the measured 3.7 khz rate is most likely higher than the actual background. Finally the Cd was removed, after which the GEM measured a rate of 4.7 khz. Apparently the cold neutron beam was dominated by a high amount of background. Assuming that the actual background is lower than the measured rate behind the Cd sheet, the measured efficiency of (4.7 khz khz) / 11 khz or 9.1% should be regarded as a lower limit of the real efficiency. Nevertheless, the measured efficiency comes reasonably close to the simulated efficiency of 12%. 3.2 Detector readout and DAQ The readout board was a Cartesian x/y strip readout [31] with 26 strips and 4 µm pitch in x and y direction. Four APV-2 [32] hybrid chips per readout were used to pre-amplify the signals. The waveforms were digitized with the Scalable Readout System (SRS) [33] and acquired with the ALICE DAQ system DATE [34] and the Automatic MOnitoRing Environment (AMORE) [3] software. Data was acquired in two different collimation settings: 9 events were acquired with a collimated beam of.3 x 1 cm 2 and 8 events with a collimated beam of 2 x 2 cm 2. As shown in figures 8b and 8c, the detector clearly recorded the shape of the beam. During the measurements with the.3 x 1 cm 2 beam, the detector was rotated by 7. With a Cd sheet in front of the detector, the beam image vanished, whereas a Pb sheet just led to a less sharp beam image due to the scattering of the neutrons. This proves in addition to the rate measurements that the detector is able to detect neutrons. About 1% of the events in both data sets were excluded from the analysis due to electrical problems like misfiring strips. The comparison of acquired and simulated spectrum (figure 4b) shows a similar shape. The conversion electron peaks at 29 kev and 39 kev are too small to be detected within the energy resolution of the detector. The 8.9 kev copper X-ray peak is caused by prompt gammas from the Gd neutron capture that interact with the copper in the GEM foils. It is visible but very small in the measured spectrum. 4. µtpc analysis The typical signals due to a neutron conversion are shown in figure. The figure depicts the raw waveforms from the y view on different strips, digitized every 2 ns. The signals on different strips are not synchronous because different drift depths require different drift times. During the drift in electric fields, electrons diffuse as a result of multiple collisions with the gas molecules, and the initially localized charge cloud spreads out. The diffusion depends on the gas mixture, the pressure and the electrical field. At normal temperature and pressure (NTP), the σ of the transversal and longitudinal diffusion amounted to around 16 µm over the 1 mm drift [36]. An ionizing particle crossing a portion of the active volume creates electrons and ions that move in opposite directions separated by the electric field. The speed of the conversion electrons is much larger than the electron drift velocity, therefore the primary charge can be assumed to be released instantaneously

7 Events [%] 3 Events [%] 1 simulated spectrum measured spectrum Deposited Energy [kev] (a) Simulated spectrum of conversion electrons Deposited Energy [kev] (b) Simulated and measured spectrum of all particles Figure 4: Simulated spectrum of conversion electrons in drift and comparison of simulated and measured spectrum in drift (bin size 1 kev). along their path. The electron cloud moves rigidly at a constant speed due to the uniform electric field. Neglecting the small amount of diffusion, the electron cloud preserves its original track shape along the drift. Knowing the instant when the primary charge is released (t ), one can extract the drift time and consequently the position along the drift of the ionizing event. In the present configuration this information is not accessible. Nevertheless, in the case of straight tracks from α particle or muons, the electrons from the end of the track always induce the signal before the electrons from the beginning of the track, unless the particle moves in a direction parallel to the cathode. This concept forms the base of the Time Projection Chamber (TPC) [37], which allows the identification of the beginning of the track even without knowing t. Conversion electrons on the other hand quite frequently propagate in a different fashion. They first move away from their origin in the converter on the cathode, but can later turn and move back towards the cathode. The x and y view of such a typical curved track as measured with the Gd-GEM setup is depicted in figure 6. Obviously the signal from the end of the conversion electron track does not always arrive first. But it still holds true that the electrons from the beginning of the track have the largest drift time. The modified µtpc analysis for gadolinium is thus only able to determine the start of a track. As shown in figure, the waveform of the signal from the gadolinium conversion electrons contains quite frequently several local maxima. This fact stems from the changes of direction and the resulting curved path that the conversion electrons travel in the drift volume. However with the α particles from 1 B 4 C [2] and their straight path, the waveforms just had one maximum. To identify the start of the track, it was thus sufficient to determine for each strip the time t where the rising flank of the pulse reaches % of its maximum amplitude (7a). The start of the track was then the strip with the largest t. For the Gd conversion electrons one reaches better results when using the time of the global maximum (7b) or the time when the falling flank falls below a set threshold (7c). The best results are obtained when the time of the last local maximum is used (7d). Interpreting the track as a point cloud composed of local maxima, an approach based on the Principal Component Analysis (PCA) [38] has been tested to reconstruct the whole track. It did not lead to an improvement in position resolution compared to the last local maximum algorithm, 6

8 (a) strip 1 (b) strip 2 (c) strip (d) strip 4 (e) strip (f) strip (g) strip 7 (h) strip 8 (i) strip 9 Figure : Typical time structure of the neutron conversion signal in different strips with the 2 ns time bins on the x-axis, and the charge in arbitrary units on the y-axis. time bin [2 ns] tstart charge [ADC counts] time bin [2 ns] tstart charge [ADC counts] start of track x strip [.4 mm] start of track y strip [.4 mm] 4 2 (a) x view (b) y view Figure 6: Typical electron track in the drift volume of the Triple-GEM detector with Gd cathode. and was discarded due to the considerable computational effort. Figure 8a shows the cathode with copper grid that was used for the measurement. The copper 7

9 charge [ADC counts] maximum charge [ADC counts] global maximum % of maximum time bin [2 ns] t (a) Discriminator rising flank time bin [2 ns] t (b) Global maximum charge [ADC counts] charge [ADC counts] local maximum 1 local maximum threshold time bin [2 ns] (c) Discriminator falling flank t time bin [2 ns] t 1 t 2 (d) Last local maximum Figure 7: Different algorithms to determine the strip of origin of the conversion electron y strip [.4 mm] counts y strip [.4 mm] counts x strip [.4 mm] x strip [.4 mm] (a) 2 µm thick Gd cathode (b) Hit distribution.3 x 1 cm 2 (c) Hit distribution 2 x 2 cm 2 Figure 8: Gd cathode with copper tape and resulting hit distributions with regions of interest. grid was added to stop the conversion electrons from escaping from the converter into the gas. It provides a sharp edge to measure the position resolution in addition to the edge of the collimated beam, which is widened by the scattering of the neutrons in readout board and GEM foils. Both types of edges can be seen in the plot of the hit distribution 8b. To determine the improvement in position resolution due to the µtpc analysis compared to the traditional center-of-mass approach, the regions with a large step in the hit distribution were studied. Figure 9 shows as example the analysis of the position resolution in region P3 of the.3 8

10 x 1 cm 2 collimated beam in figure 8b. The standard deviation σ of the position resolution is extracted from the fit of the complementary error function to the data. The reconstructed position obtained with the µtpc concept (last local maximum) is shown in red, the center-of-mass approach in blue. The center-of-mass approach represents the best possible result for a non-tpc method, and is thus shown for comparison purposes. The σ of the position resolution amounts to 32 µm for the µtpc concept and 22 µm for the center-of-mass approach. Since the copper tape did not stop all conversion electrons, an improvement of the position resolution is to be expected in presence of a sharper edge. The results of all analyzed regions can be found in table 1. The σ of the position resolution is everywhere better than 3 µm for the tape edges and the µtpc technique. The values for the beam edges are indicative of the divergence of the collimation system and the scattering of the beam. A σ of 1.7 mm for the position uncertainty due to the divergence of the beam has been determined with a Monte Carlo simulation. For the scattering of the beam in the detector under test, a σ of 17 µm has been obtained with Geant4. The position resolution measured with the edge of the beam is thus clearly dominated by the divergence of the neutron beam. Even here the µtpc technique is usually a factor two better then the center-of-mass approach. Table 2 shows the results of the measurements with the 2 x 2 cm 2 beam. Since the edges in the 2 x 2 cm 2 beam hit distribution in figure 8c are a combination of beam edge and copper tape edge, the resulting position resolution amounts to on average µm. counts 2 2 last-local-maximum fit local-maximum (σ=32 um +/- 21 um) center-of-gravity fit center-of-gravity (σ=22 um +/- 391 um) y strip [.4 mm] Figure 9: Region P3 (.3 x 1 cm 2 collimated beam): Distribution of the reconstructed y coordinate using a center-of-mass-based technique (in blue) and the µtpc analysis (in red). 9

11 last-local-maximum falling-flank global- maximum center-of-mass Position σ [µm] error [µm] σ [µm] error [µm] σ [µm] error [µm] σ [µm] error [µm] tape P1 248 (2) 229 (31) 266 (28) 887 (214) tape P2 341 (22) 43 (29) 489 (42) 313 (86) tape P3 32 (21) 371 (28) 312 (2) 22 (391) tape P4 242 (17) 16 (86) 288 (29) 11 (294) beam P 142 (46) 13 (41) 1621 (3) 472 (92) beam P (8) 1672 (1) 1847 (61) 448 (87) beam P7 23 () 1812 (48) 1869 (2) 1982 (133) beam P (72) 264 (62) 291 (74) 478 (791) beam P (4) 1746 (49) 1769 () 4126 (2) beam P (1) 226 (77) 2393 (112) - - * average tape 288 (2) 3 (44) 339 (31) 1891 (44) average beam 233 (64) 1778 () 1932 (68) 3387 (2) Table 1: Fit results for the.3 x 1 cm 2 collimated beam and the regions of interest depicted in figure 8b. * Fit not possible (no edge in data). last-local-maximum falling-flank global-maximum center-of-mass Position σ [µm] error [µm] σ [µm] error [µm] σ [µm] error [µm] σ [µm] error [µm] beam/tape P1 77 (37) 884 (68) 81 (6) 331 (81) beam/tape P2 472 (39) 697 (47) 44 (36) 31 (8) beam/tape P3 734 (6) 829 (83) 828 (6) 3341 (1) beam/tape P4 41 (42) 891 (7) 391 (84) 2617 (67) beam/tape P 396 (31) 66 (46) 467 (39) 6483 (682) * beam/tape P6 3 (2) 441 (36) 419 (3) 242 (97) beam/tape P7 387 (41) 681 (77) 88 (79) 66 (7123) * beam/tape P8 49 (72) 794 (76) 498 (6) 289 (743) average 49 (43) 723 (61) 73 (8) 386 (2318) Table 2: Fit results for the 2 x 2 cm 2 collimated beam and the regions of interest depicted in figure 8c. * Limit of fit reached.. Conclusions and Outlook After already having demonstrated the µtpc concept with thermal neutrons using 1 B 4 C as neutron converter [2], it is here for the first time applied to data from a Gd-based neutron detector. The results show that the σ of the intrinsic position resolution is better than 3 µm. Gd-GEM detectors in combination with the µtpc concept thus fulfill the required position resolution for NMX. As next step, the development of a readout made from material with a lower scattering cross section for cold neutrons is foreseen to improve the detection efficiency and reduce the impact of the scattering on position resolution. To determine the position resolution without influence of the divergence of the beam, measurements with a cadmium or gadolinium mask in front of the readout board are planned. Further, more systematic measurements of the neutron detection efficiency are needed to crosscheck the simulation results. Natural and enriched Gd converters of various thicknesses have to be evaluated to increase the detection efficiency of the detector. 1

12 6. Acknowledgements This work was partially funded by the EU Horizon 22 framework, BrightnESS project

13 References [1] European Spallation Source ESS ERIC, [2] Institut Laue-Langevin, [3] Spallation Neutron Source, [4] J-PARC, [] S. Peggs et al, ESS Technical Design Report, ESS (213). [6] O. Kirstein et al, Neutron Position Sensitive Detectors for the ESS,Proceedings from the 23rd International Workshop on Vertex Detectors, 1-19 September 214, Macha Lake, The Czech Republic, PoS(Vertex214)29. [7] A.-J. Dianoux, G. Lander (ed), Neutron Data Booklet, 2nd edition, ILL (23). [8] A. Cho, Helium-3 shortage could put freeze on low-temperature research, Science 326 (29) [9] R.T. Kouzes, The 3 He supply problem, Pacific Northwest National Laboratory, Richland, WA, U.S.A., Tech. Rep. PNNL-18388, April 29. [1] International Collaboration for the Development of Neutron Detectors, [11] K. Zeitelhack, Search for alternative techniques to helium-3 based detectors for neutron scattering applications, Neutron News 23 (212) [12] B. Guerard, R. Hall-Wilton, F. Murtas, Prospects in MPGDs development for neutron detection, Summary based on presentations during RD1 Academia-Industry Matching Event, CERN October 14-1, 213, RD1-NOTE-214-3, [13] C. Höglund et al, B 4 C thin films for neutron detection, Journal of Applied Physics 111 (212) [14] K. Andersen et al, 1 B multi-grid proportional gas counters for large area thermal neutron detectors, Nuclear Instruments and Methods in Physics Research A 72 (213) [1] T. Bigault et al, 1 B multi-grid proportional gas counters for large area thermal neutron detectors, Neutron News 23 (212) 2-2. [16] N. Niimura et al, An imaging plate neutron detector, Nuclear Instruments and Methods in Physics Research A 349 (1994) [17] T. Hosoya et al, Development of a new detector and DAQ systems for ibix, Nuclear Instruments and Methods in Physics Research A 6 (29) [18] L. Coates et al, The macromolecular neutron diffractometer (MaNDi) at the Spallation Neutron Source, Oak Ridge: enhanced optics design, high-resolution neutron detectors and simulated diffraction, Journal of Applied Crystallography 43 (21) [19] M. Titov and L. Ropelewski, Micro-pattern gaseous detector technologies and RD1 collaboration, Modern Physics Letters A 28 (213) [2] D. Pfeiffer et al., The µtpc method: improving the position resolution of neutron detectors based on MPGDs, Journal of Instrumentation 1 (21) P44. [21] G. Iakovidis, The Micromegas project for the ATLAS upgrade, Journal of Instrumentation 8 (213) C127. [22] Institute for Energy Technology IFE, 12

14 [23] C. Schulz, Entwicklung eines neuartigen Niederdruck-Detektors mit hoher Orts- und Flugzeitauflösung zum Nachweis thermischer Neutronen, PhD Thesis Freie Universität Berlin (1999). [24] B. Mindur et al., Performance of a Micro-Strip Gas Chamber for event wise, high rate thermal neutron detection with accurate 2D position determination, Journal of Instrumentation 9 (214) P124. [2] S. Biagi, DEGRAD, [26] S. Agostinelli et al, GEANT4: A Simulation toolkit, Nuclear Instruments and Methods in Physics Research A 6 (23) [27] A.A. Harms and G. McCormack, Isotopic Conversion in Gadolinium-Exposure Neutron Imaging, Nuclear Instruments and Methods in Physics Research 118 (1974) [28] G. Bruckner et al, Position sensitive detection of thermal neutrons with solid state detectors (Gd Si planar detectors), Nuclear Instruments and Methods in Physics Research A 424 (1999) [29] F. Sauli, GEM: A new concept for electron amplification in gas detectors, Nuclear Instruments and Methods in Physics Research A 386 (1997) 31. [3] A. Peisert and F. Sauli, Drift and Diffusion of Electrons in Gases: A Compilation, CERN yellow report CERN84-8 (1984). [31] S. Bachmann et al., High rate X-ray imaging using multi-gem detectors with a novel readout design, Nuclear Instruments and Methods in Physics Research A 478 (22) 14. [32] M. French et al, Design and Results from the APV2, a Deep Submicron CMOS Front-End Chip for the CMS Tracker, Nuclear Instruments and Methods in Physics Research A 466 (21) [33] S. Martoiu et al, Development of the scalable readout system for micro-pattern gas detectors and other applications, Journal of Instrumentation 8 (213) C31. [34] F. Carena et al, The ALICE data acquisition system, Nuclear Instruments and Methods in Physics Research A 741 (214) [3] F. Roukoutakis and B. von Haller, Commissioning of the ALICE-LHC online data quality monitoring framework, Nuclear Instruments and Methods in Physics Research A 63 (29) [36] F. Sauli, Gaseous Radiation Detectors. Fundamentals and Applications, Cambridge University Press (214) [37] J.N. Marx and D. Nygren, The Time Projection Chamber, Physics Today 1 (1978) [38] O. E. Ruiz et al., Parametric Curve Reconstruction from Point Clouds using Minimization Techniques, GRAPP/IVAPP, page SciTePress, (213) 13