SEARCH · Engineering Papers

Results for “Data unfolding methods”

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.

SourceResearch connectionFrom yearThrough year

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

53 records · Page 3

All sources NASA NTRS DOE OSTI

NASA NTRS1985

s-wave threshold in electron attachment - Observations and cross sections in CCl4 and SF6 at ultralow electron energies

The threshold photoionization method was used to study low-energy electron attachment phenomena in and cross sections of CCl4 and SF6 compounds, which have applications in the design of gaseous dielectrics and diffuse discharge opening switches. Measurements were made at electron energies from below threshold to 140 meV at resolutions of 6 and 8 meV. A narrow resolution-limited structure was observed in electron attachment to CCl4 and SF6 at electron energies below 10 meV, which is attributed to the divergence of the attachment cross section in the limit epsilon, l approaches zero. The results are compared with experimental collisional-ionization results, electron-swarm unfolded cross sections, and earlier threshold photoionization data.

Chutjian, A.↗

DOE OSTI2022

Neutron Leakage Spectra Sensitivities for ICSBEP Benchmarks

Neutron leakage spectra have been measured, simulated, and investigated by many groups. These spectra have many uses, including determining shielding requirements and calculating dose for radiation protection purposes, validating nuclear data, and determining material composition. It has even been proposed to use measurements of neutron leakage spectra to determine the soil composition of Mars. As part of the Los Alamos National Laboratory (LANL) Experiments Underpinned by Computational Learning for Improvements in nuclear Data (EUCLID) Laboratory Directed Research Development (LDRD) project, the authors are working on developing methods to calculate sensitivities to neutron leakage spectra. This may allow nuclear data evaluators to better use neutron leakage spectra data to constrain or adjust nuclear data. This may also be particularly useful as many neutron leakage spectra measurements do not require fissile material and can be performed with well characterized neutron sources. Neutron leakage spectra measurements are usually performed by placing a detector at the outside of a nuclear system. This may be a reactor, a neutron generator, or a neutron source. Certain detection systems can use pulse height data to infer neutron energies, other detector systems rely on other ways of determining neutron energy (for example, a Bonner sphere with multiple moderator thicknesses can be used to measure neutron spectrum). These measurements typically rely on some unfolding of the measured results. For pulsed systems like the "Livermore Pulsed Spheres," time-of-flight information can also help to constrain the neutron energy spectra data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Data Augmentation for Neutron Spectrum Unfolding with Neural Networks

Neural networks require a large quantity of training spectra and detector responses in order to learn to solve the inverse problem of neutron spectrum unfolding. In addition, due to the under-determined nature of unfolding, non-physical spectra which would not be encountered in usage should not be included in the training set. While physically realistic training spectra are commonly determined experimentally or generated through Monte Carlo simulation, this can become prohibitively expensive when considering the quantity of spectra needed to effectively train an unfolding network. In this paper, we present three algorithms for the generation of large quantities of realistic and physically motivated neutron energy spectra. Using an IAEA compendium of 251 spectra, we compare the unfolding performance of neural networks trained on spectra from these algorithms, when unfolding real-world spectra, to two baselines. We also investigate general methods for evaluating the performance of and optimizing feature engineering algorithms.

McGreivy, James (ORCID:0000000321723411)↗

DOE OSTI2020

Inclusive search for highly boosted Higgs bosons decaying to bottom quark-antiquark pairs in proton-proton collisions at $\sqrt{s} =$ 13 TeV

A search for standard model Higgs bosons (H) produced with transverse momentum (p$_{T}$) greater than 450 GeV and decaying to bottom quark-antiquark pairs ($ \mathrm{b}\overline{\mathrm{b}} $) is performed using proton-proton collision data collected by the CMS experiment at the LHC at $ \sqrt{s} $ = 13 TeV. The data sample corresponds to an integrated luminosity of 137 fb$^{−1}$. The search is inclusive in the Higgs boson production mode. Highly Lorentz-boosted Higgs bosons decaying to $ \mathrm{b}\overline{\mathrm{b}} $ are reconstructed as single large-radius jets, and are identified using jet substructure and a dedicated b tagging technique based on a deep neural network. The method is validated with Z →$ \mathrm{b}\overline{\mathrm{b}} $ decays. For a Higgs boson mass of 125 GeV, an excess of events above the background assuming no Higgs boson production is observed with a local significance of 2.5 standard deviations (σ), while the expectation is 0.7. The corresponding signal strength and local significance with respect to the standard model expectation are μ$_{H}$ = 3.7 ± 1.2(stat)$ {}_{-0.7}^{+0.8} $(syst)$ {}_{-0.5}^{+0.8} $(theo) and 1.9 σ. Additionally, an unfolded differential cross section as a function of Higgs boson p$_{T}$ for the gluon fusion production mode is presented, assuming the other production modes occur at the expected rates.[graphic not available: see fulltext]

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

