Dear Luigi
Hello, following from your help, I was able to come up with the MGDRAW routine.
It compiled okay and I ran the simulations. However, I got following two problems:
1. The files with unit 79 (the unit that I assigned in USDRAW) is not produced after the simulations. I tried to find out why it is not producing.
2. Is there any way to directly get the neutrons undergo capture process? Because my objective is to see how many neutrons undergo capture.
I can see in the list of ICODE that there is no specifically process associated with neutron capture (like nCapture in GEANT4). Hence, I put the ICODE as 101 (inelastic interactions) then look for the their particle ID number to find out.
I enclose the MGDRAW routine and my input file.
Thank you so much for your help.
From: Luigi Salvatore Esposito
Sent: 19 January 2018 12:28
To: David Sangcheol Lee <david.sangcheol.lee_at_cern.ch>
Cc: fluka-discuss <fluka-discuss_at_fluka.org>
Subject: Re: [fluka-discuss]: Possibility to filter the neutrons by the types of interactions in FLUKA
Dear David,
please keep the mail exchange on the list: this way also others can benefit or, even better, contribute.
To activate the entry USDRAW you need to USERDUMP WHAT(4) >= 1, as stated in the manual.
The region number is available as input parameter to the entry USDRAW: it is MREG.
I recommend you to identify the region number the first time the USDRAW is called, like for example with the following code snippet
SUBROUTINE MGDRAW ( ICODE, MREG )
….
SAVE MY_REG
….
ENTRY USDRAW ( ICODE, MREG, XSCO, YSCO, ZSCO )
IF ( .NOT. LFCOPE ) THEN
LFCOPE = .TRUE.
* identify the region
* xxx is your region name (it is CHARACTER*8!)
CALL GEON2R( ‘xxx ', MY_REG, IERR)
IF(IERR.EQ.1) CALL FLABRT ('USDRAW’,'xxx not found')
END IF
then you can check that you are in the wanted region, i.e. MREG == MY_REG
In case you need further help, don’t esitate to write another email and possibly include also all relevant files (FLUKA input and user routines)
Best, luigi
On 19 Jan 2018, at 10:06, David Sangcheol Lee <david.sangcheol.lee_at_cern.ch<mailto:david.sangcheol.lee_at_cern.ch>> wrote:
Thank you so much for your help. I can see that the information can be accessible. I am very new with this kind of source routine stuff so I apology if I ask few more questions.
Firstly, is it possible to limit within the region (e.g. neutrons undergo capture process within one region like target)?
Secondly, is it default parameter for userdump card that I put in fluka code? Or is there specific setup within the userdump card in fluka/flair that I have to do in order to get this code working?
I really appreciated all your help
Thank you
From: Luigi Salvatore Esposito
Sent: 18 January 2018 16:27
To: David Sangcheol Lee <david.sangcheol.lee_at_cern.ch<mailto:david.sangcheol.lee_at_cern.ch>>
Cc: fluka-discuss_at_fluka.org<mailto:fluka-discuss_at_fluka.org>
Subject: Re: [fluka-discuss]: Possibility to filter the neutrons by the types of interactions in FLUKA
Dear David,
you can access this information with the USDRAW entry of MGDRAW.
Please give a look at this post<
http://www.fluka.org/web_archive/earchive/new-fluka-discuss/10527.html> on the list that provides a code snippet that you can adapt to your need.
Best, luigi
On 18 Jan 2018, at 14:57, David Sangcheol Lee <david.sangcheol.lee_at_cern.ch<mailto:david.sangcheol.lee_at_cern.ch>> wrote:
Dear FLUKA expert
Hello, I would like to ask a question. Is it possible to filter the particle by the types of interactions (e.g. inelastic, elastic, capture and fission) in a region?
For example, I would like to see how many neutrons undergo “capture” in a region. Is it possible to do that in FLUKA?
Thank you so much for your help
Kind regards
Cordialement
<image001.jpg> David Sangcheol Lee
CERN / HSE-RP
European Organisation for Nuclear Research
1211 Geneva 23, Switzerland
Email: david.sangcheol.lee_at_cern.ch<mailto:david.sangcheol.lee_at_cern.ch>
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id=acc_info
Received on Thu Feb 01 2018 - 02:30:51 CET