Numerical integration of detector response functions via Monte Carlo simulations
Abstract
Calculations of detector response functions are complicated because they include the intricacies of signal creation from the detector itself as well as a complex interplay between the detector, the particleemitting target, and the entire experimental environment. As such, these functions are typically only accessible through timeconsuming Monte Carlo simulations. Furthermore, the output of thousands of Monte Carlo simulations can be necessary in order to extract a physics result from a single experiment. Here we describe a method to obtain a full description of the detector response function using Monte Carlo simulations. We also show that a response function calculated in this way can be used to create Monte Carlo simulation output spectra a factor of ~1000× faster than running a new Monte Carlo simulation. A detailed discussion of the proper treatment of uncertainties when using this and other similar methods is provided as well. Here, this method is demonstrated and tested using simulated data from the ChiNu experiment, which measures prompt fission neutron spectra at the Los Alamos Neutron Science Center.
 Authors:
 Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
 Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
 Publication Date:
 Research Org.:
 Los Alamos National Lab. (LANL), Los Alamos, NM (United States)
 Sponsoring Org.:
 USDOE
 OSTI Identifier:
 1364576
 Report Number(s):
 LAUR1722123
Journal ID: ISSN 01689002; TRN: US1702526
 Grant/Contract Number:
 AC5206NA25396
 Resource Type:
 Journal Article: Accepted Manuscript
 Journal Name:
 Nuclear Instruments and Methods in Physics Research. Section A, Accelerators, Spectrometers, Detectors and Associated Equipment
 Additional Journal Information:
 Journal Volume: 866; Journal ID: ISSN 01689002
 Publisher:
 Elsevier
 Country of Publication:
 United States
 Language:
 English
 Subject:
 46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND TECHNOLOGY; Detector Response, MCNP, lithiumglass detector
