search

GENIE-MC / Generator

The popular GENIE Generator product is used by nearly all accelerator neutrino experiments and it plays a key role in the exploitation of neutrino data. The Generator implements a modern software framework and it includes state-of-the-art physics modules. It captures the latest results of the GENIE global analysis of neutrino scattering data and includes several tunes that were produced using the proprietary Comparisons and Tuning products. The GENIE physics model is universal and comprehensive: It handles all neutrinos and targets, and all processes relevant from MeV to PeV energy scales. The Generator includes several tools (flux drivers, detector geometry navigators, specialized event generation apps, event reweighting engines) to simulate complex experimental setups in full detail and to support generator-related analysis tasks.
http://www.genie-mc.org
42 stars 92 forks source link

Igor/fsi fixes #355

Closed idkakorin closed 10 months ago

idkakorin commented 10 months ago

Combine updates of models HINCLCascadeIntranuke and HG4BertCascIntranuke in one branch to test them.

idkakorin commented 10 months ago

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 _G1810a/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
--------------------------------------------------------------------------------------------------------------
idkakorin commented 10 months ago

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 10 months ago

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.