Abstract. Knowledge about muon tomography has spread in recent years in the geoscientific community and several collaborations between geologists and physicists have been founded. As the data analysis is still mostly done by particle physicists, we address the need of the geoscientific community to participate in the data analysis, while not having to worry too much about the particle physics equations in the background. The result hereof is SMAUG, a toolbox consisting of several modules that cover the various aspects of data analysis in a muon tomographic experiment. In this study we show how a comprehensive geophysical model can be built from basic physics equations. The emerging uncertainties are dealt with by a probabilistic formulation of the inverse problem, which is finally solved by a Monte Carlo Markov Chain algorithm. Finally, we benchmark the SMAUG results against those of a recent study, which however, have been established with an approach that is not easily accessible to the geoscientific community. We show that they reach identical results with the same level of accuracy and precision.