Considering two continuous random variables $X$ and $Y$ with $d.f \; F_X, F_Y$ I want to fin the distribution and distribution function of the sum $Z=X+Y$.
\begin{align} P\{Z \leq z\} &= P\{X+Y \leq z\} \\ &=P\{X \leq z - Y\} \\ &=F_X(z - Y) \end{align}
I know here that I have to use somehow the joint distribution but I can't figure out how. We know that if $T(\omega) = (X(\omega), Y(\omega))$:
\begin{align} P\{T \leq (a,b)\} &= P\{X\leq a, Y \leq b\} = \int \int_{x \leq a, y \leq b} f_{XY}(a,b)\;dxdy \end{align}
Using the Radon-Nikodym Theorem if and only $T$ is absolutely continuous w.r.t the two dimensional lebesgue measure.
How to connect these two to find that the distribution function of the sum of the convolution of the densities ? Can somebody give me a formal derivation of this ? My text book and my online research provide poor explanations.
EDIT : Tentative of formulation :
Let $T(\omega) = (X(\omega),Y(\omega))$ be a random element, and $\phi(x,y) = x+y$ we have :
\begin{align} P\{Z \in B\} &= P\{X+Y \in B\} \\ &=P\{\phi(X,Y) \in B\} \\ &=P\{(X,Y) \in \phi^{-1}(B)\} \\ &=P\{T \in \phi^{-1}(B)\} \\ &= \int_{\phi^{-1}(B)} f_{XY}(x,y)d(x,y) \quad \text{If $PT^{-1}$ is absolutly continuous w.r.t 2D LM} \\ &= \int\int_{\phi^{-1}(B)} f_{XY}(x,y)dx dy \quad \text{Radon-Nikodym} \\ &= \int_{\mathbb{R}}\int_{-\infty}^{B - y} f_{XY}(x,y)dx dy \quad \text{because $(x,y) \in \phi^{-1}(B) \implies x \in B - y\quad \forall y$} \end{align}