def METHOD_NAME(self, indices):
    subA = []
    subB = []
    for i in indices:
        subA.append(self.a[i])
        subB.append(self.b[i])
    py_dist = UniformNdPy(subA, subB)
    return ot.Distribution(py_dist)