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 991 records · Page 55

A hierarchical Bayesian model for background variation in radiation source localization

In this paper, we apply a new model to account for varying background radiation in radiological source localization. We present a hierarchical Bayesian model that simultaneously infers background and source location parameters without requiring separate estimation of the background radiation at each detector location. We employ a simplified photon transport model to reduce the computational expense of Bayesian model calibration. Furthermore, we demonstrate the model accuracy by localizing a cesium-137 source in a simulated city block, and we analyze experimental field measurements with varying background. In both cases, the model provides sufficient fidelity that we can locate the source while simultaneously estimating background radiation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Low-energy 23 Al β-delayed proton decay and 22 Na destruction in novae

The radionuclide 22 Na is a target of γ-ray astronomy searches, predicted to be produced during thermonuclear runaways driving classical novae. The 22 Na (p, γ) 23 Mg reaction is the main destruction channel of 22 Na during a nova, hence, its rate is needed to accurately predict the 22 Na yield. However, experimental determinations of the resonance strengths have led to inconsistent results. In this Rapid Communication, we report a measurement of the branching ratios of the 23 Al β-delayed protons as a probe of the key 204-keV (center-of-mass) 22 Na (p, γ) 23 Mg resonance strength. We report a factor of 5 lower branching ratio compared to the most recent literature value. The variation in 22 Na yield due to nuclear data inconsistencies was assessed using a series of hydrodynamic nova outburst simulations and has increased to a factor of 3.8, corresponding to a factor of ≈ 2 uncertainty in the maximum detectability distance. Finally, this is the first reported scientific measurement using the Gaseous Detector with Germanium Tagging system.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Beam Dynamics of the Muon $g\textrm{-}2$ Experiment

The Muon $g\textrm{-}2$ Experiment (E989) at Fermilab aims to measure the muon anomalous magnetic moment $a_{\mu}$ with unprecedented precision, potentially uncovering physics beyond the Standard Model of particle physics. The result based on Runs 1-3, released in 2023, achieved a precision of 0.20 ppm. The experiment circulates muons in a storage ring, measuring $a_{\mu}$ from decay positron time and energy measurements collected with calorimeters. To achieve the required accuracy, it is crucial to measure and control the magnetic field in the ring with high precision. Beam dynamics corrections are necessary for muons not orbiting exactly in the midplane, for their oscillations, and for electric field effects. Highly accurate beam dynamics simulations are instrumental for quantifying and validating the beam dynamics corrections, ultimately improving the precision of the $a_{\mu}$ measurement and facilitating the achievement of the ambitious $70\:\mathrm{ppb}$ systematic uncertainty goal. The measured field data was incorporated into models for simulations using three codes: \texttt{gm2ringsim} (an internal Geant4-based code), \textit{COSY INFINITY}, and \textit{BMAD}. The advantages of \texttt{gm2ringsim} include using CAD-based geometry and modelling the detector effects. \textit{COSY INFINITY} is a highly accurate and efficient code that uses high-order differential-algebraic transfer maps, precise fringe field calculations, and advanced symplectification methods. Symplectification is important for maintaining the physical correctness of the muon beam behaviour with high precision over the storage time, ensuring conservation of phase space volume and preventing artificial damping or excitation of particle motion. The experiment completed its final Run 6 in July 2023, collecting 21 times more data than the previous BNL experiment. Analyses of data from Runs 4-6 are ongoing, with results planned for release in 2025, potentially resolving the current tension between experiment and theory.

43 PARTICLE ACCELERATORS↗

Jamming Detection for Low-Resolution SC-FDE Systems: A Machine Learning Approach

Jammers interfere with communication between base stations (BSs) and legitimate users, leading to degradation of wireless system performance. Our study focuses on jamming detection for wideband single-carrier frequency domain equalization (SC-FDE) systems with low-resolution analog-to digital converters (ADCs). In such systems, jamming detection is challenging because traditional analytical approaches cannot be directly applied due to the delay dispersion in wideband channels and the non-linearity induced by low-resolution ADCs. We propose a machine learning (ML)-based jamming detection method that directly uses the quantized receive signals. Significantly, our ML-based detector can be integrated into existing standard frameworks, such as unique word (UW)-based SC-FDE systems, as it uses existing pilots without requiring additional pilots for jamming detection. Through numerical simulations, we show that two or more bits provide satisfactory performance compared to unquantized scenarios. Additionally, we demonstrate that using more and well-separated pilot symbols improves performance.

99 GENERAL AND MISCELLANEOUS↗

Understanding hadronic interactions using advanced reconstruction techniques

