Re: (n,gamma) interaction from Vasilis Vlachoudis on 2011-06-06 (fluka discuss archive)

From: Vasilis Vlachoudis <Vasilis.Vlachoudis_at_cern.ch>
Date: Mon, 6 Jun 2011 13:55:22 +0200

Content-Type: multipart/mixed;
boundary="------------010506060005030703000702"

--------------010506060005030703000702
Content-Type: multipart/alternative;
boundary="------------020305040400090108040208"

--------------020305040400090108040208
Content-Type: text/plain; charset="ISO-8859-1"; format=flowed
Content-Transfer-Encoding: 7bit

Hi Julian,

photonuclear interactions are not on by default. Please use the PHOTONUC
card to enable them and possibly the LAM-BIAS to enhance them.
Please find attached an input file that I used to test
The output of the MDSTCK is printed in the .log file

Cheers
Vasilis

On 01/06/11 11:03, Julian Rassolov wrote:
> 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 <mailto: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
> <mailto:owner-fluka-discuss_at_mi.infn.it>
> [owner-fluka-discuss_at_mi.infn.it
> <mailto:owner-fluka-discuss_at_mi.infn.it>] on behalf of Julian
> Rassolov [julian.rassolov_at_gmail.com
> <mailto:julian.rassolov_at_gmail.com>]
> Sent: 16 May 2011 19:50
> To: fluka-discuss_at_fluka.org <mailto: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.
>
>
>

--------------020305040400090108040208
Content-Type: text/html; charset="ISO-8859-1"
Content-Transfer-Encoding: 7bit

<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
 <head>
 <meta content="text/html; charset=ISO-8859-1"
 http-equiv="Content-Type">
 <title></title>
 </head>
 <body bgcolor="#ffffff" text="#000000">
 Hi Julian, 
 
 photonuclear interactions are not on by default. Please use the
 PHOTONUC card to enable them and possibly the LAM-BIAS to enhance
 them. 
 Please find attached an input file that I used to test 
 The output of the MDSTCK is printed in the .log file 
 
 Cheers 
 Vasilis 
 
 On 01/06/11 11:03, Julian Rassolov wrote:
 <blockquote
 cite="mid:BANLkTi=PsqOHtS7FuhD6ERhADv7Uu_9GMg_at_mail.gmail.com"
 type="cite">
 <meta http-equiv="Content-Type" content="text/html;
 charset=ISO-8859-1">
 <div dir="ltr">
 <div>Good afternoon  , Vasilis.</div>
 <div> </div>
 <div>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.</div>
 <div> </div>
 <div>Thank you in advance , Julian. 
 
 </div>
 <div class="gmail_quote">On Wed, May 18, 2011 at 3:32 PM,
 Vasilis Vlachoudis <<a
 moz-do-not-send="true"
 href="mailto:Vasilis.Vlachoudis_at_cern.ch">Vasilis.Vlachoudis_at_cern.ch</a>>
 wrote: 
 <blockquote class="gmail_quote" style="padding-left: 1ex;
 margin: 0px 0px 0px 0.8ex; border-left: 1px solid rgb(204,
 204, 204);">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: <a moz-do-not-send="true"
 href="mailto:owner-fluka-discuss_at_mi.infn.it">owner-fluka-discuss_at_mi.infn.it</a>
 [<a moz-do-not-send="true"
 href="mailto:owner-fluka-discuss_at_mi.infn.it">owner-fluka-discuss_at_mi.infn.it</a>]
 on behalf of Julian Rassolov [<a moz-do-not-send="true"
 href="mailto:julian.rassolov_at_gmail.com">julian.rassolov_at_gmail.com</a>] 
 Sent: 16 May 2011 19:50 
 To: <a moz-do-not-send="true"
 href="mailto:fluka-discuss_at_fluka.org">fluka-discuss_at_fluka.org</a> 
 Subject: (n,gamma) interaction 
 <div>
 <div class="h5"> 
 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. 
 
 
 </div>
 </div>
 </blockquote>
 </div>
 
 </div>
 </blockquote>
 
 </body>
