Accurate ab initio modeling of spectroscopic signals in nonlinear electronic spectra, such as bidimensional (2D) spectra, requires the computation of the electronic transitions induced by the incoming pump/probe pulses, resulting in a challenging calculation of many electronic excited states. A protocol is thus required to evaluate the variations of spectral properties, like transition energies and dipole moments, with the computational level, and to estimate the sensitivity of the spectra to these variations. Such a protocol is presented here within the framework of complete and restricted active space self-consistent field (CASSCF/RASSCF) theory and its second-order perturbation theory extensions (CASPT2/RASPT2). The electronic excited-state manifolds of pyrimidine nucleobases (thymine, uracil, and cytosine) are carefully characterized in vacuo employing high-level RAS(0,0|10,8|2,12)//SS-RASPT2 calculations. The results provide a reference data set that can be used for optimizing computational efforts and costs, as required for studying computationally more demanding multichromophoric systems (e.g., di- and oligonucleotides). The spectroscopic signatures of the 2D electronic spectrum of a perfectly stacked uracil–cytosine dimer model are characterized, and experimental setups are proposed that can resolve non-covalent interchromophoric interactions in canonical pyrimidine nucleobase-stacked dimers.