Dear FLUKA experts,
I am carrying out detailed simulations of the neutron flux in the Canfranc Underground Laboratory (Canfranc, Spain). Three different neutron sources are considered: (alpha, neutron) reactions, spontaneous fission and (muon, neutron) reactions. The code works as follows:
- For the alpha-induced neutron flux primary neutrons are created within the laboratory walls with energies selected according to an (alpha, neutron) energy spectrum which has been previously calculated using specialized codes. The neutron flux in different regions of the main hall is then computed using the USTRACK card.
- For the spontaneous fission neutron flux an equivalent procedure is used.
- For the muon-induced neutron flux primary muons are created within the laboratory. The (muon, neutron) reactions are therefore computed using FLUKA. Again, the neutron flux in different regions of the main hall is then computed using the USTRACK card.
For the calculations using the neutrons as primary FLUKA particles (alpha-induced reactions and spontaneous fission) the code works fine. However, there is a problem when muons are used: at some point of the simulations the calculation freezes. I have tracked each particle using the mgdraw routine and I have concluded that the freezing happens when muons enter the black hole region. I have no idea of which could be the reason for this behavior.
- The FLUKA version is 4-3.2
- The FLAIR version 3.2-4.4
Relevant files attached:
- LSC_HallA_v3.inp (the input file)
- RunTest.exe (the FLUKA executable)
- Source_LSC.f (the source routine)
- mgdraw_user.f (the mgdraw routine)