Dear Andrea,
As I understood, 13 gammas are sum of the primary and secondary
particles. And if I used requisition (LTRACK .EQ. 1) then I will see my
primary particles in the output file.
Regards,
Viktor
2017-02-12 0:00 GMT+02:00 Andrea Fontana <andrea.fontana_at_pv.infn.it>:
> 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 - 13:05:35 CET