We investigate the fortune of young neutron stars (NS) in the whole volume of the Milky Way with new code for population synthesis. We start our modeling from the birth of massive OB stars and follow their motion in the Galaxy up to the Supernova explosion. Next we integrate the equations of motion of NS in the averaged gravitational potential of the Galaxy. We estimate the mean kick velocities from a comparison the model Z and R-distributions of radio emitting NS with that for galactic NS accordingly ATNF pulsar catalog. We follow the history of the rotational velocity and the surface magnetic field of NS taking into account the significant magnetic field decay during the first million year of a neutron star's life. The derived value for the mean time of ohmic decay is 2.3ċ105 years. We model the subsample of galactic radio pulsars which can be detected with available radio telescopes, using a radio beaming model with inhomogeneous distribution of the radio emission in the cone. The distributions functions o