Digitization modeling of a CGEM detector based on Garﬁeld++ simulation

Background The Cylindrical Gas Electron Multiplier Inner Tracker (CGEM-IT) is one of the candidates for the replacement of the current BESIII inner drift chamber, which has developed serious aging effects due to beam-related backgrounds. Purpose To build a full digitization model, it is necessary to understand the electron drift behavior and signal induction. Methods A simulation study based on Garﬁeld++ and ANSYS software has been performed. Results Parameters describing the Lorentz angle, diffusion effect, drift time, multiplication, and induction are obtained from the Garﬁeld++ simulation. Based on these results, a preliminary digitization model is implemented in the BESIII ofﬂine software system.


Introduction
The Main Drift Chamber (MDC), a key sub-detector of the Beijing Spectrometer III (BESIII), has been running for more than ten years at the Beijing Electron Positron Collider II (BEPCII). The inner part of the drift chamber, which is closest to the beam pipe, has developed serious aging effects due to the radiation from beam-related backgrounds. The gain in this region of the drift chamber has decreased year by year, causing a degradation of the hit efficiency and spatial resolution. A Cylindrical Gas Electron Multiplier (CGEM) [1] detector has been an important candidate for the upgrade of the inner tracker due to its good spatial resolution and stability in high-radiation environments. The CGEM inner tracker [2] is designed to be composed of 3 layers of cylindrical GEM detectors. The inner radii of each layer are 76.9 mm, 121.4 mm, and 161.9 mm, and the outer radii are 90.2 mm, 134.7 mm, and 175.2 mm, respectively, as shown in Fig. 1 The length of each layer is 532 mm, 690 mm, and 847 mm, respectively. Each CGEM layer consists of a cathode, 3 GEM foils, and a readout anode, as shown in Fig. 2. The four gaps are called Drift, Transfer1, Trans-fer2, and Induction, with sizes of 5 mm, 2 mm, 2 mm, and 2 mm, respectively. The gas mixture is Ar/iC4H10 with a ratio of 90/10. Electric fields in the drift, transfer, and induction gaps are 1.5 kV/cm, 3 kV/cm, and 5 kV/cm, respectively. The GEM foil has a large number of double-cone shaped holes made by a single mask process using standard photolithography technology. The diameter of each hole on the foil surface is from 50 to 70 µm, and the adjacent holes' pitch is 140 µm, as shown in Fig. 3a. The GEM foil has a three-layer structure, i.e., two layers of copper foil with a standard thickness of 5 µm, and a middle layer of a Kapton foil with a thickness of 50 µm, as shown in Fig. 3b. A high voltage difference of 270 V is applied across each GEM foil. The readout anode is segmented with 650 µm pitch XV patterned strips. The stereo angle between X and V strips changes with the layer. The electrons generated from the ionization by an incident charged particle in the gas drift toward the anode and are multiplied in the holes of the GEM foil. The moving of the  An analog readout is adopted for the CGEM-IT to record both time (T ) and charge (Q) information for each firing anode strip. The position of an incident particle as it enters the drift gap can be estimated from the reconstructed cluster using Q in a charge centroid method or using T in a micro-TPC mode [3].
A Geant4-based [4] simulation of the CGEM-IT has been developed within the BESIII offline software system (BOSS) [5]. The parameters of the incident particle, including the particle type, charge, momentum, and position in the drift gap, can be obtained from the Geant4 simulation. In addition to this information, a digitization process is needed to simulate the outputs of the detector, i.e., a set of firing strips and the T and Q values for each firing strip.
The digitization process contains three steps. The first step is to simulate the ionization, which is implemented by calling the Heed package [6] through the interface provided by the Garfield++ [7] program. The second step is to simulate the drift and multiplication of the ionized and multiplied electrons. Since the complete simulation with Garfield++ of the drift line and the multiplication is time-consuming, this step is simplified by sampling the number, position, and time of the multiplied electrons. The parameters describing the Lorentz angle, multiplication, and diffusion are obtained by a prior simulation of different parts of the CGEM with Garfield++. The third step is to simulate the signal induction which is also based on a prior simulation with Garfield++ and ANSYS.

