Parton Shower Monte Carlo Programs Research Articles

An event generator for Lepton-Hadron deep inelastic scattering at NLO+PS with POWHEG including mass effects

We present a generator for lepton nucleon collisions in the DIS regime, focusing in particular on processes with a massive lepton and/or a massive quark in the final state. We have built a full code matching NLO QCD corrections to parton shower Monte Carlo programs in the POWHEG-BOX framework. Our code can be used to compute NLO+PS accurate fully differential predictions for neutral current and charged current processes, including processes with an incoming tau neutrino, and/or including charm quarks in the final state. We also made comparisons with available data and predictions for the new neutrino experiments at CERN.

Fitting a deep generative hadronization model

Hadronization is a critical step in the simulation of high-energy particle and nuclear physics experiments. As there is no first principles understanding of this process, physically-inspired hadronization models have a large number of parameters that are fit to data. Deep generative models are a natural replacement for classical techniques, since they are more flexible and may be able to improve the overall precision. Proof of principle studies have shown how to use neural networks to emulate specific hadronization when trained using the inputs and outputs of classical methods. However, these approaches will not work with data, where we do not have a matching between observed hadrons and partons. In this paper, we develop a protocol for fitting a deep generative hadronization model in a realistic setting, where we only have access to a set of hadrons in data. Our approach uses a variation of a Generative Adversarial Network with a permutation invariant discriminator. We find that this setup is able to match the hadronization model in Herwig with multiple sets of parameters. This work represents a significant step forward in a longer term program to develop, train, and integrate machine learning-based hadronization models into parton shower Monte Carlo programs.

Exploring the universality of hadronic jet classification

The modeling of jet substructure significantly differs between Parton Shower Monte Carlo (PSMC) programs. Despite this, we observe that machine learning classifiers trained on different PSMCs learn nearly the same function. This means that when these classifiers are applied to the same PSMC for testing, they result in nearly the same performance. This classifier universality indicates that a machine learning model trained on one simulation and tested on another simulation (or data) will likely be optimal. Our observations are based on detailed studies of shallow and deep neural networks applied to simulated Lorentz boosted Higgs jet tagging at the LHC.

HELAC-Onia 2.0: An upgraded matrix-element and event generator for heavy quarkonium physics

We present an upgraded version (denoted as version 2.0) of the program HELAC-Onia for the automated computation of heavy-quarkonium helicity amplitudes within non-relativistic QCD framework. The new code has been designed to include many new and useful features for practical phenomenological simulations. It is designed for job submissions under cluster environment for parallel computations via Python scripts. We have interfaced HELAC-Onia to the parton shower Monte Carlo programs Pythia 8 and QEDPS to take into account the parton-shower effects. Moreover, the decay module guarantees that the program can perform the spin-entangled (cascade-)decay of heavy quarkonium after its generation. We have also implemented a reweighting method to automatically estimate the uncertainties from renormalization and/or factorization scales as well as parton-distribution functions to weighted or unweighted events. A further update is the possibility to generate one-dimensional or two-dimensional plots encoded in the analysis files on the fly. Some dedicated examples are given at the end of the writeup. Program summaryProgram title: HELAC-Onia 2.0Catalogue identifier: AEPR_v2_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEPR_v2_0.htmlProgram obtainable from: CPC Program Library, Queen’s University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 1513282No. of bytes in distributed program, including test data, etc.: 17036140Distribution format: tar.gzProgramming language: Python, Fortran 77, Fortran 90, C++.Operating system: Unix-like platform.Classification: 4.4, 11.1, 11.2, 11.5.Catalogue identifier of previous version: AEPR_v1_0Journal reference of previous version: Comput. Phys. Comm. 184(2013)2562External routines: LHAPDFDoes the new version supersede the previous version?: YesNature of problem: Heavy quarkonium production processes provide an important way to investigate QCD in its poorly known non-perturbative regime. Its production mechanism has attracted extensive interests from the high-energy physics community in decades. The qualitative and quantitative description of heavy-quarkonium production requires complex perturbative computations for high-multiplicity processes in the framework of the well established non-relativistic effective theory, NRQCD, and reliable Monte Carlo simulations to reproduce the collider environment.Solution method: Based on a recursion relation, the program is able to calculate the helicity amplitudes of the high-multiplicity heavy- quarkonium-production processes. Several modules are also designed for dedicated simulations: 1.The code has been interfaced with the parton shower Monte Carlo programs;2.A decay module to let heavy quarkonia decay with correct spin-correlations has been implemented;3.The code estimates the theoretical uncertainties and analyzes the generated events on the fly;4.The code is compliant with multi-threading/multi-core usage or cluster processors.Reasons for new version: Improved and expanded functionalities.Summary of revisions: Many new features were added and several important bugs were fixed. The new features extend the range of the physical applications. With the new interface, it also helps to improve the user-friendliness of the program.Running time: It depends on the process to be calculated and the required accuracy.

Squark production and decay matched with parton showers at NLO

Extending previous work on the predictions for the production of supersymmetric (SUSY) particles at the LHC, we present the fully differential calculation of the next-to-leading order (NLO) SUSY-QCD corrections to the production of squark and squark–antisquark pairs of the first two generations. The NLO cross sections are combined with the subsequent decay of the final state (anti)squarks into the lightest neutralino and (anti)quark at NLO SUSY-QCD. No assumptions on the squark masses are made, and the various subchannels are taken into account independently. In order to obtain realistic predictions for differential distributions the fixed-order calculations have to be combined with parton showers. Making use of the Powheg method we have implemented our results in the Powheg-Box framework and interfaced the NLO calculation with the parton shower Monte Carlo programs Pythia6 and Herwig++. The code is publicly available and can be downloaded from the Powheg-Box webpage. The impact of the NLO corrections on the differential distributions is studied and parton shower effects are investigated for different benchmark scenarios.

