We present 10 years of R-band monitoring data of 31 northern blazars which were either detected at very high energy (VHE) gamma rays or listed as potential VHE gamma-ray emitters. The data comprise 11820 photometric data points in the R-band obtained
in 2002-2012. We analyze the light curves by determining their power spectral density (PSD) slopes assuming a power-law dependence with a single slope $beta$ and a Gaussian probability density function (PDF). We use the multiple fragments variance function (MFVF) combined with a forward-casting approach and likelihood analysis to determine the slopes and perform extensive simulations to estimate the uncertainties of the derived slopes. We also look for periodic variations via Fourier analysis and quantify the false alarm probability through a large number of simulations. Comparing the obtained PSD slopes to values in the literature, we find the slopes in the radio band to be steeper than those in the optical and gamma rays. Our periodicity search yielded one target, Mrk 421, with a significant (p<5%) period. Finding one significant period among 31 targets is consistent with the expected false alarm rate, but the period found in Mrk~421 is very strong and deserves further consideration}.
We have derived R-band host galaxy fluxes of 20 TeV candidate BL Lacertae objects as a function of aperture radius and FWHM. The results are given as correction tables, that list the fluxes (in mJy) of all ``contaminating sources (host galaxy + signi
ficant nearby objects) as a function of aperture radius and FWHM. We found that the derived fluxes depend strongly on aperture radius, but the FWHM has only a minor effect (a few percent). We also discuss the implications of our findings to optical monitoring programs and potential sources of error in our derived fluxes. During this work we have also constructed new calibration star sequences for 9 objects and present the finding charts and calibrated magnitudes.
