New: calculate likelihood and loop
Merge branch 'feature_likelihood' into develop