From 1d95df3e7a2e5902ca9bcf2403087ff9eb554322 Mon Sep 17 00:00:00 2001
From: Felix Riehn <felix@matilda>
Date: Thu, 6 Feb 2020 18:41:19 +0100
Subject: [PATCH] set vertical eas to produce FPE due to zero tracklength
---
Documentation/Examples/vertical_EAS.cc | 11 ++++-------
1 file changed, 4 insertions(+), 7 deletions(-)
diff --git a/Documentation/Examples/vertical_EAS.cc b/Documentation/Examples/vertical_EAS.cc
index 76d78283..8c8ba1a8 100644
--- a/Documentation/Examples/vertical_EAS.cc
+++ b/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;
--
GitLab