[fluka-discuss]: Re: How to score the number and positions of protons stopping in a region

From: Ding, Qi <q.ding_at_fz-juelich.de>
Date: Sun, 13 Jun 2021 15:40:56 +0000

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

      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)')
      END IF

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 dos Santos Augusto
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,


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

This archive was generated by hypermail 2.3.0 : Sun Jun 13 2021 - 20:03:40 CEST