cn[0] = 1.0
print "Couplings to protons [GeV^{-2}]:", cp
print "Couplings to neutrons [GeV^{-2}]:", cn
# **Note that the couplings in the NREFT are *dimensionful* and should have units of GeV$^{-2}$**.
# Now let's calculate the differential spectrum for a DM particle of mass 100 GeV:
# In[37]:
m_x = 100 #DM mass in GeV
E_list = np.linspace(0, 100, 1000)
R_list_SI = E_list * 0.0
for i, E in enumerate(E_list):
R_list_SI[i] = DMU.dRdE_NREFT(E, A_I, m_x, cp, cn, FF_I)
# Now let's plot the result, which will be in units of keV$^{-1}$ kg$^{-1}$ day$^{-1}$:
# In[38]:
pl.figure()
pl.plot(E_list, R_list_SI)
pl.xlabel(r'$E_R$ [keV]')
pl.ylabel(r'$\mathrm{d}R/\mathrm{d}E_R$ [keV$^{-1}$ kg$^{-1}$ day$^{-1}$]')
pl.title('Standard SI interactions')
pl.show()
# Now let's try something more exotic. Let's just take some random values for all the couplings to protons and neutrons (setting the couplings to $\mathcal{O}_1$ to zero):
# In[16]:
