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...
Ciao
Giuseppe
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 https://www.fluka.org/fluka.php?id=acc_info
Received on Mon Sep 25 2017 - 21:09:54 CEST