Geant4 Simulation Research Articles

Nanodosimetric investigation of the track structure of therapeutic carbon ion radiation. Part 2: Detailed simulation.


A previous study reported nanodosimetric measurements of therapeutic-energy carbon ions penetrating simulated tissue. The results are incompatible with the predicted mean energy of the carbon ions in the nanodosimeter and previous experiments with lower energy monoenergetic beams. The purpose of this study is to explore the origin of these discrepancies.
Approach: Detailed simulations using the Geant4 toolkit were performed to investigate the radiation field in the nanodosimeter and provide input data for track structure simulations, which were performed with a developed version of the PTra code.

Main results: 
The Geant4 simulations show that with the narrow-beam geometry employed in the experiment, only a small fraction of the carbon ions traverse the nanodosimeter and their mean energy is between 12 % and 30 % lower than the values estimated using the SRIM software. Only about one-third or less of these carbon ions hit the trigger detector. The track structure simulations indicate that the observed enhanced ionization cluster sizes are mainly due to coincidences with events in which carbon ions miss the trigger detector. In addition, the discrepancies observed for high absorber thicknesses of carbon ions traversing the target volume could be explained by assuming an increase in thickness or interaction cross-sections in the order of 1 %.

Significance: 
The results show that even with strong collimation of the radiation field, future nanodosimetric measurements of clinical carbon ion beams will require large trigger detectors to register all events with carbon ions traversing the nanodosimeter. Energy loss calculations of the primary beam in the absorbers are insufficient and should be replaced by detailed simulations when planning such experiments. Uncertainties of the interaction cross-sections in simulation codes may shift the Bragg peak position.

High-dimensional Bayesian likelihood normalisation for CRESST's background model

Using CaWO4 crystals as cryogenic calorimeters, the CRESST experiment searches for nuclear recoils caused by the scattering of potential Dark Matter particles. A reliable identification of a potential signal crucially depends on an accurate background model. In this work, we introduce an improved normalisation method for CRESST's model of electromagnetic backgrounds, which is an important technical step towards developing a more accurate background model. Spectral templates based on Geant4 simulations are normalised via a Bayesian likelihood fit to experimental background data. Contrary to our previous work, no explicit assumption of partial secular equilibrium is required a priori, which results in a more robust and versatile applicability. This new method also naturally considers the correlation between all background components. Due to these purely technical improvements, the presented method has the potential to explain up to 82.7 % of the experimental background within [1 keV,40 keV], an improvement of at most 18.6 % compared to our previous method. The actual value is subject to ongoing validations of the included physics.

Advancing Proton Therapy: A Review of Geant4 Simulation for Enhanced Planning and Optimization in Hadron Therapy

Abstract Proton therapy is a cancer treatment method that uses high-energy proton beams to target and destroy cancer cells. In recent years, the use of proton therapy in cancer treatment has increased due to its advantages over traditional radiation methods, such as higher accuracy and reduced damage to healthy tissues. For accurate planning and delivery of proton therapy, advanced software tools are needed to model and simulate the interaction between the proton beam and the patient’s body. One of these tools is the Monte Carlo simulation software called Geant4, which provides accurate modeling of physical processes during radiation therapy. The purpose of this study is to investigate the effectiveness of the Geant4 toolbox in proton therapy in the conducted research. This review article searched for published articles between 2002 and 2023 in reputable international databases including Scopus, PubMed, Scholar, Google Web of Science, and ScienceDirect. Geant4 simulations are reliable and accurate and can be used to optimize and evaluate the performance of proton therapy systems. Obtaining some data from experiments carried out in the real world is very effective. This makes it easy to know how close the values obtained from simulations are to the behavior of ions in reality.

Unveiling the origin of XMM-Newton soft proton flares

