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)
Beispiel #2
0
    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]
Beispiel #5
0
 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]