From: Holbourn, MP \(Martin\) (firstname.lastname@example.org)
Date: Mon Nov 27 2006 - 18:27:17 CET
My geometry is similar to that found at most storage rings - a concrete ratchet wall with a beampipe through the end wall. My job runs successfully and I do indeed get high energy photons streaming through the hole in the concrete from electrons lost upstream on a collimator. What I want to do is use these photons as the source to a separate phased run as described in the section of the manual relating to collision tapes. The manual seems to suggest that one can record all particles crossing a given boundary but how do you specify which boundary crossing you are interested in? Option USERDUMP does not seem to contain any info on how to specify a boundary so am I missing something? Should WHAT(3) of USERDUMP be = 1?
Also how does one write subroutine SOURCE to read the collision tape? Page 416 of the manual says for case 3 (source particles) the first record contains
NCASE, LSTACK, LSTMAX, TKESUM, WEIPRI
and the next record
(ILO(I), ETOT(I), WT(I), XA(I), YA(I), ZA(I), TX(I), TY(I), TZ(I), I=1,LSTACK)
BUT P438 of the manual suggests reading the collision tape by the statement
READ(21,*) IPART,X ,Y ,Z, COSX, COSY, COSZ, ENERGY, WEIGHT
which is a different order of variables than written to the tape - again I might have misunderstood what is going on here.
Can anyone point me in the right direction or perhaps share a working source.f routine that does what I am trying to do.
Thanks for any help,
Radiation Protection Adviser
Tel: 01925 603266
Fax: 01925 603381
This archive was generated by hypermail 2.1.6 : Mon Nov 27 2006 - 20:41:52 CET