Monte Carlo Event Generator Research Articles

Precise measurements of W- and Z-boson transverse momentum spectra with the ATLAS detector using pp collisions at s=5.02\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\sqrt{s} = 5.02$$\\end{document} TeV and 13 TeV

This paper describes measurements of the transverse momentum spectra of W and Z bosons produced in proton–proton collisions at centre-of-mass energies of s=5.02\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\sqrt{s}=5.02$$\\end{document} TeV and s=13\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\sqrt{s}=13$$\\end{document} TeV with the ATLAS experiment at the Large Hadron Collider. Measurements are performed in the electron and muon channels, W→ℓν\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$W\\rightarrow \\ell \ u $$\\end{document} and Z→ℓℓ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$Z\\rightarrow \\ell \\ell $$\\end{document} (ℓ=e\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\ell =e$$\\end{document} or μ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\mu $$\\end{document}), and for W events further separated by charge. The data were collected in 2017 and 2018, in dedicated runs with reduced instantaneous luminosity, and correspond to 255 and 338 pb-1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\hbox {pb}^{-1}$$\\end{document} at s=5.02\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\sqrt{s}=5.02$$\\end{document} TeV and 13 TeV, respectively. These conditions optimise the reconstruction of the W-boson transverse momentum. The distributions observed in the electron and muon channels are unfolded, combined, and compared to QCD calculations based on parton shower Monte Carlo event generators and analytical resummation. The description of the transverse momentum distributions by Monte Carlo event generators is imperfect and shows significant differences largely common to W-\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$W^-$$\\end{document}, W+\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$W^+$$\\end{document} and Z production. The agreement is better at s=5.02\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\sqrt{s}=5.02$$\\end{document} TeV, especially for predictions that were tuned to Z production data at s=7\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\sqrt{s}=7$$\\end{document} TeV. Higher-order, resummed predictions based on DYTurbo generally match the data best across the spectra. Distribution ratios are also presented and test the understanding of differences between the production processes.

Quark spin effects in e+e− annihilation: A Monte Carlo event generator study

Quark spin effects in e+e− annihilation to pseudoscalar and vector mesons are implemented for the first time in the Monte Carlo event generator. The spin-dependent fragmentation of the string stretched between the produced quark-antiquark pair with correlated spin states is described by the string+P30 model implemented in the string fragmentation routine of by using the package. The simulated events are used to study the model predictions for the Collins asymmetries of mesons produced back-to-back in the e+e− center of mass system by using both the thrust axis method and the hadronic plane method. The obtained asymmetries are compared to the available data from the BELLE and experiments and the underlying Collins analyzing power from the string+P30 model is compared with phenomenological extractions. Published by the American Physical Society 2024

Cross Section Evaluation for Exclusive Channels of K+Λ and K+Σ0 Electroproduction off Protons using CLAS Detector Data

In this work, a method for evaluating the cross sections of electroproduction of K+Λ and K+Σ 0 off protons in the region of invariant masses of final hadrons MK + MY < W < 2.65 GeV (MK and MY being the masses of the kaon and hyperon, respectively) and squares of four-momentum transfers of virtual photons, i.e., photon virtualities 0 < Q2 < 5 GeV2 , is developed based on experimental data of these exclusive channels’ cross sections measured by the CLAS detector in Hall B at Jefferson Lab. A set of algorithms has been implemented to evaluate the differential cross sections of these channels, along with their statistical and systematic uncertainties. A program was developed for the evaluation of differential cross sections and structure functions using C++ and Python libraries. An interactive website was created for working with the program, enabling the analysis of one-dimensional and two-dimensional dependences of structure functions and differential cross sections. The evaluation of differential cross sections for the K+Λ and K+Σ 0 electroproduction channels is necessary for extracting the structure function σLT ′ from the data on the polarization asymmetry of electroproduction reactions of these final states with longitudinally polarized electrons. The obtained results are also important for the development of realistic Monte Carlo event generators in planning future experiments and for evaluating the efficiency of detecting final particles when extracting reaction cross sections from experimental data

Herwig 7.3 release note

