An iPython-based Monte-Carlo Glauber Model of Nuclear Collisions

An iPython-based package to simulate the Monte-Carlo Glauber model of nuclear interactions was developed and published online for broad distribution. The simulation utilizes compiled data on nuclear charge density distributions to create nuclei and cross-section data from the Particle Data Group to generate large statistics simulations across a broad range of energies and collision systems. The simulation correlates the number of nucleons participating in a collision as well as the number of binary collisions with the impact parameter for each event. Individual collisions can be visually represented, demonstrating the event-by-event variation of particular geometric overlaps, which are obscured in the ensemble data. Good agreement between the program and expected results for Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV is demonstrated. Examples of possible future RHIC experiments such as He$^3$+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV will be shown.