diff --git a/.github/workflows/test.yml b/.github/workflows/test.yml
index d028aee86..109a5e7cf 100644
--- a/.github/workflows/test.yml
+++ b/.github/workflows/test.yml
@@ -90,6 +90,7 @@ jobs:
     - uses: actions/checkout@v3
     - name: Install dependencies
       run: |
+        brew update
         brew install \
             cmake \
             coreutils \
diff --git a/CMakeLists.txt b/CMakeLists.txt
index 49dd87f7f..947a8a56f 100644
--- a/CMakeLists.txt
+++ b/CMakeLists.txt
@@ -19,6 +19,15 @@ set( Boost_USE_STATIC_RUNTIME OFF )
 
 set( CMAKE_CXX_STANDARD 14 CACHE STRING "The C++ standard whose features are requested to build this target" )
 mark_as_advanced( FORCE CMAKE_CXX_STANDARD )
+
+# ROOT may force a newer C++ standard
+# Using the same C++ standard as ROOT is recommended in https://root.cern/manual/integrate_root_into_my_cmake_project/ .
+# Older ROOT releases don't export ROOT_CXX_STANDARD at all, so only override our default when it's actually set.
+find_package(ROOT QUIET)
+if(ROOT_FOUND AND KASPER_USE_ROOT AND ROOT_CXX_STANDARD)
+    set( CMAKE_CXX_STANDARD ${ROOT_CXX_STANDARD} CACHE STRING "The C++ standard, as requested by the ROOT library. Can't be changed." FORCE)
+endif()
+
 set( CMAKE_CXX_STANDARD_REQUIRED ON )
 set( CMAKE_CXX_EXTENSIONS ON )
 set( CMAKE_EXPORT_COMPILE_COMMANDS ON )
@@ -187,7 +196,7 @@ mark_as_advanced( CLEAR CMAKE_VERBOSE_MAKEFILE )
 #option( KASPER_USE_BOOST "Build Boost dependent modules" ${Boost_FOUND} )
 set(KASPER_USE_BOOST ON)
 
-find_package(ROOT QUIET)
+#ROOT was already searched at the beginning for the C++ standard
 if(ROOT_FOUND AND NOT DEFINED KASPER_USE_ROOT)
     message("ROOT was found and enabled automatically. Set KASPER_USE_ROOT=OFF to change this.")
 endif()
