evanoconnor / NuLib

open-source neutrino interaction library
23 stars 16 forks source link

beta decay #32

Closed evanoconnor closed 5 years ago

evanoconnor commented 5 years ago

This PR includes some small updates in order to make beta decays (e-) and positron captures work in the NuLib weakrates library. Wright et al. updated NuLib to include antineutrino from weak interactions and these changes are taken from their work.

Previously, the infrastructure was in place but it halted upon usage. There were a few changes required to make it consistent. First, when looking if a rate exists for a Z,A nucleus, we must look for the Z+1,A electron capture rate entry, that is where the 'negative' direction data is stored. Also when interpolating the rate we must look there too. Note, the approximation is not working the beta decay direction. Changes will be needed for that.