HELAC-Onia: An automatic matrix element generator for heavy quarkonium physics

Abstract By the virtues of the Dyson–Schwinger equations, we upgrade the published code HELAC to be capable to calculate the heavy quarkonium helicity amplitudes in the framework of NRQCD factorization, which we dub HELAC-Onia . We rewrote the original HELAC to make the new program be able to calculate helicity amplitudes of multi P -wave quarkonium states production at hadron colliders and electron–positron colliders by including new P -wave off-shell currents. Therefore, besides the high efficiencies in computation of multi-leg processes within the Standard Model, HELAC-Onia is also sufficiently numerical stable in dealing with P -wave quarkonia (e.g. h c , b , χ c , b ) and P -wave color-octet intermediate states. To the best of our knowledge, it is a first general-purpose automatic quarkonium matrix elements generator based on recursion relations on the market. Program summary Program title: HELAC-Onia . Catalogue identifier:  AEPR_v1_0 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEPR_v1_0.html Program obtainable from:  CPC Program Library, Queen’s University, Belfast, N. Ireland Licensing provisions:  Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.:  114595 No. of bytes in distributed program, including test data, etc.:  1555550 Distribution format:  tar.gz Programming language:  Fortran 90. Computer:  Any. Operating system:  Windows, Unix. Classification:  4.4, 11.1, 11.2, 11.5. Nature of problem: An important way to explore the law of the nature is to investigate the heavy quarkonium physics at B factories and hadron colliders. However, its production mechanism is still unclear, though NRQCD can explain its decay mechanism in a sufficiently satisfactory manner. The substantial K-factor in heavy quarkonium production processes also implies that the associated production of quarkonium and a relatively large number of particles may play a crucial role in unveiling its production mechanism. Solution method: A labor-saving and efficient way is to make the tedious amplitudes calculation automatic. Based on a recursive algorithm derived from the Dyson–Schwinger equations, the goal of automatic calculation of heavy quarkonium helicity amplitudes in NRQCD has been achieved. Inheriting from the virtues of the recursion relations with the lower computational cost compared to the traditional Feynman-diagram based method, the multi-leg processes (with or without multi-quarkonia up to P -wave states) at colliders are also accessible. Running time: It depends on the process that is to be calculated. However, typically, for all of the tested processes, they take from several minutes to tens of minutes.

[1]  B. Kniehl,et al.  Complete next-to-leading-order corrections to J/psi photoproduction in nonrelativistic quantum chromodynamics. , 2009, Physical review letters.

[2]  K. Chao,et al.  Pair production of heavy quarkonium and B c ( * ) mesons at hadron colliders , 2009, 0903.2250.

[3]  J. Huston,et al.  New generation of parton distributions with uncertainties from global QCD analysis , 2002, hep-ph/0201195.

[4]  Hong-Fei Zhang,et al.  Polarization for prompt J/ψ and ψ(2s) production at the Tevatron and LHC. , 2012, Physical review letters.

[5]  G. Lepage.  A new algorithm for adaptive multidimensional integration , 1978 .

[6]  S. D. Ellis,et al.  A New Monte Carlo Treatment of Multiparticle Phase Space at High-energies , 1986 .

[7]  B. Kniehl,et al.  Probing nonrelativistic QCD factorization in polarized J/ψ photoproduction at next-to-leading order. , 2011, Physical review letters.

[8]  F. Maltoni,et al.  Upsilon production at Fermilab Tevatron and LHC energies. , 2008, Physical review letters.

[9]  W. Giele,et al.  Recursive calculations for processes with n gluons , 1988 .

[10]  H. Shao,et al.  Spin correlations in polarizations of P-wave charmonia $\chi_{cJ}$ and impact on $J/\psi$ polarization , 2012, 1209.4610.

[11]  F. Dyson.  The S Matrix in Quantum Electrodynamics , 1949 .

