Dear Alfredo,
Thanks a lot for your reply.
Best regards,
Qi Ding
________________________________
From: alfredo.ferrari_at_mi.infn.it <alfredo.ferrari_at_mi.infn.it>
Sent: 17 June 2021 08:10:17
To: Ding, Qi
Cc: alfredo.ferrari_at_mi.infn.it; Ricardo Manuel Dos Santos Augusto; fluka-discuss_at_fluka.org
Subject: Re: [fluka-discuss]: Re: How to score the number and positions of protons stopping in a region
Yes
> Dear Alfredo,
>
>
> Thanks a lot for your comments.
>
> As I mainly want to record the number and positions of stopping protons,
>
> do I need to make the minimum transport threshold to 1keV or something
> else?
>
>
> Thanks in advance.
>
>
> Best regards,
>
>
> Qi
>
> ________________________________
> From: alfredo.ferrari_at_mi.infn.it <alfredo.ferrari_at_mi.infn.it>
> Sent: 14 June 2021 14:05:14
> To: Ding, Qi
> Cc: Ricardo Manuel Dos Santos Augusto; fluka-discuss_at_fluka.org
> Subject: Re: [fluka-discuss]: Re: How to score the number and positions of
> protons stopping in a region
>
> Dear Qi
>
> your snippet of code is correct, however setting the transport threshold
> of protons to 1E-14 GeV does not make any sense. The minimum transport
> energy for ll (charged) hadrons is 1 keV, and indeed protons at that
> energy do not travel a lot at all. Transport energies down to thermal
> energies and below make sense for neutrons only
>
>> 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 Thu Jun 17 2021 - 09:56:52 CEST