Search NASA⌕ Search

SEARCH · Search NASA

Results for “Particle tracking”

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 145 records · Page 8

Real-time charged track reconstruction for CLAS12

Abstract This paper presents the results of charged particle track reconstruction in CLAS12 using artificial intelligence. In our approach, we use machine learning algorithms to reconstruct tracks, including their momentum and direction, with high accuracy from raw hits of the CLAS12 drift chambers. The reconstruction is performed in real-time, with the rate of data acquisition, and allows for the identification of event topologies in real-time. This approach revolutionizes the Nuclear Physics experiments' data processing, allowing us to identify and categorize the experimental data on the fly, and will lead to a significant reduction in experiment data processing. It can also be used in streaming readout applications leading to more efficient data acquisition and post-processing.

Instruments & Instrumentation↗

The Influence of Ocean Topography on the Upwelling of Carbon in the Southern Ocean

The physical circulation of the Southern Ocean sets the surface concentration and thus air-sea exchange of CO 2 . However, we have a limited understanding of the three-dimensional circulation that brings deep carbon-rich waters to the surface. Here, we introduce and analyze a novel high-resolution ocean model simulation with active biogeochemistry and online Lagrangian particle tracking. We focus our attention on a subset of particles with high dissolved inorganic carbon (DIC) that originate below 1,000 m and eventually upwell into the near-surface layer (upper 200 m). We find that 71% of the DIC-enriched water upwelling across 1,000 m is concentrated near topographic features, which occupy just 33% of the Antarctic Circumpolar Current. Once particles upwell to the near-surface layer, they exhibit relatively uniform pCO 2 levels and DIC decorrelation timescales, regardless of their origin. Furthermore, our results show that Southern Ocean bathymetry plays a key role in delivering carbon-rich waters to the surface.

54 ENVIRONMENTAL SCIENCES↗

Test Beam Results of Planar Pixel Sensor for the CMS Phase 2 Inner Tracker Upgrade

Results of the test beam measurements that characterise the performance of CMS Readout Chip (CROC) sensors to be used in the High Luminosity era of the Large Hadron Collider (HL-LHC) are presented. The HL-LHC peak instantaneous luminosity of $7.5 \times 10^{34} \ \text{cm}^{-2} s^{-1}$ corresponds to an average of around 200 inelastic proton-proton collisions per beam-crossing every 25 ns. In order to efficiently reconstruct and track particles in this extreme and challenging conditions, the present CMS tracking detector will be completely replaced. The new tracking detector consists of an Inner Tracker closest to the beamline and an Outer Tracker surrounding it. These are populated with modules that comprise of readout chips and silicon sensors. The test beam measurements of these modules are vital to understand the performance of the related technologies. Using a primary 120 GeV proton beam from the Main Injector at Fermilab, data was collected at the Fermilab Test Beam Facility (FTBF) using the silicon tracker telescope that provides a precision position measurement of the track impact point with less than 5 $\mu$m uncertainty. The proton beam was incident on a 1x2 planar CROC module developed by Hamamatsu. The sensor has 100x25 $\mu m^2$ standard pixels and also a smaller number of 225x25 $\mu m^2$ longer pixels at the boundary between the two ROCs. We present characterization of these modules that includes pixel efficiency, resolution, cluster size and charge distributions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Modeling betatron radiation using particle-in-cell codes for plasma wakefield accelerator diagnostics

The analysis of plasma wakefield acceleration experimental measurements, particularly in the characterization of photons emitted through the betatron radiation mechanism, requires the development of accurate numerical models. These computational models are crucial for supporting modern instrumentation designed to measure the single-shot, double-differential angular-energy radiation spectra resulting from interactions between beams and plasmas. Motivated by the needs of such applications, this paper presents detailed numerical models of betatron radiation generated in beam-plasma acceleration experiments. These models are based on the integration of the Liénard-Wiechert (LW) potentials, applied to computed particle trajectories. The particle trajectories are generated using three distinct methods: first, by tracking particles through idealized fields in the blowout regime of PWFA; second, by obtaining trajectories using the fast quasistatic particle-in-cell (PIC) code quickpic; and third, obtaining trajectories from the fully self-consistent PIC code osiris. To ensure the accuracy and reliability of these models, the paper includes various benchmark tests using analytical expressions, as well as employing the PIC code epoch, which takes an alternative approach by using a Monte Carlo quantum electrodynamics (QED)-based radiation model. Additionally, the paper presents simulations of the expected experimental betatron radiation spectra, taking into account parameters relevant to PWFA and plasma photocathode experiments at the SLAC FACET-II facility.