[12]  F. Maltoni,et al.  Hadroproduction of J/psi and Upsilon in association with a heavy-quark pair , 2007 .

[13]  Alessandro Cafarella,et al.  Helac-Phegas: A generator for all parton level processes , 2007, Comput. Phys. Commun..

[14]  Dimitri Bourilkov,et al.  The Les Houches accord PDFs (LHAPDF) and LHAGLUE , 2005, hep-ph/0508110.

[15]  Costas G. Papadopoulos,et al.  PHEGAS : A phase-space generator for automatic cross-section computation , 2000, hep-ph/0007335.

[16]  B. Kniehl,et al.  Reconciling J/ψ production at HERA, RHIC, Tevatron, and LHC with nonrelativistic QCD factorization at next-to-leading order. , 2010, Physical review letters.

[17]  Chao-Hsi Chang,et al.  BCVEGPY2.0: An upgraded version of the generator BCVEGPY with the addition of hadroproduction of the P-wave Bc states , 2006, Comput. Phys. Commun..

[18]  BCVEGPY: an event generator for hadronic production of the B-c meson , 2003, hep-ph/0309120.

[19]  M. Moretti,et al.  An algorithm to compute Born scattering amplitudes without Feynman graphs , 1995, hep-ph/9507237.

[20]  Hua-Sheng Shao,et al.  J/ψ polarization at hadron colliders in nonrelativistic QCD. , 2012, Physical review letters.

[21]  L. Lönnblad,et al.  Comparative study of various algorithms for the merging of parton showers and matrix elements in hadronic collisions , 2007, 0706.2569.

[22]  J. Kanzaki,et al.  Generic User Process Interface for Event Generators , 2001 .

[23]  G. Bodwin.  Theory of Charmonium Production , 2012, 1208.5506.

[24]  Bodwin,et al.  Rigorous QCD analysis of inclusive annihilation and production of heavy quarkonium. , 1995, Physical review. D, Particles and fields.

[25]  F. Maltoni,et al.  Upsilon production at Fermilab Tevatron and LHC energies. , 2008, Physical review letters.

[26]  M. Mangano,et al.  Matching matrix elements and shower evolution for top-pair production in hadronic collisions , 2006, hep-ph/0611129.

[27]  J. Kuhn,et al.  Rare decays of the Z0 , 1980 .

[28]  E. Berger,et al.  Inelastic Photoproduction of J/psi and Upsilon by Gluons , 1981 .

[29]  B. Kniehl,et al.  J/ψ polarization at the Tevatron and the LHC: nonrelativistic-QCD factorization at the crossroads. , 2012, Physical review letters.

[30]  Aggeliki Kanaki,et al.  HELAC: A package to compute electroweak helicity amplitudes , 2000, hep-ph/0002082.

[31]  F. Maltoni,et al.  Hadroproduction of J / ψ and Υ in association with a heavy-quark pair , 2007 .

[32]  F. Maltoni,et al.  Automatic generation of quarkonium amplitudes in NRQCD , 2007, 0712.2770.

[33]  H. K. Wöhri,et al.  Heavy quarkonium: progress, puzzles, and opportunities , 2010, 1010.5827.

[34]  J Schwinger,et al.  On the Green's Functions of Quantized Fields: I. , 1951, Proceedings of the National Academy of Sciences of the United States of America.

[35]  C. Papadopoulos,et al.  HELAC-PHEGAS: Automatic computation of helicity amplitudes and cross sections , 2000, hep-ph/0012004.

[36]  G. Hooft.  A Planar Diagram Theory for Strong Interactions , 1974 .

[37]  K. Chao,et al.  J/ψ(ψ') production at the Tevatron and LHC at O(α(s)4v4) in nonrelativistic QCD. , 2010, Physical review letters.

[38]  Jian-Xiong Wang.  Progress in FDC project , 2004 .

[39]  M. Greco,et al.  NLO production and decay of quarkonium , 1997, hep-ph/9707223.