In [1]:

    

import dmdd
import numpy as np
%matplotlib inline


    

    




WARNING:root:pymultinest not imported!
WARNING:root:DMDD_MAIN_PATH environment variable not defined, defaulting to:   ~/.dmdd

In [2]:

    

# shortcut for scattering models corresponding to rates coded in rate_UV:
anapole_model = dmdd.UV_Model('Anapole', ['mass','sigma_anapole'])
SI_model = dmdd.UV_Model('SI', ['mass','sigma_si'])

print 'model: {}, parameters: {}.'.format(anapole_model.name, anapole_model.param_names)
print 'model: {}, parameters: {}.'.format(SI_model.name, SI_model.param_names)


    

    




model: Anapole, parameters: ['mass', 'sigma_anapole'].
model: SI, parameters: ['mass', 'sigma_si'].

In [3]:

    

# intialize an Experiment with XENON target, to be passed to Simulation_AM:
xe_low = dmdd.Experiment('1xe', 'xenon', 5, 30, 1000, dmdd.eff.efficiency_unit, energy_resolution=True)


    

In [5]:

    

xe_si = dmdd.Simulation_AM('SI', xe_low, SI_model, 
                        {'mass':50.,'sigma_si':200.}, Qmin = np.asarray([5.]), 
                        Qmax = np.asarray([30.]), 
                        Tmin = 0, Tmax = 365, sigma_si = 200., 
                        element = 'xenon', force_sim = True)
#make large cross section and low Q range to get as many points as possible to fill all possible bins


    

    




Simulation data and/or pickle file does not exist. Forcing simulation.


simulated: 414 events (expected 419).






    




//anaconda/lib/python2.7/site-packages/matplotlib/image.py:359: UserWarning: Images are not supported on non-linear axes.
  warnings.warn("Images are not supported on non-linear axes.")
//anaconda/lib/python2.7/site-packages/matplotlib/transforms.py:2547: RuntimeWarning: divide by zero encountered in double_scalars
  x_scale = 1.0 / inw
//anaconda/lib/python2.7/site-packages/matplotlib/transforms.py:2548: RuntimeWarning: divide by zero encountered in double_scalars
  y_scale = 1.0 / inh






    

In [ ]: