DigiMouse Voxel problem

Advanced Features and User Routines
1

Dear FLUKA experts,
I am trying to build simulations using the DigiMouse model, and run into some issues.

  1. When I first created the .vxl file I set all materials to AIR, and then changed that using ASSGNTMAT card. I keep getting a warning that “Multiple ASSIGNMAT for region”. The simulation still runs as intended, and in the output files I see that the material assignment has been set correctly. Can this warning be ignored?
  2. The simulation runs and finishes without errors only if the beam position is not within the voxel area. If the beam position is set to any point within the voxel the run does not seem to finish at all. When I turn on the geometry debugger option, the error is “Too many errors in geometry”, although I cannot find any errors.
  3. When I use a user routine to sample a location within the voxel file, at first there seems that the simulation runs without any problems. However, when I look at the .out files I see the error that the region of the particle cannot be identified.

Please find the attached files.

Thank you kindly in advance,
hen

DigiMouse

DigiMouse3.flair (10.5 KB)
source_newgen_organs.f (19.8 KB)
DigiMouse3.inp (7.6 KB)

3

Dear @rachel.hen.shukrun

  1. The VOXEL file mici_new.vxl has embeded material definitions of the various organs (all to AIR). So flair warns you that there is a duplicate definition in the .inp. In the FLUKA run the last assignmat will be taken. So, yes you can ignore that warning. Probably the best would have been when you create the .vxl to assign the proper materials there
  2. You are hitting a precision problem, since you have set as starting position (BEAMPOS) exactly on the vertex of 8 voxels. If you displace the beam by a tiny bit (few micros) in all x,y,z ensuring your particle doesn’t start on the vertex, then the code will safely find the starting point and it finishes without any problem
  3. I could not reproduce the problem, but it seems you are looping recursively in the regionsource routine, while you should make a simple do while loop
4

Dear @vasilis,
Thank you kindly for your answer.
Regarding the first issue, I have not succeeded to recreate the .vxl with the proper materials, neither using the flair GUI nor directly using the .mat file, since the materials I need are not in the default list. The materials are defined in the .inp file but are not recognized during the .vxl creation, leading to failure.

Thank you kindly for your help,
Hen

5

You can add in the material database the ones that are missing.
From the Materials page you can import from an input or manually add the materials you want. Then you can select them in the Dicom voxel creation page.