Fix geometric distribution.
Place prior on theta_effective instead of theta=4Nu.
Report theta_herm=4Nu, as opposed to 2Nu.
Update sum_out_coals for integer t.
Switch t!!k to geometric s.
Use a C++ function for the geometric density.