Issue with setting up pm.Potential("likelihood", logl(param))

At a first glance it seems like you have a positive log likelihood for 2 of your parameters. That looks wrong, but maybe it’s just unnormalized PDFs?