The well-known Monte Carlo simulation package JETSET is developed to simula
te the hadronic final states in e(+) e(-) collision with center-of-mass ene
rgy E-cm greater than or similar to 10GeV, but not suitable for few GeV ene
rgy region. We develop a Monte Carlo program, LUARLW, which can well simula
te two-body up to six-body hadronic final states at BEPC energy region. In
LUARLW,the effects of all gluonic emissions are neglected. The LUARLW produ
ces more than 14 distributions, which agree with the BES data very well. Th
is package is used for estimating hadronic events detection efficiencies in
the measurements of total cross section of hadron production at, BES betwe
en 2-5 GeV.