gamma distribution working
plotting ncoll and npart