I am recently trying to understand dimensional regularization in the context of quantum field theory. So to solve an integral $$ \int_{\mathbb R^d} \frac{\text d ^d p}{(2 \pi)^d} \frac{1}{(p^2 + m^2)^2} $$ which is divergent for $d = 4$ it is often said that if one assumes that $d \in \mathbb C$ with sufficiently small $\text{Re}(d)$ it is possible to do some kind of analytic continuation. But in examples it is often shown that one uses $d = 4 - \varepsilon$ and then transforms to hyper spherical coordinates to get an expression like: $$ \int_{\mathbb R^+}\frac{\text d p}{(2\pi)^{4-\varepsilon}} \frac{2\pi^{(4-\varepsilon)/2}}{\Gamma\left(\frac{4 - \varepsilon}{2}\right)} \frac{p^{3-\varepsilon}}{(p^2 + m^2)^2} $$ which has a solution. We then expand for small values of $\varepsilon$.
Question:
Why is it necessary for $d$ to be complex?$\quad\rightarrow$ is there a possibility to stay real ?
Why is it legal to transform in hyper spherical coordinates when complex dimension is assumed ?
Do i get right that we expand the solution to the second integral in a laurent series for small $\varepsilon$ to extract terms which do not diverge at $\varepsilon=0$ and then when adding amplitudes the divergent terms cancel out ( dependent on the process of course ) ?