Implemented ``get_theta_distr`` function that use specified gamma distribution to compute samples from theta distribution. See Cohen 1994 A15.