Python integration with poles

Question

I am trying to integrate an expression with poles, for example:

where the poles are at z = ib-a, ib + a, -ib-a, -ib + a

Python has a range of integration tools in scipy , but I don't know if there are Python functions for this type of integration, or how to specify poles in an integrand.