One of the most important areas of application for equation-of-motion coupled-cluster (EOM-CC) theory is the prediction, simulation, and analysis of various types of electronic spectra. In this work, the EOM-CC method for ionized states, known as EOM-IP-CC, is applied to the closely lying and coupled pair of states of the ozone cation─X̃2A1 and Ã2B2─using highly accurate treatments including up to the full single, double, triple, and quadruple excitations (EOM-IP-CCSDTQ). Combined with a venerable and powerful method for calculating vibronic spectra from the Hamiltonian produced by EOM-IP-CC calculations, the simulations yield a spectrum that is in good agreement with the photoelectron spectrum of ozone. Importantly, the calculations suggest that the adiabatic gap separating these two electronic states is somewhat smaller than currently thought; an assignment of the simulated spectrum together with more precise band positions of the experimental measurements suggests T00 = 1,368 ± 65 cm-1.