Ultra-violet photoemission spectroscopy is a widely-used experimental technique to investigate the valence electronic structure of surfaces and interfaces. When detecting the intensity of the emitted electrons not only as a function of their kinetic energy, but also depending on their emission angle, as is done in angle-resolved photoemission spectroscopy (ARPES), extremely rich information about the electronic structure of the investigated sample can be extracted. For organic molecules adsorbed as well-oriented ultra-thin films on metallic surfaces, ARPES has evolved into a technique called photoemission tomography (PT). By approximating the final state of the photoemitted electron as a free electron, PT uses the angular dependence of the photocurrent, a so-called momentum map or k-map, and interprets it as the Fourier transform of the initial state's molecular orbital, thereby gains insights into the geometric and electronic structure of organic/metal interfaces.In this contribution, we present kMap.py which is a Python program that enables the user, via a PyQt-based graphical user interface, to simulate photoemission momentum maps of molecular orbitals and to perform a one-toone comparison between simulation and experiment. Based on the plane wave approximation for the final state, simulated momentum maps are computed numerically from a fast Fourier transform (FFT) of real space molecular or- * Corresponding author.