ProtoDUNE-SP served as the prototype of the future Deep Underground Neutrino Experiment (DUNE). With a total liquid argon mass of 0.77 kt, it stood as the most extensive monolithic single-phase Liquid Argon Time Projection Chamber (LArTPC) ever constructed. Strategically located and operated at CERN, ProtoDUNE-SP benefited from a specialized charged-particle test beam, offering a momentum range of 0.3 - 7 GeV/c. This setup provided a unique opportunity to conduct in-depth studies on hadronic interactions within argon. Hadronic interactions play a fundamental role in neutrino physics, particularly in the intricacies of neutrino detection. We depend on analyzing the outcomes of these interactions to accurately determine both the flavor and energy of the interacting neutrinos. Pions, as one of the primary byproducts of neutrino interactions, are crucial for understanding the overall dynamics and kinematics of these processes. Accurately characterizing and understanding the interaction of pions with argon can significantly enhance the precision of neutrino simulations and measurements.Being a Liquid Argon Time Projection Chamber (LArTPC), ProtoDUNE-SP is distinguished by its ability to produce high-definition images of charged particles as they traverse through the detector's active volume. Yet, accurately reconstructing these particle interactions and determining their kinematic attributes remains a challenging task. The successes of deep learning in diverse domains, especially in image recognition, provide a promising approach for addressing this challenge in ProtoDUNE. This thesis discusses a measurement of the Piplus - argon inelastic cross section in the energy range of 2400 to 3000 MeV, using data taken by protoDUNE during the fall of 2018. Additionally, panoptic segmentation, a machine learning technique, is introduced and validated through its application in the reconstruction of the neutral pion rest mass. The aim is to showcase how these advanced methods can enhance the quality of event reconstruction.

Sarasty Segura, Carlos Eduardo↗

Detectability of neutrino-signal fluctuations induced by the hadron-quark phase transition in failing core-collapse supernovae

Here we introduce a systematic and quantitative methodology for establishing the presence of neutrino oscillatory signals due to the hadron-quark phase transition (PT) in failing core-collapse supernovae from the observed neutrino event rate in water- or ice-based neutrino detectors. The methodology uses a likelihood ratio in the frequency domain as a test-statistic; it is employed for quantitative analysis of neutrino signals without assuming the frequency, amplitude, starting time, and duration of the PT-induced oscillations present in the neutrino events and thus it is suitable for analyzing neutrino signals from a wide variety of numerical simulations. We test the validity of this method by using a core-collapse simulation of a 17 solar-mass star by Zha et al. [Astrophys. J. 911, 74 (2021) ]. Based on this model, we further report the presence of a PT-induced oscillations quantitatively for a core-collapse supernovae out to a distance of ~10 kpc, ~5 kpc for IceCube and to a distance of ~10 kpc, ~5 kpc, and ~1 kpc for a 0.4 Mt mass water Cherenkov detector. This methodology will aid the investigation of a future galactic supernova and the study of hadron-quark phase in the core of core-collapse supernovae.

79 ASTRONOMY AND ASTROPHYSICS↗

Intelligent Experiments through Real-Time AI: Fast Data Processing and Autonomous Detector Control for High-Energy Nuclear Experiments

The aim of this project is to develop software and hardware for fast real-time data processing and autonomous detector control and calibration for the sPHENIX and the future EIC experiments. Below summarizes Georgia Tech team efforts in the past year: 1. We developed a real-time clustering algorithm and FPGA-based pipeline architecture for processing fired pixel data from ALPIDE sensors in sPHENIX experiments. Our Columnar Clustering Co-Design introduces a hardware-aware, stream-friendly approach that segments pixel data by column pairs using a Column Pair Clustering (CPC) strategy, followed by Cluster Stitching to merge adjacent subclusters. Implemented in Vitis HLS, the pipeline comprises five stages—read-in, subclustering, stitching, analysis, and write-out—connected by tagged HLS streams with custom end-of-event signaling for robust synchronization. We designed a pipelined dataflow model optimized for throughput, low latency, and minimal buffering, enabling scalable clustering across events of arbitrary size. Our system maintains spatial precision via center-of-mass and shape key extraction and efficiently handles edge cases such as fragmented or nested clusters. Compared against DBSCAN in both software and hardware, our approach demonstrates competitive performance under FPGA constraints. 2. We also conducted a comprehensive algorithm-to-hardware co-design of connected component analysis tailored for sPHENIX experiments, focusing on real-time, low-latency processing using FPGAs and High-Level Synthesis (HLS). Starting from a Python-based particle tracking pipeline, the team translated the core logic—graph traversal via DFS and Union-Find—into an HLS-compatible C++ model, replacing dynamic memory and recursion with static arrays and pipelined control flow. The final design includes a fully streamed and dataflow-compatible Union-Find kernel optimized across five iterations, incorporating loop pipelining, array partitioning, AXI/FIFO interface tuning, and function flattening. Experimental results show up to 14.8× speedup over the CPU baseline, reducing per-graph latency to 1.58 μs and demonstrating strong resource efficiency with only ~7k LUTs and zero BRAM usage. The design maintains functional correctness against the Python reference using a Python-based C-simulation framework and Mean Squared Error metrics. This work validates the potential of HLS-driven FPGA designs for edge-level HEP data acquisition, laying a scalable foundation for future integration with real-time detector pipelines and multi-graph processing systems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Development of an MCP-Based Timing Layer for the LHCb ECAL Upgrade-2

