Search NASA⌕ Search

SEARCH · Search NASA

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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Calibration of the DUNE Far Detector Using Cosmic-ray Muon Events

The Deep Underground Neutrino Experiment (DUNE) aims to set new limits on parameters associated with neutrino oscillations, neutrino astrophysics, and beyond the Standard Model (SM) searches such as nucleon decay. DUNE will quantify the magnitude of CP violation in the lepton sector, and determine the neutrino mass ordering. These benefit highly from the large target mass and excellent imaging, tracking, and particle identification capabilities of Liquid Argon Time Projection Chambers (LArTPCs). Detector calibration is essential to make precise physics measurements. For instance, accurate energy reconstruction is necessary for measuring many of the aforementioned quantities with the precision required for discovering new physics and fully exploiting the capabilities of the detector. Cosmic muons are a freely available natural source of calorimetric data and can be used for calibrating various detector parameters. This thesis provides an analysis of simulated cosmic-ray muon events generated with the Muon Simulation Underground (MUSUN) generator in the DUNE horizontal drift (HD) far detector (FD). The study focuses on analysing the energy and angular distribution of various classes of muon events, as well as characterising the different particles produced by cosmic muon interactions. The analysis of π0 → 2γ events within the cosmic-ray muon sample is presented in this thesis with a detailed study of reconstructing electromagnetic showers. The π0 mass is reconstructed within the DUNE FD, yielding a value of (136 ± 7) MeV/c2. Additionally, the thesis introduces methods for dE/dx calibration using simulated and reconstructed muon tracks. A calibration constant Ccal = (5.469 ± 0.003) × 10−3 ADC × tick/e is obtained through a model-dependent calibration process, where 1 tick corresponds to 500 ns of sampling time of an ADC. Furthermore, a calibration technique is presented, demonstrating precise translation from dQ/dx to dE/dx. This calibration method is applied to stopping muons, charged pions, and protons in the DUNE FD, addressing the measurement of energy loss in the detector volume. These are important calibrations of the DUNE FD and will contribute to achieving the exciting physics goals of the experiment.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2011

Dijet invariant mass distribution in top quark hadronic decay with QCD corrections

The dijet invariant mass distributions from the hadronic decay of unpolarized top quark (t{yields}bW{sup +} followed by W{sup +}{yields}ud) are calculated, including the next-to-leading order QCD radiative corrections. We treat the top decay in the complex mass scheme due to the existence of the intermediate state W boson. Our analytical expressions are also available in different dimensional regularization schemes and {gamma}{sub 5} strategies. Finally, in order to construct the jets, we use different jet algorithms to compare their influences on our results. The obtained dijet mass distributions from the top quark decay are useful to distinguish these dijets from those produced via other sources and to clarify the issue about the recent CDF Collaboration's Wjj anomaly.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2010

Response of lithium formate EPR dosimeters at photon energies relevant to the dosimetry of brachytherapy

Purpose: To investigate experimentally the energy dependence of the detector response of lithium formate EPR dosimeters for photon energies below 1 MeV relative to that at {sup 60}Co energies. High energy photon beams are used in calibrating dosimeters for use in brachytherapy since the absorbed dose to water can be determined with high accuracy in such beams using calibrated ion chambers and standard dosimetry protocols. In addition to any differences in mass-energy absorption properties between water and detector, variations in radiation yield (detector response) with radiation quality, caused by differences in the density of ionization in the energy imparted (LET), may exist. Knowledge of an eventual deviation in detector response with photon energy is important for attaining high accuracy in measured brachytherapy dose distributions. Methods: Lithium formate EPR dosimeters were irradiated to known levels of air kerma in 25-250 kV x-ray beams and in {sup 137}Cs and {sup 60}Co beams at the Swedish Secondary Standards Dosimetry Laboratory. Conversions from air kerma free in air into values of mean absorbed dose to the detectors were made using EGSnrc MC simulations and x-ray energy spectra measured or calculated for the actual beams. The signals from the detectors were measured using EPR spectrometry. Detector response (the EPR signal per mean absorbed dose to the detector) relative to that for {sup 60}Co was determined for each beam quality. Results: Significant decreases in the relative response ranging from 5% to 6% were seen for x-ray beams at tube voltages {<=}180 kV. No significant reduction in the relative response was seen for {sup 137}Cs and 250 kV x rays. Conclusions: When calibrated in {sup 60}Co or MV photon beams, corrections for the photon energy dependence of detector response are needed to achieve the highest accuracy when using lithium formate EPR dosimeters for measuring absorbed doses around brachytherapy sources emitting photons in the energy range of 20-150 keV such as {sup 169}Yb and electronic sources.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2022

