def get_tr(src, ana=ana7, ext=np.radians(0), llh_ext=None, flux=hyp.PowerLawFlux(2), inj=True, use_energy=True, cut_n_sigma=None, sigsub=False, _fit_null=True, mp_cpus=mp_cpus): llh_ext = ext if llh_ext is None else llh_ext llh_src = cy.utils.Sources(ra=src.ra, dec=src.dec, extension=llh_ext) if cut_n_sigma is None: cut_n_sigma = 5 if llh_ext < np.radians(3) else 3 # get_llh needs the src list llh_kw = dict(src=llh_src, cut_n_sigma=cut_n_sigma, sigsub=sigsub) # get_injs needs the src list and signal spectrum inj_kw = dict(src=src, llh_src=llh_src, flux=flux, cut_n_sigma=cut_n_sigma, inj=inj) # trial.get_trial_runner loops over sub analyses and gives a single TrialRunner return trial.get_trial_runner(ana, get_llh, get_injs, llh_kw=llh_kw, inj_kw=inj_kw, mp_cpus=mp_cpus)
def get_tr (src, ana=ana7, flux=hyp.PowerLawFlux (2), inj=True, use_energy=True, cut_n_sigma=5, sigsub=False, _fit_null=True, mp_cpus=mp_cpus): # get_llh needs the src list llh_kw = dict (src=src, sigsub=sigsub) # get_injs needs the src list and signal spectrum inj_kw = dict (src=src, flux=flux, cut_n_sigma=cut_n_sigma, inj=inj) # trial.get_trial_runner loops over sub analyses and gives a single TrialRunner return trial.get_trial_runner ( ana, get_llh, get_injs, llh_kw=llh_kw, inj_kw=inj_kw, mp_cpus=mp_cpus)