Advances in High Energy Physics
Volume 2021 (2021), Article ID 9931708, 5 pages
https://doi.org/10.1155/2021/9931708

Research Article

Determination of Position Resolution for LYSO Scintillation Crystals Using Geant4 Monte Carlo Code

M. F. O. Yahya and F. Kocak

Department of Physics, Bursa Uludag University, 16059 Bursa, Turkey

Correspondence should be addressed to F. Kocak; fkocak@uludag.edu.tr

Received 1 April 2021; Revised 14 July 2021; Accepted 26 August 2021; Published 15 September 2021

Academic Editor: Burak Bilki

Copyright © 2021 M. F. O. Yahya and F. Kocak. 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. The publication of this article was funded by SCOAP³.

Abstract.

LYSO scintillation crystals, due to their significant characteristics such as high light yield, fast decay time, small Moliére radius, and good radiation hardness, are proposed to be used for the electromagnetic calorimeter section of the Turkish Accelerator Center Particle Factory (TAC-PF) detector. In this work, the center of gravity technique was used to determine the impact coordinates of an electron initiating an electromagnetic shower in a LYSO array, in a calorimeter module containing nine crystals, each in cross-section and 200 mm in length. The response of the calorimeter module has been studied with electrons having energies in the range 0.1 GeV-2 GeV. By using the Monte Carlo simulation based on Geant4, the two-dimensional position resolution of the module is obtained as at the center of the crystal.

1. Introduction

Cerium-doped silicate-based heavy scintillation crystals were initially developed for medical applications. Later, the scintillation characteristics of lutetium oxyorthosilicate (LSO) [1] and lutetium-yttrium oxyorthosilicate (LYSO) [2, 3] were found. LYSO crystals have high stopping power (>7 g/cm³), fast decay time (40 ns) and high light yield (200 times of PWO), and superior radiation hardness against gamma rays, neutrons, and protons [4–6]. Furthermore, the crystal emits light in the wavelength region between 360 nm and 600 nm, reaching a peak at 402 nm. Due to all of the above-mentioned properties, this crystal has also attracted the attention of experimental high energy physics research groups working to improve the performance of electromagnetic calorimeters, such as the proposed SuperB forward parts of the endcap calorimeter [7], the KLOE experiment [8], the COMET experiment at J-PARC [9], and the Muon-to-Electron (Mu2e) experiment [10]. In addition to lead tungstate (PWO) and Thallium-activated Cesium Iodide (CsI(Tl)) crystals, LYSO scintillation crystals may also be studied for the electromagnetic calorimeter (ECAL) section of the proposed Turkish Accelerator Center-Particle Factory (TAC-PF) detector [11]. This paper reports the results of a simulation study carried out with Geant4 to evaluate the position resolution of LYSO scintillation crystals for incident electrons at an energy range from 0.1 GeV to 2 GeV.

2. Position Resolution of Electromagnetic Calorimeters

When a particle is sent to the electromagnetic calorimeter, its energy is deposited in the central crystal as well as the crystals around it. The impact position of the particle can be found from the weighted mean of the position of energy deposits in the crystals. The technique used to find the position of the incident particle is called the center of gravity method and calculated by

for the -coordinate and similarly,

for the -coordinate, where is the energy deposited in the crystal and and are the coordinates of the crystal [12]. For the ideal position resolution, it is preferable to take the sum over nine crystals due to the effect of large energy fluctuations in the tails of an electromagnetic cascade [13].

3. Geant4 Simulation and Results

Geant4 simulation code [14] was employed to perform the simulation process for the electrons passing the electromagnetic calorimeter module consisting of LYSO scintillation crystals. The electrons were directed perpendicular to the module in the energy range of 0.1 GeV to 2 GeV. The simulation was performed with Geant4.10.04-patch-03 with the QGSP-BERT4.5 physics list. The LYSO crystals have a length of 200 mm with a cross-section of mm². In the simulation, in order to obtain the distributions of the center of gravity of the deposited energies in the crystals, electrons were injected into the central crystal at fourteen different positions (-12.5 mm, -11mm, -9mm, -7mm, -5mm, -3mm, -1mm, 0, 2mm, 4mm, 6mm, 8mm, 10mm, and 12.5mm), respectively, to scan the entire surface of the crystal. The coordinates of the center of the central crystal were defined to be , . The relation between and in the center of the LYSO matrix for 1 GeV electrons can be seen from Figure 1.

