Search NASA⌕ Search

SEARCH · Search NASA

Results for “pythia”

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.

53 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Production cross sections of light and charmed mesons in 𝑒 + ⁢𝑒 − annihilation near 10.58 GeV

We report measurements of production cross sections for 𝜌 + , 𝜌 0 , 𝜔, 𝐾* + , 𝐾* 0 , 𝜙, 𝜂, 𝐾$^0_𝑆$, 𝑓 0 ⁡(980), 𝐷 + , 𝐷 0 , 𝐷$^+_𝑠$, 𝐷* + , 𝐷* 0 , and 𝐷*$^+_𝑠$ in 𝑒 + ⁢𝑒 − collisions at a center-of-mass energy near 10.58 GeV. The data were recorded by the Belle experiment, consisting of 571 fb −1 at 10.58 GeV and 74 fb −1 at 10.52 GeV. Production cross sections are extracted as a function of the fractional hadron momentum 𝑥 𝑝 . The measurements are compared to pythia Monte Carlo generator predictions with various fragmentation settings, including those that have increased fragmentation into vector mesons over pseudoscalar mesons. The cross sections measured for light hadrons are consistent with no additional increase of vector over pseudoscalar mesons. The charmed-meson cross sections are compared to earlier measurements—when available—including older Belle results, which they supersede. They are in agreement before application of an improved initial-state radiation correction procedure that causes slight changes in their 𝑥 𝑝 shapes.

fragmentation functions↗
DOE OSTI2025

First measurement of 𝐷 𝑠⁢1 ⁡(1 + )⁢(2536) + and 𝐷$^{*}_{𝑠⁢2}$⁡(2 + )⁢(2573) + production in proton-proton collisions at $\sqrt{s}$ = 13 TeV at the LHC

The production yields of the orbitally excited charm-strange mesons D s⁢1 ⁡(1 + )⁢(2536)+ and 𝐷$^{*}_{𝑠⁢2}$⁡(2 + )⁢(2573) + were measured for the first time in proton-proton (pp) collisions at a center-of-mass energy of $\sqrt{s}$ =13 TeV with the ALICE experiment at the LHC. The D$^{+}_{s⁢1}$ and D$^{*+}_{s⁢2}$ mesons were measured at midrapidity (|𝑦| < 0.5) in minimum-bias and high-multiplicity pp collisions in the transverse-momentum interval 2< 𝑝 T < 24 GeV/𝑐. Their production yields relative to the D$^{+}_{s}$ ground-state yield were found to be compatible between minimum-bias and high-multiplicity collisions, as well as with previous measurements in e ± ⁢p and e + ⁢e − collisions. The measured D$^{+}_{s⁢1}$/D$^{+}_{s}$ and D$^{*+}_{s⁢2}$/D$^{+}_{s}$ yield ratios are described by statistical hadronization models and can be used to tune the parameters governing the production of excited charm-strange hadrons in Monte Carlo generators, such as pythia 8.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Minimizing Selection Bias in Inclusive Jets in Heavy-Ion Collisions with Energy Correlators

The first-ever measurement of energy correlators within inclusive jets produced in heavy-ion collisions, revealed by the CMS Collaboration, shows a clear enhancement at large angles relative to the proton-proton (p-p) baseline. However, interpreting this enhancement is complicated due to selection bias from energy loss, which also distorts the energy correlator heavy-ion to p-p ratio in the hadronization region, hindering our understanding of parton/hadron dynamics in a colored medium. In this Letter, we introduce a new ratio of energy correlator observables that removes the leading effects of selection bias from the two-point energy correlator spectrum (E2C). PYTHIA and HERWIG simulations show that the impact of selection bias in the E2C is reduced by an order of magnitude, while sensitivity to any other medium modifications is retained. This quantity can be obtained directly from the experimental measurements presented by CMS, as illustrated in an accompanying note [C. Andres and J. Holguin, arXiv:2409.07526].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Soft-hard framework with exact four-momentum conservation for small systems

