All cross sections of proton induced reactions, angular distributions, energy spectra and double differential cross sections of neutron, proton, deuteron, triton, helium and alpha-particle emissions for p+209Bi reaction are consistently calculated and analyzed at incident proton energies below 250MeV. The optical model, the intra-nuclear cascade model, direct reaction theories, the unified Hauser–Feshbach and exciton model which includes the improved Iwamoto–Harada model are used. Theoretically calculated results are compared with the existing experimental data.