  Over seventy blazars identified as sources of TeV gamma ray emission have
been detected, including approximately sixty BL Lac objects and seven flat
spectrum radio quasars. The distribution in space of these objects can be
described by a luminosity function, which gives their comoving space density as
a function of their luminosity. We investigate a source selection method to be
used for determining the luminosity function of TeV gamma-ray blazars using
observations from VERITAS, a ground-based gamma-ray observatory consisting of
an array of four atmospheric Cherenkov telescopes located in southern Arizona.