Robust Group Subspace Recovery: A New Approach for Multi-Modality Data Fusion

Robust Subspace Recovery (RoSuRe) algorithm was recently introduced as a principled and numerically efficient algorithm that unfolds underlying Unions of Subspaces (UoS) structure, present in the data. The union of Subspaces (UoS) is capable of identifying more complex trends in data sets than simple linear models. In this work, we build on and extend RoSuRe to prospect the structure of different data modalities individually. We propose a novel multi-modal data fusion approach based on group sparsity which we refer to as Robust Group Subspace Recovery (RoGSuRe). Relying on a bi-sparsity pursuit paradigm and non-smooth optimization techniques, the introduced framework learns a new joint representation of the time series from different data modalities, respecting an underlying UoS model. We subsequently integrate the obtained structures to form a unified subspace structure. The proposed approach exploits the structural dependencies between the different modalities data to cluster the associated target objects. The resulting fusion of the unlabeled sensors’ data from experiments on audio and magnetic data has shown that our method is competitive with other state of the art subspace clustering methods. The resulting UoS structure is employed to classify newly observed data points, highlighting the abstraction capacity of the proposed method.

47 OTHER INSTRUMENTATION↗

NASA NTRS2017

Examination of Multiple Lithologies Within the Primitive Ordinary Chondrite NWA 5717

Northwest Africa 5717 is a primitive (subtype 3.05) ungrouped ordinary chondrite which contains two apparently distinct lithologies. In large cut meteorite slabs, the darker of these, lithology A, looks to host the second, much lighter in color, lithology B (upper left, Fig. 1). The nature of the boundary between the two is uncertain, ranging from abrupt to gradational and not always following particle boundaries. The distinction between the lithologies, beyond the obvious color differences, has been supported by a discrepancy in oxygen isotopes and an incongruity in the magnesium contents of chondrule olivine. Here, quantitative textural analysis and mineralogical methods have been used to investigate the two apparent lithologies within NWA 5717. Olivine grains contained in a thin section from NWA 7402, thought to be paired to 5717, were also measured to re-examine the distinct compositional range among the light and dark areas. Procedure: Particles from a high-resolution mosaic image of a roughly 13x15cm slice of NWA 5717 were traced in Adobe Photoshop. Due to the large size of the sample, visually representative regions of each lithology were chosen to be analyzed. The resulting layers of digitized particles were imported into ImageJ, which was used to measure their area, along with the axes, the angle from horizontal, and the centroid coordinates of ellipses fitted to each particle following the approach. Resulting 2D pixel areas were converted to spherical diameters employing the unfolding algorithm, which outputs a 3D particle size distribution based on digitized 2D size frequency data. Spatstat was used to create kernel density plots of the centroid coordinates for each region. X-ray compositional maps, microprobe analyses, and Mossbauer spectroscopy was conducted on a thin section of NWA 7402, tentatively paired to NWA 5717.

Cato, M. J.↗

NASA NTRSDate not supplied

Probing the Sun with Imaging Spectrographs

