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

GENIE G4 #67

Closed
jtenavidal opened this issue Jul 24, 2023 · 9 comments
Closed

GENIE G4 #67

jtenavidal opened this issue Jul 24, 2023 · 9 comments

Comments

@jtenavidal
Copy link
Member

jtenavidal commented Jul 24, 2023 •
edited
Loading

Some occasional crash

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
@jtenavidal
Copy link
Member Author

$GENIE/data/evgen/catalogues/pdg/genie_pdg_table.txt
Oxigen 11 isotope not available

@jtenavidal
Copy link
Member Author

Problem non existent for MEC Empirial using neutrinos

@jtenavidal
Copy link
Member Author

jtenavidal commented Jul 25, 2023 •
edited
Loading

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

@jtenavidal
Copy link
Member Author

Crash output:

1690281231 NOTICE gevgen : [n] <gEvGen.cxx::GenerateEventsAtFixedInitState (336)> :  *** Generating event............ 3568
 Process event 

|------------------------------------------------------------------------------------------------------------------|
|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 |   1 | 1000060110 |   1 |  -1 |  -1 |  -1 |  -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 | 
|------------------------------------------------------------------------------------------------------------------|
|       Fin-Init:                                                |   0.121 |  -0.122 |  -0.280 |  -1.231 |         | 
|------------------------------------------------------------------------------------------------------------------|
|       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
--------------------------------------------------------------------------------------------------------------


1690281231 NOTICE HG4BertCascIntranuke : [n] <HG4BertCascIntranuke.cxx::SetTrackingRadius (356)> : Setting tracking radius to R = 9.6156
1690281231 ERROR PDG : [n] <PDGLibrary.cxx::Find (93)> : Requested missing particle with PDG: 1000080110
1690281231 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

It seems a very unusual crash... I cannot reproduce it without the exact setup above. can we simply add O11?
Produced only after DIS scattering.

@jtenavidal
Copy link
Member Author 

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

@jtenavidal
Copy link
Member Author

Wrong things:

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

https://github.com/GENIE-MC/Generator/blob/efac7f52d1975956751df658fdbf0913e309c7d9/src/Physics/HadronTransport/HG4BertCascIntranuke.cxx#L600C81-L600C89

@jtenavidal
Copy link
Member Author

https://github.com/GENIE-MC/Generator/blob/efac7f52d1975956751df658fdbf0913e309c7d9/src/Physics/HadronTransport/HG4BertCascIntranuke.cxx#L595

@jtenavidal
Copy link
Member Author

moved the discussion here: GENIE-MC/Generator#305

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
1 participant
@jtenavidal