Novel approach for evaluating detector-related uncertainties in a LArTPC using MicroBooNE data

Primary challenges for current and future precision neutrino experiments using liquid argon time projection chambers (LArTPCs) include understanding detector effects and quantifying the associated systematic uncertainties. This paper presents a novel technique for assessing and propagating LArTPC detector-related systematic uncertainties. The technique makes modifications to simulation waveforms based on a parameterization of observed differences in ionization signals from the TPC between data and simulation, while remaining insensitive to the details of the detector model. The modifications are then used to quantify the systematic differences in low- and high-level reconstructed quantities. This approach could be applied to future LArTPC detectors, such as those used in SBN and DUNE.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Portable Programming Model Exploration for LArTPC Simulation in a Heterogeneous Computing Environment: OpenMP vs. SYCL

The evolution of the computing landscape has resulted in the proliferation of diverse hardware architectures, with different flavors of GPUs and other compute accelerators becoming more widely available. To facilitate the efficient use of these architectures in a heterogeneous computing environment, several programming models are available to enable portability and performance across different computing systems, such as Kokkos, SYCL, OpenMP and others. As part of the High Energy Physics Center for Computational Excellence (HEP-CCE) project, we investigate if and how these different programming models may be suitable for experimental HEP workflows through a few representative use cases. One of such use cases is the Liquid Argon Time Projection Chamber (LArTPC) simulation which is essential for LArTPC detector design, validation and data analysis. Following up on our previous investigations of using Kokkos to port LArTPC simulation in the Wire-Cell Toolkit (WCT) to GPUs, we have explored OpenMP and SYCL as potential portable programming models for WCT, with the goal to make diverse computing resources accessible to the LArTPC simulations. In this work, we describe how we utilize relevant features of OpenMP and SYCL for the LArTPC simulation module in WCT. We also show performance benchmark results on multi-core CPUs, NVIDIA and AMD GPUs for both the OpenMP and the SYCL implementations. Comparisons with different compilers will also be given where appropriate.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2011

Characterization of Monoenergetic Low Energy Neutron Fields with the {mu}TPC Detector

The AMANDE facility produces monoenergetic neutron fields from 2 keV to 20 MeV for metrological purposes. To be considered as a reference facility, fluence and energy distributions of neutron fields have to be determined by primary measurement standards. For this purpose, a micro Time Projection Chamber is being developed to be dedicated to measure neutron fields with energy ranging from 2 keV up to 1 MeV. We present simulations showing that such a detector, which allows the measurement of the ionization energy and the 3D reconstruction of the recoil nucleus, provides the determination of neutron energy and fluence of such low energy neutron fields.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2013

A novel method for modeling the neutron time of flight (nTOF) detector response in current mode to inertial confinement fusion experiments.

