Exemple #1
0
def test_short():
    from jetset.plot_sedfit import plt
    plt.ioff()
    test_jet()
    sed_data = data()
    print('done')
    my_shape = spectral_indices(sed_data)
    print('done')
    my_shape = sed_shaper(my_shape)
    print('done')
    prefit_jet = model_constr(my_shape)
    print('done')
    model_fit_lsb(sed_data, my_shape)
def test_full():
    from jetset.plot_sedfit import plt
    plt.ioff()
    test_jet()
    sed_data = data()
    print('done')
    my_shape = spectral_indices(sed_data)
    print('done')
    my_shape = sed_shaper(my_shape)
    print('done')
    prefit_jet = model_constr(my_shape)
    print('done')
    jet_lsb, model_minimizer_lsb, fit_model_lsb = model_fit_lsb(
        sed_data, my_shape)
    jet_minuit, model_minimizer_minuit = model_fit_minuit(sed_data, my_shape)
def test_short(plot=False):
    from jetset.plot_sedfit import plt
    try:
        plt.ioff()
        test_jet(plot)
        sed_data = data(plot)
        print('done')
        my_shape = spectral_indices(sed_data, plot)
        print('done')
        my_shape = sed_shaper(my_shape, plot)
        print('done')
        prefit_jet = model_constr(my_shape, plot)
        print('done')
        jet_lsb, model_minimizer_lsb, fit_model_lsb = model_fit_lsb(
            sed_data, my_shape, plot)
        print('TEST PASSED: OK')
    except Exception as e:
        return RuntimeError('test failed', e)
Exemple #4
0
def test_full(plot=False):
    from jetset.plot_sedfit import plt
    plt.ioff()
    test_jet(plot=plot)
    test_hadronic_jet(plot=plot)
    fit_model, model_minimizer, data = test_model_fit(minimizer='lsb',
                                                      sed_number=1)
    fit_model, model_minimizer, data = test_model_fit(minimizer='minuit',
                                                      sed_number=1)
    fit_model, model_minimizer, data = test_model_fit(minimizer='lsb',
                                                      sed_number=2)
    test_ebl(plot=plot)
    test_ebl_jet(plot=plot)
    test_ebl_jet_fit(plot=plot)
    test_emcee(plot=plot)
    test_dep_par(plot=plot)
    test_dep_par_jet(plot=plot)
    test_dep_par_composite_model(plot=plot)
    test_hadronic_jet(plot=plot)
Exemple #5
0
def test_my_foo():
    from jetset.plot_sedfit import plt
    plt.ioff()
    j = Jet()
    j.eval()
    j.energetic_report()