Context. Low-energy (<300 keV) protons entering the field of view of XMM-Newton can scatter with the X-ray mirror surface and reach the focal plane. They are observed in the form of a sudden increase in the background level, the so-called soft proton flares, affecting up to 40% of the mission observing time. Soft protons can hardly be disentangled from true X-ray events and cannot be rejected on board. Aims. All future high throughput grazing incidence X-ray telescopes operating outside the radiation belts are potentially affected by soft proton-induced contamination that must be foreseen and limited since the design phase. In-flight XMM-Newton’s observations of soft protons represent a unique laboratory to validate and improve our understanding of their interaction with the mirror, optical filters, and X-ray instruments. At the same time, such models would link the observed background flares to the primary proton population encountered by the telescope, converting XMM-Newton into a monitor for soft protons. Methods. We built a Geant4 simulation of XMM-Newton, including a verified mass model of the X-ray mirror, the focal plane assembly, and the EPIC MOS and pn-CCDs. Analytical computations and, when available, laboratory measurements collected from literature were used to verify the correct modelling of the proton scattering and transmission to the detection plane. Similarly to the instrument X-ray response, we encoded the energy redistribution and proton transmission efficiency into a redistribution matrix file (RMF), mapping the probability that a proton from 2 to 300 keV is detected in a certain detector channel, and an auxiliary response file (ARF), storing the grasp towards protons. Both files were formatted according to the standard NASA calibration database and any compliant X-ray data analysis tool can be used to simulate or analyse soft proton-induced background spectra. An overall systematic uncertainty of 30% was assumed on the basis of the estimated accuracy of the mirror geometry and transmission models. Results. For the validation, three averaged soft proton spectra, one for each filter configuration, were extracted from a collection of 13 years of MOS observations of the focused non X-ray background and analysed with Xspec. A similar power-law distribution is found for the three filter configurations, plus black-body-like emission below tens of keV used as a correction factor, based on the dedicated spectral analysis of 55 in-flight proton flares presented in Paper II. The best-fit model is in agreement with the power-law distribution predicted from independent measurements for the XMM-Newton orbit, spent mostly in the magnetosheath and nearby regions. For the first time we are able to link detected soft proton flares with the proton radiation environment in the Earth’s magnetosphere, while proving the validity of the simulation chain in predicting the background of future missions. Benefiting from this work and contributions from the Athena instrument consortia, we also present the response files for the Athena mission and updated estimates for its focused charged background.

Unveiling the origin of XMM-Newton soft proton flares

Context. Low-energy (< 300 keV) protons entering the field of view of the XMM-Newton telescope scatter with the X-ray mirror surface and might reach the X-ray detectors on the focal plane. They manifest in the form of a sudden increase in the rates, usually referred to as soft proton flares. By knowing the conversion factor between the soft proton energy and the deposited charge on the detector, it is possible to derive the incoming flux and to study the environment of the Earth magnetosphere at different distances, given the wide and elliptical XMM-Newton orbit. Thanks to detailed Geant4 simulations, we were able to build specific soft proton response matrices for MOS and PN. Aims. In this second paper, we present the results of testing these matrices with real data for the first time, while also exploring the seasonal and solar activity effect on the proton environment. The selected spectra are relative to 55 simultaneous MOS and PN observations with flares raised in four different temporal windows: December-January and July-August of 2001-2002 (solar maximum) and 2019-2020 (solar minimum). Methods. We selected and extracted the flare mean spectra and count rates in the 2–11.5 keV energy range for the four epochs. After investigating the rate variations among the MOS1, MOS2, and PN instruments, we fit the X-ray spectra using XSPEC and the proton response matrices. The best-fitting parameters derived for the three instruments were compared in order to obtain the systematic errors. Results. There is no seasonal or solar activity effect on the soft proton mean count rates, but we find large discrepancies in the instrument cross-correlations across the 20 years of satellite operations. In 2001-2002, after a few years of operation, the MOS1 and MOS2 rates are similar, and about 20% with regard to the PN ones. After 20 years, PN does not present any variation in its response, while MOS1 suffers a reduction of ∼30%, in addition to the 30% loss due to the damage of two CCDs, and MOS2 is affected by an even worse degradation (70%). The main result of the spectral analysis is that the physical model representative of the proton spectra at the input of the telescope is a power law. However, a second and phenomenological component is necessary to take into account imprecision in the generation of the matrices at softer (< 5 keV) energies. This component contributes for 21% for the MOS and 5% for the PN to the total flux in the 2–5 keV energy range. Conclusions. This study, which is the first application of the soft proton response matrices to real data, shows coherent results between detectors and allows us to estimate systematic uncertainties in the measured spectra of 3% between the two MOS detectors and of 24% between MOS and PN, together with a systematic in the input flux of about a factor of two. They are all likely due to uncertainties in the proton transmission models, with the presence of additional passive material in front of the front-illuminated MOS, and element deposition on its electrode structure across the mission life. Dedicated studies and laboratory measurements are required for improving the accuracy of the proton response files.

