Next Article in Journal
Energy Reconstruction and Calibration of the MicroBooNE LArTPC
Previous Article in Journal
Can Deviation from Maximal θ23 Be Resolved in DUNE?
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Proceeding Paper

From Inclusive to Semi-Inclusive One-Nucleon Knockout in Neutrino Event Generators †

by
Alexis Nikolakopoulos
1,*,
Steven Gardiner
1,
Afroditi Papadopoulou
2,
Stephen Dolan
3 and
Raúl González-Jiménez
4
1
Fermi National Accelerator Laboratory (Fermilab), P.O. Box 500, Batavia, IL 60510, USA
2
Argonne National Laboratory (ANL), Lemont, IL 60439, USA
3
European Organization for Nuclear Research (CERN), 1211 Geneva, Switzerland
4
Grupo de Física Nuclear, Departamento de Estructura de la Materia, Física Térmica y Electrónica, Facultad de Ciencias Físicas, Universidad Complutense de Madrid and IPARCOS, CEI Moncloa, 28040 Madrid, Spain
*
Author to whom correspondence should be addressed.
Presented at the 23rd International Workshop on Neutrinos from Accelerators, Salt Lake City, UT, USA, 30–31 July 2022.
Phys. Sci. Forum 2023, 8(1), 10; https://doi.org/10.3390/psf2023008010
Published: 13 July 2023
(This article belongs to the Proceedings of The 23rd International Workshop on Neutrinos from Accelerators)

Abstract

:
In neutrino event generators, for models for neutrino and electron scattering, only inclusive cross sections are implemented. When these models are used to describe a semi-inclusive cross section, the event generator attaches the hadron variables based on some assumptions. In this work, we compared the nucleon kinematics given by the method used in the GENIE event generator, e.g., in the implementation of the SuSAv2 model, to a fully unfactorized calculation using the relativistic distorted wave impulse approximation (RDWIA). We focused on kinematics relevant to the e 4 ν analysis and showed that observables obtained with RDWIA differ significantly from those of the approximate method used in GENIE; the latter should be considered unrealistic.

1. 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 re-interacts with the nucleus. These strong interactions are dubbed final-state interactions (FSI). Within this picture, inelastic FSI, where the energy is 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 considered the contribution of quasielastic interactions, where the exchanged boson is absorbed on a single nucleon that is excited to the continuum, to a 1-lepton, 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 that 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 compared this result to fully unfactorized calculations with the relativistic distorted-wave impulse approximation (RDWIA). We focused on the ( e , e p ) process for kinematics relevant to the e 4 ν analysis [17]. We showed that the approximate procedure used in GENIE leads to significantly different observables than the RDWIA result.

2. Results

We used 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 100 MeV ), 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.159 GeV . The events are distributed according to the ( Q 2 ) 2 -weighted cross section obtained in the RDWIA (for details, see Refs. [3,20]).
P ( E l , θ l , T N , Ω N ) = Q 2 1 GeV 2 2 M B d 4 σ ( E e , M B ) d E e d cos θ e d Ω N .
The ( Q 2 ) 2 weighting, which was introduced in the e 4 ν 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 imposed cuts for the outgoing electron 40 ° > θ e > 17 ° and E e > 400 MeV but included the full nucleon phase space.
To test the approximate treatment, we replaced 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 used the local Fermi gas (FG) obtained from the nuclear density, also taken from GENIE,
ρ ( r ) = N 1 + r a 2 α e r a 2 ,
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
E N = p m 2 + M N 2 + ω E b ( q )
if p m 2 + M N 2 M N E b ( q ) is negative; otherwise, a new p m is generated. Here, E b ( q ) is a q-dependent binding energy inspired by the energy shift from the SuSAv2 model,
E b = max ( 5 , 17.687 + 0.0564 q ) [ MeV ] q < 827 MeV ,
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 Figure 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 Figure 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 Equation (3) will not agree with momentum conservation, E N 2 M N 2 | p m + q | . To impose momentum conservation, the magnitude of the nucleon momentum is taken from Equation (3), while the direction is taken from the momentum vectors, i.e.,
k N = E N 2 M N 2 q + p m | q + p m | .
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 ‘ h N ’ INC [10]. These results are shown in Figure 2. The left panel shows the transverse momentum P T = | k e T + k N T | , 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.

