Advances in High Energy Physics

Volume 2013 (2013), Article ID 256230, 7 pages

http://dx.doi.org/10.1155/2013/256230

Research Article

A Mobile Detector for Muon Measurements Based on Two Different Techniques

B. Mitrica,¹ D. Stanca,¹ M. Petcu,¹ I. M. Brancus,¹ R. Margineanu,¹ A. Apostu,¹ C. Gomoiu,¹ A. Saftoiu,¹ G. Toma,¹ H. Rebel,² A. Haungs,² O. Sima,³ A. Gherghel-Lascu,¹ and M. Niculescu-Oglinzanu¹

¹Horia Hulubei National Institute for Physics and Nuclear Engineering, 077125 Magurele, Romania

²Institute of Experimental Nuclear Physics, Karlsruhe Institute of Technology-Campus North, 76021 Karlsruhe, Germany

³Department of Physics, University of Bucharest, 077125 Magurele, Romania

Received 25 October 2012; Accepted 17 February 2013

Academic Editor: Maury Goodman

Copyright © 2013 B. Mitrica et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Precise measurements of the muon flux are important for different practical applications, both in environmental studies and for the estimation of the water equivalent depths of underground sites. A mobile detector for cosmic muon flux measurements has been set up at IFIN-HH, Romania. The device is used to measure the muon flux on different locations at the surface and underground. Its first configuration, not used in the present, has been composed of two 1 m² scintillator plates, each viewed by wave length shifters and read out by two Photomultiplier Tubes (PMTs). A more recent configuration, consists of two 1 m² detection layers, each one including four 1 · 0,25 m² large scintillator plates. The light output in each plate is collected by twelve optical fibers and then read out by one PMT. Comparative results were obtained with both configurations.

1. Introduction

The knowledge of cosmic ray muon flux is important in many applications, from low background measurements and the estimation of water equivalent depth for different underground locations to the characterization of the background in the vicinity of mega detectors for neutrino physics, like Super-Kamiokande or the future LAGUNA project [1], also for numerous environmental experiments like the solar activity characterization or the weather and climate change surveillance [2].

Also, it has successfully been used as muon tomography technique in the search for hidden rooms in pyramids [3] or in volcanology [4]. It has a long list of possible applications, like the detection of unknown caverns, to increase the safety procedures in mining excavations, in the oil industry as an easy way to search for oil bags or at the customs checkpoints, to scan the passing vehicles.

A mobile detector for measuring the muon flux is in use at IFIN-HH Bucharest [5]. Placed in a van for mobility, it interprets the coincidence events between two overlaid active layers as passing muons. A new configuration of the active layers, using optical fiber readout of the plastic scintillator sheet, was implemented and tested, and its performances were compared with the previous one, which used waveguides for the light collection. Using this new developed configuration, the possibility of directional measurements of muon flux was observed.

2. The Apparatus

The detection principle of the muon mobile detector is to identify the traversing muons as coincidence events between two overlaid 1 m² active layers of scintillator material. Used for the current detector, two different concepts for the active layer are presented in Figure 1. One important property of this device is its mobility, the device being conceived to measure the muon flux at different locations at the surface of the ground or in underground.

(a)

(b)

Figure 1: The detection module of (a) the first detector setup. Design of KASCADE [3]; (b) the new detector setup [5].

For the former configuration of the mobile detector, two detection modules were used, each module being a scintillator plate (NE114 type) of 0.9025 m² and 3 cm thickness. Divided into four parts ( m²) [6, 7], it is readout through a wavelength shifter (NE174A type) by two photomultiplier tubes (EMI 9902 type), as we can see in Figure 1(a).

The present configuration consists of two detection modules, each composed of 4 plastic scintillator sheets (Polystyrol 80%, Methylmethacrylate 20%) having  cm³. Every sheet is crossed by 13 longitudinal strips, 12 of them being filled with an optical fiber (Figure 1(b)). The light signal of each sheet is readout by a PMT. The modules are arranged on top of each other (at 30 cm distance). The signals from the eight photomultiplier tubes are OR-ed four by four (1 + 2 + 3 + 4) and (5 + 6 + 7 + 8) and then measured in coincidence using a gate of 50 ns, so no correction due to the dead time of the detector is necessary. A counter module registers the coincidence events (see Figure 2).

