Particle production modeling using negative binomial distribution, NBD, in order to fit data from STAR and CMS. Discussion of convolution (convolving the NBD with the Ncoll distribution from the Glauber Model). Discussion of fitting the NBD parameters using chi-square test.