With the Large Hadron Collider (LHC) experiment, the particle physics community is searching for new physics beyond the standard model (SM). For this, theoretical models have to be translated to measurable quantities, observables. Monte Carlo (MC) event generators are the computational tools that simulate the experimental processes that occur in the LHC to make predictions of observables. In this project, models for SM Z bosons and BSM (beyond the SM) Z' bosons will be used for high precision calculations with a new MC event generator that has to be constructed. This new generator will make use of the recently developed parton branching (PB) method for transverse momentum dependent parton distribution functions (TMDs) and will be made available as an analysis tool for experimental and phenomenological collaborations. The TMD formalism is well suited to study non-inclusive observables such as transverse momentum and correlations. The implementation of the PB method in a new event generator is one of the main challenges because the separate parts that form the generator should be combined in a consistent and correct manner. Another challenge is to calculate hard scattering matrix elements for Z and Z' production from benchmark models. With the implementation of these matrix elements in the generator, it is expected to achieve high precision calculations that can contribute to a more detailed search for BSM physics.