A new framework, called x-scape, for the combined study of both hard and soft transverse momentum sectors in high-energy proton-proton (𝑝−𝑝) and proton-nucleus (𝑝−𝐴) collisions is set up. A dynamical initial state is set up using the 3d-Glauber model with transverse locations of hotspots within each incoming nucleon. A hard scattering that emanates from two colliding hotspots is carried out using the Pythia generator. Initial state radiation from the incoming hard partons is carried out in a new module called I-matter, which includes the longitudinal location of initial splits. The energy-momentum of both the initial hard partons and their associated beam remnants is removed from the hot spots, depleting the energy-momentum available for the formation of the bulk medium. Outgoing showers are simulated using the matter generator, and results are presented for both cases, allowing for and not allowing for energy loss. First comparisons between this hard-soft model and single inclusive hadron and jet data from 𝑝−𝑝 and minimum bias 𝑝−Pb collisions are presented. Single hadron spectra in 𝑝−𝑝 are used to carry out a limited (in number of parameters) Bayesian calibration of the model. Fair comparisons with data are indicative of the utility of this new framework. Theoretical studies of the correlation between jet 𝑝 𝑇 and event activity at mid and forward rapidity are carried out.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Measurement of the inclusive isolated-photon production cross section in pp and Pb–Pb collisions at $\mathbf {\sqrt{\textit{s}_{NN }} = 5.02}$ TeV

The ALICE Collaboration at the CERN LHC has measured the inclusive production cross section of isolated photons at midrapidity as a function of the photon transverse momentum (p$^{γ}_{T}$), in Pb–Pb collisions in different centrality intervals, and in pp collisions, at centre-of-momentum energy per nucleon pair of $\sqrt{s_{NN}}$ = 5.02 TeV. The photon transverse momentum range is between 10–14 and 40–140 GeV/c, depending on the collision system and on the Pb–Pb centrality class. The result extends to lower p$^{γ}_{T}$ than previously published results by the ATLAS and CMS experiments at the same collision energy. The covered pseudorapidity range is |η γ | < 0.67. The isolation selection is based on a charged particle isolation momentum threshold p$^{iso, ch}_{T}$ GeV/c within a cone of radii R = 0.2 and 0.4. The nuclear modification factor is calculated and found to be consistent with unity in all centrality classes, and also consistent with the HG-PYTHIA model, which describes the event selection and geometry biases that affect the centrality determination in peripheral Pb–Pb collisions. The measurement is compared to next-to-leading order perturbative QCD calculations and to the measurements of isolated photons and Z 0 bosons from the CMS experiment, which are all found to be in agreement.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Investigating late-stage particle production in pp collisions with balance functions

Balance functions have been regarded in the past as a method of investigating the late-stage hadronization found in the presence of a strongly-coupled medium. They are also used to constrain mechanisms of particle production in large and small collision systems. Charge balance functions for inclusive and identified particle pairs are reported as a function of charged particle multiplicity in proton–proton collisions simulated with the PYTHIA8 and the EPOS4 models. The charge balance functions of inclusive, pion, kaon, and proton pairs exhibit amplitudes and shapes that depend on particle species and differ significantly in the two models due to the different particle production mechanisms implemented in PYTHIA and EPOS. The shapes and amplitudes also evolve with multiplicity in both models. In addition, the evolution of the longitudinal rms width and that of balance functions integrals with multiplicity (and average transverse momentum) feature significant differences in the two models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

System size and energy dependence of the mean transverse momentum fluctuations at the LHC

Event-by-event fluctuations of the event-wise mean transverse momentum, $\langle$p T $\rangle$, of charged particles produced in proton–proton (pp) collisions at $\sqrt{s}$ = 5.02 TeV, Xe–Xe collisions at $\sqrt{s_{NN}}$ = 5.44 TeV, and Pb–Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV are studied using the ALICE detector based on the integral correlator $\langle$$\langle$Δp T Δp T $\rangle$$\rangle$. The correlator strength is found to decrease monotonically with increasing produced charged-particle multiplicity measured at midrapidity in all three systems. In Xe–Xe and Pb–Pb collisions, the multiplicity dependence of the correlator deviates significantly from a simple power-law scaling as well as from the predictions of the HIJING and AMPT models. The observed deviation from power-law scaling is expected from transverse radial flow in semicentral to central Xe–Xe and Pb–Pb collisions. In pp collisions, the correlation strength is also studied by classifying the events based on the transverse spherocity, S 0 , of the particle production at midrapidity, used as a proxy for the presence of a pronounced back-to-back jet topology. Low-spherocity (jetty) events feature a larger correlation strength than those with high spherocity (isotropic). The strength and multiplicity dependence of jetty and isotropic events are well reproduced by calculations with the PYTHIA 8 and EPOS LHC models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

