Potential Lack of Simulation of Recoil Nuclei in LArSoft
Discussing with Steve Dytman and Gianluca Petrillo, we have concluded that there seems to be a lack of vertex information associated with status 15 particles coming out of GENIE, which have been identified as final state nuclear remnants. Gianluca also confirmed that only status 1 particle are simulated by GEANT and passed to LArG4, so, naively, it seems like recoiling nuclei associated with neutrino interactions are not being simulated in our MC - at least when using GENIE. While the contribution to the neutrino interaction energy scale may be small, this seems like something we want to study in our MC and ought to be simulating (if we are indeed not currently doing so).
The GENIE authors have been notified of this, but Gianluca suggested it would be appropriate to submit a LArSoft ticket on this matter.
#2 Updated by Gianluca Petrillo about 2 years ago
- Status changed from New to Assigned
- Assignee set to Robert Hatcher
A meeting happened with Robert Hatcher, Michael Mooney, Steven Dytman, Gabriel Perdue, Erica Snider and Gianluca Petrillo.
In there, GENIE agreed to explore a solution involving to better characterise the "hardonic blob" with atomic and mass number and electric charge, and make it digestible for Geant4.
Our contact person will be Robert Hatcher.