EUV and X-ray images of the Sun have revolutionized our understanding of our closest star. With them, we can probe the morphology and temperature structure of the solar atmosphere and see how they evolve as a function of space and time. However, image data cannot be used to determine line-of-sight velocities, abundances, or densities. This information is required to calculate the energy budget of eruptive events, provide boundary conditions for global solar models, and explore fundamental processes occurring in the solar atmosphere. For those diagnostics, we require spectroscopy. Because the structures on the Sun are extended sources, most modern-day spectrometers observe the Sun through long narrow slits. Two-dimensional, spectrally pure solar images with velocity, abundance, and density information are built up by stepping the slit over regions of interest. This method implies that two-dimensional information is highly limited by cadence and the temporal evolution and spatial structure of these parameters can never be truly separated. Both spatial and spectral information can be obtained in a single snapshot with slitless spectrometers, which were often used in the 1950-1970s, but were abandoned due to the difficulty of unfolding the overlapping spatial and spectral information. Thanks to advances in computer processing speeds and machine learning algorithms, there have been several techniques developed to complete the spatial/spectral unfolding, unlocking the full capability of slitless spectrometers for solar observations. The goal of this talk is to give an overview of the capability of such instruments, including recent results from a sounding rocket flight, and demonstrate their usefulness in the next decade of solar observatories and beyond.

Amy Winebarger↗

NASA NTRSDate not supplied

Probing the Sun with Imaging Spectrographs

Amy Winebarger↗

NASA NTRSDate not supplied

Probing the Sun with Imaging Spectrographs

EUV and X-ray images of the Sun have revolutionized our understanding of our closest star. With them, we can probe the structure of the solar atmosphere and see how they evolve as a function of space and time. However, image data cannot be used to determine line-of-sight velocities, abundances, or densities. This information is required to calculate the energy budget of eruptive events, provide boundary conditions for global solar models, and explore fundamental processes occurring in the solar atmosphere. For those diagnostics, we require spectroscopy. Because the structures on the Sun are extended sources, most modern-day spectrometers observe the Sun through long narrow slits. Two-dimensional, spectrally pure solar images with velocity, abundance, and density information are built up by stepping the slit over regions of interest. This method implies that two-dimensional information is highly limited by cadence and the temporal evolution and spatial structure of these parameters can never be truly separated. Both spatial and spectral information can be obtained in a single snapshot with slitless spectrometers, which were often used in the 1950-1970s, but were abandoned due to the difficulty of unfolding the overlapping spatial and spectral information. Thanks to advances in computer processing speeds and machine learning algorithms, there have been several techniques developed to complete the spatial/spectral unfolding, unlocking the full capability of slitless spectrometers for solar observations. The goal of this talk is to give an overview of the capability of such instruments and demonstrate their usefulness in the next decade of solar observatories and beyond.

Amy Winebarger↗

DOE OSTI2025

Measurement of photonuclear jet production in ultraperipheral Pb + Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV with the ATLAS detector

In ultrarelativistic heavy ion collisions at the LHC, each nucleus acts a sources of high-energy real photons that can scatter off the opposing nucleus in ultraperipheral photonuclear (𝛾 + 𝐴) collisions. Hard scattering processes initiated by the photons in such collisions provide a novel method for probing nuclear parton distributions in a kinematic region not easily accessible to other measurements. ATLAS has measured production of dijet and multijet final states in ultraperipheral Pb + Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV using a dataset recorded in 2018 with an integrated luminosity of 1.72 nb −1 . Photonuclear final states are selected by requiring a rapidity gap in the photon direction; this selects events where one of the outgoing nuclei remains intact. Jets are reconstructed using the anti-𝑘 t algorithm with radius parameter, 𝑅 = 0.4. Triple-differential cross sections, unfolded for detector response, are measured and presented using two sets of kinematic variables. The first set consists of the total transverse momentum (𝐻 T ), rapidity, and mass of the jet system. The second set uses 𝐻 T and particle-level nuclear and photon parton momentum fractions, 𝑥 A and 𝑧 𝛾 , respectively. The results are compared with leading-order perturbative QCD calculations of photonuclear jet production cross sections, where all leading order predictions using existing fits fall below the data in the shadowing region. More detailed theoretical comparisons will allow these results to strongly constrain nuclear parton distributions, and these data provide results from the LHC directly comparable to early physics results at the planned Electron-Ion Collider.

Parton distribution functions↗

DOE OSTI2023

Development of a compact fast-neutron spectrometer for nuclear emergency response applications

