diff --git a/postpic/analyzer/particles.py b/postpic/analyzer/particles.py
index 1374a9ef..bac0c2bd 100644
--- a/postpic/analyzer/particles.py
+++ b/postpic/analyzer/particles.py
@@ -577,6 +577,25 @@ def var(self, func, weights=1.0):
         m = np.average(data, weights=w)
         return np.average((data - m)**2, weights=w)
 
+    def quantile(self, func, q, weights=1.0):
+        '''
+        The qth-quantile of the distribution.
+        '''
+        if q < 0 or q > 1:
+            raise ValueError('Quantile q ({:}) must be in range [0, 1]'.format(q))
+        w = self.weight() * weights
+        data = func(self)
+        sortidx = np.argsort(data)
+        wcs = np.cumsum(w[sortidx])
+        idx = np.searchsorted(wcs, wcs[-1]*np.array(q))
+        return data[sortidx[idx]]
+
+    def median(self, func, weights=1.0):
+        '''
+        The median
+        '''
+        return self.quantile(func, 0.5, weights=weights)
+
     # ---- Functions to create a Histogram. ---
 
     def createHistgram1d(self, scalarfx, optargsh={'bins': 300},