Figure 2: The schematic view of the electronic detection system.

The calibration has been made by comparing the measured energy deposit spectrum of the minimum ionizing particles with GEANT simulated responses. The comparison of the measured and simulated spectra, of both configurations, is shown in Figure 3. The difference in the most probable energy deposit between the two configurations comes from the different thicknesses of the plates (3 cm versus 1 cm).

(a)

(b)

Figure 3: (a) The energy calibration for one scintillator plate of the previous configuration, made by comparing the measured energy deposit spectrum of the minimum ionizing particles with GEANT 3.21 [28] simulated one. (b) The energy calibration for one scintillator sheet of the recent detector setup, made by comparison of the measured energy deposit spectrum of the minimum ionizing particles with GEANT 4.94 [29] simulated one.

For the first configuration, having the most probable energy deposit at 6.3 MeV, the signal threshold was set to 2.1 MeV [5]. For the new configuration, with the most probable energy deposit of 2.4 MeV, a signal threshold of 1.8 MeV is set.

Considering the fact that not all muons that interact with the first layer manage to pass to the second one, due to scattering loses in between, the acceptance of the detector was also investigated, using GEANT4 code (see Figure 4). The muon flux at the surface of the ground was estimated using CORSIKA code [8], using the primary cosmic ray spectrum of proton and helium, obtained from AMS results during a space shuttle mission [9].

Figure 4: The schematically view of the muons interacting with the new detection system, observed with GEANT 4.94 code (magenta lines represent muons). We can see that not all muons that cross the first scintillator layer will pass through the second one.

A correction factor has to be applied on the observed muon rate at ground level. The distance between the centers of the two active layers of the detector is 30 cm. For a distance of 8 cm, obtained if the two active layers are overlapping, a correction factor is obtained. No correction factor due to the geometry of the detector’s surface is required. The corrected muon flux is given by

where is the counting rate of the detector and   is the real flux.

3. Measurements and Results

Using the mobility of the detector, measurements of cosmic muon flux have been performed for both configurations of the active layers at different locations on the surface with different elevations. The results, presented in Table 1, are compared with those previously reported in [5] and displayed in Figure 5.

Table 1: Cosmic muon flux measured at different geographic locations with different elevations, using both configurations of the active layers. The altitude and the coordinates were measured with a GPS system.

Latitude (°)
Longitude (°)
Altitude (m a.s.l.)
Muon flux (m⁻² s⁻¹)

Former configuration
44.36
28.05

119.1 ± 3.6

44.40
26.10

122.3 ± 3.7

44.32
28.19

128.1 ± 3.8

45.24
25.94

143.2 ± 4.3

45.28
25.97

145.3 ± 4.4

45.29
25.94

146.7 ± 4.4

Recent configuration
44.97
26.02

129.1 ± 4.5

45.07
26.03

136.7 ± 4.7

45.24
25.94

142.7 ± 4.9

45.26
25.96

143.4 ± 4.9

45.26
25.97

144.6 ± 4.9

45.48
25.91

146.8 ± 5.0

45.48
25.90

150.3 ± 5.1

45.36
25.52

155.9 ± 5.3

Figure 5: Measured results of the muon flux variation with altitude, obtained with both detector configurations of the active layers of the detector. The circles represent the data obtained with the former configuration. The data obtained with the new configuration is presented with triangles. The cross point represents the results from [30].

The acquisition time for each data set was 1 h, all runs being performed at approximately the same time of the day (noon) in order to reduce the eventual influence of the solar activity and of atmospheric conditions.

Measurements of the cosmic muon flux have been performed also underground in the Unirea mine from Slanic Prahova (−208 m depth from surface), where a low-radiation level laboratory of IFIN-HH is set up, using both configurations. The detection modules were removed from the car and transported by an elevator to the observation level.

The responses are similar, 0.18 ± 0.01 m⁻² s⁻¹ for the first configuration and 0.19 ± 0.002 m⁻² s⁻¹ for the second one. The acquisition time was 1 hour for the first set of measurements and 10 hours for the second one.

The estimation of the water equivalent depth from muon flux data is made by the formula

