search

GENIE-MC / Generator

The popular GENIE Generator product is used by nearly all accelerator neutrino experiments and it plays a key role in the exploitation of neutrino data. The Generator implements a modern software framework and it includes state-of-the-art physics modules. It captures the latest results of the GENIE global analysis of neutrino scattering data and includes several tunes that were produced using the proprietary Comparisons and Tuning products. The GENIE physics model is universal and comprehensive: It handles all neutrinos and targets, and all processes relevant from MeV to PeV energy scales. The Generator includes several tools (flux drivers, detector geometry navigators, specialized event generation apps, event reweighting engines) to simulate complex experimental setups in full detail and to support generator-related analysis tasks.
http://www.genie-mc.org
42 stars 92 forks source link

Apparent energy conservation issue for MEC events when G4 FSI model is used #313

Closed wilkinson-nu closed 8 months ago

wilkinson-nu commented 1 year ago

I'll preface this by saying this is an issue I saw with v3.2.0 with G18_10*_00_000 tunes, but failed to report properly. Some studies Richie showed in the context of DUNE using the new AR23_20i_00_000 model, but playing with the FSI model prompted me to check again.

The plots below are using files Richie made and put on the DUNE gpvms and converted to GST format, which can be found here: /pnfs/dune/persistent/users/rdiurba/DUNEND_FHC_numu_gst/, just in case anybody wants to check.

Looking at the total kinetic energy of protons and neutrons (labelled Ehad here) from CCMEC events vs true q0 (Enu - Elep) for the hA2018 FSI model (AR23_20i_00_000), you get a plot that looks sensible: image This makes sense to me, there's some smearing due to the initial state nucleon momentum, but ignoring that, Ehad < q0. If anybody wants to reproduce these plots, you can probably guess the TTree draw command from the histogram title.

With the G4 FSI model (labelled AR23_20l_00_000 in Richie's files, which may be the usual convention, I'm not 100% sure): image The obvious second population doesn't make sense to me, and looks like an energy conservation issue, although presumably there's a more subtle issue, it just looks like that.

I haven't poked into this in any more detail, so apologies for this slightly unhelpful issue report... I just wanted to flag this up.

mroda88 commented 1 year ago

Hi, thanks for the report. We are investigating a number of issues related to G4 and INCL++ actually, and I think this is somehow related to those issues. I agree with you that for now those FSI should not be used for production.

mroda88 commented 10 months ago

The fix should be in #360

wilkinson-nu commented 8 months ago

I finally got around the checking the update, and it does indeed seem to have resolved this issue, so I'll close it now. Thanks!