The paper proposes a method to combine QCD matrix elements and parton showers in Monte Carlo simulations of hadronic final states in $e^+e^-$ annihilation. The approach separates the domains of applicability of matrix elements and parton showers at a jet resolution value $y_{\text{ini}}$, defined by the $k_T$-clustering algorithm. Matrix elements are modified by Sudakov form factors, while parton showers are subjected to a veto procedure to cancel the $y_{\text{ini}}$ dependence to next-to-leading logarithmic accuracy. This method provides a leading-order description of hard multi-jet configurations and jet fragmentation, avoiding double counting issues. The authors present initial results using the event generator APACIC++, showing good agreement with experimental data and highlighting the need for further improvements such as including next-to-leading order matrix elements.The paper proposes a method to combine QCD matrix elements and parton showers in Monte Carlo simulations of hadronic final states in $e^+e^-$ annihilation. The approach separates the domains of applicability of matrix elements and parton showers at a jet resolution value $y_{\text{ini}}$, defined by the $k_T$-clustering algorithm. Matrix elements are modified by Sudakov form factors, while parton showers are subjected to a veto procedure to cancel the $y_{\text{ini}}$ dependence to next-to-leading logarithmic accuracy. This method provides a leading-order description of hard multi-jet configurations and jet fragmentation, avoiding double counting issues. The authors present initial results using the event generator APACIC++, showing good agreement with experimental data and highlighting the need for further improvements such as including next-to-leading order matrix elements.