A new release of the Monte Carlo event generator Herwig (version 7.3) has been launched. This iteration encompasses several enhancements over its predecessor, version 7.2. Noteworthy upgrades include: the implementation of a process-independent electroweak angular-ordered parton shower integrated with QCD and QED radiation; a new recoil scheme for initial-state radiation improving the behaviour of the angular-ordered parton shower; the incorporation of the heavy quark effective theory to refine the hadronization and decay of excited heavy mesons and heavy baryons; a dynamic strategy to regulate the kinematic threshold of cluster splittings within the cluster hadronization model; several improvements to the structure of the cluster hadronization model allowing for refined models; the possibility to extract event-by-event hadronization corrections in a well-defined way; the possibility of using the string model, with a dedicated tune. Additionally, a new tuning of the parton shower and hadronization parameters has been executed. This article discusses the novel features introduced in version 7.3.0.

Observation of Enhanced Long-Range Elliptic Anisotropies Inside High-Multiplicity Jets in pp Collisions at sqrt[s]=13 TeV.

A search for collective effects inside jets produced in proton-proton collisions is performed via correlation measurements of charged particles using the CMS detector at the CERN LHC. The analysis uses data collected at a center-of-mass energy of sqrt[s]=13 TeV, corresponding to an integrated luminosity of 138 fb^{-1}. Jets are reconstructed with the anti-k_{T} algorithm with a distance parameter of 0.8 and are required to have transverse momentum greater than 550GeV and pseudorapidity |η^{jet}|<1.6. Two-particle correlations among the charged particles within the jets are studied as functions of the particles' azimuthal angle and pseudorapidity separations (Δϕ^{*} and Δη^{*}) in a jet coordinate basis, where particles' η^{*}, ϕ^{*} are defined relative to the direction of the jet. The correlation functions are studied in classes of in-jet charged-particle multiplicity up to N_{ch}^{j}≈100. Fourier harmonics are extracted from long-range azimuthal correlation functions to characterize azimuthal anisotropy for |Δη^{*}|>2. For low-N_{ch}^{j} jets, the long-range elliptic anisotropic harmonic, v_{2}^{*}, is observed to decrease with N_{ch}^{j}. This trend is well described by MonteCarlo event generators. However, a rising trend for v_{2}^{*} emerges at N_{ch}^{j}≳80, hinting at a possible onset of collective behavior, which is not reproduced by the models tested. This observation yields new insights into the dynamics of jet evolution in the vacuum.

Multiturn simulation of radiative Bhabha scattering in the equivalent photon approximation

In this paper, we present a Monte Carlo event generator for radiative Bhabha scattering, based on the method of equivalent photons, and optimized for multiturn tracking simulations in the framework. We demonstrate the accuracy of the event generator with successful benchmarks of the luminosity and event cross section against currently existing alternative tools. Consequently, we use it to investigate existing estimates of the radiative Bhabha scattering beam lifetime at the FCC-ee, as well as the impact of radiative Bhabha scattering on the beam dynamics, using full nonlinear tracking models. Published by the American Physical Society 2024

Observation of quantum entanglement with top quarks at the ATLAS detector

Entanglement is a key feature of quantum mechanics1–3, with applications in fields such as metrology, cryptography, quantum information and quantum computation4–8. It has been observed in a wide variety of systems and length scales, ranging from the microscopic9–13 to the macroscopic14–16. However, entanglement remains largely unexplored at the highest accessible energy scales. Here we report the highest-energy observation of entanglement, in top–antitop quark events produced at the Large Hadron Collider, using a proton–proton collision dataset with a centre-of-mass energy of √s = 13 TeV and an integrated luminosity of 140 inverse femtobarns (fb)−1 recorded with the ATLAS experiment. Spin entanglement is detected from the measurement of a single observable D, inferred from the angle between the charged leptons in their parent top- and antitop-quark rest frames. The observable is measured in a narrow interval around the top–antitop quark production threshold, at which the entanglement detection is expected to be significant. It is reported in a fiducial phase space defined with stable particles to minimize the uncertainties that stem from the limitations of the Monte Carlo event generators and the parton shower model in modelling top-quark pair production. The entanglement marker is measured to be D = −0.537 ± 0.002 (stat.) ± 0.019 (syst.) for 340GeV<mtt¯<380GeV\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$340\\,{\\rm{GeV}} < {m}_{t\\bar{t}} < 380\\,{\\rm{GeV}}$$\\end{document}. The observed result is more than five standard deviations from a scenario without entanglement and hence constitutes the first observation of entanglement in a pair of quarks and the highest-energy observation of entanglement so far.