Figure 1: Relation between and in the center of the LYSO matrix for 1 GeV electrons.

The calculated positions () versus the true positions () are shown in Figure 2 for 1 GeV electrons (S-curve). If the impact point is at the center of the crystal or near the boundary between different crystals, the particle’s position can be reconstructed correctly, as can be seen from Figure 2. In other cases, since most of the energy in the shower is deposited in the hit crystal and there is an exponential fall in the energy shared among neighboring crystals, the position of the particle is systematically miscalculated from equation (1).

Figure 2: The position versus the position for 1 GeV electrons. The solid blue line is the S-like curve.

To remove this nonlinear dependence among the true positions and the calculated positions , the S-curve fit function was utilized. This function is an empirical algorithm and given by

Here, and are presented in mm. As a result of this fit, the values of , , and at 1 GeV are found to be 2.669, 0.108, and 0.001, respectively. Table 1 shows the obtained fit results for the incident electron energies between 0.1 and 2 GeV in the direction. As the Moliére radius is loosely dependent on energy, the parameters differ slightly depending on the energy of the incident electron. By employing the values , , and which are obtained from the fit, the corrected position () is calculated by

Table 1: Fit parameters for the incident electrons in the direction.

Electron energy (GeV)
(mm)
(rad/mm)
(mm)

0.1
2.166
0.112
-0.007

0.25
2.327
0.110
0.017

0.5
2.739
0.107
-0.003

0.75
2.702
0.108
-0.004

1
2.669
0.108
0.001

1.25
2.637
0.108
-0.023

1.5
2.651
0.108
-0.003

1.75
2.654
0.108
0.003

2
2.613
0.108
-0.01

Similar calculations have been also made for the direction to obtain the corrected position (). As the corrected position distribution spectra have roughly a Gaussian shape (see Figure 3), they have been fitted using a Gaussian function to obtain the position resolution. The sigmas of these distributions shown in Figure 3 give the calorimeter position resolutions in the and directions for 1 GeV electron energy. Figure 4 shows the corrected position () versus the true position () for 1 GeV electrons in the direction. Figure 5 displays the position resolutions at the center of the central LYSO crystal (at coordinate ) for the incident electron energies from 0.1 GeV to 2 GeV. The Geant4 simulations indicate that the position resolution improves as the energy of the incident electron increases as shown in Figure 5. At the center of the central crystal, the position resolutions depending on the energy of the incident electrons can be parameterized as

for the direction and

for the direction.

Figure 3: and positions after the S-curve correction, for electron energy of 1 GeV.

Figure 4: Dependence of the corrected position in the LYSO crystals on the true coordinate .

Figure 5: Position resolutions in the and directions at the center of the LYSO matrix as a function of the energy. The solid lines represent the fits to the data. The error bars are negligible compared to the symbols shown.

The spatial resolution is not constant on the entire surface of the scintillation crystal and changes depending on the impact position of the incident electron. As can be seen in Figure 6, the simulated position resolution improves towards the edges of the crystal scintillator, and the minimum spatial resolutions are acquired on the edges. This is because the electromagnetic cascade sharing among neighboring crystals starts to become important in that position of the crystal.

Figure 6: Position resolutions versus the impact position of the electron in the and directions.

One can get the two-dimensional position resolution by using and values:

Figure 7 shows the two-dimensional position resolution as a function of energy at the center of the LYSO matrix.

Figure 7: Dependence of the two-dimensional position resolution on electron energy.

In addition, in order to check the accuracy of the simulation code, a prototype electromagnetic calorimeter developed for the COMET experiment given in Ref. [9] was simulated with the Geant4 code. The mentioned prototype consists of LYSO crystals with a dimension of  cm³. The position resolution of the prototype was calculated for the incident electrons using the center of gravity technique mentioned above. The simulation results are well-matched with the experimental data with a slight difference at low energies as can be shown in Figure 8.

Figure 8: Comparison of the position resolutions obtained from Ref. [9] with Geant4 simulation for the crystal matrix.

4. Conclusion

