MARLEYGen module fails to use loaded nuclear structure data when simulating de-excitations
The larsim MarleyGen module is used to configure the MARLEY supernova neutrino generator for use with LArSoft. In versions v06_04_00 through v06_18_00 of LArSoft, the module loads nuclear structure data (level energies, gamma-ray branching ratios, etc.) into MARLEY, but the events produced using the module fail to use these data. Instead, MARLEY falls back to using a continuum de-excitation model (the Hauser-Feshbach statistical model) for all level energies, including bound states. The command-line MARLEY executable included with LArSoft, however, gives the correct behavior when identical configuration file settings are used.
Fix bug that prevents MARLEY from using nuclear structure data loaded using the larsim MARLEYGen module. This resolves LArSoft issue #14847.
#2 Updated by Steven Gardiner almost 3 years ago
- Status changed from Assigned to Resolved
- % Done changed from 0 to 100
- Occurs In 2016-4-quarter added
- Occurs In deleted (
I've pushed a bug fix to the larsim develop branch. In MARLEY v0_9_5, this problem can occur if a marley::Generator object initializes a marley::NuclearReaction that it owns before its marley::StructureDatabase is completely initialized. MARLEY will be modified in an upcoming release to eliminate this dependence on initialization order, but LArSoft will function correctly with v0_9_5 as of larsim commit d7c9bf5d31cda419820d6508b43d33ad0a809703, which resolves this issue.