Citation Formats
Kelly, Keegan John, O'Donnell, John M., Gomez, Jaime A., Taddeucci, Terry Nicholas, Devlin, Matthew James, Haight, Robert Cameron, White, Morgan Curtis, Mosby, Shea Morgan, Neudecker, D., Buckner, Matthew Quinn, Wu, Ching Yen, and Lee, Hye Young. Numerical integration of detector response functions via Monte Carlo simulations. United States: N. p., 2017.
Web. doi:10.1016/j.nima.2017.05.048.
Kelly, Keegan John, O'Donnell, John M., Gomez, Jaime A., Taddeucci, Terry Nicholas, Devlin, Matthew James, Haight, Robert Cameron, White, Morgan Curtis, Mosby, Shea Morgan, Neudecker, D., Buckner, Matthew Quinn, Wu, Ching Yen, & Lee, Hye Young. Numerical integration of detector response functions via Monte Carlo simulations. United States. doi:10.1016/j.nima.2017.05.048.
Kelly, Keegan John, O'Donnell, John M., Gomez, Jaime A., Taddeucci, Terry Nicholas, Devlin, Matthew James, Haight, Robert Cameron, White, Morgan Curtis, Mosby, Shea Morgan, Neudecker, D., Buckner, Matthew Quinn, Wu, Ching Yen, and Lee, Hye Young. Tue .
"Numerical integration of detector response functions via Monte Carlo simulations". United States.
doi:10.1016/j.nima.2017.05.048.
@article{osti_1364576,
title = {Numerical integration of detector response functions via Monte Carlo simulations},
author = {Kelly, Keegan John and O'Donnell, John M. and Gomez, Jaime A. and Taddeucci, Terry Nicholas and Devlin, Matthew James and Haight, Robert Cameron and White, Morgan Curtis and Mosby, Shea Morgan and Neudecker, D. and Buckner, Matthew Quinn and Wu, Ching Yen and Lee, Hye Young},
abstractNote = {Calculations of detector response functions are complicated because they include the intricacies of signal creation from the detector itself as well as a complex interplay between the detector, the particleemitting target, and the entire experimental environment. As such, these functions are typically only accessible through timeconsuming Monte Carlo simulations. Furthermore, the output of thousands of Monte Carlo simulations can be necessary in order to extract a physics result from a single experiment. Here we describe a method to obtain a full description of the detector response function using Monte Carlo simulations. We also show that a response function calculated in this way can be used to create Monte Carlo simulation output spectra a factor of ~1000× faster than running a new Monte Carlo simulation. A detailed discussion of the proper treatment of uncertainties when using this and other similar methods is provided as well. Here, this method is demonstrated and tested using simulated data from the ChiNu experiment, which measures prompt fission neutron spectra at the Los Alamos Neutron Science Center.},
doi = {10.1016/j.nima.2017.05.048},
journal = {Nuclear Instruments and Methods in Physics Research. Section A, Accelerators, Spectrometers, Detectors and Associated Equipment},
number = ,
volume = 866,
place = {United States},
year = {Tue Jun 13 00:00:00 EDT 2017},
month = {Tue Jun 13 00:00:00 EDT 2017}
}

A methodology for obtaining empirical curves relating absolute measured scintillation light output to beta energy deposited is presented. Output signals were measured from thin plastic scintillator using NIST traceable beta and gamma sources and MCNP5 was used to model the energy deposition from each source. Combining the experimental and calculated results gives the desired empirical relationships. To validate, the sensitivity of a beryllium/scintillatorlayer neutron activation detector was predicted and then exposed to a known neutron fluence from a DeuteriumDeuterium fusion plasma (DD). The predicted and the measured sensitivity were in statistical agreement.

Monte Carlo path integration in real time via complex coordinates
A method is described for Monte Carlo path integration that is valid for real time propagation. More specifically, it is shown how matrix elements of the complextime propagator e/sup //sup ..beta..//sup >//sup c/H$, ..beta../sub c/ = ..beta../2+it/hdashbar, can be evaluated by straightforward Monte Carlo for values t>>..beta..hdashbar/2. The key feature is that one distorts the path of the integration variables so that the kinetic energy part of the integrand is real. This in turn means that the paths are complex valued, but it is shown that, at least for barriertype potentials, this causes no difficulties. 
Numerical integration of the Langevin equation: Monte Carlo simulation
Monte Carlo simulation techniques are derived for solving the ordinary Langevin equation of motion for a Brownian particle in the presence of an external force. These methods allow considerable freedom in selecting the size of the time step, which is restricted only by the rate of change in the external force. This approach is extended to the generalized Langevin equation which uses a memory function in the friction force term. General simulation techniques are derived which are independent of the form of the memory function. A special method requiring less storage space is presented for the case of the exponentialmore » 
Monte Carlo simulations of hadronic fragmentation functions using the NambuJonaLasiniojet model
The recently developed NambuJonaLasiniojet model is used as an effective chiral quark theory to calculate the quark fragmentation functions to pions, kaons, nucleons, and antinucleons. The effects of the vector mesons {rho}, K{sup *}, and {phi} on the production of secondary pions and kaons are included. The fragmentation processes to nucleons and antinucleons are described by using the quarkdiquark picture, which has been shown to give a reasonable description of quark distribution functions. We incorporate effects of nexttoleading order in the Q{sup 2} evolution, and compare our results with the empirical fragmentation functions. 
An improved multilevel Monte Carlo method for estimating probability distribution functions in stochastic oil reservoir simulations
In this paper, we develop an improved multilevel Monte Carlo (MLMC) method for estimating cumulative distribution functions (CDFs) of a quantity of interest, coming from numerical approximation of largescale stochastic subsurface simulations. Compared with Monte Carlo (MC) methods, that require a significantly large number of highfidelity model executions to achieve a prescribed accuracy when computing statistical expectations, MLMC methods were originally proposed to significantly reduce the computational cost with the use of multifidelity approximations. The improved performance of the MLMC methods depends strongly on the decay of the variance of the integrand as the level increases. However, the main challengemore »