corsika/detail/framework/geometry/Point.inl

@@ -67,11 +67,15 @@ namespace corsika {
     return Point(cs, getCoordinates() + pVec.getComponents(cs));
   }
 
+  inline Point Point::operator-(Vector<length_d> const& pVec) const {
+    CoordinateSystemPtr const& cs = BaseVector<length_d>::getCoordinateSystem();
+    return Point(cs, getCoordinates() - pVec.getComponents(cs));
+  }
+
   inline Vector<length_d> Point::operator-(Point const& pB) const {
     CoordinateSystemPtr const& cs = BaseVector<length_d>::getCoordinateSystem();
     return Vector<length_d>(cs, getCoordinates() - pB.getCoordinates(cs));
   }
-
   inline std::ostream& operator<<(std::ostream& os, corsika::Point const& p) {
     auto const& qv = p.getCoordinates();
     os << qv << " (ref:" << fmt::ptr(p.getCoordinateSystem()) << ")";

corsika/framework/geometry/Point.hpp

@@ -72,6 +72,7 @@ namespace corsika {
     inline void rebase(CoordinateSystemPtr const& pCS);
 
     inline Point operator+(Vector<length_d> const& pVec) const;
+    inline Point operator-(Vector<length_d> const& pVec) const;
 
     /*!
      * returns the distance Vector between two points

corsika/modules/proposal/ProposalProcessBase.hpp

@@ -46,9 +46,10 @@ namespace corsika::proposal {
          corsika::units::si::HEPEnergyType
              emCut) { //!< Stochastic losses smaller than the given cut
                       //!< will be handeled continuously.
-        auto p_cut =
-            std::make_shared<const PROPOSAL::EnergyCutSettings>(emCut / 1_MeV, 1, false);
-        return PROPOSAL::GetStdCrossSections(T(), m, p_cut, true);
+        auto particle_def = T();
+        auto p_cut = std::make_shared<const PROPOSAL::EnergyCutSettings>(
+            particle_def.mass + emCut / 1_MeV, 1, false);
+        return PROPOSAL::GetStdCrossSections(particle_def, m, p_cut, true);
       };
 
   //!

data @ 331c0aad

-Subproject commit eded833608bbc6845f9d7e2af4e151683f27818f
+Subproject commit 331c0aad7c6c17fe456c4a80f8b2575ef341ce90

tests/framework/testGeometry.cpp

@@ -60,6 +60,20 @@ TEST_CASE("Geometry CoordinateSystems") {
     CHECK(testV3.getNorm() / (tesla * meter) == Approx(6));
   }
 
+  SECTION("point") {
+    Point const p00 = Point(rootCS, {0_m, 0_m, 0_m});
+    Point const p01 = Point(rootCS, {3_m, 3_m, 3_m});
+    Point const p02 = Point(rootCS, {1_m, -5_m, 6_m});
+
+    LengthVector const d01 = p00 - p01;
+    LengthVector const d12 = p01 - p02;
+    LengthVector const d20 = p02 - p00;
+
+    CHECK(d12.getNorm() / (p01 - p02).getNorm() == Approx(1));
+    CHECK(d12.getNorm() / distance(p01, p02) == Approx(1));
+    CHECK(d12.getNorm() / distance(p02, p01) == Approx(1));
+  }
+
   SECTION("translations") {
     QuantityVector<length_d> const translationVector{0_m, 4_m, 0_m};
     CORSIKA_LOG_INFO("QuantityVector<length_d> translationVector={}", translationVector);
@@ -527,4 +541,4 @@ TEST_CASE("Path") {
     P3.removeFromEnd();
     CHECK(P3.getNSegments() == Approx(0));
   }
-}
\ No newline at end of file
+}