with  m⁻² s⁻¹ being an overall normalization constant, mwe representing an effective attenuation length for high-energy muons, and being a free parameter [9]. The variation of the muon flux with the water equivalent depth is still an open problem and one of our detector’s goals. The formula has been obtained by Monte Carlo simulation using the MUSIC code [10].

A value of mwe was obtained using the flux measured with the first configuration; from the second one, a mwe is being calculated. We can see the evolution of the muon flux as a function of mwe depth (see Figure 6). The data are compared with others over the world (see Table 2).

Table 2: A compilation of different muon flux measurements over the world.

Location (laboratory)/country
Latitude (°)
Longitude (°)
Depth (mwe)
Muon flux (m⁻² s⁻¹)

CUPP/Finland [13]
63.39 N
26.02 E
0

Belgrade/Serbia [14]
44.51 N
20.23 E
0

Slanic/Romania [5]
45.24 N
25.94 E
0

Belgrade/Serbia
44.51 N
20.23 E
25
45

Belgrade/Serbia
44.51 N
20.23 E
75
11.9

Belgrade/Serbia
44.51 N
20.23 E
100
7.24

Ogoya/Japan [15]
36.20 N
136.5 E
270

Belgrade/Serbia
44.51 N
20.23 E
300

Slanic/Romania
45.24 N
25.94 E
601

Slanic/Romania
45.24 N
25.94 E
610

Slanic/Romania
45.24 N
25.94 E
790

CUPP/Finland
63.39 N
26.02 E
980

Solotwina/Ukraine [16]
50.11 N
23.16 E
1000

LSC/Spain [17]
42.43 N
00.31 E
1200

Oto-Cosmo/Japan [18]
34.41 N
135.5 E
1400

Carlsbad/USA [19]
33.09 N
117.2 W
1526

CUPP/Finland
63.39 N
26.02 E
1900

Y2L/Korea [20]
38.04 N
128.4 E
2000

Soudan/USA [21]
47.49 N
92.14 W
2000

LSC/Spain
42.43 N
00.31 E
2400

Kamioka/Japan [22]
39.29 N
140.3 E
2700

Boulby/UK [23]
54.34 N
00.58 W
2800

CUPP/Finland
63.39 N
26.02 E
2810

Gran Sasso/Italy [24]
42.48 N
13.33 E
3200

CUPP/Finland
63.39 N
26.02 E
3960

BNO/Russia [25]
43.40 N
43.32 E
4700

Frejus/France [26]
45.08 N
06.41 E
4800

Sudbury/USA [27]
46.28 N
81.10 W
6010

Figure 6: The variation of flux as function of mwe depth for the results presented in Table 2.

In order to check if there are any preferential arrival directions of the muons, coincidence measurements have been performed in the underground between the scintillator sheets of the detector setup. The measurements have been performed in Unirea mine from Slanic Prahova salt ore, close to one of the mine’s wall (see Figure 7). The “Unirea” mine is characterised by temperature: 12.0–13.0°C, humidity: 65–70%, excavated volume:  m³, floor area:  m², average high: 52–57 m, aerosols <10 mm: 2108 part/m³, distance between walls: 32–36 m, and existing infrastructures: electricity, roads, railway, elevator, phone, Internet, and GSM networks (also inside the galleries). More details regarding the radioactive background in the mine are presented in [11].

Figure 7: Map of Unirea mine, with the mBq laboratory of IFIN-HH.

The order of the scintillator sheets of the detector is shown in Figure 8. The results are presented in Table 3. The system was positioned in the cavern, so the sheets 1 and 5 are to be placed near the wall of the mine and sheets 4 and 8 toward the empty space of the cavern.

Table 3: Rate of coincidences between sheets. Measurements were performed in the Unirea mine (−208 m depth from surface). The acquisition time for each measurement was set to 10 hours. The coincidence is represented by “” sign.

Channels
Counting rate (p/s)

(1 + 2 + 3 + 4) (5 + 6 + 7 + 8)

1 + 6

4 + 7

2 + 5

3 + 8

2 + 7

3 + 6

(1 + 2) (7 + 8)

(3 + 4) (5 + 6)

Figure 8: The orientation of the scintillator sheets of the new detector setup.

