Use gaussian likelihood for number counts, it will require a covariance matrix. In a first approach, it can just be estimated from the data.