From inclusive to semi-inclusive one-nucleon knockout in neutrino event generators

: In neutrino event generators, for models for neutrino and electron scattering only inclusive


Introduction
In recent years, accelerator-based neutrino experiments have performed measurements of the hadronic final-state in charged-current interactions.The hadron information is useful to distinguish between different interaction mechanisms and probe part of the nuclear momentum distribution through, for example, transverse kinematic imbalance [1,2].Additionally, a precise determination of outgoing nucleon kinematics, with suitable kinematic cuts leads to a more precise reconstruction of neutrino energy on an event-by-event basis [3].
The main challenge in accelerator-based neutrino experiments is that the incoming energy distribution is broad.This means that, to describe a A(ν µ , µp)X signal one has to account for a wealth of interaction mechanisms.At sufficiently high energy and momentum transfers, we could assume that the energy-momentum is absorbed by a single nucleon which reinteracts with the nucleus.These strong interactions are dubbed final-state interactions (FSI).Within this picture, inelastic FSI, where the energy gets distributed over many different final-states, is particularly important.The large phase-space covered in experiments, and the fact that the total energy of the residual system X is not strongly constrained, make a quantum-mechanical description of all the possible coupled final-state channels intractable.For this reason experiments deal with this problem using intranuclear cascade models (INC), which provide an explicit, albeit (semi-)classical, description of rescattering [4][5][6][7][8][9][10].
For the following we consider the contribution of quasielastic interactions, where the exchanged boson is absorbed on a single nucleon that is excited to the continuum, to a 1lepton 1-proton final state, e.g. to A(ν µ , µp)X or A(e, e p)X.Within the GENIE event generator the process is described by introducing this nucleon into the INC, which redistributes the strength over final-states in a unitary way.This means that after integration over all final-states the original inclusive cross section is recovered.A number of models which describe the inclusive cross sections in the vicinity of the quasielastic peak, have been implemented in GENIE [11][12][13][14][15].The caveat is that only inclusive cross sections are provided, the full kinematic and dynamic structure of the semi-inclusive cross section [16] is not available.Therefore GENIE provides a procedure to attach the outgoing nucleon kinematics given some nuclear momentum distribution, in Ref. [12] discrepancies between this method and a microscopic calculation were discussed for flux-averaged neutrino cross sections.
In the following we present the algorithm used in GENIE to compute nucleon variables.We compare this result to fully unfactorized calculations with the relativistic distorted-wave impulse approximation (RDWIA).We focus on the (e, e p) process for kinematics relevant to the e4ν analysis [17].We show that the approximate procedure used in GENIE leads to significantly different observables than the RDWIA result.

Results
We use an unfactorized RDWIA calculation where the final-state nucleon is described in the real Energy-Dependent Relativistic Mean Field (ED-RMF) potential [18].The ED-RMF model yields an inclusive (e, e ) cross section that is practically the same as the SuSAv2 model for sufficient momentum transfer [19].The main power of this approach is that it provides an exclusive cross section, unlike the SuSAv2 approach.In Ref. [20] the RDWIA with a real potential was used as input to the INC in the NEUT generator [6].The results compare well with the T2K data of Ref. [21].It was moreover found that, at sufficiently high nucleon energy (T N 100MeV), the resulting exclusive cross sections agree with optical potential calculations.
We have generated events for the process e + A → e + p + B with E e = 1.159GeV.The events are distributed according to the (Q 2 ) 2 -weighted cross section obtained in the RDWIA (for details see Refs.[3,20]) The (Q 2 ) 2 weighting, which was introduced in the e4ν analysis of Ref. [17], makes the kinematic dependence of the weighted cross section similar to the neutrino scattering case.The sum is over the invariant masses of the residual system, given by the RMF model as in Ref. [20]; we only include scattering with protons.We impose cuts for the outgoing electron 40 • > θ e > 17 •  and E e > 400 MeV, but include the full nucleon phase space.
To test the approximate treatment, we replace the nucleon variables for every event by the ones produced by the algorithm used in GENIE, which is described in the following.As the momentum distribution we use the local Fermi gas (FG) obtained from the nuclear density, also taken from GENIE, with a = 1.69 fm and α = 1.08 for carbon.Given a vector p m sampled from the momentum distribution, the outgoing nucleon energy is determined as Here E b (q) is a q-dependent binding energy inspired by the energy shift from the SuSAv2 model and E b (q) = E b (q = 827 MeV) for q > 827 MeV.The resulting distributions of the nucleon kinetic energy are shown in the left panel of Fig. 1.The procedure is found to lead to an overall shift of approximately 20 MeV of the whole distribution.This could be amended by including an additional energy shift, which GENIE allows for.On the other hand the implementation of the on-shell dispersion relation in this expression comes from the FG approach and should be considered unrealistic [22].The angular distributions (right panel of Fig. 1) exhibit significant shape differences, partly due to the more restricted missing momentum distribution.An approximation is made when determining the nucleon angle that also contributes to the narrower angular distribution.It is clear that the magnitude of the nucleon momentum generated with Eq. 3 will not agree with To impose momentum conservation, the magnitude of the nucleon momentum is taken from Eq. 3, while the direction is taken from the momentum vectors, i.e.
The residual momentum is given to the remnant nucleus.
One might expect that these discrepancies are smeared out by rescattering in the INC, and by flux-folding in neutrino experiments.Indeed, in Ref. [23] it is seen that the SuSAv2 implementation is more similar to the EDRMF than found here, for flux-folded cross sections.Some differences remain in the hadron distributions, but note that in Ref. [23] the SuSAv2 results include the INC, while the EDRMF ones do not.We have computed the effect of rescattering by propagating the nucleon, both from the full EDRMF calculation and the GENIE procedure described above, through the nucleus using the hN INC [10].This results are shown in Fig. 2. The left panel shows the transverse momentum P T = | k T e + k T N |, with T denoting the components orthogonal to the beam direction.The angular distributions (right panel) are smeared out, but the effect of the approximation remains visible.For the P T distribution the differences are significant, and they affect the interpretation of experiments at fixed electron energy if this approximate treatment is used.

