where \(n_1\)  and \(n_2\)  are integers and \(d= 4-2\epsilon \)  is the space-time dimension.
For such integrals one can derive a closed formula for arbitrary tensor rank \(N\) 
(see e.g. [1]} and Eq. (REF ) in Appendix ).
After performing the \(p_\ell \)  integration, we obtain integrals with an effective propagator \(1/q^2\) 
raised to an \(\epsilon \) -dependent power.
