When computing a recovered model from gravity survey, I’m struggling to quantitatively analyse my results.

For example in the Sparse Norm Gravity Inversion example on the SimPEG documentation, the recovered model looks like a good fit to the true model but I am not sure how to quantitatively process the data

Is it possible to find the index or coordinates of the maximum density values on this recovered model, in the example above, and compare it to the true model to analyse the quality of the inversion?

Thanks