The increase in instantaneous luminosity during the high-luminosity phase of the LHC represents a significant challenge for future detectors. A strategy to cope with high-pileup conditions is to add a fourth dimension to the measurements of the hits, by exploiting the time separation of the various proton–proton primary collisions. According to LHCb simulation studies, resolutions of about 10–20 picoseconds, at least an order of magnitude shorter than the average time span between primary interactions, would be greatly beneficial for the physics reach of the experiment. Microchannel plate (MCP) photomultipliers are compact devices capable of measuring the arrival time of charged particles with the required resolution. The technology of large-area picosecond photodetectors (LAPPDs) is under investigation to implement a timing layer that can be placed within a sampling calorimeter module with the purpose of measuring the arrival time of electromagnetic showers. LAPPD performances, using a Gen-I tile with a delay-line anode and a Gen-II with a capacitively coupled anode, have been measured thoroughly both with laser (wavelength of 405 nm and pulse width of 27.5 ps FWHM) and high-energy electron (1–5.8 GeV) beams. Time resolutions of the order of 30 ps for single photoelectrons and 15 ps for electromagnetic showers initiated by 5-GeV electrons, as measured at the shower maximum, are obtained.

47 OTHER INSTRUMENTATION↗

Boosted decision tree reweighting of simulated neutrino interactions for O ( 1 ) GeV neutrino cross-section measurements

This paper illustrates a generic method for multidimensional reweighting of O ( 1 ) GeV neutrino interaction Monte Carlo samples. The reweighting is based on a boosted decision tree algorithm trained on high-dimensional space in detector final-state observables. This enables one generator’s events to be reweighted so that its reconstructed particle content and kinematics distributions, as well as detector efficiency, match those of a target model. The approach establishes an efficient way to reuse legacy Monte Carlo data, avoiding regeneration. As an example, we test its use in a measurement of transverse kinematic imbalance of the μ - and proton in charged-current quasielastic like ν μ events from the MINERvA experiment.

Lin, Z. [Rochester U.] (ORCID:0009000188903698)↗

Geant4 Monte-Carlo (GEMC) A database-driven simulation program

GEMC[1] is an application that harnesses the power of databases to execute Geant4 Monte-Carlo simulations. The databases (MYSQL, CSQL, TEXT) define the geometry, materials, digitization algorithms, readout electronics and output formats. Implemented in C++, GEMC also boasts a user-friendly Python API that facilitates detector construction and database population. GEMC can handle real-life scenarios such as geometry variations and the run number-dependent calibration constants and digitization parameters. This abstract provides an overview of GEMC, accompanied by examples that showcase its versatility. We delve into the practical application of GEMC within the the CLAS12 experimental program at Jefferson Lab.

Ungaro, Maurizio↗

Data-driven Alignment of the Tracker for the Mu2e Experiment

