diff --git a/Documentation/Examples/cascade_example.cc b/Documentation/Examples/cascade_example.cc index 8b11f9ddfb050bc8c9f91c6f9c0a4ee2a6f522b1..13997ef13121518205128e3c8a25fd6294e0f325 100644 --- a/Documentation/Examples/cascade_example.cc +++ b/Documentation/Examples/cascade_example.cc @@ -18,10 +18,11 @@ #include <corsika/setup/SetupTrajectory.h> #include <corsika/random/RNGManager.h> + #include <corsika/cascade/sibyll2.3c.h> #include <corsika/cascade/SibStack.h> +#include <corsika/process/sibyll/ParticleConversion.h> -//#include <corsika/units/PhysicalConstants.h> #include <corsika/units/PhysicalUnits.h> using namespace corsika; using namespace corsika::process; @@ -41,6 +42,11 @@ public: template <typename Particle> double MinStepLength(Particle& p) const { // beam particles for sibyll : 1, 2, 3 for p, pi, k + Code c_id = p.GetPID(); + //sibyll::SibyllCodeIntType s_id = process::sibyll::ConvertToSibyll( p.GetPID() ); + std::cout << "MinStepLength: particle input " << "corsika id: " << c_id << std::endl; + //std::cout << "MinStepLength: particle input " << "sibyll id: " << s_id << std::endl; + // std::cout << "MinStepLength: particle input " << "sibyll id: " << process::sibyll::ConvertToSibyll( p.GetPID() ) << std::endl; int kBeam = 1; /* diff --git a/Framework/Cascade/Cascade.cc b/Framework/Cascade/Cascade.cc index aeab1de22fa08617bbdda403958dca8152fdfdeb..5844785262187001c913927b3fd7bbc7e87b5d44 100644 --- a/Framework/Cascade/Cascade.cc +++ b/Framework/Cascade/Cascade.cc @@ -38,5 +38,6 @@ void Cascade::Step(Particle& particle) { corsika::geometry::LineTrajectory trajectory = fProcesseList.Transport(particle, nextStep); sequence.DoContinuous(particle, trajectory); + // whats going on here? Everywhere else DoDiscrete is passed a Stack reference as well sequence.DoDiscrete(particle); }