BCSH has been benchmarked in a variety of model systems and can recover almost the exact branching ratios as well as the time-dependent populations compared with the exact quantum dynamics. In the present study, we further investigate whether BCSH can produce accurate time-dependent density matrix through properly analyzing the surface hopping trajectories.