We have developed a Compact Fast Neutron Spectrometer (CFNS) for passive assay of special nuclear material (SNM) through the observation of fast neutrons. The CFNS consists of eight organic glass scintillators (OGS) coupled to silicon photomultipliers and a waveform digitizer, which are integrated within a human-portable box. The CFNS determines the neutron energy profile by spectrum unfolding using the Maximum-Likelihood Expectation Maximization method. The detector acquisition system was optimized to have a dynamic range of up to 10 MeV neutron energy. Bulk special nuclear material (SNM) measurements from the National Criticality Experiments Research Center were analyzed for SNM validation/examination. Additionally, the results show that the CFNS can be used to distinguish between fission and (α, n) neutron emitters, regardless of intervening material type (Cu and polyethylene) and thickness, by taking the ratio of neutron counts at different regions in the unfolded energy spectrum. Additionally, by fitting an exponential curve to the unfolded energy spectrum of PuO 2 and Pu neutron emitters, the CFNS showed the ability of distinguishing between pure Pu oxide, pure Pu metal and mixed oxide-metal configurations.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2024

Adjusting the Energy Profile for CH–O Interactions Leads to Improved Stability of RNA Stem-Loop Structures in MD Simulations

The role of ribonucleic acid (RNA) in biology continues to grow, but insight into important aspects of RNA behavior is lacking, such as dynamic structural ensembles in different environments, how flexibility is coupled to function, and how function might be modulated by small molecule binding. In the case of proteins, much progress in these areas has been made by complementing experiments with atomistic simulations, but RNA simulation methods and force fields are less mature. It remains challenging to generate stable RNA simulations, even for small systems where well-defined, thermostable structures have been established by experiments. Further many different aspects of RNA energetics have been adjusted in force fields, seeking improvements that are transferable across a variety of RNA structural motifs. In this work, the role of weak CH···O interactions is explored, which are ubiquitous in RNA structure but have received less attention in RNA force field development. By comparing data extracted from high-resolution RNA crystal structures to energy profiles from quantum mechanics and force field calculations, it is shown that CH···O interactions are overly repulsive in the widely used Amber RNA force fields. A simple, targeted adjustment of CH···O repulsion that leaves the remainder of the force field unchanged was developed. Then, the standard and modified force fields were tested using molecular dynamics (MD) simulations with explicit water and salt, amassing over 300 μs of data for multiple RNA systems containing important features such as the presence of loops, base stacking interactions as well as canonical and noncanonical base pairing. In this work and others, standard force fields lead to reproducible unfolding of the NMR-based structures. Including a targeted CH···O adjustment in an otherwise identical protocol dramatically improves the outcome, leading to stable simulations for all RNA systems tested.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Differential transcriptomic alterations in nasal versus lung tissue of acrolein-exposed rats

Introduction: Acrolein is a significant component of anthropogenic and wildfire emissions, as well as cigarette smoke. Although acrolein primarily deposits in the upper respiratory tract upon inhalation, patterns of site-specific injury in nasal versus pulmonary tissues are not well characterized. This assessment is critical in the design of in vitro and in vivo studies performed for assessing health risk of irritant air pollutants. Methods: In this study, male and female Wistar-Kyoto rats were exposed nose-only to air or acrolein. Rats in the acrolein exposure group were exposed to incremental concentrations of acrolein (0, 0.1, 0.316, 1 ppm) for the first 30 min, followed by a 3.5 h exposure at 3.16 ppm. In the first cohort of male and female rats, nasal and bronchoalveolar lavage fluids were analyzed for markers of inflammation, and in a second cohort of males, nasal airway and left lung tissues were used for mRNA sequencing. Results: Protein leakage in nasal airways of acrolein-exposed rats was similar in both sexes; however, inflammatory cells and cytokine increases were more pronounced in males when compared to females. No consistent changes were noted in bronchoalveolar lavage fluid of males or females except for increases in total cells and IL-6. Acrolein-exposed male rats had 452 differentially expressed genes (DEGs) in nasal tissue versus only 95 in the lung. Pathway analysis of DEGs in the nose indicated acute phase response signaling, Nrf2-mediated oxidative stress, unfolded protein response, and other inflammatory pathways, whereas in the lung, xenobiotic metabolism pathways were changed. Genes associated with glucocorticoid and GPCR signaling were also changed in the nose but not in the lung. Discussion: These data provide insights into inhaled acrolein-mediated sex-specific injury/inflammation in the nasal and pulmonary airways. The transcriptional response in the nose reflects acrolein-induced acute oxidative and cytokine signaling changes, which might have implications for upper airway inflammatory disease susceptibility.

