Test of HINCLCascadeIntranuke and HG4BertCascIntranuke models.
Opened this issue · 7 comments
idkakorin commented
I begin to test the HINCLCascadeIntranuke model. The following event seems strange to me, because in it a proton is mother of proton, neutron, deuterium and helium.
|------------------------------------------------------------------------------------------------------------------|
|GENIE GHEP Event Record [print level: 3] |
|------------------------------------------------------------------------------------------------------------------|
| Idx | Name | Ist | PDG | Mother | Daughter | Px | Py | Pz | E | m |
|------------------------------------------------------------------------------------------------------------------|
| 0 | e- | 0 | 11 | -1 | -1 | 4 | 4 | 0.000 | 0.000 | 3.000 | 3.000 | 0.001 |
| 1 | C12 | 0 | 1000060120 | -1 | -1 | 2 | 3 | 0.000 | 0.000 | 0.000 | 11.175 | 11.175 |
| 2 | proton | 11 | 2212 | 1 | -1 | 5 | 5 | -0.064 | 0.112 | -0.102 | 0.921 | **0.938 | M = 0.906
| 3 | B11 | 2 | 1000050110 | 1 | -1 | 10 | 10 | 0.064 | -0.112 | 0.102 | 10.254 | 10.253 |
| 4 | e- | 1 | 11 | 0 | -1 | -1 | -1 | -0.172 | -0.185 | 2.920 | 2.931 | 0.001 | P = (0.059,0.063,-0.996)
| 5 | proton | 14 | 2212 | 2 | -1 | 6 | 9 | 0.108 | 0.296 | -0.022 | 0.990 | 0.938 |
| 6 | H2 | 1 | 1000010020 | 5 | -1 | -1 | -1 | 0.162 | 0.216 | -0.121 | 1.899 | 1.876 |
| 7 | proton | 1 | 2212 | 5 | -1 | -1 | -1 | -0.015 | -0.064 | -0.067 | 0.943 | 0.938 |
| 8 | neutron | 1 | 2112 | 5 | -1 | -1 | -1 | 0.017 | 0.031 | 0.057 | 0.942 | 0.940 |
| 9 | He4 | 1 | 1000020040 | 5 | -1 | -1 | -1 | 0.051 | -0.082 | 0.120 | 3.731 | 3.727 |
| 10 | He4 | 1 | 1000020040 | 3 | -1 | -1 | -1 | -0.044 | 0.083 | 0.090 | 3.730 | 3.727 |
|------------------------------------------------------------------------------------------------------------------|
| Fin-Init: | 0.000 | -0.000 | -0.000 | -0.000 | |
|------------------------------------------------------------------------------------------------------------------|
| Vertex: e- @ (x = 0.00000 m, y = 0.00000 m, z = 0.00000 m, t = 0.000000e+00 s) |
|------------------------------------------------------------------------------------------------------------------|
| Err flag [bits:15->0] : 0000000000000000 | 1st set: none |
| Err mask [bits:15->0] : 1111111111111111 | Is unphysical: NO | Accepted: YES |
|------------------------------------------------------------------------------------------------------------------|
| sig(Ev) = 3.37581e-29 cm^2 | dsig(Q2;E)/dQ2 = 4.52020e-29 cm^2/GeV^2 | Weight = 1.00000 |
|------------------------------------------------------------------------------------------------------------------|
--------------------------------------------------------------------------------------------------------------
GENIE Interaction Summary
--------------------------------------------------------------------------------------------------------------
[-] [Init-State]
|--> probe : PDG-code = 11 (e-)
|--> nucl. target : Z = 6, A = 12, PDG-Code = 1000060120 (C12)
|--> hit nucleon : PDC-Code = 2212 (proton)
|--> hit quark : no set
|--> probe 4P : (E = 3.000000, Px = 0.000000, Py = 0.000000, Pz = 3.000000)
|--> target 4P : (E = 11.174863, Px = 0.000000, Py = 0.000000, Pz = 0.000000)
|--> nucleon 4P : (E = 0.921003, Px = -0.063703, Py = 0.111874, Pz = -0.101640)
[-] [Process-Info]
|--> Interaction : EM
|--> Scattering : QES
[-] [Kinematics]
|--> *Selected* Bjorken x = 0.525450
|--> *Selected* Inelasticity y = 0.020267
|--> *Selected* Momentum transfer Q2 (>0) = 0.065342
|--> *Selected* Hadronic invariant mass W = 0.938272
[-] [Exclusive Process Info]
|--> charm prod. : false |--> strange prod. : false
|--> f/s nucleons : N(p) = 0 N(n) = 0
|--> f/s pions : N(pi^0) = 0 N(pi^+) = 0 N(pi^-) = 0
|--> f/s Other : N(gamma) = 0 N(Rho^0) = 0 N(Rho^+) = 0 N(Rho^-) = 0
|--> resonance : [not set]
|--> final quark prod. : false
|--> final lepton prod. : false
--------------------------------------------------------------------------------------------------------------
The above event has number 2 and were generated with command
gevgen -n 200 -p 11 -e 3 -t 1000060120 --tune G18_10a_00_000 --cross-sections xsec_splines_G18_10a_00_000.xml --event-generator-list EMQE --seed 65539
I set <param type="alg" name="HadronTransp-Model"> genie::HINCLCascadeIntranuke/Default </param> in the configuration file G18_10a/ModelConfiguration.xml.
Another example of such event has number 5:
|------------------------------------------------------------------------------------------------------------------|
|GENIE GHEP Event Record [print level: 3] |
|------------------------------------------------------------------------------------------------------------------|
| Idx | Name | Ist | PDG | Mother | Daughter | Px | Py | Pz | E | m |
|------------------------------------------------------------------------------------------------------------------|
| 0 | e- | 0 | 11 | -1 | -1 | 4 | 4 | 0.000 | 0.000 | 3.000 | 3.000 | 0.001 |
| 1 | C12 | 0 | 1000060120 | -1 | -1 | 2 | 3 | 0.000 | 0.000 | 0.000 | 11.175 | 11.175 |
| 2 | proton | 11 | 2212 | 1 | -1 | 5 | 5 | -0.020 | -0.029 | 0.075 | 0.922 | **0.938 | M = 0.918
| 3 | B11 | 2 | 1000050110 | 1 | -1 | 11 | 11 | 0.020 | 0.029 | -0.075 | 10.253 | 10.253 |
| 4 | e- | 1 | 11 | 0 | -1 | -1 | -1 | 0.367 | -0.098 | 2.856 | 2.881 | 0.001 | P = (-0.128,0.034,-0.991)
| 5 | proton | 14 | 2212 | 2 | -1 | 6 | 10 | -0.388 | 0.069 | 0.219 | 1.041 | 0.938 |
| 6 | neutron | 1 | 2112 | 5 | -1 | -1 | -1 | 0.079 | 0.153 | 0.120 | 0.963 | 0.940 |
| 7 | proton | 1 | 2212 | 5 | -1 | -1 | -1 | -0.057 | -0.084 | -0.098 | 0.949 | 0.938 |
| 8 | proton | 1 | 2212 | 5 | -1 | -1 | -1 | -0.066 | 0.008 | -0.019 | 0.941 | 0.938 |
| 9 | H2 | 1 | 1000010020 | 5 | -1 | -1 | -1 | -0.019 | 0.201 | 0.076 | 1.888 | 1.876 |
| 10 | H3 | 1 | 1000010030 | 5 | -1 | -1 | -1 | -0.220 | -0.075 | 0.116 | 2.821 | 2.809 |
| 11 | He4 | 1 | 1000020040 | 3 | -1 | -1 | -1 | -0.084 | -0.105 | -0.051 | 3.730 | 3.727 |
|------------------------------------------------------------------------------------------------------------------|
| Fin-Init: | -0.000 | 0.000 | -0.000 | -0.002 | |
|------------------------------------------------------------------------------------------------------------------|
| Vertex: e- @ (x = 0.00000 m, y = 0.00000 m, z = 0.00000 m, t = 0.000000e+00 s) |
|------------------------------------------------------------------------------------------------------------------|
| Err flag [bits:15->0] : 0000000000000000 | 1st set: none |
| Err mask [bits:15->0] : 1111111111111111 | Is unphysical: NO | Accepted: YES |
|------------------------------------------------------------------------------------------------------------------|
| sig(Ev) = 3.37581e-29 cm^2 | dsig(Q2;E)/dQ2 = 6.09729e-30 cm^2/GeV^2 | Weight = 1.00000 |
|------------------------------------------------------------------------------------------------------------------|
--------------------------------------------------------------------------------------------------------------
GENIE Interaction Summary
--------------------------------------------------------------------------------------------------------------
[-] [Init-State]
|--> probe : PDG-code = 11 (e-)
|--> nucl. target : Z = 6, A = 12, PDG-Code = 1000060120 (C12)
|--> hit nucleon : PDC-Code = 2212 (proton)
|--> hit quark : no set
|--> probe 4P : (E = 3.000000, Px = 0.000000, Py = 0.000000, Pz = 3.000000)
|--> target 4P : (E = 11.174863, Px = 0.000000, Py = 0.000000, Pz = 0.000000)
|--> nucleon 4P : (E = 0.921982, Px = -0.020381, Py = -0.028634, Pz = 0.074807)
[-] [Process-Info]
|--> Interaction : EM
|--> Scattering : QES
[-] [Kinematics]
|--> *Selected* Bjorken x = 0.802749
|--> *Selected* Inelasticity y = 0.037037
|--> *Selected* Momentum transfer Q2 (>0) = 0.151128
|--> *Selected* Hadronic invariant mass W = 0.938272
[-] [Exclusive Process Info]
|--> charm prod. : false |--> strange prod. : false
|--> f/s nucleons : N(p) = 0 N(n) = 0
|--> f/s pions : N(pi^0) = 0 N(pi^+) = 0 N(pi^-) = 0
|--> f/s Other : N(gamma) = 0 N(Rho^0) = 0 N(Rho^+) = 0 N(Rho^-) = 0
|--> resonance : [not set]
|--> final quark prod. : false
|--> final lepton prod. : false
--------------------------------------------------------------------------------------------------------------
Originally posted by @idkakorin in #355 (comment)
idkakorin commented
This event is generated with command gevgen -n 200 -p 11 -e 3 -t 1000060120 --tune G18_10a_00_000 --cross-sections xsec_splines_G18_10a_00_000.xml --event-generator-list EMRES --seed 65539 look strange too:
|------------------------------------------------------------------------------------------------------------------|
|GENIE GHEP Event Record [print level: 3] |
|------------------------------------------------------------------------------------------------------------------|
| Idx | Name | Ist | PDG | Mother | Daughter | Px | Py | Pz | E | m |
|------------------------------------------------------------------------------------------------------------------|
| 0 | e- | 0 | 11 | -1 | -1 | 4 | 4 | 0.000 | 0.000 | 3.000 | 3.000 | 0.001 |
| 1 | C12 | 0 | 1000060120 | -1 | -1 | 2 | 3 | 0.000 | 0.000 | 0.000 | 11.175 | 11.175 |
| 2 | proton | 11 | 2212 | 1 | -1 | 5 | 5 | 0.020 | -0.144 | -0.100 | 0.921 | **0.938 | M = 0.904
| 3 | B11 | 2 | 1000050110 | 1 | -1 | 11 | 12 | -0.020 | 0.144 | 0.100 | 10.254 | 10.253 |
| 4 | e- | 1 | 11 | 0 | -1 | -1 | -1 | 0.287 | -0.187 | 1.956 | 1.986 | 0.001 | P = (-0.145,0.094,-0.985)
| 5 | N+(1720) | 3 | 202214 | 2 | -1 | 6 | 7 | -0.267 | 0.044 | 0.943 | 1.935 | **1.720 | M = 1.667
| 6 | proton | 14 | 2212 | 5 | -1 | 8 | 8 | -0.716 | 0.178 | 0.748 | 1.409 | 0.938 | FSI = 1
| 7 | pi0 | 14 | 111 | 5 | -1 | 9 | 10 | 0.449 | -0.135 | 0.195 | 0.526 | 0.135 |
| 8 | proton | 1 | 2212 | 6 | -1 | -1 | -1 | -0.716 | 0.178 | 0.748 | 1.409 | 0.938 |
| 9 | pi0 | 1 | 111 | 7 | -1 | -1 | -1 | 0.171 | 0.002 | 0.403 | 0.458 | 0.135 |
| 10 | neutron | 1 | 2112 | 7 | -1 | -1 | -1 | 0.213 | -0.196 | -0.107 | 0.989 | 0.940 |
| 11 | He4 | 1 | 1000020040 | 3 | -1 | -1 | -1 | -0.025 | 0.109 | -0.070 | 3.730 | 3.727 |
| 12 | Li6 | 1 | 1000030060 | 3 | -1 | -1 | -1 | 0.070 | 0.094 | 0.070 | 5.603 | 5.602 |
|------------------------------------------------------------------------------------------------------------------|
| Fin-Init: | 0.000 | -0.000 | 0.000 | -0.000 | |
|------------------------------------------------------------------------------------------------------------------|
| Vertex: e- @ (x = 0.00000 m, y = 0.00000 m, z = 0.00000 m, t = 0.000000e+00 s) |
|------------------------------------------------------------------------------------------------------------------|
| Err flag [bits:15->0] : 0000000000000000 | 1st set: none |
| Err mask [bits:15->0] : 1111111111111111 | Is unphysical: NO | Accepted: YES |
|------------------------------------------------------------------------------------------------------------------|
| sig(Ev) = 1.72623e-31 cm^2 | d2sig(W,Q2;E)/dWdQ2 = 2.39707e-30 cm^2/GeV^3 | Weight = 1.00000 |
|------------------------------------------------------------------------------------------------------------------|
--------------------------------------------------------------------------------------------------------------
GENIE Interaction Summary
--------------------------------------------------------------------------------------------------------------
[-] [Init-State]
|--> probe : PDG-code = 11 (e-)
|--> nucl. target : Z = 6, A = 12, PDG-Code = 1000060120 (C12)
|--> hit nucleon : PDC-Code = 2212 (proton)
|--> hit quark : no set
|--> probe 4P : (E = 3.000000, Px = 0.000000, Py = 0.000000, Pz = 3.000000)
|--> target 4P : (E = 11.174863, Px = 0.000000, Py = 0.000000, Pz = 0.000000)
|--> nucleon 4P : (E = 0.920798, Px = 0.020237, Py = -0.143581, Pz = -0.100412)
[-] [Process-Info]
|--> Interaction : EM
|--> Scattering : RES
[-] [Kinematics]
|--> *Selected* Bjorken x = 0.083466
|--> *Selected* Inelasticity y = 0.349569
|--> *Selected* Momentum transfer Q2 (>0) = 0.178776
|--> *Selected* Hadronic invariant mass W = 1.667299
[-] [Exclusive Process Info]
|--> charm prod. : false |--> strange prod. : false
|--> f/s nucleons : N(p) = 0 N(n) = 0
|--> f/s pions : N(pi^0) = 0 N(pi^+) = 0 N(pi^-) = 0
|--> f/s Other : N(gamma) = 0 N(Rho^0) = 0 N(Rho^+) = 0 N(Rho^-) = 0
|--> resonance : P13(1720)
|--> final quark prod. : false
|--> final lepton prod. : false
--------------------------------------------------------------------------------------------------------------
My opinion is that it should look like:
|------------------------------------------------------------------------------------------------------------------|
|GENIE GHEP Event Record [print level: 3] |
|------------------------------------------------------------------------------------------------------------------|
| Idx | Name | Ist | PDG | Mother | Daughter | Px | Py | Pz | E | m |
|------------------------------------------------------------------------------------------------------------------|
| 0 | e- | 0 | 11 | -1 | -1 | 4 | 4 | 0.000 | 0.000 | 3.000 | 3.000 | 0.001 |
| 1 | C12 | 0 | 1000060120 | -1 | -1 | 2 | 3 | 0.000 | 0.000 | 0.000 | 11.175 | 11.175 |
| 2 | proton | 11 | 2212 | 1 | -1 | 5 | 5 | 0.020 | -0.144 | -0.100 | 0.921 | **0.938 | M = 0.904
| 3 | B11 | 2 | 1000050110 | 1 | -1 | 11 | 12 | -0.020 | 0.144 | 0.100 | 10.254 | 10.253 |
| 4 | e- | 1 | 11 | 0 | -1 | -1 | -1 | 0.287 | -0.187 | 1.956 | 1.986 | 0.001 | P = (-0.145,0.094,-0.985)
| 5 | N+(1720) | 3 | 202214 | 2 | -1 | 6 | 7 | -0.267 | 0.044 | 0.943 | 1.935 | **1.720 | M = 1.667
| 6 | proton | 14 | 2212 | 5 | -1 | 8 | 10 | -0.716 | 0.178 | 0.748 | 1.409 | 0.938 | FSI = 1
| 7 | pi0 | 14 | 111 | 5 | -1 | 8 | 10 | 0.449 | -0.135 | 0.195 | 0.526 | 0.135 |
| 8 | proton | 1 | 2212 | 6 | 7 | -1 | -1 | -0.716 | 0.178 | 0.748 | 1.409 | 0.938 |
| 9 | pi0 | 1 | 111 | 6 | 7 | -1 | -1 | 0.171 | 0.002 | 0.403 | 0.458 | 0.135 |
| 10 | neutron | 1 | 2112 | 6 | 7 | -1 | -1 | 0.213 | -0.196 | -0.107 | 0.989 | 0.940 |
| 11 | He4 | 1 | 1000020040 | 3 | -1 | -1 | -1 | -0.025 | 0.109 | -0.070 | 3.730 | 3.727 |
| 12 | Li6 | 1 | 1000030060 | 3 | -1 | -1 | -1 | 0.070 | 0.094 | 0.070 | 5.603 | 5.602 |
|------------------------------------------------------------------------------------------------------------------|
| Fin-Init: | 0.000 | -0.000 | 0.000 | -0.000 | |
|------------------------------------------------------------------------------------------------------------------|
| Vertex: e- @ (x = 0.00000 m, y = 0.00000 m, z = 0.00000 m, t = 0.000000e+00 s) |
|------------------------------------------------------------------------------------------------------------------|
| Err flag [bits:15->0] : 0000000000000000 | 1st set: none |
| Err mask [bits:15->0] : 1111111111111111 | Is unphysical: NO | Accepted: YES |
|------------------------------------------------------------------------------------------------------------------|
| sig(Ev) = 1.72623e-31 cm^2 | d2sig(W,Q2;E)/dWdQ2 = 2.39707e-30 cm^2/GeV^3 | Weight = 1.00000 |
|------------------------------------------------------------------------------------------------------------------|
--------------------------------------------------------------------------------------------------------------
GENIE Interaction Summary
--------------------------------------------------------------------------------------------------------------
[-] [Init-State]
|--> probe : PDG-code = 11 (e-)
|--> nucl. target : Z = 6, A = 12, PDG-Code = 1000060120 (C12)
|--> hit nucleon : PDC-Code = 2212 (proton)
|--> hit quark : no set
|--> probe 4P : (E = 3.000000, Px = 0.000000, Py = 0.000000, Pz = 3.000000)
|--> target 4P : (E = 11.174863, Px = 0.000000, Py = 0.000000, Pz = 0.000000)
|--> nucleon 4P : (E = 0.920798, Px = 0.020237, Py = -0.143581, Pz = -0.100412)
[-] [Process-Info]
|--> Interaction : EM
|--> Scattering : RES
[-] [Kinematics]
|--> *Selected* Bjorken x = 0.083466
|--> *Selected* Inelasticity y = 0.349569
|--> *Selected* Momentum transfer Q2 (>0) = 0.178776
|--> *Selected* Hadronic invariant mass W = 1.667299
[-] [Exclusive Process Info]
|--> charm prod. : false |--> strange prod. : false
|--> f/s nucleons : N(p) = 0 N(n) = 0
|--> f/s pions : N(pi^0) = 0 N(pi^+) = 0 N(pi^-) = 0
|--> f/s Other : N(gamma) = 0 N(Rho^0) = 0 N(Rho^+) = 0 N(Rho^-) = 0
|--> resonance : P13(1720)
|--> final quark prod. : false
|--> final lepton prod. : false
--------------------------------------------------------------------------------------------------------------
similar to this MEC event
1699479032 NOTICE gevdump : [n] <gEvDump.cxx::main (177)> : ** Event: 1
|------------------------------------------------------------------------------------------------------------------|
|GENIE GHEP Event Record [print level: 3] |
|------------------------------------------------------------------------------------------------------------------|
| Idx | Name | Ist | PDG | Mother | Daughter | Px | Py | Pz | E | m |
|------------------------------------------------------------------------------------------------------------------|
| 0 | e- | 0 | 11 | -1 | -1 | 4 | 4 | 0.000 | 0.000 | 3.000 | 3.000 | 0.001 |
| 1 | C12 | 0 | 1000060120 | -1 | -1 | 2 | 3 | 0.000 | 0.000 | 0.000 | 11.175 | 11.175 |
| 2 | n+n | 11 | 2000000200 | 1 | -1 | 5 | 5 | 0.194 | -0.056 | -0.268 | 1.909 | 1.879 |
| 3 | C10 | 2 | 1000060100 | 1 | -1 | 14 | 15 | -0.194 | 0.056 | 0.268 | 9.266 | **9.328 | M = 9.260
| 4 | e- | 1 | 11 | 0 | -1 | -1 | -1 | 0.137 | -0.344 | 2.638 | 2.664 | 0.001 | P = (-0.052,0.129,-0.990)
| 5 | n+n | 3 | 2000000200 | 2 | -1 | 6 | 7 | 0.057 | 0.288 | 0.094 | 2.245 | **1.879 | M = 2.224
| 6 | neutron | 14 | 2112 | 5 | -1 | 8 | 13 | -0.180 | -0.415 | -0.017 | 1.043 | 0.940 |
| 7 | neutron | 14 | 2112 | 5 | -1 | 8 | 13 | 0.237 | 0.703 | 0.111 | 1.202 | 0.940 |
| 8 | neutron | 1 | 2112 | 6 | 7 | -1 | -1 | -0.010 | -0.251 | -0.184 | 0.990 | 0.940 |
| 9 | neutron | 1 | 2112 | 6 | 7 | -1 | -1 | 0.073 | 0.561 | 0.387 | 1.163 | 0.940 |
| 10 | neutron | 1 | 2112 | 6 | 7 | -1 | -1 | 0.123 | 0.070 | -0.050 | 0.951 | 0.940 |
| 11 | proton | 1 | 2212 | 6 | 7 | -1 | -1 | -0.070 | -0.035 | 0.168 | 0.956 | 0.938 |
| 12 | proton | 1 | 2212 | 6 | 7 | -1 | -1 | -0.033 | 0.062 | 0.025 | 0.941 | 0.938 |
| 13 | proton | 1 | 2212 | 6 | 7 | -1 | -1 | 0.004 | 0.045 | -0.064 | 0.942 | 0.938 |
| 14 | H2 | 1 | 1000010020 | 3 | -1 | -1 | -1 | -0.175 | -0.125 | 0.129 | 1.892 | 1.876 |
| 15 | He4 | 1 | 1000020040 | 3 | -1 | -1 | -1 | -0.048 | 0.017 | -0.047 | 3.728 | 3.727 |
|------------------------------------------------------------------------------------------------------------------|
| Fin-Init: | -0.000 | -0.000 | -0.000 | 0.052 | |
|------------------------------------------------------------------------------------------------------------------|
| Vertex: e- @ (x = 0.00000 m, y = 0.00000 m, z = 0.00000 m, t = 0.000000e+00 s) |
|------------------------------------------------------------------------------------------------------------------|
| Err flag [bits:15->0] : 0000000000000000 | 1st set: none |
| Err mask [bits:15->0] : 1111111111111111 | Is unphysical: NO | Accepted: YES |
|------------------------------------------------------------------------------------------------------------------|
| sig(Ev) = 8.92144e-32 cm^2 | dsig(Ev;{K_s})/dK = 0.00000e+00 cm^2/{K} | Weight = 1.00000 |
|------------------------------------------------------------------------------------------------------------------|
--------------------------------------------------------------------------------------------------------------
GENIE Interaction Summary
--------------------------------------------------------------------------------------------------------------
[-] [Init-State]
|--> probe : PDG-code = 11 (e-)
|--> nucl. target : Z = 6, A = 12, PDG-Code = 1000060120 (C12)
|--> hit nucleon : no set
|--> hit quark : no set
|--> probe 4P : (E = 3.000000, Px = 0.000000, Py = 0.000000, Pz = 3.000000)
|--> target 4P : (E = 11.174863, Px = 0.000000, Py = 0.000000, Pz = 0.000000)
[-] [Process-Info]
|--> Interaction : EM
|--> Scattering : MEC
[-] [Kinematics]
|--> *Selected* Bjorken x = 0.098832
|--> *Selected* Inelasticity y = 0.120152
|--> *Selected* Momentum transfer Q2 (>0) = 0.155110
|--> *Selected* Hadronic invariant mass W = 2.223811
[-] [Exclusive Process Info]
|--> charm prod. : false |--> strange prod. : false
|--> f/s nucleons : N(p) = 0 N(n) = 0
|--> f/s pions : N(pi^0) = 0 N(pi^+) = 0 N(pi^-) = 0
|--> f/s Other : N(gamma) = 0 N(Rho^0) = 0 N(Rho^+) = 0 N(Rho^-) = 0
|--> resonance : [not set]
|--> final quark prod. : false
|--> final lepton prod. : false
--------------------------------------------------------------------------------------------------------------
idkakorin commented
Another strange event 2 (gevgen -n 200 -p 11 -e 3 -t 1000060120 --tune G18_10a_00_000 --cross-sections xsec_splines_G18_10a_00_000.xml --event-generator-list EMDIS --seed 65539):
|------------------------------------------------------------------------------------------------------------------|
|GENIE GHEP Event Record [print level: 3] |
|------------------------------------------------------------------------------------------------------------------|
| Idx | Name | Ist | PDG | Mother | Daughter | Px | Py | Pz | E | m |
|------------------------------------------------------------------------------------------------------------------|
| 0 | e- | 0 | 11 | -1 | -1 | 4 | 4 | 0.000 | 0.000 | 3.000 | 3.000 | 0.001 |
| 1 | C12 | 0 | 1000060120 | -1 | -1 | 2 | 3 | 0.000 | 0.000 | 0.000 | 11.175 | 11.175 |
| 2 | neutron | 11 | 2112 | 1 | -1 | 5 | 5 | 0.119 | -0.001 | 0.038 | 0.920 | **0.940 | M = 0.911
| 3 | C11 | 2 | 1000060110 | 1 | -1 | 21 | 22 | -0.119 | 0.001 | -0.038 | 10.255 | 10.254 |
| 4 | e- | 1 | 11 | 0 | -1 | -1 | -1 | 0.075 | 0.110 | 1.059 | 1.067 | 0.001 | P = (-0.070,-0.103,-0.992)
| 5 | HadrSyst | 12 | 2000000001 | 2 | -1 | 6 | 11 | 0.045 | -0.111 | 1.979 | 2.853 | **0.000 | M = 2.051
| 6 | neutron | 14 | 2112 | 5 | -1 | 12 | 12 | -0.220 | -0.044 | 0.605 | 1.140 | 0.940 | FSI = 1
| 7 | pi0 | 14 | 111 | 5 | -1 | 16 | 20 | 0.137 | 0.022 | 0.190 | 0.271 | 0.135 |
| 8 | pi0 | 14 | 111 | 5 | -1 | 13 | 13 | 0.051 | -0.005 | 0.311 | 0.343 | 0.135 | FSI = 1
| 9 | pi0 | 14 | 111 | 5 | -1 | 14 | 14 | 0.053 | 0.099 | 0.400 | 0.437 | 0.135 | FSI = 1
| 10 | pi+ | 14 | 211 | 5 | -1 | 16 | 20 | -0.057 | 0.040 | 0.063 | 0.168 | 0.140 |
| 11 | pi- | 14 | -211 | 5 | -1 | 15 | 15 | 0.080 | -0.223 | 0.409 | 0.493 | 0.140 | FSI = 1
| 12 | neutron | 1 | 2112 | 6 | -1 | -1 | -1 | -0.220 | -0.044 | 0.605 | 1.140 | 0.940 |
| 13 | pi0 | 1 | 111 | 8 | -1 | -1 | -1 | 0.051 | -0.005 | 0.311 | 0.343 | 0.135 |
| 14 | pi0 | 1 | 111 | 9 | -1 | -1 | -1 | 0.053 | 0.099 | 0.400 | 0.437 | 0.135 |
| 15 | pi- | 1 | -211 | 11 | -1 | -1 | -1 | 0.080 | -0.223 | 0.409 | 0.493 | 0.140 |
| 16 | pi- | 1 | -211 | 7 | 10 | -1 | -1 | 0.026 | 0.151 | 0.074 | 0.220 | 0.140 |
| 17 | proton | 1 | 2212 | 7 | 10 | -1 | -1 | 0.086 | 0.077 | 0.138 | 0.955 | 0.938 |
| 18 | proton | 1 | 2212 | 7 | 10 | -1 | -1 | 0.028 | -0.009 | -0.294 | 0.984 | 0.938 |
| 19 | proton | 1 | 2212 | 7 | 10 | -1 | -1 | -0.107 | -0.012 | 0.029 | 0.945 | 0.938 |
| 20 | proton | 1 | 2212 | 7 | 10 | -1 | -1 | 0.003 | -0.020 | -0.119 | 0.946 | 0.938 |
| 21 | He3 | 1 | 1000020030 | 3 | -1 | -1 | -1 | 0.002 | -0.253 | 0.713 | 2.909 | 2.808 |
| 22 | He4 | 1 | 1000020040 | 3 | -1 | -1 | -1 | -0.078 | 0.129 | -0.326 | 3.745 | 3.727 |
|------------------------------------------------------------------------------------------------------------------|
| Fin-Init: | 0.000 | -0.000 | 0.000 | 0.009 | |
|------------------------------------------------------------------------------------------------------------------|
| Vertex: e- @ (x = 0.00000 m, y = 0.00000 m, z = 0.00000 m, t = 0.000000e+00 s) |
|------------------------------------------------------------------------------------------------------------------|
| Err flag [bits:15->0] : 0000000000000000 | 1st set: none |
| Err mask [bits:15->0] : 1111111111111111 | Is unphysical: NO | Accepted: YES |
|------------------------------------------------------------------------------------------------------------------|
| sig(Ev) = 3.02243e-30 cm^2 | d2sig(x,y;E)/dxdy = 3.88497e-29 cm^2 | Weight = 1.00000 |
|------------------------------------------------------------------------------------------------------------------|
--------------------------------------------------------------------------------------------------------------
GENIE Interaction Summary
--------------------------------------------------------------------------------------------------------------
[-] [Init-State]
|--> probe : PDG-code = 11 (e-)
|--> nucl. target : Z = 6, A = 12, PDG-Code = 1000060120 (C12)
|--> hit nucleon : PDC-Code = 2112 (neutron)
|--> hit quark : PDC-Code = 2 (u) [valence]
|--> probe 4P : (E = 3.000000, Px = 0.000000, Py = 0.000000, Pz = 3.000000)
|--> target 4P : (E = 11.174863, Px = 0.000000, Py = 0.000000, Pz = 0.000000)
|--> nucleon 4P : (E = 0.920077, Px = 0.119469, Py = -0.001260, Pz = 0.038037)
[-] [Process-Info]
|--> Interaction : EM
|--> Scattering : DIS
[-] [Kinematics]
|--> *Running* Hadronic invariant mass W = 2.051104
|--> *Selected* Bjorken x = 0.014575
|--> *Selected* Inelasticity y = 0.647390
|--> *Selected* Momentum transfer Q2 (>0) = 0.049934
|--> *Selected* Hadronic invariant mass W = 2.051104
[-] [Exclusive Process Info]
|--> charm prod. : false |--> strange prod. : false
|--> f/s nucleons : N(p) = 0 N(n) = 0
|--> f/s pions : N(pi^0) = 0 N(pi^+) = 0 N(pi^-) = 0
|--> f/s Other : N(gamma) = 0 N(Rho^0) = 0 N(Rho^+) = 0 N(Rho^-) = 0
|--> resonance : [not set]
|--> final quark prod. : false
|--> final lepton prod. : false
--------------------------------------------------------------------------------------------------------------
Particle with index = 16 has mothers from 7 to 10, but particles with indexes 8 and 9 couldn't be its mothers, because each of them have only one child: 8 has child 13 and 9 has child 14.
idkakorin commented
Let's compare the outputs of HINCLCascadeIntranuke and HG4BertCascIntranuke models.
Two first events are generated with HINCLCascadeIntranuke:
|------------------------------------------------------------------------------------------------------------------|
|GENIE GHEP Event Record [print level: 3] |
|------------------------------------------------------------------------------------------------------------------|
| Idx | Name | Ist | PDG | Mother | Daughter | Px | Py | Pz | E | m |
|------------------------------------------------------------------------------------------------------------------|
| 0 | e- | 0 | 11 | -1 | -1 | 4 | 4 | 0.000 | 0.000 | 3.000 | 3.000 | 0.001 |
| 1 | C12 | 0 | 1000060120 | -1 | -1 | 2 | 3 | 0.000 | 0.000 | 0.000 | 11.175 | 11.175 |
| 2 | proton | 11 | 2212 | 1 | -1 | 5 | 5 | -0.064 | 0.112 | -0.102 | 0.921 | **0.938 | M = 0.906
| 3 | B11 | 2 | 1000050110 | 1 | -1 | 10 | 10 | 0.064 | -0.112 | 0.102 | 10.254 | 10.253 |
| 4 | e- | 1 | 11 | 0 | -1 | -1 | -1 | -0.172 | -0.185 | 2.920 | 2.931 | 0.001 | P = (0.059,0.063,-0.996)
| 5 | proton | 14 | 2212 | 2 | -1 | 6 | 9 | 0.108 | 0.296 | -0.022 | 0.990 | 0.938 |
| 6 | H2 | 1 | 1000010020 | 5 | -1 | -1 | -1 | 0.162 | 0.216 | -0.121 | 1.899 | 1.876 |
| 7 | proton | 1 | 2212 | 5 | -1 | -1 | -1 | -0.015 | -0.064 | -0.067 | 0.943 | 0.938 |
| 8 | neutron | 1 | 2112 | 5 | -1 | -1 | -1 | 0.017 | 0.031 | 0.057 | 0.942 | 0.940 |
| 9 | He4 | 1 | 1000020040 | 5 | -1 | -1 | -1 | 0.051 | -0.082 | 0.120 | 3.731 | 3.727 |
| 10 | He4 | 1 | 1000020040 | 3 | -1 | -1 | -1 | -0.044 | 0.083 | 0.090 | 3.730 | 3.727 |
|------------------------------------------------------------------------------------------------------------------|
| Fin-Init: | 0.000 | -0.000 | -0.000 | -0.000 | |
|------------------------------------------------------------------------------------------------------------------|
| Vertex: e- @ (x = 0.00000 m, y = 0.00000 m, z = 0.00000 m, t = 0.000000e+00 s) |
|------------------------------------------------------------------------------------------------------------------|
| Err flag [bits:15->0] : 0000000000000000 | 1st set: none |
| Err mask [bits:15->0] : 1111111111111111 | Is unphysical: NO | Accepted: YES |
|------------------------------------------------------------------------------------------------------------------|
| sig(Ev) = 3.37581e-29 cm^2 | dsig(Q2;E)/dQ2 = 4.52020e-29 cm^2/GeV^2 | Weight = 1.00000 |
|------------------------------------------------------------------------------------------------------------------|
--------------------------------------------------------------------------------------------------------------
GENIE Interaction Summary
--------------------------------------------------------------------------------------------------------------
[-] [Init-State]
|--> probe : PDG-code = 11 (e-)
|--> nucl. target : Z = 6, A = 12, PDG-Code = 1000060120 (C12)
|--> hit nucleon : PDC-Code = 2212 (proton)
|--> hit quark : no set
|--> probe 4P : (E = 3.000000, Px = 0.000000, Py = 0.000000, Pz = 3.000000)
|--> target 4P : (E = 11.174863, Px = 0.000000, Py = 0.000000, Pz = 0.000000)
|--> nucleon 4P : (E = 0.921003, Px = -0.063703, Py = 0.111874, Pz = -0.101640)
[-] [Process-Info]
|--> Interaction : EM
|--> Scattering : QES
[-] [Kinematics]
|--> *Selected* Bjorken x = 0.525450
|--> *Selected* Inelasticity y = 0.020267
|--> *Selected* Momentum transfer Q2 (>0) = 0.065342
|--> *Selected* Hadronic invariant mass W = 0.938272
[-] [Exclusive Process Info]
|--> charm prod. : false |--> strange prod. : false
|--> f/s nucleons : N(p) = 0 N(n) = 0
|--> f/s pions : N(pi^0) = 0 N(pi^+) = 0 N(pi^-) = 0
|--> f/s Other : N(gamma) = 0 N(Rho^0) = 0 N(Rho^+) = 0 N(Rho^-) = 0
|--> resonance : [not set]
|--> final quark prod. : false
|--> final lepton prod. : false
--------------------------------------------------------------------------------------------------------------
|------------------------------------------------------------------------------------------------------------------|
|GENIE GHEP Event Record [print level: 3] |
|------------------------------------------------------------------------------------------------------------------|
| Idx | Name | Ist | PDG | Mother | Daughter | Px | Py | Pz | E | m |
|------------------------------------------------------------------------------------------------------------------|
| 0 | e- | 0 | 11 | -1 | -1 | 4 | 4 | 0.000 | 0.000 | 3.000 | 3.000 | 0.001 |
| 1 | C12 | 0 | 1000060120 | -1 | -1 | 2 | 3 | 0.000 | 0.000 | 0.000 | 11.175 | 11.175 |
| 2 | proton | 11 | 2212 | 1 | -1 | 5 | 5 | -0.020 | -0.029 | 0.075 | 0.922 | **0.938 | M = 0.918
| 3 | B11 | 2 | 1000050110 | 1 | -1 | 11 | 11 | 0.020 | 0.029 | -0.075 | 10.253 | 10.253 |
| 4 | e- | 1 | 11 | 0 | -1 | -1 | -1 | 0.367 | -0.098 | 2.856 | 2.881 | 0.001 | P = (-0.128,0.034,-0.991)
| 5 | proton | 14 | 2212 | 2 | -1 | 6 | 10 | -0.388 | 0.069 | 0.219 | 1.041 | 0.938 |
| 6 | neutron | 1 | 2112 | 5 | -1 | -1 | -1 | 0.079 | 0.153 | 0.120 | 0.963 | 0.940 |
| 7 | proton | 1 | 2212 | 5 | -1 | -1 | -1 | -0.057 | -0.084 | -0.098 | 0.949 | 0.938 |
| 8 | proton | 1 | 2212 | 5 | -1 | -1 | -1 | -0.066 | 0.008 | -0.019 | 0.941 | 0.938 |
| 9 | H2 | 1 | 1000010020 | 5 | -1 | -1 | -1 | -0.019 | 0.201 | 0.076 | 1.888 | 1.876 |
| 10 | H3 | 1 | 1000010030 | 5 | -1 | -1 | -1 | -0.220 | -0.075 | 0.116 | 2.821 | 2.809 |
| 11 | He4 | 1 | 1000020040 | 3 | -1 | -1 | -1 | -0.084 | -0.105 | -0.051 | 3.730 | 3.727 |
|------------------------------------------------------------------------------------------------------------------|
| Fin-Init: | -0.000 | 0.000 | -0.000 | -0.002 | |
|------------------------------------------------------------------------------------------------------------------|
| Vertex: e- @ (x = 0.00000 m, y = 0.00000 m, z = 0.00000 m, t = 0.000000e+00 s) |
|------------------------------------------------------------------------------------------------------------------|
| Err flag [bits:15->0] : 0000000000000000 | 1st set: none |
| Err mask [bits:15->0] : 1111111111111111 | Is unphysical: NO | Accepted: YES |
|------------------------------------------------------------------------------------------------------------------|
| sig(Ev) = 3.37581e-29 cm^2 | dsig(Q2;E)/dQ2 = 6.09729e-30 cm^2/GeV^2 | Weight = 1.00000 |
|------------------------------------------------------------------------------------------------------------------|
--------------------------------------------------------------------------------------------------------------
GENIE Interaction Summary
--------------------------------------------------------------------------------------------------------------
[-] [Init-State]
|--> probe : PDG-code = 11 (e-)
|--> nucl. target : Z = 6, A = 12, PDG-Code = 1000060120 (C12)
|--> hit nucleon : PDC-Code = 2212 (proton)
|--> hit quark : no set
|--> probe 4P : (E = 3.000000, Px = 0.000000, Py = 0.000000, Pz = 3.000000)
|--> target 4P : (E = 11.174863, Px = 0.000000, Py = 0.000000, Pz = 0.000000)
|--> nucleon 4P : (E = 0.921982, Px = -0.020381, Py = -0.028634, Pz = 0.074807)
[-] [Process-Info]
|--> Interaction : EM
|--> Scattering : QES
[-] [Kinematics]
|--> *Selected* Bjorken x = 0.802749
|--> *Selected* Inelasticity y = 0.037037
|--> *Selected* Momentum transfer Q2 (>0) = 0.151128
|--> *Selected* Hadronic invariant mass W = 0.938272
[-] [Exclusive Process Info]
|--> charm prod. : false |--> strange prod. : false
|--> f/s nucleons : N(p) = 0 N(n) = 0
|--> f/s pions : N(pi^0) = 0 N(pi^+) = 0 N(pi^-) = 0
|--> f/s Other : N(gamma) = 0 N(Rho^0) = 0 N(Rho^+) = 0 N(Rho^-) = 0
|--> resonance : [not set]
|--> final quark prod. : false
|--> final lepton prod. : false
--------------------------------------------------------------------------------------------------------------
The following one with HG4BertCascIntranuke:
|------------------------------------------------------------------------------------------------------------------|
|GENIE GHEP Event Record [print level: 3] |
|------------------------------------------------------------------------------------------------------------------|
| Idx | Name | Ist | PDG | Mother | Daughter | Px | Py | Pz | E | m |
|------------------------------------------------------------------------------------------------------------------|
| 0 | e- | 0 | 11 | -1 | -1 | 4 | 4 | 0.000 | 0.000 | 3.000 | 3.000 | 0.001 |
| 1 | C12 | 0 | 1000060120 | -1 | -1 | 2 | 3 | 0.000 | 0.000 | 0.000 | 11.175 | 11.175 |
| 2 | proton | 11 | 2212 | 1 | -1 | 5 | 5 | 0.044 | -0.038 | 0.008 | 0.922 | **0.938 | M = 0.920
| 3 | B11 | 2 | 1000050110 | 1 | -1 | 9 | 10 | -0.044 | 0.038 | -0.008 | 10.253 | 10.253 |
| 4 | e- | 1 | 11 | 0 | -1 | -1 | -1 | 0.057 | -0.133 | 2.975 | 2.978 | 0.001 | P = (-0.019,0.045,-0.999)
| 5 | proton | 14 | 2212 | 2 | -1 | 6 | 8 | -0.013 | 0.095 | 0.033 | 0.944 | 0.938 |
| 6 | gamma | 1 | 22 | 5 | -1 | -1 | -1 | -0.001 | 0.000 | -0.001 | 0.002 | 0.000 |
| 7 | gamma | 1 | 22 | 5 | -1 | -1 | -1 | 0.001 | -0.000 | -0.000 | 0.001 | 0.000 |
| 8 | gamma | 1 | 22 | 5 | -1 | -1 | -1 | -0.000 | 0.000 | 0.000 | 0.000 | 0.000 |
| 9 | He4 | 1 | 1000020040 | 3 | -1 | -1 | -1 | -0.032 | 0.225 | -0.103 | 3.736 | 3.727 |
| 10 | Be8 | 1 | 1000040080 | 3 | -1 | -1 | -1 | -0.024 | -0.093 | 0.130 | 7.457 | 7.455 |
|------------------------------------------------------------------------------------------------------------------|
| Fin-Init: | 0.000 | -0.000 | 0.000 | -0.000 | |
|------------------------------------------------------------------------------------------------------------------|
| Vertex: e- @ (x = 0.00000 m, y = 0.00000 m, z = 0.00000 m, t = 0.000000e+00 s) |
|------------------------------------------------------------------------------------------------------------------|
| Err flag [bits:15->0] : 0000000000000000 | 1st set: none |
| Err mask [bits:15->0] : 1111111111111111 | Is unphysical: NO | Accepted: YES |
|------------------------------------------------------------------------------------------------------------------|
| sig(Ev) = 3.37582e-29 cm^2 | dsig(Q2;E)/dQ2 = 5.31281e-28 cm^2/GeV^2 | Weight = 1.00000 |
|------------------------------------------------------------------------------------------------------------------|
--------------------------------------------------------------------------------------------------------------
GENIE Interaction Summary
--------------------------------------------------------------------------------------------------------------
[-] [Init-State]
|--> probe : PDG-code = 11 (e-)
|--> nucl. target : Z = 6, A = 12, PDG-Code = 1000060120 (C12)
|--> hit nucleon : PDC-Code = 2212 (proton)
|--> hit quark : no set
|--> probe 4P : (E = 3.000000, Px = 0.000000, Py = 0.000000, Pz = 3.000000)
|--> target 4P : (E = 11.174863, Px = 0.000000, Py = 0.000000, Pz = 0.000000)
|--> nucleon 4P : (E = 0.922149, Px = 0.043767, Py = -0.037856, Pz = 0.007689)
[-] [Process-Info]
|--> Interaction : EM
|--> Scattering : QES
[-] [Kinematics]
|--> *Selected* Bjorken x = 0.386032
|--> *Selected* Inelasticity y = 0.009916
|--> *Selected* Momentum transfer Q2 (>0) = 0.021002
|--> *Selected* Hadronic invariant mass W = 0.938272
[-] [Exclusive Process Info]
|--> charm prod. : false |--> strange prod. : false
|--> f/s nucleons : N(p) = 0 N(n) = 0
|--> f/s pions : N(pi^0) = 0 N(pi^+) = 0 N(pi^-) = 0
|--> f/s Other : N(gamma) = 0 N(Rho^0) = 0 N(Rho^+) = 0 N(Rho^-) = 0
|--> resonance : [not set]
|--> final quark prod. : false
|--> final lepton prod. : false
--------------------------------------------------------------------------------------------------------------
The output of HG4BertCascIntranuke seems more logigal for me.
idkakorin commented
Resonance events:
HINCLCascadeIntranuke:
|------------------------------------------------------------------------------------------------------------------|
|GENIE GHEP Event Record [print level: 3] |
|------------------------------------------------------------------------------------------------------------------|
| Idx | Name | Ist | PDG | Mother | Daughter | Px | Py | Pz | E | m |
|------------------------------------------------------------------------------------------------------------------|
| 0 | e- | 0 | 11 | -1 | -1 | 4 | 4 | 0.000 | 0.000 | 3.000 | 3.000 | 0.001 |
| 1 | C12 | 0 | 1000060120 | -1 | -1 | 2 | 3 | 0.000 | 0.000 | 0.000 | 11.175 | 11.175 |
| 2 | proton | 11 | 2212 | 1 | -1 | 5 | 5 | 0.020 | -0.144 | -0.100 | 0.921 | **0.938 | M = 0.904
| 3 | B11 | 2 | 1000050110 | 1 | -1 | 11 | 12 | -0.020 | 0.144 | 0.100 | 10.254 | 10.253 |
| 4 | e- | 1 | 11 | 0 | -1 | -1 | -1 | 0.287 | -0.187 | 1.956 | 1.986 | 0.001 | P = (-0.145,0.094,-0.985)
| 5 | N+(1720) | 3 | 202214 | 2 | -1 | 6 | 7 | -0.267 | 0.044 | 0.943 | 1.935 | **1.720 | M = 1.667
| 6 | proton | 14 | 2212 | 5 | -1 | 8 | 8 | -0.716 | 0.178 | 0.748 | 1.409 | 0.938 | FSI = 1
| 7 | pi0 | 14 | 111 | 5 | -1 | 9 | 10 | 0.449 | -0.135 | 0.195 | 0.526 | 0.135 |
| 8 | proton | 1 | 2212 | 6 | -1 | -1 | -1 | -0.716 | 0.178 | 0.748 | 1.409 | 0.938 |
| 9 | pi0 | 1 | 111 | 7 | -1 | -1 | -1 | 0.171 | 0.002 | 0.403 | 0.458 | 0.135 |
| 10 | neutron | 1 | 2112 | 7 | -1 | -1 | -1 | 0.213 | -0.196 | -0.107 | 0.989 | 0.940 |
| 11 | He4 | 1 | 1000020040 | 3 | -1 | -1 | -1 | -0.025 | 0.109 | -0.070 | 3.730 | 3.727 |
| 12 | Li6 | 1 | 1000030060 | 3 | -1 | -1 | -1 | 0.070 | 0.094 | 0.070 | 5.603 | 5.602 |
|------------------------------------------------------------------------------------------------------------------|
| Fin-Init: | 0.000 | -0.000 | 0.000 | -0.000 | |
|------------------------------------------------------------------------------------------------------------------|
| Vertex: e- @ (x = 0.00000 m, y = 0.00000 m, z = 0.00000 m, t = 0.000000e+00 s) |
|------------------------------------------------------------------------------------------------------------------|
| Err flag [bits:15->0] : 0000000000000000 | 1st set: none |
| Err mask [bits:15->0] : 1111111111111111 | Is unphysical: NO | Accepted: YES |
|------------------------------------------------------------------------------------------------------------------|
| sig(Ev) = 1.72623e-31 cm^2 | d2sig(W,Q2;E)/dWdQ2 = 2.39707e-30 cm^2/GeV^3 | Weight = 1.00000 |
|------------------------------------------------------------------------------------------------------------------|
--------------------------------------------------------------------------------------------------------------
GENIE Interaction Summary
--------------------------------------------------------------------------------------------------------------
[-] [Init-State]
|--> probe : PDG-code = 11 (e-)
|--> nucl. target : Z = 6, A = 12, PDG-Code = 1000060120 (C12)
|--> hit nucleon : PDC-Code = 2212 (proton)
|--> hit quark : no set
|--> probe 4P : (E = 3.000000, Px = 0.000000, Py = 0.000000, Pz = 3.000000)
|--> target 4P : (E = 11.174863, Px = 0.000000, Py = 0.000000, Pz = 0.000000)
|--> nucleon 4P : (E = 0.920798, Px = 0.020237, Py = -0.143581, Pz = -0.100412)
[-] [Process-Info]
|--> Interaction : EM
|--> Scattering : RES
[-] [Kinematics]
|--> *Selected* Bjorken x = 0.083466
|--> *Selected* Inelasticity y = 0.349569
|--> *Selected* Momentum transfer Q2 (>0) = 0.178776
|--> *Selected* Hadronic invariant mass W = 1.667299
[-] [Exclusive Process Info]
|--> charm prod. : false |--> strange prod. : false
|--> f/s nucleons : N(p) = 0 N(n) = 0
|--> f/s pions : N(pi^0) = 0 N(pi^+) = 0 N(pi^-) = 0
|--> f/s Other : N(gamma) = 0 N(Rho^0) = 0 N(Rho^+) = 0 N(Rho^-) = 0
|--> resonance : P13(1720)
|--> final quark prod. : false
|--> final lepton prod. : false
--------------------------------------------------------------------------------------------------------------
HG4BertCascIntranuke:
|------------------------------------------------------------------------------------------------------------------|
|GENIE GHEP Event Record [print level: 3] |
|------------------------------------------------------------------------------------------------------------------|
| Idx | Name | Ist | PDG | Mother | Daughter | Px | Py | Pz | E | m |
|------------------------------------------------------------------------------------------------------------------|
| 0 | e- | 0 | 11 | -1 | -1 | 4 | 4 | 0.000 | 0.000 | 3.000 | 3.000 | 0.001 |
| 1 | C12 | 0 | 1000060120 | -1 | -1 | 2 | 3 | 0.000 | 0.000 | 0.000 | 11.175 | 11.175 |
| 2 | neutron | 11 | 2112 | 1 | -1 | 5 | 5 | 0.188 | 0.036 | 0.109 | 0.918 | **0.940 | M = 0.892
| 3 | C11 | 2 | 1000060110 | 1 | -1 | 19 | 19 | -0.188 | -0.036 | -0.109 | 10.256 | 10.254 |
| 4 | e- | 1 | 11 | 0 | -1 | -1 | -1 | 0.018 | -0.137 | 2.519 | 2.522 | 0.001 | P = (-0.007,0.054,-0.998)
| 5 | Delta0 | 3 | 2114 | 2 | -1 | 6 | 7 | 0.171 | 0.173 | 0.591 | 1.396 | **1.233 | M = 1.241
| 6 | neutron | 14 | 2112 | 5 | -1 | 8 | 18 | 0.367 | 0.092 | 0.489 | 1.124 | 0.940 |
| 7 | pi0 | 14 | 111 | 5 | -1 | 8 | 18 | -0.196 | 0.082 | 0.102 | 0.272 | 0.135 |
| 8 | neutron | 1 | 2112 | 6 | 7 | -1 | -1 | 0.413 | -0.071 | 0.055 | 1.030 | 0.940 |
| 9 | pi+ | 1 | 211 | 6 | 7 | -1 | -1 | -0.056 | -0.157 | 0.017 | 0.218 | 0.140 |
| 10 | neutron | 1 | 2112 | 6 | 7 | -1 | -1 | -0.024 | 0.135 | 0.159 | 0.963 | 0.940 |
| 11 | neutron | 1 | 2112 | 6 | 7 | -1 | -1 | 0.067 | 0.069 | 0.014 | 0.945 | 0.940 |
| 12 | proton | 1 | 2212 | 6 | 7 | -1 | -1 | -0.019 | 0.185 | 0.207 | 0.979 | 0.938 |
| 13 | proton | 1 | 2212 | 6 | 7 | -1 | -1 | -0.039 | 0.071 | 0.006 | 0.942 | 0.938 |
| 14 | neutron | 1 | 2112 | 6 | 7 | -1 | -1 | 0.052 | -0.068 | 0.028 | 0.944 | 0.940 |
| 15 | neutron | 1 | 2112 | 6 | 7 | -1 | -1 | -0.058 | -0.082 | 0.037 | 0.946 | 0.940 |
| 16 | proton | 1 | 2212 | 6 | 7 | -1 | -1 | 0.001 | -0.046 | 0.070 | 0.942 | 0.938 |
| 17 | gamma | 1 | 22 | 6 | 7 | -1 | -1 | 0.003 | -0.004 | 0.006 | 0.008 | 0.000 |
| 18 | gamma | 1 | 22 | 6 | 7 | -1 | -1 | -0.001 | 0.001 | 0.000 | 0.001 | 0.000 |
| 19 | He4 | 1 | 1000020040 | 3 | -1 | -1 | -1 | -0.358 | 0.104 | -0.119 | 3.734 | **3.727 | M = 3.713
|------------------------------------------------------------------------------------------------------------------|
| Fin-Init: | 0.000 | -0.000 | 0.000 | -0.002 | |
|------------------------------------------------------------------------------------------------------------------|
| Vertex: e- @ (x = 0.00000 m, y = 0.00000 m, z = 0.00000 m, t = 0.000000e+00 s) |
|------------------------------------------------------------------------------------------------------------------|
| Err flag [bits:15->0] : 0000000000000000 | 1st set: none |
| Err mask [bits:15->0] : 1111111111111111 | Is unphysical: NO | Accepted: YES |
|------------------------------------------------------------------------------------------------------------------|
| sig(Ev) = 3.41717e-30 cm^2 | d2sig(W,Q2;E)/dWdQ2 = 3.40509e-28 cm^2/GeV^3 | Weight = 1.00000 |
|------------------------------------------------------------------------------------------------------------------|
--------------------------------------------------------------------------------------------------------------
GENIE Interaction Summary
--------------------------------------------------------------------------------------------------------------
[-] [Init-State]
|--> probe : PDG-code = 11 (e-)
|--> nucl. target : Z = 6, A = 12, PDG-Code = 1000060120 (C12)
|--> hit nucleon : PDC-Code = 2112 (neutron)
|--> hit quark : no set
|--> probe 4P : (E = 3.000000, Px = 0.000000, Py = 0.000000, Pz = 3.000000)
|--> target 4P : (E = 11.174863, Px = 0.000000, Py = 0.000000, Pz = 0.000000)
|--> nucleon 4P : (E = 0.918470, Px = 0.188142, Py = 0.036337, Pz = 0.109403)
[-] [Process-Info]
|--> Interaction : EM
|--> Scattering : RES
[-] [Kinematics]
|--> *Selected* Bjorken x = 0.029574
|--> *Selected* Inelasticity y = 0.158322
|--> *Selected* Momentum transfer Q2 (>0) = 0.022729
|--> *Selected* Hadronic invariant mass W = 1.241262
[-] [Exclusive Process Info]
|--> charm prod. : false |--> strange prod. : false
|--> f/s nucleons : N(p) = 0 N(n) = 0
|--> f/s pions : N(pi^0) = 0 N(pi^+) = 0 N(pi^-) = 0
|--> f/s Other : N(gamma) = 0 N(Rho^0) = 0 N(Rho^+) = 0 N(Rho^-) = 0
|--> resonance : P33(1232)
|--> final quark prod. : false
|--> final lepton prod. : false
--------------------------------------------------------------------------------------------------------------
Again, the output of HG4BertCascIntranuke seems more logigal for me.
idkakorin commented
MEC events:
HINCLCascadeIntranuke:
|------------------------------------------------------------------------------------------------------------------|
|GENIE GHEP Event Record [print level: 3] |
|------------------------------------------------------------------------------------------------------------------|
| Idx | Name | Ist | PDG | Mother | Daughter | Px | Py | Pz | E | m |
|------------------------------------------------------------------------------------------------------------------|
| 0 | e- | 0 | 11 | -1 | -1 | 4 | 4 | 0.000 | 0.000 | 3.000 | 3.000 | 0.001 |
| 1 | C12 | 0 | 1000060120 | -1 | -1 | 2 | 3 | 0.000 | 0.000 | 0.000 | 11.175 | 11.175 |
| 2 | n+n | 11 | 2000000200 | 1 | -1 | 5 | 5 | 0.194 | -0.056 | -0.268 | 1.909 | 1.879 |
| 3 | C10 | 2 | 1000060100 | 1 | -1 | 14 | 15 | -0.194 | 0.056 | 0.268 | 9.266 | **9.328 | M = 9.260
| 4 | e- | 1 | 11 | 0 | -1 | -1 | -1 | 0.137 | -0.344 | 2.638 | 2.664 | 0.001 | P = (-0.052,0.129,-0.990)
| 5 | n+n | 3 | 2000000200 | 2 | -1 | 6 | 7 | 0.057 | 0.288 | 0.094 | 2.245 | **1.879 | M = 2.224
| 6 | neutron | 14 | 2112 | 5 | -1 | 8 | 13 | -0.180 | -0.415 | -0.017 | 1.043 | 0.940 |
| 7 | neutron | 14 | 2112 | 5 | -1 | 8 | 13 | 0.237 | 0.703 | 0.111 | 1.202 | 0.940 |
| 8 | neutron | 1 | 2112 | 6 | 7 | -1 | -1 | -0.010 | -0.251 | -0.184 | 0.990 | 0.940 |
| 9 | neutron | 1 | 2112 | 6 | 7 | -1 | -1 | 0.073 | 0.561 | 0.387 | 1.163 | 0.940 |
| 10 | neutron | 1 | 2112 | 6 | 7 | -1 | -1 | 0.123 | 0.070 | -0.050 | 0.951 | 0.940 |
| 11 | proton | 1 | 2212 | 6 | 7 | -1 | -1 | -0.070 | -0.035 | 0.168 | 0.956 | 0.938 |
| 12 | proton | 1 | 2212 | 6 | 7 | -1 | -1 | -0.033 | 0.062 | 0.025 | 0.941 | 0.938 |
| 13 | proton | 1 | 2212 | 6 | 7 | -1 | -1 | 0.004 | 0.045 | -0.064 | 0.942 | 0.938 |
| 14 | H2 | 1 | 1000010020 | 3 | -1 | -1 | -1 | -0.175 | -0.125 | 0.129 | 1.892 | 1.876 |
| 15 | He4 | 1 | 1000020040 | 3 | -1 | -1 | -1 | -0.048 | 0.017 | -0.047 | 3.728 | 3.727 |
|------------------------------------------------------------------------------------------------------------------|
| Fin-Init: | -0.000 | -0.000 | -0.000 | 0.052 | |
|------------------------------------------------------------------------------------------------------------------|
| Vertex: e- @ (x = 0.00000 m, y = 0.00000 m, z = 0.00000 m, t = 0.000000e+00 s) |
|------------------------------------------------------------------------------------------------------------------|
| Err flag [bits:15->0] : 0000000000000000 | 1st set: none |
| Err mask [bits:15->0] : 1111111111111111 | Is unphysical: NO | Accepted: YES |
|------------------------------------------------------------------------------------------------------------------|
| sig(Ev) = 8.92144e-32 cm^2 | dsig(Ev;{K_s})/dK = 0.00000e+00 cm^2/{K} | Weight = 1.00000 |
|------------------------------------------------------------------------------------------------------------------|
--------------------------------------------------------------------------------------------------------------
GENIE Interaction Summary
--------------------------------------------------------------------------------------------------------------
[-] [Init-State]
|--> probe : PDG-code = 11 (e-)
|--> nucl. target : Z = 6, A = 12, PDG-Code = 1000060120 (C12)
|--> hit nucleon : no set
|--> hit quark : no set
|--> probe 4P : (E = 3.000000, Px = 0.000000, Py = 0.000000, Pz = 3.000000)
|--> target 4P : (E = 11.174863, Px = 0.000000, Py = 0.000000, Pz = 0.000000)
[-] [Process-Info]
|--> Interaction : EM
|--> Scattering : MEC
[-] [Kinematics]
|--> *Selected* Bjorken x = 0.098832
|--> *Selected* Inelasticity y = 0.120152
|--> *Selected* Momentum transfer Q2 (>0) = 0.155110
|--> *Selected* Hadronic invariant mass W = 2.223811
[-] [Exclusive Process Info]
|--> charm prod. : false |--> strange prod. : false
|--> f/s nucleons : N(p) = 0 N(n) = 0
|--> f/s pions : N(pi^0) = 0 N(pi^+) = 0 N(pi^-) = 0
|--> f/s Other : N(gamma) = 0 N(Rho^0) = 0 N(Rho^+) = 0 N(Rho^-) = 0
|--> resonance : [not set]
|--> final quark prod. : false
|--> final lepton prod. : false
--------------------------------------------------------------------------------------------------------------
HG4BertCascIntranuke:
|------------------------------------------------------------------------------------------------------------------|
|GENIE GHEP Event Record [print level: 3] |
|------------------------------------------------------------------------------------------------------------------|
| Idx | Name | Ist | PDG | Mother | Daughter | Px | Py | Pz | E | m |
|------------------------------------------------------------------------------------------------------------------|
| 0 | e- | 0 | 11 | -1 | -1 | 4 | 4 | 0.000 | 0.000 | 3.000 | 3.000 | 0.001 |
| 1 | C12 | 0 | 1000060120 | -1 | -1 | 2 | 3 | 0.000 | 0.000 | 0.000 | 11.175 | 11.175 |
| 2 | n+p | 11 | 2000000201 | 1 | -1 | 5 | 5 | -0.056 | -0.096 | 0.096 | 1.884 | 1.878 |
| 3 | B10 | 2 | 1000050100 | 1 | -1 | 15 | 16 | 0.056 | 0.096 | -0.096 | 9.291 | **9.324 | M = 9.290
| 4 | e- | 1 | 11 | 0 | -1 | -1 | -1 | -0.286 | -0.462 | 2.537 | 2.594 | 0.001 | P = (0.110,0.178,-0.978)
| 5 | n+p | 3 | 2000000201 | 2 | -1 | 6 | 7 | 0.230 | 0.366 | 0.559 | 2.289 | **1.878 | M = 2.177
| 6 | neutron | 14 | 2112 | 5 | -1 | 8 | 14 | 0.362 | 0.659 | 0.475 | 1.294 | 0.940 |
| 7 | proton | 14 | 2212 | 5 | -1 | 8 | 14 | -0.132 | -0.293 | 0.084 | 0.995 | 0.938 |
| 8 | neutron | 1 | 2112 | 6 | 7 | -1 | -1 | 0.342 | 0.620 | 0.438 | 1.255 | 0.940 |
| 9 | proton | 1 | 2212 | 6 | 7 | -1 | -1 | -0.073 | -0.176 | 0.047 | 0.959 | 0.938 |
| 10 | proton | 1 | 2212 | 6 | 7 | -1 | -1 | -0.032 | 0.035 | 0.113 | 0.946 | 0.938 |
| 11 | neutron | 1 | 2112 | 6 | 7 | -1 | -1 | -0.106 | 0.149 | -0.064 | 0.959 | 0.940 |
| 12 | neutron | 1 | 2112 | 6 | 7 | -1 | -1 | 0.050 | -0.019 | 0.027 | 0.941 | 0.940 |
| 13 | proton | 1 | 2212 | 6 | 7 | -1 | -1 | 0.031 | -0.007 | -0.006 | 0.939 | 0.938 |
| 14 | gamma | 1 | 22 | 6 | 7 | -1 | -1 | -0.000 | -0.000 | -0.000 | 0.000 | 0.000 |
| 15 | H2 | 1 | 1000010020 | 3 | -1 | -1 | -1 | 0.043 | 0.012 | 0.071 | 1.878 | 1.876 |
| 16 | He4 | 1 | 1000020040 | 3 | -1 | -1 | -1 | 0.031 | -0.151 | -0.162 | 3.736 | **3.727 | M = 3.730
|------------------------------------------------------------------------------------------------------------------|
| Fin-Init: | -0.000 | -0.000 | -0.000 | 0.033 | |
|------------------------------------------------------------------------------------------------------------------|
| Vertex: e- @ (x = 0.00000 m, y = 0.00000 m, z = 0.00000 m, t = 0.000000e+00 s) |
|------------------------------------------------------------------------------------------------------------------|
| Err flag [bits:15->0] : 0000000000000000 | 1st set: none |
| Err mask [bits:15->0] : 1111111111111111 | Is unphysical: NO | Accepted: YES |
|------------------------------------------------------------------------------------------------------------------|
| sig(Ev) = 7.13715e-31 cm^2 | dsig(Ev;{K_s})/dK = 0.00000e+00 cm^2/{K} | Weight = 1.00000 |
|------------------------------------------------------------------------------------------------------------------|
--------------------------------------------------------------------------------------------------------------
GENIE Interaction Summary
--------------------------------------------------------------------------------------------------------------
[-] [Init-State]
|--> probe : PDG-code = 11 (e-)
|--> nucl. target : Z = 6, A = 12, PDG-Code = 1000060120 (C12)
|--> hit nucleon : no set
|--> hit quark : no set
|--> probe 4P : (E = 3.000000, Px = 0.000000, Py = 0.000000, Pz = 3.000000)
|--> target 4P : (E = 11.174863, Px = 0.000000, Py = 0.000000, Pz = 0.000000)
[-] [Process-Info]
|--> Interaction : EM
|--> Scattering : MEC
[-] [Kinematics]
|--> *Selected* Bjorken x = 0.221276
|--> *Selected* Inelasticity y = 0.145499
|--> *Selected* Momentum transfer Q2 (>0) = 0.345331
|--> *Selected* Hadronic invariant mass W = 2.177485
[-] [Exclusive Process Info]
|--> charm prod. : false |--> strange prod. : false
|--> f/s nucleons : N(p) = 0 N(n) = 0
|--> f/s pions : N(pi^0) = 0 N(pi^+) = 0 N(pi^-) = 0
|--> f/s Other : N(gamma) = 0 N(Rho^0) = 0 N(Rho^+) = 0 N(Rho^-) = 0
|--> resonance : [not set]
|--> final quark prod. : false
|--> final lepton prod. : false
--------------------------------------------------------------------------------------------------------------
The both looks good for me.
idkakorin commented
DIS events:
HINCLCascadeIntranuke:
|------------------------------------------------------------------------------------------------------------------|
|GENIE GHEP Event Record [print level: 3] |
|------------------------------------------------------------------------------------------------------------------|
| Idx | Name | Ist | PDG | Mother | Daughter | Px | Py | Pz | E | m |
|------------------------------------------------------------------------------------------------------------------|
| 0 | e- | 0 | 11 | -1 | -1 | 4 | 4 | 0.000 | 0.000 | 3.000 | 3.000 | 0.001 |
| 1 | C12 | 0 | 1000060120 | -1 | -1 | 2 | 3 | 0.000 | 0.000 | 0.000 | 11.175 | 11.175 |
| 2 | neutron | 11 | 2112 | 1 | -1 | 5 | 5 | 0.119 | -0.001 | 0.038 | 0.920 | **0.940 | M = 0.911
| 3 | C11 | 2 | 1000060110 | 1 | -1 | 21 | 22 | -0.119 | 0.001 | -0.038 | 10.255 | 10.254 |
| 4 | e- | 1 | 11 | 0 | -1 | -1 | -1 | 0.075 | 0.110 | 1.059 | 1.067 | 0.001 | P = (-0.070,-0.103,-0.992)
| 5 | HadrSyst | 12 | 2000000001 | 2 | -1 | 6 | 11 | 0.045 | -0.111 | 1.979 | 2.853 | **0.000 | M = 2.051
| 6 | neutron | 14 | 2112 | 5 | -1 | 12 | 12 | -0.220 | -0.044 | 0.605 | 1.140 | 0.940 | FSI = 1
| 7 | pi0 | 14 | 111 | 5 | -1 | 16 | 20 | 0.137 | 0.022 | 0.190 | 0.271 | 0.135 |
| 8 | pi0 | 14 | 111 | 5 | -1 | 13 | 13 | 0.051 | -0.005 | 0.311 | 0.343 | 0.135 | FSI = 1
| 9 | pi0 | 14 | 111 | 5 | -1 | 14 | 14 | 0.053 | 0.099 | 0.400 | 0.437 | 0.135 | FSI = 1
| 10 | pi+ | 14 | 211 | 5 | -1 | 16 | 20 | -0.057 | 0.040 | 0.063 | 0.168 | 0.140 |
| 11 | pi- | 14 | -211 | 5 | -1 | 15 | 15 | 0.080 | -0.223 | 0.409 | 0.493 | 0.140 | FSI = 1
| 12 | neutron | 1 | 2112 | 6 | -1 | -1 | -1 | -0.220 | -0.044 | 0.605 | 1.140 | 0.940 |
| 13 | pi0 | 1 | 111 | 8 | -1 | -1 | -1 | 0.051 | -0.005 | 0.311 | 0.343 | 0.135 |
| 14 | pi0 | 1 | 111 | 9 | -1 | -1 | -1 | 0.053 | 0.099 | 0.400 | 0.437 | 0.135 |
| 15 | pi- | 1 | -211 | 11 | -1 | -1 | -1 | 0.080 | -0.223 | 0.409 | 0.493 | 0.140 |
| 16 | pi- | 1 | -211 | 7 | 10 | -1 | -1 | 0.026 | 0.151 | 0.074 | 0.220 | 0.140 |
| 17 | proton | 1 | 2212 | 7 | 10 | -1 | -1 | 0.086 | 0.077 | 0.138 | 0.955 | 0.938 |
| 18 | proton | 1 | 2212 | 7 | 10 | -1 | -1 | 0.028 | -0.009 | -0.294 | 0.984 | 0.938 |
| 19 | proton | 1 | 2212 | 7 | 10 | -1 | -1 | -0.107 | -0.012 | 0.029 | 0.945 | 0.938 |
| 20 | proton | 1 | 2212 | 7 | 10 | -1 | -1 | 0.003 | -0.020 | -0.119 | 0.946 | 0.938 |
| 21 | He3 | 1 | 1000020030 | 3 | -1 | -1 | -1 | 0.002 | -0.253 | 0.713 | 2.909 | 2.808 |
| 22 | He4 | 1 | 1000020040 | 3 | -1 | -1 | -1 | -0.078 | 0.129 | -0.326 | 3.745 | 3.727 |
|------------------------------------------------------------------------------------------------------------------|
| Fin-Init: | 0.000 | -0.000 | 0.000 | 0.009 | |
|------------------------------------------------------------------------------------------------------------------|
| Vertex: e- @ (x = 0.00000 m, y = 0.00000 m, z = 0.00000 m, t = 0.000000e+00 s) |
|------------------------------------------------------------------------------------------------------------------|
| Err flag [bits:15->0] : 0000000000000000 | 1st set: none |
| Err mask [bits:15->0] : 1111111111111111 | Is unphysical: NO | Accepted: YES |
|------------------------------------------------------------------------------------------------------------------|
| sig(Ev) = 3.02243e-30 cm^2 | d2sig(x,y;E)/dxdy = 3.88497e-29 cm^2 | Weight = 1.00000 |
|------------------------------------------------------------------------------------------------------------------|
--------------------------------------------------------------------------------------------------------------
GENIE Interaction Summary
--------------------------------------------------------------------------------------------------------------
[-] [Init-State]
|--> probe : PDG-code = 11 (e-)
|--> nucl. target : Z = 6, A = 12, PDG-Code = 1000060120 (C12)
|--> hit nucleon : PDC-Code = 2112 (neutron)
|--> hit quark : PDC-Code = 2 (u) [valence]
|--> probe 4P : (E = 3.000000, Px = 0.000000, Py = 0.000000, Pz = 3.000000)
|--> target 4P : (E = 11.174863, Px = 0.000000, Py = 0.000000, Pz = 0.000000)
|--> nucleon 4P : (E = 0.920077, Px = 0.119469, Py = -0.001260, Pz = 0.038037)
[-] [Process-Info]
|--> Interaction : EM
|--> Scattering : DIS
[-] [Kinematics]
|--> *Running* Hadronic invariant mass W = 2.051104
|--> *Selected* Bjorken x = 0.014575
|--> *Selected* Inelasticity y = 0.647390
|--> *Selected* Momentum transfer Q2 (>0) = 0.049934
|--> *Selected* Hadronic invariant mass W = 2.051104
[-] [Exclusive Process Info]
|--> charm prod. : false |--> strange prod. : false
|--> f/s nucleons : N(p) = 0 N(n) = 0
|--> f/s pions : N(pi^0) = 0 N(pi^+) = 0 N(pi^-) = 0
|--> f/s Other : N(gamma) = 0 N(Rho^0) = 0 N(Rho^+) = 0 N(Rho^-) = 0
|--> resonance : [not set]
|--> final quark prod. : false
|--> final lepton prod. : false
--------------------------------------------------------------------------------------------------------------
HG4BertCascIntranuke:
|------------------------------------------------------------------------------------------------------------------|
|GENIE GHEP Event Record [print level: 3] |
|------------------------------------------------------------------------------------------------------------------|
| Idx | Name | Ist | PDG | Mother | Daughter | Px | Py | Pz | E | m |
|------------------------------------------------------------------------------------------------------------------|
| 0 | e- | 0 | 11 | -1 | -1 | 4 | 4 | 0.000 | 0.000 | 3.000 | 3.000 | 0.001 |
| 1 | C12 | 0 | 1000060120 | -1 | -1 | 2 | 3 | 0.000 | 0.000 | 0.000 | 11.175 | 11.175 |
| 2 | neutron | 11 | 2112 | 1 | -1 | 5 | 5 | 0.062 | 0.150 | -0.087 | 0.919 | **0.940 | M = 0.900
| 3 | C11 | 2 | 1000060110 | 1 | -1 | 25 | 26 | -0.062 | -0.150 | 0.087 | 10.256 | 10.254 |
| 4 | e- | 1 | 11 | 0 | -1 | -1 | -1 | 0.070 | 0.082 | 0.595 | 0.605 | 0.001 | P = (-0.116,-0.136,-0.984)
| 5 | HadrSyst | 12 | 2000000001 | 2 | -1 | 6 | 11 | -0.008 | 0.068 | 2.318 | 3.314 | **0.000 | M = 2.368
| 6 | proton | 14 | 2212 | 5 | -1 | 12 | 12 | 0.018 | -0.178 | 1.451 | 1.737 | 0.938 | FSI = 1
| 7 | pi- | 14 | -211 | 5 | -1 | 15 | 24 | 0.022 | 0.132 | 0.289 | 0.347 | 0.140 |
| 8 | pi0 | 14 | 111 | 5 | -1 | 13 | 13 | -0.030 | 0.100 | 0.088 | 0.192 | 0.135 | FSI = 1
| 9 | pi0 | 14 | 111 | 5 | -1 | 14 | 14 | 0.190 | -0.220 | 0.366 | 0.487 | 0.135 | FSI = 1
| 10 | pi+ | 14 | 211 | 5 | -1 | 15 | 24 | -0.189 | 0.007 | -0.037 | 0.238 | 0.140 |
| 11 | pi- | 14 | -211 | 5 | -1 | 15 | 24 | -0.019 | 0.228 | 0.160 | 0.312 | 0.140 |
| 12 | proton | 1 | 2212 | 6 | -1 | -1 | -1 | 0.018 | -0.178 | 1.451 | 1.737 | 0.938 |
| 13 | pi0 | 1 | 111 | 8 | -1 | -1 | -1 | -0.030 | 0.100 | 0.088 | 0.192 | 0.135 |
| 14 | pi0 | 1 | 111 | 9 | -1 | -1 | -1 | 0.190 | -0.220 | 0.366 | 0.487 | 0.135 |
| 15 | pi- | 1 | -211 | 7 | 11 | -1 | -1 | 0.022 | 0.128 | 0.282 | 0.340 | 0.140 |
| 16 | proton | 1 | 2212 | 7 | 11 | -1 | -1 | -0.091 | 0.271 | 0.285 | 1.021 | 0.938 |
| 17 | pi- | 1 | -211 | 7 | 11 | -1 | -1 | 0.030 | -0.019 | -0.124 | 0.190 | 0.140 |
| 18 | pi+ | 1 | 211 | 7 | 11 | -1 | -1 | 0.017 | -0.038 | 0.100 | 0.177 | 0.140 |
| 19 | proton | 1 | 2212 | 7 | 11 | -1 | -1 | -0.140 | 0.035 | -0.121 | 0.957 | 0.938 |
| 20 | proton | 1 | 2212 | 7 | 11 | -1 | -1 | 0.025 | 0.054 | 0.016 | 0.940 | 0.938 |
| 21 | neutron | 1 | 2112 | 7 | 11 | -1 | -1 | 0.016 | 0.078 | 0.111 | 0.949 | 0.940 |
| 22 | gamma | 1 | 22 | 7 | 11 | -1 | -1 | 0.003 | -0.001 | 0.004 | 0.005 | 0.000 |
| 23 | gamma | 1 | 22 | 7 | 11 | -1 | -1 | 0.001 | 0.001 | -0.001 | 0.002 | 0.000 |
| 24 | gamma | 1 | 22 | 7 | 11 | -1 | -1 | -0.002 | 0.000 | -0.000 | 0.002 | 0.000 |
| 25 | H3 | 1 | 1000010030 | 3 | -1 | -1 | -1 | -0.310 | -0.166 | -0.051 | 2.831 | 2.809 |
| 26 | He4 | 1 | 1000020040 | 3 | -1 | -1 | -1 | 0.181 | -0.128 | -0.000 | 3.736 | **3.727 | M = 3.730
|------------------------------------------------------------------------------------------------------------------|
| Fin-Init: | -0.000 | -0.000 | 0.000 | -0.002 | |
|------------------------------------------------------------------------------------------------------------------|
| Vertex: e- @ (x = 0.00000 m, y = 0.00000 m, z = 0.00000 m, t = 0.000000e+00 s) |
|------------------------------------------------------------------------------------------------------------------|
| Err flag [bits:15->0] : 0000000000000000 | 1st set: none |
| Err mask [bits:15->0] : 1111111111111111 | Is unphysical: NO | Accepted: YES |
|------------------------------------------------------------------------------------------------------------------|
| sig(Ev) = 3.02243e-30 cm^2 | d2sig(x,y;E)/dxdy = 2.04566e-29 cm^2 | Weight = 1.00000 |
|------------------------------------------------------------------------------------------------------------------|
--------------------------------------------------------------------------------------------------------------
GENIE Interaction Summary
--------------------------------------------------------------------------------------------------------------
[-] [Init-State]
|--> probe : PDG-code = 11 (e-)
|--> nucl. target : Z = 6, A = 12, PDG-Code = 1000060120 (C12)
|--> hit nucleon : PDC-Code = 2112 (neutron)
|--> hit quark : PDC-Code = 2 (u) [valence]
|--> probe 4P : (E = 3.000000, Px = 0.000000, Py = 0.000000, Pz = 3.000000)
|--> target 4P : (E = 11.174863, Px = 0.000000, Py = 0.000000, Pz = 0.000000)
|--> nucleon 4P : (E = 0.919185, Px = 0.062292, Py = 0.150192, Pz = -0.087051)
[-] [Process-Info]
|--> Interaction : EM
|--> Scattering : DIS
[-] [Kinematics]
|--> *Running* Hadronic invariant mass W = 2.367954
|--> *Selected* Bjorken x = 0.012078
|--> *Selected* Inelasticity y = 0.804144
|--> *Selected* Momentum transfer Q2 (>0) = 0.058637
|--> *Selected* Hadronic invariant mass W = 2.367954
[-] [Exclusive Process Info]
|--> charm prod. : false |--> strange prod. : false
|--> f/s nucleons : N(p) = 0 N(n) = 0
|--> f/s pions : N(pi^0) = 0 N(pi^+) = 0 N(pi^-) = 0
|--> f/s Other : N(gamma) = 0 N(Rho^0) = 0 N(Rho^+) = 0 N(Rho^-) = 0
|--> resonance : [not set]
|--> final quark prod. : false
|--> final lepton prod. : false
--------------------------------------------------------------------------------------------------------------
The both outputs look similar.
idkakorin commented
My conclusion
is that events in the table should be rearranged in the proper way. The particle with status =14 and FSI = 1 should go first, after them should go particles with status = 14, then particles which mothers are ones with status =14 and FSI = 1, the particles which mothers are ones with status = 14 and in the end nuclei, which mother is nuclei with index 3.
I attached all outputs for the both models here.
HG4BertCascIntranuke_DIS.txt
HG4BertCascIntranuke_MEC.txt
HG4BertCascIntranuke_QE.txt
HG4BertCascIntranuke_RES.txt
HINCLCascadeIntranuke_DIS.txt
HINCLCascadeIntranuke_MEC.txt
HINCLCascadeIntranuke_QE.txt
HINCLCascadeIntranuke_RES.txt