I am reading about amplitude techniques, such as the spinor-helicity formalism: more specifically, the BCFW recursion relations. The lecture notes I am reading are titled as Scattering Amplitudes, and are written by Elvang and Huang. At some point in Chapter 3, the authors claim that I can shift some of the momenta of the outgoing particles, $$p_i^{\mu}\rightarrow \hat{p}_i^{\mu}=p_i^{\mu}+zr_i^{\mu}\tag{3.1}$$ such that the scattering amplitude at hand, $A_n$, becomes holomorphic, i.e. $A_n\rightarrow A_n(z)$. By applying some complex analysis techniques, we are allowed to write the physical amplitude $A_n$ (obtained by setting $z=0$) as $$A_n=-\sum_{z_{\mathcal{I}}}\text{Res}_{z=z_{\mathcal{I}}} \frac{A_n(z)}{z}+B_n\tag{3.4}$$ where $B_n$ is the residue of $A_n(z)/z$ at $z\rightarrow\infty$.
Later on, they argue that each of the residues in the same takes the following form $$\text{Res}_{z=z_{\mathcal{I}}}=-\hat{A}_L(z_{\mathcal{I}}) \frac{1}{P_I^2}\hat{A}_R(z_{\mathcal{I}})\tag{3.5}$$ where $\hat{A}_L$ and $\hat{A}_R$ are two on-shell sub-amplitudes of $A_n$, and $P_I$ is a scalar propagator connecting the Feynman diagrams of the two sub-amplitudes.
I have three questions that are related to that:
Why would the amplitude factorize like that?
Why are the sub-amplitudes on-shell? Is it guaranteed beforehand, or is it an empirical statement?
Later on in eq. (3.12), when they apply the BCFW recursion relations to an $n$ point amplitude $\hat{A}_L$ takes arguments $(1,\hat{P}_I^{h_I},k,...,n)$, whereas $\hat{A}_R$ takes arguments $(-\hat{P}_I^{-h_I},2,...,k-1)$. Why have the arguments corresponding to the propagator lines have been identified as the momentum carried by the scalar propagator?
Any help will be appreciated.