Dear Ricardo,
Thanks a lot for your response.
I am trying to record the positions of stopping protons through the whole target(including water channels).
And I set the transport threshold for protons in the PARA-THR at 1e-14 GeV.
Then I use ENDRAW ENTRY in MGDRAW subroutine, the code is
ENTRY ENDRAW ( ICODE, MREG, RULL, XSCO, YSCO, ZSCO )
IF ( .NOT. LFCOPE ) THEN
LFCOPE = .TRUE.
OPEN ( UNIT = 78, FILE = "collfile", STATUS = 'UNKNOWN')
END IF
IF(ICODE .EQ. 12) THEN
IF( JTRACK .EQ. 1 ) WRITE(78,'(1P,4G25.15)')
& XSCO, YSCO, ZSCO
END IF
RETURN
Could you please help me check if this is a correct way of recording the positions of stopping protons ?
I am attaching the input file and mgdraw file.
Thanks in advance.
Best regards,
Qi Ding
________________________________
From: Ricardo Manuel Dos Santos Augusto <r.s.augusto_at_cern.ch>
Sent: 20 May 2021 21:59:45
To: Ding, Qi; fluka-discuss_at_fluka.org
Subject: RE:How to score the number and positions of protons stopping in a region
Dear Qi,
You could score proton fluence using a USRBIN inside the target, and see how the proton beam attenuates. If you want to know about the beam protons only then you could always use BEAMPART.
A more complex approach would be to use mgdraw, since ENDRAW's call from Kaskad has in Icode 12 an option for stopping particle. Note that to capture the whole picture you will need adequate transport cut-offs.
Best regards,
Ricardo
//Ricardo dos Santos Augusto
TRIUMF-EH&S
Office: +1 604 222 10 47 (6375)
4004 Wesbrook Mall,
Vancouver, BC V6T 2A3, Canada//
________________________________________
De : owner-fluka-discuss_at_mi.infn.it [owner-fluka-discuss_at_mi.infn.it] de la part de Ding, Qi [q.ding_at_fz-juelich.de]
Envoyé : mardi 18 mai 2021 9:42
À : fluka-discuss_at_fluka.org
Objet : [fluka-discuss]: How to score the number and positions of protons stopping in a region
Dear FLUKA experts,
In my simulation there is a proton beam(MeV) impinging a metal target, a water channel is inside there. I want to know the number of protons which are stopping inside the target and the positions where protons stopped.
Could you please tell me how to get these information?
Thanks in advance.
Best regards,
Qi
------------------------------------------------------------------------------------------------
------------------------------------------------------------------------------------------------
Forschungszentrum Juelich GmbH
52425 Juelich
Sitz der Gesellschaft: Juelich
Eingetragen im Handelsregister des Amtsgerichts Dueren Nr. HR B 3498
Vorsitzender des Aufsichtsrats: MinDir Volker Rieke
Geschaeftsfuehrung: Prof. Dr.-Ing. Wolfgang Marquardt (Vorsitzender),
Karsten Beneke (stellv. Vorsitzender), Prof. Dr.-Ing. Harald Bolt,
Prof. Dr. Frauke Melchior
------------------------------------------------------------------------------------------------
------------------------------------------------------------------------------------------------
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id=acc_info
Received on Sun Jun 13 2021 - 20:03:32 CEST