SUMMARY
A pulse height tally response expansion (PHRE) method is developed for detectors. By expanding the incident flux of the detector, a set of response functions is constructed via Monte Carlo estimators for pulse height tallies. B-spline functions are selected for the expansion of the incident flux in the energy phase space. Response functions are generated for a CsI(Na) crystal using an energy dependent B-spline form of the PHRE method. The method is verified for several incident flux spectra on the crystal. Results of the energy dependent B-spline form of the PHRE method are compared to the solutions generated using MCNP5. The method is shown to be several orders faster than MCNP5 while maintaining paralleled accuracy.