Z not conserving in few reaction channels for proton on 7Li reaction

Advanced Features and User Routines
1

Dear experts,

In case of proton on 7Li, while analyzing the mgdraw output for different reactions, I noticed that some reaction channels are not conserving the Z value. Although the number of such reactions are very small, still I want to understand the output.

The output file size is large, hence I am adding those output lines below: I gave a run with 1 cycle with 1E+8 history.

For e.g. Line no. 2081834 : Here Z is not conserved

MGDRAW
Reaction with ICODE 300 in region 3
Primary history: 64079842
JTRACK, LTRACK: 8 2
Generated 3 particles of ID: 8 8 7
Ep 0.15E-02 0.67E-03 0.79E-04
Generated 1 fragments with
A = 6
Z = 3
E = 0.20E-02
Generated a residual with A = 0 , Z = 0

For e.g. Line no. 2035011: Here Z is not conserved

MGDRAW
Reaction with ICODE 300 in region 3
Primary history: 62627997
JTRACK, LTRACK: 8 2
Generated 3 particles of ID: 8 8 -6
Ep 0.71E-03 0.16E-02 0.44E-03
Generated 1 fragments with
A = 2
Z = 1
E = 0.37E-02
Generated a residual with A = 0 , Z = 0

For e.g. Line no. 2003438 Here Z is not conserved

MGDRAW
Reaction with ICODE 300 in region 3
Primary history: 61665215
JTRACK, LTRACK: 8 2
Generated 2 particles of ID: 8 8
Ep 0.36E-02 0.40E-02
Generated 2 fragments with
A = 2 4
Z = 1 2
E = 0.12E-02 0.33E-03
Generated a residual with A = 0 , Z = 0

mgdraw.f (12.4 KB)
p_onLi7_20MeV.inp (2.2 KB)

Regards,

Riya

2

As a general rule, charge and baryon numbers are obviously conserved in any reaction.
Specifically, the three cases you mention refer to n+7Li reactions (and not p+7Li), as the JTRACK value you report indicates. Therefore, 3 protons are always found in the final state as expected (as part of 6Li or alpha plus deuteron).

1 Like
3

Thank you @ceruttif for the clarification. I totally overlooked the JTRACK values.

Regards,

Riya