Multiple Energies in sample_discrete_momentum_energy File

Scoring
1

Dear FLUKA experts,

I calculated the efficiency of a detector using the DETECT card in combination with the source_newgen.f user routine. Everything works well when I use a single energy value—either directly by setting momentum_energy = 0.00066158 or by using a file (DiscreteEnergy1.dat) that contains only one energy. The results are consistent and correct.

However, when I use a file (DiscreteEnergyAll.dat) containing multiple (13) discrete energies in the sample_discrete_momentum_energy function, the calculated efficiency appears to be divided by the number of energy values used. For example, the resulting efficiency is roughly 1/13 of what it should be.

Could someone please advise how to correctly handle multiple discrete energies in this context? Is there a normalization step I am missing ?

Best regards
DiscreteEnergy1.dat (56 Bytes)
DiscreteEnergyAll.dat (368 Bytes)
source_newgen.f (21.3 KB)

image
image1398×386 41.8 KB

2

Dear @m.assalmi ,

The source routine will normalize the total porbability to 1, so you need to normalize manually the results. In this case multiplying it by 13.

Cheers,
David