However, when I launch my simulation with 2.5 MeV neutrons,I ended with more In116[m] or In113[m] than In115[m] ( with a 25 : 2 and a 3:2 ratio) whereas cross section for theses reactions are normally in favor of In115[m] production (and even more for In113[m] given its weak proportion) .
And when I try to add _HP, there is no metastable production.
Hi ,
I have already tried the Hadr03 and the result is the same,
even in your output file: for 10⁵ incident neutron, In115[336.244]: 3533
Cross section for this reaction is 0.343b
The total cross section is about 1.87b without elastic scattering and 4.96b with it.
So, like all primary particle will interact in material, the final distribution should be according the cross section distribution.
And if we follow this point with 10⁵ primary , we should have for In115m (0.343/1.87)*10⁵=18342 or (0.343/4.96)*10⁵=6915 activated atoms which is about 5 times more than what we can observe.
I agree with your analyse. Below, a mail sent to Alberto and Vladimir :
Neutron (2.5 MeV) on In115. QGSP_BIC constructor selected.
The number of created ions (or excited) do not match the cross sections given by G4HadronicProcessStore.
This does not happen with QGSP_BIC_HP. But there is no metastable states created !
JEFF data for metastable level production may be converted to HP files for In target. You may add these new files to HP data set.
Concerning non-HP model: there are no partial cross sections but transition probabilities for gamma de-excitation and analytical probabilities of other reactions. Here one may check levels and probabilities inside G4LEVELGAMMADATA for In.
But this change also kill the metastable production.
Should I just use QGSP_BIC_HP, retrieving number and energy of incident neutrons on my target, and after the run use tabulated cross section for determining the number of metastable created ?