3. 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 performed in the relativistic distorted wave impulse approximation for kinematics relevant to the e 4 ν 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 computed the effect of rescattering by using the ‘ h N ’ 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 in Equation (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 of 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 undertaken 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.

Author Contributions

Conceptualization, A.N. and S.D.; Methodology, A.N., S.G. and R.G.-J.; Software, A.N., S.G., A.P. and R.G.-J.; Writing—original draft, A.N.; Writing—review & editing, A.N., S.G., A.P., S.D. and R.G.-J. All authors have read and agreed to the published version of the manuscript.

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. This material is based upon work supported by Laboratory Directed Research and Development (LDRD) funding from Argonne National Laboratory, provided by the Director, Office of Science, of the U.S. Department of Energy under Contract No. DE-AC02-06CH11357.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The event distributions used for this work, and a C++ program to compute alternative nucleon variables are available from A. Nikolakopoulos upon request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lu, X.G.; Pickering, L.; Dolan, S.; Barr, G.; Coplowe, D.; Uchida, Y.; Wark, D.; Wascko, M.O.; Weber, A.; Yuan, T. Measurement of nuclear effects in neutrino interactions with minimal dependence on neutrino energy. Phys. Rev. C 2016, 94, 015503. [Google Scholar] [CrossRef] [Green Version]
  2. Dolan, S. Exploring Nuclear Effects in Neutrino-Nucleus Interactions Using Measurements of Transverse Kinematic Imbalance from T2K and MINERvA. arXiv 2018, arXiv:hep-ex/1810.06043. [Google Scholar]
  3. González-Jiménez, R.; Barbaro, M.B.; Caballero, J.A.; Donnelly, T.W.; Jachowicz, N.; Megias, G.D.; Niewczas, K.; Nikolakopoulos, A.; Van Orden, J.W.; Udías, J.M. Neutrino energy reconstruction from semi-inclusive samples. Phys. Rev. C 2022, 105, 025502. [Google Scholar] [CrossRef]
  4. Golan, T.; Juszczak, C.; Sobczyk, J.T. Effects of final-state interactions in neutrino-nucleus interactions. Phys. Rev. C 2012, 86, 015505. [Google Scholar] [CrossRef] [Green Version]
  5. Hayato, Y. A neutrino interaction simulation program library NEUT. Acta Phys. Polon. 2009, B40, 2477–2489. [Google Scholar] [CrossRef]
  6. Hayato, Y.; Pickering, L. The NEUT neutrino interaction simulation program library. Eur. Phys. J. Spec. Top. 2021, 230, 4469–4481. [Google Scholar] [CrossRef]
  7. Salcedo, L.; Oset, E.; Vicente-Vacas, M.; Garcia-Recio, C. Computer simulation of inclusive pion nuclear reactions. Nucl. Phys. A 1988, 484, 557–592. [Google Scholar] [CrossRef]
  8. Andreopoulos, C.; Bell, A.; Bhattacharya, D.; Cavanna, F.; Dobson, J.; Dytman, S.; Gallagher, H.; Guzowski, P.; Hatcher, R.; Kehayias, P.; et al. The GENIE Neutrino Monte Carlo Generator. Nucl. Instrum. Meth. 2010, A614, 87–104. [Google Scholar] [CrossRef] [Green Version]
  9. Isaacson, J.; Jay, W.I.; Lovato, A.; Machado, P.A.N.; Rocco, N. ACHILLES: A novel event generator for electron- and neutrino-nucleus scattering. arXiv 2022, arXiv:hep-ph/2205.06378. [Google Scholar]
  10. Dytman, S.; Hayato, Y.; Raboanary, R.; Sobczyk, J.T.; Tena Vidal, J.; Vololoniaina, N. Comparison of validation methods of simulations for final state interactions in hadron production experiments. Phys. Rev. D 2021, 104, 053006. [Google Scholar] [CrossRef]
  11. Dolan, S.; Nikolakopoulos, A.; Page, O.; Gardiner, S.; Jachowicz, N.; Pandey, V. Implementation of the continuum random phase approximation model in the GENIE generator and an analysis of nuclear effects in low-energy transfer neutrino interactions. Phys. Rev. D 2022, 106, 073001. [Google Scholar] [CrossRef]
  12. Dolan, S.; Megias, G.D.; Bolognesi, S. Implementation of the SuSAv2-MEC 1p1h and 2p2h models in GENIE and analysis of nuclear effects in T2K measurements. arXiv 2019, arXiv:hep-ex/1905.08556. [Google Scholar]
  13. Nieves, J.; Simo, I.R.; Vacas, M.J.V. Inclusive charged-current neutrino-nucleus reactions. Phys. Rev. C 2011, 83, 045501. [Google Scholar] [CrossRef] [Green Version]
  14. Nieves, J.; Simo, I.R.; Vacas, M.V. Two particle–hole excitations in charged current quasielastic antineutrino-nucleus scattering. Phys. Lett. B 2013, 721, 90–93. [Google Scholar] [CrossRef] [Green Version]
  15. Martini, M.; Ericson, M. Quasielastic and multinucleon excitations in antineutrino-nucleus interactions. Phys. Rev. C 2013, 87, 065501. [Google Scholar] [CrossRef] [Green Version]
  16. Moreno, O.; Donnelly, T.W.; Van Orden, J.W.; Ford, W.P. Semi-inclusive charged-current neutrino-nucleus reactions. Phys. Rev. D 2014, 90, 013014. [Google Scholar] [CrossRef] [Green Version]
  17. Khachatryan, M.; Papadopoulou, A.; Ashkenazi, A.; Hauenstein, F.; Nambrath, A.; Hrnjic, A.; Weinstein, L.B.; Hen, O.; Piasetzky, E.; Betancourt, M.; et al. Electron-beam energy reconstruction for neutrino oscillation measurements. Nature 2021, 599, 565–570. [Google Scholar] [CrossRef]
  18. González-Jiménez, R.; Nikolakopoulos, A.; Jachowicz, N.; Udías, J.M. Nuclear effects in electron-nucleus and neutrino-nucleus scattering within a relativistic quantum mechanical framework. Phys. Rev. C 2019, 100, 045501. [Google Scholar] [CrossRef]
  19. González-Jiménez, R.; Barbaro, M.; Caballero, J.; Donnelly, T.; Jachowicz, N.; Megias, G.; Niewczas, K.; Nikolakopoulos, A.; Udías, J. Constraints in modeling the quasielastic response in inclusive lepton-nucleus scattering. Phys. Rev. C 2020, 101, 015503. [Google Scholar] [CrossRef] [Green Version]
  20. Nikolakopoulos, A.; González-Jiménez, R.; Jachowicz, N.; Niewczas, K.; Sánchez, F.; Udías, J.M. Benchmarking intranuclear cascade models for neutrino scattering with relativistic optical potentials. Phys. Rev. C 2022, 105, 054603. [Google Scholar] [CrossRef]
  21. Abe, K.; Amey, J.; Andreopoulos, C.; Anthony, L.; Antonova, M.; Aoki, S.; Ariga, A.; Ashida, Y.; Azuma, Y.; Ban, S.; et al. Characterization of nuclear effects in muon-neutrino scattering on hydrocarbon with a measurement of final-state kinematics and correlations in charged-current pionless interactions at T2K. Phys. Rev. D 2018, 98, 032003. [Google Scholar] [CrossRef] [Green Version]
  22. Van Orden, J.W.; Donnelly, T.W. Nuclear theory and event generators for charge-changing neutrino reactions. Phys. Rev. C 2019, 100, 044620. [Google Scholar] [CrossRef] [Green Version]
  23. Franco-Patino, J.M.; González-Jiménez, R.; Dolan, S.; Barbaro, M.B.; Caballero, J.A.; Megias, G.D.; Udias, J.M. Final state interactions in semi-inclusive neutrino-nucleus scattering: Application to T2K and MINERνA experiments. arXiv 2022, arXiv:nucl-th/2207.02086. [Google Scholar] [CrossRef]
  24. Bourguille, B.; Nieves, J.; Sánchez, F. Inclusive and exclusive neutrino-nucleus cross sections and the reconstruction of the interaction kinematics. JHEP 2021, 4, 4. [Google Scholar] [CrossRef]
  25. Benhar, O.; Farina, N.; Nakamura, H.; Sakuda, M.; Seki, R. Electron- and neutrino-nucleus scattering in the impulse approximation regime. Phys. Rev. D 2005, 72, 053005. [Google Scholar] [CrossRef] [Green Version]
  26. Rocco, N.; Barbieri, C.; Benhar, O.; De Pace, A.; Lovato, A. Neutrino-Nucleus Cross Section within the Extended Factorization Scheme. Phys. Rev. C 2019, 99, 025502. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Distributions of nucleon kinetic energy (left) and scattering angle with respect to the electron beam (right). Results for scattering of 1.159 GeV electrons of 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 1. Distributions of nucleon kinetic energy (left) and scattering angle with respect to the electron beam (right). Results for scattering of 1.159 GeV electrons of 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.
Psf 08 00010 g001
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 of carbon. The solid lines correspond to the same calculations as in Figure 1. The corresponding dashed lines use in addition the GENIE hN INC model.
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 of carbon. The solid lines correspond to the same calculations as in Figure 1. The corresponding dashed lines use in addition the GENIE hN INC model.
Psf 08 00010 g002
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Nikolakopoulos, A.; Gardiner, S.; Papadopoulou, A.; Dolan, S.; González-Jiménez, R. From Inclusive to Semi-Inclusive One-Nucleon Knockout in Neutrino Event Generators. Phys. Sci. Forum 2023, 8, 10. https://doi.org/10.3390/psf2023008010

AMA Style

Nikolakopoulos A, Gardiner S, Papadopoulou A, Dolan S, González-Jiménez R. From Inclusive to Semi-Inclusive One-Nucleon Knockout in Neutrino Event Generators. Physical Sciences Forum. 2023; 8(1):10. https://doi.org/10.3390/psf2023008010

Chicago/Turabian Style

Nikolakopoulos, Alexis, Steven Gardiner, Afroditi Papadopoulou, Stephen Dolan, and Raúl González-Jiménez. 2023. "From Inclusive to Semi-Inclusive One-Nucleon Knockout in Neutrino Event Generators" Physical Sciences Forum 8, no. 1: 10. https://doi.org/10.3390/psf2023008010

APA Style

Nikolakopoulos, A., Gardiner, S., Papadopoulou, A., Dolan, S., & González-Jiménez, R. (2023). From Inclusive to Semi-Inclusive One-Nucleon Knockout in Neutrino Event Generators. Physical Sciences Forum, 8(1), 10. https://doi.org/10.3390/psf2023008010

Article Metrics

Back to TopTop