From: Chris Theis (Christian.Theis@cern.ch)
Date: Tue Nov 28 2006 - 01:19:24 CET
Hi Martin,
I think that some time ago a colleague and I have already solved the problem that you are tackling at the moment. I have a bunch of routines (modified mgdraw, source + parser extensions) that allow you to flag geometry regions as particle sinks in the input file. All particles entering those regions are stopped and their type, weight, energy, position & direction is saved to a file. In a second step these files can be loaded and the particles are re-started for further transport and also the normalization with respect to the original run is taken care of. In case you're interested send me a private message and I will dig those routines up.
Cheers
Chris
----------------------------------------------------------------------------
Chris Theis
CERN/SC-RP - European Organization for Nuclear Research
1211 Geneva 23, Switzerland
Phone: +41 22 767 8069 Office: 892-2A-015
e-mail: Christian.Theis@cern.ch
----------------------------------------------------------------------------
________________________________
Von: owner-fluka-discuss@fisica.unimi.it im Auftrag von Holbourn, MP (Martin)
Gesendet: Mo 27.11.2006 18:27
An: Fluka-Discuss (E-mail)
Betreff: USERDUMP and User written SOURCE routine
Experts,
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,
Regards
Martin Holbourn
Radiation Protection Adviser
Daresbury Laboratory
Daresbury
Warrington
Cheshire
WA4 4AD
Tel: 01925 603266
Fax: 01925 603381
mailto:m.p.holbourn@dl.ac.uk
This archive was generated by hypermail 2.1.6 : Tue Nov 28 2006 - 14:50:59 CET