Interactions with Argon not supported by SIBYLL
The corsika7Atmosphere
created by create_5layer_atmosphere
which is currently used in the examples includes Argon as does the real atmosphere. Unfortunately SIBYLL is only made for nuclei up to A=20
. This means currently high energy showers in which enough interactions occur for Argon to be sampled will crash with SIBYLL.
The crash looks like this in case you run into that:
terminate called after throwing an instance of 'std::runtime_error'
what(): Invalid target/projectile/energy combination
The immediate solution is to replace Argon with Oxygen in the SIBYLL interface whenever interactions occur (corsika7 does the same). In the long run SIBYLL may be extended to allow Argon interactions but that will not be soon and it won't solve the problem for SIBYLL 2.1 which we want to add eventually.