Re: [fluka-discuss]: muon + mgdraw

From: Giuseppe Battistoni <>
Date: Mon, 25 Sep 2017 20:01:33 +0200

Hi Andrea
some partial answers.
I leave to others the question of EMF-BIAS because there are more
important things and I will concentrate on the use of mgdraw.

Such a topic is not easy at all and I'll try to give you some hints.

You are placing the write statements in the USDRAW section (actually
these sections technically are called ENTRY), and you are not using the
proper variables.
In this entry JTRACK is the code of the particle interacting at
coordinates XSCO, YSCO, ZSCO. It has nothing to do with particles
created in that interaction.
Looping on NP you must scroll the arrays of common Genstk
[ see in $FLUPRO/flukapro the content of (GENSTK) ] and you will select
muon by asking for kpart(ip) = 10 or 11. Look carefully at the other
arrays of GENSTK common for all kinematical infos.

JTRACK, ETRACK, XTRACK etc are variables of the common (TRACKR) in which
are copied, one at a time, the variables of transported particles,
taking them from GENSTK. You will use them in other Entries of Mgdraw.
For instance in the BXDRAW entry, those variables will identify the
running particle during transport at the moment in which a boundary
between regions is crossed. USDRAW Entry instead "captures" the moment
in which particles are created, not the action of transport.
Individual particle transport can be followed in the main initial part
of mgdraw.f using the variables of (TRACKR) common block.

I don't know if this is clear enough...


On 25/09/2017 16:47, andrea russomando wrote:
> Dear FLUKA experts
> Iím studying the interaction of an electron beam with a lead target
> I would like to record all the muons created, and save the variables needed to reconstruct their kinematic
> For this purpose Iím using a modified version of mgdraw
> I never used it before, so I have some doubts
> 1) Am I using the right cards?
> I donít know if Iím using in the right way the card EMF-BIAS (the manual suggests to be very careful in the mixing of bias and dump) and I donít know if I should activate the card photonuc (type electrnuc).
> If this card is enable seems to me that the muons are produced at rest, while otherwise have very high energy (I attached the two different dumps)
> 2) I canít understand the results of the dump
> It seems that the same particle is written more than once in the stack
> Also I always obtained Ntrack=0 (from what I understand it means a local energy deposition), even though Iím using a low threshold in the card PART-THR
> Thanks in advance
> andrea

INFN Milano
via Celoria 16, 20133 Milano
tel: +39 02 50317307
fax: +39 02 50317617
You can manage unsubscription from this mailing list at
Received on Mon Sep 25 2017 - 21:09:54 CEST

This archive was generated by hypermail 2.3.0 : Mon Sep 25 2017 - 21:09:56 CEST