This paper describes the concept, technical realisation and validation of a largely datadriven method to model events with Z→tt decays. In Z→μμ events selected from proton-proton collision data recorded at √s=8 TeV with the ATLAS experiment at the LHC in 2012, the Z decay muons are replaced by t leptons from simulated Z→tt decays at the level of reconstructed tracks and calorimeter cells. The t lepton kinematics are derived from the kinematics of the original muons. Thus, only the well-understood decays of the Z boson and t leptons as well as the detector response to the t decay products are obtained from simulation. All other aspects of the event, such as the Z boson and jet kinematics as well as effects from multiple interactions, are given by the actual data. This so-called t-embedding method is particularly relevant for Higgs boson searches and analyses in tt final states, where Z → tt decays constitute a large irreducible background that cannot be obtained directly from data control samples. In this paper, the relevant concepts are discussed based on the implementation used in the ATLAS Standard Model H ! tt analysis of the full datataset recorded during 2011 and 2012.