diff --git a/PeakOTron.py b/PeakOTron.py index 9baa771be0e870ed2d9c39f2a3c97ffb4b47cc84..514f8ea7d097795f9c23f65dd165a5de3bb2d989 100644 --- a/PeakOTron.py +++ b/PeakOTron.py @@ -70,15 +70,13 @@ class BandWidthOptimiser: - def __init__(self, data, kernel="gaussian", alpha = 0.99, n_kde_samples=2**13): + def __init__(self, data, kernel="gaussian", alpha = 0.95, n_kde_samples=2**13): if(alpha<=0.5): raise ValueError("alpha must be above 0.5") self.data = data self.PPF = self._PPF(self.data) - self.x_low = self.PPF(0.5-alpha/2) - self.x_hi = self.PPF(0.5+alpha/2) - self.data = self.data[(self.data>self.x_low) & (self.data<self.x_hi)] + self.data = self.data[(self.data<self.PPF(alpha))] self.kernel=kernel self.n_kde_samples=n_kde_samples self.last_arg = None @@ -94,6 +92,7 @@ class BandWidthOptimiser: return loss + class BinnedLH: