"Mo", "Tc", "Mo", "Ru", "Rh", "Pd", "Ag", "Cd", "Sn", "Te",
"Ba", "La", "Ce", "Pr", "Nd", "Sm", "Eu", "Gd", "Tb", "Dy",
"Ho", "Er", "Tm", "Yb", "Lu", "Hf", "W", "Os", "Ir", "Pt",
"Au", "Hg", "Pb", "Bi", "Th", "U")
# Start logging our progress
logging.basicConfig(
level=logging.INFO,
format='[%(asctime)s] %(levelname)s:%(filename)s:%(message)s',
datefmt='%d/%m/%Y %H:%M:%S')
logger = logging.getLogger(__name__)
logger.info("Synthesizing SAGA spectra")
# Instantiate base synthesizer
synthesizer = Synthesizer(library_name="saga_sample",
logger=logger,
docstring=__doc__)
def valid_float(n):
try:
test = re.match("<?(.*)+-.*", n)
if test is not None:
n = test.group(1)
n = float(n)
except ValueError:
return None
except TypeError:
return None
return n
# List of elements whose abundances we pass to TurboSpectrum
element_list = ("[C/H]", "[N/H]", "[O/H]", "[Na/H]", "[Mg/H]", "[Al/H]",
"[Si/H]", "[Ca/H]", "[Ti/H]", "[Mn/H]", "[Co/H]", "[Ni/H]",
"[Ba/H]", "[Sr/H]")
# Start logging our progress
logging.basicConfig(
level=logging.INFO,
format='[%(asctime)s] %(levelname)s:%(filename)s:%(message)s',
datefmt='%d/%m/%Y %H:%M:%S')
logger = logging.getLogger(__name__)
logger.info("Synthesizing APOKASC grid of spectra")
# Instantiate base synthesizer
synthesizer = Synthesizer(library_name="apokasc_training_set",
logger=logger,
docstring=__doc__)
# Table supplies list of stars in the APOKASC training set, giving the stellar labels for each star in the training set
stellar_data = Table.read("../../4MOST_testspectra/trainingset_param.tab",
format="ascii")
# Iterate over the spectra we're supposed to be synthesizing
star_list = []
for star in stellar_data:
# Look up stellar parameters of the star we're about to synthesize
metadata = synthesizer.astropy_row_to_dict(star)
star_data = {
"name": metadata["Starname"],
"stellar_mass": metadata['Mass'] if 'Mass' in metadata else 1,
"""
import logging
from lib.base_synthesizer import Synthesizer
from fourgp_speclib import SpectrumLibrarySqlite
# Start logging our progress
logging.basicConfig(level=logging.INFO, format='[%(asctime)s] %(levelname)s:%(filename)s:%(message)s',
datefmt='%d/%m/%Y %H:%M:%S')
logger = logging.getLogger(__name__)
logger.info("Synthesizing spectra of pepsi")
# Instantiate base synthesizer
synthesizer = Synthesizer(library_name="turbo_pepsi_replica_3label",
logger=logger,
docstring=__doc__)
spectra = SpectrumLibrarySqlite.open_and_search(
library_spec='pepsi_4fs_hrs/',
workspace='/home/travegre/Projects/4GP/4most-4gp-scripts/workspace/',
extra_constraints={"continuum_normalised": True}
)
pepsi_library, pepsi_library_items = [spectra[i] for i in ("library", "items")]
# Load test set
pepsi_library_ids = [i["specId"] for i in pepsi_library_items]
spectra = [pepsi_library.open(ids=i).extract_item(0) for i in pepsi_library_ids]
print(spectra)
star_list = []
import logging
from lib.base_synthesizer import Synthesizer
# Start logging our progress
logging.basicConfig(
level=logging.INFO,
format='[%(asctime)s] %(levelname)s:%(filename)s:%(message)s',
datefmt='%d/%m/%Y %H:%M:%S')
logger = logging.getLogger(__name__)
logger.info(
"Synthesizing spectra for some pollution sources for cross-talk tests")
# Instantiate base synthesizer
synthesizer = Synthesizer(library_name="pollution_sources_turbospec",
logger=logger,
docstring=__doc__)
star_list = [
{
'name': "rc,fe_h=-3",
'Teff': 4500,
'logg': 3,
'[Fe/H]': -3,
'extra_metadata': {}
},
{
'name': "to,fe_h=-3",
'Teff': 6250,
'logg': 4,
'[Fe/H]': -3,