CROCUS Air Quality Data at Argonne National Laboratory Prairie Site

The AQT (Vaisala AQT530) instrument provides observations on meteorological conditions, including particulate matter (PM2.5, PM10), gas species concentrations (NO, NO2, O3, CO), and environment temperature and moisture. These measurements are critical for understanding air quality. These measurements are useful for understanding changes in aerosol properties, air quality research, and comparing to model experiments especially in urban environments. These measurements are collected at the Argonne Testbed for Multiscale Observational Science (ATMOS), a prairie field site at Argonne National Laboratory in Lemont, Illinois. Data is available in the netCDF data format, we encourage data users review documentation through Project Pythia to understand how to work with netCDF data https://foundations.projectpythia.org/core/data-formats/netcdf-cf.html. Each file contains one day's worth of data (24 hours, starting at 0000 UTC). The data is aggregated into daily frequency to make it easier to process multiple days, and compress the higher-resolution fields. File naming convention includes the project (CROCUS), location (atmos), data level (raw, a1), date (year, month, day), and hour (0000).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

CROCUS Weather Data at Chicago State University Prairie Site

Vaisala WXT sensor is an all-in-one weather instrument that provides 6 of the most important weather parameters: barometric pressure, temperature, relative humidity, rainfall, wind speed and direction. Temperature, pressure, relative humidity, and rainfall are sampled at 1 second frequency, while wind speed/direction is measured at ten per second (10Hz) frequency. These measurements are useful for looking at characterizing local weather, identifying unique weather events, and studying local turbulence, especially given the high temporal resolution of the wind measurements. These measurements are collected at Chicago State University (CSU) at a prarie site on campus in Chicago, Illinois. Data is available in the netCDF data format, we encourage data users review documentation through Project Pythia to understand how to work with netCDF data https://foundations.projectpythia.org/core/data-formats/netcdf-cf.html. The data is aggregated into daily frequency to make it easier to process multiple days, and compress the higher-resolution fields. Each file contains one day's worth of data (24 hours, starting at 0000 UTC). File naming convention includes the project (CROCUS), location (CSU), data level (raw, a1), date (year, month, day), and hour (0000).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

CROCUS Weather Data at Argonne National Laboratory Prairie Site

Vaisala WXT sensor is an all-in-one weather instrument that provides 6 of the most important weather parameters: barometric pressure, temperature, relative humidity, rainfall, wind speed and direction. Temperature, pressure, relative humidity, and rainfall are sampled at 1 second frequency, while wind speed/direction is measured at ten per second (10Hz) frequency. These measurements are useful for looking at characterizing local weather, identifying unique weather events, and studying local turbulence, especially given the high temporal resolution of the wind measurements. These measurements are collected at the Argonne Testbed for Multiscale Observational Science (ATMOS), a prairie field site at Argonne National Laboratory in Lemont, Illinois. Data is available in the netCDF data format, we encourage data users review documentation through Project Pythia to understand how to work with netCDF data https://foundations.projectpythia.org/core/data-formats/netcdf-cf.html. The data is aggregated into daily frequency to make it easier to process multiple days, and compress the higher-resolution fields. Each file contains one day's worth of data (24 hours, starting at 0000 UTC). File naming convention includes the project (CROCUS), location (atmos), data level (raw, a1), date (year, month, day), and hour (0000).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

CROCUS Weather Data at University of Illinois - Chicago Tower

