Skip to content

Test of HINCLCascadeIntranuke and HG4BertCascIntranuke models. #359

Open
@idkakorin

Description

@idkakorin

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

Type

No type

Projects

No projects

Milestone

No milestone

Relationships

None yet

Development

No branches or pull requests

Issue actions