From 1f94c2b1863e4d05d4d93d5fc939e907f2394059 Mon Sep 17 00:00:00 2001
From: ralfulrich <ralf.ulrich@kit.edu>
Date: Thu, 28 Oct 2021 22:37:41 +0200
Subject: [PATCH] examples
---
examples/corsika.cpp | 8 ++++----
examples/em_shower.cpp | 4 ++--
examples/vertical_EAS.cpp | 2 +-
3 files changed, 7 insertions(+), 7 deletions(-)
diff --git a/examples/corsika.cpp b/examples/corsika.cpp
index e42008e80..a3db1ccfa 100644
--- a/examples/corsika.cpp
+++ b/examples/corsika.cpp
@@ -37,7 +37,7 @@
#include <corsika/media/MediumPropertyModel.hpp>
#include <corsika/media/UniformMagneticField.hpp>
#include <corsika/media/ShowerAxis.hpp>
-#include <corsika/media/SlidingPlanarExponential.hpp>
+#include <corsika/media/CORSIKA7Atmospheres.hpp>
#include <corsika/modules/BetheBlochPDG.hpp>
#include <corsika/modules/LongitudinalProfile.hpp>
@@ -212,7 +212,7 @@ int main(int argc, char** argv) {
ofstream atmout("earth.dat");
for (LengthType h = 0_m; h < 110_km; h += 100_m) {
- Point const ptest{rootCS, 0_m, 0_m, builder.getPlanetRadius() + h};
+ Point const ptest{rootCS, 0_m, 0_m, constants::EarthRadius::Mean + h};
auto rho =
env.getUniverse()->getContainingNode(ptest)->getModelProperties().getMassDensity(
ptest);
@@ -257,8 +257,8 @@ int main(int argc, char** argv) {
/* === END: CONSTRUCT PRIMARY PARTICLE === */
/* === START: CONSTRUCT GEOMETRY === */
- auto const observationHeight = 0_km + builder.getPlanetRadius();
- auto const injectionHeight = 111.75_km + builder.getPlanetRadius();
+ auto const observationHeight = 0_km + constants::EarthRadius::Mean;
+ auto const injectionHeight = 111.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/em_shower.cpp b/examples/em_shower.cpp
index 07a2b6391..8dd068d77 100644
--- a/examples/em_shower.cpp
+++ b/examples/em_shower.cpp
@@ -118,8 +118,8 @@ int main(int argc, char** argv) {
cout << "input momentum: " << plab.getComponents() / 1_GeV
<< ", norm = " << plab.getNorm() << endl;
- auto const observationHeight = 1.4_km + builder.getPlanetRadius();
- auto const injectionHeight = 112.75_km + builder.getPlanetRadius();
+ auto const observationHeight = 1.4_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 925824a8e..5fde4edc8 100644
--- a/examples/vertical_EAS.cpp
+++ b/examples/vertical_EAS.cpp
@@ -39,7 +39,7 @@
#include <corsika/media/MediumPropertyModel.hpp>
#include <corsika/media/UniformMagneticField.hpp>
#include <corsika/media/ShowerAxis.hpp>
-#include <corsika/media/SlidingPlanarExponential.hpp>
+#include <corsika/media/CORSIKA7Atmospheres.hpp>
#include <corsika/modules/BetheBlochPDG.hpp>
#include <corsika/modules/LongitudinalProfile.hpp>
--
GitLab