Search NASA⌕ Search

SEARCH · Search NASA

Results for “Methods:statistical”

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 703 records · Page 39

Model independent approach for calculating galaxy rotation curves for low S/N MaNGA galaxies

Internal kinematics of galaxies, traced through the stellar rotation curve or two dimensional velocity map, carry important information on galactic structure and dark matter. With upcoming surveys, the velocity map may play a key role in the development of kinematic lensing as an astrophysical probe. Here, we improve techniques for extracting velocity information from integral field spectroscopy at low signal-to-noise (S/N), without a template, and demonstrate substantial advantages over the standard Penalized PiXel-Fitting method (pPXF) approach. Robust rotation curves can be derived down to S/N ≈ 2 using our method.

79 ASTRONOMY AND ASTROPHYSICS↗

PyCBC Inference: A Python-based Parameter Estimation Toolkit for Compact Binary Coalescence Signals

We introduce new modules in the open-source PyCBC gravitational-wave astronomy toolkit that implement Bayesian inference for compact-object binary mergers. We review the Bayesian inference methods implemented and describe the structure of the modules. Here, we demonstrate that the PyCBC Inference modules produce unbiased estimates of the parameters of a simulated population of binary black hole mergers. We show that the parameters' posterior distributions obtained using our new code agree well with the published estimates for binary black holes in the first Advanced LIGO–Virgo observing run.

79 ASTRONOMY AND ASTROPHYSICS↗

A new algorithm for identifying the flavour of B 0 s mesons at LHCb

A new algorithm for the determination of the initial flavour of B 0 s mesons is presented. The algorithm is based on two neural networks and exploits the b hadron production mechanism at a hadron collider. The first network is trained to select charged kaons produced in association with the B 0 s meson. The second network combines the kaon charges to assign the B 0 s flavour and estimates the probability of a wrong assignment. The algorithm is calibrated using data corresponding to an integrated luminosity of 3 fb –1 collected by the LHCb experiment in proton-proton collisions at 7 and 8 TeV centre-of-mass energies. The calibration is performed in two ways: by resolving the B 0 s – B ¯ 0 s flavour oscillations in B 0 s → D – s π + decays, and by analysing flavour-specific B * s 2 (5840) 0 → B + K – decays. The tagging power measured in B 0 s → D – s π + decays is found to be (1.80 ± 0.19 (stat) ± 0.18 (syst))%, which is an improvement of about 50% compared to a similar algorithm previously used in the LHCb experiment.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

The use of minimal spanning trees in particle physics

Minimal spanning trees (MSTs) have been used in cosmology and astronomy to distinguish distributions of points in a multi-dimensional space. They are essentially unknown in particle physics, however. We briefly define MSTs and illustrate their properties through a series of examples. We show how they might be applied to study a typical event sample from a collider experiment and conclude that MSTs may prove useful in distinguishing different classes of events.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Deuterium charge exchange recombination spectroscopy from the top of the pedestal to the scrape off layer in H-mode plasmas

Recent completion of the thirty two channel main-ion (deuterium) charge exchange recombination spectroscopy (CER) diagnostic on DIII-D enables detailed comparisons between impurity and main-ion temperature, density, and toroidal rotation. Sixteen sightlines cover the core of the plasma and another sixteen are densely packed towards the edge, providing high resolution measurements of the pedestal and steep gradient edge region of H-mode plasmas. The complexities of the D α spectrum require fitting with a comprehensive model, as well as using iterative collisional radiative modeling to determine the underlying thermal deuterium ion properties. Large differences in the structure and magnitude of impurity (C 6+ ) and main-ion (D + ) toroidal rotation profiles are seen in the H-mode pedestal. Additionally the D + temperature can be half the value of the C 6+ temperature at the separatrix and shows more of a pedestal structure. Typically only the impurity properties are measured and the main-ion properties are either assumed to be the same, or inferred using neoclassical models, which require validation in the steep gradient region. Furthermore, these measured differences have implications for transport model validation, intrinsic rotation studies, pedestal stability, and the boundary conditions for scrape off layer and plasma material interactions studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Deep neural networks for energy and position reconstruction in EXO-200

Here, we apply deep neural networks (DNN) to data from the EXO-200 experiment. In the studied cases, the DNN is able to reconstruct the relevant parameters—total energy and position—directly from raw digitized waveforms, with minimal exceptions. For the first time, the developed algorithms are evaluated on real detector calibration data. The accuracy of reconstruction either reaches or exceeds what was achieved by the conventional approaches developed by EXO-200 over the course of the experiment. Most existing DNN approaches to event reconstruction and classification in particle physics are trained on Monte Carlo simulated events. Such algorithms are inherently limited by the accuracy of the simulation. We describe a unique approach that, in an experiment such as EXO-200, allows to successfully perform certain reconstruction and analysis tasks by training the network on waveforms from experimental data, either reducing or eliminating the reliance on the Monte Carlo.

97 MATHEMATICS AND COMPUTING↗

