EDITED
John is playing a game on $n$ days, each day being independent.
On each day $i$, his probability of success is $p_i$.
We have $\frac{1}{n} \sum_{i=1}^n p_i = p$, and typically, the standard deviation $\sigma$ of these $p_i$ is small. $\sigma$ is known.
So I have a succession of $Bernoulli(p_i)$. I want to model the probability of having k success after n trials, using $p$ and $\sigma$ for instance.
I can approximate this with a Binomial distribution, but even though it's close it still leads to biases. Any thoughts on how to improve this?