A simple geometric model which was motivated by the interesting property of negative binomial multiplicity distribution is used to calculate the pseudo-rapidity distribution of hadron-hadron collisions for all charged multiplicity and for various multiplicity intervals at collider energies. The model with no adjustable parameter yields compatible results with the experimental data at the centre-mass-system energy of square root s=200-900 GeV and can be easily extrapolated to higher energy regions and used in super-high energy cosmic ray studies.