Thanks John,
It's useful to know I can generate the material files in PfE too.
The Penelope.in file is very similar to the Penepma one (Penepma is just a version of Penelope configured specifically for EPMA setups whereas Penelope is the full generic version that can be used for e.g. model the radiation dose received by different body tissues - I think this is what it was originally created for). In the Penelope.in file the main difference is how the detectors need to be defined using IMPDET, IDBODY and IDKPAR:
TITLE AWE shuttle loaded with 8 4x10x10mm Pu98Ga samples
>>>>>>>> Source definition.
SKPAR 1 [Primary particles: 1=electron, 2=photon, 3=positron]
SENERG 2.5e4 [Initial energy (monoenergetic sources only)]
SPOSIT 0.0 0.0 1.0 [Coordinates of the source]
SCONE 180 0 0 [Conical beam; angles in deg]
.
>>>>>>>> Material data and simulation parameters.
Up to MAXMAT materials; 2 lines for each material.
MFNAME Pu98Ga.mat [Material file, up to 20 chars]
MSIMPA 1e3 1e3 1e3 0.2 0.2 1e3 1e3 [EABS(1:3),C1,C2,WCC,WCR]
MFNAME epoxy.mat [Material file, up to 20 chars]
MSIMPA 1e3 1e3 1e3 0.2 0.02 1e3 1e3 [EABS(1:3),C1,C2,WCC,WCR]
MFNAME Fe.mat [Material file, up to 20 chars]
MSIMPA 1e3 1e3 1e3 0.2 0.02 1e3 1e3 [EABS(1:3),C1,C2,WCC,WCR]
.
>>>>>>>> Geometry and local simulation parameters.
GEOMFN shuttle.geo [Geometry file, up to 20 chars]
DSMAX 1 1.0e-4 [IB, Maximum step length (cm) in body IB]
DSMAX 2 1.0e-4 [IB, Maximum step length (cm) in body IB]
.
>>>>>>>> Interaction forcing.
IFORCE 1 1 4 -5 0.9 1.0 [KB,KPAR,ICOL,FORCER,WLOW,WHIG]
IFORCE 1 1 5 -5 0.9 1.0 [KB,KPAR,ICOL,FORCER,WLOW,WHIG]
IFORCE 1 2 2 10 0.9 1.0 [KB,KPAR,ICOL,FORCER,WLOW,WHIG]
IFORCE 1 2 2 10 0.9 1.0 [KB,KPAR,ICOL,FORCER,WLOW,WHIG]
IFORCE 2 1 4 -5 0.9 1.0 [KB,KPAR,ICOL,FORCER,WLOW,WHIG]
IFORCE 2 1 5 -5 0.9 1.0 [KB,KPAR,ICOL,FORCER,WLOW,WHIG]
IFORCE 2 2 2 10 0.9 1.0 [KB,KPAR,ICOL,FORCER,WLOW,WHIG]
IFORCE 2 2 2 10 0.9 1.0 [KB,KPAR,ICOL,FORCER,WLOW,WHIG]
.
>>>>>>>> Bremsstrahlung splitting.
IBRSPL 1 5 [KB,splitting factor]
IBRSPL 2 5 [KB,splitting factor]
.
>>>>>>>> X-ray splitting.
IXRSPL 1 5 [KB,splitting factor]
IXRSPL 2 5 [KB,splitting factor]
.
>>>>>>>> Emerging particles. Energy and angular distributions.
NBE 1e3 1e5 200 [Energy window and no. of bins]
NBANGL 60 60 [No. of bins for the angles THETA and PHI]
.
>>>>>>>> Impact detectors (up to 25 different detectors).
IPSF=0; no psf is created.
IPSF=1; a psf is created (for only one detector).
IDCUT=0; tracking is discontinued at the detector entrance.
IDCUT=1; the detector does not affect the tracking.
IDCUT=2; the detector does not affect tracking, the energy
distribution of particle fluence (integrated over the
volume of the detector) is calculated.
IMPDET 0 2e4 1000 0 0 [E-window, no. of bins, IPSF, IDCUT]
IDBODY 11 [Active body]
IDKPAR 2 [Kind of detected particles]
.
>>>>>>>> Job properties.
RESUME dump1.dmp [Resume from this dump file, 20 chars]
DUMPTO dump1.dmp [Generate this dump file, 20 chars]
DUMPP 60 [Dumping period, in sec]
.
NSIMSH 1e9 [Desired number of simulated showers]
TIME 1e9 [Allotted simulation time, in sec]
END [Ends the reading of input data]
The IDBODY links the detector to a body/module in the .geo file that specifies the shape and position of the detector so, for example, the standard annular Penepma detector is modelled as the intersection between two concentric cone surfaces, two concentric spheres, and the z=0 plane:
SURFACE ( CD1) Plane Z=0
INDICES=( 0, 0, 0, 1, 0)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( CD2) Cone, interior
INDICES=( 1, 1,-1, 0, 0)
X-SCALE=( 1.732050000000000E+00, 0)
Y-SCALE=( 1.732050000000000E+00, 0)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( CD3) Cone, exterior
INDICES=( 1, 1,-1, 0, 0)
X-SCALE=( 5.773502000000000E-01, 0)
Y-SCALE=( 5.773502000000000E-01, 0)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( CD4) Inner sphere
INDICES=( 1, 1, 1, 0,-1)
X-SCALE=( 1.000000000000000E+01, 0)
Y-SCALE=( 1.000000000000000E+01, 0)
Z-SCALE=( 1.000000000000000E+01, 0)
0000000000000000000000000000000000000000000000000000000000000000
SURFACE ( CD5) Inner sphere
INDICES=( 1, 1, 1, 0,-1)
X-SCALE=( 1.010000000000000E+01, 0)
Y-SCALE=( 1.010000000000000E+01, 0)
Z-SCALE=( 1.010000000000000E+01, 0)
0000000000000000000000000000000000000000000000000000000000000000
BODY ( CD1) Cylindrical detector 1
MATERIAL( 3)
SURFACE ( CD1), SIDE POINTER=( 1)
SURFACE ( CD2), SIDE POINTER=(-1)
SURFACE ( CD3), SIDE POINTER=( 1)
SURFACE ( CD4), SIDE POINTER=( 1)
SURFACE ( CD5), SIDE POINTER=(-1)
(Thanks to Xavier for this - It would have taken me hours to work this out). All the other bodies that make up the sample are described exactly as for Penepma.
Mike