Hadron therapy centres are evolving towards reduced footprint layouts, often featuring a single treatment room. The evaluation of beam properties, radiation protection quantities, and concrete shielding activation via numerical simulations poses new challenges that can be tackled using the numerical beam transport and Monte-Carlo code Beam Delivery Simulation (BDSIM), allowing a seamless simulation of the dynamics as a whole. Specific developments have been carried out in BDSIM to advance its efficiency toward such applications, and a detailed 4D Monte-Carlo scoring mechanism has been implemented. It produces tal- lies such as the spatial-energy differential fluence in arbitrary scoring meshes. The feature makes use of the generic boost::histogram library and allows an event-by-event serialisation and storage in the ROOT data format. The pyg4ometry library is extended to improve the visualisa- tion of critical features such as the complex geometries of BDSIM models, the beam tracks, and the scored quantities. Data are converted from Geant4 and ROOT to a 3D visualisation using the VTK framework. These features are applied to a complete IBA Proteus One model.