MARTINI for heavy quarks: initialization, parton evolution, and hadronization with event generation

We present additions to the MARTINI event generator for examining heavy quarks and quarkonia in heavy-ion collisions. All stages of a heavy-ion collision affect the observables associated with heavy quarks: the initial phase space of the heavy quarks are sampled with PYTHIA8.1, the heavy quarks are evolved using Langevin dynamics and a 3+1-dimensional hydrodynamical description of the heavy-ion collision, and are fragmented and hadronized using a modified version of the color evaporation model that takes into account non-trivial evolution in position space, as well as the possibility of recombinant quarkonium production in heavy-ion collisions. We use this to re-examine the production of quarkonium at RHIC, and anticipating vertex detection we predict yields of B_c mesons at RHIC and the LHC.