The main objective of the thesis is the implementation of Monte Carlo simulation in ROOT framework which is provided by CERN to use the Glauber Model for making geometrical calculations and estimations of quantities such as the number of participating nucleons N part and binary nucleon-nucleon collisions N coll in High Energy Nuclear Collisions, in our case between Gold (Au) ions. I reached the main objective of this project by programming the Glauber MC simulation to produce results based on theoretical background. This method is more sufficient and applicable than the Optical limit model of Glauber. After running the simulation, I was able to estimate the dataset for values of the quantities N part and N coll under different circumstances (different centrality levels), which are non-measurable data during the experiments of nuclei collisions, compared to the fixed-target experiments. By providing different input parameters to Wood-Saxon distribution and plotting dependence between obtained values and centrality level, I tried to present the best fit which can produce similar results to the theoretical calculations. It is agreed that the given relations between centrality/impact parameter and number of participants/collisions are correct and supported by theory.