We perform an exact diagonalization study of the topological order in topological flat band models through calculating entanglement entropy and spectra of low energy states. We identify multiple independent minimal entangled states, which form a set of orthogonal basis states for the ground-state manifold. We extract the modular transformation matrices S (U) which contains the information of mutual (self) statistics, quantum dimensions and fusion rule of quasi-particles. Moreover, we demonstrate that these matrices are robust and universal in the whole topological phase against different perturbations until the quantum phase transition takes place.