Hello FLUKA experts,
I am trying to create a muon source with the source_newgen.f template. These muons are to have a zenith angle distribution and an energy distribution. I have an error, and a few questions:
-
The error:
Region of particle in position -1.136407157E+27 -2.886724341E+26 3.000000000E+03 cannot be determined, error count: 1
I believe the issue is related to my use of thesample_annular_distribution
subroutine. The muon beam I’m trying to define must have particles coming from various positions uniformly distributed on a disk of radius 50 m. With a lower max radius, the simulation seems to run fine. However, above a certain threshold, the annular subroutine could be failing point. This point is so far out of the geometry which otherwise has no error. -
I’d like to use the
sample_histogram_momentum_energy
subroutine. Are there any peculiarities for making a histogram.txt file for this function? Must the probabilities be normalized? -
How can I sample the zenith angles in the same way as I’m doing it for the energy? Can I repurpose the
sample_histogram_momentum_energy
subroutine?
Thank you for your helpful advice.
Regan
nEXO-mu-n.inp (19.1 KB)
source_muon.f (19.3 KB)