History-by-History Variance in Monte Carlo Simulation of Radiation Interactions with Matter

• 作者: Mary Pik Wai Chin
• 来源期刊: 应用数学(英文)
• 学科: 医学
• 关键词: Monte Carlo FLUKA Error STATISTICS Variance
• 年，卷（期）: yysxyw_2017,（3）
• 页码范围: 293-306
• 页数: 14页
• 分类号: R73

摘要

Radiation interaction with matter is by nature stochastic. Monte Carlo simulation makes possible, practical, economical and ethical, many experiments otherwise not. In-silico simulation of random samples of radiation histories places into our hands data which may be treated and analysed in radically different ways. This work reports history-by-history computation of the variance in simulation output from FLUKA, a general-purpose Monte Carlo code. Variance computed history by history is compared with variance estimated from varying batch sizes. This work also addresses the issue of under-sampling where, against conventional expectations, the variance spikes as we sample more histories. The background to the spikes is traced by reconstructing single histories interaction by interaction—a novel level of details atypical of Monte Carlo simulations in the field, where statistical convergence of averaged quantities has been the focus.