Implementing primordial binaries in simulations of star cluster formation with a hybrid MHD and direct N-Body method