Gradient Hölder regularity for weak solutions of nonlinear sub-elliptic systems with drift terms in the Heisenberg group