Example #1
0
def multi_init(n, Nexp, NMC, **kwargs):
    energy = kwargs.pop("energy", False)
    seed = kwargs.pop("seed", 0)
    llh = MultiPointSourceLLH(seed=seed, **kwargs)

    P = np.random.RandomState(seed)

    for i in range(n):
        llh.add_sample("S{0:02d}".format(i), init(P.poisson(Nexp // n),
                                                  NMC, energy=energy,
                                                  seed=seed + i + 1))

    return llh
Example #2
0
def multi_init(n, Nexp, NMC, **kwargs):
    energy = kwargs.pop("energy", False)

    llh = MultiPointSourceLLH(nsource=25,
                              nsource_bounds=(-Nexp / 2.,
                                              Nexp / 2.) if not energy else
                              (0., Nexp / 2.),
                              seed=np.random.randint(2**32),
                              **kwargs)

    for i in xrange(n):
        llh.add_sample(str(i), init(Nexp, NMC, energy=energy))

    return llh
Example #3
0
def multi_init(n, Nexp, NMC, **kwargs):
    energy = kwargs.pop("energy", False)

    llh = MultiPointSourceLLH(hemispheres=dict(Full=[-np.inf, np.inf]),
                              nsource=25,
                              nsource_bounds=(-Nexp / 2., Nexp / 2.)
                                             if not energy else (0., Nexp / 2.),
                              seed=np.random.randint(2**32),
                              **kwargs)

    for i in xrange(n):
        llh.add_sample(str(i), init(Nexp, NMC, energy=energy))

    return llh