From 5f34a3590a2a794ec4a1e8cf87bde6321e262ab1 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 231f6030..aa70f3d6 100644
--- a/Documentation/Examples/cascade_example.cc
+++ b/Documentation/Examples/cascade_example.cc
@@ -236,6 +236,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");
@@ -283,7 +286,7 @@ int main() {
process::sibyll::NuclearInteraction sibyllNuc(env, sibyll);
process::sibyll::Decay decay(trackedHadrons);
ProcessCut cut(20_GeV);
- ObservationLevel obsLevel(1400_m);
+ ObservationLevel obsLevel(height_atmosphere - 2000_m);//1400_m);
process::TrackWriter::TrackWriter trackWriter("tracks.dat");
process::EnergyLoss::EnergyLoss eLoss;
@@ -318,7 +321,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