Hi Ana
complementing Mikhail's answer, which already hit the main point, the
proton mass is more precisely 938.27231 MeV, as given by the AMPRTN
parameter in $FLUPRO/flukapro/\(DBLPRC\).
Concerning ICODE = 101 in mgdraw, it refers to any non-elastic nuclear
reaction, including what you previously indicated as inelastic, which in
that meaning is just a specific subset of non-elastic channels. Note
however that often in real life inelastic is used (as in the USDRAW
header) as synonym of non-elastic, and not in the restricted meaning you
reported
Cheers
Francesco
**************************************************
Francesco Cerutti
CERN-EN/STI
CH-1211 Geneva 23
Switzerland
tel. ++41 22 7678962
fax ++41 22 7668854
On Mon, 16 Nov 2015, 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.
>
>
>
>
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id=acc_info
Received on Mon Nov 16 2015 - 22:56:51 CET