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

Error in G4 and INCL++ cascades due to wrong physics #305

Open jtenavidal opened 1 year ago

jtenavidal commented 1 year ago

I am running e-+C12 events using G18_10d (G4 FSI). The event generation crashes due to an unexpected isotope in the final state. The crash observed is the following:

1690205249 NOTICE gevgen : [n] <gEvGen.cxx::GenerateEventsAtFixedInitState (336)> :  *** Generating event............ 3567
1690205249 NOTICE HG4BertCascIntranuke : [n] <HG4BertCascIntranuke.cxx::SetTrackingRadius (354)> : Setting tracking radius to R = 9.6156
1690205249 NOTICE gevgen : [n] <gEvGen.cxx::GenerateEventsAtFixedInitState (348)> : Generated Event GHEP Record: 

|------------------------------------------------------------------------------------------------------------------|
|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 |   2.261 |   2.261 |   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.071 |   0.041 |  -0.070 |   0.922 | **0.938 | M = 0.915 
|   3 |           B11 |   2 | 1000050110 |   1 |  -1 |  10 |  12 |  -0.071 |  -0.041 |   0.070 |  10.253 |  10.253 | 
|   4 |            e- |   1 |         11 |   0 |  -1 |  -1 |  -1 |  -0.215 |  -0.051 |   2.187 |   2.198 |   0.001 | P = (0.098,0.023,-0.995)
|   5 |        proton |  14 |       2212 |   2 |  -1 |  -1 |  -1 |   0.286 |   0.092 |   0.005 |   0.985 |   0.938 | 
|   6 |        proton |   1 |       2212 |  -1 |  -1 |  -1 |  -1 |   0.105 |   0.050 |  -0.016 |   0.946 |   0.938 | 
|   7 |       neutron |   1 |       2112 |  -1 |  -1 |  -1 |  -1 |   0.079 |   0.084 |   0.056 |   0.948 |   0.940 | 
|   8 |         gamma |   1 |         22 |  -1 |  -1 |  -1 |  -1 |   0.000 |   0.001 |  -0.000 |   0.001 |   0.000 | 
|   9 |         gamma |   1 |         22 |  -1 |  -1 |  -1 |  -1 |  -0.000 |  -0.000 |   0.000 |   0.000 |   0.000 | 
|  10 |            H2 |   1 | 1000010020 |   3 |   0 |  -1 |  -1 |   0.004 |   0.124 |   0.072 |   1.881 |   1.876 | 
|  11 |           He4 |   1 | 1000020040 |   3 |   0 |  -1 |  -1 |  -0.029 |  -0.015 |  -0.021 |   3.728 |   3.727 | 
|  12 |           He4 |  15 | 1000020040 |   3 |  -1 |  -1 |  -1 |   0.056 |  -0.192 |  -0.017 |   3.734 |   3.727 | 
|------------------------------------------------------------------------------------------------------------------|
|       Fin-Init:                                                |  -0.000 |  -0.000 |  -0.000 |  -0.001 |         | 
|------------------------------------------------------------------------------------------------------------------|
|       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.35944e-29 cm^2  | dsig(Q2;E)/dQ2 =          8.07790e-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 =     2.261000, Px =     0.000000, Py =     0.000000, Pz =     2.261000)
 |--> target 4P    : (E =    11.174863, Px =     0.000000, Py =     0.000000, Pz =     0.000000)
 |--> nucleon 4P   : (E =     0.921755, Px =     0.070841, Py =     0.040529, Pz =    -0.069544)
[-] [Process-Info]  
 |--> Interaction : EM
 |--> Scattering  : QES
[-] [Kinematics]
 |--> *Selected* Bjorken x = 0.543826
 |--> *Selected* Inelasticity y = 0.020647
 |--> *Selected* Momentum transfer Q2 (>0) = 0.050333
 |--> *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
--------------------------------------------------------------------------------------------------------------

