The sum I am looking for is the following sum as $M \to \infty$: 

$$  L(\omega)   = \sum_{m = 1}^{M} \frac{\sin\left( N \frac{\omega_m - \omega}{2} \right)}{\sin\left( \frac{\omega_m - \omega}{2} \right)} \cos\left(N \frac{\omega_m - \omega}{2} + \beta_m \right) $$

where
* $\omega_m$ is a random number from a Gaussian distribution having the parameters (mean $\mu$ and variance $\sigma^2$). 
$$ \omega_m \sim \mathcal{N}(\mu, \sigma^2) $$
* the $\beta_m$ are random numbers drawn from a uniform distribution from $-\pi$ to $+\pi$, 

$$ \beta_m \sim \mathcal{U}[-\pi, +\pi] $$