StarPy derives the quenching star formation history (SFH) of a single galaxy through the Bayesian Markov Chain Monte Carlo method code emcee (ascl:1303.002). The sample function implements the emcee EnsembleSampler function for the galaxy colors input. Burn-in is run and calculated for the length specified before the sampler is reset and then run for the length of steps specified. StarPy provides the ability to use the look-up tables provided or creating your own.
http://adsabs.harvard.edu/abs/2015MNRAS.450..435S and acknowledgement: This work has made use of the StarPy star formation history Python code described in Smethurst et al. (2015).