In this work, we investigate the potential connection between entanglement and thermalization, as well as the dynamics of hadronization in QCD physics. Through simulations of quantum dynamics on classical hardware, we analyze the real-time response of the massive Schwinger model coupled to external sources, which serves as an analog for the production and fragmentation of quark jets. This analogy is justified by the shared properties of confinement and chiral symmetry breaking between the Schwinger model and QCD. We observe the growth of time-dependent entanglement between the produced jets, which is driven by the increasing number of contributing eigenstates of the reduced density matrix with sufficiently large and closely spaced eigenvalues. This growth is intrinsically linked to the process of thermalization. The emergent eigenstates are identified as meson-like bound states, and by tracking their temporal evolution, we are able to observe real-time hadronization. At mid-rapidity, the long-time values of local observables approach approximately constant values, indicating the onset of equilibrium and the approach to thermalization.