Ejemplo n.º 1
0
res_tables = []
for regime in jd_tuples:
    print('\nStarting new fitting regime={}\n'.format(regime['label']))

    # Select subset of the dataset.
    handler.select_fits_subset(pre_processing=None,
                               binning=None,
                               exposure=(None, 'max'),
                               reduction=None,
                               observatory=None,
                               jd=regime['jd_rule'])

    # Define fitting routines.
    fitter = SpectralFits(GeminiTelescopeInstruments.gmos)
    fitter.add(
        etacar_templates.h_i_3771(version=regime['template'],
                                  solver=regime['solver']))

    # Ready fitting routine and pre-processing options.
    fitter.rcs_edge_size = (2, 10)
    fitter.compile(helio_correction=False,
                   continuum_normalisation=True,
                   re_bin=None,
                   refine_continuum=True,
                   bad_jds=fetch_bad_jds(db_path=DB_PATH,
                                         fit='h_iota_master',
                                         comp='Any'))

    # Execute fitting routines.
    fitter.fit(handler.data_subset,
               diagnostics=False,
Ejemplo n.º 2
0
res_tables = []
for regime in jd_tuples:
    print('\nStarting new fitting regime={}\n'.format(regime['label']))

    # Select subset of the dataset.
    handler.select_fits_subset(pre_processing=None,
                               binning=None,
                               exposure=(None, 'max'),
                               reduction='endeavour',
                               observatory=regime['label'],
                               jd=regime['jd_rule'])

    # Define fitting routines.
    fitter = SpectralFits(regime['telescope'])
    fitter.add(
        novacar_templates.h_i_6563(version=regime['template'],
                                   solver=regime['solver']))

    # Ready fitting routine and pre-processing options.
    fitter.compile(helio_correction=False,
                   continuum_normalisation=True,
                   re_bin=None,
                   refine_continuum=True,
                   bad_jds=fetch_bad_jds())

    # Execute fitting routines.
    fitter.fit(handler.data_subset,
               diagnostics=False,
               draw=False,
               db_path=DB_PATH)