We can observe that the coincidence rate between 3 + 4 and 5 + 6 is 3 times higher than the coincidence rate of 1 + 2 and 7 + 8, indicating the fact that the muons which are arriving at the detector from the direction of the cavern are more numerous than those which are coming from the direction of the wall. This is an important observation because in this way we can determine the variations of the density in different materials or in the terrestrial crust, with many practical applications, like the detection of oil deposits, hidden rooms in pyramids, or at customs checkpoints in the search of contraband products.

4. Conclusions

A versatile detector for cosmic muon flux measurements was set up and studied at IFIN-HH, Bucharest. It can be used in many applications where the information about the flux of atmospheric muons is of interest. One of them that were demonstrated in this paper is the determination of water equivalent depths in Unirea mine of Slanic [9] and could be used to measure the water equivalent depth on any site in order to check the possibility to install a large underground experiment like LAGUNA, MINOS [12], and so forth. The results show that the mobile detector could be used to measure the directional variation of the muon flux in order to study hadronic interaction models or muon multiplicities.

In this work, the technique of collecting light signals by PMTs through optical fibers was investigated. Its reliability was demonstrated by the good agreement between the results obtained with this configuration and another one, where the light is transferred from the scintillator plate to a PMT through a light guide, in the determination of the cosmic muon flux in underground locations (Unirea mine from Slanic Prahova) and the variation of muon flux with the altitude. This type of reading technique combined with the mobility of the detector could be used easily in various applications, like archaeology and volcanology, or as a scanner used in customs control.

Acknowledgment

The Romanian authors would like to thank the support of the Romanian Authority for Scientific Research by the Projects Parteneriate 194/2012—CORONA and PN 09 37 01 05.

References

A. Rubbia, “The LAGUNA Design Study: towards giant liquid based underground detectors for neutrino physics and astrophysics and proton decay searches,” Acta Physica Polonica B, vol. 41, pp. 1727–1732, 2010.

A. Saftoiu, A. Bercuci, I. M. Brancus, et al., “Measurements of the cosmic muon flux with the willi detector as a source of information about solar events,” Romanian Journal of Physics, vol. 56, no. 5-6, pp. 664–672, 2011.

J. Marteau, D. Gibert, N. Lesparre, F. Nicollin, P. Noli, and F. Giacoppo, “Muons tomography applied to geosciences and volcanology,” Nuclear Instruments and Methods in Physics Research, vol. 695, no. 11, pp. 23–28, 2012.

