  The method to simulate the rescattering and energy loss of hard partons in
ultrarelativistic heavy ion collisions has been developed. The model is a fast
Monte-Carlo tool introduced to modify a standard PYTHIA jet event. The full
heavy ion event is obtained as a superposition of a soft hydro-type state and
hard multi-jets. The model is applied to analysis of the jet quenching pattern
at RHIC.