LO, NLO, and NNLO parton distributions for LHC event generators

We present NNPDF4.0MC, a variant of the NNPDF4.0 set of parton distributions (PDFs) at LO, NLO and NNLO, with and without inclusion of the photon PDF, suitable for use with Monte Carlo (MC) event generators, which require PDFs to satisfy additional constraints in comparison to standard PDF sets. These requirements include PDF positivity down to a low scale Q ~ 1 GeV, smooth extrapolation in the very small and large x regions, and numerically stable results even in extreme regions of phase space for all PDFs. We compare the NNPDF4.0MC PDFs to their baseline NNPDF4.0 counterparts, and to the NNPDF2.3LO set entering the Monash tune of the Pythia8 event generator. We briefly assess the phenomenological impact of these PDFs on the cross-sections for hard and soft QCD processes at the LHC.

Resolving the ultracollinear paradox with effective field theory

Naive intuition about scale decoupling breaks down in the presence of fermion masses. Kinematic enhancements can greatly extend the range where one needs to keep a finite mass in calculations to obtain a correct result at even the O(1) level. Treating a light fermion as massive though, leads to a known but somewhat obscure paradox, a seeming leading-order sensitivity to an arbitrarily small mass. We show how a proper formulation in effective field theory not only resolves the physical conundrum but also answers the very practical question of when fermion masses are required. This has important implications for the development of shower Monte Carlo event generators above the weak scale. Published by the American Physical Society 2024

Measurement of the 1-jettiness event shape observable in deep-inelastic electron-proton scattering at HERA

The H1 Collaboration reports the first measurement of the 1-jettiness event shape observable τ1b\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\ au _1^b$$\\end{document} in neutral-current deep-inelastic electron-proton scattering (DIS). The observable τ1b\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\ au _1^b$$\\end{document} is equivalent to a thrust observable defined in the Breit frame. The data sample was collected at the HERA ep collider in the years 2003–2007 with center-of-mass energy of s=319GeV\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\sqrt{s}=319\\,\ extrm{GeV} $$\\end{document}, corresponding to an integrated luminosity of 351.1 pb-1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\ extrm{pb}^{-1}$$\\end{document}. Triple differential cross sections are provided as a function of τ1b\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\ au _1^b$$\\end{document}, event virtuality Q2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$Q^{2}$$\\end{document}, and inelasticity y, in the kinematic region Q2>150GeV2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$Q^{2} >150\\,\ extrm{GeV}^2 $$\\end{document}. Single differential cross sections are provided as a function of τ1b\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\ au _1^b$$\\end{document} in a limited kinematic range. Double differential cross sections are measured, in contrast, integrated over τ1b\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\ au _1^b$$\\end{document} and represent the inclusive neutral-current DIS cross section measured as a function of Q2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$Q^{2}$$\\end{document} and y. The data are compared to a variety of predictions and include long-standing and more recent Monte Carlo event generators, predictions in fixed-order perturbative QCD where calculations up to O(αs3)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\mathcal {O}(\\alpha _\ extrm{s} ^3)$$\\end{document} are available for τ1b\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\ au _1^b$$\\end{document} or inclusive DIS, and resummed predictions at next-to-leading logarithmic accuracy matched to fixed order predictions at O(αs2)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\mathcal {O}(\\alpha _\ extrm{s} ^2)$$\\end{document}. These comparisons reveal sensitivity of the 1-jettiness observable to QCD parton shower and resummation effects, as well as the modeling of hadronization and fragmentation. Within their range of validity, the fixed-order predictions provide a good description of the data. Monte Carlo event generators are predictive over the full measured range and hence their underlying models and parameters can be constrained by comparing to the presented data.

