Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Detector simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4

Accurate contrast determination for X-ray speckle visibility spectroscopy

X-ray speckle visibility spectroscopy using X-ray free-electron lasers has long been proposed as a probe of fast dynamics in noncrystalline materials. In this paper, numerical modeling is presented to show how the data interpretation of visibility spectroscopy can be impacted by the nonidealities of real-life X-ray detectors. Using simulated detector data, this work provides a detailed analysis of the systematic errors of several contrast extraction algorithms in the context of low-count-rate X-ray speckle visibility spectroscopy and their origins are discussed. Here, it was found that the finite detector charge cloud and pixel size lead to an unavoidable `degeneracy' in photon position determination, and that the contrasts extracted using different algorithms can all be corrected by a simple linear model. The results suggest that experimental calibration of the correction coefficient at the count rate of interest is possible and essential. This allows computationally lightweight algorithms to be implemented for on-the-fly analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

GeantV: Results from the Prototype of Concurrent Vector Particle Transport Simulation in HEP

Full detector simulation was among the largest CPU consumers in all CERN experiment software stacks for the first two runs of the Large Hadron Collider. In the early 2010s, it was projected that simulation demands would scale linearly with increasing luminosity, with only partial compensation from increasing computing resources. The extension of fast simulation approaches to cover more use cases that represent a larger fraction of the simulation budget is only part of the solution, because of intrinsic precision limitations. The remainder corresponds to speeding up the simulation software by several factors, which is not achievable by just applying simple optimizations to the current code base. In this context, the GeantV R&D project was launched, aiming to redesign the legacy particle transport code in order to benefit from features of fine-grained parallelism, including vectorization and increased locality of both instruction and data. This paper provides an extensive presentation of the results and achievements of this R&D project, as well as the conclusions and lessons learned from the beta version prototype.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Advancing set-conditional set generation: Diffusion models for fast simulation of reconstructed particles

The computational intensity of detector simulation and event reconstruction poses a significant difficulty for data analysis in collider experiments. This challenge inspires the continued development of machine learning techniques to serve as efficient surrogate models. We propose a fast emulation approach that combines simulation and reconstruction. In other words, a neural network generates a set of reconstructed objects conditioned on input particle sets. To make this possible, we advance set-conditional set generation with diffusion models. Using a realistic, generic, and public detector simulation and reconstruction package (COCOA), we show how diffusion models can accurately model the complex spectrum of reconstructed particles inside jets.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Organic Scintillator Detector Response Simulations with DRiFT

Here, this work presents the organic scintillation simulation capabilities of DRiFT, a post-processing Detector Response Function Toolkit for MCNPR output. DRiFT is used to create realistic scintillation detector response functions to incident neutron and gamma mixed- field radiation. As a post-processing tool, DRiFT leverages the extensively validated radiation transport capabilities of MCNPR ® 6, which also provides the ability to simulate complex sources and geometries. DRiFT is designed to be flexible, it allows the user to specify scintillator material, PMT type, applied PMT voltage, and quenching data used in simulations. The toolkit's capabilities, which include the generation of pulse shape discrimination plots and full-energy detector spectra, are demonstrated in a comparison of measured and simulated neutron contributions from 252 Cf and PuBe, and photon spectra from 22 Na and 228 Th sources. DRiFT reproduced energy resolution effects observed in EJ-301 measurements through the inclusion of scintillation yield variances, photon transport noise, and PMT photocathode and multiplication noise.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Using standard calibrated geometries to characterize a coaxial high purity germanium gamma detector for Monte Carlo simulations

A detector model optimization procedure based on matching Monte Carlo simulations with measurements for two experimentally calibrated sample geometries which are frequently used in radioactivity measurement laboratories results in relative agreement within 5% between simulated and measured efficiencies for a high purity germanium detector. The optimization procedure indicated that the increase in dead layer thickness is largely responsible for a detector efficiency decrease in time. The optimized detector model allows Monte Carlo efficiency calibration for all other samples of which the geometry and bulk composition is known. The presented method is a competitive and economic alternative to more elaborate detector scanning methods and results in a comparable accuracy.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

SU-E-I-67: X-Ray Fluorescence for Energy Response Calibration of a Photon Counting Detector: A Simulation Study

Purpose: To investigate the feasibility of energy calibration and energy response characterization of a photon counting detector using x-ray fluorescence. Methods: A comprehensive Monte Carlo simulation study was done to investigate the influence of various geometric components on the x-ray fluorescence measurement. Different materials, sizes, and detection angles were simulated using Geant4 Application for Tomographic Emission (GATE) Monte Carlo package. Simulations were conducted using 100 kVp tungsten-anode spectra with 2 mm Al filter for a single pixel cadmium telluride (CdTe) detector with 3 × 3 mm2 in detection area. The fluorescence material was placed 300 mm away from both the x-ray source and the detector. For angular dependence measurement, the distance was decreased to 30 mm to reduce the simulation time. Compound materials, containing silver, barium, gadolinium, hafnium, and gold in cylindrical shape, were simulated. The object size varied from 5 to 100 mm in diameter. The angular dependence of fluorescence and scatter were simulated from 20° to 170° with an incremental step of 10° to optimize the fluorescence to scatter ratio. Furthermore, the angular dependence was also experimentally measured using a spectrometer (X-123CdTe, Amptek Inc., MA) to validate the simulation results. Results: The detection angle between 120° to 160° resulted in more optimal x-ray fluorescence to scatter ratio. At a detection angle of 120°, the object size did not have a significant effect on the fluorescence to scatter ratio. The experimental results of fluorescence angular dependence are in good agreement with the simulation results. The Kα and Kβ peaks of five materials could be identified. Conclusion: The simulation results show that the x-ray fluorescence procedure has the potential to be used for detector energy calibration and detector response characteristics by using the optimal system geometry.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Nuclear signal simulation applied to gas ionizing chambers

Particle transport codes used in detector simulation allow the calculation of the energy deposited by charged particles produced following an interaction. The pulses temporal shaping is more and more used in nuclear measurement into pulse shape analysis techniques. A model is proposed in this paper to simulate the pulse temporal shaping and the associated noise level thanks to the output track file PTRAC provides by Monte-Carlo particle transport codes. The model has been dedicated to ion chambers and more especially for High Pressure Xenon chambers HPXe where the pulse shape analysis can resolve some issues regarding with this technology as the ballistic deficit phenomenon. The model is fully described and an example is presented as a validation of such full detector simulation. (authors)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Physics validation studies for muon collider detector background simulations

Within the broad discipline of physics, the study of the fundamental forces of nature and the most basic constituents of the universe belongs to the field of particle physics. While frequently referred to as 'high-energy physics,' or by the acronym 'HEP,' particle physics is not driven just by the quest for ever-greater energies in particle accelerators. Rather, particle physics is seen as having three distinct areas of focus: the cosmic, intensity, and energy frontiers. These three frontiers all provide different, but complementary, views of the basic building blocks of the universe. Currently, the energy frontier is the realm of hadron colliders like the Tevatron at Fermi National Accelerator Laboratory (Fermilab) or the Large Hadron Collider (LHC) at CERN. While the LHC is expected to be adequate for explorations up to 14 TeV for the next decade, the long development lead time for modern colliders necessitates research and development efforts in the present for the next generation of colliders. This paper focuses on one such next-generation machine: a muon collider. Specifically, this paper focuses on Monte Carlo simulations of beam-induced backgrounds vis-a-vis detector region contamination. Initial validation studies of a few muon collider physics background processes using G4beamline have been undertaken and results presented. While these investigations have revealed a number of hurdles to getting G4beamline up to the level of more established simulation suites, such as MARS, the close communication between us, as users, and the G4beamline developer, Tom Roberts, has allowed for rapid implementation of user-desired features. The main example of user-desired feature implementation, as it applies to this project, is Bethe-Heitler muon production. Regarding the neutron interaction issues, we continue to study the specifics of how GEANT4 implements nuclear interactions. The GEANT4 collaboration has been contacted regarding the minor discrepancies in the neutron interaction cross sections for boron. While corrections to the data files themselves are simple to implement and distribute, it is quite possible, however, that coding changes may be required in G4beamline or even in GEANT4 to fully correct nuclear interactions. Regardless, these studies are ongoing and future results will be reflected in updated releases of G4beamline.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Rapid Monte Carlo simulation of detector DQE(f)

Purpose: Performance optimization of indirect x-ray detectors requires proper characterization of both ionizing (gamma) and optical photon transport in a heterogeneous medium. As the tool of choice for modeling detector physics, Monte Carlo methods have failed to gain traction as a design utility, due mostly to excessive simulation times and a lack of convenient simulation packages. The most important figure-of-merit in assessing detector performance is the detective quantum efficiency (DQE), for which most of the computational burden has traditionally been associated with the determination of the noise power spectrum (NPS) from an ensemble of flood images, each conventionally having 10{sup 7} − 10{sup 9} detected gamma photons. In this work, the authors show that the idealized conditions inherent in a numerical simulation allow for a dramatic reduction in the number of gamma and optical photons required to accurately predict the NPS. Methods: The authors derived an expression for the mean squared error (MSE) of a simulated NPS when computed using the International Electrotechnical Commission-recommended technique based on taking the 2D Fourier transform of flood images. It is shown that the MSE is inversely proportional to the number of flood images, and is independent of the input fluence provided that the input fluence is above a minimal value that avoids biasing the estimate. The authors then propose to further lower the input fluence so that each event creates a point-spread function rather than a flood field. The authors use this finding as the foundation for a novel algorithm in which the characteristic MTF(f), NPS(f), and DQE(f) curves are simultaneously generated from the results of a single run. The authors also investigate lowering the number of optical photons used in a scintillator simulation to further increase efficiency. Simulation results are compared with measurements performed on a Varian AS1000 portal imager, and with a previously published simulation performed using clinical fluence levels. Results: On the order of only 10–100 gamma photons per flood image were required to be detected to avoid biasing the NPS estimate. This allowed for a factor of 10{sup 7} reduction in fluence compared to clinical levels with no loss of accuracy. An optimal signal-to-noise ratio (SNR) was achieved by increasing the number of flood images from a typical value of 100 up to 500, thereby illustrating the importance of flood image quantity over the number of gammas per flood. For the point-spread ensemble technique, an additional 2× reduction in the number of incident gammas was realized. As a result, when modeling gamma transport in a thick pixelated array, the simulation time was reduced from 2.5 × 10{sup 6} CPU min if using clinical fluence levels to 3.1 CPU min if using optimized fluence levels while also producing a higher SNR. The AS1000 DQE(f) simulation entailing both optical and radiative transport matched experimental results to within 11%, and required 14.5 min to complete on a single CPU. Conclusions: The authors demonstrate the feasibility of accurately modeling x-ray detector DQE(f) with completion times on the order of several minutes using a single CPU. Convenience of simulation can be achieved using GEANT4 which offers both gamma and optical photon transport capabilities.

61 RADIATION PROTECTION AND DOSIMETRY↗

Evaluation of simulated HPGe detector efficiencies in OpenMC compared to MCNP

In late 2023, the open-source radiation transport code OpenMC introduced a pulse-height tally (PHT) feature, enabling users to track the total energy deposited by individual photons in OpenMC cells. This function represents an important improvement of the OpenMC code because it allows users to simulate the response function of a gamma detector without relying on closed-source alternatives. Despite this, limited work has been published evaluating OpenMC as a radiation transport code for simulating gamma spectroscopy experiments. This study attempts to demonstrate the usefulness of OpenMC in this space by directly comparing its PHT output to MCNP (a trusted industry-standard Monte Carlo code). In the first half of this study, a radiation transport experiment is described in which a detector with a complex internal geometry is exposed to various gamma-emitting isotopes ( 133 Ba, 137 Cs, and 60 Co) over a range of distances. The setup of this experiment was modeled in MCNP with sufficient detail to capture the efficiency characteristics of a high purity germanium (HPGe) detector for the primary gammas of each isotope. After good agreement between the MCNP model and experiment was achieved, an identical model was produced in OpenMC to allow direct comparison between the full energy peak (FEP) values produced in OpenMC and MCNP. The results show strong agreement between OpenMC and MCNP across the full range of tested energies, with each model’s FEP values typically within 2% of each other and most FEP areas within 4% of experimental data. Worse agreement was seen between the Monte Carlo codes and experiment below 300 keV (an expected result). For the 662 keV line of 137 Cs, both codes were found in poor agreement with experiment and each other over the full range of distances tested (possibly indicating an error with the 137 Cs experimental data). Ignoring the anomalous results of the 662 keV line, all other data show good qualitative and quantitative agreement between MCNP and OpenMC. This result demonstrates the accuracy of OpenMC’s PHT feature for spectroscopic applications in which detector efficiency is a primary concern.

07 - ISOTOPES AND RADIATION SOURCES↗

Porting HEP Parameterized Calorimeter Simulation Code to GPUs

The High Energy Physics (HEP) experiments, such as those at the Large Hadron Collider (LHC), traditionally consume large amounts of CPU cycles for detector simulations and data analysis, but rarely use compute accelerators such as GPUs. As the LHC is upgraded to allow for higher luminosity, resulting in much higher data rates, purely relying on CPUs may not provide enough computing power to support the simulation and data analysis needs. As a proof of concept, we investigate the feasibility of porting a HEP parameterized calorimeter simulation code to GPUs. We have chosen to use FastCaloSim, the ATLAS fast parametrized calorimeter simulation. While FastCaloSim is sufficiently fast such that it does not impose a bottleneck in detector simulations overall, significant speed-ups in the processing of large samples can be achieved from GPU parallelization at both the particle (intra-event) and event levels; this is especially beneficial in conditions expected at the high-luminosity LHC, where extremely high per-event particle multiplicities will result from the many simultaneous proton-proton collisions. We report our experience with porting FastCaloSim to NVIDIA GPUs using CUDA. A preliminary Kokkos implementation of FastCaloSim for portability to other parallel architectures is also described.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Verification of Electromagnetic Physics Models for Parallel Computing Architectures in the GeantV Project

An intensive R&D and programming effort is required to accomplish new challenges posed by future experimental high-energy particle physics (HEP) programs. The GeantV project aims to narrow the gap between the performance of the existing HEP detector simulation software and the ideal performance achievable, exploiting latest advances in computing technology. The project has developed a particle detector simulation prototype capable of transporting in parallel particles in complex geometries exploiting instruction level microparallelism (SIMD and SIMT), task-level parallelism (multithreading) and high-level parallelism (MPI), leveraging both the multi-core and the many-core opportunities. We present preliminary verification results concerning the electromagnetic (EM) physics models developed for parallel computing architectures within the GeantV project. In order to exploit the potential of vectorization and accelerators and to make the physics model effectively parallelizable, advanced sampling techniques have been implemented and tested. In this paper we introduce a set of automated statistical tests in order to verify the vectorized models by checking their consistency with the corresponding Geant4 models and to validate them against experimental data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Exploring the Great Pyramid: Stand-Alone Monte Carlo Simulations for Detector Optimization

Cosmic-ray muon imaging has been used to non-destructively examine the Pyramids of Khufu and Khafre on the Giza Plateau; the EGP project will continue this line of research by undertaking a full tomographic scan of the former and, in doing so, will increase the sensitivity of the technique by upwards of two orders of magnitude. This paper details the development of stand-alone Monte Carlo simulations intended to aid in optimizing the EGP detector system. While the rectangular scintillator model considered here does provide sucient position resolution in a few specic test cases, it is not the optimal telescope design. An attempted detector modication to permit tagging of muons by their incident momenta is also discussed.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

GEANT4 code for DANCE with NEUANCE detector-response simulations

A new GEANT4 model for DANCE with NEUANCE has been developed to account for the change in the DANCE configuration when NEUANCE is installed. In this model, the energy resolution and the shape of the threshold of the individual DANCE detectors are used as an input to provide realistic representation of the measured γ-ray spectra. In addition, the distance from each detector to the target is defined in an input file as well as the physical presence of the a detector in the DANCE and NEUANCE arrays. The azimuthal angle of NEUANCE is also an external parameter to be tuned to represent the measurement. Simulated spectra are compared with measurements with standard γ-ray calibrated sources.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Simulating the Phonon Collection Efficiency in Kinetic Inductance Phonon-Mediated Detectors

Kinetic inductance phonon-mediated (KIPM) detectors use microwave kinetic inductance detectors (MKIDs) to read out phonon signals in the substrate. They are a promising class of detectors to be used in light dark matter (DM) searches due to their potential eV-scale sensitivity and native frequency-domain multiplexability. In order to improve upon the design of these detectors, simulations are needed to understand the effects of detector design on measurable physical parameters including the phonon collection efficiency, $\eta_{ph}$, which is defined as the ratio of phonon energy detected by the sensitive target element to the incident energy deposited in the substrate. This work simulates the phonon collection efficiency for two existing detector architectures to study phonon absorption at substrate-sensor interfaces and phonon loss to non-sensitive detector elements.

Dang, Stella Q.↗

Detector and Beamline Simulation for Next-Generation High Energy Physics Experiments

The success of high energy physics programs relies heavily on accurate detector simulations and beam interaction modeling. The increasingly complex detector geometries and beam dynamics require sophisticated techniques in order to meet the demands of current and future experiments. Common software tools used today are unable to fully utilize modern computational resources, while data-recording rates are often orders of magnitude larger than what can be produced via simulation. In this paper, we describe the state, current and future needs of high energy physics detector and beamline simulations and related challenges, and we propose a number of possible ways to address them.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Optimizing Geant4 Hadronic Models

Geant4, the leading detector simulation toolkit used in high energy physics, employs a set of physics models to simulate interactions of particles with matter across a wide range of energies. These models, especially the hadronic ones, rely largely on directly measured cross-sections and inclusive characteristics, and use physically motivated parameters. However, they generally aim to cover a broad range of possible simulation tasks and may not always be optimized for a particular process or a given material. The Geant4 collaboration recently made many parameters of the models accessible via a configuration interface. This opens a possibility to fit simulated distributions to the thin target experimental datasets and extract optimal values of the model parameters and the associated uncertainties. Such efforts are currently undertaken by the Geant4 collaboration with the goal of offering alternative sets of model parameters, also known as "tunes", for certain applications. The effort should subsequently lead to more accurate estimates of the systematic errors in physics measurements given the detector simulation role in performing the physics measurements. Results of the study are presented to illustrate how Geant4 model parameters can be optimized through applying fitting techniques, to improve the agreement between the Geant4 and the experimental data.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