From the Gaia mission, kinematic and star count data, together with the physical parameters of the stars - ages and metallicities-, will allow to characterise our Galaxy populations and, from that, the overall Galactic gravitational potential. One of the promising procedures to reach such goal will be to optimise the present Population Synthesis models by fitting, through robust statistical techniques, the large and small scale structure and kinematics parameters that best will reproduce Gaia data. As a first step we present here the comparison between the data from the Tycho-2 catalogue and the Besançon Galaxy Model simulations, as well as our present and future work on the model optimization.