diff --git a/Documentation/Examples/vertical_EAS.cc b/Documentation/Examples/vertical_EAS.cc
index c50c74ae10e378f9647f6feb126aaf6f04f12d0c..f49bc4be92c3175e16da3c42be2837d3a995aca6 100644
--- a/Documentation/Examples/vertical_EAS.cc
+++ b/Documentation/Examples/vertical_EAS.cc
@@ -13,6 +13,7 @@
 #include <corsika/environment/FlatExponential.h>
 #include <corsika/environment/LayeredSphericalAtmosphereBuilder.h>
 #include <corsika/environment/NuclearComposition.h>
+#include <corsika/environment/ShowerAxis.h>
 #include <corsika/geometry/Plane.h>
 #include <corsika/geometry/Sphere.h>
 #include <corsika/process/ProcessSequence.h>
@@ -140,9 +141,8 @@ int main(int argc, char** argv) {
             particles::Code::Proton, E0, plab, injectionPos, 0_ns});
   }
 
-  Line const line(injectionPos, plab.normalized() * 1_m * 1_Hz);
-  auto const velocity = line.GetV0().norm();
-  setup::Trajectory const showerAxis(line, (injectionPos - showerCore).norm() / velocity);
+  environment::ShowerAxis const showerAxis{injectionPos, (injectionPos - showerCore) * 2,
+                                           env};
 
   // setup processes, decays and interactions