IAP GITLAB

Skip to content
Snippets Groups Projects
Commit 71270188 authored by ralfulrich's avatar ralfulrich
Browse files

just some cleanup

parent cb9eef08
No related branches found
No related tags found
No related merge requests found
...@@ -189,7 +189,7 @@ namespace corsika::process { ...@@ -189,7 +189,7 @@ namespace corsika::process {
}; };
/// the +operator assembles many BaseProcess, ContinuousProcess, and /// the +operator assembles many BaseProcess, ContinuousProcess, and
/// DiscreteProcess objects into a ProcessSequence, all combinatoris /// DiscreteProcess objects into a ProcessSequence, all combinatorics
/// must be allowed, this is why we define a macro to define all /// must be allowed, this is why we define a macro to define all
/// combinations here: /// combinations here:
......
...@@ -12,6 +12,8 @@ ...@@ -12,6 +12,8 @@
#include <corsika/process/stack_inspector/StackInspector.h> #include <corsika/process/stack_inspector/StackInspector.h>
#include <corsika/units/PhysicalUnits.h> #include <corsika/units/PhysicalUnits.h>
#include <corsika/logging/Logger.h>
#include <iostream> #include <iostream>
using namespace std; using namespace std;
...@@ -31,20 +33,15 @@ process::EProcessReturn StackInspector<Stack, Trajectory>::DoContinuous(Particle ...@@ -31,20 +33,15 @@ process::EProcessReturn StackInspector<Stack, Trajectory>::DoContinuous(Particle
Trajectory&, Trajectory&,
Stack& s) const { Stack& s) const {
// using namespace corsika::particles::io;
static int countStep = 0; static int countStep = 0;
if (!fReport) return EProcessReturn::eOk; if (!fReport) return EProcessReturn::eOk;
// std::cout << "generation " << countStep << std::endl;
[[maybe_unused]] int i = 0; [[maybe_unused]] int i = 0;
EnergyType Etot = 0_GeV; EnergyType Etot = 0_GeV;
for (auto& iterP : s) { for (auto& iterP : s) {
EnergyType E = iterP.GetEnergy(); EnergyType E = iterP.GetEnergy();
Etot += E; Etot += E;
// std::cout << " particle data: " << i++ << ", id=" << iterP << " | " << std::endl;
} }
countStep++; countStep++;
// cout << "#=" << countStep << " " << s.GetSize() << " " << Etot/1_GeV << endl;
cout << countStep << " " << s.GetSize() << " " << Etot / 1_GeV << " " << endl; cout << countStep << " " << s.GetSize() << " " << Etot / 1_GeV << " " << endl;
return EProcessReturn::eOk; return EProcessReturn::eOk;
} }
...@@ -54,9 +51,6 @@ double StackInspector<Stack, Trajectory>::MinStepLength(Particle&) const { ...@@ -54,9 +51,6 @@ double StackInspector<Stack, Trajectory>::MinStepLength(Particle&) const {
return 0; return 0;
} }
template <typename Stack, typename Trajectory>
void StackInspector<Stack, Trajectory>::DoDiscrete(Particle&, Stack&) const {}
template <typename Stack, typename Trajectory> template <typename Stack, typename Trajectory>
void StackInspector<Stack, Trajectory>::Init() {} void StackInspector<Stack, Trajectory>::Init() {}
......
...@@ -36,9 +36,6 @@ namespace corsika::process { ...@@ -36,9 +36,6 @@ namespace corsika::process {
// template <typename Particle> // template <typename Particle>
double MinStepLength(Particle&) const; double MinStepLength(Particle&) const;
// template <typename Particle, typename Stack>
void DoDiscrete(Particle&, Stack&) const;
private: private:
bool fReport; bool fReport;
}; };
......
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