The proper treatment of hadronic resonances plays an important role in many aspects of heavy ion collisions. This is expected to be the case also for hadronization, due to the large degeneracies of excited states, and the abundant production of hadrons from their decays. We first show how a comprehensive treatment of excited meson states can be incorporated into quark recombination, and in extension, into Hybrid Hadronization. We then discuss the quantum mechanics of forming excited states, utilizing the Wigner distribution functions of angular momentum eigenstates of isotropic 3-D harmonic oscillators. We further describe how resonance decays can be handled, based on a set of minimal assumptions, by creating an extension of hadron decays in PYTHIA 8. Finally, we present first results by simulating $e^++e^-$
collisions using PYTHIA and Hybrid Hadronization with excited mesons up to orbital angular momentum $L=4$ and radial quantum number 2. We find that states up to $L=2$ are produced profusely by quark recombination.