Abstract: The accurate determination of the response function of organic scintillation neutron detectors complements their experimental characterization. Monte Carlo simulations with GEANT4 can reduce the effort and cost implied, especially for complex detection systems for which the characterization is more challenging. Previous studies have reported on the inaccuracy of GEANT4 in the calculation of the neutron response of organic scintillation detectors above 6 MeV, due to an incomplete description of the neutron-induced alpha production reactions on carbon. We have improved GEANT4 in this direction by incorporating models and data from NRESP, an excellent Monte Carlo simulation tool developed at the Physikalisch-Technische Bundesanstalt (PTB), Germany, for the specific purpose of calculating the neutron response function of organic scintillation detectors. The results have been verified against simulations with NRESP and validated against Time-Of-Flight measurements with an NE213 detector at PTB. This work has potential applications beyond organic scintillation detectors, to other types of detectors where reactions induced by fast neutrons on carbon require an accurate description.