One of the most remarkable characteristics of blazars is their highly variable and luminous emission across the entire electromagnetic spectrum, including very high energy (VHE) gamma rays. Specifically, the TeV gamma-ray variability is a powerful probe of the physical processes in the innermost regions of the jet that constrains the models of particle acceleration and radiation in these extreme environments. Despite the discovery of dozens of blazars with ground-based TeV gamma-ray telescopes in recent decades, the origin of VHE gamma-ray emissions remains unclear.

In this work, we investigate the TeV variability properties of three BL Lac objects: PKS 2155-304, Mrk 421, and Mrk 501, by analyzing the data collected by the MAGIC, H.E.S.S., and VERITAS telescopes over the last two decades. In contrast to previous studies, we employ the same methodology using the Bayesian blocks and the Eisenstein-Hut HOP algorithm to investigate the TeV variability of these BL Lac objects. This consistent framework is intended to enable a comprehensive analysis of the temporal properties and yield deeper insights into the physical mechanisms that may underlie the observed flux variability.