Created a PoissonGammaGaussianLikelihood function