New error control methods for integration of trajectories


Geant4 and GeantV use Runge-Kutta methods to integrate the motion of charged particles in a non-uniform electromagnetic field. Methods must provide good integration accuracy for the integration and to cost a minimum of computation time. Integration is used to identify the intersection point between the curved track and the volume boundaries. Due to the large number of steps and the cost of the evaluations of the field, the integration and intersection are a performance critical part of detector simulation. Recent work has introduced new RK methods which reduce the number of field evaluations required, and has the potential to decreased the computation time.

Task ideas

Expected results

Working implementation of:


This project requires prior exposure to Numerical Analysis and familiarity with either C++, C or Java programming. Exposure to either numerical methods for solving Ordinary Differential equations (ODEs) or tools for analysing data such as ROOT or R will be valuable. Both programming skill and knowledge of numerical methods for ODEs will be improved by undertaking this project.


John Apostolakis

Web site

Source code

Corresponding Project

Participating Organizations