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

Full-Text HTML XML Download Download as PDF (Size:1182KB) PP. 293-306
DOI: 10.4236/am.2017.83024    433 Downloads   562 Views  


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.

Cite this paper

Chin, M. (2017) History-by-History Variance in Monte Carlo Simulation of Radiation Interactions with Matter. Applied Mathematics, 8, 293-306. doi: 10.4236/am.2017.83024.


[1] Ferrari, A., Sala, P.R., Fassò, A. and Ranft, J. (2005) FLUKA: A Multi-Particle Transport Code. Technical Report CERN-2005-10, INFN/TC 05/11, SLAC-R-773, CERN, INFN, SLAC.
[2] Goorley, T., James, M., Booth, T., Brown, F., Bull, J., Cox, L.J., Durkee, J., Elson, J., Fensin, M., Forster, R.A., Hendricks, J., Hughes, H.G., Johns, R., Kiedrowski, B., Martz, R., Mashnik, S., McKinney, G., Pelowitz, D., Prael, R., Sweezy, J., Waters, L., Wilcox, T. and Zukaitis, T. (2016) Features of MCNP6. Annals of Nuclear Energy, 87, 772-783.
[3] Kawrakow, I. (2000) Accurate Condensed History Monte Carlo Simulation of Electron Transport. I. EGSnrc, The New EGS4 Version. Medical Physics, 27, 485-498.
[4] Sempau, J., Fernández-Varea, J.M., Acosta, E. and Salvat, F. (2003) Experimental Benchmarks of the Monte Carlo Code PENELOPE. Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms, 207, 107-123.
[5] Agostinelli, S., Allison, J., Amako, K., Apostolakis, J., Araujo, H., Arce, P., Asai, M., Axen, D., Banerjee, S., Barrand, G., Behner, F., Bellagamba, L., Boudreau, J., Broglia, L., Brunengo, A., Burkhardt, H., Chauvie, S., Chuma, J., Chytracek, R., Cooperman, G., Cosmo, G., Degtyarenko, P., Dell’Acqua, A., Depaola, G., Dietrich, D., Enami, R., Feliciello, A., Ferguson, C., Fesefeldt, H., Folger, G., Foppiano, F., Forti, A., Garelli, S., Giani, S., Giannitrapani, R., Gibin, D., Gómez Cadenas, J.J., González, I., Gracia Abril, G., Greeniaus, G., Greiner, W., Grichine, V., Grossheim, A., Guatelli, S., Gumplinger, P., Hamatsu, R., Hashimoto, K., Hasui, H., Heikkinen, A., Howard, A., Ivanchenko, V., Johnson, A., Jones, F.W., Kallenbach, J., Kanaya, N., Kawabata, M., Kawabata, Y., Kawaguti, M., Kelner, S., Kent, P., Kimura, A., Kodama, T., Kokoulin, R., Kossov, M., Kurashige, H., Lamanna, E., Lampén, T., Lara, V., Lefebure, V., Lei, F., Liendl, M., Lockman, W., Longo, F., Magni, S., Maire, M., Medernach, E., Minamimoto, K., Mora de Freitas, P., Morita, Y., Murakami, K., Nagamatu, M., Nartallo, R., Nieminen, P., Nishimura, T., Ohtsubo, K., Okamura, M., O’Neale, S., Oohata, Y., Paech, K., Perl, J., Pfeiffer, A., Pia, M.G., Ranjard, F., Rybin, A., Sadilov, S., Di Salvo, E., Santin, G., Sasaki, T., Savvas, N., Sawada, Y., Scherer, S., Sei, S., Sirotenko, V., Smith, D., Starkov, N., Stoecker, H., Sulkimo, J., Takahata, M., Tanaka, S., Tcherniaev, E., Safai Tehrani, E., Tropeano, M., Truscott, P., Uno, H., Urban, L., Urban, P., Verderi, M., Walkden, A., Wander, W., Weber, H., Wellisch, J.P., Wenaus, T., Williams, D.C., Wright, D., Yamada, T., Yoshida, H. and Zschiesche, D. (2003) Geant4—A Simulation Toolkit. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment, 506, 250-303.
[6] Battistoni, G., Bauer, J., Boehlen, T.T., Cerutti, F., Chin, M.P.W., Dos Santos Augusto, R., Ferrari, A., Garcia Ortega, P., Koz lowska, W., Magro, G., Mairani, A., Parodi, K., Sala, P.R., Schoofs, P., Tessonnier, T. and Vlachoudis, V. (2016) The FLUKA Code: An Accurate Simulation Tool for Particle Therapy. Frontiers in Oncology, 6, 116.
[7] Battistoni, G., Bohlen, T., Cerutti, F., Chin, M.P.W., Luigi, S.E., Fassò, A., Ferrari, A., Lechner, A., Empl, A., Mairani, A., Mereghetti, A., Garcia Ortega, P., Ranft, J., Roesler, S., Sala, P.R., Vlachoudis, V. and Smirnov, G. (2015) Overview of the FLUKA Code. Annals of Nuclear Energy, 82, 10-18.
[8] Bohlen, T., Cerutti, F., Chin, M.P.W., Fassò, A., Ferrari, A., Garcia Ortega, P., Mairani, A., Sala, P.R., Smirnov, G. and Vlachoudis, V. (2014) The FLUKA Code: Developments and Challenges for High Energy and Medical Applications. Nuclear Data Sheets, 120, 211-214.
[9] Chin, M.P.W. (2015) FTREE: Single-History Monte Carlo Analysis for Radiation Detection and Measurement. Journal of Radioanalytical & Nuclear Chemistry, 306, 565.
[10] Chin, M.P.W., Cerutti, F., Ferrari, A., Mairani, A. and Sala, P.R. (2011) Sample Histories from Carbon Therapy of a Human Brain. Transactions American Nuclear Society, 105, 55-56.
[11] Chin, M.P.W. and Spyrou, N.M. (2007) Event-by-Event Monte Carlo Tracking of Neutron-Nucleus Collisions in Neutron Detectors. Transactions American Nuclear Society, 97, 288-289.
[12] Chin, M.P.W., Cerutti, F., Ferrari, A., Garcia Ortega, P. and Sala, P.R. (2014) Monte Carlo Latching Studies of Prompt-Gamma Detection during Hadron Therapy. IEEE Transactions on Nuclear Science, 61, 2540-2546.
[13] Sempau, J., Sánchez-Reyes, A., Salvat, F., Tahar, H.O., Jiang, S.B. and Fernández-Varea, J.M. (2001) Monte Carlo Simulation of Electron Beams from an Accelerator Head Using PENELOPE. Physics in Medicine and Biology, 46, 1163.
[14] Walters, B.R.B., Kawrakow, I. and Rogers, D.W.O. (2002) History by History Statistical Estimators in the BEAM Code System. Medical Physics, 29, 2745-2752.
[15] Chiaveri, E., Altstadt, S., Andrzejewski, J., Audouin, L., Barbagallo, M., Bécares, V., Becvár, F., Belloni, F., Berthoumieux, E., Billowes, J., Boccone, V., Bosnar, D., Brugger, M., Calviani, M., Calvino, F., Cano-Ott, D., Carrapico, C., Cerutti, F., Chin, M.P.W., Colonna, N., Cortés, G., Cortés-Giraldo, M.A., Diakaki, M., Domingo-Pardo, C., Duran, I., Dressler, R., Dzysiuk, N., Eleftheriadis, C., Ferrari, A., Fraval, K., Ganesan, S., García, A.R., Giubrone, G., Gómez-Hornillos, M.B., Gonialves, I.F., González-Romero, E., Griesmayer, E., Guerrero, C., Gunsing, F., Gurusamy, P., Hernández-Prieto, A., Jenkins, D.G., Jericha, E., Kadi, Y., Kappeler, F., Karadimos, D., Kivel, N., Koehler, P., Kokkoris, M., Krticka, M., Kroll, J., Lampoudis, C., Langer, C., Leal-Cidoncha, E., Lederer, C., Leeb, H., Leong, L.S., Losito, R., Mallick, A., Manousos, A., Marganiec, J., Martínez, T., Massimi, C., Mastinu, P.F., Mastromarco, M., Meaze, M., Mendoza, E., Mengoni, A., Milazzo, P.M., Mingrone, F., Mirea, M., Mondalaers, W., Paradela, C., Pavlik, A., Perkowski, J., Plompen, A., Praena, J., Quesada, J.M., Rauscher, T., Reifarth, R., Riego, A., Robles, M.S., Roman, F., Rubbia, C., Sabaté-Gilarte, M., Sarmento, R., Saxena, A., Schillebeeckx, P., Schmidt, S., Schumann, D., Tagliente, G., Tain, J.L., Tarrío, D., Tassan-Got, L., Tsinganis, A., Valenta, S., Vannini, G., Variale, V., Vaz, P., Ventura, A., Versaci, R., Vermeulen, M.J., Vlachoudis, V., Vlastou, R., Wallner, A., Ware, T., Weigand, M., Weiss, C., Wright, T. and Zugec, P. (2014) The CERN n TOF Facility: Neutron Beams Performances for Cross Section Measurements. Nuclear Data Sheets, 119, 1-4.

comments powered by Disqus

Copyright © 2017 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.