Here is an illustration on a Hierarchical Poisson failure rates from Clark and Gelfand, using Python and the PyMC package.

The Python code is as simple as the R code, although it is obviously more object-oriented. The main part are highlighted below. These are in fact where we specify the prior distributions:

Running the model:

Plot one
Plot two