Vaisala WXT sensor is an all-in-one weather instrument that provides 6 of the most important weather parameters: barometric pressure, temperature, relative humidity, rainfall, wind speed and direction. Temperature, pressure, relative humidity, and rainfall are sampled at 1 second frequency, while wind speed/direction is measured at ten per second (10Hz) frequency. These measurements are useful for looking at characterizing local weather, identifying unique weather events, and studying local turbulence, especially given the high temporal resolution of the wind measurements. These measurements are collected at the University of Illinois in Chicago, Illinois, on the meteorological tower near the greenhouse on campus. Data is available in the netCDF data format, we encourage data users review documentation through Project Pythia to understand how to work with netCDF data https://foundations.projectpythia.org/core/data-formats/netcdf-cf.html. File naming convention includes the project (CROCUS), location (UIC), data level (raw, a1), date (year, month, day), and hour (0000).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

CROCUS Air Quality Data at University of Illinois - Chicago Tower

The AQT (Vaisala AQT530) instrument provides observations on meteorological conditions, including particulate matter (PM2.5, PM10), gas species concentrations (NO, NO2, O3, CO), and environment temperature and moisture. These measurements are critical for understanding air quality. These measurements are useful for understanding changes in aerosol properties, air quality research, and comparing to model experiments especially in urban environments. These measurements are collected at the University of Illinois in Chicago, Illinois, on the meteorological tower near the greenhouse on campus. Data is available in the netCDF data format, we encourage data users review documentation through Project Pythia to understand how to work with netCDF data https://foundations.projectpythia.org/core/data-formats/netcdf-cf.html. Each file contains one day's worth of data (24 hours, starting at 0000 UTC). File naming convention includes the project (CROCUS), location (UIC), data level (raw, a1), date (year, month, day), and hour (0000).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

CROCUS Dual Polarization Ceilometer Data at Northeastern Illinois University Rooftop

