2D histogram interpolating/smoothing in unbinned fit

Is there a function or suggested way to do this? It is quite common to do unbinned fitting > 1D and it seems strange that I cannot find such a function anywhere.

Thanks

Have a look at class TGraph2D and examples in $ROOTSYS/tutorials/fit

Rene