  The design and analysis of time-domain sky surveys requires the ability to
simulate accurately realistic populations of core collapse supernova (SN)
events. We present a set of spectral time-series templates designed for this
purpose, for both hydrogen-rich (type II, IIn, IIb) and stripped envelope
(types Ib, Ic, Ic-BL) core collapse supernovae. We use photometric and
spectroscopic data for 67 core collapse supernovae from the literature, and for
each generate a time-series spectral template. The techniques used to build the
templates are fully data-driven with no assumption of any parametric form or
model for the light curves. The template-building code is open-source, and can
be applied to any transient for which well-sampled multi-band photometry and
multiple spectroscopic observations are available. We extend these spectral
templates into the near-ultraviolet to $\lambda \lambda \sim 1600 \AA$ using
observer-frame ultraviolet photometry. We also provide a set of templates
corrected for host galaxy dust extinction, and provide a set of luminosity
functions that can be used with our spectral templates in simulations. We give
an example of how these templates can be used by integrating them within the
popular SN simulation package $\textsc{snana}$, and simulating core collapse
supernovae in photometrically-selected cosmological type Ia supernova samples,
prone to contamination from core collapse events.
