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 enviroment for parallel computations via Python scripts. We have interfaced HELACOnia 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 futher update is the possiblity 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.
Programming language:Python, Fortran 77, Fortran 90, C++ Keywords: heavy quarkonium, NRQCD, Monte Carlo simulation Nature of physical problem: Heavy quarkonium production processes provide an important way to investigate QCD in its poorly known non-perturbative regime. Its production mechanism has been 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 effecitive theory, NRQCD, and reliable Monte Carlo simulations to repreduce the collider enviroment.
Method of solution:Based on a recursion relation, the program is able to calculate the helicity ampltiudes of the high-multiplicity heavy-qurkonium-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 compilant with multi-threading/multi-core usage or cluster processors.