python code for Bayesian inference of Star Formation Histories, using MCMC