There are several machines in this country that produce short bursts of neutrons for various applications. A few examples are the Z machine, operated by Sandia National Laboratories in Albuquerque, NM ; the OMEGA Laser Facility at the University of Rochester in Rochester, NY ; and the National Ignition Facility (NIF) operated by the Department of Energy at Lawrence Livermore National Laboratory in Livermore, California . They all incorporate neutron time of flight (nTOF) detectors which measure neutron yield, and the shapes of the waveforms from these detectors contain germane information about the plasma conditions that produce the neutrons. However, the signals can also be %E2%80%9Cclouded%E2%80%9D by a certain fraction of neutrons that scatter off structural components and also arrive at the detectors, thereby making analysis of the plasma conditions more difficult. These detectors operate in current mode - i.e., they have no discrimination, and all the photomultiplier anode charges are integrated rather than counted individually as they are in single event counting. Up to now, there has not been a method for modeling an nTOF detector operating in current mode. MCNP PoliMi was developed in 2002 to simulate neutron and gamma ray detection in a plastic scintillator, which produces a collision data output table about each neutron and photon interaction occurring within the scintillator; however, the post processing code which accompanies MCNP PoliMi assumes a detector operating in single event counting mode and not current mode. Therefore, the idea for this work had been born: could a new post processing code be written to simulate an nTOF detector operating in current mode? And if so, could this process be used to address such issues as the impact of neutron scattering on the primary signal? Also, could it possibly even identify sources of scattering (i.e., structural materials) that could be removed or modified to produce %E2%80%9Ccleaner%E2%80%9D neutron signals? This process was first developed and then applied to the axial neutron time of flight detectors at the Z Facility mentioned above. First, MCNP PoliMi was used to model relevant portions of the facility between the source and the detector locations. To obtain useful statistics, variance reduction was utilized. Then, the resulting collision output table produced by MCNP PoliMi was further analyzed by a MATLAB postprocessing code. This converted the energy deposited by neutron and photon interactions in the plastic scintillator (i.e., nTOF detector) into light output, in units of MeVee%D1%84 (electron equivalent) vs time. The time response of the detector was then folded into the signal via another MATLAB code. The simulated response was then compared with experimental data and shown to be in good agreement. To address the issue of neutron scattering, an %E2%80%9CIdeal Case,%E2%80%9D (i.e., a plastic scintillator was placed at the same distance from the source for each detector location) with no structural components in the problem. This was done to produce as %E2%80%9Cpure%E2%80%9D a neutron signal as possible. The simulated waveform from this %E2%80%9CIdeal Case%E2%80%9D was then compared with the simulated data from the %E2%80%9CFull Scale%E2%80%9D geometry (i.e., the detector at the same location, but with all the structural materials now included). The %E2%80%9CIdeal Case%E2%80%9D was subtracted from the %E2%80%9CFull Scale%E2%80%9D geometry case, and this was determined to be the contribution due to scattering. The time response was deconvolved out of the empirical data, and the contribution due to scattering was then subtracted out of it. A transformation was then made from dN/dt to dN/dE to obtain neutron spectra at two different detector locations.

Nelson, Alan J.↗
DOE OSTI2018

Comparison between simulated and observed LHC beam backgrounds in the ATLAS experiment at E beam =4 TeV

Results of dedicated Monte Carlo simulations of beam-induced background (BIB) in the ATLAS experiment at the Large Hadron Collider (LHC) are presented and compared with data recorded in 2012. During normal physics operation this background arises mainly from scattering of the 4 TeV protons on residual gas in the beam pipe. Methods of reconstructing the BIB signals in the ATLAS detector, developed and implemented in the simulation chain based on the \textsc{Fluka} Monte Carlo simulation package, are described. The interaction rates are determined from the residual gas pressure distribution in the LHC ring in order to set an absolute scale on the predicted rates of BIB so that they can be compared quantitatively with data. Through these comparisons the origins of the BIB leading to different observables in the ATLAS detectors are analysed. The level of agreement between simulation results and BIB measurements by ATLAS in 2012 demonstrates that a good understanding of the origin of BIB has been reached.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

A new track finding algorithm based on a multi-dimensional extension of the Hough Transform

We introduce a new pattern recognition algorithm for track finding in High Energy Physics Experiments based on an extension of the Hough Transform to multiple dimensions. A remarkable property of this algorithm is that the execution time is simply proportional to the total number of the hits to be processed, making it particularly attractive for high occupancy situations. The algorithm needs to be trained using a sufficiently large set of simulated tracks. The same track finding algorithm can be used for very different detector geometries and only the set of simulated tracks used for training needs to be changed. The particular structure of the algorithm also lends itself naturally to parallel hardware implementations which, combined with its intrinsic flexibility, should provide a most powerful tool for triggering at future colliders.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2018

Lasnex simulated neutron production rates of MagLIF driven by 17.5 MA

Time resolved neutron detection will provide important information of MagLIF implosions. To aid the design of such detectors Lasnex has been used to simulate the neutron production rates for two MagLIF configurations. The results are shown in Fig 1 a) B=15 T, preheat=1.1 kJ, D 2 density=1.05 mg/cc and b) B=10 T, preheat=0.8 kJ, D 2 density=0.7 mg/cc.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1997

Simulation results of Pulse Shape Discrimination (PSD) for background reduction in INTEGRAL Spectrometer (SPI) germanium detectors