Measurement of the 12C(e,e′)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\mathrm {{}^{12}C}(e,e')$$\\end{document} Cross Sections at Q2=0.8GeV2/c2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$Q^2=0.8\\,\ extrm{GeV}^2/c^2$$\\end{document}

We present the findings of a study based on a new inelastic electron-scattering experiment on the 12\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${}^{12}$$\\end{document}C nucleus focusing on the kinematic region of Q2=0.8GeV2/c2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$Q^2=0.8\\,\ extrm{GeV}^2/{c}^2$$\\end{document}. The measured cross section is sensitive to the transverse response function and provides a stringent test of theoretical models, as well as of the theoretical assumptions made in Monte-Carlo event-generator codes developed for the interpretation of neutrino-nucleus experiments, such as DUNE and HyperK. We find that modern generators such as GENIE and GiBUU reproduce our new experimental data within 10%\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\%$$\\end{document}.

Observation and differential cross section measurement of neutral current DIS events with an empty hemisphere in the Breit frame

The Breit frame provides a natural frame to analyze lepton–proton scattering events. In this reference frame, the parton model hard interactions between a quark and an exchanged boson defines the coordinate system such that the struck quark is back-scattered along the virtual photon momentum direction. In Quantum Chromodynamics (QCD), higher order perturbative or non-perturbative effects can change this picture drastically. As Bjorken-x decreases below one half, a rather peculiar event signature is predicted with increasing probability, where no radiation is present in one of the two Breit-frame hemispheres and all emissions are to be found in the other hemisphere. At higher orders in αs\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\alpha _{s}$$\\end{document} or in the presence of soft QCD effects, predictions of the rate of these events are far from trivial, and that motivates measurements with real data. We report on the first observation of the empty current hemisphere events in electron–proton collisions at the HERA collider using data recorded with the H1 detector at a center-of-mass energy of 319 GeV. The fraction of inclusive neutral-current DIS events with an empty hemisphere is found to be 0.0112±3.9%stat±4.5%syst±1.6%mod\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$0.0112 \\pm 3.9\\%_\ ext {stat} \\pm 4.5\\%_{\ ext {syst}} \\pm 1.6\\%_{\ ext {mod}}$$\\end{document} in the selected kinematic region of 150<Q2<1500GeV2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$150<Q^{2} <1500\\, \ extrm{GeV}^2 $$\\end{document} and inelasticity 0.14<y<0.7\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$0.14<y<0.7$$\\end{document}. The data sample corresponds to an integrated luminosity of 351.1 pb-1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$${}^{-1}$$\\end{document}, sufficient to enable differential cross section measurements of these events. The results show an enhanced discriminating power at lower Bjorken-x among different Monte Carlo event generator predictions.

Production of hidden-heavy and double-heavy hadronic molecules at the Z factory of CEPC

With a clean environment and high collision energy, the Circular Electron Positron Collider (CEPC) would be an excellent facility for heavy flavor physics. Using the Monte Carlo event generator ythia, we simulate the production of the charmed (bottom) hadron pairs in the electron-positron collisions at the Z factory of CEPC, and the inclusive production rates for typical candidates of the hidden/double-charm and hidden/double-bottom S-wave hadronic molecules are estimated at an order-of-magnitude level with the final state interactions after the hadron pair production. The predicted cross sections for the hidden-charm meson-meson molecules X(3872) and Zc(3900) are at pb level, which are about two to three orders of magnitude larger than the production cross sections for the double-charm meson-meson molecules Tcc and Tcc*, as the double-charmed ones require the production of two pairs of cc¯ from the Z boson decay. The production cross sections for the hidden-charm pentaquark states Pc and Pcs as meson-baryon molecules are a few to tens of fb, which are about one magnitude larger than those of the possible hidden-charm baryon-antibaryon and double-charm meson-baryon molecules. In the bottom sector, the production cross sections for the Zb states as B(*)B¯* molecules are about tens to hundreds of fb, indicating 106–107 events from a two-year operation of CEPC, and the expected events from the double-bottom molecules are about 2–5 orders of magnitude smaller than the Zb states. Our results shows great prospects of probing heavy exotic hadrons at CEPC. Published by the American Physical Society 2024

Isolating perturbative QCD splittings in heavy-ion collisions

We define a new strategy to scan jet substructure in heavy-ion collisions. The scope is multifold: (i) test the dominance of vacuum jet dynamics at early times, (ii) capture the transition from coherent to incoherent jet energy loss, and (iii) study elastic scatterings in the medium, which are either hard and perturbative or soft and responsible for jet thermalization. To achieve that, we analyze the angular distribution of the hardest splitting, θhard, above a transverse momentum scale, ktmin, in high-pt jets. Sufficiently high values of ktmin target the regime in which the observable is uniquely determined by vacuumlike splittings and energy loss, leaving the jet substructure unmodified compared to proton-proton collisions. Decreasing ktmin enhances the sensitivity to the relation between energy loss and the intrajet structure and, in particular, to observe signatures of color decoherence at small angles. At wider angles it also becomes sensitive to hard elastic scatterings with the medium and, therefore, the perturbative regime of medium response. Choosing ktmin≈0 leads to order one effects of nonperturbative origin such as hadronization and, potentially, soft scatterings responsible for jet thermalization. We perform a comprehensive analysis of this observable with three state-of-the-art jet-quenching Monte Carlo event generators. Our study paves the way for defining jet observables in heavy-ion collisions dominated by perturbative QCD and thus calculable from first principles. Published by the American Physical Society 2024

Exploring the time axis within medium-modified jets

In this manuscript, we illustrate how to use the newly proposed τ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\ au $$\\end{document} re-clustering algorithm to select jets with different degrees of quenching without biasing their initial transverse momentum spectrum. Our study is based on Z+jet simulated events using the JEWEL Monte Carlo event generator to account for jet quenching effects. We apply the τ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\ au $$\\end{document} re-clustering algorithm to extract a proxy for a time axis (formation time) within the evolving medium. This information allows us to label jets according to their fragmentation pattern and select populations with enhanced sensitivity to quenching effects. Our results illustrate the potential of jets as precision tools for QGP tomography. Further, we show that the discussed method minimizes the biases stemming from pT\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$p_{T}$$\\end{document}-, ΔR\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\Delta R$$\\end{document}- or mass-based jet selection.

Comparative analysis of jet and underlying event properties across various models as a function of charged particle multiplicity at 7 TeV* *Supported by Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2024R106), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia

