Simulation study on performance optimization of a prototype scintillation detector for the GRANDProto35 experiment

Simulation study on performance optimization of a prototype scintillation detector for the GRANDProto35 experiment 1 corresp first-author qianxl@sdmu.edu.cn qianxl@sdmu.edu.cn qianxl@sdmu.edu.cn 1 corresp sunhy@sdmu.edu.cn sunhy@sdmu.edu.cn sunhy@sdmu.edu.cn 2 corresp liuc@ihep.ac.cn liuc@ihep.ac.cn liuc@ihep.ac.cn 1 3 School of Intelligent Engineering, Shandong Management University, Jinan 250357, China School of Intelligent Engineering, Shandong Management University, Jinan 250357, China Key Laboratory of Particle Astrophysics, Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 100049, China Key Laboratory of Particle Astrophysics, Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 100049, China Sorbonne Universite, Universite Paris Diderot, Sorbonne Paris Cite, CNRS/IN2P3, LPNHE, Paris Sorbonne Universite, Universite Paris Diderot, Sorbonne Paris Cite, CNRS/IN2P3, LPNHE, Paris As a proposed detector, the giant radio array for neutrino detection (GRAND) is primarily designed to discover and study the origin of ultra-high-energy cosmic rays, with ultra-high-energy neutrinos presenting the main method for detecting ultra-high-energy cosmic rays and their sources. The main principle is to detect radio emissions generated by ultra-high-energy neutrinos interacting with the atmosphere as they travel. GRAND is the largest neutrino detection array to be built in China. GRANDProto35, as the first stage of the GRAND experiment, is a coincidence array composed of radio antennas and a scintillation detector, the latter of which, as a traditional detector, is used to perform cross-validation with radio detection, thus verifying the radio detection efficiency and enabling study of the background exclusion method. This study focused on the implementation of the optimization simulation and experimental testing of the performance of the prototype scintillation detector used in GRANDProto35. A package based on GEANT4 was used to simulate the details of the scintillation detector, including the optical properties of its materials, the height of the light guide box, and position inhomogeneity. The surface of the scintillator and the reflective materials used in the detector were optimized, and the influence of light guide heights and position inhomogeneity on the energy and time resolutions of the detector were studied. According to the simulation study, the number of scintillator photoelectrons increased when changing from the polished surface to the ground surface, with the appropriate design height for the light guide box being 50 cm and the appropriate design area for the scintillator being 0.5 m 2 . The performance of the detector was tested in detail through a coincidence experiment, and the test results showed that the number of photoelectrons collected in the detector was ∼ 84 with a time resolution of ∼ 1 ns, indicating good performance. The simulation results were consistent with those obtained from the tests, which also verified the reliability of the simulation software. These studies provided a full understanding of the performance of the scintillation detector and guidance for the subsequent operation and analysis of the GRANDProto35 experimental array. GRANDProto35 GEANT4 Scintillation detector Light guide height Photoelectrons 1 Introduction 1 Ultra-high-energy cosmic rays (UHECRs), with energies of ≥ e18eV, are the most energetically charged particles ever observed, although their origin remains unclear [ 1 1 ]. They are probably extragalactic in origin and purportedly made in powerful cosmic accelerators; yet, neither is identified. The direct strategy to discover UHECR sources is to detect localized excesses from the arrival direction of UHECRs. However, cosmic-ray charged particles tend to be easily changed in terms of direction because of the effect of galactic and intergalactic magnetic fields at the time of traveling, which prevents us from precisely retracing their trajectories back to their sources. UHECRs interact with cosmic microwave background photons, thus lowering their energy compared to the original energy. Consequently, only a few UHECRs above 40 EeV can reach the Earth from distances beyond 100 Mpc—the Greisen–Zatsepin–Kuzmin cutoff [ 2 2 , 3 3 ]. In contrast, the indirect strategy entails looking for EeV gamma rays and neutrinos produced by UHECRs. UHECRs interact with their accelerating sources and surrounding matter to produce ultra-high-energy (UHE) gamma rays and neutrinos, which can point back at their sources because they are not charged nor affected by cosmic magnetic fields during propagation. However, similar to UHECRs, interactions with the cosmic microwave background cause the UHE gamma rays to lose both energy and their original information. Gamma rays do not reach the Earth beyond 10 Mpc. Instead, they cascade down to GeV–TeV, which makes it difficult to disentangle the gamma rays produced in unrelated phenomena. However, UHE neutrinos barely interact with intergalactic matter during intergalactic propagation because of their small cross section; therefore, they present the best means for the observation of UHECRs. Figure 1 Figure 1 shows the propagation of UHECRs and the associated secondary particles from their sources to the Earth. Because the cosmogenic neutrino flux is not constant and tiny fluxes are involved, an extremely sensitive neutrino detector array is needed for their discovery. The giant radio array for neutrino detection (GRAND), a proposed large-scale neutrino observation array planned in China, was designed for this purpose with the ultimate goal of discovering and studying the sources of UHECRs [ 4 4 - 6 6 4-6 ]. With GRAND, the sources of UHECRs may be better observed by combining cosmic rays, neutrinos, and gamma rays at different energies. Upon arrival at the Earth, UHECRs, gamma rays, and neutrinos initiate large particle showers in the atmosphere. The propagation of charged particles through the geomagnetic field may result in radio emissions. GRAND follows its principle to detect the ground footprint of the radio emission, using an array of 200,000 radio antennas distributed over an area of 200,000 km 2 and operating in the 50–200 MHz band. Because radio detectors are under novel development, and their performance parameters have to be tested by mature traditional detectors, the GRAND Cooperative Group established the GRANDProto35 experiment for testing. The GRANDProto35 prototype coincidence array, as the first construction stage of the GRAND experiment, will lay the foundation for future stages. It is built based on the Tianshan Radio Experiment for Neutrino Detection project [ 7 7 ] located in the Tianshan Mountains in the Xinjiang province of China. It consists of 35 radio detectors and 24 scintillation detectors. A scintillation detector with high sensitivity and detection efficiency for cosmic-ray observation was used to perform cross-validation with radio detection [ 8 8 ]. The goal is to achieve an efficiency of > 80% for the radio detection of showers, with a background rejection that maintains the ratio of false positives to true positives to < 10%. The GRANDProto35 coincidence array is arranged in a rectangle, and the entire array is 800 m long in the east–west axis and 2,400 m long along the north–south axis. Because of the geomagnetic field, radio emission signals generated by cosmic rays perpendicular to the geomagnetic field direction are more concentrated; therefore, such an array arrangement is more conducive to the detection of air showers coming from the north. To improve the efficiency of cosmic-ray detection, scintillation detectors are inclined to the north. The plastic scintillator is a unique detector that is used to detect and measure the properties of energetic charged particles in high-energy accelerators and nuclear and cosmic-ray physics experiments. Large-area plastic scintillators are commonly used as basic detector elements in cosmic-ray-induced extensive air shower experiments to measure the particle densities and their relative arrival times at the observational site. It is ideally suited for this purpose because of its rugged nature, fast response time ( ∼ 1 ns), and reasonable cost. For example, plastic scintillation detectors have been selected for MINOS [ 9 9 ], KASCADE-Grande [ 10 10 ], GRAPES-3 [ 11 11 ], AugerPrime [ 12 12 ], Tibet AS γ [ 13 13 ], and LHAASO [ 14 14 ]. There are several different designs for light collection in scintillation detectors: (1) The photons are collected by direct coupling of a photomultiplier tube (PMT) to one edge or one face of the scintillator [ 15 15 ]. Because the PMT cathode is much smaller than the scintillator area, only a small fraction of the emitted photons is collected. Furthermore, the nonuniform collection of photons is the primary factor influencing the variation in the detector response. (2) The scintillation light is collected using wavelength-shifting fibers, which have been designed for ATHLET [ 16 16 ], AugerPrime, etc. (3) By coupling the tapered plastic (adiabatic acrylic) light guide to the scintillator and PMT to improve detector uniformity, the addition of a plastic light guide to one end of the scintillator results in a tremendous loss of light and an increase in cost. In addition, the Cherenkov photons caused by charged particles in the light guide may cause spurious signals. To collect a large fraction of photon signals uniformly and with a good time resolution, a simple design of photon collection by the PMT mounted at a height below the scintillator is employed in GRANDProto35. The main physical performance requirements for scintillation detectors used in GRANDProto35 are listed in Table 1 Table 1 [ 17 17 ]. The experiment requires that the scintillation detector have high energy and time resolutions. Therefore, in this study, GEANT4 was used to simulate the light emitting, transmitting, and collecting processes of the scintillator. In addition, the simulation also accounted for how the shape of the light guide box and the position inhomogeneity of the detector influence the energy and time resolutions of the detector, thus optimizing the performance of the detector and meeting the requirements for physical indexes of the experiment. MIPs Minimum ionizing particles 1 1 Dynamic range (MIPs[ 1 1 ]) 1 1 1–2000 1 1 Energy resolution (at MIP) 1 1 <30% 1 1 Time resolution (ns) 1 1 <2 1 1 Detection efficiency 1 1 >95% 1 1 Threshold (at MIP) 1 1 0.3 1 1 Detector normal orientation (zenith angle) 1 1 40°–60° The paper is constructed as follows: Sect. 2 describes the GEANT4 simulation process for the detector in detail, Sect. 3 compares the simulation and experimental results, and Sect. 4 details the optimization simulation and test results for the detector, i.e., the height of the light guide box and the position inhomogeneity of detector. Section 5 presents the conclusions of this study. 2 GEANT4 simulation of the detector 1 A plastic scintillation detector uses the fluorescence generated by ionization and excitation of radiation in substances to detect ionizing radiation. The main working process is described as follows: (1) When the charged particles enter the scintillator, they interact with the scintillation crystal, thereby exciting and ionizing the atoms and molecules. (2) the excited atom is then deexcited to emit fluorescent photons with a wavelength in the visible band. (3) The light guide is used to collect as many scintillation photons to the PMT photocathode as possible, after which photoelectrons are emitted from the photocathode. (4) The photoelectrons multiply in the PMT to generate pulse signals. (5) The pulse signals are recorded and analyzed by electronics. The best energy resolution was obtained by collecting only the maximum number of photons. For this reason, very good light collection is necessary; in this case, the parameters affecting light collection are an important criterion. For instance, light collection is largely influenced by different factors, including the type and reflectivity coefficient of reflective materials, the response of reflection and refraction of the photon at the boundary media, the roughness of the scintillator surface, the shape of the light guide, and the area of the scintillator. Therefore, these parameters need to be considered comprehensively in the detector simulation and design. In the past decades, various simulation studies on parameters influencing the light collection of scintillation detectors have been conducted [ 18 18 , 19 19 ]. In terms of the studies on the surface reflectance properties, Janecek and Moses [ 20 20 ] proposed the use of measured reflectance data based on the experimental characterization of the angular reflection distribution of scintillators. Roncali and Cherry [ 21 21 ] proposed a different approach based on the three-dimensional measurements of crystal surfaces. Surface samples were scanned using three-dimensional atomic force microscopy and used to compute the reflectance properties of the crystal surfaces. These studies provided the reference methods and models for the design of the scintillation detector. Optical Monte Carlo (MC) simulations are frequently used to predict the light distribution in scintillation crystals, together with light collection by photodetectors. Current MC simulation software such as DETECT [ 22 22 ], Litrani [ 23 23 ], GEANT4 [ 24 24 , 25 25 ], and GATE [ 26 26 ] have been developed to perform such tasks. GEANT4 is a toolkit for simulating the passage of particles through matter. It can simulate the processes of particle ionizing in the scintillator and the transmission of scintillation light [ 27 27 , 28 28 ]. All physical processes are simulated in detail based on the time when charged particles enter the detector at the time of detecting scintillation photons in photosensitive areas. This section introduces the setting of parameters in the GEANT4 simulation and the entire simulation process. 2.1 Geometric structural design of the detector 2 In this design, the scintillator is placed at the top of an inverted trapezoidal shaped, light-tight container with reflective internal surfaces. The reflective surface guides the scintillation photons to the PMT after multiple diffuse reflections. Therefore, the detector mainly consists of a plastic scintillator, reflective coating, light guide device, light guide box (air light guide), and light collecting device (PMT). A stainless-steel housing was adopted for the detector. The dimensions of the optimized plastic scintillator were 707 × 707 × 20 mm 3 (for a total area of 0.5 m 2 ), and the upper and lower dimensions of the light guide box housing were 730 × 730 mm 2 and 150 x150 mm, respectively, with a vertical height of 500 mm. In view of the requirements for a large dynamic range of GRANDProto35, a dynode readout design was adopted for the PMT, for which the Hamamatsu R7725 model was selected. Because this experiment mainly focuses on cosmic rays with large zenith angles, an angle-adjustable baseplate was designed below the detector. Figure 2 Figure 2 shows a schematic of the scintillation detector. 2.2 Optical interface design for the scintillator 2 When the scintillation light produced in the scintillator encounters the edges of the scintillator, photons may be reflected or transmitted, and some of them may escape from the crystal. To strengthen the collection of light, Tyvek, a reflective material, was affixed to the upper and side surfaces of the scintillator, and the scintillator and the reflective material were tightly fitted by air. According to the refractive indexes of the scintillator and air, light at an emission angle within the range of the critical angle ( θ c ) can escape from the scintillator and enter the light guide to reach the PMT. The remaining light escapes from the scintillator after some total internal reflections at the ends of the scintillator or is lost during the total internal reflections between the upper and lower interfaces of the scintillator, as shown in Fig. 3a Fig. 3a . θ c , the critical angle (39.2°) is calculated as follows: <math display="block" id="M1"><mrow><msub><mi>θ</mi><mtext>c</mtext></msub><mo>=</mo><msup><mstyle mathsize="140%" displaystyle="true"><mrow><mi>sin</mi></mrow></mstyle><mrow><mo>−</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><mrow><mfrac><mrow><msub><mi>n</mi><mrow><mtext>a</mtext><mi>i</mi><mi>r</mi></mrow></msub></mrow><mrow><msub><mi>n</mi><mrow><mtext>s</mtext><mi>c</mi><mi>i</mi><mi>n</mi><mi>t</mi></mrow></msub></mrow></mfrac></mrow><mo>)</mo></mrow><mn>,</mn></mrow></math> where n air and n scint refer to the refractive indexes of the air and scintillator, respectively. To reduce the number of total internal reflections in the scintillator and increase the output of the light, the lower surface of the scintillator was ground, so that the photons to be collected may hit the ground surface at a large incident angle, thus escaping from the scintillator more easily. The process is illustrated in Fig. 3b Fig. 3b . In practical applications, the upper surface is polished, the lower surface is ground, and the side surface is cut for the design of the scintillator. 2.3 Scintillator emission spectrum, PMT quantum efficiency, and Tyvek reflectivity 2 An EJ-200 scintillator (Eljen Technology, USA) was used as the plastic scintillator. It is characterized by a long optical attenuation length and fast time response, with a light yield of 10,000 photons/MeV. The light emission spectrum of the plastic scintillator and the quantum efficiency of the PMT photocathode are shown in Fig. 4 Fig. 4 . In the figure, the black line indicates the EJ-200 emission spectrum, which is proportional to the probability of light emission at a given wavelength when the particles deposit energy, while the red line shows the quantum efficiency of the PMT at the corresponding wavelength. Within the scope of the emission spectrum of the scintillator, on average, the quantum efficiency of the PMT photocathode was ∼ 0.21. The performance parameters of the scintillator and PMT used in the simulation are listed in Table 2 Table 2 . 2 1 EJ-200 material properties 1 1 Light output (% anthracene) 1 1 64 1 1 Scintillation efficiency (photons/MeV) 1 1 10,000 1 1 Wavelength of maximum emission (nm) 1 1 425 1 1 Light attenuation length (cm) 1 1 380 1 1 Density (g/cm 3 ) 1 1 1.023 1 1 Refractive index 1 1 1.58 2 1 PMT (R7725) 1 1 Effective photocathode diameter (mm) 1 1 51 1 1 Transit time spread (ns) 1 1 1.85 1 1 Quantum efficiency 1 1 ∼ 0.21 1 1 Electron collection factor 1 1 0.6 1 1 Rise time (ns) 1 1 2.5 1 1 Fall time (ns) 1 1 3.4 Tyvek, as the reflective material, was also affixed to the inner wall of the light guide box increase the detector’s light collection efficiency. Figure 5 Figure 5 shows the reflectivity of different Tyvek models in air measured by the National Institute of Metrology, China. Because Tyvek 1082D clearly has the highest reflectivity among the three models, it was chosen. Its performance parameters are listed in Table 3 Table 3 . When the wavelength exceeded 380 nm, the reflectivity was ∼ 95%, and the reflectivity was high within the sensitive wavelength range of the scintillator and PMT. All of these performance parameters were digitized and added to the GEANT4 code to perform the simulation properly. 1 1 Unit mass (g/m 2 ) 1 1 105 1 1 Thickness ( μ m) 1 1 275 1 1 Opacity 1 1 98.4% 1 1 Longitudinal elongation (N/(25.4 mm) 1 1 290 1 1 Transverse elongation (N/(25.4 mm)) 1 1 330 1 1 Bursting strength (kPa) 1 1 1,700 2.4 Determination of optical simulation parameters 2 In GEANT4, there are two optical reflection models to be selected by the user: GLISUR and UNIFIED. These are used to describe the optical properties of medium surfaces. The GLISUR model, originating from GEANT3, defines the interface between two media as polished or ground, with only an adjustable parameter describing the roughness of the interface. The UNIFIED model is available for an accurate description of the optical process of the ground surface of various surface shapes over a long wavelength range. In the UNIFIED model, four types of surface reflections are possible: specular spike, specular lobe, backscatter, and Lambertian, as shown in Fig. 6 Fig. 6 . For specular reflection on an average surface, the reflected photons are reflected at the average surface normal. For backscatter reflection, the photons are reflected back in the direction from which they came. In Lambertian reflection, the photons are reflected with a Lambertian probability distribution, that is, into a cosine distribution around the average surface normal [ 20 20 , 31 31 ]. It is assumed in the optical UNIFIED model that the surface of a medium is composed of an infinite number of “micro-facets". The entire surface is defined as an average surface, the micro-facets have their own normal directions, and the directions of all the micro-facets are averaged to obtain a direction. The included angle of the normal direction for the micro-facet and the normal direction for the average plane ( α ) follows a Gaussian distribution. The standard deviation of the distribution ( σ α ) is used to describe the roughness of the surface. Each time a specular lobe interaction occurs, a micro-facet is randomly selected from this distribution, and a specular reflection is then calculated through the summation of all of the micro-facet orientations, as shown in Fig. 7 Fig. 7 . When a beam of photons with a certain momentum is incident on the surface of the medium, UNIFIED uses four parameters, namely, C sl , the specular reflection probability of the micro-facet; C ss , the specular reflection probability of the average surface; C bs , the backscattering probability; and C dl , the Lambertian reflection probability, to describe the radiant intensity of the surface. According to the model, C sl = 1, and σ α = 0 corresponds to the case when the surface is completely polished. Figure 8 Figure 8 shows 240× magnified images of the polished surface, the ground surface, and the cut surface of EJ-200 scintillator measured by using a VMS-1510F Video Measuring Instrument. Obviously, the three surfaces are different in roughness, and the polished surface of the scintillator is not completely polished; in that case, roughness was added in the simulation to a certain extent. The probabilities of specular spike reflection and backscatter reflection for these surfaces were very small; therefore, they were set to 0. Tyvek reflective material is made of high-density polyethylene fiber, which is opaque and not smooth. Therefore, both specular and diffuse reflection can take place on the Tyvek surface, and diffuse reflection is the main component [ 32 32 , 33 33 ]. In addition, the probability of backward scattering is very small ( C bs = 0). Tyvek exhibited a micro-faceted structure; therefore, the probability of specular spike reflection off it is also very small ( C ss = 0). The optical parameters of the scintillator surfaces and Tyvek surface used in the simulation are listed in Table 4 Table 4 . 1 1   1 1 C sl 1 1 C ss 1 1 C dl 1 1 C bs 1 1 σ α 1 1 Polished surface 1 1 98.9% 1 1 0 1 1 1.1% 1 1 0 1 1 0.04 1 1 Ground surface 1 1 89% 1 1 0 1 1 11% 1 1 0 1 1 0.1 1 1 Cut surface 1 1 96.7% 1 1 0 1 1 3.3% 1 1 0 1 1 0.12 1 1 Tyvek surface 1 1 20% 1 1 0 1 1 80% 1 1 0 1 1 0.2 2.5 Simulation of the PMT output signal 2 The number of photons emitted from the scintillator is proportional to the energy loss of the charged particles at a rate assumed to be one photon per 100 eV [ 34 34 ]. The angular distribution of the emitted light is isotropic. The time profile of scintillation emission is described by the following emission time ( t emit ) probability function [ 35 35 ]: <math display="block" id="M2"><mrow><mi>n</mi><mo stretchy="false">(</mo><msub><mi>t</mi><mrow><mtext>e</mtext><mi>m</mi><mi>i</mi><mi>t</mi></mrow></msub><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><msub><mi>n</mi><mtext>f</mtext></msub></mrow><mrow><msub><mi>τ</mi><mtext>f</mtext></msub></mrow></mfrac><msup><mi>e</mi><mrow><mo>−</mo><mfrac><mrow><msub><mi>t</mi><mrow><mtext>e</mtext><mi>m</mi><mi>i</mi><mi>t</mi></mrow></msub></mrow><mrow><msub><mi>τ</mi><mtext>f</mtext></msub></mrow></mfrac></mrow></msup><mo>+</mo><mfrac><mrow><msub><mi>n</mi><mtext>s</mtext></msub></mrow><mrow><msub><mi>τ</mi><mtext>s</mtext></msub></mrow></mfrac><msup><mi>e</mi><mrow><mo>−</mo><mfrac><mrow><msub><mi>t</mi><mrow><mtext>e</mtext><mi>m</mi><mi>i</mi><mi>t</mi></mrow></msub></mrow><mrow><msub><mi>τ</mi><mtext>s</mtext></msub></mrow></mfrac></mrow></msup><mn>,</mn></mrow></math> where τ f refers to the fast decay time constants, τ s represents the decay time of the slow component, and n f and n s denote the numbers of fast and slow components, respectively. Because the signal timing is determined at the leading edge in the discriminator, τ f is the most important parameter for the simulation. In this simulation, τ f is 2.1 ns and τ s is 14.2 ns, where the emission time probability function can produce the same rise time as the scintillator light output. The PMT produces a pulse signal if a single photoelectron is generated. The pulse waveform of a single photoelectron was simulated using the following time response function [ 36 36 ]: <math display="block" id="M3"><mrow><msub><mi>v</mi><mtext>i</mtext></msub><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mo>=</mo><mi>G</mi><msub><mi>C</mi><mtext>e</mtext></msub><mfrac><mrow><msup><mi>t</mi><mn>2</mn></msup><msup><mi>e</mi><mrow><mo>−</mo><msup><mi>t</mi><mn>2</mn></msup><mo>/</mo><msup><mi>τ</mi><mn>2</mn></msup></mrow></msup></mrow><mrow><mstyle displaystyle="true"><mrow><mo>∫</mo></mrow></mstyle><msup><mi>t</mi><mn>2</mn></msup><msup><mi>e</mi><mrow><mo>−</mo><msup><mi>t</mi><mn>2</mn></msup><mo>/</mo><msup><mi>τ</mi><mn>2</mn></msup></mrow></msup><mtext>d</mtext><mi>t</mi></mrow></mfrac><mn>,</mn></mrow></math> where G is the gain of the PMT, C e represents the charge-to-voltage conversion factor, and τ denotes the time constant of 2.5 ns, that is, the rise time of the PMT R7725 used for the detector. When a charged particle deposits its energy in the scintillator, optical photons are emitted, and multiple photoelectrons can be produced from the PMT. Because each photoelectron pulse arrives at the PMT at a different time, the output signal of the PMT can be obtained through the summation of the pulse waveforms of individual photoelectrons, which can be calculated using the following formula: <math display="block" id="M4"><mrow><msub><mi>V</mi><mrow><mtext>P</mtext><mi>M</mi><mi>T</mi></mrow></msub><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mo>=</mo><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mtext>i</mtext><mo>=</mo><mn>1</mn></mrow><mrow><msub><mi>n</mi><mrow><mtext>p</mtext><mi>e</mi></mrow></msub></mrow></munderover></mstyle><msub><mi>v</mi><mtext>i</mtext></msub><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mn>,</mn></mrow></math> where n pe represents the number of photoelectrons. Figure 9a, b Figure 9a, b show a typical single photoelectron pulse waveform and the total PMT output signal waveform, respectively, in the simulation. The red dots in Fig. 9b Fig. 9b represent the waveform of the simulated output pulse signal. To make a comparison with the simulated output pulse, the experimental output pulse waveform of the optimized prototype detector was added. It can be clearly seen that the simulation data were consistent with the output signal obtained in the experiment, which proved the validity of the simulation for signal amplitude and time information. In the actual operation of the detector, the PMT output analog signal was amplified and converted into a digital signal by the front-end electronics, after which the charge and time information of the signal were obtained. Given the influence of noise, such as electronic noise, threshold discrimination was applied to the pulse in the GRANDProto35 electronics. In the simulation, the threshold voltage was set to -50 mV, corresponding to ∼ 0.3 MIPs. In fact, there was a time-walk effect when measuring time in this way. The effect was corrected in both the experiment and the simulation. 3 Comparison between simulation and experiment 1 The photoelectron energy spectrum and arrival time distribution of the detector can be obtained based on the GEANT4 simulation of the material composition, parameter setting, and the pulse signal output of the detector. To verify the final simulation results, a prototype detector was constructed using the geometric structure described in Sect. 2.1. (The height and area optimization of the detector are described in Sect. 4.) The details were then tested and analyzed. To eliminate the interference of noise and accurately measure the charge and time information of different positions on the detector surface, a 5 × 5 × 5 cm 3 scintillator was placed above the prototype detector for the cosmic-ray coincidence event test. The testing principle is shown in Fig 10 Fig 10 . When passing through the 1/2 divider, the output signals of the coincidence scintillator and prototype scintillator are divided into two channels, one of which enters the analog-to-digital converters (CAEN V965) for charge measurement and the other is discriminated by a constant fraction discriminator (CAEN N843) and used as the start and stop signals of the time-to-digital converter (CAEN V775N) for time measurement. To quantify the contribution of the Tyvek reflective materials to the increase in the light output, a comparative experiment was performed. The Tyvek material on the inner walls of the scintillator and light guide were changed to black paper. Four different packaging methods were tested for simulation: Scin tyvek - LG tyvek , Scin tyvek - LG black , Scin black - LG tyvek , and Scin black - LG black . Figure 11 Figure 11 illustrates the comparison of the photoelectron spectra between the simulation and experimental data collected by the PMT under Scin tyvek - LG tyvek packaging. It is clear that, when the largest number of photoelectrons, ∼ 84 photoelectrons, are collected by the PMT under this packaging, the spectra of the simulation and experimental data are in good agreement. The energy resolution of the detector is defined as <math display="block" id="M5"><mrow><mi>σ</mi><mo>=</mo><mfrac><mi>τ</mi><mrow><mn>2.355</mn></mrow></mfrac><mn>,</mn></mrow></math> where τ refers to the ratio of the full width at half maximum (Δ E ) to the most probable value ( E ) of the photoelectron distribution, from which the energy resolution is calculated to be ∼ 25%. Table 5 Table 5 shows the test results of all four packaging types. The experiments and MC simulations indicate that the design of Tyvek packaging around the scintillator and on the inside wall of the light guide box significantly improved the collection efficiency of the detector. Consequently, Tyvek packaging was used for subsequent simulations. To compare the difference of light output between the ground and the polished surfaces of the scintillator, a comparative test was also performed. The test results showed that the number of photoelectrons increased by ∼ 20% when the lower surface was changed to ground. Figure 12 Figure 12 shows the time difference distribution information detected by the prototype detector and the MC simulation data. It can be clearly seen that the experimental and simulation results are also consistent, and the time resolution is ∼ 1 ns. 2 1 Packaging type 1 2 N pe (data) 1 2 N pe (MC) 1 2 <math id="M6"><mrow><mfrac><mrow><mtext>M</mtext><mi>C</mi><mo>−</mo><mi>d</mi><mi>a</mi><mi>t</mi><mi>a</mi></mrow><mrow><mtext>d</mtext><mi>a</mi><mi>t</mi><mi>a</mi></mrow></mfrac></mrow></math> 1 1 Scintillator 1 1 Light guide 1 1 Tyvek 1 1 Tyvek 1 1 84.0 1 1 86.1 1 1 2.5% 1 1 Tyvek 1 1 Black paper 1 1 8.5 1 1 8.3 1 1 -2.3% 1 1 Black paper 1 1 Tyvek 1 1 11.5 1 1 12.9 1 1 12.1% 1 1 Black paper 1 1 Black paper 1 1 3.4 1 1 4.1 1 1 20.6% 4 Simulation of detector performance optimization 1 4.1 Influence of light guide heights on energy and time resolutions of the detector 2 The energy and time resolutions of the detector are closely related to the structure of the light guide box. Therefore, we simulated the optimization for the light guide heights under the condition of a fixed scintillator area (0.5 m 2 ). Muons are the dominant component of charged particles at sea level. Muon samples with energies of > 500 MeV were simulated and distributed over the entire surface of the detector. Figure 13 Figure 13 shows the distributions of photoelectrons collected by the PMT for different light guide heights ( LG H ), that is, 10, 30, 50, 70, and 90 cm. It can be seen from the figure that the number of photoelectrons gradually decreased with the increase of the light guide height, with the most probable value decreasing from 134.7 to 56.6, and the energy resolution gradually becoming better to a certain extent. Table 6 Table 6 lists the specific values of the number of photoelectrons and the energy resolutions at different light guide heights. Owing to the short propagation distance, more photons can be collected by the PMT for a low light guide height; however, the photoelectron distribution can be widened at the same time. This is largely due to the obvious position inhomogeneity at a low light guide height. For a high light guide height, the opposite is true. 1 1 LG H (cm) 1 1 10 1 1 30 1 1 50 1 1 70 1 1 90 1 1 N pe 1 1 134.7 1 1 112.8 1 1 84.3 1 1 70.7 1 1 56.6 1 1 Energy resolution (%) 1 1 30.7 1 1 27.2 1 1 26.5 1 1 26.3 1 1 25.1 1 1 Time resolution (ns) 1 1 0.86 1 1 0.91 1 1 0.94 1 1 1.12 1 1 1.28 Figure 14a Figure 14a shows the time resolution of the photons generated by the scintillator arriving at the PMT at different light guide heights, and Fig. 14b Fig. 14b shows the total time distribution of all the photons generated by the scintillator after transmission through the light guide box. As can be seen from the two figures, the lower the light guide height, the less the time it took for the photons to reach PMT, and the better the time resolution. The reason for this is that, when the PMT was placed close to the scintillator, the photons had less distance to travel in the light guide box, and more photons arrived at approximately the same time; therefore, the time resolution was better, which was consistent with our expectation. The specific values of the time resolution are given in Table 6 Table 6 . Based on the above results, we concluded that the lower the light guide height, the greater the number of collected photoelectrons. However, the energy resolution was poor, whereas the time performance was relatively good. Based on this, we selected the middle height (50 cm) as the light guide height. However, we still need to further consider the position response of the scintillator to determine the light guide height. 4.2 Influence of position inhomogeneity on energy and time resolutions of the detector 2 In the case in which the cosmic-ray particles hit different positions on the scintillator, the number of photoelectrons and time response will be different, which will affect the energy and time resolutions of the detector. Therefore, specific simulation and experimental studies were conducted on the position homogeneity of the detector. 4.2.1 Simulation of position inhomogeneity 3 In the simulation, plumbing of different positions on the scintillator was used, different positions on the diagonal of the scintillator were selected, and a total of 18 points were tested. At each position, the photoelectron spectra at different light guide heights were simulated. Each photoelectron spectrum was fitted, and Landau and exponential convolution functions were adopted to obtain the most probable number of photoelectrons. The number of photoelectrons collected at the central point of the scintillator was regarded as 1, and the ratio of photoelectrons at other positions is defined as follows: <math display="block" id="M7"><mrow><mtext>r</mtext><mi>a</mi><mi>t</mi><mi>i</mi><mi>o</mi><mo>=</mo><mfrac><mrow><msub><mi>n</mi><mrow><mtext>p</mtext><mi>o</mi><mi>s</mi><mi>i</mi><mi>t</mi><mi>i</mi><mi>o</mi><mi>n</mi></mrow></msub></mrow><mrow><msub><mi>n</mi><mrow><mtext>c</mtext><mi>e</mi><mi>n</mi><mi>t</mi><mi>e</mi><mi>r</mi></mrow></msub></mrow></mfrac><mn>.</mn></mrow></math> According to the simulation results shown in Fig. 15 Fig. 15 , it is clear that, the closer to the edge of the scintillator one gets, the fewer the collected photoelectrons. This is because the photons generated at the edge of the scintillator were more easily lost during the propagation in the light guide box, and they cannot be received by the PMT. However, the higher the light guide height, the better the position homogeneity of the scintillator. When the height was 90 cm, the difference between the number of photoelectrons at the edge and that at the center of the scintillator was ∼ 15%. Furthermore, the position inhomogeneity changed more significantly as the light guide height decreased. These simulation results are consistent with expectations. Based on the results of the number of photoelectrons and energy and time resolutions of the scintillator under different light guide heights simulated in Sect. 4.1, a 50 cm light guide height was selected as optimal. 4.2.2 Influence of position inhomogeneity on the energy resolution of the detector 3 To verify the simulation, the coincidence test method used in Sect. 3 was adopted to test the position inhomogeneity of the scintillation detector. Points A – G , that is, a total of seven position points on the scintillator, were tested, as shown in Fig. 16a Fig. 16a , and the light guide height was set to the optimum value (50 cm). Figure 16b Figure 16b shows the results of the comparison between the measured experimental data and the simulation. When the distance was < 35 cm (the “ F " position) from the central position of the scintillator, the number of photoelectrons decreases slowly; in contrast, when the distance is > 35 cm, the number of photoelectrons decreases more quickly and the edge effect is more significant. The maximum inhomogeneity of the scintillation detector is ∼ 15%. Table 7 Table 7 lists the number of photoelectrons, energy resolution, and position inhomogeneity results for all test points measured in both the experiment and simulation. 1 1 Performance 7 1 Position 1 1   1 1 A 1 1 B 1 1 C 1 1 D 1 1 E 1 1 F 1 1 G 1 1 Distance from center (cm) 1 1 0 1 1 7.07 1 1 14.14 1 1 21.21 1 1 28.28 1 1 35.35 1 1 42.42 1 1 N pe (data) 1 1 84.0 1 1 83.58 1 1 82.94 1 1 81.96 1 1 80.40 1 1 78.39 1 1 73.75 1 1 N pe (MC) 1 1 86.15 1 1 85.52 1 1 85.11 1 1 84.05 1 1 82.41 1 1 80.0 1 1 74.45 1 1 Homogeneity (data) (%) 1 1 100 1 1 99.50 1 1 98.75 1 1 97.58 1 1 95.73 1 1 93.33 1 1 87.80 1 1 Homogeneity (MC) (%) 1 1 100 1 1 99.28 1 1 98.80 1 1 97.56 1 1 95.67 1 1 92.87 1 1 86.43 1 1 Energy resolution (data) (%) 1 1 27.38 1 1 26.42 1 1 27.38 1 1 26.00 1 1 29.42 1 1 26.58 1 1 29.06 1 1 Energy resolution (MC) (%) 1 1 26.88 1 1 24.97 1 1 26.87 1 1 26.31 1 1 26.02 1 1 26.35 1 1 26.18 4.2.3 Influence of position inhomogeneity on the time resolution of the detector 3 Before a pulse signal enters the time-to-digital converter, it is discriminated by the discriminator, and the time-walk effect may appear [ 37 37 ]. When the leading edge of the pulse signal drops rapidly, the relation between the amplitude of the output signal and time can be expressed as follows: <math display="block" id="M8"><mrow><mi>V</mi><mo stretchy="false">(</mo><mi>T</mi><mo stretchy="false">)</mo><mo>=</mo><msub><mi>V</mi><mtext>0</mtext></msub><msup><mrow><mrow><mo>(</mo><mrow><mfrac><mrow><mi>T</mi><mo>−</mo><msub><mi>T</mi><mtext>0</mtext></msub></mrow><mrow><msub><mi>T</mi><mtext>R</mtext></msub></mrow></mfrac></mrow><mo>)</mo></mrow></mrow><mn>2</mn></msup><mn>,</mn></mrow></math> where V 0 refers to the maximum value of the pulse signal, T R represents the rise time of the signal, and T 0 denotes the start time of the signal. The time at the discrimination threshold, T th , can be described as follows: <math display="block" id="M9"><mrow><msub><mi>T</mi><mrow><mtext>t</mtext><mi>h</mi></mrow></msub><mo>=</mo><msub><mi>T</mi><mtext>R</mtext></msub><msqrt><mrow><mfrac><mrow><msub><mi>V</mi><mrow><mtext>t</mtext><mi>h</mi></mrow></msub></mrow><mrow><msub><mi>V</mi><mtext>0</mtext></msub></mrow></mfrac></mrow></msqrt><mo>+</mo><msub><mi>T</mi><mtext>0</mtext></msub><mn>.</mn></mrow></math> Because the charge of the pulse signal is proportional to the amplitude, that is, Q ∝ V 0 , the relation between T th and Q can be expressed as follows: <math display="block" id="M10"><mrow><msub><mi>T</mi><mrow><mtext>t</mtext><mi>h</mi></mrow></msub><mo>=</mo><mfrac><mi>k</mi><mrow><msqrt><mi>Q</mi></msqrt></mrow></mfrac><mo>+</mo><msub><mi>T</mi><mtext>0</mtext></msub><mn>,</mn></mrow></math> where k is the time-walk factor. To accurately obtain the real-time resolution of the prototype scintillation detector, the charge and time were corrected, and the results showed that the time resolution of the detector can be improved by ∼ 20% through charge and time correction. Figure 17 Figure 17 shows the time resolution at different positions of the detector measured in the experiment. With an increase in the distance from the center of the detector, the time resolution at each point gradually deteriorated. However, the time resolutions were relatively good overall. The overall average time resolution is ∼ 1 ns, and the time-walk factor k obtained by fitting Eq. 9 is ∼ 78 ns·pC 1/2 . The experiment and simulation of position inhomogeneity showed that, with the increase in distance from the center, the number of photoelectrons gradually decreased and the time resolution also gradually deteriorated. Given the number of photoelectrons, energy resolution, and time resolution, the optimal area of the scintillation detector is 0.5 m 2 (i.e., the length of the side is 70.7 cm). In addition, the influence of scintillator thickness was also tested with different values, and the investigation showed that a scintillator with 2 cm thickness could fulfill the requirements of the experiment. 5 Conclusion 1 As a coincidence array, the GRANDProto35 experimental array consists of radio antennas and a scintillation detector. The scintillation detector is mainly used for cross-validation with radio observations to improve the efficiency of radio detection for cosmic rays. The performance of the scintillation detector has an impact on the ability to detect UHECRs; therefore, it is necessary to conduct simulations and experiments to optimize the performance of the scintillation detector. A specific GEANT4 simulation was conducted in which the optical properties of the materials used for the detector, height of the light guide, and position inhomogeneity of the detector were simulated. Based on the simulation study, to increase the number of photoelectrons collected in the detector, a ground surface of the scintillator was designed, and a Tyvek reflective material was added around the scintillator and on the inside wall of the light guide box. The light guide height and position inhomogeneity of the detector were studied in terms of the energy and time resolutions of the detector. According to the results, a light guide height of 50 cm and an area of 0.5 m 2 were optimal for the structural design of the detector. A series of verification experiments were performed, and the simulation results were found to be consistent with the experimental results, which verified the validity of our simulation. These results are very important for subsequent GRANDProto35 experiments, as well as for analyses of relevant physical targets. References K. Kotera , A.V. Olinto , The astrophysics of ultrahigh-energy cosmic rays . Annu. Rev. Astron. Astr . 49 , 119 ( 2011 ). doi: 10.1146/annurev-astro-081710-102620 http://doi.org/10.1146/annurev-astro-081710-102620 The astrophysics of ultrahigh-energy cosmic rays K. Greisen , End to the cosmic-ray spectrum? Phys. Rev. Lett . 16 , 748 ( 1966 ). doi: 10.1103/PhysRevLett.16.748 http://doi.org/10.1103/PhysRevLett.16.748 End to the cosmic-ray spectrum? G.T. Zatsepin , V.A. Kuz'min , Upper limit of the spectrum of cosmic rays . Jetp L . 4 , 78 ( 1966 ). Upper limit of the spectrum of cosmic rays J. Álvarez-Muñiz , R.A. Batista , J. Bolmont et al ., The giant radio array for neutrino detection (GRAND): science and design . Sci. China. Phys. Chem . 63 , 219501 ( 2020 ). doi: 10.1007/s11433-018-9385-7 http://doi.org/10.1007/s11433-018-9385-7 The giant radio array for neutrino detection (GRAND): science and design O. Martineau-Huynh , The Giant Radio Array for Neutrino Detection . Paper Presented at the 36th International Cosmic Ray Conference ( Madison, Wisconsin, USA , July 24 -August 1. 2019 ) The Giant Radio Array for Neutrino Detection M. Tueros , Grand, a giant radio array for neutrino detection: objectives, design and current status . EPJ Web. Conf . 216 , 01006 ( 2019 ). doi: 10.1051/epjconf/201921601006 http://doi.org/10.1051/epjconf/201921601006 Grand, a giant radio array for neutrino detection: objectives, design and current status D. Charrier , K.D. de Vries , Q.B. Gou et al ., Autonomous radio detection of air showers with the TREND50 antenna array . Astropart. Phys . 110 , 15 ( 2019 ). doi: 10.1016/j.astropartphys.2019.03.002 http://doi.org/10.1016/j.astropartphys.2019.03.002 Autonomous radio detection of air showers with the TREND50 antenna array Q.B. Gou , O. Martineau-Huynh , J. David et al ., The GRANDproto35 experiment . Paper Presented at the 35th International Cosmic Ray Conference ( Bexco, Busan, Korea 10-20-July. 2017 ) The GRANDproto35 experiment D.G. Michael , P. Adamson , T. Alexopoulos et al ., The magnetized steel and scintillator calorimeters of the MINOS experiment . Nucl. Instrum. Meth. A . 596 , 190 - 228 ( 2008 ). doi: 10.1016/j.nima.2008.08.003 http://doi.org/10.1016/j.nima.2008.08.003 The magnetized steel and scintillator calorimeters of the MINOS experiment E.T. Antoni , W. Apel , F. Badea et al ., The cosmic-ray experiment KASCADE . Nucl. Instrum. Meth. A . 513 , 490 ( 2003 ). doi: 10.1016/S0168-9002(03)02076-X http://doi.org/10.1016/S0168-9002(03)02076-X The cosmic-ray experiment KASCADE S.K. Gupta , Y. Aikawa , N.V. Gopalakrishnan et al ., GRAPES-3—A high-density air shower array for studies on the structure in the cosmic-ray energy spectrum near the knee . Nucl. Instrum. Meth. A . 540 , 311 - 323 ( 2005 ). doi: 10.1016/j.nima.2004.11.025 http://doi.org/10.1016/j.nima.2004.11.025 GRAPES-3—A high-density air shower array for studies on the structure in the cosmic-ray energy spectrum near the knee A. Aab , P. Abreu , M. Aglietta et al ., The pierre auger observatory upgrade-preliminary design report . arXiv:1604.03637 ( 2016 ). The pierre auger observatory upgrade-preliminary design report L.K. Ding , X.H. Ding , C. Fan et al ., Tibet air shower array: results and future plan . J. Phys. Conf. Ser . 120 , 062024 ( 2008 ). doi: 10.1088/1742-6596/120/6/062024 http://doi.org/10.1088/1742-6596/120/6/062024 Tibet air shower array: results and future plan J. Zhao , J. Liu , X.D. Sheng et al ., Design and performances of electromagnetic particle detector for LHAASO-KM2A . Chinese Phys. C . 38 , 036002 ( 2014 ). doi: 10.1088/1674-1137/38/3/036002 http://doi.org/10.1088/1674-1137/38/3/036002 Design and performances of electromagnetic particle detector for LHAASO-KM2A P. Aguiar , E. Casarejos , J.A. Vilan et al ., Geant4-GATE simulation of a large plastic scintillator for muon radiography . IEEE T. Nucl. Sci . 62 , 1233 - 1238 ( 2015 ). doi: 10.1109/TNS.2015.2431297 http://doi.org/10.1109/TNS.2015.2431297 Geant4-GATE simulation of a large plastic scintillator for muon radiography G.I. Britvich , S.K. Chernichenko , A.P. Chubenko et al ., The large scintillation charged particles detector of the Tien-Shan complex “ATHLET” . Nucl. Instrum. Meth. A . 564 , 225 - 234 ( 2006 ). doi: 10.1016/j.nima.2006.03.042 http://doi.org/10.1016/j.nima.2006.03.042 The large scintillation charged particles detector of the Tien-Shan complex “ATHLET” X.L. Qian , X. Wang , H.Y. Sun et al ., Simulation study on cosmic ray background at large zenith angle based on GRANDProto35 coincidence array experiment . Nucl. Sci. Tech . 32 , 1 - 11 ( 2021 ). doi: 10.1007/s41365-020-00841-3 http://doi.org/10.1007/s41365-020-00841-3 Simulation study on cosmic ray background at large zenith angle based on GRANDProto35 coincidence array experiment M. Kandemir , A. Cakir , Simulation and efficiency studies of optical photon transportation and detection with plastic antineutrino detector modules . Nucl. Instrum. Meth. A . 898 , 30 ( 2018 ). doi: 10.1016/j.nima.2018.04.059 http://doi.org/10.1016/j.nima.2018.04.059 Simulation and efficiency studies of optical photon transportation and detection with plastic antineutrino detector modules S. Scheu , H. Kaspar , P. Robmann et al ., Studies on wrapping materials and light collection geometries in plastic scintillators . Nucl. Instrum. Meth. A . 567 , 345 - 349 ( 2006 ). doi: 10.1016/j.nima.2006.05.153 http://doi.org/10.1016/j.nima.2006.05.153 Studies on wrapping materials and light collection geometries in plastic scintillators M. Janecek , W.W. Moses , Simulating scintillator light collection using measured optical reflectance . IEEE T. Nucl. Sci . 57 , 964 ( 2010 ). doi: 10.1109/TNS.2010.2042731 http://doi.org/10.1109/TNS.2010.2042731 Simulating scintillator light collection using measured optical reflectance E. Roncali , S.R. Cherry , Simulation of light transport in scintillators based on 3D characterization of crystal surfaces . PHYS. MED. BIOL . 58 , 2185 ( 2013 ) doi: 10.1088/0031-9155/58/7/2185 http://doi.org/10.1088/0031-9155/58/7/2185 Simulation of light transport in scintillators based on 3D characterization of crystal surfaces G.F. Knoll , T.F. Knoll , T.M. Henderson , Light collection in scintillation detector composites for neutron detection . IEEE T. Nucl. Sci . 35 , 872 - 875 ( 1988 ). doi: 10.1109/23.12850 http://doi.org/10.1109/23.12850 Light collection in scintillation detector composites for neutron detection F.-X. Gentit , Litrani: a general-purpose Monte Carlo program simulating light propagation in isotropic or anisotropic media . Nucl. Instrum. Meth. A . 486 , 35 - 39 ( 2002 ). doi: 10.1016/S0168-9002(02)00671-X http://doi.org/10.1016/S0168-9002(02)00671-X Litrani: a general-purpose Monte Carlo program simulating light propagation in isotropic or anisotropic media S. Agostinelli , J. Allison , K. Amako et al ., GEANT4—a simulation toolkit . Nucl. Instrum. Meth. A . 506 , 250 ( 2003 ). doi: 10.1016/S0168-9002(03)01368-8 http://doi.org/10.1016/S0168-9002(03)01368-8 GEANT4—a simulation toolkit J. Allison , K. Amako , J. Apostolakis et al ., Geant4 developments and applications . IEEE T. Nucl. Sci . 53 , 270 ( 2006 ). doi: 10.1109/TNS.2006.869826 http://doi.org/10.1109/TNS.2006.869826 Geant4 developments and applications G. Santin , D. Strul , D. Lazaro et al ., GATE: A Geant4-based simulation platform for PET and SPECT integrating movement and time management . IEEE T. Nucl. Sci . 50 , 1516 - 1521 ( 2003 ). doi: 10.1109/TNS.2003.817974 http://doi.org/10.1109/TNS.2003.817974 GATE: A Geant4-based simulation platform for PET and SPECT integrating movement and time management W.Z. Xu , Y.F. Liu , Z.Q. Tan et al ., Geant4 simulation of plastic scintillators for a prototype μSR spectrometer . Nucl. Sci. Tech . 24 , 040401 ( 2013 ). doi: 10.13538/j.1001-8042/nst.2013.04.011 http://doi.org/10.13538/j.1001-8042/nst.2013.04.011 Geant4 simulation of plastic scintillators for a prototype μSR spectrometer S. Riggi , R. La Rocca , E. Leonora et al ., GEANT4 simulation of plastic scintillator strips with embedded optical fibers for a prototype of tomographic system . Nucl. Instrum. Meth. A . 624 , 583 - 590 ( 2010 ). doi: 10.1016/j.nima.2010.10.012 http://doi.org/10.1016/j.nima.2010.10.012 GEANT4 simulation of plastic scintillator strips with embedded optical fibers for a prototype of tomographic system Eljen Technology , EJ-200 Plastic Scintillator , http://www.eljentechnology.com/images/products/data_sheets/EJ-200_EJ-204_EJ-208_EJ-212.pdf http://www.eljentechnology.com/images/products/data_sheets/EJ-200_EJ-204_EJ-208_EJ-212.pdf EJ-200 Plastic Scintillator Hamamatsu Photonics , 51 mm (2") Photomultiplier tube assembly R7725 , http://www.hamamatsu.com/resources/pdf/etd/R7723_R7724_R7725_TPMH1315E.pdf http://www.hamamatsu.com/resources/pdf/etd/R7723_R7724_R7725_TPMH1315E.pdf 51 mm (2") Photomultiplier tube assembly R7725 A. Levin , C. Moisan , A more physical approach to model the surface treatment of scintillation counters and its implementation into DETECT . IEEE Nucl. Sci. Conf. R . 2 , 702 ( 1996 ). doi: 10.1109/NSSMIC.1996.591410 http://doi.org/10.1109/NSSMIC.1996.591410 A more physical approach to model the surface treatment of scintillation counters and its implementation into DETECT A. Filevich , P. Bauleo , H. Bianchi et al ., Spectral-directional reflectivity of Tyvek immersed in water . Nucl. Instrum. Meth. A . 423 , 108 ( 1999 ). doi: 10.1016/S0168-9002(98)01194-2 http://doi.org/10.1016/S0168-9002(98)01194-2 Spectral-directional reflectivity of Tyvek immersed in water J. Arteaga-Velázquez , C. Vazquez-Lopez , A. Zepeda , A measurement of the diffuse reflectivity of 1056 Tyvek in air and water . Nucl. Instrum. Meth. A . 553 , 312 ( 2005 ). doi: 10.1016/j.nima.2005.08.030 http://doi.org/10.1016/j.nima.2005.08.030 A measurement of the diffuse reflectivity of 1056 Tyvek in air and water D. Clark , The intrinsic scintillation efficiency of plastic scintillators for 60Co gamma excitation . Nucl. Instrum. Methods . 117 , 295 ( 1974 ). doi: 10.1016/0029-554X(74)90412-1 http://doi.org/10.1016/0029-554X(74)90412-1 The intrinsic scintillation efficiency of plastic scintillators for 60Co gamma excitation C. Hawkes , M. Kuhlen , B. Milliken et al ., Decay time and light yield measurements for plastic scintillating fibers . Nucl. Instrum. Meth. A . 292 , 329 ( 1990 ). doi: 10.1016/0168-9002(90)90388-M http://doi.org/10.1016/0168-9002(90)90388-M Decay time and light yield measurements for plastic scintillating fibers J. Nam , Y. Choi , D. Kim et al ., A detailed Monte Carlo simulation for the Belle TOF system . Nucl. Instrum. Meth. A . 491 , 54 ( 2002 ). doi: 10.1016/S0168-9002(02)01231-7 http://doi.org/10.1016/S0168-9002(02)01231-7 A detailed Monte Carlo simulation for the Belle TOF system X.L. Qian , J.F. Chang , C.F. Feng et al ., Calibration of pulse transit time through a cable for EAS experiments . Chinese Phys. C . 38 , 066202 ( 2014 ). doi: 10.1088/1674-1137/38/6/066202 http://doi.org/10.1088/1674-1137/38/6/066202 Calibration of pulse transit time through a cable for EAS experiments All authors contributed to the study conception and design. Material preparation, data collection and analysis were performed by Xiang-Li Qian, Hui-Ying Sun, Cheng Liu, Xu Wang, and Olivier Martineau-Huynh. The first draft of the manuscript was written by Xiang-Li Qian and all authors commented on previous versions of the manuscript. All authors read and approved the final manuscript. 10.1007/s41365-021-00882-2 2020-12-24 2021-03-19 2021-03-25 2021-05-18 2021-05-01 2021-12-20 © China Science Publishing & Media Ltd. (Science Press), Shanghai Institute of Applied Physics, the Chinese Academy of Sciences, ChineseNuclear Society and Springer Nature Singapore Pte Ltd. 2021 This is a post-peer-review, pre-copyedit version of an article published in Nuclear Science and Techniques. The final authenticated version is available online at: http://dx.doi.org/10.1007/s41365-021-00882-2 http://dx.doi.org/10.1007/s41365-021-00882-2 . 2021 Nuclear Science and Techniques 05 32