From d7aa7265182aa0b6791de3feb1a8deec8e379ea7 Mon Sep 17 00:00:00 2001
From: Remy Prechelt <prechelt@hawaii.edu>
Date: Fri, 21 May 2021 11:04:01 -1000
Subject: [PATCH] Add helper to create US standard atmosphere.
---
.../LayeredSphericalAtmosphereBuilder.inl | 2 +-
corsika/framework/utility/ImplementsMixin.hpp | 40 +++++++++++++++
corsika/media/USStandardAtmosphere.hpp | 51 +++++++++++++++++++
examples/vertical_EAS.cpp | 6 +--
4 files changed, 95 insertions(+), 4 deletions(-)
create mode 100644 corsika/framework/utility/ImplementsMixin.hpp
create mode 100644 corsika/media/USStandardAtmosphere.hpp
diff --git a/corsika/detail/media/LayeredSphericalAtmosphereBuilder.inl b/corsika/detail/media/LayeredSphericalAtmosphereBuilder.inl
index 0f17cf7c8..88a49ee76 100644
--- a/corsika/detail/media/LayeredSphericalAtmosphereBuilder.inl
+++ b/corsika/detail/media/LayeredSphericalAtmosphereBuilder.inl
@@ -167,7 +167,7 @@ namespace corsika {
template <typename TMediumInterface, template <typename> typename MExtraEnvirnoment>
struct make_layered_spherical_atmosphere_builder {
template <typename... TArgs>
- static auto create(Point const& center, LengthType planetRadius, TArgs... args) {
+ static auto create(Point const& center, LengthType const earthRadius, TArgs... args) {
return LayeredSphericalAtmosphereBuilder<TMediumInterface, MExtraEnvirnoment,
TArgs...>{std::forward<TArgs>(args)...,
center, planetRadius};
diff --git a/corsika/framework/utility/ImplementsMixin.hpp b/corsika/framework/utility/ImplementsMixin.hpp
new file mode 100644
index 000000000..45eec946e
--- /dev/null
+++ b/corsika/framework/utility/ImplementsMixin.hpp
@@ -0,0 +1,40 @@
+/*
+ * (c) Copyright 2021 CORSIKA Project, corsika-project@lists.kit.edu
+ *
+ * This software is distributed under the terms of the GNU General Public
+ * Licence version 3 (GPL Version 3). See file LICENSE for a full version of
+ * the license.
+ */
+#pragma once
+
+/**
+ * @file ImplementsMixin.hpp
+ */
+
+#include <type_traits>
+
+namespace corsika {
+
+ namespace detail {
+
+ /**
+ Helper traits class (partial) for static compile time checking.
+
+ This method checks whether a given class implements a particular
+ mixin - this is particularly useful in the media heirarchy that utilizes
+ mixins for the interface definition.
+ */
+ template <template <typename> typename Mixin, typename T>
+ struct implements_mixin {
+
+ template <typename U>
+ static std::true_type test(Mixin<U>&);
+ static std::false_type test(...);
+
+ public:
+ using type = decltype(test(std::declval<T&>()));
+ static constexpr bool value = type::value;
+ };
+
+ } // namespace detail
+} // namespace corsika
diff --git a/corsika/media/USStandardAtmosphere.hpp b/corsika/media/USStandardAtmosphere.hpp
new file mode 100644
index 000000000..7194d79f1
--- /dev/null
+++ b/corsika/media/USStandardAtmosphere.hpp
@@ -0,0 +1,51 @@
+/*
+ * (c) Copyright 2021 CORSIKA Project, corsika-project@lists.kit.edu
+ *
+ * This software is distributed under the terms of the GNU General Public
+ * Licence version 3 (GPL Version 3). See file LICENSE for a full version of
+ * the license.
+ */
+
+#pragma once
+
+#include <corsika/media/IRefractiveIndexModel.hpp>
+#include <corsika/media/LayeredSphericalAtmosphereBuilder.hpp>
+#include <corsika/framework/utility/ImplementsMixin.hpp>
+
+namespace corsika {
+
+ template <typename TEnvironmentInterface, template <typename> typename TExtraEnv,
+ typename TEnvironment, typename... TArgs>
+ auto create_us_standard_atmosphere(TEnvironment& env, Point const& center,
+ LengthType const& earthRadius, TArgs... args) {
+
+ // construct the atmosphere builder
+ auto builder = make_layered_spherical_atmosphere_builder<
+ TEnvironmentInterface, TExtraEnv>::create(center, earthRadius,
+ std::forward<TArgs>(args)...);
+
+ // as per the vertical_EAS reference, we do not include Ar for now.
+ // TODO: This is not a US standard atmosphere
+ builder.setNuclearComposition(
+ {{Code::Nitrogen, Code::Oxygen}, {0.7847f, 1.f - 0.7847f}});
+
+ // add the standard atmosphere layers
+ builder.addExponentialLayer(1222.6562_g / (1_cm * 1_cm), 994186.38_cm, 2_km);
+ builder.addExponentialLayer(1222.6562_g / (1_cm * 1_cm), 994186.38_cm, 4_km);
+ builder.addExponentialLayer(1144.9069_g / (1_cm * 1_cm), 878153.55_cm, 10_km);
+ builder.addExponentialLayer(1305.5948_g / (1_cm * 1_cm), 636143.04_cm, 40_km);
+ builder.addExponentialLayer(540.1778_g / (1_cm * 1_cm), 772170.16_cm, 100_km);
+ builder.addLinearLayer(1e9_cm, 112.8_km + constants::EarthRadius::Mean);
+
+ // check if we want to also add the US standard refractivity
+ if constexpr (detail::implements_mixin<IRefractiveIndexModel,
+ TEnvironmentInterface>::value) {
+
+ // TODO: Add US Standard refractivity
+ }
+
+ // and assemble the environment
+ builder.assemble(env);
+ };
+
+} // namespace corsika
diff --git a/examples/vertical_EAS.cpp b/examples/vertical_EAS.cpp
index 09e229bca..5bc1d5fcd 100644
--- a/examples/vertical_EAS.cpp
+++ b/examples/vertical_EAS.cpp
@@ -34,7 +34,7 @@
#include <corsika/media/FlatExponential.hpp>
#include <corsika/media/HomogeneousMedium.hpp>
#include <corsika/media/IMagneticFieldModel.hpp>
-#include <corsika/media/LayeredSphericalAtmosphereBuilder.hpp>
+#include <corsika/media/USStandardAtmosphere.hpp>
#include <corsika/media/NuclearComposition.hpp>
#include <corsika/media/MediumPropertyModel.hpp>
#include <corsika/media/UniformMagneticField.hpp>
@@ -180,8 +180,8 @@ int main(int argc, char** argv) {
cout << "input momentum: " << plab.getComponents() / 1_GeV
<< ", norm = " << plab.getNorm() << endl;
- 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));
--
GitLab