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

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

Open
jtenavidal opened this issue Jul 25, 2023 · 17 comments
Open

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

jtenavidal opened this issue Jul 25, 2023 · 17 comments

Comments

@jtenavidal
Copy link
Contributor

jtenavidal commented Jul 25, 2023 •
edited
Loading

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:

  • Missing links between daughters and mothers pre and post FSI
  • Assuming that the interaction is e-+C12->e- + pi- + n, the remanent nucleus should be nitrogen. We suspect that oxygen is generated assuming that two pi- are created due to charge conservation and the missing links between the particles.
  • Usually, in GENIE, the remaining nucleus is a generic HadronicBlob. Nitrogen should not even show in the event record.
@jtenavidal jtenavidal mentioned this issue Jul 25, 2023
Closed
@jtenavidal jtenavidal changed the title G4 error due to wrong physics Error in G4 and INCL++ cascades due to wrong physics Jul 26, 2023
@jtenavidal
Copy link
Contributor Author

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

dytman commented Jul 26, 2023

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

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:

  • Why is the code adding O11? Where in the code does this happen? Is it a bug in the interface or the G4/INCL codes somehow behave like this?
  • Why do we have the exact same problem for G4 and INCL?
  • Why do we not keep the mother/daughter relation for either of these codes? This information is important for users.

@jtenavidal
Copy link
Contributor Author

INCL:

  • e-C12 at 2 GeV
  • G18_10c_00_000
  • 400k events generated
  • 196 of then contain O11

@idkakorin
Copy link
Contributor

Hi Julia,

My opinion is that if we add some isotopes in genie_pdg_table.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
Copy link
Member

mroda88 commented Jul 27, 2023

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

dytman commented Jul 27, 2023

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

dytman commented Jul 27, 2023

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

@idkakorin idkakorin mentioned this issue Jul 29, 2023
Closed
@jtenavidal
Copy link
Contributor Author

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

rdiurba commented Aug 2, 2023

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
Copy link
Member

mroda88 commented Aug 2, 2023

@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
Copy link

rdiurba commented Aug 2, 2023 •
edited
Loading

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
Copy link
Member

mroda88 commented Aug 2, 2023

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

rdiurba commented Aug 2, 2023

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

jtenavidal commented Aug 2, 2023 •
edited
Loading

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

rdiurba commented Aug 3, 2023

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
Copy link
Member

mroda88 commented Nov 9, 2023

THe fix is in #360

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Assignees
No one assigned
Labels
None yet
Projects
None yet
Milestone
No milestone
Development

No branches or pull requests
5 participants
@dytman @idkakorin @mroda88 @jtenavidal @rdiurba