def test_blh_with_parts():
np.random.seed(0)
data = np.random.randn(10000)
shifted = data + 3.
data = np.append(data, [shifted])
g1 = rename(gaussian, ['x', 'lmu', 'lsigma'])
g2 = rename(gaussian, ['x', 'rmu', 'rsigma'])
allpdf = AddPdfNorm(g1, g2)
blh = BinnedLH(allpdf, data)
blh.draw(args=(0, 1, 3, 1, 0.5), parts=True)
def test_ulh_with_parts():
np.random.seed(0)
data = np.random.randn(10000)
shifted = data + 3.0
data = np.append(data, [shifted])
g1 = rename(gaussian, ["x", "lmu", "lsigma"])
g2 = rename(gaussian, ["x", "rmu", "rsigma"])
allpdf = AddPdfNorm(g1, g2)
ulh = UnbinnedLH(allpdf, data)
ulh.draw(args=(0, 1, 3, 1, 0.5), parts=True)