GQLink: an implementation of Quantized State Systems (QSS) methods in Geant4

• Published in: Journal of physics. Conference series Vol. 1085; no. 5; pp. 52015 - 52020
• Main Authors: Santi, Lucio, Bergero, Federico, Jun, Soon Yung, Genser, Krzysztof, Elvira, Daniel, Castro, Rodrigo
• Format: Journal Article
• Language: English
• Published: Bristol IOP Publishing 01.09.2018
• Subjects: Accuracy; Charged particles; Differential equations; Equations of motion; Iterative algorithms; Iterative methods; Mathematical analysis; Numerical analysis; Numerical methods; Ordinary differential equations; Physics; PHYSICS OF ELEMENTARY PARTICLES AND FIELDS; Polynomials; Simulation
• ISSN: 1742-6588; 1742-6596; 1742-6596
• DOI: 10.1088/1742-6596/1085/5/052015

Simulations in high energy physics (HEP) often require the numerical solution of ordinary differential equations (ODE) to determine the trajectories of charged particles in a magnetic field when particles move throughout detector volumes. Each crossing of a volume interrupts the underlying numerical method that solves the equations of motion, triggering iterative algorithms to estimate the intersection point within a given accuracy. The computational cost of this procedure can grow significantly depending on the application at hand. Quantized State System (QSS) is a recent family of discrete-event driven numerical methods exhibiting attractive features for this type of problems, such as native dense output (sequences of polynomial segments updated only by accuracy-driven events) and lightweight detection and handling of volume crossings. In this work we present GQLink, a proof-of-concept integration of QSS with the Geant4 simulation toolkit which stands as an interface for co-simulation that orchestrates robustly and transparently the interaction between the QSS simulation engine and aspects such as geometry definition and physics processes controlled by Geant4. We validate the accuracy and study the performance of the method in simple geometries (subject to intense volume crossing activity) and then in a realistic HEP application using a full CMS detector configuration.