We present a method for parametric modelling of the physical components of the Galaxys magnetised interstellar medium, simulating the observables, and mapping out the likelihood space using a Markov Chain Monte-Carlo analysis. We then demonstrate it using total and polarised synchrotron emission data as well as rotation measures of extragalactic sources. With these three datasets, we define and study three components of the magnetic field: the large-scale coherent field, the small-scale isotropic random field, and the ordered field. In this first paper, we use only data along the Galactic plane and test a simple 2D logarithmic spiral model for the magnetic field that includes a compression and a shearing of the random component giving rise to an ordered component. We demonstrate with simulations that the method can indeed constrain multiple parameters yielding measures of, for example, the ratios of the magnetic field components. Though subject to uncertainties in thermal and cosmic ray electron densities and depending on our particular model parametrisation, our preliminary analysis shows that the coherent component is a small fraction of the total magnetic field and that an ordered component comparable in strength to the isotropic random component is required to explain the polarisation fraction of synchrotron emission. We outline further work to extend this type of analysis to study the magnetic spiral arm structure, the details of the turbulence as well as the 3D structure of the magnetic field.