1690205249 NOTICE gevgen : [n] <gEvGen.cxx::GenerateEventsAtFixedInitState (336)> :  *** Generating event............ 3568
1690205249 NOTICE HG4BertCascIntranuke : [n] <HG4BertCascIntranuke.cxx::SetTrackingRadius (354)> : Setting tracking radius to R = 9.6156
1690205249 ERROR PDG : [n] <PDGLibrary.cxx::Find (93)> : Requested missing particle with PDG: 1000080110
1690205249 FATAL GHepParticle : [n] <GHepParticle.cxx::AssertIsKnownParticle (539)> : 
** You are attempting to insert particle with PDG code = 1000080110 into the event record.
** This particle can not be found in $GENIE/data/evgen/catalogues/pdg/genie_pdg_table.txt

The event 3568 is not successfully generated due to this error: You are attempting to insert particle with PDG code = 1000080110 into the event record.. As a consequence, you cannot see the event record (notice the one you can see is from the previous event).

You can recover the crash with the following command:

gevgen -n 10000 -p 11 -e 2.261 -t 1000060120 --tune G18_10d_00_000 --cross-sections G18_10d_nuclei_Q2_def_total_xsec.xml --event-generator-list EM --seed 65539

The trick we followed to debug this further is to add the Oxigen 11 isotope in the corresponding GENIE file: $GENIE/data/evgen/catalogues/pdg/genie_pdg_table.txt. In that way we can recover the event record for this specific event:

1690283540 NOTICE gevgen : [n] <gEvGen.cxx::GenerateEventsAtFixedInitState (348)> : Generated Event GHEP Record: 

|------------------------------------------------------------------------------------------------------------------|
|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 |   2.261 |   2.261 |   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.016 |   0.112 |  -0.036 |   0.920 | **0.940 | M = 0.913 
|   3 |           C11 |   2 | 1000060110 |   1 |  -1 |  10 |  10 |  -0.016 |  -0.112 |   0.036 |  10.255 |  10.254 | 
|   4 |            e- |   1 |         11 |   0 |  -1 |  -1 |  -1 |   0.138 |  -0.011 |   1.945 |   1.950 |   0.001 | P = (-0.071,0.006,-0.997)
|   5 |      HadrSyst |  12 | 2000000001 |   2 |  -1 |   6 |   7 |  -0.121 |   0.122 |   0.280 |   1.231 | **0.000 | M = 1.186 
|   6 |        proton |  14 |       2212 |   5 |  -1 |  -1 |  -1 |  -0.139 |   0.200 |   0.077 |   0.972 |   0.938 | 
|   7 |           pi- |  14 |       -211 |   5 |  -1 |  -1 |  -1 |   0.017 |  -0.077 |   0.203 |   0.259 |   0.140 | 
|   8 |       neutron |   1 |       2112 |  -1 |  -1 |  -1 |  -1 |   0.076 |  -0.195 |   0.210 |   0.985 |   0.940 | 
|   9 |           pi- |   1 |       -211 |  -1 |  -1 |  -1 |  -1 |   0.047 |   0.094 |  -0.049 |   0.181 |   0.140 | 
|  10 |           O11 |  15 | 1000080110 |   3 |  -1 |  -1 |  -1 |  -0.261 |   0.112 |   0.155 |  10.317 | *11.206 | M = 10.312 
|------------------------------------------------------------------------------------------------------------------|
|       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) =       2.43153e-30 cm^2  | d2sig(x,y;E)/dxdy =       7.19287e-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 =     2.261000, Px =     0.000000, Py =     0.000000, Pz =     2.261000)
 |--> target 4P    : (E =    11.174863, Px =     0.000000, Py =     0.000000, Pz =     0.000000)
 |--> nucleon 4P   : (E =     0.920162, Px =     0.016454, Py =     0.111531, Pz =    -0.035750)
[-] [Process-Info]  
 |--> Interaction : EM
 |--> Scattering  : DIS
[-] [Kinematics]
 |-->  *Running* Hadronic invariant mass W = 1.186371
 |--> *Selected* Bjorken x = 0.037070
 |--> *Selected* Inelasticity y = 0.138086
 |--> *Selected* Momentum transfer Q2 (>0) = 0.022127
 |--> *Selected* Hadronic invariant mass W = 1.186371