Alewel, Devin I.↗

DOE OSTI2025

Tools for unbinned unfolding

Machine learning has enabled differential cross section measurements that are not discretized. Going beyond the traditional histogram-based paradigm, these unbinned unfolding methods are rapidly being integrated into experimental workflows. Here, in order to enable widespread adaptation and standardization, we develop methods, benchmarks, and software for unbinned unfolding. For methodology, we demonstrate the utility of boosted decision trees for unfolding with a relatively small number of high-level features. This complements state-of-the-art deep learning models capable of unfolding the full phase space. To benchmark unbinned unfolding methods, we develop an extension of existing dataset to include acceptance effects, a necessary challenge for real measurements. Additionally, we directly compare binned and unbinned methods using discretized inputs for the latter in order to control for the binning itself. Lastly, we have assembled two software packages for the OmniFold unbinned unfolding method that should serve as the starting point for any future analyses using this technique. One package is based on the widely-used RooUnfold framework and the other is a standalone package available through the Python Package Index (PyPI).

47 OTHER INSTRUMENTATION↗

NASA NTRSDate not supplied

Functional Task Tests in Partial Gravity During Parabolic Flight

BACKGROUND Understanding how critical mission tasks are performed in partial gravity such as on the moon or Mars is necessary to define effective and comprehensive countermeasure strategies for preserving crew performance during exploration missions. We studied the performance of tasks such as standing, walking, and jumping during the partial gravity phases of parabolic flight. We hypothesized that the acute effects of partial gravity on vestibular, proprioceptive, and sensorimotor functions would negatively impact performance. METHODS Twelve subjects were tested over three flights of 30 parabolas each, including 10 parabolas at 0.25g, 10 parabolas at 0.5g, and 10 parabolas at 0.75g. Subjects also performed tests in 1g between parabolas. During the sit-to-stand with obstacle walk task, subjects rose from a seated position and walked as quickly as possible straight ahead towards a cone (4 m distance), walked around the cone making a 180° left turn, returned, and sat in the chair. On the way to and from the cone, subjects stepped over a 30 cm high obstacle. For the recovery from fall task, subjects lay prone for the pull-up phase and initial 10 sec of the parabola. Then they were asked to rise as quickly as possible and maintain a quiet stance for 10 sec. The tandem rail balance task involved standing with both feet on a 4.5 cm wide rail. The time ended when subjects either stepped off the rail or grabbed on to support straps. The jump down task started with the subjects standing on a 30 cm high platform, then the subjects were instructed to step off the platform, land with both feet simultaneously, and settle in a quiet stance. The cone of stability task involved subjects leaning about the ankles as far as they could in the anterior, posterior, and lateral directions without unfolding their arms or taking a step. The center of pressure distance between endpoints was calculated. Data were collected using inertial measurement units (Opal V2, APDM, Portland, OR) worn on the head and trunk, heart rate monitors (RS800CX, Polar, Kempele, Finland), and a force plate (Bertec, Columbus, OH). RESULTS Gravity level had a significant effect on performance, with the greatest changes from 1g tending to be at the 0.25g level (Table 1). Lower gravity levels were associated with increased times to complete the sit-to-stand obstacle walk task and the recovery from fall task, decreased change in heart rate during the recovery from fall task, decreased rail balance times, and increased cone of stability distance in the anterior-posterior direction. Table 1. Functional task performance at different gravity levels during parabolic flight. Measure0.25g0.5g0.75g1gp-valueSit-to-stand with obstacle walk time (sec)9.8 ±1.4*7.2 ±0.76.9 ±0.8*7.4 ±0.9<0.001Recovery from fall time to settle (sec)5.3 ± 0.9*4.6 ± 0.84.2 ± 0.64.3 ± 0.70.002Recovery from fall change in heart rate (bpm)6.0 ± 7.3*11.2 ± 6.5*14.9 ± 4.915.9 ± 6.5<0.001Rail balance with eyes open time (sec)2.7 ±0.8*4.8 ±2.16.6 ±4.78.4 ±6.70.031Rail balance with eyes closed time (sec)1.6 ±0.4*2.1 ±0.42.4 ±0.62.6 ±0.8<0.001Jump down time to settle (sec)2.1 ± 0.41.9 ± 0.42.0 ± 0.31.8 ± 0.30.328Cone of stability –anterior-posterior (cm)20.8 ±2.7*18.8 ±2.218.6 ±1.717.9 ±2.00.039Cone of stability –lateral (cm)26.8 ±6.624.4 ±2.123.0 ±2.123.8 ±2.30.215p-value: one-way repeated measures analysis of variance; *Significant pairwise difference from 1g (p<0.05). DISCUSSION These data suggest that there is a dose-response relationship between gravity level and functional task performance. The largest changes in performance were expected at the lowest gravity level (0.25g) because subjects would no longer be able to use the gravitational reference for the perception of upright. Understanding the extent of performance deficits informs the risks and design of countermeasures for exploration spaceflight missions. ACKNOWLEDGEMENT This work is supported by NASA’s Human Research Program Human Health Countermeasures Element.

