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.
.