Electroweak W+W− jj prodution at NLO in QCD matched with parton shower in the POWHEG-BOX

We present an implementation of electroweak W+W-jj production at hadron colliders in the POWHEG framework, a method that allows the interfacing of a next-to-leading order QCD calculation with parton shower Monte Carlo programs. We provide results for both, fully and semi-leptonic decay modes of the weak bosons, taking resonant and non-resonant contributions and spin correlations of the final-state particles into account. To illustrate the versatility of our implementation, we provide phenomenological results for two representative scenarios with a light and with a heavy Higgs boson, respectively, and in a kinematic regime of highly boosted gauge bosons. The impact of the parton shower is found to depend on the setup and the observable under investigation. In particular, distributions related to a central-jet veto are more sensitive to these effects. Therefore the impact of radiation by the parton shower on next-to-leading order predictions should be assessed carefully on a case-by-case basis.

Slepton pair production in the POWHEG BOX

We present an implementation for slepton pair production at hadron colliders in the POWHEG BOX, a framework for combining next-to-leading order QCD calculations with parton-shower Monte-Carlo programs. Our code provides a SUSY Les Houches Accord interface for setting the supersymmetric input parameters. Decays of the sleptons and parton-shower effects are simulated with PYTHIA. Focussing on a representative point in the supersymmetric parameter space we show results for kinematic distributions that can be observed experimentally. While next-to-leading order QCD corrections are sizable for all distributions, the parton shower affects the color-neutral particles only marginally. Pronounced parton-shower effects are found for jet distributions.

Next-to-leading order QCD corrections to electroweak Zjj production in the POWHEG BOX

We present an implementation of electroweak Z-boson production in association with two jets at hadron colliders in the POWHEG framework, a method that allows the interfacing of NLO-QCD calculations with parton-shower Monte Carlo programs. We focus on the leptonic decays of the weak gauge boson, and take photonic and non-resonant contributions to the matrix elements fully into account. We provide results for observables of particular importance for the suppression of QCD backgrounds to vector-boson fusion processes by means of central-jet-veto techniques. While parton-shower effects are small for most observables associated with the two hardest jets, they can be more pronounced for distributions that are employed in central-jet-veto studies.

Single-top t-channel hadroproduction in the four-flavour scheme with POWHEG and aMC@NLO

Abstract We present results for the QCD next-to-leading order (NLO) calculation of single-top t-channel production in the 4-flavour scheme, interfaced to Parton Shower (PS) Monte Carlo programs according to the POWHEG and MC@NLO methods. Comparisons between the two methods, as well as with the corresponding process in the 5-flavour scheme are presented. For the first time results for typical kinematic distributions of the spectator-b jet are presented in an NLO + PS approach.

Hadronic top-quark pair-production with one jet and parton showering

Abstract We present a calculation of heavy-flavor production in hadronic collisions in association with one jet matched to parton shower Monte Carlo programs at next-to-leading order in perturbative QCD. Top-quark decays are included and spin correlations in the decay products are taken into account. The calculation builds on existing results for the radiative corrections to heavy-quark plus one jet production and uses the POWHEG BOX for the interface to the parton shower programs PYTHIA or HERWIG. A broad phenomenological study for the Large Hadron Collider and the Tevatron is presented. In particular we study — as one important sample application — the impact of the parton shower on the top-quark charge asymmetry.

Next-to-leading order QCD effects and the top quark mass measurements at the LHC

It is anticipated that a number of techniques to measure the top quark mass at the LHC will yield m_top with uncertainties of about 0.5-1 percent. These uncertainties are mostly theoretical; they are usually estimated using parton shower Monte Carlo programs whose reliability at this level of precision is difficult to assess. The goal of this paper is to contrast those estimates with the results of NLO QCD computations for a few observables, often discussed in the context of high-precision top quark mass measurements at the LHC. In particular, we study the NLO QCD corrections to the invariant mass distribution of a charged lepton and a B-meson in lepton+jets channels. In the dilepton channel we investigate the invariant mass distribution of a charged lepton and a b-jet, the average energy of the two leptons and the average energy of the b-jets from top decays.

Angular multiplicity fluctuations in hadronic Z decays and comparison to QCD models and analytical calculations

Abstract Local multiplicity fluctuations in angular phase space intervals are studied using factorial moments measured in hadronic events at s ≃91.2 GeV , which were collected by the L3 detector at LEP. Parton shower Monte Carlo programs agree well with the data. On the other hand, first-order QCD calculations in the Double Leading Log Approximation and the Modified Leading Log Approximation are found to deviate significantly from the data.

Measurement of inclusive η production in hadronic decays of the Z 0

We present a study of the inclusive η production based on 300 000 hadronic Z 0 decays. The measured inclusive momentum distribution can be reproduced by parton shower Monte Carlo programs and also by an analytical QCD calculation. Comparing our results with low energy e +e − data, we find that QCD describes both the shape and the energy evolution of the η spectrum. The comparison of η production rates in quark- and gluon-enriched jet samples does not show statistically significant evidence for more abundant production of η mesons in gluon fragmentation.

Measurement of the inclusive production of neutral pions and charged particles on the Z 0 resonance

We present a study of the inclusive production of neutral pions and charged particles from 112 000 hadronic Z 0 decays. The measured inclusive momentum distributions can be reproduced by parton shower Monte Carlo programs and also by an analytical QCD calculation. Comparing our results to e +e − data between √ s = 9 and 91 GeV, we findfind that the evolution of the spectra with center of mass energy is consistent with the QCD predictions.