Higher order tensor renormalization group for relativistic fermion systems