Create a Poisson-gamma conjugate sampler/step method

[brandonwillard]

In cases for which we are only estimating Poisson rate parameters, it would be nice to have a simple conjugate sampler when we set gamma priors (on the rate parameters). Such a conjugate sampler would follow essentially the same pattern as the conjugate categorical-Dirichlet step method TransMatConjugateStep.