The tempered Lefschetz thimble method [14] is a parallel-tempering algorithm to solve the sign problem in Monte Carlo simulations. It uses the flow time of the antiholomorphic gradient flow as a tempering parameter and is expected to tame both the sign and multimodal problems simultaneously. In this letter, we further develop the algorithm so that the expectation values can be estimated precisely with a criterion ensuring global equilibrium and the sufficiency of the sample size. The key is the use of the fact that the expectation values should be the same for all replicas due to Cauchy's theorem. To demonstrate that this algorithm works well, we apply it to the quantum Monte Carlo simulation of the Hubbard model away from half-filling on a two-dimensional lattice of small size, and show that the numerical results agree nicely with exact values.