Question about counting pions created in a specific region using MGDRAW

Scoring
1

Versions

  • FLUKA version: 4-5.0

FLUKA:
Flair:

Description

Dear FLUKA experts,

I am a FLUKA user and I have a question regarding the correct way to count the number of positive pions (π⁺) created (produced) within a specific geometry region, not just those passing through it.

My setup:

  • FLUKA version: 4-5.0

  • Input file: Xe_W_3082MeV_0.inp(attached)

  • Geometry: I have a region named TARGET (made of tungsten) where a heavy ion beam (Xe-132) interacts.

  • I want to count the π⁺ produced exclusively inside the TARGET region.

What I have tried so far:

  1. Using USERDUMP -1.0 magdraw TARGET to write only particles created in that region, but I could not do it properly.

My specific questions are:

  1. What is the most straightforward and recommended method to obtain the total number of π⁺ produced in a specific region?

  2. If I use USERDUMP with magdraw, how can I do this? and how can I interpret the output file?

  3. Alternatively, would USRTRACK or USRBDX be more appropriate for this purpose?

I would greatly appreciate any guidance or examples you can provide.

Thank you very much for your time and help.

Best regards,

Jilberto

Input files

Please upload all rele

Xe_W_3082MeV_0.flair (1.6 KB)

Xe_W_3082MeV_0.inp (2.1 KB)

3

Hi Jilberto,

I believe you are on the right track with mgdraw, this routine indeed being the most appropriate method for your problem. You say you also tried to use the mgdraw routine already, could you also provide this file in addition to the .inp and .flair file?

If you get back to me on this I can work on answering your questions.

Best regards,

Andreas

4

Dear Andreas, first of all thank you very much for your answer. I have atached the files.

regards,

Jilberto

Xe_W_3082MeV.inp (2.0 KB)

Xe_W_3082MeV.flair (2.7 KB)

mgdraw.f (3.7 KB)

5

Hello Jilberto,

Sorry for the late reply and thanks for sharing the files. From the looks of it, you are using a highly customised version of the standard mgdraw.f routine. Using only this, I am not able to compile this into an executable to check your results. Could you tell me if any other routines or commons are needed and share those as well?

You also mentioned you used the command “USERDUMP -1.0 magdraw TARGET”, however, in your .inp and /flair files you do not specify the different WHAT entries, it is only “USERDUMP”. For a highly customized version of mgdraw.f like you wrote it is not guaranteed that the calls to the different entries are properly made.

Cheers,

Andreas