def stems(x, y, ax=None, label=None, **kwargs):
    """
    A faster replacement for stem plot using vlines (= LineCollection)
    LineCollection keywords are supported.
    """
    bottom = kwargs.pop('bottom', 0)
    ax.axhline(bottom, **kwargs)
    ax.vlines(x, y, bottom, label=label, **kwargs)
