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.
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id=acc_info
Received on Mon Nov 16 2015 - 20:54:06 CET