Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Test of HINCLCascadeIntranuke and HG4BertCascIntranuke models. #359

Open
idkakorin opened this issue Nov 8, 2023 · 7 comments
Open

Test of HINCLCascadeIntranuke and HG4BertCascIntranuke models. #359

idkakorin opened this issue Nov 8, 2023 · 7 comments
Assignees

Comments

@idkakorin
Copy link
Contributor

idkakorin commented Nov 8, 2023

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 idkakorin changed the title 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. Test of HINCLCascadeIntranuke and HG4BertCascIntranuke models. Nov 8, 2023
@idkakorin
Copy link
Contributor Author

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
Copy link
Contributor Author

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
Copy link
Contributor Author

idkakorin commented Nov 12, 2023

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
Copy link
Contributor Author

idkakorin commented Nov 12, 2023

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
Copy link
Contributor Author

idkakorin commented Nov 12, 2023

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
Copy link
Contributor Author

idkakorin commented Nov 12, 2023

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
Copy link
Contributor Author

idkakorin commented Nov 12, 2023

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

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

7 participants