Monte Carlo simulation of a plastic scintillator response function in β–γ coincidence measurement