add functions to calculate and output binary collision density
add functions to output nuclear thickness function and the code is passed a test
add output functions for rho_binary and TA, TB in the ebe_jet mode
fix a bug in the collision probability for gaussian nucleon, width = sqrt(2*B)