Methods of the simulation study
The construction of the geometry of the detector and the calculation of the electric field are implemented with ANSYS Coordinate system used in the Garfield++ simulation. The electric field is in the Y direction, and the magnetic field is in the Z direction which is a finite element analysis software [8]. Since the drift and the avalanche of an ionized electron usually occur within a few millimeters, which is very small compared to the size of the whole detector, a planar approximation to the cylindrical shape is adopted to simplify the construction of the detector within ANSYS. The definition of the coordinates is shown in Fig. 4, where the GEM foil is perpendicular to the Y axis and the magnetic field is parallel to the Z axis. Figure 5 shows a unit of the GEM foil and a unit of the anode plane constructed within ANSYS. The field map is solved using ANSYS, which is then used as an input for the simulation with Garfield++. Due to the periodicity of the detector structure, the CGEM detector is realized in Garfield++ by an extension of the unit structure in the X and Z axes.
The drift and multiplication of a single electron are studied with Garfield++. The method is similar to the one in reference [9], but a different parameterization method is adopted in this paper for the needs of the digitization model. As a full simulation of all drift and multiplication processes for a particle passing through the CGEM with Garfield++ is very time-consuming and hard to parameterize, it is divided into three parts. Figure 4 shows the three planes which are used in the Garfield++ simulation. A single electron is simulated as it travels from its generated position to a reference plane after the next GEM foil. This path contains the drift in one gap region and the multiplication in one layer of the foil. The distance between the reference plane and the GEM foil cannot be too small because some of the multiplied electrons may be captured by the electrode on the foil under the influence of the electric field in the region near the foil, as shown in Fig. 11. There are three reference planes, each of which is 150 microns from the GEM foil, as shown in Fig. 4. For the first part, the generated position is the location of ionization in the drift region. The electron is generated uniformly in the X -Z plane with zero initial velocity. In the other two parts, the electron starts from 150 microns after the previous foil and is also generated uniformly in the X -Z plane with zero initial velocity. The Penning coefficient parameter is set to be 0.44 [10] in the simulation. The parameters describing the Lorentz angle, diffusion, gain, and drift time, which will be used in the digitization model, are extracted from the simulation results of a large number of electrons.
After the drift and multiplication in the drift and transfer regions, the electron goes into the induction region. In order to make a model to simulate the induction process, a simulation study of the induced signal on anode strips was performed using ANSYS and Garfield++. The induced signal on an anode strip is described by the Shockley-Ramo theorem: where i(t) is the induced current, q is the particle charge, v is the particle velocity, and E w (r) is the weighting field for a given anode at particle position r. The weighting field is defined as the field produced when a unit potential is applied to the readout electrode and zero potential is applied to all other electrodes [11]. The weighting field of each anode strip was calculated with ANSYS, and then it was imported into Garfield++ with the corresponding standard class to simulate the induced current on the strips. Figure 6 shows the weighting potential on the central X and V strips.

Results of Garfield++ simulation
From the Garfield++ simulation, we can obtain the drift and multiplication behavior of a single electron, including the Lorentz angle, transverse diffusion, drift time, and electron multiplication distribution.

Lorentz angle
The Lorentz force causes a displacement of the electron in the X direction (δ X ), as shown in Fig. 4. To get the Lorentz angle, the electron is generated at different Y . The distance between the electron and the top surface of the GEM foil is δY . The Lorentz angle can be obtained from the line fit to the graph of δ X versus δY . Figure 7 shows the result from the drift region, and Table 1 shows the estimation of the Lorentz angle for each gap. The Lorentz angle in the drift region is bigger because of the lower electric field in this region.

Transverse diffusion
The transverse diffusion has an impact on the cluster size of the multiplied electrons. It can be extracted from fitting the displacement distribution of the multiplied electrons, which is a Gaussian distribution. For the Transfer1 and Trans-fer2 regions, the diffusion effect can simply be described by the width of the Gaussian distribution since the drift distance is almost constant, as shown in Fig. 9b, c. In the drift region, the diffusion effect is a function of δY , as shown in Fig. 8. Figure 9a shows the distribution of multiplied electrons with ionized electrons generated in the middle of the drift region. Table 1 shows the transverse diffusion of each region. The diffusion in the X direction is a little bigger than that in the Z direction because the drift line is not perpendicular to the X direction due to the Lorentz force.

Drift time
As previously mentioned in Section 2, the drift time and its resolution are functions of the drift distance in the drift  Table 1. Also, Table 1 shows the value of the drift time for each region.