Initial results on energy resolution of the NEXT-White detector

One of the major goals of the NEXT-White (NEW) detector is to demonstrate the energy resolution that an electroluminescent high pressure xenon TPC can achieve for high energy tracks. For this purpose, energy calibrations with 137Cs and 232Th sources have been carried out as a part of the long run taken with the detector during most of 2017. This paper describes the initial results obtained with those calibrations, showing excellent linearity and an energy resolution that extrapolates to approximately 1% FWHM at Q$_{\beta\beta}$.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Characterization of the ePix100a and the FastCCD Semiconductor Detectors for the European XFEL

Here, the European X-ray Free Electron Laser (EuXFEL) is a research facility providing spatially coherent X-ray flashes in the energy range from 0 . 25 keV to 25 keV of unprecedented brilliance and with unique time structure: X-ray pulses with a 4 . 5 MHz repetition rate arranged in trains with 2700 pulses every 100 ms. The facility operates three photon beamlines called SASE 1, SASE 2 and SASE 3. Each of the beamlines is hosting two scientific experiments. The SASE 1 beamline started its user operation in September 2017, followed by successful first lasing at the SASE 2 beamline in May 2018. Early user experiments are planned to start in 2019 at this beamline, while early user experiments for the SASE 3 beamline are scheduled for the end of 2018. The quality of the experimental data will gain substantial benefits from an accurate characterization and calibration of the X-ray detectors. Supplementing high repetition rate detectors at MHz speeds, slower detectors such as the ePix100a and the FastCCD will be operated at the train repetition rate of 10 Hz. These 2D silicon pixelized detectors use fast parallel column-wise readout implemented as a CCD or as a hybrid pixel detector. In the following, characterization and analysis approaches for the FastCCD and the ePix100a detectors are discussed and the performance of the detectors is evaluated using appropriate state-of-the-art analysis techniques.

47 OTHER INSTRUMENTATION↗

Parametrizing the detector response with neural networks

In high energy physics, characterizing the response of a detector to radiation is one of the most important and fundamental experimental tasks. In many cases, this task is accomplished by parameterizing summary statistics of the full detector response probability density. The parameterized detector response can then be used for calibration as well as for directly improving physics analysis sensitivity. This paper reports on how to parameterize summary statistics of the detector response using neural networks. Specifically, neural networks are powerful tools for incorporating multidimensional data and the loss function used during training determines which summary statistic is learned. One common summary statistic that has not been combined with deep learning (as far as the authors are aware) is the mode. A neural network-based approach to mode learning is proposed and empirically demonstrated in the context of high energy jet calibrations. Altogether, the neural network-based toolkit for detector response parameterization can enhance the utility of data collected at high energy physics experiments and beyond.

47 OTHER INSTRUMENTATION↗

Automated detector simulation and reconstruction parametrization using machine learning

Rapidly applying the effects of detector response to physics objects (e.g. electrons, muons, showers of particles) is essential in high energy physics. Presently available tools for the transformation from truth-level physics objects to reconstructed detector-level physics objects involve manually defining resolution functions. These resolution functions are typically derived in bins of variables that are correlated with the resolution (e.g. pseudorapidity and transverse momentum). This process is time consuming, requires manual updates when detector conditions change, and can miss important correlations. Machine learning offers a way to automate the process of building these truth-to-reconstructed object transformations and can capture complex correlation for any given set of input variables. Such machine learning algorithms, with sufficient optimization, could have a wide range of applications: improving phenomenological studies by using a better detector representation, allowing for more efficient production of Geant4 simulation by only simulating events within an interesting part of phase space, and studies on future experimental sensitivity to new physics.

47 OTHER INSTRUMENTATION↗

Fast inference of Boosted Decision Trees in FPGAs for particle physics

We describe the implementation of Boosted Decision Trees in the hls4ml library, which allows the translation of a trained model into FPGA firmware through an automated conversion process. Thanks to its fully on-chip implementation, hls4ml performs inference of Boosted Decision Tree models with extremely low latency. With a typical latency less than 100 ns, this solution is suitable for FPGA-based real-time processing, such as in the Level-1 Trigger system of a collider experiment. These developments open up prospects for physicists to deploy BDTs in FPGAs for identifying the origin of jets, better reconstructing the energies of muons, and enabling better selection of rare signal processes.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Geant4Reweight: a framework for evaluating and propagating hadronic interaction uncertainties in Geant4

Geant4Reweight is an open-source C++ framework that allows users to weight tracks produced by the Geant4 particle transport Monte Carlo simulation according to hadron interaction cross section variations and estimate uncertainties in Geant4 interaction models by comparing the simulation’s hadron interaction cross section predictions to data. The ability to weight hadron transport as simulated by Geant4 is crucial to the propagation of systematic uncertainties related to secondary hadronic interactions in current and upcoming neutrino oscillation experiments, including MicroBooNE, NOvA, and DUNE, aswell as hadron test beam experiments such as ProtoDUNE. Here, we provide motivation for weighting hadron tracks in Geant4 in the context of systematic uncertainty propagation, a description of Geant4’s transport simulation technique, and a description of our weighting technique and fitting framework in the momentum range 0–10 GeV/c, which is typical for the hadrons produced by neutrino interactions in these experiments.

