  We combine recent advances in excited state variational principles, fast
multi-Slater Jastrow methods, and selective configuration interaction to create
multi-Slater Jastrow wave function approximations that are optimized for
individual excited states. In addition to the Jastrow variables and linear
expansion coefficients, this optimization includes state-specific orbital
relaxations in order to avoid the compromises necessary in state-averaged
approaches. We demonstrate that, when combined with variance matching to help
balance the quality of the approximation across different states, this approach
delivers accurate excitation energies even when using very modest multi-Slater
expansions. Intriguingly, this accuracy is maintained even when studying a
difficult chlorine-anion-to-$\pi^{*}$ charge transfer in which traditional
state-averaged multi-reference methods must contend with different states that
require drastically different orbital relaxations.