Pulse Shape Discrimination (PSD) for background reduction will be used in the INTErnational Gamma Ray Astrophysics Laboratory (INTEGRAL) imaging spectrometer (SPI) to improve the sensitivity from 200 keV to 2 MeV. The observation of significant astrophysical gamma ray lines in this energy range is expected, where the dominant component of the background is the beta(sup -) decay in the Ge detectors due to the activation of Ge nuclei by cosmic rays. The sensitivity of the SPI will be improved by rejecting beta(sup -) decay events while retaining photon events. The PSD technique will distinguish between single and multiple site events. Simulation results of PSD for INTEGRAL-type Ge detectors using a numerical model for pulse shape generation are presented. The model was shown to agree with the experimental results for a narrow inner bore closed end cylindrical detector. Using PSD, a sensitivity improvement factor of the order of 2.4 at 0.8 MeV is expected.

Slassi-Sennou, S. A.↗
DOE OSTI2018

Electromagnetic Simulations of Newly Designed Semiconductor Bolometers for Submillimeter Observations

In the frame of the SPICA mission, new semiconductor bolometers have been designed for the POL instrument which is directly part of the HERSCHEL/PACS legacy. These detectors are cooled down to 50 mK and are sensitive to polarization for the 100, 200 and 350 μm wavelength bands. In this paper, 3D simulations have been carried out with the COMSOL Multiphysics software to design the two suspended interlaced silicon spirals which support the orthogonal absorbers (sensitive to the two polarizations) of the POL detectors. These simulations allow to optimize the design of the absorbers and to quantify the cross-polarization. In the last part of this paper, we explain how it is possible to reach the 200 and 350 μm bands from the layout of the 100-μm detector.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Measurement of the muon neutrino charged-current mesonless cross section in the NOvA near detector

NOvA is a long-baseline accelerator neutrino experiment at Fermilab. Its physics goals include precision neutrino oscillation measurements, neutrino interaction cross-section measurements and beyond Standard Model explorations. We present a measurement of muon neutrino charged-current cross section with zero mesons in the final state at the NOvA near detector. This measurement is performed as a function of the kinematics of the final state muon. Our chosen interaction channel is especially sensitive to quasielastic and meson exchange current interactions and it provides handles for constraining the cross section systematic uncertainties in oscillation analyses in present and future experiments. For particle identification, we use a convolutional neural network (CNN) trained on individual particles simulated in the NOvA Near detector. This allows us to select the desired signal while reducing the potential bias from neutrino interaction modeling. We study strategies for constraining the remaining charged-pion background via Michel electron information in a template fitting approach. The main experimental result is a two-dimensional differential cross section as a function of final-state muon kinetic energy and polar angle. The parameters of this measurement, including binning and unfolding, were optimized to reduce the expected systematic uncertainty in the total cross section. The final result shows good agreement with the main GENIE-based simulation framework that was independently fine-tuned in NOvA. We finally propose improvements and subsequent steps that build on this analysis and further dissect the final states of neutrino interactions. This work has been supported by US DOE grant DE-SC0015684.

Sánchez Falero, Sebastián Jesús [Iowa State U.]↗
DOE OSTI2021

Novel features and GPU performance analysis for EM particle transport in the Celeritas code

Celeritas is a new computational transport code designed for high-performance simulation of high-energy physics detectors. This work describes some of its current capabilities and the design choices that enable the rapid development of efficient on-device physics. The abstractions that underpin the code design facilitate low-level performance tweaks that require no changes to the higher-level physics code. We evaluate a set of independent changes that together yield an almost 40% speedup over the original GPU code for a net performance increase of 220× for a single GPU over a single CPU running 8.4M tracks on a small demonstration physics app.

Johnson, Seth R.↗
DOE OSTI2018

Construction and performance of the barrel electromagnetic calorimeter for the Gluex experiment

The barrel calorimeter is part of the new spectrometer installed in Hall D at Jefferson Lab for the GlueX experiment. The calorimeter was installed in 2013, commissioned in 2014 and has been operating routinely since early 2015. The detector configuration, associated Monte Carlo simulations, calibration and operational performance are described in this paper. The calorimeter records the time and energy deposited by charged and neutral particles created by a multi-GeV photon beam. It is constructed as a lead and scintillating-fiber calorimeter and read out with 3840 large-area silicon photomultiplier arrays. Particles impinge on the detector over a wide range of angles, from normal incidence at 90 degrees down to 11.5 degrees, which defines a geometry that is fairly unique among calorimeters. The response of the calorimeter has been measured during a running experiment and performs as expected for electromagnetic showers below 2.5 GeV. Finally, we characterize the performance of the BCAL using the energy resolution integrated over typical angular distributions for π 0 and η production of σ E / E = 5 . 2 % / E ( GeV ) ⊕ 3 . 6 % and a timing resolution of σ = 150 ps at 1 GeV.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Modelling CN Zeeman effect observations of the envelopes of a low-mass protostellar disc and a massive protostar