</html>

--------------020305040400090108040208--

--------------010506060005030703000702
Content-Type: application/x-fluka; name="photonuc.inp"
Content-Transfer-Encoding: base64
Content-Disposition: attachment; filename="photonuc.inp"

VElUTEUKUGhvdG8gTnVjbGVhciBzY29yaW5nCkdMT0JBTCAgICAgICAgICAgICAgICAgICAg
ICAgICAgICAgICAgICAgICAgICAgMS4wICAgICAgIDEuMApERUZBVUxUUyAgICAgICAgICAg
ICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgUFJF
Q0lTSU8KQkVBTSAgICAgICAgICAgIC0xLjAgICAgICAgICAgICAgICAgICAgICAgICAgICAg
ICAgICAgICAgICAgICAgICAgICAgIFBIT1RPTgpCRUFNUE9TCkRJU0NBUkQgICAgIEVMRUNU
Uk9OICBQT1NJVFJPTiAgICAgICAgICAgICAgICAgICAgICAgICBQSU9OKyAgICAgUElPTi0K
UEFSVC1USFIgICAgICAgLTAuMDIgICBORVVUUk9OICAgICAgICAgICAgICAgICAgICAgICAg
ICAgMC4wClBBUlQtVEhSICAgICAgIC0wLjAyICAgIFBST1RPTiAgICAgICAgICAgICAgICAg
ICAgICAgICAgIDAuMApQSE9UT05VQyAgICAgICAgICAxLiAgICAgICAgICAgICAgICAgICAg
ICAgICAgTEVBRApHRU9CRUdJTiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAg
ICAgICAgICAgICAgICAgICAgICAgICAgICAgQ09NQk5BTUUKICAgIDAgICAgMCAgICAgICAg
ICAKKiBCbGFjayBib2R5ClNQSCBibGtib2R5ICAgIDAuMCAwLjAgMC4wIDEwMDAwMC4wCiog
Vm9pZCBzcGhlcmUKU1BIIHZvaWQgICAgICAgMC4wIDAuMCAwLjAgMTAwMDAuMAoqIEN5bGlu
ZHJpY2FsIHRhcmdldApSQ0MgdGFyZ2V0ICAgICAwLjAgMC4wIDAuMCAwLjAgMC4wIDEwMDAu
IDEwMC4KRU5ECiogQmxhY2sgaG9sZQpCTEtCT0RZICAgICAgNSArYmxrYm9keSAtdm9pZAoq
IFZvaWQgYXJvdW5kClZPSUQgICAgICAgICA1ICt2b2lkIC10YXJnZXQKKiBUYXJnZXQKVEFS
R0VUICAgICAgIDUgK3RhcmdldApFTkQKR0VPRU5ECiogLi4rLi4uLjEuLi4uKy4uLi4yLi4u
LisuLi4uMy4uLi4rLi4uLjQuLi4uKy4uLi41Li4uLisuLi4uNi4uLi4rLi4uLjcuLgpBU1NJ
R05NQSAgICBCTENLSE9MRSAgIEJMS0JPRFkKQVNTSUdOTUEgICAgICBWQUNVVU0gICAgICBW
T0lECkFTU0lHTk1BICAgICAgICBMRUFEICAgIFRBUkdFVAoqTEFNLUJJQVMgICAgICAgICAw
LjAgICAgICAwLjAyICAgICAgTEVBRCAgICBQSE9UT04KKlVTRVJEVU1QICAgICAgICAxMDAu
ICAgICAgIDUwLiAgICAgICAwLjAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBkdW1w
ClJBTkRPTUlaICAgICAgICAgMS4wClNUQVJUICAgICAgICAgMTAwMDAuClNUT1AK
--------------010506060005030703000702--
Received on Mon Jun 06 2011 - 17:40:45 CEST

This archive was generated by hypermail 2.2.0 : Mon Jun 06 2011 - 17:40:45 CEST