Let $\{f_k\}$ be a sequence of entire functions that are real valued on the real-line. We study the expected number of real zeros of random sums of the form $P_n(z)=\sum_{k=0}^nη_k f_k(z)$, where $\{η_k\}$ are real valued i.i.d.~random variables. We establish a formula for the density function $ρ_n$ for the expected number of real zeros of $P_n$. As a corollary, taking the random variables $\{η_k\}$ to be i.i.d.~standard Gaussian, appealing to Fourier inversion we recover the representation for the density function previously given by Vanderbei through means of a different proof. Placing the restrictions on the common characteristic function $φ$ of $\{η_k\}$ that $|φ(s)|\leq (1+as^2)^{-q}$, with $a>0$ and $q\geq 1$, as well as that $φ$ is three times differentiable with each the second and third derivatives being uniformly bounded, we achieve an upper bound on the density function $ρ_n$ with explicit constants that depend only on the restrictions on $φ$. As an application we considered the limiting value of $ρ_n$ when the spanning functions $f_k(z)=p_k(z)$, $k=0,1,\dots, n$, where $\{p_k\}$ are Bergman polynomials on the unit disk.