#import WimPyDD as WD
import numpy as np
import matplotlib.pyplot as pl

# plots the H(v) response function for halo-independent calculations
vmin=np.linspace(1e-3,782,100)
pl.clf()
energy,h=WD.wimp_dd_rate(WD.XENON_1T_2018, WD.SI, vmin, 1/vmin, 100,sum_over_streams=False,sigma_p=9.066504208570347e-47)
pl.plot(vmin,h[0])
pl.text(500,15,r'$m_{\chi}$ = 100 GeV',fontsize=12)
pl.text(500,13,r'${\sigma}^p = 9.06 \times 10^{-47} cm^2$',fontsize=12)
pl.xlabel(r'$v_{min}$ (km/s)',size=12)
pl.ylabel('H(v) (events)',size=12)
pl.show()