---------------------------------------------------------------------------
KeyboardInterrupt Traceback (most recent call last)
<ipython-input-4-39e078dcf6ce> in <module>()
3 Qmax = np.asarray([80.]),
4 Tmin = 0, Tmax = 365, sigma_si = 75.5,
----> 5 element = 'xenon', force_sim = True)
6
7 #maybe use a histogram to see what it looks like?
/Users/katelynneese/dmddACT/dmdd/dmdd.py in __init__(self, name, experiment, model, parvals, Qmin, Qmax, element, sigma_si, sigma_anapole, mass, v_amplitude, Tmin, Tmax, path, force_sim, asimov, nbins_asimov, plot_nbins, plot_theory, silent)
1068 raise ValueError('Asimov simulations not yet implemented!')
1069 else:
-> 1070 Q = self.simulate_data()
1071 np.savetxt(self.datafile,Q)
1072 fout = open(self.picklefile,'wb')
/Users/katelynneese/dmddACT/dmdd/dmdd.py in simulate_data(self)
1134 pdf_value = PDF(Q_rand, T_rand, element = self.element, mass = self.mass,
1135 sigma_si= self.sigma_si, sigma_anapole = self.sigma_anapole,
-> 1136 Qmin = np.asarray([self.Qmin]), Qmax = np.asarray([self.Qmax]), Tmin = self.Tmin, Tmax = self.Tmax)/env
1137 if U < pdf_value:
1138 #increment matches
/Users/katelynneese/dmddACT/dmdd/dmdd.py in PDF(Q, time, element, mass, sigma_si, sigma_anapole, Qmin, Qmax, Tmin, Tmax)
117 Tmin = Tmin, Tmax = Tmax,
118 element = element, sigma_si = sigma_si,
--> 119 sigma_anapole = sigma_anapole, mass = mass)
120
121 return drdq/norm[0] #for now removed efficiency and made it 1
/Users/katelynneese/dmddACT/dmdd/dmdd.py in integral(Qmin, Qmax, Tmin, Tmax, Qpoints, Tpoints, function, sigma_si, sigma_anapole, mass, element, v_amplitude)
102 for j,T in enumerate(T_box):
103 a_sum = np.sum(function(Q = np.asarray(Q), time = T, sigma_si = sigma_si, sigma_anapole = sigma_anapole,
--> 104 mass = mass, element = element, v_amplitude = v_amplitude))
105 total_sum.append(a_sum)
106
/Users/katelynneese/dmddACT/dmdd/dmdd.py in dRdQ_AM(mass, sigma_si, sigma_anapole, Q, time, element, vlag_mean, v_amplitude)
81 #print type(v_lag) #vlag must be a number not an array
82
---> 83 rate_QT = rate_UV.dRdQ(Q = energy, v_lag = v_lag, mass = mass, sigma_si = sigma_si, sigma_anapole = sigma_anapole, element = element)
84
85 "Return a 1D array with the rate based on the time and energy given"
KeyboardInterrupt: