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

Authors: M. F. O. Yahya; F. Kocak
Journal: Advances in High Energy Physics
 (2021)

---

## 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, each25mm×25mm 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 σRmm=3.95±0.08/E⊕1.91±0.11 at the center of the crystal.


---

## Body

## 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/cm3), 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(1)xgravity=∑ixiEi∑iEi,for the x-coordinate and similarly,
(2)ygravity=∑iyiEi∑iEi,for they-coordinate, where Ei is the energy deposited in the ith crystal and xi and yi are the coordinates of the ith 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 3×3 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 17.5X0 with a cross-section of 25×25mm21.2RM. 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 x=0, y=0. The relation between xgravity and ygravity in the center of the LYSO matrix for 1 GeV electrons can be seen from Figure 1.Figure 1
Relation betweenxgravity and ygravity in the center of the LYSO matrix for 1 GeV electrons.The calculated positions (xgravity) versus the true positions (xtrue) 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
Thexgravity position versus the xtrue position for 1 GeV electrons. The solid blue line is the S-like curve.To remove this nonlinear dependence among the true positionsxtrue and the calculated positions xgravity, the S-curve fit function was utilized. This function is an empirical algorithm and given by
(3)xgravity=ctandxtrue−e.Here,xgravity and xtrue are presented in mm. As a result of this fit, the values of c, d, and e 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 X 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 c, d, and e which are obtained from the fit, the corrected X position (xcorrected) is calculated by
(4)xcorrected=1dtan−1xgravityc+e.Table 1
Fit parameters for the incident electrons in theX direction.
Electron energy (GeV)c (mm)d (rad/mm)e (mm)0.12.1660.112-0.0070.252.3270.1100.0170.52.7390.107-0.0030.752.7020.108-0.00412.6690.1080.0011.252.6370.108-0.0231.52.6510.108-0.0031.752.6540.1080.00322.6130.108-0.01Similar calculations have been also made for theY direction to obtain the corrected Y position (ycorrected). 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 X and Y directions for 1 GeV electron energy. Figure 4 shows the corrected position (xcorrected) versus the true position (xtrue) for 1 GeV electrons in the X direction. Figure 5 displays the position resolutions at the center of the central LYSO crystal (at coordinate x=y=0) 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
(5)σxmm=2.77±0.07E⊕1.46±0.10for the X direction and
(6)σymm=2.77±0.05E⊕1.31±0.07for the Y direction.Figure 3
xcorrected and ycorrected positions after the S-curve correction, for electron energy of 1 GeV.Figure 4
Dependence of the corrected positionxcorrected in the LYSO crystals on the true coordinate xtrue.Figure 5
Position resolutions in theX and Y directions at the center of the 3×3 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 Figure6, 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 theX and Y directions.One can get the two-dimensional position resolutionσR by using σx and σy values:
(7)σRmm=σx2+σy2.Figure7 shows the two-dimensional position resolution as a function of energy at the center of the 3×3 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 7×7 LYSO crystals with a dimension of 2×2×12 cm3. 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 7×7 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 Figure6. 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 X (Y) 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 σxmm=2.77±0.07/E⊕1.46±0.10 for the X direction, and a similar result was also obtained for the Y direction. The two-dimensional position resolution was calculated as σRmm=3.95±0.08/E⊕1.91±0.11. 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.



---

*Source: 9931708-2021-09-15.xml*