[-] [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
--------------------------------------------------------------------------------------------------------------

Issues with the event:

jtenavidal commented 1 year ago

INCL++ shows the same behaviour.

Output before adding O11:

1690379613 NOTICE gevgen : [n] <gEvGen.cxx::GenerateEventsAtFixedInitState (336)> :  *** Generating event............ 216
1690379613 NOTICE PauliBlock : [n] <PauliBlocker.cxx::ProcessEventRecord (106)> :  *** The generated event is Pauli-blocked (|p_{nucleon}| = 0.156324 GeV < Fermi momentum = 0.192431 GeV) ***
1690379613 WARN EventGenerator : [n] <EventGenerator.cxx::ProcessEventRecord (142)> : An unphysical event was generated and was rejected
1690379613 NOTICE HINCLCascadeIntranuke : [n] <HINCLCascadeIntranuke.cxx::ProcessEventRecord (459)> : ************ Running HINCLCascadeIntranuke MODE INTRANUKE ************
1690379613 NOTICE HINCLCascadeIntranuke : [n] <HINCLCascadeIntranuke.cxx::TransportHadrons (500)> : Propagating hadrons within nucleus found in position = 3
1690379613 NOTICE HINCLCascadeIntranuke : [n] <HINCLCascadeIntranuke.cxx::TransportHadrons (516)> : Nucleus (A,Z) = (11, 5)
1690379613 NOTICE gevgen : [n] <gEvGen.cxx::GenerateEventsAtFixedInitState (348)> : Generated Event GHEP Record: 

|------------------------------------------------------------------------------------------------------------------|
|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 |   2.261 |   2.261 |   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.006 |   0.170 |  -0.085 |   0.921 | **0.938 | M = 0.901 
|   3 |           B11 |   2 | 1000050110 |   1 |  -1 |  -1 |  -1 |  -0.006 |  -0.170 |   0.085 |  10.254 |  10.253 | 
|   4 |            e- |   1 |         11 |   0 |  -1 |  -1 |  -1 |   0.243 |  -0.068 |   2.170 |   2.185 |   0.001 | P = (-0.111,0.031,-0.993)
|   5 |        proton |  14 |       2212 |   2 |  -1 |   6 |   9 |  -0.237 |   0.238 |   0.006 |   0.996 |   0.938 | 
|   6 |       neutron |   1 |       2112 |   5 |   3 |  -1 |  -1 |  -0.180 |   0.109 |   0.138 |   0.973 |   0.940 | 
|   7 |        proton |   1 |       2212 |   5 |   3 |  -1 |  -1 |  -0.051 |   0.057 |   0.007 |   0.941 |   0.938 | 
|   8 |           Li6 |   1 | 1000030060 |   5 |   3 |  -1 |  -1 |  -0.098 |  -0.148 |  -0.001 |   5.604 |   5.602 | 
|   9 |           He4 |   1 | 1000020040 |   5 |   3 |  -1 |  -1 |   0.086 |   0.050 |  -0.053 |   3.729 |   3.727 | 
|------------------------------------------------------------------------------------------------------------------|
|       Fin-Init:                                                |  -0.000 |   0.000 |   0.000 |  -0.003 |         | 
|------------------------------------------------------------------------------------------------------------------|
|       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.89137e-31 cm^2  | dsig(Q2;E)/dQ2 =          4.37022e-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 =     2.261000, Px =     0.000000, Py =     0.000000, Pz =     2.261000)
 |--> target 4P    : (E =    11.174863, Px =     0.000000, Py =     0.000000, Pz =     0.000000)
 |--> nucleon 4P   : (E =     0.920562, Px =     0.006256, Py =     0.169538, Pz =    -0.084653)
[-] [Process-Info]  
 |--> Interaction : EM
 |--> Scattering  : QES
[-] [Kinematics]
 |--> *Selected* Bjorken x = 0.490088
 |--> *Selected* Inelasticity y = 0.029712
 |--> *Selected* Momentum transfer Q2 (>0) = 0.066191
 |--> *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
--------------------------------------------------------------------------------------------------------------

