For simulating more accurately neutron or proton production from photonuclear reactions, a data-based photonuclear reaction simulation algorithm has been developed. Reliable photonuclear cross sections from evaluated or experimental database are chosen as input data. For checking the validity of the use of the data-based photonuclear algorithm, benchmarking simulations are presented in detail. We calculate photonuclear cross sections or reaction yield for ~9Be, ~(48)Ti, ~(133)Cs, and ~(197)Au and compare them with experimental data in the region of incident photon energy below *30 MeV. While G...