def findSecondary(time, flux, period, epoch, duration):
	t2, f2 = removeTransits(time, flux, period, epoch, duration)
	minp, maxp = period - 0.1, period + 0.1
	if t2[-1] - t2[0] == 0 or 1./maxp < 1./(t2[-1] - t2[0]):
		return (np.nan,)*5
	if minp < 0.5:
		minp = 0.5
	planetInfo = do_bls_and_fit(t2, f2, minp, maxp)
	return (t2,) + planetInfo[0:4] + (planetInfo[5],)