Add in poisson likelihood

Get this from the xray PE project