Recently we derived a method for constructing a quasi diabatic Hamiltonian, H(Q), in the vicinity of a two or three state conical intersection based on ab initio electronic structure data. Here Q denotes the N^int internal degrees of freedom. By using the derivative couplings, energy difference gradients and average energy gradients the method requires electronic structure data at only ~ 2N^int points. This enables very flexible multireference configuration interaction expansions to be used in determining the electronic structure data ( ~130 million CSFs in one example reported ). The Hamiltonian contains all terms through second order. This Hamiltonian can be transformed into the form used in the multimode vibronic coupling model to describe vibronic spectra. Thus when the minimum energy conical intersection is near the minima of the states in question this transformed Hamiltonian can be used to determine photoelectron spectra. This approach is described and results presented.