IAP GITLAB

Pythia 8 crash during decay of particles::DsMinus and particles::DsPlus

I've encountered two crashes happening very rarely (so far 2 out of several 100 showers) at high energies (E0 = 4.41e+19 eV and 1.65e+20 eV with vertical_EAS) which I traced back to Pythia. The crashes happened during the decay of particles::DsMinus and particles::DsPlus, respectively.

Here is the gdb backtrace:

#0  0x00000000006524e2 in Pythia8::HelicityMatrixElement::pBreitWigner(double, double, double, double, double) ()
#1  0x0000000000662116 in Pythia8::HMETau2TwoMesonsViaVector::initHadronicCurrent(std::vector<Pythia8::HelicityParticle, std::allocator<Pythia8::HelicityParticle> >&) ()
#2  0x00000000006565b2 in Pythia8::HelicityMatrixElement::decayWeight(std::vector<Pythia8::HelicityParticle, std::allocator<Pythia8::HelicityParticle> >&) ()
#3  0x000000000078470b in Pythia8::TauDecays::decay(int, Pythia8::Event&) ()
#4  0x00000000006d7478 in Pythia8::ParticleDecays::decay(int, Pythia8::Event&) ()
#5  0x000000000061c91d in Pythia8::HadronLevel::next(Pythia8::Event&) ()
#6  0x00000000004e1a07 in Pythia8::Pythia::forceHadronLevel(bool) ()
#7  0x00000000004e2320 in Pythia8::Pythia::next() ()
#8  0x000000000046d939 in corsika::process::pythia::Decay::DoDecay<corsika::stack::StackIteratorInterface<corsika::stack::CombinedStackImpl<corsika::stack::nuclear_extension::NuclearStackExtensionImpl<corsika::stack::super_stupid::SuperStupidStackImpl>, GeometryData<corsika::environment::Environment<corsika::environment::IMediumModel> > >, corsika::setup::detail::StackWithGeometryInterface, corsika::stack::SecondaryView<corsika::stack::CombinedStackImpl<corsika::stack::nuclear_extension::NuclearStackExtensionImpl<corsika::stack::super_stupid::SuperStupidStackImpl>, GeometryData<corsika::environment::Environment<corsika::environment::IMediumModel> > >, corsika::setup::detail::StackWithGeometryInterface> > > (this=0x7ffd93822190, vP=...)
    at /home/ka/ka_ikp/ka_zv2315/corsika/Processes/Pythia/Decay.cc:226

The error message is Program terminated with signal 8, Arithmetic exception..

To upload designs, you'll need to enable LFS and have an admin enable hashed storage. More information