We benchmark here a Monte Carlo model simulating an inelastic neutron scattering (INS) system for quantitative analysis of carbon in soil. Specifically, we compare the simulations with experimental results of copper foils activations, INS system calibration, INS system optimization of the height above the ground and comparing pulse height distributions due to 137Cs and 60Co sources. Most of the simulations and the measurements agree better than 10%, although some of them registered discrepancies larger than 20%.