gSeaGen code by KM3NeT: An efficient tool to propagate muons simulated with CORSIKA

Kód výsledku v IS VaVaI

<a href="https://www.isvavai.cz/riv?ss=detail&h=RIV%2F68407700%3A21670%2F25%3A00384174" target="_blank" >RIV/68407700:21670/25:00384174 - isvavai.cz</a>

Výsledek na webu

<a href="https://doi.org/10.1016/j.cpc.2025.109660" target="_blank" >https://doi.org/10.1016/j.cpc.2025.109660</a>

DOI - Digital Object Identifier

<a href="http://dx.doi.org/10.1016/j.cpc.2025.109660" target="_blank" >10.1016/j.cpc.2025.109660</a>

Jazyk výsledku

angličtina

Název v původním jazyce

gSeaGen code by KM3NeT: An efficient tool to propagate muons simulated with CORSIKA

Popis výsledku v původním jazyce

The KM3NeT Collaboration has tackled a common challenge faced by the astroparticle physics community, namely adapting the experiment-specific simulation software to work with the CORSIKA air shower simulation output. The proposed solution is an extension of the open source code gSeaGen, which allows the transport of muons generated by CORSIKA to a detector of any size at an arbitrary depth. The gSeaGen code was not only extended in terms of functionality but also underwent a thorough redesign of the muon propagation routine, resulting in a more accurate and efficient simulation. This paper presents the capabilities of the new gSeaGen code as well as prospects for further developments. Program summary: Program title: gSeaGen CPC Library link to program files: https://doi.org/10.17632/ymgxvy2br4.2 Developer's respository link: git.km3net.de/opensource/gseagen Licensing provisions: BSD 3-Clause Programming language: C++ Nature of problem: Integration of the state-of-the-art extensive air shower Monte Carlo event generator CORSIKA [1] into the atmospheric muon simulation for water Cherenkov neutrino telescopes. The primary use case considered is the KM3NeT experiment [2], however, the code should be able to cover other similar experiments as well. The challenges in this work included interfacing the CORSIKA binary output, efficient handling of already generated events to reduce the overall computational cost, and preserving all the additional available information, which can be invaluable in physics analyses. Solution method: The readout of CORSIKA simulation was adapted from the base script provided together with CORSIKA and implemented as a standalone flux driver in gSeaGen. The propagation routine has been redesigned to support the geometry of extensive air shower simulations and to improve its efficiency in propagating particles to the detector. To ensure a reliable modelling of muon energy loss and scattering, PROPOSAL [3] was set as the default internal code for muon transport. PROPOSAL is an open-source software developed and maintained by the IceCube collaboration [4] and is a well-established solution used by the neutrino physics community. Additional comments including restrictions and unusual features: The code was tested with GENIE [5] version 3.4.0 and PROPOSAL 6.1.5. Currently, linking of gSeaGen to GENIE is mandatory, even in the case of a muon-only simulation using CORSIKA. References: [1] D. Heck et al., FZKA-6019 (1998). [2] S. Adrián-Martínez et al., J. Phys. G, Nucl. Part. Phys. 43 (2016) 084001. [3] J.H. Koehne et al., Comput. Phys. Commun. 184 (2013) 2070. [4] M.G. Aartsen et al., J. Instrum. 12 (2017) P03012. [5] C. Andreopoulos et al., Nucl. Instrum. Methods A 614 (2010) 87.

Název v anglickém jazyce

gSeaGen code by KM3NeT: An efficient tool to propagate muons simulated with CORSIKA

Popis výsledku anglicky

The KM3NeT Collaboration has tackled a common challenge faced by the astroparticle physics community, namely adapting the experiment-specific simulation software to work with the CORSIKA air shower simulation output. The proposed solution is an extension of the open source code gSeaGen, which allows the transport of muons generated by CORSIKA to a detector of any size at an arbitrary depth. The gSeaGen code was not only extended in terms of functionality but also underwent a thorough redesign of the muon propagation routine, resulting in a more accurate and efficient simulation. This paper presents the capabilities of the new gSeaGen code as well as prospects for further developments. Program summary: Program title: gSeaGen CPC Library link to program files: https://doi.org/10.17632/ymgxvy2br4.2 Developer's respository link: git.km3net.de/opensource/gseagen Licensing provisions: BSD 3-Clause Programming language: C++ Nature of problem: Integration of the state-of-the-art extensive air shower Monte Carlo event generator CORSIKA [1] into the atmospheric muon simulation for water Cherenkov neutrino telescopes. The primary use case considered is the KM3NeT experiment [2], however, the code should be able to cover other similar experiments as well. The challenges in this work included interfacing the CORSIKA binary output, efficient handling of already generated events to reduce the overall computational cost, and preserving all the additional available information, which can be invaluable in physics analyses. Solution method: The readout of CORSIKA simulation was adapted from the base script provided together with CORSIKA and implemented as a standalone flux driver in gSeaGen. The propagation routine has been redesigned to support the geometry of extensive air shower simulations and to improve its efficiency in propagating particles to the detector. To ensure a reliable modelling of muon energy loss and scattering, PROPOSAL [3] was set as the default internal code for muon transport. PROPOSAL is an open-source software developed and maintained by the IceCube collaboration [4] and is a well-established solution used by the neutrino physics community. Additional comments including restrictions and unusual features: The code was tested with GENIE [5] version 3.4.0 and PROPOSAL 6.1.5. Currently, linking of gSeaGen to GENIE is mandatory, even in the case of a muon-only simulation using CORSIKA. References: [1] D. Heck et al., FZKA-6019 (1998). [2] S. Adrián-Martínez et al., J. Phys. G, Nucl. Part. Phys. 43 (2016) 084001. [3] J.H. Koehne et al., Comput. Phys. Commun. 184 (2013) 2070. [4] M.G. Aartsen et al., J. Instrum. 12 (2017) P03012. [5] C. Andreopoulos et al., Nucl. Instrum. Methods A 614 (2010) 87.

Druh

J_imp - Článek v periodiku v databázi Web of Science

CEP obor

—

OECD FORD obor

10308 - Astronomy (including astrophysics,space science)

Projekt

<a href="/cs/project/GA24-12702S" target="_blank" >GA24-12702S: Experimentální a teoretický výzkum astrofyzikálních vysokoenergetických neutrin v neutrinovém teleskopu KM3NeT</a><br>

Návaznosti

P - Projekt vyzkumu a vyvoje financovany z verejnych zdroju (s odkazem do CEP)

Rok uplatnění

2025

Kód důvěrnosti údajů

S - Úplné a pravdivé údaje o projektu nepodléhají ochraně podle zvláštních právních předpisů

Název periodika

Computer Physics Communications

ISSN

0010-4655

e-ISSN

1879-2944

Svazek periodika

314

Číslo periodika v rámci svazku

109660

Stát vydavatele periodika

NL - Nizozemsko

Počet stran výsledku

15

Strana od-do

1-15

Kód UT WoS článku

001511648500001

EID výsledku v databázi Scopus

2-s2.0-105005518957