In this study, a comprehensive analysis of jets and underlying events as a function of charged particle multiplicity in proton-proton (pp) collisions at a center-of-mass energy of TeV is conducted. Various Monte Carlo (MC) event generators, including Pythia8.308, EPOS1.99, EPOSLHC, EPOS4 , and EPOS4 , are employed to predict particle production. The predictions from these models are compared with experimental data from the CMS collaboration. The charged particles are categorized into those associated with underlying events and those linked to jets, and the analysis is restricted to charged particles with and GeV/c. By comparing the MC predictions with CMS data, we find that EPOS , EPOSLHC, and Pythia8 consistently reproduce the experimental results for all charged particles, underlying events, intrajets, and leading charged particles. For charged jet rates with GeV/c, EPOS4 and Pythia8 perform exceptionally well. In the case of charged jet rates with GeV/c, EPOSLHC reproduces satisfactorily good results, whereas EPOS4 exhibits good agreement with the data at higher charged particle multiplicities compared to the other models. This can be attributed to the conversion of energy into flow when "Hydro=on," leading to an increase in multiplicity. The EPOSLHC model describes the data better owing to the new collective flow effects, correlated flow treatment, and parameterization compared to EPOS1.99. However, the examination of the jet spectrum and normalized charged density reveals that EPOS4 , EPOS4 , and EPOSLHC exhibit good agreement with the experimental results, whereas Pythia8 and EPOS1.99 do not perform as well owing to the lack of correlated flow treatment.