Record high counting rate of positron annihilation lifetime spectrometer achieved by [formula omitted] coincidence

Conventional positron annihilation lifetime (PAL) spectrometers, which employ the γ-γ coincidence with two perpendicularly-positioned detectors, have a typical counting rate of 100–300 counts per second (cps). To increase the counting rate, according to the optimized structural parameters using Geant4 simulation, a 22Na positron source (∼ 1.665 MBq), a silicon photomultiplier (SiPM), and digital waveform technology are utilized for the first time to a β+-γ coincidence PAL spectrometer. After the optimization and verification of accuracy and stability, this spectrometer can achieve a time resolution of about 206 ps, and a record high effective counting rate of approximately 19000 cps (∼ 11000 cps/MBq), which is two orders of magnitude greater than those of conventional PAL spectrometers. The significant shortening of measurement time for a single PAL spectrum enables us to observe sub-minute-scale evolution of microstructure during rapid physical and chemical processes in the future.

Significance of the proton energy loss mechanism to gold nanoparticles in proton therapy: a Geant4 simulation

The biological effectiveness in proton therapy is only slightly greater than of the treatment by X-ray and hence, many researches have suggested the use of gold nanoparticles for increasing the ionization interactions to produce more secondary electrons and elevate the yield of DNA damage. But the ionization interactions also lead to protons energy loss inside the nanoparticles. The present study shows that, the protons slowed-down by High-Z nanoparticles are responsible for dose enhancement rather than the produced secondary electrons. To this purpose, using Geant4 Monte Carlo tool, one million nanoparticles distributed in a proton irradiated volume and variation of the proton’s spectra and the dose related to this variation has been demonstrated. It was found that, the elevation in proton’s LET values when passing through the gold nanoparticles will lead to a more significant dose enhancement than the increased dose due to the extra secondary electrons. Also, it was found that, the mechanism of protons slowing-down by gold nanoparticles has another useful aspect in proton therapy in which, the dose leakage to surrounding healthy tissues will be reduced which must be considered in future investigations more precisely.

Contributions of cosmic-ray components to the HPGe gamma spectrometer background spectrum within the 0°–45° Zenith angle range

This study investigates the contributions of various cosmic-ray-induced components to the energy response of a High-Purity Germanium (HPGe) gamma spectrometer within the 0°–45° zenith angle range. The analysis focuses on muons, neutrons, protons, electrons, positrons, and photons. It also examines the impact of particle showers induced by muons interacting with the lead chamber and reaching the HPGe detector. Utilizing the GEANT4 simulation toolkit, we provide a detailed examination of these components' influence on the spectrum. The results confirm that within this zenith angle range, muons account for approximately 91% of the recorded cosmic-ray induced energy spectrum and contribute 54.6% to the overall cosmic-ray-induced background spectrum. The study also highlights the significant role of showers, especially those resulting from muon interactions with the lead shielding, in shaping the low-energy spectrum below 3 MeV.

Effect of Hydrogen on Albedo Neutrons in the Lunar Surface

