Event generator for<mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mi>J</mml:mi><mml:mo>/</mml:mo><mml:mi>ψ</mml:mi></mml:math>and<mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mi>ψ</mml:mi><mml:mn/><mml:mo>(</mml:mo><mml:mn>2</mml:mn><mml:mi>S</mml:mi><mml:mo>)</mml:mo><mml:mn/></mml:math>decay

We have developed a Monte Carlo generator for simulating charmonium $J/\ensuremath{\psi}$ and $\ensuremath{\psi}(2S)$ inclusive decay. In the model, charmonium decay via gluons is described by the QCD partonic theory, and the partonic hadronization is handled by the LUND model. Extended C- and G-parity conservation are assumed and abnormal suppression effects of charmonium decay are included. This model reproduces the properties of hadronic events in the charmonium inclusive decay, such as the branching ratios of hadronic resonance, the ratios of stable hadrons and the radiative products, and as the global properties of hadronic events.