I have a question regarding energy deposition in a 15-crystal setup. When I simulate energy deposition using, for example, 40 keV photons, I notice that the energy deposited in the later crystals is sometimes higher than in the earlier ones. For instance, here is a sample output:
Ideally, the energy deposition should decrease as we move from the first to the last crystal, but that’s not happening consistently.
I’m generating an energy deposition response matrix from 20 to 200 keV in 10 keV steps and using it to unfold experimental spectra. However, due to these fluctuations, the chi-square function doesn’t perform well.
Is there a known reason for these non-monotonic energy deposition patterns? And is there a recommended method to smooth or correct these fluctuations to improve the response matrix?
Perhaps this is due to statistical fluctuation. There are no errors, only Energy deposition fluctuations. But I would like to avoid it as it is crucial for me. I am taking 1e7 primaries with 5 cycles.
Please find attached the screenshot of the LIS file showing the energy deposition with errors. Since I’m having trouble understanding this, I am also including my inp and flair files for your reference.
From this file, I notice that some of the values you scored have errors of up to 60%. As a general guideline, keep in mind that scored values with errors larger than 50% are completely unreliable. You should aim for errors smaller than 10% ideally and sometimes even smaller for anything serious.
I advise you find a strategy to reduce these errors. This could be run more primaries, use biasing or use a more targeted source.
I tried implementing biasing (see image)
but it didn’t help, could you please see if I applied it correctly. Then I tried to increase primaries to 1e8 which reduces errors significantly compare to 1e7 primaries now the errors are 1e-2 to 32 percent.
But I would like to reduce it more as you mentioned less than 10 percent, what are the other steps should I take?
Regarding the source configuration, I am currently using a point-like, collimated source. If there are any adjustments you would recommend for a more “targeted” source setup, I would appreciate your guidance.
Increasing the primaries by a factor 10 will reduce the error by a factor 1/sqrt(10)=0.3 in the best case which is consistent with what you observe here. To reduce them by a factor 6 you would need 36 times more primaries. See the relevant lecture on Monte Carlo techniques, also in the latest FLUKA course.
As for the source, you may want to consult the relevant page in the manual. In general you just want to make sure that all primaries you simulate actually contribute in a meaningful way to the result you want to score. I do not see an obvious problem in your simulation but simply be aware of that.