The remote sensing observations of the lunar soil by using neutron spectroscopy as well as infrared spectroscopy in various studies suggest the presence of hydrogen in the lunar surface, particularly around the lunar polar regions. Additionally, the cosmic ray protons incident on the lunar surface induce the generation of secondary neutrons depending on the composition of the surface soil and the associated density, a certain number of which, called albedo neutrons, leak from the lunar regolith. Motivated by the interaction of the cosmic ray protons with the lunar surface in the absence and presence of hydrogen, a series of GEANT4 simulations are employed in the current study to unveil the influence of hydrogen on the energy spectrum of the albedo neutrons that escape from the lunar surface. Initially, a discrete 69-bin proton energy spectrum between 0.4 and 115 GeV based on the PAMELA spectrometer is implemented into GEANT4. Subsequently, a single-volume lunar surface of 2-m thickness is constructed where it is assumed that a single layer of either 1.6 or 1.93 g/cm^3 constitutes the lunar regolith. The material composition in the present study includes 43.7 wt.% oxygen, 0.3 wt.% sodium, 5.6 wt.% magnesium, 9 wt.% aluminium, 21.1 wt.% silicon, 8.5 wt.% calcium, 1.5 wt.% titanium, 0.1 wt.% manganese, and 10.2 wt.% iron in accordance with another study. Then, 0.1 wt.% hydrogen is introduced by replacing oxygen in order to assess the impact of hydrogen on the secondary neutrons. In the wake of irradiating a lunar surface of 3*2*3 m^3 with a planar vertical PAMELA proton beam of 20*20 cm^2, the depth profile of the generated neutrons as well as the corresponding initial energy spectrum is primarily obtained in the absence and presence of 0.1 wt.% hydrogen by voxelizing the entire geometry with 100 cells of 2-cm thickness. Next, a surface detector is placed at the top of the lunar surface in order to collect the albedo neutrons in both cases, and the energy spectra of the albedo neutrons are acquired in terms of thermal, epithermal, and fast neutrons. From the GEANT4 simulations in this study, it is shown that the presence of 0.1 wt.% hydrogen is observable in each energy regime of the albedo neutrons at the lunar surface, thereby providing an indication about the elemental variation of the lunar soil.

Research on B4C/PEEK Composite Material Radiation Shielding

There are various types of charged particles in the space environment, which can cause different types of radiation damage to materials and devices, leading to on-orbit failures and even accidents for spacecraft. Developing lightweight and efficient radiation-shielding materials is an effective approach to improving the inherent protection of spacecraft. The protective performance of different materials against proton and electron spectra in the Earth’s radiation belts is evaluated using a Geant4 simulation. Based on the simulation results, suitable hardening components were selected to design composite materials, and B4C/PEEK composites with different B4C contents were successfully prepared. The experimental results demonstrate that the simulated and experimental results for the electron, proton and neutron shielding performance of the B4C/PEEK composites are consistent. These composites exhibit excellent radiation shielding capabilities against electrons, protons and neutrons, and the radiation protection performance improves with increasing B4C content in the B4C/PEEK composite materials.

Radiation shielding transmission data for modern digital radiography.

Radiography is one of the most widely used x-ray imaging modalities. In National Council on Radiation Protection and Measurements (NCRP) Report No. 147, transmission data for radiographic systems were evaluated on those installed before 2000. For x-ray systems (except intraoral dental) manufactured on or after June 10, 2006, the U.S. required minimum half-value layer (HVL) was increased; for example, 2.9 (2.3) mm Al at 80kV, where the value in parenthesis denotes the earlier requirement before the above date. To calculate the transmission of the broad x-ray beam of modern digital radiography (DR) through shielding materials. X-ray beam HVLs on two DR systems (Agfa DR 600, GE Revolution XR/d) were measured in 10kV increments between 60 and 120kV, with a calibrated ionization chamber (Radcal model 10×5-60) and varying thickness of aluminum 1100 plates. Monte Carlo (Geant4) simulation was performed to calculate the transmission of broad x-ray beams through lead, concrete, gypsum, and steel, with x-ray HVLs matching those of the DR 600 at two beam filtrations (default, 1mm Al plus 0.1mm Cu added filtration). The transmission data were fitted to the Archer equation. HVLs on two DR systems with default beam filtration were consistent (median difference, 2.1%; maximum difference, 5.5%). An additional beam filtration option (1mm Al plus 0.1mm Cu) on the DR 600 substantially increased HVLs by 45.2%-61.2%. Transmission fitting parameters were provided for seven tube voltages (60-120kV) at two beam filtrations. This work presents transmission data for modern DR systems, indicating increased x-ray beam filtration compared to the primary x-ray beam in NCRP Report No. 147. The updated transmission data can enhance structural shielding evaluations at individual tube voltages or with a workload distribution.

