def __init__(self, statics_log_filename, energies, which=-1): CanonicalBase.__init__(self, statics_log_filename, which) # Save the energies self.energies = energies # Make a histogram of the energies self.histogram = numpy.histogram(energies, bins=self.binning)[0] # In some versions of numpy, there will be an extra bin, with # values faling outside the histogram if len(self.histogram)==len(self.lnG_support)+1: self.histogram = self.histogram[:-1] # Calculate the support between lnG and the histogram self.support = self.lnG_support & (self.histogram>0) if (numpy.any(self.support)): # Calculate the bin number for each energy self.bin_number_for_energies = [] for energy in self.energies: bin_number = self.calc_bin(energy) if 0 <= bin_number < len(self.bin_centers): self.bin_number_for_energies.append(bin_number) else: self.bin_number_for_energies.append(-1)
def __init__(self, statics_log_filename, which=-1): CanonicalBase.__init__(self, statics_log_filename, which)
def __init__(self, statics_log_filename, which=-1): CanonicalBase.__init__(self, statics_log_filename, which) # Save lnG only where there is support self.lnGs = self.lnG[self.lnG_support] self.Es = self.bin_centers[self.lnG_support]