1690379613 NOTICE gevgen : [n] <gEvGen.cxx::GenerateEventsAtFixedInitState (336)> :  *** Generating event............ 217
1690379613 NOTICE HINCLCascadeIntranuke : [n] <HINCLCascadeIntranuke.cxx::ProcessEventRecord (459)> : ************ Running HINCLCascadeIntranuke MODE INTRANUKE ************
1690379613 NOTICE HINCLCascadeIntranuke : [n] <HINCLCascadeIntranuke.cxx::TransportHadrons (500)> : Propagating hadrons within nucleus found in position = 3
1690379613 NOTICE HINCLCascadeIntranuke : [n] <HINCLCascadeIntranuke.cxx::TransportHadrons (516)> : Nucleus (A,Z) = (11, 6)
1690379613 ERROR PDG : [n] <PDGLibrary.cxx::Find (93)> : Requested missing particle with PDG: 1000080110
Segmentation fault

after adding oxigen 11:

1690380108 NOTICE gevgen : [n] <gEvGen.cxx::GenerateEventsAtFixedInitState (336)> :  *** Generating event............ 217
1690380108 NOTICE HINCLCascadeIntranuke : [n] <HINCLCascadeIntranuke.cxx::ProcessEventRecord (459)> : ************ Running HINCLCascadeIntranuke MODE INTRANUKE ************
1690380108 NOTICE HINCLCascadeIntranuke : [n] <HINCLCascadeIntranuke.cxx::TransportHadrons (500)> : Propagating hadrons within nucleus found in position = 3
1690380108 NOTICE HINCLCascadeIntranuke : [n] <HINCLCascadeIntranuke.cxx::TransportHadrons (516)> : Nucleus (A,Z) = (11, 6)
1690380108 NOTICE gevgen : [n] <gEvGen.cxx::GenerateEventsAtFixedInitState (348)> : Generated Event GHEP Record: 

|------------------------------------------------------------------------------------------------------------------|
|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 |   2.261 |   2.261 |   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.019 |  -0.164 |  -0.135 |   0.919 | **0.940 | M = 0.894 
|   3 |           C11 |   2 | 1000060110 |   1 |  -1 |  -1 |  -1 |   0.019 |   0.164 |   0.135 |  10.256 |  10.254 | 
|   4 |            e- |   1 |         11 |   0 |  -1 |  -1 |  -1 |   0.456 |   0.131 |   1.723 |   1.788 |   0.001 | P = (-0.255,-0.073,-0.964)
|   5 |        Delta0 |   3 |       2114 |   2 |  -1 |   6 |   7 |  -0.475 |  -0.294 |   0.402 |   1.392 | **1.233 | M = 1.210 
|   6 |       neutron |  14 |       2112 |   5 |  -1 |   8 |   8 |  -0.592 |  -0.289 |   0.301 |   1.186 |   0.940 | 
|   7 |           pi0 |  14 |        111 |   5 |  -1 |   9 |  10 |   0.117 |  -0.005 |   0.102 |   0.206 |   0.135 | 
|   8 |       neutron |   1 |       2112 |   6 |  -1 |  -1 |  -1 |  -0.592 |  -0.289 |   0.301 |   1.186 |   0.940 | 
|   9 |           pi- |   1 |       -211 |   7 |   3 |  -1 |  -1 |   0.096 |   0.040 |   0.051 |   0.181 |   0.140 | 
|  10 |           O11 |   1 | 1000080110 |   7 |   3 |  -1 |  -1 |   0.040 |   0.118 |   0.186 |  11.208 |  11.206 | 
|------------------------------------------------------------------------------------------------------------------|
|       Fin-Init:                                                |  -0.000 |  -0.000 |  -0.000 |   0.928 |         | 
|------------------------------------------------------------------------------------------------------------------|
|       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.47523e-31 cm^2  | d2sig(W,Q2;E)/dWdQ2 =     1.84473e-29 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 =     2.261000, Px =     0.000000, Py =     0.000000, Pz =     2.261000)
 |--> target 4P    : (E =    11.174863, Px =     0.000000, Py =     0.000000, Pz =     0.000000)
 |--> nucleon 4P   : (E =     0.918629, Px =    -0.018918, Py =    -0.163703, Pz =    -0.135161)
