Loading peakintegration/peakfitter.py +21 −17 Original line number Diff line number Diff line Loading @@ -1905,32 +1905,36 @@ class PeaksIntegrator(object): # exit() # plt.close('all') # plt.subplot(211) mantdI = self.peak_mantidI[peak_ids_main] bayesI = np.array([peak.getIntensity() for peak in self.peaks])[peak_ids_main] if self.save_path is not None: with open(make_path(self.save_path,'main_I.npy'), 'wb') as f: np.save(f, mantdI) np.save(f, bayesI) # mantdI = self.peak_mantidI[peak_ids_main] # bayesI = np.array([peak.getIntensity() for peak in self.peaks])[peak_ids_main] # if self.save_path is not None: # with open(make_path(self.save_path,'main_I.npy'), 'wb') as f: # np.save(f, mantdI) # np.save(f, bayesI) # ind = np.argsort(bayesI) # ind2 = np.argsort(mantdI) # plt.semilogy(mantdI[ind], label='Mantid') # plt.semilogy(bayesI[ind], label='Bayesian') # plt.legend() # plt.subplot(212) mantdI = self.peak_mantidI[peak_ids_sat] bayesI = np.array([peak.getIntensity() for peak in self.peaks])[peak_ids_sat] # print(np.array([peak.getIntensity() for peak in self.peaks])[313]) # print(self.peak_mantidI[313]) # exit() if self.save_path is not None: with open(make_path(self.save_path,'sat_I.npy'), 'wb') as f: np.save(f, mantdI) np.save(f, bayesI) # mantdI = self.peak_mantidI[peak_ids_sat] # bayesI = np.array([peak.getIntensity() for peak in self.peaks])[peak_ids_sat] # print(np.array([peak.getIntensity() for peak in self.peaks])[385]) # print(self.peak_mantidI[385]) # print(self.main_ids[385]) # # exit() # if self.save_path is not None: # with open(make_path(self.save_path,'sat_I.npy'), 'wb') as f: # np.save(f, mantdI) # np.save(f, bayesI) # ind = np.argsort(bayesI) # ind2 = np.argsort(mantdI) # # pritn(ind2) # plt.semilogy(mantdI[ind2], label='Mantid') # k = np.argmax(mantdI[ind]>100) # print(k, np.array(peak_ids_sat)[ind][k], np.array(peak_ids_sat)[ind]) # plt.semilogy(mantdI[ind2], '.', label='Mantid') # plt.semilogy(bayesI[ind], label='Bayesian') # plt.legend() # plt.show() Loading Loading
peakintegration/peakfitter.py +21 −17 Original line number Diff line number Diff line Loading @@ -1905,32 +1905,36 @@ class PeaksIntegrator(object): # exit() # plt.close('all') # plt.subplot(211) mantdI = self.peak_mantidI[peak_ids_main] bayesI = np.array([peak.getIntensity() for peak in self.peaks])[peak_ids_main] if self.save_path is not None: with open(make_path(self.save_path,'main_I.npy'), 'wb') as f: np.save(f, mantdI) np.save(f, bayesI) # mantdI = self.peak_mantidI[peak_ids_main] # bayesI = np.array([peak.getIntensity() for peak in self.peaks])[peak_ids_main] # if self.save_path is not None: # with open(make_path(self.save_path,'main_I.npy'), 'wb') as f: # np.save(f, mantdI) # np.save(f, bayesI) # ind = np.argsort(bayesI) # ind2 = np.argsort(mantdI) # plt.semilogy(mantdI[ind], label='Mantid') # plt.semilogy(bayesI[ind], label='Bayesian') # plt.legend() # plt.subplot(212) mantdI = self.peak_mantidI[peak_ids_sat] bayesI = np.array([peak.getIntensity() for peak in self.peaks])[peak_ids_sat] # print(np.array([peak.getIntensity() for peak in self.peaks])[313]) # print(self.peak_mantidI[313]) # exit() if self.save_path is not None: with open(make_path(self.save_path,'sat_I.npy'), 'wb') as f: np.save(f, mantdI) np.save(f, bayesI) # mantdI = self.peak_mantidI[peak_ids_sat] # bayesI = np.array([peak.getIntensity() for peak in self.peaks])[peak_ids_sat] # print(np.array([peak.getIntensity() for peak in self.peaks])[385]) # print(self.peak_mantidI[385]) # print(self.main_ids[385]) # # exit() # if self.save_path is not None: # with open(make_path(self.save_path,'sat_I.npy'), 'wb') as f: # np.save(f, mantdI) # np.save(f, bayesI) # ind = np.argsort(bayesI) # ind2 = np.argsort(mantdI) # # pritn(ind2) # plt.semilogy(mantdI[ind2], label='Mantid') # k = np.argmax(mantdI[ind]>100) # print(k, np.array(peak_ids_sat)[ind][k], np.array(peak_ids_sat)[ind]) # plt.semilogy(mantdI[ind2], '.', label='Mantid') # plt.semilogy(bayesI[ind], label='Bayesian') # plt.legend() # plt.show() Loading
