The application of Monte Carlo (MC) simulation in the context of Medical Linear Accelerator (LinAc) treatment planning has increased in recent decades owing to its accuracy in dose calculations. The Multileaf Collimator (MLC), a component of the LinAc heads, plays an essential role by allowing the precise customization of the treatment beam’s shape, effectively targeting the tumor area while minimizing radiation exposure to surrounding healthy structures. During a jaws static treatment, this MLC configuration varies for each angular position of the accelerator gantry, while the rest of the LinAc geometry remains static. To address the operation of the accelerator during the treatment planning avoiding the need of computationally intensive and repetitive simulations, phase-space files (PSF) are widely employed in LinAc MC calculations. This study aims to introduce and validate a methodology for the computational generation of PSFs specifically defined just below the MLC. These PSFs will assume the appropriate shape corresponding to the angle at which the accelerator gantry is positioned during treatment. The core of this methodology involves the creation of a comprehensive database containing precalculated probability functions derived from MC. In order to validate the methodology, simulated dose values of a LinAc emitting a 6 MV photon beam have been employed. The results have been tested against experimental data, and compared with the results from complete simulations, yielding PSFs of varying shapes. Finally, this study reveals that the results obtained by comparing organ dose values computed from the generated PSFs to those from the reference PSF are statistically compatible. This innovative method not only offers a valuable tool for optimizing LinAc treatment planning but also underscores the importance of accurate dose calculations in the field of radiation therapy.