We use the POLARIS radiative transfer code to produce simulated circular polarization Zeeman emission maps of the cyanide (CN) J = 1–0 molecular line transition for two types of protostellar envelope magnetohydrodynamic simulations. Our first model is a low-mass disc envelope system (box length L = 200 au), and our second model is the envelope of a massive protostar (L = 104 au) with a protostellar wind and a CN-enhanced outflow shell. We compute the velocity-integrated Stokes I and V, as well as the implied V/I polarization percentage, for each detector pixel location in our simulated emission maps. Our results show that both types of protostellar environments are in principle accessible with current circular polarization instruments, with each containing swaths of envelope area that yield percentage polarizations that exceed the 1.8 percent nominal sensitivity limit for circular polarization experiments with the Atacama Large Millimeter/submillimeter Array. In both systems, high-polarization (≳1.8 percent) pixels tend to lie at an intermediate distance away from the central star and where the line-centre opacity of the CN emission is moderately optically thin (τLC ~ 0.1–1). Furthermore, our computed V/I values scale roughly with the density-weighted mean line-of-sight magnetic field strength, indicating that Zeeman observations can effectively diagnose the strength of envelope-scale magnetic fields. We also find that pixels with large V/I are preferentially co-located where the absolute value of the velocity-integrated V is also large, suggesting that locations with favourable percentage polarization are also favourable in terms of raw signal.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Test beam study of SiPM-on-tile configurations

Light yield and spatial uniformity for a large variety of configurations of scintillator tiles were studied. The light from each scintillator was collected by a Silicon Photomultiplier (SiPM) directly viewing the produced scintillation light (SiPM-on-tile technique). The varied parameters included tile transverse size, tile thickness, tile wrapping material, scintillator composition, and SiPM model. Furthermore, these studies were performed using 120 GeV protons at the Fermilab Test Beam Facility. External tracking allowed the position of each proton penetrating a tile to be measured. The results were compared to a GEANT4 simulation of each configuration of scinitillator, wrapping, and SiPM.

47 OTHER INSTRUMENTATION↗
DOE OSTI2015

Muon Neutrino to Electron Neutrino Oscillation in NO$\nu$A

NOvA is a long-baseline neutrino oscillation experiment optimized for electron neu- trino ( e) appearance in the NuMI beam, a muon neutrino ( $\nu$ μ ) source at Fermilab. It consists of two functionally identical, nearly fully-active liquid-scintillator tracking calorimeters. The near detector (ND) at Fermilab is used to study the neutrino beam spectrum and composition before oscillation, and measure background rate to the $\nu$ e appearance search. The far detector, 810 km away in Northern Minnesota, observes the oscillated beam and is used to extract oscillation parameters from the data. NO$\nu$ A's long baseline, combined with the ability of the NuMI beam to operate in the anti- neutrino mode, makes NO$\nu$ A sensitive to the last unmeasured parameters in neutrino oscillations- mass hierarchy, CP violation and the octant of mixing angle 23. This thesis presents the search for e appearance in the first data collected by the NO A detectors from October 2013 till May 2015. Studies of the NuMI neutrino data collected in the NO A near detector are also presented, which show large discrepancies between the ND simulation and data. Muon- removed electron (MRE) events, constructed by replacing the muon in $\nu$ μ charged cur- rent interactions by a simulated electron, are used to correct the far detector $\nu$e appearance prediction for these discrepancies. In the analysis of the first data, a total of 6 e candidate events are observed in the far detector on a background of 1, a 3.46 σ excess, which is interpreted as strong evidence for e appearance. The results are consistent with our expectation, based on constraints from other neutrino oscillation experiments. The result presented here differs from the officially published $\nu$ e appearance result from the NO A experiment where the systematic error is assumed to cover the MRE correction.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