In summary, we have proposed a consistent method to construct the density matrix for surface hopping simulations. Both population and coherence distributions have been obtained. In particular, the occupation of active states and the electronic wavefunction have been utilized to define the amplitude and phase of coherence, respectively. It has been shown that the separation of WPs on different adiabatic PESs can be properly captured by the adiabatic nuclear dynamics on the active PES, and thus the occupation of active states is more reliable than the electronic wavefunction when analyzing the adiabatic populations. Utilizing this advantage, our consistent construction method has significantly improved the performance even using traditional FSSH trajectories as demonstrations. Most importantly, the consistently constructed density matrix from BCSH trajectories has shown highly accurate results in representative one-and twodimensional scattering models.