Discussion
In neutrino event generators, when only the inclusive cross section is known, an approximation is used to attach hadron variables to every event in order to describe semi-inclusive signals.We have provided an overview of the algorithm used in GENIE, e.g. in the implementations of Refs.[11,12].We have compared the nucleon observables that result from this algorithm using a local FG in quasielastic electron scattering, with calculations done in the relativistic distorted wave impulse approximation for kinematics relevant to the e4ν analysis [17].We find that the distribution of proton angles with respect to the beam and of the transverse momentum P T are significantly affected by this approximate treatment.We compute the effect of rescattering by using the hN cascade model from GENIE, and find that differences remain large, in particular for P T .
One can attribute part of these discrepancies to the use of the LFG momentum distribution.Indeed, while the LFG has been relatively successful in describing semi-inclusive flux-averaged neutrino cross sections [24], the LFG spectral function should be considered unrealistic [22] and hence fails under more restricted kinematic conditions.Additionally, due to the approximation made in selecting the nucleon angles Eq. ( 5), the missing momentum distribution of this procedure will not correspond to the LFG momentum distribution used as input.These issues can be amended by using instead a more realistic spectral function, e.g.Refs.[3,25,26], and a modification to the selection of the angle which will be presented in future work.Nonetheless, even in this case one does not obtain a fully consistent treatment suitable for every model, and in any case one loses the full kinematic and dynamical structure of the cross section [16].
The RDWIA treatment used here, possibly with the inclusion of a spectral function beyond the mean field (see Ref. [3]), does retain this structure and provides a way to study and benchmark such approximate treatments.Alternatively, the RDWIA events can be used directly in conjunction with a cascade model as done here and in Ref. [20].This approach will be used in future work to study the influence of nucleon distortion on the semi-inclusive cross section, and to establish the applicability and limitations of the cascade models used in neutrino event generators.

Figure 1 .
Figure1.Distributions of nucleon kinetic energy (left) and scattering angle with respect to the electron beam (right).Results for scattering of 1.159 GeV electrons off carbon.The red lines show the unfactorized RDWIA calculation using the ED-RMF potential, the blue lines show the result of the algorithm used in GENIE that uses the same inclusive cross section as input.

Figure 2 .
Figure 2. Distributions of transverse momentum (left) and scattering angle with respect to the electron beam (right).Results for scattering of 1.159 GeV electrons off carbon.The solid lines correspond to the same calculations as in Fig. 1.The corresponding dashed lines use in addition the GENIE hN INC model.

Funding:
Fermilab is operated by the Fermi Research Alliance, LLC under contract No. DE-AC02-07CH11359 with the United States Department of Energy.R.G.-J was supported by the government of Madrid and Complutense University under project PR65/19-22430