Optimization of the Electron Spectrometer Telescope geometry for the PRESET satellite through Monte Carlo simulation

The Electron Spectrometer Telescope (EST) aboard the PRESET satellite mission aims to measure the pitch angle distribution of 0.3–4 MeV electrons in the outer Van Allen belts. The PRESET satellite is planned to launch into a sun-synchronous low Earth orbit in Q2, 2026. We present comprehensive Monte Carlo simulations to optimize the design and response of the EST instrument. The EST consists of a stack of silicon strip detectors and a collimator to take pitch-angle dependent electron spectral measurements with a target angular resolution of 6° while rejecting the proton, heavy ion and gamma detection events. Various collimator designs and detector stacking configurations are investigated using the Geant4 code to deduce an optimal design and configuration. For validation of the Geant4 simulations, a benchmark test spectrometer was set up using a stack of four silicon detectors and a good agreement between the measured and simulated responses was found for a beta spectrum from a 90Sr/90Y source. The collimator design was optimized by adjusting total length, aperture size, number of view-ports and collimator materials. The optimum aperture size and the number of view-ports were determined by varying them and selecting the best cases that meet the angular resolution and the geometric factor requirements. Moreover, to address the under-utilization of the front position-sensitive detector encountered in a typical pin-hole collimator, two offset collimators with tilted axes were employed. Extensive simulations were carried out by varying the number of silicon detectors and the thickness of each detector to optimize the configuration of the silicon detector stack. An optimum thickness of the front detector was found to be 140 µm, which can minimize the perturbation counts caused by the proton detection events. For the other detectors, a stack of four 1.5 mm thick detectors was chosen to achieve a good sensitivity to low energy electrons at a tolerable complexity of the system. The angular response simulated for the optimum design of the EST showed a resolution of 5.5°, which is slightly better than the target resolution of 6°. The response matrices of the final design simulated for isotropic electron and proton fields showed a high geometric factor, which is expected to produce an average electron counting rate of 230 cps within the trapped region with a negligible contamination of the electron spectrum by protons except for a mild contamination by the 1.0–1.1 MeV protons.

Scintillator-SiPM detector for tracking and energy deposition measurements

An innovative particle detector that offers a compelling combination of cost-effectiveness and high accuracy is introduced. The detector features plastic scintillators paired with a sparse arrangement of SiPMs, strategically positioned within a unique opto-mechanical framework. This configuration delivers precise measurements of spatial impact position and energy deposition of impinging particles. The manuscript describes the detector’s physical model complemented by an analytical representation. These calculations underpin a numerical algorithm, facilitating the estimation of particle impingement position and energy deposition. The results of the numerical calculations are compared with the output of GEANT4 simulations and evaluated by rigorous laboratory testing. An array of these detectors, intended for deployment in a spaceborne experiment, underwent detailed design, manufacturing, and testing. Their performance and alignment with the physical model were validated through meticulously conducted ground-based laboratory experiments, conclusively affirming the detector’s properties.

(G4-SEOPTIM): A GEANT4 application for shielded enclosure, designing and shielding optimizing: The case of a new Moroccan shielded enclosure

