From 3ad101ec4dd7daaa5702f70c784b11da542e3a57 Mon Sep 17 00:00:00 2001
From: Maximilian Reininghaus <maximilian.reininghaus@tu-dortmund.de>
Date: Mon, 13 Jul 2020 00:08:06 +0200
Subject: [PATCH] setup environment w/ magnetic field
---
Documentation/Examples/vertical_EAS.cc | 7 +++++--
Setup/SetupEnvironment.h | 1 +
2 files changed, 6 insertions(+), 2 deletions(-)
diff --git a/Documentation/Examples/vertical_EAS.cc b/Documentation/Examples/vertical_EAS.cc
index c6f2dbe61..7068bb26a 100644
--- a/Documentation/Examples/vertical_EAS.cc
+++ b/Documentation/Examples/vertical_EAS.cc
@@ -15,9 +15,12 @@
#include <corsika/cascade/Cascade.h>
#include <corsika/environment/Environment.h>
#include <corsika/environment/FlatExponential.h>
+#include <corsika/environment/HomogeneousMedium.h>
+#include <corsika/environment/IMagneticFieldModel.h>
#include <corsika/environment/LayeredSphericalAtmosphereBuilder.h>
#include <corsika/environment/NuclearComposition.h>
#include <corsika/environment/ShowerAxis.h>
+#include <corsika/environment/UniformMagneticField.h>
#include <corsika/geometry/Plane.h>
#include <corsika/geometry/Sphere.h>
#include <corsika/logging/Logging.h>
@@ -143,8 +146,8 @@ int main(int argc, char** argv) {
cout << "input momentum: " << plab.GetComponents() / 1_GeV << ", norm = " << plab.norm()
<< endl;
- auto const observationHeight = 0_km + builder.getEarthRadius();
- auto const injectionHeight = 112.75_km + builder.getEarthRadius();
+ auto const observationHeight = 1.4_km + seaLevel;
+ auto const injectionHeight = 112.75_km + seaLevel;
auto const t = -observationHeight * cos(thetaRad) +
sqrt(-units::static_pow<2>(sin(thetaRad) * observationHeight) +
units::static_pow<2>(injectionHeight));
diff --git a/Setup/SetupEnvironment.h b/Setup/SetupEnvironment.h
index f4265e3dd..9c893beb2 100644
--- a/Setup/SetupEnvironment.h
+++ b/Setup/SetupEnvironment.h
@@ -13,6 +13,7 @@
#include <corsika/environment/IMediumModel.h>
#include <corsika/environment/IMediumPropertyModel.h>
#include <corsika/environment/IRefractiveIndexModel.h>
+#include <corsika/environment/IMagneticFieldModel.h>
namespace corsika::setup {
--
GitLab