Hi Steve,
I simulated 10000 nu_e (90-100 MeV) with carbon 12 target.
About 7000 are CC events, almost all of them have at least one free proton.
In them, there are 2200 events with a free neutron in the final states
produced
with the same mechanism. So the fraction is indeed high and raised my
attention.
So I am wondering what is the best (reasonable) treatment for them?
Shall I just considering the proton kinetic energy and ignore these
neutrons?
Do you have some suggestions? Thank you.
--
Cheers,
Zhe Wang
On 10/27/2013 10:43 PM, Steve Dytman wrote:
> It is the result of a final state interaction (FSI) of type 4=inelastic.
> This event shows other problems that we are working on.
>
> Steve
>
> On 10/27/2013 6:55 AM, Zhe Wang wrote:
>> Dear expert,
>>
>> I am trying to understand a snip of the generation log file.
>> I don't quite understand how the track no. 7, neutron is generated.
>> Do you know why? Thank you in advance.
>>
>> --
>> Cheers,
>> Zhe Wang
>>
>> |------------------------------------------------------------------------------------------------------------------|
>>
>> |GENIE GHEP Event Record [print level: 3] |
>> |------------------------------------------------------------------------------------------------------------------|
>>
>> | Idx | Name | Ist | PDG | Mother | Daughter
>> | Px | Py | Pz | E | m |
>> |------------------------------------------------------------------------------------------------------------------|
>>
>> | 0 | nu_e | 0 | 12 | -1 | -1 | 4 | 4 |
>> 0.000 | 0.000 | 0.100 | 0.100 | 0.000 |
>> | 1 | C12 | 0 | 1000060120 | -1 | -1 | 2 | 3 |
>> 0.000 | 0.000 | 0.000 | 11.179 | 11.179 |
>> | 2 | neutron | 11 | 2112 | 1 | -1 | 5 | 5 |
>> -0.046 | 0.148 | 0.072 | 0.934 | **0.940 | M = 0.919
>> | 3 | C11 | 2 | 1000060110 | 1 | -1 | 8 | 8 |
>> 0.046 | -0.148 | -0.072 | 10.245 | 10.243 |
>> | 4 | e- | 1 | 11 | 0 | -1 | -1 | -1 |
>> 0.025 | -0.043 | -0.023 | 0.054 | 0.001 | P = (-0.461,0.783,0.418)
>> | 5 | proton | 14 | 2212 | 2 | -1 | 6 | 7 |
>> -0.071 | 0.190 | 0.194 | 0.979 | 0.938 | FSI = 4
>> | 6 | proton | 1 | 2212 | 5 | -1 | -1 | -1 |
>> -0.000 | 0.000 | -0.000 | 0.938 | 0.938 |
>> | 7 | neutron | 1 | 2112 | 5 | -1 | -1 | -1 |
>> -0.006 | -0.003 | 0.036 | 0.940 | 0.940 |
>> | 8 | HadrBlob | 15 | 2000000002 | 3 | -1 | -1 | -1 |
>> 0.032 | 0.009 | -0.008 | 9.315 | **0.000 | M = 9.315
>> | 9 | NucBindE | 1 | 2000000101 | -1 | -1 | -1 | -1 |
>> -0.019 | 0.053 | -0.087 | 0.025 | **0.000 | M = -0.100
>> | 10 | NucBindE | 1 | 2000000101 | -1 | -1 | -1 | -1 |
>> -0.031 | -0.016 | 0.181 | 0.025 | **0.000 | M = -0.183
>> |------------------------------------------------------------------------------------------------------------------|
>>
>> | Fin-Init: | -0.000 | 0.000 | 0.000 | 0.019 | |
>> |------------------------------------------------------------------------------------------------------------------|
>>
>> | Vertex: nu_e @ (x = 0.00000 m, y = 0.00000 m, z
>> = 0.00000 m, t = 0 s) |
>> |------------------------------------------------------------------------------------------------------------------|
>>
>> | Err flag [bits:15->0] : 0000000000000000 | 1st
>> set: none |
>> | Err mask [bits:15->0] : 1111111111111111 | Is unphysical: NO
>> | Accepted: YES |
>> |------------------------------------------------------------------------------------------------------------------|
>>
>> | sig(Ev) = 1.7352e-39 cm^2 | dsig(Q2;E)/dQ2 = 2.8019e-38
>> cm^2/GeV^2 | Weight = 1 |
>> |------------------------------------------------------------------------------------------------------------------|
>>
>>
>>
>> --------------------------------------------------------------------------------------------------------------
>>
>> GENIE Interaction Summary
>> --------------------------------------------------------------------------------------------------------------
>>
>> [-] [Init-State]
>> |--> probe : PDG-code = 12 (nu_e)
>> |--> nucl. target : Z = 6, A = 12, PDG-Code = 1000060120 (C12)
>> |--> hit nucleon : PDC-Code = 2112 (neutron)
>> |--> hit quark : no set
>> |--> probe 4P : (E = 0.0995661, Px = 0, Py
>> = 0, Pz = 0.0995661)
>> |--> target 4P : (E = 11.179, Px = 0, Py
>> = 0, Pz = 0)
>> |--> nucleon 4P : (E = 0.9342, Px = -0.0456271, Py =
>> 0.147645, Pz = 0.0720215)
>> [-] [Process-Info]
>> |--> Interaction : Weak[CC]
>> |--> Scattering : QES
>> [-] [Kinematics]
>> |--> *Selected* Bjorken x = 0.294936
>> |--> *Selected* Inelasticity y = 0.303082
>> |--> *Selected* Momentum transfer Q2 (>0) = 0.0153471
>> |--> *Selected* Hadronic invariant mass W = 0.93827
>> [-] [Exclusive Process Info]
>> |--> charm prod. : false
>> |--> f/s nucleons : N(p) = 0 N(n) = 0
>> |--> f/s pions : N(pi^0) = 0 N(pi^+) = 0 N(pi^-) = 0
>> |--> resonance : [not set]
>> --------------------------------------------------------------------------------------------------------------
>>
|