I have $N$ Independent Random variables Laplacian distributions with $\mu=0$ and positive $b=\sigma^2/2$. I also have dominant random variable $(X_s)$ with Laplacian distribution with $\mu=0$ and $b=\sigma^2_s/2.$ $(\sigma^2_s \gg \sigma^2).$ I want to find the joint distribution for $$X= \sum_1^N X_i +X_s$$
My work: I took the product of the CF of the random variables on the right to get
$$ \varphi_X(t)= \frac{2}{2+\sigma^2_s t^2} \left( \frac{2}{2+\sigma^2t^2} \right)^N$$
The denominator of this expression can be simplified to
\begin{align} & (2+\sigma^2_s t^2) (2+\sigma^2t^2)^N \\[6pt] = {} & (2+\sigma^2_s t^2)(\sigma^2_s)^N(2/\sigma^2_s + \sigma^2/\sigma^2_s t^2)^N \\[6pt] \approx {} & (2+\sigma_s^2 t^2)2^N \end{align}
Thus the CF will be that of the $X_s$.
I would like to get some approximate form of the distribution of $X$ based on $N.$ In this approach, I am not able to. Can someone please help?