diff --git a/DDSegmentation/include/DDSegmentation/Segmentation.h b/DDSegmentation/include/DDSegmentation/Segmentation.h
index b20b462ffae64fbc6d46fa195808ac6786c20cff..68de568d62a00b4ab27736dd9bec39a3534ff3d5 100644
--- a/DDSegmentation/include/DDSegmentation/Segmentation.h
+++ b/DDSegmentation/include/DDSegmentation/Segmentation.h
@@ -133,6 +133,11 @@ protected:
 	/// Helper method to convert a 1D position to a cell ID
 	static int positionToBin(double position, double cellSize, double offset = 0.);
 
+	/// Helper method to convert a bin number to a 1D position given a vector of binBoundaries
+        static double binToPosition(CellID bin, std::vector<double> const& cellBoundaries, double offset = 0.);
+	/// Helper method to convert a 1D position to a cell ID given a vector of binBoundaries
+	static int positionToBin(double position, std::vector<double> const& cellBoundaries, double offset = 0.);
+
 	/// The segmentation name
 	std::string _name;
 	/// The segmentation type
diff --git a/DDSegmentation/src/Segmentation.cpp b/DDSegmentation/src/Segmentation.cpp
index 2995828cd995a199ae2fc34872311c8be52cee4b..8430bb59c00d4ef26b808a01710f8eb15adb3e52 100644
--- a/DDSegmentation/src/Segmentation.cpp
+++ b/DDSegmentation/src/Segmentation.cpp
@@ -11,6 +11,7 @@
 #include <sstream>
 #include <stdexcept>
 #include <cmath>
+#include <algorithm>
 
 namespace DD4hep {
 namespace DDSegmentation {
@@ -144,5 +145,32 @@ int Segmentation::positionToBin(double position, double cellSize, double offset)
 	return int(floor((position + 0.5 * cellSize - offset) / cellSize));
 }
 
+/// Helper method to convert a bin number to a 1D position given a vector of binBoundaries
+double Segmentation::binToPosition(CellID bin, std::vector<double> const& cellBoundaries, double offset) {
+  return (cellBoundaries[bin+1] + cellBoundaries[bin])*0.5 + offset;
+}
+/// Helper method to convert a 1D position to a cell ID given a vector of binBoundaries
+int Segmentation::positionToBin(double position, std::vector<double> const& cellBoundaries, double offset) {
+
+  // include the lower edge to the segmentation
+  if(fabs(position - cellBoundaries.front()) < 1e-12) return 0;
+
+  // include the upper edge of the last bin to the segmentation
+  if(fabs(position - cellBoundaries.back())  < 1e-12) return int(cellBoundaries.size()-2);
+
+  // hits outside cannot be treated
+  if(position < cellBoundaries.front()) throw std::runtime_error("Hit Position is outside of segmentation");
+  if(position > cellBoundaries.back() ) throw std::runtime_error("Hit Position is outside of segmentation");
+
+
+  std::vector<double>::const_iterator bin = std::upper_bound(cellBoundaries.begin(),
+							     cellBoundaries.end(),
+							     position-offset);
+
+  // need to reduce found bin by one, because upper_bound works that way, lower_bound would not help
+  return bin - cellBoundaries.begin() - 1 ;
+
+}
+
 } /* namespace DDSegmentation */
 } /* namespace DD4hep */