PENELOPE 2003 - OECD Nuclear Energy Agency
PENELOPE 2003 - OECD Nuclear Energy Agency
PENELOPE 2003 - OECD Nuclear Energy Agency
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
6.2. Examples of MAIN programs 197<br />
Note: The labels KL, KC that identify each body are defined by the ordering in the<br />
input geometry list. These labels are written on the output geometry report.<br />
STHICK ... For an external source, thickness (height) of the active volume of the<br />
source (cylinder or ring).<br />
-- Default: STHICK=0.0<br />
SRADII ... inner and outer radii of the active source volume.<br />
-- Defaults: SRIN=0.0, SROUT=0.0<br />
SPOSIT ... coordinates of the centre of the source volume.<br />
-- Defaults: SX0=SY0=0, SZ0=-1.0E15<br />
SDIREC ... polar angle θ and azimuthal angle φ of the source beam axis direction, in<br />
deg.<br />
-- Defaults: STHETA=0.0, SPHI=0.0<br />
SAPERT ... angular aperture α of the source beam, in deg.<br />
-- Default: SALPHA=0.0<br />
NMAT ... number of different materials (up to 10 with the original program dimensions).<br />
Materials are identified by their ordering in penelope’s input material<br />
data file.<br />
-- Default: NMAT=1<br />
SIMPAR ... set of simulation parameters for the M-th material; absorption energies,<br />
EABS(1:3,M), elastic scattering parameters, C1(M) and C2(M), and cutoff energy<br />
losses for inelastic collisions and bremsstrahlung emission, WCC(M) and WCR(M).<br />
One line for each material.<br />
-- Defaults: EABS(1,M)=EABS(3,M)=0.01*EPMAX, EABS(2,M)=0.001*EPMAX<br />
C1(M)=C2(M)=0.1, WCC=EABS(1,M), WCR=EABS(2,M)<br />
EPMAX is the maximum energy of all particles found in the simulation (depends on<br />
the source energies).<br />
PFNAME ... name of penelope’s input material data file (18 characters).<br />
-- Default: ’material.mat’<br />
NBE<br />
... number of energy bins and limits of the interval where energy distributions<br />
are tallied.<br />
-- Defaults: NBE=100, EMIN=0.0, EMAX=EPMAX<br />
NBTH ... number of bins for the polar angle θ.<br />
-- Default: NBTH=90<br />
NBPH ... number of bins for the azimuthal angle φ.<br />
-- Default: NBPH=1 (azimuthal average).<br />
NBZ<br />
... number of bins for the z-coordinate.<br />
-- Default: NBZ=100<br />
NBR ... number of bins for the radial variable, r = (x 2 + y 2 ) 1/2 .<br />
-- Default: NBR=100