started refactoring verticalEAS

Documentation/Examples/vertical_EAS.cc

@@ -222,26 +222,25 @@ int main() {
   // setup environment, geometry
   using EnvType = Environment<setup::IEnvironmentModel>;
   EnvType env;
+  const CoordinateSystem& rootCS = env.GetCoordinateSystem();
   auto& universe = *(env.GetUniverse());
 
   auto theMedium =
-      EnvType::CreateNode<Sphere>(Point{env.GetCoordinateSystem(), 0_m, 0_m, 0_m},
+      EnvType::CreateNode<Sphere>(Point{rootCS, 0_m, 0_m, 0_m},
                                   1_km * std::numeric_limits<double>::infinity());
 
   // fraction of oxygen
-  const float fox = 0.20946;
-  using MyHomogeneousModel = HomogeneousMedium<environment::IMediumModel>;
-  theMedium->SetModelProperties<MyHomogeneousModel>(
-      1_kg / (1_m * 1_m * 1_m),
+  using FlatExp = environment::FlatExponential<environment::IMediumModel>;
+  theMedium->SetModelProperties<FlatExp>(
+      Point{rootCS, 0_m, 0_m, },
       environment::NuclearComposition(
           std::vector<particles::Code>{particles::Code::Nitrogen,
-                                       particles::Code::Oxygen},
-          std::vector<float>{(float)1. - fox, fox}));
+                                       particles::Code::Oxygen,
+                                       particles::Code::Argon},
+          std::vector<float>{0.7847,0.2105,0.005})); // values taken from AIRES manual
 
   universe.AddChild(std::move(theMedium));
 
-  const CoordinateSystem& rootCS = env.GetCoordinateSystem();
-
   // setup processes, decays and interactions
   tracking_line::TrackingLine tracking;