Create CCSNe database¶
This notebook creates a database file containing 18 different CCSNe models from 6 different datasets.
Running this notebook requires the user has the raw data from the different models on their computer. The raw data is not available as part of SIMPLE. It may be available upon request from the original authors.
The database file this notebook generates contains the abundance for each mass coordinate in the CCSNe that is ejected. See the CCSNe tutorial on how to use this database file.
List of models
La22Lawson et al. (2021)LC18Limongi & Chieffi (2018)Pi16Pignatari et al. (2016)Ra02Rauscher et al. (2002)Ri18Ritter et al. (2018)Si18Sieverding et al. (2018)
In [1]:
Copied!
import simple
from simple import models, ccsne
import logging
simple.set_logging_level('INFO')
logger = logging.getLogger('SIMPLE.CreateDataBase')
modelcollection = simple.new_collection()
datadir = 'data/'
savefilename = 'SIMPLE_CCSNe_v3.hdf5'
#########################
### Import ref item ###
#########################
name = 'W17'
logger.info(f'Generating Ref: {name}')
values, keys = models.load_csv_h(f'{datadir}ref/isotope_mass_W17.csv')
model = modelcollection.new_ref('IsoRef', name,
type='ref_isomass', citation='',
data_values=values, data_keys=keys, data_unit='Da')
name = 'iniab2.0E-02GN93'
logger.info(f'Generating Ref: {name}')
values, keys = models.load_ppn(f'{datadir}ref/{name}.ppn')
model = modelcollection.new_ref('IsoRef', name,
type='ref_isoabu', citation='',
data_values=values, data_keys=keys,
data_unit='mass')
name = 'iniab2.0E-02GN93_Lo03'
logger.info(f'Generating Ref: {name}')
values, keys = models.load_ppn(f'{datadir}ref/{name}.ppn')
model = modelcollection.new_ref('IsoRef', name,
type='ref_isoabu', citation='',
data_values=values, data_keys=keys,
data_unit='mass')
name = 'iniab1.3E-02Lo03_Lo03_merged'
logger.info(f'Generating Ref: {name}')
values, keys = models.load_ppn(f'{datadir}ref/{name}.ppn')
model = modelcollection.new_ref('IsoRef', name,
type='ref_isoabu', citation='',
data_values=values, data_keys=keys,
data_unit='mass')
name = 'rau_solar_ref'
logger.info(f'Generating Ref: {name}')
values, keys = models.load_ppn(f'{datadir}ref/{name}.ppn')
model = modelcollection.new_ref('IsoRef', name,
type='ref_isoabu', citation='',
data_values=values, data_keys=keys,
data_unit='mass')
name = 'iniab1.3E-02As09_As09_merged'
logger.info(f'Generating Ref: {name}')
values, keys = models.load_ppn(f'{datadir}ref/{name}.ppn')
model = modelcollection.new_ref('IsoRef', name,
type='ref_isoabu', citation='',
data_values=values, data_keys=keys,
data_unit='mass')
###########################
### Import CCSNe models ###
###########################
logger.info(f'Generating CCSNe dataset: Ri18')
models = ccsne.load_Ri18(datadir+'R18/', 'iniab2.0E-02GN93', 'W17')
for name, attrs in models.items(): model = modelcollection.new_model('CCSNe', name, **attrs)
logger.info(f'Generating CCSNe dataset: Pi16')
models = ccsne.load_Pi16(datadir+'P16/', 'iniab2.0E-02GN93_Lo03', 'W17')
for name, attrs in models.items(): model = modelcollection.new_model('CCSNe', name, **attrs)
logger.info(f'Generating CCSNe dataset: La22')
models = ccsne.load_La22(datadir+'LAW22/', 'iniab2.0E-02GN93_Lo03', 'W17', remove_Mrem=True)
for name, attrs in models.items(): model = modelcollection.new_model('CCSNe', name, **attrs)
logger.info(f'Generating CCSNe dataset: Si18')
models = ccsne.load_Si18(datadir+'SIE18/', 'iniab1.3E-02Lo03_Lo03_merged', 'W17')
for name, attrs in models.items(): model = modelcollection.new_model('CCSNe', name, **attrs)
logger.info(f'Generating CCSNe dataset: Ra02')
models = ccsne.load_Ra02(datadir+'R02/', 'rau_solar_ref', 'W17')
for name, attrs in models.items(): model = modelcollection.new_model('CCSNe', name, **attrs)
logger.info(f'Generating CCSNe dataset: LC18')
models = ccsne.load_LC18(datadir+'LC18/', 'iniab1.3E-02As09_As09_merged', 'W17')
for name, attrs in models.items(): model = modelcollection.new_model('CCSNe', name, **attrs)
#################
### Save file ###
#################
modelcollection.save(f'{datadir}{savefilename}')
import simple
from simple import models, ccsne
import logging
simple.set_logging_level('INFO')
logger = logging.getLogger('SIMPLE.CreateDataBase')
modelcollection = simple.new_collection()
datadir = 'data/'
savefilename = 'SIMPLE_CCSNe_v3.hdf5'
#########################
### Import ref item ###
#########################
name = 'W17'
logger.info(f'Generating Ref: {name}')
values, keys = models.load_csv_h(f'{datadir}ref/isotope_mass_W17.csv')
model = modelcollection.new_ref('IsoRef', name,
type='ref_isomass', citation='',
data_values=values, data_keys=keys, data_unit='Da')
name = 'iniab2.0E-02GN93'
logger.info(f'Generating Ref: {name}')
values, keys = models.load_ppn(f'{datadir}ref/{name}.ppn')
model = modelcollection.new_ref('IsoRef', name,
type='ref_isoabu', citation='',
data_values=values, data_keys=keys,
data_unit='mass')
name = 'iniab2.0E-02GN93_Lo03'
logger.info(f'Generating Ref: {name}')
values, keys = models.load_ppn(f'{datadir}ref/{name}.ppn')
model = modelcollection.new_ref('IsoRef', name,
type='ref_isoabu', citation='',
data_values=values, data_keys=keys,
data_unit='mass')
name = 'iniab1.3E-02Lo03_Lo03_merged'
logger.info(f'Generating Ref: {name}')
values, keys = models.load_ppn(f'{datadir}ref/{name}.ppn')
model = modelcollection.new_ref('IsoRef', name,
type='ref_isoabu', citation='',
data_values=values, data_keys=keys,
data_unit='mass')
name = 'rau_solar_ref'
logger.info(f'Generating Ref: {name}')
values, keys = models.load_ppn(f'{datadir}ref/{name}.ppn')
model = modelcollection.new_ref('IsoRef', name,
type='ref_isoabu', citation='',
data_values=values, data_keys=keys,
data_unit='mass')
name = 'iniab1.3E-02As09_As09_merged'
logger.info(f'Generating Ref: {name}')
values, keys = models.load_ppn(f'{datadir}ref/{name}.ppn')
model = modelcollection.new_ref('IsoRef', name,
type='ref_isoabu', citation='',
data_values=values, data_keys=keys,
data_unit='mass')
###########################
### Import CCSNe models ###
###########################
logger.info(f'Generating CCSNe dataset: Ri18')
models = ccsne.load_Ri18(datadir+'R18/', 'iniab2.0E-02GN93', 'W17')
for name, attrs in models.items(): model = modelcollection.new_model('CCSNe', name, **attrs)
logger.info(f'Generating CCSNe dataset: Pi16')
models = ccsne.load_Pi16(datadir+'P16/', 'iniab2.0E-02GN93_Lo03', 'W17')
for name, attrs in models.items(): model = modelcollection.new_model('CCSNe', name, **attrs)
logger.info(f'Generating CCSNe dataset: La22')
models = ccsne.load_La22(datadir+'LAW22/', 'iniab2.0E-02GN93_Lo03', 'W17', remove_Mrem=True)
for name, attrs in models.items(): model = modelcollection.new_model('CCSNe', name, **attrs)
logger.info(f'Generating CCSNe dataset: Si18')
models = ccsne.load_Si18(datadir+'SIE18/', 'iniab1.3E-02Lo03_Lo03_merged', 'W17')
for name, attrs in models.items(): model = modelcollection.new_model('CCSNe', name, **attrs)
logger.info(f'Generating CCSNe dataset: Ra02')
models = ccsne.load_Ra02(datadir+'R02/', 'rau_solar_ref', 'W17')
for name, attrs in models.items(): model = modelcollection.new_model('CCSNe', name, **attrs)
logger.info(f'Generating CCSNe dataset: LC18')
models = ccsne.load_LC18(datadir+'LC18/', 'iniab1.3E-02As09_As09_merged', 'W17')
for name, attrs in models.items(): model = modelcollection.new_model('CCSNe', name, **attrs)
#################
### Save file ###
#################
modelcollection.save(f'{datadir}{savefilename}')
SIMPLE.CreateDataBase - INFO - Generating Ref: W17 SIMPLE.CreateDataBase - INFO - Generating Ref: iniab2.0E-02GN93 SIMPLE.CreateDataBase - INFO - Generating Ref: iniab2.0E-02GN93_Lo03 SIMPLE.CreateDataBase - INFO - Generating Ref: iniab1.3E-02Lo03_Lo03_merged SIMPLE.CreateDataBase - INFO - Generating Ref: rau_solar_ref SIMPLE.CreateDataBase - INFO - Generating Ref: iniab1.3E-02As09_As09_merged SIMPLE.CreateDataBase - INFO - Generating CCSNe dataset: Ri18 SIMPLE.ccsne - INFO - Fudging the masscoord mass from masscoord SIMPLE.ccsne - INFO - Calculating Default Onion Structure SIMPLE.ccsne - INFO - Fudging the masscoord mass from masscoord SIMPLE.ccsne - INFO - Calculating Default Onion Structure SIMPLE.ccsne - INFO - Fudging the masscoord mass from masscoord SIMPLE.ccsne - INFO - Calculating Default Onion Structure SIMPLE.CreateDataBase - INFO - Generating CCSNe dataset: Pi16 SIMPLE.ccsne - INFO - Fudging the masscoord mass from masscoord SIMPLE.ccsne - INFO - Calculating Default Onion Structure SIMPLE.ccsne - INFO - Fudging the masscoord mass from masscoord SIMPLE.ccsne - INFO - Calculating Default Onion Structure SIMPLE.ccsne - INFO - Fudging the masscoord mass from masscoord SIMPLE.ccsne - INFO - Calculating Default Onion Structure SIMPLE.CreateDataBase - INFO - Generating CCSNe dataset: La22 SIMPLE.ccsne - INFO - Fudging the masscoord mass from masscoord SIMPLE.ccsne - INFO - Calculating Default Onion Structure SIMPLE.ccsne - INFO - Removing datapoints from the Mrem zone SIMPLE.ccsne - INFO - Fudging the masscoord mass from masscoord SIMPLE.ccsne - INFO - Calculating Default Onion Structure SIMPLE.ccsne - INFO - Fudging the masscoord mass from masscoord SIMPLE.ccsne - INFO - Calculating Default Onion Structure SIMPLE.ccsne - INFO - Removing datapoints from the Mrem zone SIMPLE.CreateDataBase - INFO - Generating CCSNe dataset: Si18 SIMPLE.ccsne - INFO - Fudging the masscoord mass from masscoord SIMPLE.ccsne - INFO - Calculating Default Onion Structure SIMPLE.ccsne - INFO - Fudging the masscoord mass from masscoord SIMPLE.ccsne - INFO - Calculating Default Onion Structure SIMPLE.ccsne - INFO - Fudging the masscoord mass from masscoord SIMPLE.ccsne - INFO - Calculating Default Onion Structure SIMPLE.CreateDataBase - INFO - Generating CCSNe dataset: Ra02 SIMPLE.ccsne - INFO - Calculating Default Onion Structure SIMPLE.ccsne - INFO - Calculating Default Onion Structure SIMPLE.ccsne - INFO - Calculating Default Onion Structure SIMPLE.CreateDataBase - INFO - Generating CCSNe dataset: LC18 SIMPLE.ccsne - INFO - Fudging the masscoord mass from masscoord SIMPLE.ccsne - INFO - Calculating Default Onion Structure SIMPLE.ccsne - INFO - Fudging the masscoord mass from masscoord SIMPLE.ccsne - INFO - Calculating Default Onion Structure SIMPLE.ccsne - INFO - Fudging the masscoord mass from masscoord SIMPLE.ccsne - INFO - Calculating Default Onion Structure SIMPLE.models - INFO - Saving ModelCollection(v3.0) as: data/SIMPLE_CCSNe_v3.hdf5 SIMPLE.models - INFO - saving ref: W17 SIMPLE.models - INFO - saving ref: iniab2.0E-02GN93 SIMPLE.models - INFO - saving ref: iniab2.0E-02GN93_Lo03 SIMPLE.models - INFO - saving ref: iniab1.3E-02Lo03_Lo03_merged SIMPLE.models - INFO - saving ref: rau_solar_ref SIMPLE.models - INFO - saving ref: iniab1.3E-02As09_As09_merged SIMPLE.models - INFO - saving model: Ri18_m15 SIMPLE.models - INFO - saving model: Ri18_m20 SIMPLE.models - INFO - saving model: Ri18_m25 SIMPLE.models - INFO - saving model: Pi16_m15 SIMPLE.models - INFO - saving model: Pi16_m20 SIMPLE.models - INFO - saving model: Pi16_m25 SIMPLE.models - INFO - saving model: La22_m15 SIMPLE.models - INFO - saving model: La22_m20 SIMPLE.models - INFO - saving model: La22_m25 SIMPLE.models - INFO - saving model: Si18_m15 SIMPLE.models - INFO - saving model: Si18_m20 SIMPLE.models - INFO - saving model: Si18_m25 SIMPLE.models - INFO - saving model: Ra02_m15 SIMPLE.models - INFO - saving model: Ra02_m20 SIMPLE.models - INFO - saving model: Ra02_m25 SIMPLE.models - INFO - saving model: LC18_m15 SIMPLE.models - INFO - saving model: LC18_m20 SIMPLE.models - INFO - saving model: LC18_m25 SIMPLE.models - INFO - Filesize: 95.79 MB SIMPLE.models - INFO - Time to save file: 0:00:08.933632
In [2]:
Copied!
modelcollection
modelcollection
Out[2]:
Models in collection:
- [0]
Ri18_m15(CCSNe) - [1]
Ri18_m20(CCSNe) - [2]
Ri18_m25(CCSNe) - [3]
Pi16_m15(CCSNe) - [4]
Pi16_m20(CCSNe) - [5]
Pi16_m25(CCSNe) - [6]
La22_m15(CCSNe) - [7]
La22_m20(CCSNe) - [8]
La22_m25(CCSNe) - [9]
Si18_m15(CCSNe) - [10]
Si18_m20(CCSNe) - [11]
Si18_m25(CCSNe) - [12]
Ra02_m15(CCSNe) - [13]
Ra02_m20(CCSNe) - [14]
Ra02_m25(CCSNe) - [15]
LC18_m15(CCSNe) - [16]
LC18_m20(CCSNe) - [17]
LC18_m25(CCSNe)
References in collection:
W17(IsoRef)iniab2.0E-02GN93(IsoRef)iniab2.0E-02GN93_Lo03(IsoRef)iniab1.3E-02Lo03_Lo03_merged(IsoRef)rau_solar_ref(IsoRef)iniab1.3E-02As09_As09_merged(IsoRef)
In [3]:
Copied!
modelcollection[6]
modelcollection[6]
Out[3]:
Name: La22_m15
Type: CCSNe
Dataset: La22
Mass: 15
Attributes: name, type, dataset, citation, refid_isoabu, refid_isomass, mass, masscoord, masscoord_mass, abundance_values, abundance_keys, abundance_unit, zone, clsname, abundance, ref_isoabu, ref_isomass