+ outcommented code that is not fully equivalent to lower_bound & upper_bound

Resample/Processed/ParticleSpans.cpp

@@ -38,22 +38,30 @@ std::vector<double> bottomLayerCoordinates(const MultiLayer& sample)
     return result;
 }
 
-size_t iLayerTop(double top, const std::vector<double> ZBottoms)
+size_t iLayerTop(double top, const std::vector<double>& ZBottoms)
 {
     if (top <= ZBottoms.back())
         return ZBottoms.size();
-    for (size_t i = 0; i <= ZBottoms.size(); ++i)
+    /* ???
+    for (size_t i = 0; i < ZBottoms.size(); ++i)
         if (ZBottoms.at(ZBottoms.size()-1-i) < top)
-            return ZBottoms.size()-1-i;
-    ASSERT(0);
-//    auto index_above = std::lower_bound(ZBottoms.rbegin(), ZBottoms.rend(), top);
-//    return static_cast<size_t>(ZBottoms.rend() - index_above);
+            return ZBottoms.size()-i;
+    return 0;
+    /* */
+    auto index_above = std::lower_bound(ZBottoms.rbegin(), ZBottoms.rend(), top);
+    return static_cast<size_t>(ZBottoms.rend() - index_above);
 }
 
-size_t iLayerBottom(double bottom, const std::vector<double> ZBottoms)
+size_t iLayerBottom(double bottom, const std::vector<double>& ZBottoms)
 {
     if (bottom < ZBottoms.back())
         return ZBottoms.size();
+    /* ???
+    for (size_t i = 0; i < ZBottoms.size(); ++i)
+        if (ZBottoms.at(ZBottoms.size()-1-i) > bottom)
+            return ZBottoms.size()-i;
+    return 0;
+    /* */
     auto index_below = std::upper_bound(ZBottoms.rbegin(), ZBottoms.rend(), bottom);
     return static_cast<size_t>(ZBottoms.rend() - index_below);
 }