I have given a function
$$G=p_1q_1 - p_2q_2$$
on a 4-dimensional phase-space. This function $G$ commutes with the Hamiltonian
$$H= \frac{p_1p_2}{m} + m\omega^2q_1q_2.$$
It generates a flow
$$(\vec{q},\vec{p}) \rightarrow (\vec{Q}(\vec{q},\vec{p},\alpha), \vec{P}(\vec{q},\vec{p},\alpha)) $$
defined by
$$\frac{dQ_i}{d\alpha} = \frac{\partial G}{\partial P_i}; \quad \vec{Q}(\vec{q},\vec{p},0) = \vec{q}$$
$$\frac{dP_i}{d\alpha} = -\frac{\partial G}{\partial Q_i}; \quad \vec{P}(\vec{q},\vec{p},0) = \vec{p}.$$
The calculation is no problem, but I do not understand its background. How does $G$ generate the transformation and why is it first defined in old coordinates $\vec{q}, \vec{p}$ but then used with the new ones $\vec{Q}, \vec{P}$? What is the connection to (Lie-) group theory? And can the transformation be understood in terms of vector fields on a manifold?
Since the whole concept is not really clear to me I would be very grateful for a comprehensive explanation. Book-suggestion are also highly welcome.