The Gamma Poisson distribution (GaP) is a mixture model with two positive parameters, α and β.
This hierarchical distribution is used to model a variety of data including failure rates, RNA-Sequencing data  and random distribution of micro-organisms in a food matrix . When a Gamma distribution doesn’t fit data because the overall distribution from multiple samples is too spread out (i.e., the variance is greater than the mean), the Gamma Poisson may be a good fit.
Gamma Poisson Distribution PDF
The probability mass function is :
Where Γ denotes the gamma function.
Connection to the Negative Binomial Distribution
The negative binomial distribution can arise from a gamma mixture of Poisson distributions. Specifically, this can happen when the density distribution across sample plots follows a gamma distribution with scale dependent parameters and a constant mean . It is perhaps for this reason that the GaP process is sometimes called the negative binomial process  although it is more technically a form of the negative binomial rather than its equivalent .
 Holmes, S. & Huber, W. 4 Mixture Models.
 Toft, N. et al. The Gamma–Poisson model as a statistical method to determine if micro-organisms are randomly distributed in a food matrix. Food Microbiology. Volume 23, Issue 1, February 2006, Pages 90-94.
 Gamma-Poisson Distribution. Retrieved December 8, 2021 from: http://www.math.wm.edu/~leemis/chart/UDR/PDFs/Gammapoisson.pdf
 Magnussen, S. (2008). A Gamma-Poisson Distribution of Point to k Nearest Event Distance. Forest Science 54(4): 429-441.
 Ewing, S. & Kunz, M. Fitting of Failure Rate Data to Gamma-Poisson Distribution Utilizing Method of Moments. Retrieved January 8, 2021 from: https://inldigitallibrary.inl.gov/sites/sti/sti/8059210.pdf
 Engelhardt, M. (1994). Events in Time: Basic Analysis of Poisson Data.