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

From: BENLIU <ben.liu_at_pi.infn.it>
Date: Thu, 06 Feb 2014 22:40:03 +0100

Dear Anton,

     Thank you very much!
     I solved the problem with my colleague's help. I made a mistake
when compiling the object files source.o. Now everything works fine.
     Your answers make me much more clear about the time issue now.

  Ciao, Ben

On 02/06/2014 04:02 PM, Anton Lechner wrote:
> 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 - 23:52:56 CET

This archive was generated by hypermail 2.3.0 : Thu Feb 06 2014 - 23:52:57 CET