On 16.11.2015 21:07, Ferreira De Almeida Lourenco, Ana wrote:
> Dear Mikhail and All,
>
> Many thanks for your e-mail and for the example of the mgdraw.f subrotine!
>
> Regarding the mgdraw subroutine I still have some questions:
> 1- When using Icode = 101, does it include non-elastic nuclear interactions as well? If not, what is the Icode that I should use for non-elastic nuclear interactions?
> 2- What are the units of ETRACK? I'm simulating a 200 MeV proton beam on a carbon target, and scoring the energy from inelastic interactions. The file "carbon_200001_fort.98" is an example of my output file and the energy of the particles that undergo an inelastic interaction is about 1 GeV which does not make sense since my projectile has a maximum energy of 200 MeV.
>
> Many thanks,
> Ana
>
> ________________________________________
> De: owner-fluka-discuss_at_mi.infn.it <owner-fluka-discuss_at_mi.infn.it> em nome de Mikhail Polkovnikov <Михаил Полковников> <pmk_at_ihep.ru>
> Enviado: 5 de novembro de 2015 03:17
> Para: fluka-discuss_at_fluka.org
> Assunto: [fluka-discuss]: Re: Primary protons absorbed in a nuclear interaction vs elastic large angle scattered
>
> Dear Ana,
>> Unfortunately, I am not familiar with the use of subrotines. From my understanding I need to use USERDUMP and link fluka with a modified version of mgdraw.f, where for my application I need to use the ENTRY USDRAW. Is this correct?
> Correct.
>
> USERDUMP should look like:
>
> *...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....
> USERDUMP 100. 98. 2. 1. dump
>
>
>> Also, how can I have information regarding the residual energy and angle of the protons that undergo an elastic nuclear interaction?
> Look at the (TRACKR) header file, it stores current track parameters as:
> energy, directions etc.
>
> In attachment USDRAW writes interaction events in the TARGET region.
>
> Best regards.
>
>
>
Dear Ana,
2- Units of ETRACK is GeV, but it's a TOTAL energy of the particle which
is a sum of kinetic plus the proton rest mass (938 MeV).
Best regards,
Mikhail.
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id=acc_info
Received on Mon Nov 16 2015 - 21:12:01 CET