Monte Carlo calculation of the energy loss spectra for gamma rays in cylindrical NaI(Tl) crystals