Open
Description
I begin to test the HINCLCascadeIntranuke model. The following event seems strange to me, because in it a proton is mother of proton, neutron, deuterium and helium.
|------------------------------------------------------------------------------------------------------------------|
|GENIE GHEP Event Record [print level: 3] |
|------------------------------------------------------------------------------------------------------------------|
| Idx | Name | Ist | PDG | Mother | Daughter | Px | Py | Pz | E | m |
|------------------------------------------------------------------------------------------------------------------|
| 0 | e- | 0 | 11 | -1 | -1 | 4 | 4 | 0.000 | 0.000 | 3.000 | 3.000 | 0.001 |
| 1 | C12 | 0 | 1000060120 | -1 | -1 | 2 | 3 | 0.000 | 0.000 | 0.000 | 11.175 | 11.175 |
| 2 | proton | 11 | 2212 | 1 | -1 | 5 | 5 | -0.064 | 0.112 | -0.102 | 0.921 | **0.938 | M = 0.906
| 3 | B11 | 2 | 1000050110 | 1 | -1 | 10 | 10 | 0.064 | -0.112 | 0.102 | 10.254 | 10.253 |
| 4 | e- | 1 | 11 | 0 | -1 | -1 | -1 | -0.172 | -0.185 | 2.920 | 2.931 | 0.001 | P = (0.059,0.063,-0.996)
| 5 | proton | 14 | 2212 | 2 | -1 | 6 | 9 | 0.108 | 0.296 | -0.022 | 0.990 | 0.938 |
| 6 | H2 | 1 | 1000010020 | 5 | -1 | -1 | -1 | 0.162 | 0.216 | -0.121 | 1.899 | 1.876 |
| 7 | proton | 1 | 2212 | 5 | -1 | -1 | -1 | -0.015 | -0.064 | -0.067 | 0.943 | 0.938 |
| 8 | neutron | 1 | 2112 | 5 | -1 | -1 | -1 | 0.017 | 0.031 | 0.057 | 0.942 | 0.940 |
| 9 | He4 | 1 | 1000020040 | 5 | -1 | -1 | -1 | 0.051 | -0.082 | 0.120 | 3.731 | 3.727 |
| 10 | He4 | 1 | 1000020040 | 3 | -1 | -1 | -1 | -0.044 | 0.083 | 0.090 | 3.730 | 3.727 |
|------------------------------------------------------------------------------------------------------------------|
| Fin-Init: | 0.000 | -0.000 | -0.000 | -0.000 | |
|------------------------------------------------------------------------------------------------------------------|
| Vertex: e- @ (x = 0.00000 m, y = 0.00000 m, z = 0.00000 m, t = 0.000000e+00 s) |
|------------------------------------------------------------------------------------------------------------------|
| Err flag [bits:15->0] : 0000000000000000 | 1st set: none |
| Err mask [bits:15->0] : 1111111111111111 | Is unphysical: NO | Accepted: YES |
|------------------------------------------------------------------------------------------------------------------|
| sig(Ev) = 3.37581e-29 cm^2 | dsig(Q2;E)/dQ2 = 4.52020e-29 cm^2/GeV^2 | Weight = 1.00000 |
|------------------------------------------------------------------------------------------------------------------|
--------------------------------------------------------------------------------------------------------------
GENIE Interaction Summary
--------------------------------------------------------------------------------------------------------------
[-] [Init-State]
|--> probe : PDG-code = 11 (e-)
|--> nucl. target : Z = 6, A = 12, PDG-Code = 1000060120 (C12)
|--> hit nucleon : PDC-Code = 2212 (proton)
|--> hit quark : no set
|--> probe 4P : (E = 3.000000, Px = 0.000000, Py = 0.000000, Pz = 3.000000)
|--> target 4P : (E = 11.174863, Px = 0.000000, Py = 0.000000, Pz = 0.000000)
|--> nucleon 4P : (E = 0.921003, Px = -0.063703, Py = 0.111874, Pz = -0.101640)
[-] [Process-Info]
|--> Interaction : EM
|--> Scattering : QES
[-] [Kinematics]
|--> *Selected* Bjorken x = 0.525450
|--> *Selected* Inelasticity y = 0.020267
|--> *Selected* Momentum transfer Q2 (>0) = 0.065342
|--> *Selected* Hadronic invariant mass W = 0.938272
[-] [Exclusive Process Info]
|--> charm prod. : false |--> strange prod. : false
|--> f/s nucleons : N(p) = 0 N(n) = 0
|--> f/s pions : N(pi^0) = 0 N(pi^+) = 0 N(pi^-) = 0
|--> f/s Other : N(gamma) = 0 N(Rho^0) = 0 N(Rho^+) = 0 N(Rho^-) = 0
|--> resonance : [not set]
|--> final quark prod. : false
|--> final lepton prod. : false
--------------------------------------------------------------------------------------------------------------
The above event has number 2 and were generated with command
gevgen -n 200 -p 11 -e 3 -t 1000060120 --tune G18_10a_00_000 --cross-sections xsec_splines_G18_10a_00_000.xml --event-generator-list EMQE --seed 65539
I set <param type="alg" name="HadronTransp-Model"> genie::HINCLCascadeIntranuke/Default </param>
in the configuration file G18_10a/ModelConfiguration.xml.
Another example of such event has number 5:
|------------------------------------------------------------------------------------------------------------------|
|GENIE GHEP Event Record [print level: 3] |
|------------------------------------------------------------------------------------------------------------------|
| Idx | Name | Ist | PDG | Mother | Daughter | Px | Py | Pz | E | m |
|------------------------------------------------------------------------------------------------------------------|
| 0 | e- | 0 | 11 | -1 | -1 | 4 | 4 | 0.000 | 0.000 | 3.000 | 3.000 | 0.001 |
| 1 | C12 | 0 | 1000060120 | -1 | -1 | 2 | 3 | 0.000 | 0.000 | 0.000 | 11.175 | 11.175 |
| 2 | proton | 11 | 2212 | 1 | -1 | 5 | 5 | -0.020 | -0.029 | 0.075 | 0.922 | **0.938 | M = 0.918
| 3 | B11 | 2 | 1000050110 | 1 | -1 | 11 | 11 | 0.020 | 0.029 | -0.075 | 10.253 | 10.253 |
| 4 | e- | 1 | 11 | 0 | -1 | -1 | -1 | 0.367 | -0.098 | 2.856 | 2.881 | 0.001 | P = (-0.128,0.034,-0.991)
| 5 | proton | 14 | 2212 | 2 | -1 | 6 | 10 | -0.388 | 0.069 | 0.219 | 1.041 | 0.938 |
| 6 | neutron | 1 | 2112 | 5 | -1 | -1 | -1 | 0.079 | 0.153 | 0.120 | 0.963 | 0.940 |
| 7 | proton | 1 | 2212 | 5 | -1 | -1 | -1 | -0.057 | -0.084 | -0.098 | 0.949 | 0.938 |
| 8 | proton | 1 | 2212 | 5 | -1 | -1 | -1 | -0.066 | 0.008 | -0.019 | 0.941 | 0.938 |
| 9 | H2 | 1 | 1000010020 | 5 | -1 | -1 | -1 | -0.019 | 0.201 | 0.076 | 1.888 | 1.876 |
| 10 | H3 | 1 | 1000010030 | 5 | -1 | -1 | -1 | -0.220 | -0.075 | 0.116 | 2.821 | 2.809 |
| 11 | He4 | 1 | 1000020040 | 3 | -1 | -1 | -1 | -0.084 | -0.105 | -0.051 | 3.730 | 3.727 |
|------------------------------------------------------------------------------------------------------------------|
| Fin-Init: | -0.000 | 0.000 | -0.000 | -0.002 | |
|------------------------------------------------------------------------------------------------------------------|
| Vertex: e- @ (x = 0.00000 m, y = 0.00000 m, z = 0.00000 m, t = 0.000000e+00 s) |
|------------------------------------------------------------------------------------------------------------------|
| Err flag [bits:15->0] : 0000000000000000 | 1st set: none |
| Err mask [bits:15->0] : 1111111111111111 | Is unphysical: NO | Accepted: YES |
|------------------------------------------------------------------------------------------------------------------|
| sig(Ev) = 3.37581e-29 cm^2 | dsig(Q2;E)/dQ2 = 6.09729e-30 cm^2/GeV^2 | Weight = 1.00000 |
|------------------------------------------------------------------------------------------------------------------|
--------------------------------------------------------------------------------------------------------------
GENIE Interaction Summary
--------------------------------------------------------------------------------------------------------------
[-] [Init-State]
|--> probe : PDG-code = 11 (e-)
|--> nucl. target : Z = 6, A = 12, PDG-Code = 1000060120 (C12)
|--> hit nucleon : PDC-Code = 2212 (proton)
|--> hit quark : no set
|--> probe 4P : (E = 3.000000, Px = 0.000000, Py = 0.000000, Pz = 3.000000)
|--> target 4P : (E = 11.174863, Px = 0.000000, Py = 0.000000, Pz = 0.000000)
|--> nucleon 4P : (E = 0.921982, Px = -0.020381, Py = -0.028634, Pz = 0.074807)
[-] [Process-Info]
|--> Interaction : EM
|--> Scattering : QES
[-] [Kinematics]
|--> *Selected* Bjorken x = 0.802749
|--> *Selected* Inelasticity y = 0.037037
|--> *Selected* Momentum transfer Q2 (>0) = 0.151128
|--> *Selected* Hadronic invariant mass W = 0.938272
[-] [Exclusive Process Info]
|--> charm prod. : false |--> strange prod. : false
|--> f/s nucleons : N(p) = 0 N(n) = 0
|--> f/s pions : N(pi^0) = 0 N(pi^+) = 0 N(pi^-) = 0
|--> f/s Other : N(gamma) = 0 N(Rho^0) = 0 N(Rho^+) = 0 N(Rho^-) = 0
|--> resonance : [not set]
|--> final quark prod. : false
|--> final lepton prod. : false
--------------------------------------------------------------------------------------------------------------
Originally posted by @idkakorin in #355 (comment)
Metadata
Metadata
Labels
No labels