A C++ Cherenkov photons simulation in CORSIKA 8

¹ LUPM : Laboratoire Univers et Particules de Montpellier, France
² DALI : Digits, Architectures et Logiciels Informatiques, Perpignan, France
³ LIRMM : Laboratoire d’Informatique, de Robotique et de Microélectronique de Montpellier, France
⁴ LPSC : Laboratoire de Physique Subatomique et de Cosmologie, Grenoble, France

^* e-mail: matthieu.carrere@lupm.in2p3.fr

Published online: 23 August 2021

Abstract

CORSIKA is a standard software for simulations of air showers induced by cosmic rays. It has been developed mainly in Fortran 77 continuously over the last thirty years. It has become very difficult to add new physics features to CORSIKA 7. CORSIKA 8 aims to be the future of the CORSIKA project. It is a framework in C++17 which uses modern concepts in object oriented programming for an efficient modularity and flexibility. The CORSIKA 8 project aims to attain high performance by exploiting techniques such as vectorization, gpu/cpu parallelization, extended use of static polymorphism and the most precise physical models available. In this paper, we focus on the Cherenkov photon propagation module of CORSIKA, which is of particular interest for gamma-ray experiments, like the Cherenkov Telescope Array. First, we present the optimizations that we have applied to the Cherenkov module thanks to the results of detailed profiling using performance counters. Then, we report our preliminary work to develop the Cherenkov Module in the CORSIKA 8 framework. Finally, we will demonstrate the first performance comparison with the current CORSIKA software as well as physics validation.