I'm looking for the distribution of the sum of independent beta-Bernoulli random variables, that is Bernoulli random variables with beta-distributed parameters, x_i ~ Bern(p_i) with p_i ~ beta(alpha,beta) (same alpha and beta, of course). I guess one can work with convolutions or characteristic...