diff --git "a/data/slim_points/slim_points_toy_0_range\050211-4781\051.pkl" "b/data/slim_points/slim_points_toy_0_range\050211-4781\051.pkl" index dda04e4..e0d7e05 100644 --- "a/data/slim_points/slim_points_toy_0_range\050211-4781\051.pkl" +++ "b/data/slim_points/slim_points_toy_0_range\050211-4781\051.pkl" Binary files differ diff --git a/plots/points/histo.png b/plots/points/histo.png index 1b8415f..8830c35 100644 --- a/plots/points/histo.png +++ b/plots/points/histo.png Binary files differ diff --git a/plots/points/histo_raw.png b/plots/points/histo_raw.png index cb8ad04..68a3044 100644 --- a/plots/points/histo_raw.png +++ b/plots/points/histo_raw.png Binary files differ diff --git a/raremodel.py b/raremodel.py index 94cb0d8..b459a91 100644 --- a/raremodel.py +++ b/raremodel.py @@ -277,11 +277,12 @@ def total_pdf_neg(q2): return -1*self.total_pdf(q2) - maxi = fminbound(total_pdf_neg, x_min**2, x_max**2) + min_x = fminbound(total_pdf_neg, x_min**2, x_max**2) - maxi = self.total_pdf(np.array(maxi)) + maxi = self.total_pdf(np.array(min_x)) - maxi /= 1000 + if np.abs(min_x-3096) < 10 or np.abs(min_x-3686): + maxi /= 1000 for toy in range(nr_of_toys): diff --git a/test.py b/test.py index 3212990..b1a4139 100644 --- a/test.py +++ b/test.py @@ -20,7 +20,7 @@ modl.mode = mode -set_size = 1e5 +set_size = 4e5 nonres_set_size = 44000 # nonres_set_size = 1000 nr_of_toys = 1