97 MATHEMATICS AND COMPUTING↗

Optical calibration of the SNO+ detector in the water phase with deployed sources

SNO+ is a large-scale liquid scintillator experiment with the primary goal of searching for neutrinoless double beta decay, and is located approximately 2 km underground in SNOLAB, Sudbury, Canada. The detector acquired data for two years as a pure water Cherenkov detector, starting in May 2017. During this period, the optical properties of the detector were measured in situ using a deployed light diffusing sphere, with the goal of improving the detector model and the energy response systematic uncertainties. The measured parameters included the water attenuation coefficients, effective attenuation coefficients for the acrylic vessel, and the angular response of the photomultiplier tubes and their surrounding light concentrators, all across different wavelengths. Finally, the calibrated detector model was validated using a deployed tagged gamma source, which showed a 0.6% variation in energy scale across the primary target volume.

47 OTHER INSTRUMENTATION↗

Single electron in a storage ring: a probe into the fundamental properties of synchrotron radiation and a powerful diagnostic tool

An experimental study of undulator radiation generated by a single electron was carried out at the Integrable Optics Test Accelerator (IOTA) storage ring at Fermilab. Photons were detected with a single-photon avalanche diode (SPAD) at an average rate of 20 kHz, corresponding to one detection every few hundred revolutions. Events were continuously recorded by a picosecond event timer for as long as 5 minutes at a time. The collected data were used to directly test the predictions of classical Poisson photostatistics. The study was motivated in part by a previous observation of sub-Poissonian statistics in a similar experiment with synchrotron radiation. In addition, we show that the time series of recorded events can be used to study the synchrotron motion of a single electron, yielding dynamical parameters such as the phase jitter of the radiofrequency cavity and the dependence of the synchrotron period on amplitude. Research plans for further experiments are also discussed.

47 OTHER INSTRUMENTATION↗

Muon track reconstruction in a segmented bolometric array using multi-objective optimization

Recent advances in segmented solid-state detector arrays for rare-event searches have allowed the technology to approach the ton-scale in detector mass and the scale of meters in size. Often focused around searches for neutrinoless double-beta decay or direct dark matter detection, such experiments also have the capability to search for exotic particles that leave track-like signatures across their volume. However, the segmented nature of such detector arrays often sets the spatial resolution and makes the problem of reconstructing track-like paths non-trivial. Here, in this paper, we present an algorithm that improves reconstruction of track-like events in segmented detectors using multi-objective optimization — a computational technique that optimizes more than one cost function at a time without specifying a quantitative weighting between them. Such a technique allows the reconstruction of tracks through a detector and the determination of path-lengths through individual elements. When combined with the reconstructed energy depositions in each element this allows for a calculation of the stopping power of track-like particles and opens the door to searches for particles with abnormal stopping power like monopoles or lightly-ionizing particles (LIPs). Results are presented which evaluate the precision of the reconstruction tools as they currently stand against Monte Carlo generated data. The algorithm is presented in the context of the CUORE experiment, but has applications to other segmented calorimeter detectors.

47 OTHER INSTRUMENTATION↗

FlameNEST: explicit profile likelihoods with the Noble Element Simulation Technique

We present FlameNEST, a framework providing explicit likelihood evaluations in noble element particle detectors using data-driven models from the Noble Element Simulation Technique. FlameNEST provides a way to perform statistical analyses on real data with no dependence on large, computationally expensive Monte Carlo simulations by evaluating the likelihood on an event-by-event basis using analytic probability elements convolved together in a single TensorFlow multiplication. Furthermore, this robust framework creates opportunities for simple inter-collaboration analyses which will be fundamental for the future of experimental dark matter physics.

47 OTHER INSTRUMENTATION↗

Calomplification — the power of generative calorimeter models

Motivated by the high computational costs of classical simulations, machine-learned generative models can be extremely useful in particle physics and elsewhere. They become especially attractive when surrogate models can efficiently learn the underlying distribution, such that a generated sample outperforms a training sample of limited size. This kind of GANplification has been observed for simple Gaussian models. We show the same effect for a physics simulation, specifically photon showers in an electromagnetic calorimeter.

47 OTHER INSTRUMENTATION↗

Characterisation of ion temperature and toroidal rotation on the ST40 tokamak

ST40 is a high field low-aspect ratio spherical tokamak built and operated by Tokamak Energy Ltd. Recent plasma operations were aimed at exploring operational scenarios to maximise the central plasma temperature and have culminated in the achievement of thermal ion temperatures of over 9 keV. Here this paper presents ion temperature and toroidal rotation measurements performed on ST40 during the 2021–22 campaign for a range of different scenarios. Several independent diagnostic systems are used, analysing their correlation and interpreting their differences using new diagnostic forward models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