feat: detray ranges

core/CMakeLists.txt

@@ -113,7 +113,15 @@ detray_add_library( detray_core core
    "include/detray/tracks/tracks.hpp"
    # utils include(s)
    "include/detray/utils/column_wise_operator.hpp"
-   "include/detray/utils/enumerate.hpp"
+   "include/detray/utils/ranges/detail/iterator_functions.hpp"
+   "include/detray/utils/ranges/empty.hpp"
+   "include/detray/utils/ranges/enumerate.hpp"
+   "include/detray/utils/ranges/iota.hpp"
+   "include/detray/utils/ranges/join.hpp"
+   "include/detray/utils/ranges/pick.hpp"
+   "include/detray/utils/ranges/ranges.hpp"
+   "include/detray/utils/ranges/single.hpp"
+   "include/detray/utils/ranges/subrange.hpp"
    "include/detray/utils/invalid_values.hpp"
    "include/detray/utils/quadratic_equation.hpp"
    "include/detray/utils/ratio.hpp"

core/include/detray/core/detail/tuple_container.hpp

@@ -9,6 +9,7 @@
 
 // Project include(s)
 #include "detray/definitions/detail/accessor.hpp"
+#include "detray/definitions/indexing.hpp"
 #include "detray/definitions/qualifiers.hpp"
 
 namespace detray {

core/include/detray/core/detector.hpp

@@ -17,6 +17,7 @@
 #include "detray/geometry/surface.hpp"
 #include "detray/geometry/volume.hpp"
 #include "detray/tools/bin_association.hpp"
+#include "detray/utils/ranges.hpp"
 
 // Vecmem include(s)
 #include <vecmem/memory/memory_resource.hpp>
@@ -479,7 +480,7 @@ class detector {
            << " volumes." << std::endl;
         ss << " local surface finders." << std::endl;
 
-        for (const auto &[i, v] : enumerate(_volumes)) {
+        for (const auto [i, v] : detray::views::enumerate(_volumes)) {
             ss << "[>>] Volume at index " << i << ": " << std::endl;
             ss << " - name: '" << v.name(names) << "'" << std::endl;
 

core/include/detray/core/transform_store.hpp

@@ -10,7 +10,7 @@
 
 #include "detray/definitions/indexing.hpp"
 #include "detray/definitions/qualifiers.hpp"
-#include "detray/utils/enumerate.hpp"
+#include "detray/utils/ranges.hpp"
 
 namespace detray {
 
@@ -71,7 +71,7 @@ class static_transform_store {
     DETRAY_HOST_DEVICE
     inline auto range(const size_t begin, const size_t end,
                       const context & /*ctx*/) const {
-        return iterator_range(_data, dindex_range{begin, end});
+        return detray::ranges::subrange(_data, dindex_range{begin, end});
     }
 
     /** Access to a predefined range of elements
@@ -229,4 +229,4 @@ inline static_transform_store_data<static_transform_store<vector_t> > get_data(
     return static_transform_store_data(store);
 }
 
-}  // namespace detray
+}  // namespace detray

core/include/detray/definitions/detail/accessor.hpp

@@ -17,12 +17,11 @@
 
 #include <algorithm>
 #include <array>
-#include <climits>
+#include <limits>
 #include <tuple>
 #include <type_traits>
 #include <utility>
 
-#include "detray/definitions/indexing.hpp"
 #include "detray/definitions/qualifiers.hpp"
 
 namespace detray {

core/include/detray/definitions/units.hpp

@@ -21,10 +21,21 @@ constexpr double m = 1000.0;
 constexpr double mm3 = mm * mm * mm;
 constexpr double cm3 = cm * cm * cm;
 
+// Time, native unit mm = [speed-of-light * time] = mm/s * s
+constexpr double s = 299792458000.0;
+constexpr double fs = 1e-15 * s;
+constexpr double ps = 1e-12 * s;
+constexpr double ns = 1e-9 * s;
+constexpr double us = 1e-6 * s;
+constexpr double ms = 1e-3 * s;
+constexpr double min = 60.0 * s;
+constexpr double h = 3600.0 * s;
+
 // Energy, native unit GeV
 constexpr double keV = 0.000001;
 constexpr double MeV = 0.001;
 constexpr double GeV = 1.0;
+constexpr double TeV = 1e3;
 
 // Atomic mass unit u
 // 1u == 0.93149410242 GeV/c
@@ -43,6 +54,9 @@ constexpr double mol = 1.0;
 // Avogadro constant
 constexpr double kAvogadro = 6.02214076e23 / unit_constants::mol;
 
+// Charge, native unit e (elementary charge)
+constexpr double e = 1.0;
+
 // Magnetic field, native unit GeV/(e*mm)
 constexpr double T = 0.000299792458;  // equivalent to c in appropriate SI units