Implement (x,y) scans of NLL(x,y,theta) for each theta slice
Implement (x,y) marginalization at each theta