We extend our earlier "virtual detector" method [X. Wang, J. Tian, and J. H. Eberly, Phys. Rev. Lett. 110, 243001 (2013)], a hybrid quantum mechanical and classical trajectory method, to include phases in the classical trajectories. Effects of quantum interferences, lost in the earlier method, are restored. The obtained photoelectron momentum distributions agree well with the corresponding numerical solutions of the time-dependent Schrödinger equation.