From 773410b8367830a5d4b491a7c775ec380be50a89 Mon Sep 17 00:00:00 2001
From: Felix Riehn <felix@matilda>
Date: Thu, 29 Nov 2018 11:26:02 +0000
Subject: [PATCH] using particle conversion
---
Documentation/Examples/cascade_example.cc | 8 +++++++-
Framework/Cascade/Cascade.cc | 1 +
2 files changed, 8 insertions(+), 1 deletion(-)
diff --git a/Documentation/Examples/cascade_example.cc b/Documentation/Examples/cascade_example.cc
index 8b11f9ddf..13997ef13 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 aeab1de22..584478526 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);
}
--
GitLab