T R Macaulay↗

DOE OSTI2022

Photosynthesis phenology, as defined by solar-induced chlorophyll fluorescence, is overestimated by vegetation indices in the extratropical Northern Hemisphere

Vegetation phenology is highly sensitive to climate change, although the data and methods used to estimate key phenological states can influence this sensitivity. Because of its direct relation to leaf photosynthetic carbon uptake, remotely sensed solar-induced chlorophyll fluorescence (SIF) can provide new insight assessing changes in vegetation phenology. In this work, we investigated the potential of using a SIF time series product named contiguous SIF (CSIF) to estimate spring, summer, and autumn phenology in the extratropical Northern Hemisphere (>30°N) and compared the results with those based on Moderate Resolution Imaging Spectroradiometer (MODIS) Normalized Difference Vegetation Index (NDVI) for the period 2001–2018. Overall, we found similar spatial patterns in phenological states. However, specific dates of key phenological events differed when using CSIF vs. MODIS NDVI data. NDVI data indicated that the growing season started earlier (by 10.1 days on average) and ended later (11.5 days on average) relative to CSIF data. This implies that actual periods for photosynthetic activity are shorter (by 21.6 days on average) than those estimated from vegetation indices more directly related to changes in canopy structure. These large differences between results from NDVI and that from CSIF suggest that vegetation indices such as NDVI seem to overestimate the period for active photosynthesis over the extratropical Northern Hemisphere. Furthermore, while phenology of the early growing season is dominated by temperature for both NDVI and CSIF data, phenology of the late growing season is mainly controlled by temperature for NDVI but by precipitation for CSIF. Our findings were further confirmed by other SIF (GOME-2 SIF) and vegetation index (MODIS EVI) datasets. Phenology modes in Earth system modelling are often parameterized using leaf unfolding and senescence from either station or satellite observations. Our results imply that canopy structure-based parameterization schemes may have overestimated photosynthesis active period, and thus productivity responses. We conclude that SIF data offers a novel and unique approach for assessing phenological change - one that is more directly tied to the carbon cycle and how it is being influenced by climate change.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Neural posterior unfolding

Differential cross section measurements are the currency of scientific exchange in particle and nuclear physics. A key challenge for these analyses is the correction for detector distortions, known as deconvolution or unfolding. Binned unfolding of cross section measurements traditionally rely on the regularized inversion of the response matrix that represents the detector response, mapping pre-detector (`particle level') observables to post-detector (`detector level') observables. In this paper we introduce Neural Posterior Unfolding, a modern, Bayesian approach that leverages normalizing flows for unfolding. By using normalizing flows for neural posterior estimation, NPU offers several key advantages including implicit regularization through the neural network architecture, fast amortized inference that eliminates the need for repeated retraining, and direct access to the full uncertainty in the unfolded result. In addition to introducing NPU, we implement a classical Bayesian unfolding method called Fully Bayesian Unfolding (FBU) in modern Python so it can also be studied. These tools are validated on simple Gaussian examples and then tested on simulated jet substructure examples from the Large Hadron Collider (LHC). We find that the Bayesian methods are effective and worth additional development to be analysis ready for cross section measurements at the LHC and beyond.

Analysis and statistical methods↗