Pythia8 MC Tuning Validation Using the Professor2 Package

This article presents a method for multi-parameter, simultaneous tuning of the Monte Carlo event generator. It is validated on the Pythia8 Monte Carlo event generator widely used in High Energy Physics (HEP). The obtained results show that the method can be used to constrain the free parameters of phenomenological models while properly taking into account the correlations existing among the parameters.

Event generators for high-energy physics experiments

We provide an overview of the status of Monte-Carlo event generators for high-energy particle physics. Guided by the experimental needs and requirements, we highlight areas of active development, and opportunities for future improvements. Particular emphasis is given to physics models and algorithms that are employed across a variety of experiments. These common themes in event generator development lead to a more comprehensive understanding of physics at the highest energies and intensities, and allow models to be tested against a wealth of data that have been accumulated over the past decades. A cohesive approach to event generator development will allow these models to be further improved and systematic uncertainties to be reduced, directly contributing to future experimental success. Event generators are part of a much larger ecosystem of computational tools. They typically involve a number of unknown model parameters that must be tuned to experimental data, while maintaining the integrity of the underlying physics models. Making both these data, and the analyses with which they have been obtained accessible to future users is an essential aspect of open science and data preservation. It ensures the consistency of physics models across a variety of experiments.

Lepton pair production in muon-nucleus scattering

The MUonE experiment aims at providing a novel determination of the leading hadronic contribution to the muon anomalous magnetic moment through the study of elastic muon-electron scattering. Since the initial-state electrons are bound in a low-Z atomic target, the interaction between the incoming muons and the nuclei is expected to be the main source of experimental background. In this article, we study the production of a real lepton pair from the muon-nucleus scattering, discussing its numerical impact in the MUonE kinematic configuration. The process is described as a scattering of a muon in an external Coulomb field with the addition of a form factor to describe the nuclear charge distribution. The calculation is implemented in the fully differential Monte Carlo event generator Mesmer, without introducing any approximation on the angular variables.

Light-flavor particle production in high-multiplicity pp collisions at s\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\sqrt{\ extrm{s}} $$\\end{document}= 13 TeV as a function of transverse spherocity

Results on the transverse spherocity dependence of light-flavor particle production (π, K, p, ϕ, K*0, KS0\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ {\ extrm{K}}_{\ extrm{S}}^0 $$\\end{document}, Λ, Ξ) at midrapidity in high-multiplicity pp collisions at s\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\sqrt{s} $$\\end{document} = 13 TeV were obtained with the ALICE apparatus. The transverse spherocity estimator SOpT=1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\left({S}_{\ extrm{O}}^{p_{\ extrm{T}}=1}\\right) $$\\end{document} categorizes events by their azimuthal topology. Utilizing narrow selections on SOpT=1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ {S}_{\ extrm{O}}^{p_{\ extrm{T}}=1} $$\\end{document}, it is possible to contrast particle production in collisions dominated by many soft initial interactions with that observed in collisions dominated by one or more hard scatterings. Results are reported for two multiplicity estimators covering different pseudorapidity regions. The SOpT=1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ {S}_{\ extrm{O}}^{p_{\ extrm{T}}=1} $$\\end{document} estimator is found to effectively constrain the hardness of the events when the midrapidity (|η| < 0.8) estimator is used.The production rates of strange particles are found to be slightly higher for soft isotropic topologies, and severely suppressed in hard jet-like topologies. These effects are more pronounced for hadrons with larger mass and strangeness content, and observed when the topological selection is done within a narrow multiplicity interval. This demonstrates that an important aspect of the universal scaling of strangeness enhancement with final-state multiplicity is that high-multiplicity collisions are dominated by soft, isotropic processes. On the contrary, strangeness production in events with jet-like processes is significantly reduced.The results presented in this article are compared with several QCD-inspired Monte Carlo event generators. Models that incorporate a two-component phenomenology, either through mechanisms accounting for string density, or thermal production, are able to describe the observed strangeness enhancement as a function of SOpT=1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ {S}_{\ extrm{O}}^{p_{\ extrm{T}}=1} $$\\end{document}.