def METHOD_NAME(system, lmax, average=False, wl=False):
    """Returns Minkowski Structure Metrics up to a maximum l value."""
    voro = freud.locality.Voronoi().compute(system=system)
    op = freud.order.Steinhardt(
        l=list(range(lmax + 1)), average=average, weighted=True, wl=wl, wl_normalize=wl
    )
    op.compute(system, neighbors=voro.nlist)
    return op.particle_order