S. Aguillar, et al., “Searching for cavities in the Teotihuacan Pyramid of the Sun using cosmic muons,” in Proceedings of the 32nd International Cosmic Ray Conference (ICRC '11), Beijing, China, 2011.

B. Mitrica, R. Margineanu, S. Stoica, et al., “A mobile detector for measurements of the atmospheric muon flux in underground sites,” Nuclear Instruments and Methods in Physics Research A, vol. 654, pp. 176–183, 2011.

H. Bozdog, M. Elzerb, H.J. Gils, et al., “The detector system for measurement of multiple cosmic muons in the central detector of KASCADE,” Nuclear Instruments and Methods in Physics Research A, vol. 465, pp. 455–471, 2001.

T. Antoni, et al., “The cosmic-ray experiment KASCADE,” Nuclear Instruments and Methods in Physics Research A, vol. 513, pp. 490–510, 2003.

D. Heck, J. Knapp, J. N. Capdevielle, G. Schatz, and T. Thouw, “CORSIKA: a Monte Carlo code to simulate extensive air showers,” Report FZKA 6019, Forschungszentrum, Karlsruhe, Germany, 1998.

J. A. Formaggio and C. J. Martoff, “Backgrounds to sensitive experiments underground,” Annual Review of Nuclear and Particle Science, vol. 54, pp. 361–412, 2004.

V. A. Kudryavtsev, “Muon simulation codes MUSIC and MUSUN for underground physics,” Computer Physics Communications, vol. 180, no. 3, pp. 339–346, 2009.

R. Margineanu, C. Simion, S. Bercea et al., “The Slanic-Prahova (ROMANIA) underground low-background radiation laboratory,” Applied Radiation and Isotopes, vol. 66, no. 10, pp. 1501–1506, 2008.

P. Schreiner, et al., “Interpretation of the atmospheric muon charge ratio in MINOS,” in Proceedings of the30th International Cosmic Ray Conference (ICRC '07), vol. 9, pp. 97–98, Merida, Mexico, 2007.

T. Enqvist, A. Mattila, V. Föhr, et al., “Measurements of muon flux in the Pyhäsalmi underground laboratory,” Nuclear Instruments and Methods in Physics Research A, vol. 554, pp. 286–290, 2005.

A. Dragic, D. Jokovića, R. Banjanaca, et al., “Measurement of cosmic ray muon flux in the Belgrade ground level and underground laboratories,” Nuclear Instruments and Methods in Physics Research A, vol. 591, pp. 470–475, 2008.

K. Komura and Y. Hamajima, “Ogoya underground laboratory for the measurement of extremely low levels of environmental radioactivity: review of recent projects carried out at OUL,” Applied Radiation and Isotopes, vol. 61, pp. 185–189, 2004.

Y. G. Zdesenko, et al., “The study of the background of the detectors in the Solotvina underground laboratory,” in Proceedings of the 2nd International Symposium on Underground Physics, vol. 1987, p. 291, Baksan Valley, Russia, 1988.

A. Bettini, “The Canfranc Underground Laboratory (LSC),” The European Physical Journal Plus, vol. 127, article 112, 2012.

H. Ohsumi, H. Ejiri, M. Fujiwara, et al., “The new underground laboratory, "Oto Cosmo Observatory" and ELEGANT Detectors,” in Proceedings of the XV RCNP Osaka International Symposium on Nuclear Physics Frontiers with Electro-Weak Probes, pp. 274–279, World Scientific, March 1996.

E.-I. Esch, T.J. Bowles, A. Hime, A. Pichlmaier, R. Reifarth, and H. Wollnik, “The cosmic ray muon flux at WIPP,” Nuclear Instruments and Methods in Physics Research A, vol. 538, pp. 516–525, 2005.

J. J. Zhu, K. J. Kang, Y. J. Li, et al., “Study on the muon background in the underground laboratory of KIMS,” High Energy Physics and Nuclear Physics, vol. 29, pp. 721–726, 2005.

W. W. M. Allison, et al., “Cosmic ray sun shadow in Soudan 2 underground muon flux,” in Proceedings of the 26th International Cosmic Ray Conference, August 1999.

Y. Fukuda, T. Hayakawa, E. Ichihara, et al., “Measurement of the flux and zenith-angle distribution of upward throughgoing muons by super-kamiokande,” Physical Review Letters, vol. 82, pp. 2644–2648, 1999.

M. Robinson, V. A. Kudryavtsev, R. Lüscher, et al., “Measurements of muon flux at 1070 m vertical depth in the Boulby underground laboratory,” Nuclear Instruments and Methods in Physics Research A, vol. 511, pp. 347–353, 2003.

G. Bellini, G. J. Benzigerf, D. Bick, et al., “Cosmic-muon flux and annual modulation in Borexino at 3800 m water-equivalent depth,” Journal of Cosmology and Astroparticle Physics, vol. 2012, no. 15, article 5, 2012.

Y. M. Andreyev, V. I. Gurentsov, and I. M. Kogai, “Seasonal and diurnal variations of cosmic ray intensity based on measurements at depth of approximately 8.5 * 10⁴ g * cm⁻²,” in Proceedings of the 20th International Cosmic Ray Conference (ICRC '87), vol. 6, p. 200, Moscow, Russia, 1987.

W. Rhode, “Measurements of the muon-flux with the Fréjus-detector,” Nuclear Physics B, vol. 35, pp. 250–253, 1994.

B. Aharmim, S. Ahmed, T. Andersen, et al., “Measurement of the cosmic ray and neutrino-induced muon flux at the Sudbury neutrino observatory,” Physical Review D, vol. 80, Article ID 012001, 15 pages, 2009.

R. Brun and F. Carminati, GEANT Detector Description and Simulation Tool, CERN Programming Library Long Writeup W5013, 1993.

S. Agostinelliae, J. Allison, and K. Amako, “Geant4-a simulation toolkit,” Nuclear Instruments and Methods in Physics Research A, vol. 506, pp. 250–303, 2003.

K. Greisen, et al., “The intensities of the hard and soft components of cosmic rays as functions of altitude and zenith angle,” Physical Review, vol. 61, pp. 212–221, 1942.