Consider a multiplicative function $f(n)$ that we write in the form $f(n) = h(n)n^a$ for a certain $a>0$ and $h(n)$ a multiplicative function such that $h(n) \asymp 1$ (more precisely, we can take $h(p^k) = (1-p^{-2})^3$ for all $k>0$, even if in fact I am interested in this function with two different values of $h(p)$ and $h(p^2)$). I denote the Dirichlet series $D_f(s)$ associated to $f$.
My question is: what is known about the poles and residues of $D_f$?
I can write what I believe to be true, but maybe extra hypotheses are needed on the function $h(n)$. At least $D_f$ converges for $\Re(s)>a+1$ by the bound on $f(n)$. I think it has a pole at $s=a+1$, and meromorphic continuation up to $\Re(s)>a$. I am particularly interested in the residue of $D_f(s)$ at $s=a$. I guess it is $$R := \prod_p \left(1 - \frac 1p \right) \left( 1 + \frac{h(p)}{p} + \frac{h(p^2)}{p^2} + \cdots \right)$$
I would like to know if this is true, even up to adding some hypothesis, and why it is so. I thought of it as follows: we know the asymptotics $$\sum_{n<x} f(n) = \sum_{n<x} h(n)n^a \sim R \frac{x^{a+1}}{a+1}$$ and we know that the Dirichlet series is related by Mellin transform to this partial sum, more precisely $$D_f(s) = s \int_1^\infty x^{-s-1} \left(\sum_{n<x} f(n) \right) dx \sim \frac{Rs}{a+1} \int_1^\infty x^{a - s} dx \sim \frac{Rs}{(a+1)(a + 1 - s)} $$
Here the residue at $s=a+1$ is obviously $R$, but of course I cannot substitute that roughly the equivalent.