def dRdE(ER_keV, m_x, cp, cn, **kwargs):
#Load in the list of nuclear spins, atomic masses and mass fractions
nuclei_Xe = ["Xe128", "Xe129", "Xe130", "Xe131", "Xe132", "Xe134", "Xe136"]
nuclei_list = np.loadtxt("Nuclei.txt", usecols=(0,), dtype='string')
frac_list = np.loadtxt("Nuclei.txt", usecols=(3,))
frac_vals = dict(zip(nuclei_list, frac_list))
dRdE = np.zeros_like(ER_keV)
for nuc in nuclei_Xe:
dRdE += frac_vals[nuc]*DMU.dRdE_NREFT(ER_keV, m_x, cp, cn, nuc, **kwargs)
return dRdE
def dRdEAr(m, E, c, eff, Nevents=True, **kwargs):
s = eff(E) * DMU.dRdE_NREFT(E, m, c, c, "Ar40", **kwargs) + 1.e-30
if Nevents:
s, sum(s * ER_width)
else:
return s
def dRdE(m, E, c, eff):
s = eff(E) * DMU.dRdE_NREFT(E, m, c, c, "Ar40") + 1.e-30
return s