TY - JOUR
T1 - Fast computation of MadGraph amplitudes on graphics processing unit (GPU)
AU - Hagiwara, K.
AU - Kanzaki, J.
AU - Li, Q.
AU - Okamura, N.
AU - Stelzer, T.
N1 - Funding Information:
This work is supported by the Grant-in-Aid for Scientific Research from the Japan Society for the Promotion of Science (Nos. 20340064 and 22740155), in part by the National Natural Science Foundation of China under Grants No. 11205008 as well as a grant from the National Science Foundation (NSF PHY-1068326). We wish to thank Dave Rainwater for inspiring us with the idea of using GPU for fast LHC physics simulation and his initial contribution to the project.
PY - 2013
Y1 - 2013
N2 - Continuing our previous studies on QED and QCD processes, we use the graphics processing unit (GPU) for fast calculations of helicity amplitudes for general Standard Model (SM) processes. Additional HEGET codes to handle all SM interactions are introduced, as well as the program MG2CUDA that converts arbitrary MadGraph generated HELAS amplitudes (FORTRAN) into HEGET codes in CUDA. We test all the codes by comparing amplitudes and cross sections for multi-jet processes at the LHC associated with production of single and double weak bosons, a top-quark pair, Higgs boson plus a weak boson or a top-quark pair, and multiple Higgs bosons via weak-boson fusion, where all the heavy particles are allowed to decay into light quarks and leptons with full spin correlations. All the helicity amplitudes computed by HEGET are found to agree with those computed by HELAS within the expected numerical accuracy, and the cross sections obtained by gBASES, a GPU version of the Monte Carlo integration program, agree with those obtained by BASES (FORTRAN), as well as those obtained by MadGraph. The performance of GPU was over a factor of 10 faster than CPU for all processes except those with the highest number of jets.
AB - Continuing our previous studies on QED and QCD processes, we use the graphics processing unit (GPU) for fast calculations of helicity amplitudes for general Standard Model (SM) processes. Additional HEGET codes to handle all SM interactions are introduced, as well as the program MG2CUDA that converts arbitrary MadGraph generated HELAS amplitudes (FORTRAN) into HEGET codes in CUDA. We test all the codes by comparing amplitudes and cross sections for multi-jet processes at the LHC associated with production of single and double weak bosons, a top-quark pair, Higgs boson plus a weak boson or a top-quark pair, and multiple Higgs bosons via weak-boson fusion, where all the heavy particles are allowed to decay into light quarks and leptons with full spin correlations. All the helicity amplitudes computed by HEGET are found to agree with those computed by HELAS within the expected numerical accuracy, and the cross sections obtained by gBASES, a GPU version of the Monte Carlo integration program, agree with those obtained by BASES (FORTRAN), as well as those obtained by MadGraph. The performance of GPU was over a factor of 10 faster than CPU for all processes except those with the highest number of jets.
UR - http://www.scopus.com/inward/record.url?scp=84886742617&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=84886742617&partnerID=8YFLogxK
U2 - 10.1140/epjc/s10052-013-2608-2
DO - 10.1140/epjc/s10052-013-2608-2
M3 - Article
AN - SCOPUS:84886742617
VL - 73
SP - 1
EP - 39
JO - European Physical Journal C
JF - European Physical Journal C
SN - 1434-6044
IS - 11
M1 - 2608
ER -