An alignment of the Mu2e Tracker is investigated using the Mu2e Offline Software and Millepede-II. Simulated misalignments of Tracker planes were inserted between Monte-Carlo event generation and before reconstruction using randomly generated $(x,y,z)$-translations and rotations of the rigid-body planes in the detector model. Cosmic ray muons, with no magnetic field present in the Detector Solenoid, will follow a straight-line trajectory. A maximum-likelihood fit is used to reconstruct a track and determine its parameters. Tracks traversing at least four planes, with ten registered straw hits, and a maximum time residual of $20$ ns are selected to be used in an alignment sample. The cuts were chosen to discard tracks with low information content, such as those with high levels of multiple scattering or those with low numbers of measurements contributing to a single track fit. Millepede-II produces the resulting alignment constants according to a simultaneous least-squares fit of the alignment parameters and the track parameters for all input tracks. The improved alignment constants are used to repeat the process and recover an increased number of suitable tracks, further improving the estimate of the alignment constants in subsequent iterations. A performance study of the process was conducted on a misaligned Tracker in the case of translation-only, and translation plus rotational misalignments. In both cases, the correct alignment constants were recovered to within statistical error after three global fit iterations. Additionally, a continuous integration system was developed for the \texttt{Mu2e/Offline} GitHub repository, introducing a convenient automated testing process for all proposed changes (``pull requests"). The tests developed included a full build of the software, followed by tests of the validation routines, geometry overlap checking, and multi-threading mode. Additionally, code quality checks and low-statistics physics comparison checks a re also available.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Impact of Time Dependent Reactor and Sensor Physics on Core Power Synthesis

Online synthesis of the power distribution is critical in the operation and control of nuclear power reactors to ensure that the core is operating within safety margins, and to provide essential knowledge associated with the burnup of the fuel. In light water reactors (LWRs), power synthesis is achieved by using some a priori knowledge of the state of the reactor core and updating based on the signals coming from in-core sensors—namely, self-powered neutron detectors (SPNDs). This report aims to study the effects of fuel burnup and sensor degradation on the ability to accurately synthesize the power distribution in a LWR. Several modeling tools were used to simulate power synthesis based on the responses of SPNDs, with emitters made out of Rh or V. A representative pressurized water reactor low-enriched uranium (LEU) core was modeled using the Polaris/Purdue Advanced Reactor Core Simulator (PARCS) approach. The Monte Carlo N-Particle Transport 6 (MCNP6) code was used, as well, to calculate response functions between different segments of fuel to individual SPNDs; this is a crucial parameter for power synthesis. The Oak Ridge Isotope GENeration (ORIGEN) package in the Standardized Computer Analyses for Licensing Evaluation (SCALE) code was used to model the time-dependent isotopic transmutation in the SPND emitters. All these data were fed into a custom code that enacted the point-based iterative (PBI) method to simulate power synthesis. Developmental work was also performed on high-fidelity SPND models in the GEometry ANd Tracking 4 (Geant4) code, which enables higher-accuracy modeling of the current responses from SPNDs. In this work, five sets of time-dependent power synthesis test cases were conducted. In these test cases, systematic changes in the input conditions enabled an analysis of the effect of (1) slightly inaccurate a priori power distribution assumptions with respect to fuel burnup, (2) highly inaccurate a priori assumptions with respect to fuel burnup (such that burnup is not included in the a priori assumed distribution), and (3) differences between Rh and V SPNDs in terms of downstream consequences of the transmutation in the emitters. The authors discovered that one may permissibly have slightly inaccurate a priori assumptions of the fuel burnup (such that the level of burnup may be slightly under- or over-approximated by the accumulated burnup in approximately 9.3 full power days), but to not account for burnup at all in the a priori assumption leads to severe levels of error, approaching 25% at maximum. The authors also discovered that V SPNDs are extraordinarily robust in the low-enriched uranium fuel cycle considered in this modeling work, whereas Rh SPNDs undergo significant transmutation that can result in large errors in the synthesized power distribution.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Fermilab's controls development with virtual accelerator

Control Systems development is often the last thing considered when designing and building new equipment, e.g. a new detector or superconducting RF LINAC; however when the new equipment is installed, it is the first thing desired to be operational for testing. Due to frequent delays in building new equipment and project deadlines, control system development and testing is often curtailed. A way to alleviate this problem is to simulate the control system, though this will be challenging for complex systems.The Fermilab PIP-II (proton improvement plan - II) project is being constructed at Fermilab to deliver $800\,MeV$ protons of $>1\,MW$ beam power to replace the present LINAC for the remainder of the existing accelerator complex. The new LINAC consists of a warm front end (WFE), 23 superconducting RF cryomodules (of 5 types), and a beam transfer line (BTL) to the existing complex.The accelerator physics group has a parallel project to create a digital twin (DT) of the PIP-II accelerator. We have coupled the EPICS controls to this DT and are developing both the DT and EPICS software in parallel. This will allow us to develop the EPICS software framework, the HMIs, sequences, high level physics applications, and other services for use in a fully functional control system.This presentation will detail the work that we have performed to date and show demonstrations of controlling and monitoring the status of the accelerator, as well as future plans for this work.

Hanlet, Pierrick [Fermilab]↗

Detectability of standing accretion shock instabilities activity in supernova neutrino signals

We introduce a novel methodology for establishing the presence of standing accretion shock instabilities (SASI) in the dynamics of a core collapse supernova from the observed neutrino event rate at water- or ice- based neutrino detectors. The methodology uses a likelihood ratio in the frequency domain as a test- statistics; it is also employed to assess the potential to estimate the frequency and the amplitude of the SASI modulations of the neutrino signal. The parameter estimation errors are consistent with the minimum possible errors as evaluated from the inverse of the Fisher information matrix, and close to the theoretical minimum for the SASI amplitude. Using results from a core-collapse simulation of a 15 solar-mass star by Kuroda et al. [Astrophys. J. 851, 62 (2017).] as a test bed for the method, we find that SASI can be identified with high confidence for a distance to the supernova of up to ~6 kpc for IceCube and up to ~3 kpc for a 0.4 Mt mass water Cherenkov detector. This methodology will aid the investigation of a future galactic supernova.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Compressive sensing optical design and simulation tool

A system or method for designing optical systems includes generating a mathematical model of an imaging system; tracing a plurality of light rays from a source through the model of the imaging system in a forward direction and a backward direction incident upon an arbitrarily defined detector location; calculating the total light rays incident upon the arbitrarily defined detector location in response to the plurality of light rays; tracing a series of calibration images first through the computational or compressive imager mathematical model; measuring an output set of parameters based on a plurality of detected rays calibration images at the defined detector location; storing the output set of parameters in software, and determining an estimated system transfer function.

42 ENGINEERING↗

Measurements of the ionization efficiency of protons in methane

The amount of energy released by a nuclear recoil ionizing the atoms of the active volume of detection appears “quenched” compared to an electron of the same kinetic energy. This different behavior in ionization between electrons and nuclei is described by the Ionization Quenching Factor (IQF) and it plays a crucial role in direct dark matter searches. For low kinetic energies (below 50 keV ), IQF measurements deviate significantly from common models used for theoretical predictions and simulations. We report measurements of the IQF for proton, an appropriate target for searches of Dark Matter candidates with a mass of approximately 1 GeV , with kinetic energies in between 2 keV and 13 keV in 100 mbar of methane. We used the Comimac facility in order to produce the motion of nuclei and electrons of controlled kinetic energy in the active volume, and a NEWS-G SPC to measure the deposited energy. The Comimac electrons are used as a reference to calibrate the detector with 7 energy points. A detailed study of systematic effects led to the final results well fitted by IQF (E K )=E$^α_K$ / (β+E$^α_K$) with α=0.70±0.08 and β=1.32±0.17 . In agreement with some previous works in other gas mixtures, we measured less ionization energy than predicted from SRIM simulations, the difference reaching 33% at 2 keV .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Adiabatic light guide with S-shaped strips

A light guide is an essential part of many scintillator counters and light collection systems. Our main interest is a light guide for a thin, wide scintillator that has high light transmission while converting the area of the light source to the shape of a photo-detector. Here, we propose a variation of the light guide that avoids a 90° twist of the strips, reduces the length of the light pipe, and reduces the complexity of production. Detailed Monte Carlo simulation studies have been performed for a three-strip S-shaped light-guide system.

47 OTHER INSTRUMENTATION↗

Search for the Anomalous Events Detected by ANITA Using the Pierre Auger Observatory

A dedicated search for upward-going air showers at zenith angles exceeding 110° and energies E > 0.1 EeV has been performed using the Fluorescence Detector of the Pierre Auger Observatory. The search is motivated by two “anomalous” radio pulses observed by the ANITA flights I and III that appear inconsistent with the standard model of particle physics. Using simulations of both regular cosmic-ray showers and upward-going events, a selection procedure has been defined to separate potential upward-going candidate events and the corresponding exposure has been calculated in the energy range [0.1–33] EeV. One event has been found in the search period between January 1, 2004, and December 31, 2018, consistent with an expected background of 0.27 ± 0.12 events from misreconstructed cosmic-ray showers. This translates to an upper bound on the integral flux of ( 7.2 ± 0.2 ) × 10 − 21 cm − 2 sr − 1 y − 1 and ( 3.6 ± 0.2 ) × 10 − 20 cm − 2 sr − 1 y − 1 for an E − 1 and E − 2 spectrum, respectively. An upward-going flux of showers normalized to the ANITA observations is shown to predict over 34 events for an E − 3 spectrum and over 8.1 events for a conservative E − 5 spectrum, in strong disagreement with the interpretation of the anomalous events as upward-going showers. Published by the American Physical Society 2025

Abdul Halim, A.↗