Hi,
I think the problem you observed could be caused by problem generating the 4-D PDF. I would not be surprised if some regions are missed when generating. You could try to increase the foam samples,
see A bug ?, extract a 2-D histogram based on Roofit
createProjection() computes instead the numerical integral first, so I would think it should work better
Lorenzo