import lhctodd
# Listing all the available DD limits 
lhctodd.list()

# Geeting the XENON1T limits
limit = dd.DD(arxiv="1805.12562")

# cross-section limit values at 100 GeV
limit.sigma(100)

# plot the limit 
limit.plot()

# get the data
data = limit.data() 
# or
data = limit.to_pandas()

axial_model = lhctodd.SD(g_chi=1.0, g_g=0.25, g_l=0.0, label="CMS MonoZ")
axial_dd = axial_model.from_csv("limit-EXO-19-003-SD-90CL.csv")
axial_model.plot()

# or

import matplotlib.pyplot as plt
plt.plot(axial_dd[:,0], axial_dd[:,1])
plt.xlabel("DM mass")
plt.ylabel("DM-Nucleon cross-section (cm2)")