ret = np.sum(-1 * np.log(pdf(p, x, fixed_parameters))) - poisson_term return ret ################################################################################ # Importing the dark matter data ################################################################################ # # This is the ``time" of the first event, in some digitial units used # by the detector first_event = 2750361.2 # Full path to the directory infile_name = '/Users/lm27apic/Documents/Fall_2012/Dark_Matter_Research/dark_matter_data/low_gain.txt' tdays, energies = cu.get_cogent_data(infile_name, first_event=first_event, calibration=999) xlo = 0.133 xhi = 0.1470 mu = (xlo + xhi) / 2 sigma = 0.0085 slope = 0.5 # Index the range of energies index0 = energies > xlo index1 = energies < xhi index = index0 * index1 x = energies[index]
ret = np.sum(-1*np.log(pdf(p,x,fixed_parameters))) - poisson_term return ret ################################################################################ # Importing the dark matter data ################################################################################ # # This is the ``time" of the first event, in some digitial units used # by the detector first_event = 2750361.2 # Full path to the directory infile_name = '/Users/lm27apic/Documents/Fall_2012/Dark_Matter_Research/dark_matter_data/low_gain.txt' tdays,energies = cu.get_cogent_data(infile_name,first_event=first_event,calibration=999) xlo = 0.133 xhi = 0.1470 mu = (xlo + xhi) / 2 sigma = 0.0085 slope = 0.5 # Index the range of energies index0 = energies>xlo index1 = energies<xhi index = index0*index1 x = energies[index] plt.figure()