Inverse Mellin transform of 3 gamma functions product I want to calculate the inverse Mellin transform  of products  of 3 gamma functions.
$$F\left ( s \right )=\frac{1}{2i\pi}\int \Gamma(s)\Gamma (2s+a)\Gamma( 2s+b)x^{-s}ds$$
Above contour integral has 3 poles.
$$s_{1}=-n$$
$$s_{2}=-\left ( \frac{n+a}{2} \right )$$
$$s_{3}=-\left ( \frac{n+b}{2} \right )$$
Which pole is used in contour integral?
$$F\left ( s \right )=\sum_{n=0}^{\infty }\frac{\left ( -1 \right )^{n}}{n!}\Gamma \left ( a-2n \right )\Gamma \left ( -2n+b \right )x^{n}+\sum_{n=0}^{\infty}\frac{\left ( -1 \right )^{n}}{n!}\Gamma \left ( -\frac{n}{2}-\frac{a}{2} \right)\Gamma \left ( -n-a+b\right )x^{\frac{n}{2}+\frac{a}{2}}+\sum_{n=0}^{\infty}\frac{(-1)^{n}}{n!}\Gamma \left (-\frac{n}{2}-\frac{b}{2} \right )\Gamma (-n-b+a)x^{\frac{n}{2}+\frac{b}{2}}$$
Is this correct??
 A: To avoid all poles in the Mellin inversion formula you want to integrate along the line $\int_{\gamma-i\infty}^{\gamma+i\infty}ds$ where $\gamma>\max(0,-a/2,-b/2)$; then Mathematica gives the result in terms of the Meijer G-function:
\begin{align}
F( x )=&\frac{1}{2i\pi}\int_{\gamma-i\infty}^{\gamma+i\infty} \Gamma(s)\Gamma (2s+a)\Gamma( 2s+b)x^{-s}\,ds\\
=&\frac{1}{\pi}2^{a+b-2} G_{0,5}^{5,0}\left(\frac{x}{16}|
\begin{array}{c}
 0,\frac{a}{2},\frac{a+1}{2},\frac{b}{2},\frac{b+1}{2} \\
\end{array}
\right).
\end{align}
The integral can alternatively be evaluated by contour integration, closing the contour in the left-half complex plane and picking up the poles at $-n$, $-(n+a)/2$, $-(n+b)/2$, $n=0,1,2,\ldots$. For this we assume that $a\neq b$ and $a,b$ are both non-integer --- so that all poles are simple. I then find that
\begin{align}
F(x)=&\sum_{n=0}^{\infty }\frac{\left ( -1 \right )^{n}}{n!}\Gamma \left ( a-2n \right )\Gamma \left ( -2n+b \right )x^{n}\\
&+\frac{1}{2}\sum_{n=0}^{\infty}\frac{\left ( -1 \right )^{n}}{n!}\Gamma \left ( -\frac{n}{2}-\frac{a}{2} \right)\Gamma \left ( -n-a+b\right )x^{\frac{n}{2}+\frac{a}{2}}\\
&+\frac{1}{2}\sum_{n=0}^{\infty}\frac{(-1)^{n}}{n!}\Gamma \left (-\frac{n}{2}-\frac{b}{2} \right )\Gamma (-n-b+a)x^{\frac{n}{2}+\frac{b}{2}}.
\end{align}
I checked numerically that this agrees with the Meijer G-function result -- it differs from the result in the OP by factors 1/2 in the second and third sum.
