Dear Viktor,
you have to distinguish between primary particles (that you generate
with your source) and secondary particles (that are generated by the
interactions in your geometry). In your example you are seeing 13 gammas
in the output because they are generated by the interactions of some of
your 10 primaries.
An easy way to check this will be, for instance, to replace all your
materials with vacuum: you will get exactly the 10 photons that you
have generated. Or, for instance, you can add the printout of
relevant properties: there are many in the FLUKA common blocks and
you could try LTRACK in (TRACKR). LTRACK is a flag recording the
generation number of the particle and is 1 for the primary and >1
for all the secondaries generated in an event.
JTRACK is the identity number of the particle as it is stored in the
(TRACKR) common block and is shared by the routines that access the
common block, while IJ is a variable in the USRMED routine:
they have the same content, but are not linked directly.
Note that the call to USRMED follows the call to MGDRAW and the print
in your collision file is done even if the kill occurs.
Hope this helps.
Regards,
Andrea
PS
Please reply to the mailing list too: the discussion might be of
interest for others.
Il 11/02/2017 20:40, Віктор Родін ha scritto:
> Dear Andrea,
> Thank you for your help and clear answers. Now i can imagine
> about format of output file. But still I have problems with killing of
> the particles. As you have seen I send 10 gammas. but in output file
> we have 13. How is it possible?I used the proposed methods and this
> don't have effect.
> IJ is alternative of JTRACK?Because if IJ is in the requsition, I
> receive empty files.
> Cheers,
> Viktor
>
> 2017-02-10 23:49 GMT+02:00 Andrea Fontana <andrea.fontana_at_pv.infn.it
> <mailto:andrea.fontana_at_pv.infn.it>>:
>
> Hello again Viktor,
> Ihave double checked the behaviour of the format specifiers
> used in the previous example, '(1P,4G10.4)' for instance, and I write
> with a correction on the meaning of 1P.
>
> 1P is used in connection with the G format specifier: G is used for
> floating point numbers and is an extension of the F format specifier
> to optionally print the result in exponential notation (depending on
> the value). 1P is meant to alter the result of the WRITE statement
> byshifting the location of the decimal point by 1 position,
> transforming for example0.123E10 in 1.234E9.
>
> You can find more details in a Fortran referencebook.
>
> Best,
> Andrea
>
>
> Il 10/02/2017 18:38, Віктор Родін ha scritto:
>
> Dear FLUKA users and experts,
> I'm looking for ways of detecting quantities of particles in
> the defined region and their coordinates of boundary crossing.
> I saw that MGDRAW subroutine gives a lot of usefull
> information about particles
> I searched how the MGDRAW subroutine works and found nice
> answer and MGDRAW file in the discuss :
> http://www.fluka.org/web_archive/earchive/new-fluka-discuss/6806.html
> <http://www.fluka.org/web_archive/earchive/new-fluka-discuss/6806.html>
> <http://www.fluka.org/web_archive/earchive/new-fluka-discuss/6806.html
> <http://www.fluka.org/web_archive/earchive/new-fluka-discuss/6806.html>>
> I'm newbie in FORTRAN77 and ask you to answer me on few
> questions from this topic:
> 1. Everything shifts in output file if I add third coordinate
> in this MGDRAW , How can this be repaired?
> 2. Can I kill each particle which cross boundary? Because
> some particles are registered twice. I read that it could be
> done only with help of the USRMED. Is it true?
> Thank you for attention. INP file and MGDRAW are attached.
> Best regards,
> Viktor
>
>
> --
> ========================================================================
> Dr. Andrea Fontana tel: +39 0382 987991
> <tel:%2B39%200382%20987991>
> Istituto Nazionale fax: +39 0382 423241
> <tel:%2B39%200382%20423241>
> di Fisica Nucleare
> Sezione di Pavia e-mail:
> andrea.fontana_at_pv.infn.it <mailto:andrea.fontana_at_pv.infn.it>
> Via Bassi 6 web :
> www.pv.infn.it/~fontana <http://www.pv.infn.it/%7Efontana>
> 27100 PAVIA, Italy
> ========================================================================
>
>
--
========================================================================
Dr. Andrea Fontana tel: +39 0382 987991
Istituto Nazionale fax: +39 0382 423241
di Fisica Nucleare
Sezione di Pavia e-mail: andrea.fontana_at_pv.infn.it
Via Bassi 6 web : www.pv.infn.it/~fontana
27100 PAVIA, Italy
========================================================================
__________________________________________________________________________
You can manage unsubscription from this mailing list at https://www.fluka.org/fluka.php?id=acc_info
Received on Sun Feb 12 2017 - 00:08:48 CET