The simulation of the position resolution of the electromagnetic calorimeter module made of LYSO crystals for the proposed TAC-PF detector has been performed using Geant4 simulation code. The position resolution was calculated every 2 mm step for the incident electron energies from 0.1 to 2 GeV, as can be shown in Figure 6. The resolution has improved significantly with S-shape correction. For electrons entering the calorimeter module in the center of the LYSO crystal, the resolution in the () direction is 7.76 mm (7.81 mm) at 0.1 GeV and 2.41 mm (2.32 mm) at 2 GeV. The position resolution was calculated as for the direction, and a similar result was also obtained for the direction. The two-dimensional position resolution was calculated as . Also, by using Function (4) with the fit parameters shown in Table 1 with a minor modification, the impact position of an electron or photon in similar LYSO calorimeters can be deduced from the center of gravity of energy deposits in the crystals.

Data Availability

All figures in the paper represent the results obtained by Monte Carlo simulation, and the data used to support the findings of this study are cited at relevant places within the text as references.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The numerical calculations reported in this paper were partially performed at TUBITAK ULAKBIM, High Performance and Grid Computing Center (TRUBA resources).

References

C. Melcher and J. Schweitzer, “Cerium-doped lutetium oxyorthosilicate: a fast, efficient new scintillator,” IEEE Transactions on Nuclear Science, vol. 39, no. 4, pp. 502–505, 1992.

D. Cooke, K. McClellan, B. Bennett et al., “Crystal growth and optical characterization of cerium-doped Lu1.8Y0.2SiO₅,” Journal of Applied Physics, vol. 88, no. 12, pp. 7360–7362, 2000.

T. Kimble, M. Chou, and B. H. T. Chai, “Scintillation properties of LYSO crystals,” in 2002 IEEE Nuclear Science Symposium Conference Record, pp. 1434–1437, Norfolk, VA, USA, November 2002.

R. Mao, L. Zhang, and R. Y. Zhu, “Gamma ray induced radiation damage in PWO and LSO/LYSO crystals,” in 2009 IEEE Nuclear Science Symposium Conference Record (NSS/MIC), pp. 2045–2049, Orlando, FL, USA, October 2009.

L. Zhang, R. Mao, and R. Y. Zhu, “Effects of neutron irradiations in various crystal samples of large size for future crystal calorimeter,” in 2009 IEEE Nuclear Science Symposium Conference Record (NSS/MIC), pp. 2041–2044, Orlando, FL, USA, October 2009.

G. Dissertori, D. Luckey, F. Nessi-Tedaldi et al., “Results on damage induced by high-energy protons in LYSO calorimeter crystals,” Nuclear Instruments and Methods in Physics Research Section A, vol. 745, pp. 1–6, 2014.

G. Eigen, Z. Zhou, D. Chao et al., “A LYSO calorimeter for the SuperB factory,” Nuclear Instruments and Methods in Physics Research Section A, vol. 718, pp. 107–109, 2013.

M. Cordelli, F. Happacher, M. Martini et al., “CCALT: a crystal calorimeter for the KLOE-2 experiment,” Journal of Physics: Conference Series, vol. 293, no. 1, article 012010, 2011.

K. Oishi, “Development of electromagnetic calorimeter using GSO and LYSO crystals for the J-PARC muon-to-electron conversion search experiment,” in 2014 IEEE Nuclear Science Symposium and Medical Imaging Conference (NSS/MIC), pp. 1–6, Seattle, WA, USA, November 2014.

R. Tschirhart, “The Mu2e experiment at Fermilab,” Nuclear Physics B-Proceedings Supplements, vol. 210-211, pp. 245–248, 2011.

F. Kocak and I. Tapan, “Simulation of LYSO crystal for the TAC-PF electromagnetic calorimeter,” Acta Physica Polonica A, vol. 131, no. 3, pp. 527–529, 2017.

R. Wigmans, Calorimetry: Energy Measurement in Particle Physics, Oxford University Press, 2000.

V. A. Batarin, T. Brennan, J. Butler et al., “Precision measurement of energy and position resolutions of the BTeV electromagnetic calorimeter prototype,” Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, vol. 510, no. 3, pp. 248–261, 2003.

S. Agostinelli, J. Allison, K. Amako et al., “Geant4--a simulation toolkit,” Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, vol. 506, no. 3, pp. 250–303, 2003.