This dataset is from the Department of Energy Office of Science funded project, CROCUS Urban Integrated Field Laboratory (https://crocus-urban.org/). The dual-polarization ceilometer (Vaisala CL61) is an autonomous lidar system operating at 910 nm wavelength, providing valuable measurements for understanding atmospheric boundary layer evolution, air quality, and cloud-aerosol interactions. The CL61 measures the backscattered signal intensity alternating between parallel- and cross-polarization signals. The unique depolarization measurement capability improves discrimination between different particle types, such as liquid droplets, ice crystals, and aerosols. The depolarization is highly dependent on the scatterer shape and orientation (spherical vs non-spherical particles), with the linear depolarization ratio providing a measure of dominant backscatter signal component from atmospheric particles at various heights, essentially allowing discrimination between liquid and solid particles. With its efficient optical system, CL61’s improved signal-to-noise ratio compared to traditional ceilometers allows studying detailed vertical profiles of aerosols and clouds up to 15 km height.Datasets are stored in a netCDF data format, and we we encourage users to make use of the associated toolkits available from Unidata (https://www.unidata.ucar.edu/software/netcdf/), Project Pythia (https://foundations.projectpythia.org/core/data-formats/netcdf-cf.html), and our “Instrument Cookbooks” (https://crocus-urban.github.io/instrument-cookbooks) for more information on how to process the metadata-rich datasets.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Post-hoc reweighting of hadron production in the Lund string model

We present a method for reweighting flavor selection in the Lund string fragmentation model. This is the process of calculating and applying event weights enabling fast and exact variation of hadronization parameters on pre-generated event samples. The procedure is post hoc, requiring only a small amount of additional information stored per event, and allowing for efficient estimation of hadronization uncertainties without repeated simulation. Weight expressions are derived from the hadronization algorithm itself, and validated against direct simulation for a wide range of observables and parameter shifts. The hadronization algorithm can be viewed as a hierarchical Markov process with stochastic rejections, a structure common to many complex simulations outside of high-energy physics. This perspective makes the method modular, extensible, and potentially transferable to other domains. We demonstrate the approach in Pythia, including both coverage considerations and timing benefits. For the purpose of this paper, our goal is to develop and demonstrate the the formalism, and we therefore exclude several model variations for baryon production (popcorn model, junction production) needed for proton collisions. These will be the topic of a future paper.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Characterizing the hadronization of parton showers using the HOMER method

We update the HOMER method, a technique to solve a restricted version of the inverse problem of hadronization – extracting the Lund string fragmentation function f(z) from data using only observable information. Here, we demonstrate its utility by extracting f(z) from synthetic Pythia simulations using high-level observables constructed on an event-by-event basis, such as multiplicities and shape variables. Four cases of increasing complexity are considered, corresponding to e + e − collisions at a center-of-mass energy of 90 ≥ v producing either a string stretched between a $q$ and $\bar{q}$ containing no gluons; the same string containing one gluon g with fixed kinematics; the same but the gluon has varying kinematics; and the most realistic case, strings with an unrestricted number of gluons that is the end-result of a parton shower. We demonstrate the extraction of f(z) in each case, with the result of only a relatively modest degradation in performance of the HOMER method with the increased complexity of the string system.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Pythia8 Quark and Gluon Jets (float8 e4m3FN)

A float8 (e4m3FN) quantized version of the quark and gluon jet dataset originally published by Komiske, Metodiev, and Thaler (Zenodo record 3164691). Only the 20-file subset without charm and bottom quark jets is included here. All simulation parameters and jet selection criteria are identical to the original: Pythia 8.226, √s = 14 TeV Quarks from WeakBosonAndParton:qg2gmZq, gluons from WeakBosonAndParton:qqbar2gmZg with the Z decaying to neutrinos FastJet 3.3.0, anti-k_t jets with R = 0.4 p_T^jet ∈ [500, 550] GeV, |y^jet| < 1.7 There are 20 files, each in compressed NumPy format (QG_jets_fp8e4m3fn_0.npz through QG_jets_fp8e4m3fn_19.npz). Each file contains two arrays: X: (100000, M, 4) — 50k quark and 50k gluon jets, randomly sorted, padded to max multiplicity M, with particle features (pt, rapidity, azimuthal angle, pdgid) y: (100000,) — jet labels, gluon = 0, quark = 1 Since NumPy has no native fp8 dtype, X is stored as float32, but the values have been quantized through TensorFlow's float8_e4m3fn type and carry only fp8 precision. The quantization procedure is as follows: a global per-channel scale factor is computed from the absolute maximum value across all 20 chunks (with FP8_MAX = 448.0, the maximum representable value of e4m3FN). Each chunk is then scaled into the fp8 dynamic range, round-tripped through tf.experimental.float8_e4m3fn, and scaled back. This global scaling ensures a consistent quantization grid across the full dataset. The y labels are unchanged. Users should be aware that e4m3FN has limited dynamic range and precision. We recommend verifying this format is appropriate for your application; for a less aggressive reduction see the float16 and float32 versions linked below. If you use this dataset, please cite the original Zenodo record and its associated paper: Komiske, Metodiev, Thaler, Energy Flow Networks: Deep Sets for Particle Jets, JHEP 01 (2019) 121, arXiv:1810.05165

DiLullo, Nicholas [Brown University] (ORCID:000000↗
DOE OSTI2026

HDSense: An efficient method for ranking observable sensitivity

Identifying which observables most effectively constrain model parameters can be computationally prohibitive when considering full likelihoods of many correlated observables. This is especially important for, e.g., hadronization models, where high precision is required to interpret the results of collider experiments. We introduce the High-Dimensional Sensitivity (HDSense) score, a computationally efficient metric for ranking observable sets using only one-dimensional histograms. Derived by profiling over unknown correlations in the Fisher information framework, the score balances total information content against redundancy between observables. We apply HDSense to rank a set observables in terms of their constraining power with respect to five parameters of the Lund string model of hadronization implemented in Pythia using simulated leptonic collider events at the $Z$ pole. Validation against machine-learning--based full-likelihood approximations demonstrates that HDSense successfully identifies near-optimal observable subsets. The framework naturally handles data from multiple experiments with different acceptances and incorporates detector effects. While demonstrated on hadronization models, the methodology applies broadly to generic parameter estimation problems where correlations are unknown or difficult to model.

Assi, Benoît [Cincinnati U.] (ORCID:00000003092433↗