Monte Carlo Simulation for Elastic Energy Loss of Hard Partons in a Hydrodynamical Background

We have developed a Monte Carlo simulation describing the $2\to 2$ scatterings of perturbatively produced, non-eikonally propagating high-energy partons with the quarks and gluons of the expanding QCD medium created in ultrarelativistic heavy ion collisions. The partonic scattering rates are computed in leading-order perturbative QCD (pQCD), while three different hydrodynamical scenarios are used to model the strongly interacting medium. We compare our results and tune the model with the neutral pion suppression observed in $\sqrt{s_{NN}}=200$ GeV Au+Au collisions at the BNL-RHIC. We find the incoherent nature of elastic energy loss incompatible with the measured angular dependence of the suppression. The effects of the initial state density fluctuations of the bulk medium are found to be small. Also the extrapolation from RHIC to the LHC is discussed.