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) = n k=0 η 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)| ≤ (1 + as 2 ) −q , with a > 0 and q ≥ 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, . . . , n, where {p k } are Bergman polynomials on the unit disk.