Re: (n,gamma) interaction

Good afternoon , Vasilis.

I have triedto use the subroutine that you sent me and have successfully
compiled it , but unfortunately due to luck of my experience I want able to
use this function properly ( actually I wasn't able to use it at all ) . If
it is possible for you to send me some small example of usage of any
function (and the way to collect its results) in the fluka - it would be
great.

Thank you in advance , Julian.

On Wed, May 18, 2011 at 3:32 PM, Vasilis Vlachoudis <
Vasilis.Vlachoudis_at_cern.ch> wrote:

> Dear Julian,
>
> there is no direct solution to identify the reaction, however by inspecting
> the particles stack you can find out if a photon pushes neutrons (after a
> photonuclear reaction). The way I found was to use the routine MDSTCK,
> which is called just after an interaction and before pushing the
> secondaries in the main FLUKA stack (FLKSTK).
> Note the particle that made the interaction is removed from the stack,
> actually only the stack pointer NPFLKA is decreased, and all its
information
> is found in the main stack with pointer NPFLKA+1
>
> Therefore you inspect the if the "last" particle (NPFLKA+1) in the stack to
> check if it was a photon ILOFLK(NPFLKA+1)=7
> and scan the secondary particles from (GENSTK) if any of the them
> KPART(#)=8 is a neutron
>
> SUBROUTINE MDSTCK ( IFLAG, NPSECN )
>
> INCLUDE '(DBLPRC)'
> INCLUDE '(DIMPAR)'
> INCLUDE '(IOUNIT)'
> *
> *----------------------------------------------------------------------*
> * *
> * Copyright (C) 1998-2007 by Alfredo Ferrari & Paola Sala *
> * All Rights Reserved. *
> * *
> * *
> * Created on 19 may 1998 by Alfredo Ferrari & Paola Sala *
> * Infn - Milan *
> * *
> * Last change on 18-jan-07 by Alfredo Ferrari *
> * *
> * *
> * Iflag = 1: standard Kaskad call *
> * = 2: Kaskad call after elastic (but hydrogen) *
> * = 3: Kasneu call *
> * = 4: Emfsco call *
> * *
> *----------------------------------------------------------------------*
> *
> INCLUDE '(EMFSTK)'
> INCLUDE '(FHEAVY)'
> INCLUDE '(FLKSTK)'
> INCLUDE '(GENSTK)'
> INCLUDE '(TRACKR)'
>
> LOGICAL LFIRST
>
> * Peak the last removed particle from the stack
> IF (ILOFLK(NPFLKA+1) .EQ. 7) THEN
> * Scan the secondary stack
> LFIRST = .TRUE.
> DO I=1,NP,1
> * Do we have a neutron there
> IF (KPART(I).EQ.8) THEN
> IF (LFIRST) THEN DO
> LFIRST = .FALSE.
> PRINT *,'Photonnuclear: E=',TKEFLK(NPFLKA+1)
> PRINT *,' Weight=',WTFLK(NPFLKA+1)
> PRINT *,' Pos=',XFLK(NPFLKA+1),
> & YFLK(NPFLKA+1),
> & ZFLK(NPFLKA+1)
> END IF
> PRINT *,' neutron E=',TKI(I), ' wgt=',WEI(I)
> END IF
> END DO
> END IF
> *
> RETURN
> *=== End of subroutine Mdstck =========================================*
> END
>
>
>
> Vasilis Vlachoudis
> Dep EN, CERN
> CH-1211 GENEVA 23
> SWITZERLAND
>
> Phone: +41-22 767 9851
> GSM: +41-76 487 4378
> Fax: +41-22 766 9644
>
>
> ________________________________________
> From: owner-fluka-discuss_at_mi.infn.it [owner-fluka-discuss_at_mi.infn.it] on
> behalf of Julian Rassolov [julian.rassolov_at_gmail.com]
> Sent: 16 May 2011 19:50
> To: fluka-discuss_at_fluka.org
> Subject: (n,gamma) interaction
>
> Good afternoon, everyone.
>
> I am facing some problem here, I need to track amount of the gamma-n
> reaction that took place in some material. How is it possible to define?
>
>
> Thank you in advance,
>
> Julian.
>
>
>

--001636833f825a784204a4a2cc02--
Received on Wed Jun 01 2011 - 15:26:43 CEST