RE: [fluka-discuss]: Age of the particles

From: Anton Lechner <Anton.Lechner_at_cern.ch>
Date: Thu, 6 Feb 2014 15:02:00 +0000

Dear Ben,

1) I tried to reproduce your problem with AGESTK, but for me everything works fine, i.e. the random time assigned in the SOURCE routine is correctly printed out in SODRAW. Might it be that you get only zeros because you didn't activate the SOURCE routine by issuing the SOURCE card in your input file? (which you didn't attach) This would explain why you only get zeros ...

2) By default, for any secondary particle step, ATRACK gives you the time passed since the creation of the corresponding primary particle (since AGESTK of a primary is by default zero). If you set AGESTK in the SOURCE routine to a value different from zero, then ATRACK is counted with respect to this time for any secondary particle created by *this* primary. I hope this answers your question.

3) Note that the purpose of SOUEVT is to store the primaries of *one* event, and not to accumulate primaries from different events. If you indeed intend to have more than one primary per event, note that it is not recommended to have excessive numbers in order to leave a maximum amount of space in the stack for secondary particles. However if you want to indeed increase MXSOEV, you can do this in a local copy of (SOUEVT).

Cheers, Anton





________________________________________
From: owner-fluka-discuss_at_mi.infn.it [owner-fluka-discuss_at_mi.infn.it] on behalf of BenLIU [ben.liu_at_pi.infn.it]
Sent: 06 February 2014 13:04
To: fluka-discuss_at_fluka.org
Subject: [fluka-discuss]: Age of the particles

Dear FLUKA users,

        I am trying to add random random time structure in fluka
source.f(AGESTK(NPFLKA)=FLRNDM (XXX)), and of course we need the
primaries' age for the further use. I tried to read the AGESTK in the
mgdraw.f SODRAW entry, but returns zero always. AGESTK is the age of
particle respect to real time zero. Am I doing this in the right way?

       I try it with simple COMMON block to pass the value out in
source.f, but not successful. I am so sure about the users' COMMON BLOCK.

       Another question is the ATRACK, if I have two primaries at
different time, a photon(daughter particle of the second primary)'s
ATRACK value is respect to the start point (loaded into the main
stack)of his father primary particle(the second primary) not the first
one. So if I want the real timestamp of the photon I should also add the
primary age respect to real time zero. Am I correct?

       At last, I also try to read the COMMON block /SOUEVT/, which is
the saved bank of source particles. But there is the limitation of the
numbers of the particles, MXSOEV = 100. Is this parameter able to be
modified?


Ben LIU
Received on Thu Feb 06 2014 - 16:58:50 CET