From 7fc69a281d48d8f9d0baf52ba1f51c0fbf012157 Mon Sep 17 00:00:00 2001
From: Maximilian Reininghaus <maximilian.reininghaus@kit.edu>
Date: Fri, 5 Oct 2018 13:38:31 +0200
Subject: [PATCH] refined HydrogenSphere

---
 Environment/HydrogenSphere/HydrogenSphere.h | 5 ++++-
 1 file changed, 4 insertions(+), 1 deletion(-)

diff --git a/Environment/HydrogenSphere/HydrogenSphere.h b/Environment/HydrogenSphere/HydrogenSphere.h
index aaf15e89..ca5b0f31 100644
--- a/Environment/HydrogenSphere/HydrogenSphere.h
+++ b/Environment/HydrogenSphere/HydrogenSphere.h
@@ -9,6 +9,8 @@
  * a fSphere homogeneously filled with hydrogen
  */
 
+namespace corsika::environment {
+
 class HydrogenSphere {
   CoordinateSystem const& fCS;
   corsika::geometry::Sphere const fSphere;
@@ -25,7 +27,7 @@ public:
                                : corsika::particles::Code::Unknown;
   }
 
-  MassDensityType GetDensity(Point const& p) const { return density; };
+  MassDensityType GetDensity(Point const& p) const { return fSphere.isInside(p) ? density : 0_kg / (meter*meter*meter); };
 
   GetMagneticField(Point const& p) {
     QuantityVector<magnetic_flux_density_d> components{0 * corsika::units::tesla,
@@ -35,4 +37,5 @@ public:
   }
 };
 
+}
 #endif
-- 
GitLab