In [1]:

    

import matplotlib.pylab as plt
import numpy as np

import dm_models as dm

from time import time


    

In [2]:

    

target_atom = dm.AGe
target_mass = 1.0 # Mass of target in kg

massDM = 10.0   # GeV
sigma_n = 7e-42 # WIMP-nucleon cross section cm^2
elo = 0.5       # keVee
ehi = 3.2       # keVee 

max_days = 365 # Number of days that you collect data

efficiency = lambda x: 1.0

model = 'shm' # Velocity distribution model: 
              # "shm": Simple Halo Model. Maxwell-Boltzman distribution
              # "stream": Dark matter stream. Sagittarius stream parameters by default.

start = time()
num_wimps = dm.expected_number_of_WIMPs(wimp_model=model,
                                        loE=elo,hiE=ehi,
                                        loDay=1,hiDay=max_days,
                                        target_atom=target_atom,
                                        target_mass=target_mass,
                                        massDM=massDM,
                                        sigma_n=sigma_n,
                                        efficiency=efficiency)
stop = time()

print "# expected WIMPs: %f" % (num_wimps)

print "time to run: %f ms" % (stop-start)


    

    




# expected WIMPs: 341.339780
time to run: 0.242491 ms

In [ ]: