set vertical eas to produce FPE due to zero tracklength

Documentation/Examples/vertical_EAS.cc

@@ -31,8 +31,6 @@
 #include <corsika/process/sibyll/Interaction.h>
 #include <corsika/process/sibyll/NuclearInteraction.h>
 
-//#include <corsika/process/pythia/Decay.h>
-
 #include <corsika/process/urqmd/UrQMD.h>
 
 #include <corsika/process/particle_cut/ParticleCut.h>
@@ -64,7 +62,6 @@ using namespace corsika::units::si;
 void registerRandomStreams() {
   random::RNGManager::GetInstance().RegisterRandomStream("cascade");
   random::RNGManager::GetInstance().RegisterRandomStream("s_rndm");
-  //  random::RNGManager::GetInstance().RegisterRandomStream("pythia");
   random::RNGManager::GetInstance().RegisterRandomStream("UrQMD");
 
   random::RNGManager::GetInstance().SeedAll();
@@ -103,9 +100,9 @@ int main() {
   // setup particle stack, and add primary particle
   setup::Stack stack;
   stack.Clear();
-  const Code beamCode = Code::Proton;
+  const Code beamCode = Code::RhoMinus; //Code::KStarMinus;
   auto const mass = particles::GetMass(beamCode);
-  const HEPEnergyType E0 = 0.1_PeV;
+  const HEPEnergyType E0 = 0.1_TeV;
   double theta = 0.;
   double phi = 0.;
 
@@ -155,7 +152,7 @@ int main() {
       particles::Code::KMinus, particles::Code::K0Long,  particles::Code::K0Short};
 
   process::sibyll::Interaction sibyll;
-  process::sibyll::NuclearInteraction sibyllNuc(sibyll, env);
+  //process::sibyll::NuclearInteraction sibyllNuc(sibyll, env);
   process::sibyll::Decay decay;
 
   //process::pythia::Decay decay(trackedHadrons);
@@ -173,7 +170,7 @@ int main() {
 
   process::UrQMD::UrQMD urqmd;
 
-  auto sibyllSequence = sibyll << sibyllNuc;
+  auto sibyllSequence = sibyll; // << sibyllNuc;
   process::switch_process::SwitchProcess switchProcess(urqmd, sibyllSequence, 55_GeV);
   auto sequence = switchProcess << decay << eLoss << cut << observationLevel
                                 << trackWriter;