gamma_dist function

Create a gamma distribution object for energy testing