Multiplication
Due to the diffusion effect and the impact of the electric field on the drift line, not all the electrons can be effectively multiplied in the GEM foil. Some electrons are captured by the upper electrode before entering the hole, as shown in Fig. 11a. Some electrons are multiplied in the hole, but without multiplied electrons leaving the hole, as shown in Fig. 11b. Only the electron shown in Fig. 11c is multiplied effectively. So the multiplication of the electrons can be described with two parameters, i.e., the transparency and gain. The transparency (τ ) is defined as: where N orign is the number of generated electrons, as shown in Fig. 11d, and N T is the number of generated electrons with at least one multiplied electron reaching the reference plane below the GEM foil. The gain of electrons can be described by the Polya distribution [12]:   where C 0 is a constant, G is the number of multiplied electrons reaching the reference plane, G 0 is the average gain of a single electron, and θ is the parameter determining the variance of the Polya distribution. The value of θ is related to the electronic proportion of energy exceeding the ionization threshold. Figure 12 shows the distribution of G and the fit to the Polya distribution. Figure 13 shows the τ and G 0 parameters as a function of δY in the drift region, from which we see that τ and G 0 are independent of the position of the ionized electron. The τ and G 0 of the drift region are 80.9% and 20.0, respectively. Table 2 shows values of τ and G 0 for each plane. The τ for the first GEM foil is larger than the others, because the electric field in the drift region is smaller than that in the transfer regions [13]. Both the high voltage on the GEM foil and the electric fields in the gap have an impact on G 0 . Different combinations of these two factors cause different values of G 0 in the different regions.  The digitization of multiplication for each GEM foil is divided into two steps. The first step is the sampling for the number of electrons that can be multiplied effectively using the transparency. The second step is the sampling of the gain with a Polya distribution.

Signal induction
The moving of each multiplied electron in the induction region, i.e., from the reference plane below the third GEM foil in Fig. 4 to the anode strip, causes induced currents in several strips. The reference plane and the anode plane with both X and V strips can be divided into periodic areas. The simulation of the induction process is simplified to a small basic area. To make a fast digitization model, the basic area is divided into several grids. The induction process for electrons in each grid is simulated. The strips with induced currents and the shapes of the signals in these strips are recorded, which will be used in the digitization model. In the digitization algorithm of the induction, the final signal in each strip is the sum of the induced currents caused by all multiplied electrons. The current from each multiplied electron is obtained from the prepared record of simulated signal shapes. Figure 14 shows the final induced current on the central X and V strip when a 1.5 GeV vertical incident electron passes through the first layer of the CGEM detector.

Digitization model
Based on the results of the Garfield++ simulation described in the previous section, a preliminary full digitization algorithm has been implemented. The distribution of all the multiplied electrons after passing through the third GEM foil can be obtained from the simulation of the previous processes. Figure 15 shows the cluster size obtained from the digitization algorithm, from which we see that the mean value of the cluster size is 4.1 and 3.4 for X and V , respectively (corresponding to a stereo angle of 43 degrees between the X and V strips). The first version of the digitization algorithm has been successfully integrated into the BOSS framework. The effect of the electronics response [14] and the threshold will also be added. The procedure will be further optimized and tuned once data are collected to make the MC simulation consistent with data.
Currently, it takes a long time to run the digitization algorithm. Simulation of 100 track segments would take 1 min. The main reason is that the sampling of large number of multiplied electrons is very time-consuming. So, how to simplify the sampling process to speed up the algorithm will be considered in the next version.

Conclusion
A simulation study using Garfield++ tool for the CGEM inner tracker has been performed according to the needs of the digitization model. Preliminary results including the Lorentz angle, diffusion, drift time, multiplication, and induction are obtained from the Garfield++ simulation and have been implemented into the digitization algorithm. Comparison with the result from the beam test [15] of the prototype shows that the cluster size from the algorithm is bigger. This difference may be due to the lack of electronic simulation and consideration of threshold effect. What is more, all parameters and settings in this work are the same as the future BESIII CGEM-IT, which are different from the prototype of beam test. Considering those factors, the cluster size from the digitization model is reasonable. The digitization model will be improved in the next step. The samplings in the process of drift and multiplication may be optimized. The difference between the experiment and simulation should be reduced after improvement and tuning with experimental data in the future.