Re: [fluka-discuss]: Re: Primary protons absorbed in a nuclear interaction vs elastic large angle scattered

From: Francesco Cerutti <>
Date: Mon, 16 Nov 2015 21:41:00 +0100

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



Francesco Cerutti
CH-1211 Geneva 23
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: <> em nome de Mikhail Polkovnikov <Михаил Полковников> <>
> Enviado: 5 de novembro de 2015 03:17
> Para:
> 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
Received on Mon Nov 16 2015 - 22:56:51 CET

This archive was generated by hypermail 2.3.0 : Mon Nov 16 2015 - 22:56:57 CET