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.