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.

Sign up for GitHub

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

Jump to bottom

Igor/fsi fixes #355

Merged
merged 2 commits into from
Nov 6, 2023
Merged

Igor/fsi fixes #355

merged 2 commits into from
Nov 6, 2023

Conversation

idkakorin
Copy link
Contributor

@idkakorin idkakorin commented Nov 6, 2023 •
edited
Loading

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

@idkakorin idkakorin merged commit f5abe45 into GENIE-MC:igor/FSIFixes Nov 6, 2023
idkakorin added a commit to idkakorin/Generator that referenced this pull request Nov 6, 2023
@idkakorin
Merge pull request #13 from GENIE-MC/igor/FSIFixes
d2b446b 
Merge pull request GENIE-MC#355 from idkakorin/igor/FSIFixes
@idkakorin
Copy link
Contributor Author

idkakorin commented Nov 8, 2023 •
edited
Loading

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
--------------------------------------------------------------------------------------------------------------

@idkakorin
Copy link
Contributor Author

idkakorin commented Nov 8, 2023 •
edited
Loading

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

idkakorin commented Nov 8, 2023 •
edited
Loading

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 idkakorin mentioned this pull request Nov 8, 2023
Open
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Reviewers
No reviews
Assignees
No one assigned
Labels
None yet
Projects
None yet
Milestone
No milestone
Development

Successfully merging this pull request may close these issues.
1 participant
@idkakorin