PurposeThis study aims to develop and validate a Geant4 simulation application using the latest G4RadioactiveDecayPhysics v5.1 library to optimize the shielding design of a Moroccan shielded enclosure. The goal is to ensure effective radiation protection by comparing simulation results with initial geometric calculations based on the attenuation law (Beer-Lambert Law) and implementing design optimizations to enhance safety and reduce costs. MethodsThe study began with a prototype design based on linear attenuation calculations, which were costly and offered limited patient safety. Using the Geant4 Monte Carlo toolkit, we simulated the decay of a Fluorine-18 source within the shielded enclosure to evaluate dose rate distributions. The simulation's accuracy was validated through comparison with experimental measurements, and the results informed the geometric optimization of the shielding enclosure. ResultsThe Geant4 simulations demonstrated that the current shielding design effectively reduces radiation doses to acceptable levels. However, the simulations also identified opportunities for further optimization of the enclosure's geometry, allowing for more precise and efficient use of shielding materials while maintaining safety standards. ConclusionThe study confirms the utility of Geant4 for optimizing radiation shielding designs. While the initial dimensions provided adequate protection, the simulation enabled more precise geometric optimization. This approach not only enhances radiation protection but also informs the design and construction of more efficient and cost-effective shielding enclosures. The adjustments made to the enclosure's faces significantly improved radiation protection, reducing the dose rate surrounding the shielded enclosure to acceptable levels.

Particle identification capability of a homogeneous calorimeter composed of oriented crystals

Recent studies have shown that the electromagnetic shower induced by a high-energy electron, positron or photon incident along the axis of an oriented crystal develops in a space more compact than the ordinary. On the other hand, the properties of the hadronic interactions are not affected by the lattice structure. This means that, inside an oriented crystal, the natural difference between the hadronic and the electromagnetic shower profile is strongly accentuated. Thus, a calorimeter composed of oriented crystals could be intrinsically capable of identifying more accurately the nature of the incident particles, with respect to a detector composed only of non-aligned crystals. Since no oriented calorimeter has ever been developed, this possibility remains largely unexplored and can be investigated only by means of numerical simulations. In this work, we report the first quantitative evaluation of the particle identification capability of such a calorimeter, focusing on the case of neutron-gamma discrimination. We demonstrate through Geant4 simulations that the use of oriented crystals significantly improves the performance of a Random Forest classifier trained on the detector data. This work is a proof that oriented calorimeters could be a viable option for all the environments where particle identification must be performed with a very high accuracy, such as future high-intensity particle physics experiments and satellite-based γ-ray telescopes.

Development of a large NaI(Tl) detection system

A large NaI(Tl) detection system was developed to address the increasing demand for high detection efficiency and high energy resolution detectors across various applications. Through Geant4 simulation and experimental research, critical characteristics of the detector prototype were examined, encompassing energy response, energy nonlinearity, and energy resolution. Given the substantial impact of variations in Tl concentration and the complex process of photon propagation on the detection performance of large NaI(Tl) crystals, the detection efficiency and energy resolution with the radioactive source placed at different incident spots were studied to evaluate its spatial response uniformity. The results demonstrate that the detector's energy response range spans from 50 keV to 3 MeV with an energy nonlinearity of 1.5%. The developed large NaI(Tl) detection system achieved an overall energy resolution of 7.93%, along with excellent uniformity in its resolution and detection efficiency, satisfying the application requirements for large dimensions and high detection efficiency.

High brightness betatron x-ray source driven by chirped laser pulses

We demonstrate high brightness betatron x-ray generation from a chirped laser pulses driven plasma accelerator. It is shown that positively chirped laser pulse leads to the initiation and enhancement of collective oscillations of electrons inside plasma bubbles, due to associated pulse front tilt (PFT). The PFT causes transverse drift of the bubbles with respect to the laser axis, which results in high brightness x-ray generation. At an optimum chirp, enhanced x-ray emission of >108 photons/pulse/sr in 0.1% BW with a critical energy of ∼18 keV was observed by a factor >2 in comparison to the case of no chirp. The role of collective oscillation in enhancing x-ray emission is also validated in the Geant4 simulations.

Optimization of the Pixel Design for Large Gamma Cameras Based on Silicon Photomultipliers.

