From fa24e6fef9e5165c9828eaf5fc99824bfd7c898c Mon Sep 17 00:00:00 2001
From: ralfulrich <ralf.ulrich@kit.edu>
Date: Fri, 29 Oct 2021 08:46:43 +0200
Subject: [PATCH] fixed examples
---
examples/em_shower.cpp | 1 +
examples/hybrid_MC.cpp | 5 +++--
examples/vertical_EAS.cpp | 1 -
3 files changed, 4 insertions(+), 3 deletions(-)
diff --git a/examples/em_shower.cpp b/examples/em_shower.cpp
index 8dd068d77..62aa91e56 100644
--- a/examples/em_shower.cpp
+++ b/examples/em_shower.cpp
@@ -26,6 +26,7 @@
#include <corsika/media/ShowerAxis.hpp>
#include <corsika/media/MediumPropertyModel.hpp>
#include <corsika/media/UniformMagneticField.hpp>
+#include <corsika/media/CORSIKA7Atmospheres.hpp>
#include <corsika/modules/LongitudinalProfile.hpp>
#include <corsika/modules/ObservationPlane.hpp>
diff --git a/examples/hybrid_MC.cpp b/examples/hybrid_MC.cpp
index 718fdbca2..4a85351b0 100644
--- a/examples/hybrid_MC.cpp
+++ b/examples/hybrid_MC.cpp
@@ -34,6 +34,7 @@
#include <corsika/media/MediumPropertyModel.hpp>
#include <corsika/media/UniformMagneticField.hpp>
#include <corsika/media/ShowerAxis.hpp>
+#include <corsika/media/CORSIKA7Atmospheres.hpp>
#include <corsika/modules/BetheBlochPDG.hpp>
#include <corsika/modules/LongitudinalProfile.hpp>
@@ -170,8 +171,8 @@ int main(int argc, char** argv) {
cout << "input momentum: " << plab.getComponents() / 1_GeV
<< ", norm = " << plab.getNorm() << endl;
- auto const observationHeight = 0_km + builder.getPlanetRadius();
- auto const injectionHeight = 112.75_km + builder.getPlanetRadius();
+ auto const observationHeight = 0_km + constants::EarthRadius::Mean;
+ auto const injectionHeight = 112.75_km + constants::EarthRadius::Mean;
auto const t = -observationHeight * cos(thetaRad) +
sqrt(-static_pow<2>(sin(thetaRad) * observationHeight) +
static_pow<2>(injectionHeight));
diff --git a/examples/vertical_EAS.cpp b/examples/vertical_EAS.cpp
index 5fde4edc8..911f45c52 100644
--- a/examples/vertical_EAS.cpp
+++ b/examples/vertical_EAS.cpp
@@ -34,7 +34,6 @@
#include <corsika/media/FlatExponential.hpp>
#include <corsika/media/HomogeneousMedium.hpp>
#include <corsika/media/IMagneticFieldModel.hpp>
-#include <corsika/media/USStandardAtmosphere.hpp>
#include <corsika/media/NuclearComposition.hpp>
#include <corsika/media/MediumPropertyModel.hpp>
#include <corsika/media/UniformMagneticField.hpp>
--
GitLab