diff --git a/KGeoBag/Source/Extensions/Mesh/Complex/Source/KGExtrudedSurfaceMesher.cc b/KGeoBag/Source/Extensions/Mesh/Complex/Source/KGExtrudedSurfaceMesher.cc
index 10a92df37..50495553b 100644
--- a/KGeoBag/Source/Extensions/Mesh/Complex/Source/KGExtrudedSurfaceMesher.cc
+++ b/KGeoBag/Source/Extensions/Mesh/Complex/Source/KGExtrudedSurfaceMesher.cc
@@ -261,7 +261,7 @@ void KGExtrudedSurfaceMesher::DiscretizeSegment(const KGExtrudedObject::Arc* arc
     double n1[3] = {};
     n1[2] = 0;
 
-    double p0[3];
+    double p0[3] = {};
     p0[2] = fExtrudedObject->GetZMin();
 
     double xy_len = 0;
diff --git a/KGeoBag/Source/Extensions/Mesh/Complex/Source/KGPortHousingSurfaceMesher.cc b/KGeoBag/Source/Extensions/Mesh/Complex/Source/KGPortHousingSurfaceMesher.cc
index 4a7bf0e69..73e0c7a5d 100644
--- a/KGeoBag/Source/Extensions/Mesh/Complex/Source/KGPortHousingSurfaceMesher.cc
+++ b/KGeoBag/Source/Extensions/Mesh/Complex/Source/KGPortHousingSurfaceMesher.cc
@@ -40,8 +40,8 @@ void KGPortHousingSurfaceMesher::VisitWrappedSurface(KGWrappedSurface<KGPortHous
     double x_loc[3];
     double y_loc[3];
     double z_loc[3];
-    double x_min;
-    double x_max;
+    double x_min = 0.;
+    double x_max = 0.;
 
     KGCoordinateTransform* coordTransform = nullptr;
 
diff --git a/Kassiopeia/Trajectories/Source/KSTrajTermPropagation.cxx b/Kassiopeia/Trajectories/Source/KSTrajTermPropagation.cxx
index 914a4abe9..4e30327c1 100644
--- a/Kassiopeia/Trajectories/Source/KSTrajTermPropagation.cxx
+++ b/Kassiopeia/Trajectories/Source/KSTrajTermPropagation.cxx
@@ -25,7 +25,7 @@ KSTrajTermPropagation::~KSTrajTermPropagation() = default;
 void KSTrajTermPropagation::Differentiate(double /*aTime*/, const KSTrajExactParticle& aParticle,
                                           KSTrajExactDerivative& aDerivative) const
 {
-    KThreeVector tVelocity = fDirection * aParticle.GetVelocity();
+    KThreeVector tVelocity = (+fDirection) * aParticle.GetVelocity();
     KThreeVector tForce =
         aParticle.GetCharge() * (aParticle.GetElectricField() + tVelocity.Cross(aParticle.GetMagneticField()));
 
@@ -39,7 +39,7 @@ void KSTrajTermPropagation::Differentiate(double /*aTime*/, const KSTrajExactSpi
 {
     double TempSN = std::sqrt(aParticle.GetSpin().MagnitudeSquared() - aParticle.GetSpin0() * aParticle.GetSpin0());
 
-    KThreeVector tVelocity = fDirection * aParticle.GetVelocity();
+    KThreeVector tVelocity = (+fDirection) * aParticle.GetVelocity();
 
     KThreeVector MagneticMoment = aParticle.GetGyromagneticRatio() * aParticle.GetSpin() /
                                   aParticle.GetSpin().Magnitude() * aParticle.GetSpinMagnitude() *
@@ -157,7 +157,7 @@ void KSTrajTermPropagation::Differentiate(double aTime, const KSTrajExactTrapped
 void KSTrajTermPropagation::Differentiate(double /*aTime*/, const KSTrajAdiabaticParticle& aParticle,
                                           KSTrajAdiabaticDerivative& aDerivative) const
 {
-    double tLongVelocity = fDirection * aParticle.GetLongVelocity();
+    double tLongVelocity = (+fDirection) * aParticle.GetLongVelocity();
     double tLongitudinalMomentum = aParticle.GetLongMomentum();
     double tTransverseMomentum = aParticle.GetTransMomentum();
     double tLorentzFactor = aParticle.GetLorentzFactor();
@@ -219,7 +219,7 @@ void KSTrajTermPropagation::Differentiate(double /*aTime*/, const KSTrajAdiabati
         GradE1 / E1.Magnitude() - KThreeMatrix::OuterProduct(E1, GradE1Magnitude) / E1.Magnitude() / E1.Magnitude();
     KThreeVector A = Grade1 * e2;
 
-    KThreeVector tVelocity = fDirection * aParticle.GetVelocity();
+    KThreeVector tVelocity = (+fDirection) * aParticle.GetVelocity();
     KThreeVector tForce =
         aParticle.GetCharge() * (aParticle.GetElectricField() + tVelocity.Cross(aParticle.GetMagneticField())) +
         aParticle.GetGyromagneticRatio() * aParticle.GetSpinMagnitude() * katrin::KConst::Hbar() *
@@ -244,7 +244,7 @@ void KSTrajTermPropagation::Differentiate(double /*aTime*/, const KSTrajAdiabati
 void KSTrajTermPropagation::Differentiate(double /*aTime*/, const KSTrajMagneticParticle& aParticle,
                                           KSTrajMagneticDerivative& aDerivative) const
 {
-    KThreeVector tVelocity = fDirection * aParticle.GetMagneticField().Unit();
+    KThreeVector tVelocity = (+fDirection) * aParticle.GetMagneticField().Unit();
 
     aDerivative.AddToVelocity(tVelocity);
 
@@ -254,7 +254,7 @@ void KSTrajTermPropagation::Differentiate(double /*aTime*/, const KSTrajMagnetic
 void KSTrajTermPropagation::Differentiate(double /*aTime*/, const KSTrajElectricParticle& aParticle,
                                           KSTrajElectricDerivative& aDerivative) const
 {
-    KThreeVector tVelocity = fDirection * aParticle.GetElectricField().Unit();
+    KThreeVector tVelocity = (+fDirection) * aParticle.GetElectricField().Unit();
 
     aDerivative.AddToVelocity(tVelocity);
 
diff --git a/Kommon/Base/CMakeLists.txt b/Kommon/Base/CMakeLists.txt
index 54b0cbd4e..5f4b839af 100644
--- a/Kommon/Base/CMakeLists.txt
+++ b/Kommon/Base/CMakeLists.txt
@@ -38,7 +38,6 @@ set(KOMMON_BASE_HEADER_FILES
     Utility/Gnuplot.hpp
 
     Utility/KUtilityMessage.h
-    Utility/KAlgorithm.h
     Utility/KException.h
     Utility/KCast.h
     Utility/KConst_2006.h
diff --git a/Kommon/Base/Utility/KAlgorithm.h b/Kommon/Base/Utility/KAlgorithm.h
deleted file mode 100644
index 0df161508..000000000
--- a/Kommon/Base/Utility/KAlgorithm.h
+++ /dev/null
@@ -1,53 +0,0 @@
-/**
- * @file KAlgorithm.h
- *
- * @date 05.10.2017
- * @author Valerian Sibille <vsibille@mit.edu>
- */
-#ifndef K_ALGORITHM_H_
-#define K_ALGORITHM_H_
-
-#include <algorithm>
-#include <iterator>
-
-namespace katrin
-{
-
-namespace Algorithm
-{
-// the following supports maps implemented as std::vector<std::pair> which lack map::at
-template<class Map, class Key> const auto& FindValue(const Map& map, const Key& key)
-{
-
-    auto it = std::find_if(std::cbegin(map), std::cend(map), [&](const auto& pair) { return pair.first == key; });
-    if (it == std::cend(map))
-        throw std::out_of_range("Utility::Algorithm::FindValue(map, key): key not in map");
-
-    return it->second;
-}
-
-template<class Map, class Value> const auto& FindKey(const Map& map, const Value& value)
-{
-
-    auto it = std::find_if(std::cbegin(map), std::cend(map), [&](const auto& pair) { return pair.second == value; });
-    if (it == std::cend(map))
-        throw std::out_of_range("Utility::Algorithm::FindKey(map, value): value not in map");
-
-    return it->first;
-}
-
-template<class RandomAccessIt, class Distance>
-RandomAccessIt SumStride(RandomAccessIt begin, RandomAccessIt end, Distance stride)
-{
-    while (begin != end) {
-        *begin += *(begin + stride);
-        ++begin;
-    }
-    return begin;
-}
-
-}  // namespace Algorithm
-
-}  // namespace katrin
-
-#endif
diff --git a/Kommon/Math/include/KMathIntegrator.h b/Kommon/Math/include/KMathIntegrator.h
index 1cb6422e6..7725ed7b1 100644
--- a/Kommon/Math/include/KMathIntegrator.h
+++ b/Kommon/Math/include/KMathIntegrator.h
@@ -516,7 +516,7 @@ template<class XIntegrandType>
 inline XFloatT KMathIntegrator<XFloatT, XSamplingPolicy>::QAGS(XIntegrandType&& integrand)
 {
     /*Use gsl implementation for integral int_a^b. Implementation is experimental. */
-	const uint ws_size = 100;
+	const size_t ws_size = 100;
     gsl_integration_workspace* workspace = gsl_integration_workspace_alloc(ws_size);
     FType proxy = integrand;
     gsl_function F=GSLFunction(proxy);
@@ -561,7 +561,7 @@ template<class XIntegrandType>
 inline XFloatT KMathIntegrator<XFloatT, XSamplingPolicy>::QAGIU(XIntegrandType&& integrand)
 {
     /*Use gsl implementation for integral int_xmin^infty. Use case is the krypton spectrum. Implementation is experimental. */
-	const uint ws_size = 100;
+	const size_t ws_size = 100;
     gsl_integration_workspace* workspace = gsl_integration_workspace_alloc(ws_size);
     FType proxy = integrand;
     gsl_function F=GSLFunction(proxy);