def METHOD_NAME(self, p,
          rmin = None,
          rmax = None):
    weightData = self.computeSliceWeighting(self.positionData, rmin, rmax)
    if p == "inf":
        Ln = la.norm(weightData*self.vectorData, np.inf)
    else:
        Ln = la.norm(weightData*self.vectorData, p)/max(1e-30, sum(weightData))**(1.0/p)
    return Ln