As scipy is open source, you can actually read the code for integrate.quad, which says:
For finite integration limits, the integration is performed using a Clenshaw-Curtis method which uses Chebyshev moments. ... If one of the integration limits is infinite, then a Fourier integral is computed (assuming w neq 0).
The call in _quad for an infinite bound is to _quadpack._qagie.
