I previously treated the sum of r.v. as the trivial sum as functions: $$ (X+Y)(\omega):=X(\omega)+Y(\omega).$$ But a confusion arises when I encountered the Strong Large Number thm.:
Let $(X_n)_(n\geq 1)$ be independent and identically distributed (i.i.d.) and defined on the same space. Let $\mu=E{X_i}$ and $\sigma^2=\sigma_{X_j}^2 \lt \infty$. Let $S_n=\sum^n_{j=1} X_j$. Then $\lim_{n \to \infty}\frac{S_n}{n} = \lim_{n \to \infty} \frac{1}{n}\sum_{j=1}^{n}X_j=\mu$ a.s and in $L^2$.
Since it says $X_i$ are i.i.d, they must be defined on the same state sapce. But I think the sum $S_n$ here cannot be viewed as the ordinary sum of functions. For example, let the state space $\Omega$ be $\{H, L\}$ where H means head and L means tail in experiments of tossing coin and $X_i(H) := 1, X_i(T) :=0$ for $i\geq 1$. If we treat the sum as ordinary sum of functions, there will be only two values in the range of $S_n$: $S_n(H) = \sum_{i=1}^{n}X_i(H)=n$ and $S_n(T)=\sum_{i=1}^{n}X_i=0$. Then, there is nothing interesting to study.
So, the proper way to define the sum may be extending the domain of $S_n$ from $\Omega$ to $\prod_{i=1}^{\infty}\Omega$, the sigma algebra from $\mathcal{A}$ to $\otimes_{i=1}^{\infty}\mathcal{A}$ and the measure from $P$ to $\otimes_{i=1}^{\infty}P$. In this case, $S_n(\omega)=S_n((\omega_1,\dots,\omega_n,\omega_{n+1},\dots)):=\sum_{i=1}^{n}X_i(\omega_i)$.
I've seen that the well known properties concerning expectation (i.e. integration) such as $E\{X+Y\}=E\{X\}+E\{Y\}$ also hold by the Fubini thm. from this perspective.
But I failed to find similar contents in all resources accessible to me. So I'm wondering if I'm correct?