HepLib: A C++ library for high energy physics (version 1.1)

• Published in: Computer physics communications Vol. 285; p. 108631
• Main Authors: Feng, Feng, Tang, Shan-Rong, Gao, Ya-Di
• Format: Journal Article
• Language: English
• Published: Elsevier B.V 01.04.2023
• Subjects: IBP; Loop corrections; Sector decomposition; Loop corrections; IBP; Sector decomposition
• ISSN: 0010-4655; 1879-2944
• DOI: 10.1016/j.cpc.2022.108631

HepLib is a C++ Library for computations in High Energy Physics, it builds on top of other well-established libraries or programs, including GiNaC, Fermat, Form, Fire, etc., its first version has been released in Comput. Phys. Commun. 265 (2021) 107982. Here we provide a minor upgraded version with parts of obligated aspects removed and many new features added, the main feature in this new version is to use Flint for the simplification on multi-variate polynomial, so one can replace the process-based parallelism by the thread-based parallelism for the related polynomial evaluations to achieve a higher performance. Program title:HepLib 1.1 CPC Library link to program files:https://doi.org/10.17632/jspjxy3z66.2 Developer's repository link:https://github.com/HepLib/src Licensing provisions: : GPLv3 Programming language:C++ Journal reference of previous version: Comput. Phys. Commun. 265 (2021) 107982 Does the new version supersede the previous version?: Yes External libraries and programs: : GiNaC [2], Qgraf [3], Flint [4], Qhull [5], Form [6], Fermat [7], Fire [8] Nature of problem: There are lots of independent packages or programs in high energy physics, which are developed in different programming languages, here we provide an integrated interface in C++ language, to generate Feynman diagrams/amplitudes, perform algebraic simplifications on Dirac/Color matrixes and reduce scalar integrals to master integrals. The evaluation of master integrals is the core part for the high-order calculations, one has to resolve ultraviolet/infrared divergence within dimensional regularization and the singularities inside the integration domain while performing multiple dimensional numerical integrations. Solution method:HepLib uses Qgraf to generate the Feynman diagrams/aptitudes, Form to perform the related computations on Dirac/Color matrix, Fire for the IBP reduction. A C++ implementation of sector decomposition method is used to extract both ultraviolet and infrared singularity of the integral, the contour deformation is adopted to avoid the singularity inside the integration domain, a parallelized numerical integration is chosen to achieve high efficiency or performance. Reasons for the new version: Fix bugs from GiNaC and implement new features. Summary of revisions: First, we have improved the installation procedure, now one can install the whole package with a single script install.sh, and furthermore one can install all libraries into a distinct directory, so the libraries will not affect the others which are installed on the user's computer. We have upgraded the GiNaC to the latest version with a few bugs fixed in the official site and removed the dependence on the libraries MinUits [9] and Cuba [10], which are used in performing sector decomposition in the previous version. To get high performance on the evaluation of multi-variate polynomials, we introduce Flint in the new version, and one can convert the expression from GiNaC to Flint, and vice versa. Finally, we also made adaptations to the codes and examples to match the new version. [1]F. Feng, Y.F. Xie, Q.C. Zhou, S.R. Tang, Comput. Phys. Commun. 265 (2021) 107982, https://doi.org/10.1016/j.cpc.2021.107982.[2]https://www.ginac.de.[3]http://cfif.ist.utl.pt/~paulo/qgraf.html.[4]https://flintlib.org.[5]http://qhull.org.[6]https://github.com/vermaseren/form.[7]http://home.bway.net/lewis.[8]https://bitbucket.org/feynmanintegrals.[9]https://root.cern.ch/doc/master/md_math_minuit2_doc_Minuit2.html.[10]https://feynarts.de/cuba/.