Yadav, M. [University of California, Los Angeles, ↗

Processing MCNP Elemental Edit Outputs

The Los Alamos National Laboratory's (LANL) Monte Carlo N-Particle® (MCNP®) transport code version 6 has the capability for tracking particles on unstructured mesh (UM) geometry models. The MCNP UM feature was originally designed for the UM models created by the Abaqus/CAE software product. The MCNP version 6.2.0 and later recognizes the UM models read from Abaqus INP files or MCNPUM files converted from Abaqus INP files. Other finite element analysis software packages may generate UM models and then convert these models into Abaqus INP formats.

97 MATHEMATICS AND COMPUTING↗

Modelling Beam Loss Within the LANSCE Proton Storage Ring

Several upgrades are being considered for the proton storage ring (PSR) at the Los Alamos Neutron Science Center (LANSCE) to reduce beam loss and thereby reduce the cooldown period of the PSR following a beam run. First, we have considered an increased beam pipe diameter would reduce beam loss due to beam scraping caused by misalignments and tuning errors. However, this would require increased pole-to-pole gap height within the dipoles and quadrupoles, which would change their effective length and alter their fringe fields. The effect of different magnet gaps on the beam optical parameters and on beam loss was studied using the simulation codes MAD-X and PyORBIT. Second, we are developing a detailed particle tracking model within the framework of the simulation package General Particle Tracer (GPT) for the PSR H - stripping system. This model will be used to study the effect of the stripper foil parameters (position, composition, areal density, depth, etc.) on first-turn losses, where most of the observed beam loss and emittance growth occurs due to foil scattering, foil stripping, and Lorentz stripping. The model implements C++ custom elements created to model foil scattering, foil stripping, and Lorentz stripping, as these are not built-in features of GPT. The preliminary results of the MAD-X and PyORBIT simulations, as well as the GPT simulation model, is described in this tech note.

43 PARTICLE ACCELERATORS↗

Development of a pattern recognition algorithm for reconstructing multi-particle events in the Mu2e detector

Mu2e is an upcoming experiment at Fermilab and its main goal is to search for the Charged Lepton Flavor Violation (CLFV) in the coherent transition of a muon into an electron on an Al target. In Mu2e, multi-particle events can occur simultaneously within the same time region and it is crucial to accurately identify each particle track, including signals, to improve the robustness of track finding methods and enhance reconstruction efficiency. ¯p annihilation is one of the background events and produces multiple particles that can mimic signal events. Additionally, photons from radiative pion captures can produce a γ → $e+e−$ pair, which can be used to calibrate the Mu2e momentum scale and the resolution. The Mu2e track reconstruction sequence begins by grouping hits produced in the tracker based on time and z coordinate information, called TimeCluster, and selected hits are processed to reconstruct helices and determine their momentum. The current pattern recognition algorithms identify a single helix per TimeCluster for single track events. A new pattern recognition algorithm is being developed to reconstruct multi-particle events and its features for finding multiple tracks and the current evaluation results are reported.

Kitagawa, H. [Pisa U.]↗

Energy Reconstruction and Calibration of the MicroBooNE LArTPC

The Liquid Argon Time Projection Chamber (LArTPC) is increasingly becoming the chosen technology for current and future precision neutrino oscillation experiments due to its superior capability in particle tracking and energy calorimetry. In LArTPCs, calorimetric information is critical for particle identification, which is the foundation for the neutrino cross-section and oscillation measurements as well as searches for beyond standard model physics. One of the primary challenges in employing LArTPC technology is characterizing its performance and quantifying the associated systematic uncertainties. MicroBooNE, the longest operating LArTPC to date, has performed numerous such measurements, including studies of detector physics and electromagnetic shower reconstruction. Here we present results on the operation and performance of the detector during its data taking, highlighting accomplishments toward calorimetric reconstruction, calibration, and detector physics.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Jet energy scale and resolution measured in proton–proton collisions at $\sqrt{s}=13$ TeV with the ATLAS detector

Jet energy scale and resolution measurements with their associated uncertainties are reported for jets using 36–81 fb -1 of proton–proton collision data with a centre-of-mass energy of √s=13 TeV collected by the ATLAS detector at the LHC. Jets are reconstructed using two different input types: topo-clusters formed from energy deposits in calorimeter cells, as well as an algorithmic combination of charged-particle tracks with those topo-clusters, referred to as the ATLAS particle-flow reconstruction method. The anti-k t jet algorithm with radius parameter R=0.4 is the primary jet definition used for both jet types. This result presents new jet energy scale and resolution measurements in the high pile-up conditions of late LHC Run 2 as well as a full calibration of particle-flow jets in ATLAS. Jets are initially calibrated using a sequence of simulation-based corrections. Next, several in situ techniques are employed to correct for differences between data and simulation and to measure the resolution of jets. The systematic uncertainties in the jet energy scale for central jets (|η|<1.2) vary from 1% for a wide range of high-p T jets (250 T <2000 GeV), to 5% at very low p T (20 GeV) and 3.5% at very high p T (>2.5 TeV). The relative jet energy resolution is measured and ranges from (24±1.5)% at 20 GeV to (6±0.5)% at 300 GeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Improvement of the NOvA Near Detector Event Reconstruction and Primary Vertexing through the Application of Machine Learning Methods

The purpose of this work is to examine the application of a deep learning model in event reconstruction of neutrino interactions. The challenges faced in event reconstruction include the placement of an accurate primary neutrino interaction vertex which is used to support the particle track and prong algorithms. The result of accurate primary vertex ensures all particles involved in a neutrino interaction are included. We propose a regression-based Convolutional Neural Network (CNN) method to predict the primary vertex of a particle interaction. We show that with raw two-dimensional pixel map views as input, the regression-based CNN can predict the primary vertex in all three coordinates. This work is applied as part of the NOvA (NuMI Off-axis $\nu_e$ Appearance) near detector reconstruction efforts. The primary vertex predicted by the regression-based CNN model shows promising results for future applications. This deep learning method can be extended to secondary vertexing through a Kernel Density Estimate algorithm discussed in this work.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Quantum Enhanced Tracker

We report the study and development of a novel new charged particle detector, a "Quantum Enhanced Tracker" for 3D imaging of charged particle tracks. The proposed detector relies on, and benefits from, the extremely high sensitivity of atoms, prepared in a specific quantum superposition, to external perturbations. Atoms are prepared in a so-called "dark superposition" state, a passing charged particle will disturb this quantum state, allowing atoms to change polarization of the probe light or absorb probe light and to fluoresce, enabling direct imaging of the particle trajectory with high resolution. This method has similarities with Cherenkov detectors where a charged particle produces light along its track, but with the major advantages of enhanced sensitivity due to the low energy required to perturb the engineered quantum state, the ability to amplify the signal using a drive laser passing through the medium, and the ability to provide a 3D track.

Zhang, Shukui↗

Dark Matter and Track Triggering with the CMS Experiment (Final Report)

This report summarizes the progress from DOE Early Career Award entitled “Dark Matter and Track Triggering with the CMS Experiment”. The first major goal of the project was to establish a sensitive Dark Matter (DM) search program for the CMS experiment in Run-2 of the LHC. Two related strategies were developed to accomplish this goal. First, we designed and executed searches for DM produced in association with heavy flavor quark pairs, top/anti-top (ttbar) and bottom/anti-bottom (bbbar). In addition, we sought to maximize the power of Run-2 DM searches through the development of a statistical combination of all major search channels within a consistent theoretical framework. The results of these aspects of the project are detailed in Section 1. The second major goal of the project was to develop a real-time “Level-1” tracking trigger system for the high-luminosity LHC (HL-LHC) CMS upgrade. Charged particle tracking in the first stage of the CMS trigger will be crucial for surviving the high-pileup environment expected at the HL-LHC. The L1 tracking trigger must process all front-end hit data sent from the inner Tracker and will have just 4us to output track primitives to the downstream L1 trigger without data loss. Our development of the CMS tracking trigger is described in Section 2. The team supported by the award is given in Section 3. The project resulted in 10 peer reviewed publications, a graduate student thesis, and contributions to the CMS Phase-2 Tracker Technical Design Report, as is detailed in Section 4. Our development work for the CMS tracking trigger has been integrated in the backend system architecture for the Tracker upgrade.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Measurement of the primary Lund jet plane density in proton-proton collisions at $ \sqrt{\textrm{s}} $ = 13 TeV

A measurement is presented of the primary Lund jet plane (LJP) density in inclusive jet production in proton-proton collisions. The analysis uses 138 fb$^{−1}$ of data collected by the CMS experiment at $ \sqrt{s} $ = 13 TeV. The LJP, a representation of the phase space of emissions inside jets, is constructed using iterative jet declustering. The transverse momentum k$_{T}$ and the splitting angle ∆R of an emission relative to its emitter are measured at each step of the jet declustering process. The average density of emissions as function of ln(k$_{T}$/GeV) and ln(R/∆R) is measured for jets with distance parameters R = 0.4 or 0.8, transverse momentum p$_{T}$> 700 GeV, and rapidity |y| < 1.7. The jet substructure is measured using the charged-particle tracks of the jet. The measured distributions, unfolded to the level of stable charged particles, are compared with theoretical predictions from simulations and with perturbative quantum chromodynamics calculations. Due to the ability of the LJP to factorize physical effects, these measurements can be used to improve different aspects of the physics modeling in event generators.[graphic not available: see fulltext]

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Numerical investigation of inertial sphere impaction on a cylinder in subcritical flow

We investigated the accuracy of a Lagrangian scheme for simulating both the motion and wall impaction of neutrally buoyant inertial spheres (8.85 ≤ St ≤ 18.04) moving in the subcritical regime (Re = 10,972 and 22,366) of flows past a circular cylinder. The accuracy of the eddy-resolving simulations of the flow field was verified based on available data at nearly comparable flow conditions. The accuracy of the eddy-resolving simulations of the flow field was verified based on available data at nearly comparable flow conditions. The accuracy of the particle tracking scheme was verified with respect to corresponding experimental data in a laboratory flume. Particle dispersion en route to the circular cylinder is mostly driven by the flow disturbance from the particle injector tube. The trends in collision rates are largely explained by the rate of particle entertainment into the boundary layer around the cylinder and the inertial condition of the spheres in this proximity region—both of which can be better assessed by post-processing of the numerical outcomes. Here, the primary incentive to undertake this work was the need for validating a particle-based method for quantifying the hydraulic stressors leading to bodily damage of fish passing flows influenced by hydraulic infrastructure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Using CUBIT to Create Unstructured Mesh Models for MCNP Simulations

The Monte Carlo N-Particle (MCNP) transport code version 6 (also known as MCNP6) has the capability for tracking particles on unstructured mesh (UM) geometry models embedded into constructive solid geometry (CSG) cells. This feature has been developed for performing calculations of complex geometry models because creating CSG models is a time-consuming and error-prone process as the complexities of geometries increase. An MCNP UM calculation requires UM geometry input files. The UM capability was originally designed to work with UM models created with the Abaqus/CAE software and ASCII input files that it generates. The Abaqus-formatted input files needed for MCNP UM calculations must have the correct Abaqus syntax and meet the additional MCNP requirements. Several software packages can generate a UM model formatted as an Abaqus input file. Cubit, the Sandia National Laboratory automated mesh generation toolkit, can generate a UM model formatted as an Abaqus input file. However, the Abaqus input files created by Cubit cannot be used for MCNP simulations. A Python script has been developed to convert an Abaqus file created by Cubit to an Abaqus file that the MCNP code can process. This report describes the process of using Cubit to create UM models for MCNP calculations.

97 MATHEMATICS AND COMPUTING↗

Processing MCNP Elemental Edit Outputs

The Monte Carlo N-Particle (MCNP) transport code version 6 (also known as MCNP6) has the capability for tracking particles on unstructured mesh (UM) geometry models embedded into constructive solid geometry (CSG) cells. A UM geometry is a collection of elements representing a solid geometry. The first step of MCNP UM modeling is using other software packages to create a finite element mesh representation of a solid 3D geometry. Computer-aided design (CAD) or computer-aided manufacturing (CAM) software is typically used to create a solid geometry model, which is later imported into mesh generation software to create a UM model. The MCNP UM feature was originally designed for models generated by the Abaqus/CAE software. The MCNP code version 6.0 and later can process UM models formatted as Abaqus input files. MCNP can process a UM model consisting of several different element types including linear tetrahedral or hexahedral elements and calculate quantities of interest such as flux and energy deposition at elements. An MCNP UM simulation provides high-fidelity elemental edit (i.e., tally) outputs, which can be further used in multiphysics calculations. The MCNP UM feature was used for multiphysics simulations where quantities of interest calculated by MCNP are used as inputs for heat transfer calculations in Abaqus. MCNP6.3 can produce two types of elemental edit output (EEOUT) file formats: ASCII and HDF5. An EEOUT file type must be requested on an EMBED card while output type (flux or energy deposition) must be requested on an EMBEE card. We wrote Python3 scripts to extract energy deposition values in an ASCII or HDF5 EEOUT file and compute a heat flux profile for an Abaqus heat transfer calculation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Athena-I CUBIT Journal Files

The Monte Carlo N-Particle (MCNP)1 transport code version 6 (also known as MCNP6) has the capability for tracking particles on unstructured mesh (UM) geometry models embedded into con structive solid geometry (CSG) cells. This feature has been developed for performing calculations of complex geometry models because manually creating CSG models is a time-consuming and error prone process as the complexities of geometries increase. The UM capability was originally designed to work with UM models created with the Abaqus software and ASCII input files that it generates. The MCNP code version 6.0 and later can process UM models formatted as Abaqus input files. Starting with version 6.3, MCNP can also process HDF5 mesh input files. External codes must be used to generate Abaqus input files for MCNP UM calculations. CUBIT, the Sandia National Laboratory automated mesh generation toolkit, can generate a UM model formatted as an Abaqus input file. However, the Abaqus input files exported from CUBIT cannot be used for MCNP simulations because it lacks the proper syntax. A Python script was developed to convert an Abaqus file created by CUBIT to an Abaqus file format that MCNP can process. Creating UM models for complex geometries is not an easy task. The process of creating UM models in CUBIT for MCNP simulations is detailed in. CUBIT provides several user interface options including a graphical user interface (GUI) and a command line interface. A GUI provides an easy way to use CUBIT without learning the CUBIT command syntax. When using CUBIT with either interface option, command lines are written into an ASCII file known as a journal file; this journal file can be edited and archived so that it can be played back in CUBIT to automatically generate a UM model. This report describes the CUBIT journal files of the UM models developed for Athena-I. The Athena platform, an energy-tuning assembly, was developed to spectrally shape the National Ignition Facility (NIF) deuterium-tritium fusion neutron source to a thermonuclear (fusion) plus prompt fission neutron spectrum with capability to act as a short-pulse neutron source. MCNP6 was used for the Athena experiment design analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Generating MCNP Input Files for Unstructured Mesh Geometries

The Los Alamos National Laboratory’s (LANL) Monte Carlo N-Particle (MCNP)1 transport code version 6.3 (also known as MCNP6.3) has the capability for tracking particles on unstructured mesh (UM) geometry models embedded into constructive solid geometry (CSG) cells. This feature has been developed for performing calculations of complex geometry models because manually creating CSG models is time-consuming and error-prone as the complexities of geometries increase. A UM geometry model is a collection of finite elements representing a solid geometry. The first step of the MCNP UM calculation is using other software packages to create a finite element mesh representation of a solid 3D geometry because the MCNP code cannot be used to generate a UM model. Computer-aided design (CAD) software is typically used to create a solid geometry model, which is later imported into mesh generation software to create a UM model. Some mesh generation software packages may also be used to create solid geometries and thus CAD files are not needed. The MCNP UM feature was originally designed for models generated by the Abaqus/CAE software suite. The MCNP code version 6.0 and later can process UM models formatted as Abaqus input files. Starting with a 6.3 version, the MCNP code can process HDF5 mesh input files. We only focus on the UM models formatted as Abaqus input files in this report since currently no external software can be used to generate HDF5 mesh input files for MCNP UM calculations. The MCNP code version 6.3 can be used to convert the Abaqus mesh input files into the HDF5 mesh input files, but this option is typically used by the MCNP code development team to test the HDF5 mesh input file feature. Several software packages (such as Abaqus, Attila4MC, or Cubit) can be used to create the Abaqus input files for MCNP UM calculations. An MCNP UM calculation using an Abaqus model requires two input file types: MCNP and Abaqus input files. The Abaqus input files needed for MCNP UM calcu lations must have the correct Abaqus syntax and meet the additional requirements by the MCNP code. The MCNP code can process only Abaqus input files that make use of part and assembly definitions, where elements in each part must be grouped into one or more element sets (i.e., elset) using *Elset keyword lines with specified naming formats. The MCNP and Abaqus input files required for MCNP UM simulations must be related; pseudo-cells in an MCNP input file must be constructed from mesh model data from an Abaqus input file. For large complex UM models, it is tedious to manually create MCNP UM input files. The um pre op (unstructured mesh pre operations) program with the -m option can be used to create a skeleton MCNP input file from an Abaqus input file [6]. Since the um pre op program was written in Fortran and was not written for optimized performance, this program is a deprecated feature in the MCNP code version 6.3 and may be removed in the next release of the code. To improve calculation flow of multiphysics calculations, a Python3 code called write mcnp um input has been developed to generate an MCNP input file instead of using the um_pre_op -m option. This Python code was initially released to the public in 2020. We have updated this Python code for MCNP6.3 and it was used to generate the MCNP input files used to verify the MCNP6.3 code. The write_mcnp_um_input code is included with the MCNP6.3 code package which will be released to the public through the Radiation Safety Information Computational Center (RSICC) at Oak Ridge National Laboratory. This report is a revision of LA-UR-20-27139 report.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