What is the distribution of the dot product of a Dirichlet vector with a fixed vector? I am trying to get the distribution of a weighted sum when the weights are uncertain:
$S = \sum\limits_{i=1}^N w_iC_i = \mathbf{w}\cdot \mathbf{C}$ where vector $\mathbf{w}$ is random with components having an N-dimensional Dirichlet distribution,: $\mathbf{w} \sim \mathcal{D}_{\theta_1,\theta_2...\theta_n}$ such that $\sum\limits_{i=1}^N w_i = 1$
The vector $\mathbf{C}$ is an N-dimensional fixed vector, whose components are the terms being randomly weighted.
I think that I can approximate the dirichlet by a multivariate gaussian, with variance-covaraince matrix determined from the dirichlet variances-covariances. Then the weighted sum would could be modeled as the a truncated normal distribution.
However, is there any theory out there about the actual distribution of the above operation?

I'm very interested in this distribution as well, unfortunately I believe TenaliRaman's answer is only an approximation. Perhaps it only holds when the distribution is sufficiently concentrated to be approximately Gaussian..?
I tried the following in R:
require(gtools)
library(MASS)

n <- 7000
C <- c(0,1,0.2,0.5)
alpha <- rep(0.2,length(C))
x <- rdirichlet(n, alpha)%*%C


fits <- fitdistr(x,"beta",list(shape1=1,shape2=1))

y <- rbeta(n,shape1=fits$estimate[1],shape2=fits$estimate[2])

plot(sort(x),sort(y))
abline(a=0,b=1)

Here's the result:

Clearly not a match.
