This paper presents a method for combining QCD matrix elements and parton showers in Monte Carlo simulations of hadronic final states in $ e^{+}e^{-} $ annihilation. The method separates the domains of applicability of matrix elements and parton showers at a jet resolution parameter $ y_{ini} $, defined according to the $ k_T $-clustering algorithm. Matrix elements are modified by Sudakov form factors, and parton showers are subjected to a veto procedure to cancel dependence on $ y_{ini} $ to next-to-leading logarithmic (NLL) accuracy. This approach provides a leading-order description of hard multi-jet configurations and jet fragmentation while avoiding double-counting issues. The method is implemented in the event generator APACIC++ and tested against experimental data. The results show that APACIC++ provides a good description of multi-jet observables and reproduces the overall shape of events. The approach is extended to lepton-hadron and hadron-hadron collisions, and further improvements such as including next-to-leading order matrix elements are suggested to enhance agreement with data. The method is effective in combining matrix elements and parton showers, merging their benefits in a simple way.This paper presents a method for combining QCD matrix elements and parton showers in Monte Carlo simulations of hadronic final states in $ e^{+}e^{-} $ annihilation. The method separates the domains of applicability of matrix elements and parton showers at a jet resolution parameter $ y_{ini} $, defined according to the $ k_T $-clustering algorithm. Matrix elements are modified by Sudakov form factors, and parton showers are subjected to a veto procedure to cancel dependence on $ y_{ini} $ to next-to-leading logarithmic (NLL) accuracy. This approach provides a leading-order description of hard multi-jet configurations and jet fragmentation while avoiding double-counting issues. The method is implemented in the event generator APACIC++ and tested against experimental data. The results show that APACIC++ provides a good description of multi-jet observables and reproduces the overall shape of events. The approach is extended to lepton-hadron and hadron-hadron collisions, and further improvements such as including next-to-leading order matrix elements are suggested to enhance agreement with data. The method is effective in combining matrix elements and parton showers, merging their benefits in a simple way.