Our study aims to investigate the dynamics of conventional memory T cells (Tconv) and regulatory T cells (Treg) following activation, and to explore potential differences between these two cell types. To achieve this, we developed advanced statistical mixed models and mathematical models of ordinary differential equations (ODE), which allowed us to transform post-vaccination immunological processes into mathematical formulas. These models were applied on in-house data from a de novo Hepatitis B vaccination trial. By accounting for inter- and intra-individual variability, our models provided good fits for both antigen-specific Tconv and Treg cells, overcoming the challenge of studying these complex processes. Our modeling approach provided a deeper understanding of the immunological processes underlying T cell development after vaccination. Specifically, our analysis revealed several important findings regarding the dynamics of Tconv and Treg cells, as well as their relationship to seropositivity for HSV-1 and EBV, and the dynamics of antibody response to vaccination. Firstly, our modeling indicated that Tconv dynamics suggest the existence of two T cell types, in contrast to Treg dynamics where only one T cell type is predicted. Secondly, we found that individuals who converted to a positive antibody response to the vaccine earlier had lower decay rates for both Tregs and Tconv cells, which may have important implications for the development of more effective vaccination strategies. Additionally, our modeling showed that HSV-1 seropositivity negatively influenced Tconv cell expansion after the second vaccination, while EBV seropositivity was associated with higher Treg expansion rates after vaccination. Overall, this study provides a critical foundation for understanding the dynamic processes underlying T cell development after vaccination.