Purpose: The need for four-dimensional (4D) treatment planning becomes indispensable when it comes to radiation therapy for moving tumors in the thoracic and abdominal regions. The primary purpose of this study is to combine the actual breathing trace during each individual treatment fraction with the Linac's log file information and Monte Carlo 4D dose calculations. We investigated this workflow on multiple computed tomography (CT) datasets in a clinical environment for stereotactic body radiation therapy (SBRT) treatment planning. Methods: We have developed a workflow, which allows us to recalculate absorbed dose to a 4DCT dataset using Monte Carlo calculation methods and accumulate all 4D doses in order to compare them to the planned dose using the Linac's log file, a 4DCT dataset, and the patient's actual breathing curve for each individual fraction. For five lung patients, three-dimensional-conformal radiation therapy (3D-CRT) and volumetric modulated arc treatment (VMAT) treatment plans were generated on four different CT image datasets: a native free-breathing 3DCT, an average intensity projection (AIP) and a maximum intensity projection (MIP) CT both obtained from a 4DCT, and a 3DCT with density overrides based on the 3DCT (DO). The Monte Carlo 4D dose has been calculated on each 4DCT phase using the Linac's log file and the patient's breathing trace as a surrogate for tumor motion and dose was accumulated to the gross tumor volume (GTV) at the 50% breathing phase (end of exhale) using deformable image registration. Results: DD 98% and DD 2% between 4D dose and planned dose differed largely for 3DCT-based planning and also for DO in three patients. Least dose differences between planned and recalculated dose have been found for AIP and MIP treatment planning which both tend to be superior to DO, but the results indicate a dependency on the breathing variability, tumor motion, and size. An interplay effect has not been observed in the small patient cohort. Conclusions: We have developed a workflow which, to our best knowledge, is the first incorporation of the patient breathing trace over the course of all individual treatment fractions with the Linac's log file information and 4D Monte Carlo recalculations of the actual treated dose. Due to the small patient cohort, no clear recommendation on which CT can be used for SBRT treatment planning can be given, but the developed workflow, after adaption for clinical use, could be used to enhance a priori 4D Monte Carlo treatment planning in the future and help with the decision on which CT dataset treatment planning should be carried out.