The computer code system PENELOPE (version 2008) performs Monte Carlo simulation of coupled
electron-photon transport in arbitrary materials for a wide energy range, from a few hundred eV to
about 1 GeV. Photon transport is simulated by means of the standard, detailed simulation scheme.
Electron and positron histories are generated on the basis of a mixed procedure, which combines
detailed simulation of hard events with condensed simulation of soft interactions. A geometry package
called PENGEOM permits the generation of random electron-photon showers in material systems
consisting of homogeneous bodies limited by quadric surfaces, i.e., planes, spheres, cylinders, etc. This
report is intended not only to serve as a manual of the PENELOPE code system, but also to provide the
user with the necessary information to understand the details of the Monte Carlo algorithm.

/pdf/penelope-2006-a-code-system-for-monte-carlo-simulation-of-43nx3zmw6c.pdf

PENELOPE-2006: A Code System for Monte Carlo Simulation of Electron and Photon Transport

An algorithm for Monte Carlo simulation of coupled electron-photon transport is described. Electron and positron tracks are generated by means of PENELOPE, a mixed procedure developed by Baro et al. [Nucl. Instr. and Meth. B 100 (1995) 31]. The simulation of photon transport follows the conventional, detailed method. Photons are assumed to interact via coherent and incoherent scattering, photoelectric absorption and electron-positron pair production. Photon interactions are simulated through analytical differential cross sections, derived from simple physical models and renormalized to reproduce accurate attenuation coefficients available from the literature. The combined algorithm has been implemented in a FORTRAN 77 computer code that generates electron-photon showers in arbitrary materials for the energy range from ∼1 GeV down to 1 keV or the binding energy of the L-shell of the heaviest element in the medium, whichever is the largest. The code is capable of following secondary particles that are generated within this energy range. The reliability of the algorithm and computer code is demonstrated by comparing simulation results with experimental data and with results from other Monte Carlo codes.

An algorithm for Monte Carlo simulation of coupled electron-photon transport

The physical algorithms implemented in the latest release of the general-purpose Monte Carlo code penelope for the simulation of coupled electron–photon transport are briefly described. We discuss the mixed (class II) scheme used to transport intermediate- and high-energy electrons and positrons and, in particular, the approximations adopted to account for the energy dependence of the interaction cross-sections. The reliability of the simulation code, i.e. of the adopted interaction models and tracking algorithms, is analyzed by means of a comprehensive comparison of simulation results with experimental data available from the literature. The present analysis demonstrates that penelope yields a consistent description of electron transport processes in the energy range from a few keV up to about 1 GeV.

Experimental benchmarks of the Monte Carlo code penelope

Electron-photon Monte Carlo calculations: The ETRAN code

This work involves a validation of the photon and electron transport of the GEANT4 particle simulation toolkit for radiotherapy physics applications. We examine the cross sections and sampling algorithms of the three electromagnetic physics models in version 4.6.1 of the toolkit: Standard, Low-energy, and Penelope. The depth dose distributions in water for incident monoenergetic and clinical beams are compared to the EGSNRC results. In photon beam simulations, all three models agree with EGSNRC to within 2%, except for the buildup region. Larger deviations are found for incident electron beams, and the differences are affected by user-imposed electron step limitations. Particle distributions through thin layers of clinical target materials, and perturbation effects near high-Z and low-Z interfaces are also investigated. The electron step size artifacts observed in our studies indicate potential problems with the condensed history algorithm. A careful selection of physics processes and transport parameters is needed for optimum efficiency and accuracy.

C.A. Quarles

Papers

Shape functions for atomic-field bremsstrahlung from electrons of kinetic energy 1–500 keV on selected neutral atoms 1 ≤ Z ≤ 92