  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.
