From 8fb7bfe529e48785ee5205053d992b8b527530bc Mon Sep 17 00:00:00 2001
From: ralfulrich <ralf.ulrich@kit.edu>
Date: Tue, 2 Apr 2019 11:43:20 +0200
Subject: [PATCH] new setting, so that ObservationLevel process sees a few
particles
---
Documentation/Examples/cascade_example.cc | 7 +++++--
1 file changed, 5 insertions(+), 2 deletions(-)
diff --git a/Documentation/Examples/cascade_example.cc b/Documentation/Examples/cascade_example.cc
index 6c4036deb..65d835c60 100644
--- a/Documentation/Examples/cascade_example.cc
+++ b/Documentation/Examples/cascade_example.cc
@@ -240,6 +240,9 @@ public:
// The example main program for a particle cascade
//
int main() {
+
+ const LengthType height_atmosphere = 112.8_km;
+
feenableexcept(FE_INVALID);
// initialize random number sequence(s)
random::RNGManager::GetInstance().RegisterRandomStream("cascade");
@@ -281,7 +284,7 @@ int main() {
// random::RNGManager::GetInstance().RegisterRandomStream("pythia");
// process::pythia::Decay decay(trackedHadrons);
ProcessCut cut(20_GeV);
- ObservationLevel obsLevel(1400_m);
+ ObservationLevel obsLevel(height_atmosphere - 2000_m);//1400_m);
// random::RNGManager::GetInstance().RegisterRandomStream("HadronicElasticModel");
// process::HadronicElasticModel::HadronicElasticInteraction
@@ -327,7 +330,7 @@ int main() {
cout << "input angles: theta=" << theta << " phi=" << phi << endl;
cout << "input momentum: " << plab.GetComponents() / 1_GeV << endl;
Point pos(rootCS, 0_m, 0_m,
- 112.8_km); // this is the CORSIKA 7 start of atmosphere/universe
+ height_atmosphere); // this is the CORSIKA 7 start of atmosphere/universe
stack.AddParticle(std::tuple<particles::Code, units::si::HEPEnergyType,
corsika::stack::MomentumVector, geometry::Point,
units::si::TimeType, unsigned short, unsigned short>{
--
GitLab