IAP GITLAB

Skip to content
Snippets Groups Projects
Commit 03e4d693 authored by Felix Riehn's avatar Felix Riehn
Browse files

using particle conversion

parent c8e1c4a6
No related branches found
No related tags found
No related merge requests found
...@@ -18,10 +18,11 @@ ...@@ -18,10 +18,11 @@
#include <corsika/setup/SetupTrajectory.h> #include <corsika/setup/SetupTrajectory.h>
#include <corsika/random/RNGManager.h> #include <corsika/random/RNGManager.h>
#include <corsika/cascade/sibyll2.3c.h> #include <corsika/cascade/sibyll2.3c.h>
#include <corsika/cascade/SibStack.h> #include <corsika/cascade/SibStack.h>
#include <corsika/process/sibyll/ParticleConversion.h>
//#include <corsika/units/PhysicalConstants.h>
#include <corsika/units/PhysicalUnits.h> #include <corsika/units/PhysicalUnits.h>
using namespace corsika; using namespace corsika;
using namespace corsika::process; using namespace corsika::process;
...@@ -41,6 +42,11 @@ public: ...@@ -41,6 +42,11 @@ public:
template <typename Particle> template <typename Particle>
double MinStepLength(Particle& p) const { double MinStepLength(Particle& p) const {
// beam particles for sibyll : 1, 2, 3 for p, pi, k // 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; int kBeam = 1;
/* /*
......
...@@ -38,5 +38,6 @@ void Cascade::Step(Particle& particle) { ...@@ -38,5 +38,6 @@ void Cascade::Step(Particle& particle) {
corsika::geometry::LineTrajectory trajectory = corsika::geometry::LineTrajectory trajectory =
fProcesseList.Transport(particle, nextStep); fProcesseList.Transport(particle, nextStep);
sequence.DoContinuous(particle, trajectory); sequence.DoContinuous(particle, trajectory);
// whats going on here? Everywhere else DoDiscrete is passed a Stack reference as well
sequence.DoDiscrete(particle); sequence.DoDiscrete(particle);
} }
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment