import os from config import REPO_PATH, DB_PATH from extract.core.lines import SpectralFits from extract.helpers import convert_vacuum_to_air_wavelengths, \ calculate_vacuum_wavelength_rydberg_formula # Pickle combination alg fits. SpectralFits.save(db_path=DB_PATH, table='Hepsilon_RegimeConcat', comp_number=[1, 2, 3, 4], lab_wavelength=convert_vacuum_to_air_wavelengths( calculate_vacuum_wavelength_rydberg_formula(2, 7)), velocity_error='pooled_intra_night_variance', pickle_path=os.path.join( REPO_PATH, 'Pickles', 'EtaCar_v3', 'fits_gmos_hepsilon_combi_[1,2,3,4].p')) # Pickle benchmark bisector half max fits. SpectralFits.save(db_path=DB_PATH, table='Hepsilon_BHM', comp_number=None, lab_wavelength=convert_vacuum_to_air_wavelengths( calculate_vacuum_wavelength_rydberg_formula(2, 7)), velocity_error=10.0, pickle_path=os.path.join(REPO_PATH, 'Pickles', 'EtaCar_v3', 'fits_gmos_hepsilon_bhm.p'))
from config import DB_PATH from extract.core.lines import SpectralFits from extract.observatories.gemini import GeminiTelescopeInstruments from extract.helpers import fetch_bad_jds # Define viewing routines. viewer = SpectralFits(GeminiTelescopeInstruments.gmos) viewer.view(db_path=DB_PATH, table='Hzeta_RegimeConcat', start_at=0, comp_number=[1, 2, 3, 4], bad_jds=fetch_bad_jds( db_path=DB_PATH, fit='h_zeta_master', comp='Any'), fold=False, save=False, save_loc='/Users/grantd/ScreenGrabs/SpectralDynamics/' 'Videos/Hzeta/Hzeta_GMOS_Fit_comp{}_{}.jpg')
# Iterate different fitting regimes. 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_3889(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(db_path=DB_PATH, fit='h_zeta_master', comp='Any')) # Execute fitting routines. fitter.fit(handler.data_subset,
from config import DB_PATH from extract.core.lines import SpectralFits from extract.observatories.gemini import GeminiTelescopeInstruments from extract.helpers import fetch_bad_jds # Define viewing routines. viewer = SpectralFits(GeminiTelescopeInstruments.gmos) viewer.rcs_edge_size = (2, 10) viewer.view(db_path=DB_PATH, table='Hiota_RegimeConcat', start_at=0, comp_number=[1, 2, 3, 4], bad_jds=fetch_bad_jds( db_path=DB_PATH, fit='h_iota_master', comp='Any'), fold=False, save=False, save_loc='/Data/ScreenGrabs/AstroDynamics Screenshots/' 'Videos/Hiota/Hiota_GMOS_Fit_comp{}_{}.jpg')