[-] [Process-Info]  
 |--> Interaction : EM
 |--> Scattering  : RES
[-] [Kinematics]
 |--> *Selected* Bjorken x = 0.304039
 |--> *Selected* Inelasticity y = 0.200414
 |--> *Selected* Momentum transfer Q2 (>0) = 0.290363
 |--> *Selected* Hadronic invariant mass W = 1.209591
[-] [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
--------------------------------------------------------------------------------------------------------------
dytman commented 1 year ago

This is exotic event and exotic final state. Is it just luck? see https://www.nndc.bnl.gov/nudat3/. O12 has very short halflife. O11 is proposed but not seen. It would be 8p and 3n, very exotic. I think you did right thing. We are using other people's code in new situation. Do we need an error condition for such an event?

jtenavidal commented 1 year ago

Hi Steve, as mentioned above the event is wrong because O11 should not be there. Also notice that the same error is obtained using both INCL and G4 (different codes). I am investigating this further for now, I want to understand how often these kind of events occur. For now I have some questions unsolved:

jtenavidal commented 1 year ago

INCL:

idkakorin commented 1 year ago

Hi Julia,

My opinion is that if we add some isotopes in _genie_pdgtable.txt we should do it in consistent way, namely, check whether that happen with other nuclei, and if yes, add corresponding isotopes for them too.

mroda88 commented 1 year ago

Hi all, I'm chipping in as I think things are getting confused. The problem is not that O11 is created. We know it is a valid isotope that it is documented in nuclide tables. And of course if we think we need, we can add it to the table. Indeed it is what she did in order to generate a complete file that we can analise.

The problem is that in events like the one Julia reported, simply for pure charge and baryon conservation, that Oxigen should be a Nitrogen. The missing PDG error simply makes the problem visible as it generates an error, but the problem is in the physics. The error is just a weird symptom. Obviously the concern here is that problems like these can happen in other situations which we haven't noticed because the final nucleus is available and no errors are thrown.

I hope to have a discussion about this on the next Core developer meeting (14 August) assuming the necessary players will be able to connect, despite summer and possible vacations. If not, please send some notes, either here or with an email.

dytman commented 1 year ago

even 196 out of 400k still seems like a lot to me.

Important question - does INCL or GENIE calculate final state? As you say, event printout is apparently wrong. I see final nucleus is C12 missing a neutron which is C11. Nothing there says it should be O11.

dytman commented 1 year ago

According to data, O11 has not been seen! Surely, someone could calculate its mass so we can include in our tables if needed.

jtenavidal commented 1 year ago

Richie Diurba also recovered the same error running gevgen_hadron with a pion on lead 10000822080 with INCL++ and Geant4. 1690877562 ERROR PDG : [n] <PDGLibrary.cxx::Find (93)> : Requested missing particle with PDG: 1000751990

rdiurba commented 1 year ago

I was generating events for studies to tune hA2018 and develop systematics for it. This was a replication of the studies documented here: https://github.com/GENIE-MC/Generator/blob/master/src/Physics/HadronTransport/HAIntranuke2018.cxx#L739. I was generating HINCL and hG4BertCasc events for piplus and piminus. I then would get errors for carbon and iron. The carbon error was for the same isotope above (O11) and this command gives me an error for Rhenium-199:

gevgen_hadron -n 1000000 -p -211 -t 1000822080 -k 0.5 -m HINCL -o /dune/data/users/rdiurba/genieFSIStudies/piMinus_1000822080_0.5GeV_HINCL_500k_0

mroda88 commented 1 year ago

@rdiurba can you to find a way to put the whole event on here? I think it's clear that the problem is not the missing nucleus in our table (that an always be added if necessary), but rather the wrong accounting that the external code does when more than one particle is propagated.

rdiurba commented 1 year ago

Unfortunately, it always exits ungracefully and spits out an empty ghep tree. Here is the printout though. Let me know if you meant something different. 

 1690970762 NOTICE gevgen_hadron : [n] <gEvGenHadronNucleus.cxx::main (200)> :  *** Generating event............ 658778
1690970762 NOTICE HINCLCascadeIntranuke : [n] <HINCLCascadeIntranuke.cxx::ProcessEventRecord (459)> : ************ Running HINCLCascadeIntranuke MODE INTRANUKE ************
1690970762 NOTICE gevgen_hadron : [n] <gEvGenHadronNucleus.cxx::main (210)> : 

|------------------------------------------------------------------------------------------------------------------|
|GENIE GHEP Event Record [print level:   3]                                                                        |
|------------------------------------------------------------------------------------------------------------------|
| Idx |          Name | Ist |        PDG |   Mother  | Daughter  |      Px |      Py |      Pz |       E |      m  | 
|------------------------------------------------------------------------------------------------------------------|
|   0 |           pi- |   0 |       -211 |  -1 |  -1 |   2 |   2 |   0.000 |   0.000 |   0.624 |   0.640 |   0.140 | 
|   1 |         Pb208 |   0 | 1000822080 |  -1 |  -1 |   3 |   3 |   0.000 |   0.000 |   0.000 | 193.687 | 193.687 | 
|   2 |           pi- |   1 |       -211 |   0 |  -1 |  -1 |  -1 |   0.000 |   0.000 |   0.624 |   0.640 |   0.140 | 
|   3 |         Pb208 |  15 | 1000822080 |   1 |  -1 |  -1 |  -1 |   0.000 |   0.000 |   0.000 | 193.687 | 193.687 | 
|------------------------------------------------------------------------------------------------------------------|
|       Fin-Init:                                                |   0.000 |   0.000 |   0.000 |  -0.000 |         | 
|------------------------------------------------------------------------------------------------------------------|
|       Vertex:            pi- @ (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) =       0.00000e+00 cm^2  | dsig(Ev;{K_s})/dK   =     0.00000e+00 cm^2/{K}   | Weight =          1.00000 |
|------------------------------------------------------------------------------------------------------------------|

--------------------------------------------------------------------------------------------------------------
GENIE Interaction Summary
--------------------------------------------------------------------------------------------------------------
[-] [Init-State] 
 |--> probe        : PDG-code = 0 (Rootino)
 |--> nucl. target : Z = 0, A = 0, PDG-Code = 0 (Rootino)
 |--> hit nucleon  : no set
 |--> hit quark    : no set
 |--> probe 4P     : (E =     0.000000, Px =     0.000000, Py =     0.000000, Pz =     0.000000)
 |--> target 4P    : (E =     0.000000, Px =     0.000000, Py =     0.000000, Pz =     0.000000)

[-] [Process-Info]  
 |--> Interaction : Unknown
 |--> Scattering  : Unknown
[-] [Kinematics]
[-] [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
--------------------------------------------------------------------------------------------------------------

1690970762 NOTICE gevgen_hadron : [n] <gEvGenHadronNucleus.cxx::main (200)> :  *** Generating event............ 658779
1690970762 NOTICE HINCLCascadeIntranuke : [n] <HINCLCascadeIntranuke.cxx::ProcessEventRecord (459)> : ************ Running HINCLCascadeIntranuke MODE INTRANUKE ************
1690970763 ERROR PDG : [n] <PDGLibrary.cxx::Find (93)> : Requested missing particle with PDG: 1000751990
mroda88 commented 1 year ago

in order to have the full print, Julia locally added the missing isotope in order for the event and the file to be completed. Then, she printed the full event using gevdump.

rdiurba commented 1 year ago

Ok, thanks. I am running this as a dunegpvm ups product and have never actually built with Geant4 or INCL++. If you can give me a couple days, I will try to do the same.

jtenavidal commented 1 year ago

You can use this script to build GENIE locally using G4: https://github.com/GENIE-MC/Generator/blob/master/src/scripts/production/python/setup_GENIE.sh

Simply do: source https://github.com/GENIE-MC/Generator.git master true false; for INCL build source https://github.com/GENIE-MC/Generator.git master false true; for G4 build

from the gpvms

rdiurba commented 1 year ago

Using: gevgen_hadron -n 1000000 -p -211 -t 1000822080 -k 0.5 -m HINCL

1691007074 NOTICE gevdump : [n] <gEvDump.cxx::main (177)> :  ** Event: 658779

|------------------------------------------------------------------------------------------------------------------|
|GENIE GHEP Event Record [print level:   3]                                                                        |
|------------------------------------------------------------------------------------------------------------------|
| Idx |          Name | Ist |        PDG |   Mother  | Daughter  |      Px |      Py |      Pz |       E |      m  | 
|------------------------------------------------------------------------------------------------------------------|
|   0 |           pi- |   0 |       -211 |  -1 |  -1 |   2 |  11 |   0.000 |   0.000 |   0.624 |   0.640 |   0.140 | 
|   1 |         Pb208 |   0 | 1000822080 |  -1 |  -1 |  12 |  12 |   0.000 |   0.000 |   0.000 | 193.687 | 193.687 | 
|   2 |       neutron |   1 |       2112 |   0 |  -1 |  -1 |  -1 |  -0.509 |  -0.165 |   0.297 |   1.121 |   0.940 | 
|   3 |        proton |   1 |       2212 |   0 |  -1 |  -1 |  -1 |  -0.041 |  -0.114 |   0.289 |   0.989 |   0.938 | 
|   4 |        proton |   1 |       2212 |   0 |  -1 |  -1 |  -1 |   0.120 |  -0.162 |   0.101 |   0.965 |   0.938 | 
|   5 |           pi0 |   1 |        111 |   0 |  -1 |  -1 |  -1 |   0.008 |  -0.010 |  -0.018 |   0.137 |   0.135 | 
|   6 |        proton |   1 |       2212 |   0 |  -1 |  -1 |  -1 |   0.142 |   0.137 |   0.192 |   0.978 |   0.938 | 
|   7 |        proton |   1 |       2212 |   0 |  -1 |  -1 |  -1 |   0.344 |  -0.049 |   0.020 |   1.001 |   0.938 | 
|   8 |        proton |   1 |       2212 |   0 |  -1 |  -1 |  -1 |   0.178 |   0.177 |   0.002 |   0.971 |   0.938 | 
|   9 |            H2 |   1 | 1000010020 |   0 |  -1 |  -1 |  -1 |   0.347 |   0.137 |  -0.123 |   1.916 |   1.876 | 
|  10 |       neutron |   1 |       2112 |   0 |  -1 |  -1 |  -1 |   0.034 |   0.006 |  -0.050 |   0.942 |   0.940 | 
|  11 |         gamma |   1 |         22 |   0 |  -1 |  -1 |  -1 |   0.000 |   0.001 |   0.002 |   0.002 |   0.000 | 
|  12 |         Re199 |  15 | 1000751990 |   1 |  -1 |  -1 |  -1 |  -0.623 |   0.042 |  -0.088 | 185.270 | 185.269 | 
|------------------------------------------------------------------------------------------------------------------|
|       Fin-Init:                                                |   0.000 |   0.000 |  -0.000 |  -0.035 |         | 
|------------------------------------------------------------------------------------------------------------------|
|       Vertex:            pi- @ (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) =       0.00000e+00 cm^2  | dsig(Ev;{K_s})/dK   =     0.00000e+00 cm^2/{K}   | Weight =          1.00000 |
|------------------------------------------------------------------------------------------------------------------|

--------------------------------------------------------------------------------------------------------------
GENIE Interaction Summary
--------------------------------------------------------------------------------------------------------------
[-] [Init-State] 
 |--> probe        : PDG-code = 0 (Rootino)
 |--> nucl. target : Z = 0, A = 0, PDG-Code = 0 (Rootino)
 |--> hit nucleon  : no set
 |--> hit quark    : no set
 |--> probe 4P     : (E =     0.000000, Px =     0.000000, Py =     0.000000, Pz =     0.000000)
 |--> target 4P    : (E =     0.000000, Px =     0.000000, Py =     0.000000, Pz =     0.000000)

[-] [Process-Info]  
 |--> Interaction : Unknown
 |--> Scattering  : Unknown
[-] [Kinematics]
[-] [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
--------------------------------------------------------------------------------------------------------------
mroda88 commented 10 months ago

THe fix is in #360