Most single-photon emission computed tomography (SPECT) scanners employ a gamma camera with a large scintillator crystal and 50-100 large photomultiplier tubes (PMTs). In the past, we proposed that the weight, size and cost of a scanner could be reduced by replacing the PMTs with large-area silicon photomultiplier (SiPM) pixels in which commercial SiPMs are summed to reduce the number of readout channels. We studied the feasibility of that solution with a small homemade camera, but the question on how it could be implemented in a large camera remained open. In this work, we try to answer this question by performing Geant4 simulations of a full-body SPECT camera. We studied how the pixel size, shape and noise could affect its energy and spatial resolution. Our results suggest that it would be possible to obtain an intrinsic spatial resolution of a few mm FWHM and an energy resolution at 140 keV close to 10%, even if using pixels more than 20 times larger than standard commercial SiPMs of 6 × 6 mm2. We have also found that if SiPMs are distributed following a honeycomb structure, the spatial resolution is significantly better than if using square pixels distributed in a square grid.

3D printed microcrystalline CsI:Tl composite scintillating thin films for X-ray imaging

The utilization of additive manufacturing techniques, especially Digital Light Printing (DLP), in fabricating CsI:Tl scintillator films demonstrates considerable potential for streamlining the production of scintillators tailored for X-ray imaging applications. This research focuses on the fabrication of CsI:Tl-based composite plastic scintillator thin films. In this study, circular films measuring 1-inch in diameter and 0.1 mm & 0.2 mm in thickness are being produced and tested for gamma photon counts under alpha and gamma radiation. To establish the stopping power range of the films, a Monte Carlo based GEANT4 simulation has been carried out. Additionally, investigations into their suitability for X-ray imaging applications are being conducted, revealing the spatial resolution of the films (0.2 and 0.1 mm) between 100 and 130 μm and 1.26 lp/mm with a contrast range of 4.0–12.3 %. The observed decrease in spatial resolution and contrast for the 0.2 mm thick film is attributed to the thickness increase exacerbating the scattering phenomenon while simultaneously enhancing the X-ray stopping power. This highlights the significance of inherent trade-off between maximizing spatial resolution and compromising light yield of 0.1 mm films compared to the 0.2 mm thick film. By utilizing 3D printing, this approach offers a cost-effective and time-efficient method for producing thin-film scintillators with enhanced flexibility and customization options compared to conventional methods.

SSLG4: A novel scintillator simulation library for Geant4

This study introduces a new Scintillator Simulation Library called SSLG4 for the Geant4 Monte Carlo simulation package. With SSLG4, we aim to enhance efficiency and accelerate progress in optical simulations within the Geant4 framework by simplifying scintillator handling and providing a rich repository of scintillators. The SSLG4 enables users to quickly include predefined scintillator materials in their simulations without requiring manual definition. The library initially contains 68 scintillators, consisting of 58 organic and 10 inorganic types. Most of these scintillators are selected from the catalogs of several scintillator manufacturers, notably Eljen and Luxium. Other scintillators are included based on their widespread use across various physics domains. The library stores optical data of scintillators in ASCII files with .mac and .txt extensions, enabling users to add, remove, or modify properties of scintillators at runtime of their applications. In addition, we made all the scintillator data available in the library on a dedicated page of our website to ensure convenient access for all users. Program summaryProgram title: SSLG4CPC Library link to program files:https://doi.org/10.17632/3zbwr5wf7z.1Developer's repository link:https://github.com/mkandemirr/SSLG4, https://neutrino.erciyes.edu.tr/SSLG4/Licensing provisions: GNU General Public License 3Programming language: C++External routines/libraries: Geant4, CMake, OPSimNature of problem: Defining a new scintillator in Geant4 is a cumbersome process for some users due to three main reasons: (1) It requires a lot of data input from users, (2) collecting the scintillator data requires an extensive literature review, and (3) the collected data needs to be converted into the desired format. In addition, the interfaces provided to define a scintillator direct users to embed scintillator data into their source code, resulting in increased code complexity, reduced code readability, and an inefficient working environment.Solution method: To solve the problems mentioned above, developing and introducing a new library consisting of fully parameterized and ready-to-use scintillators would greatly increase the useability of the Geant4 simulation package for scintillator studies and interest a wide range of scientific communities.