Choosing a cosmology model

Choosing a cosmology model#

Since the new implementation allows to switch on the fly from the src to the obs frame, we need to define a cosmolgoy model in order to evaluate the luminosity distance also for dataset independently from the Jet model. This is implemented in usign the astropy cosmology subpackaged https://docs.astropy.org/en/stable/cosmology/ Moreover, this allows also to change easilyt the cosmology used for the Jet model, passing a specifc cosmology model form astropy

import jetset
print('tested on jetset',jetset.__version__)

tested on jetset 1.3.0rc7

from jetset.cosmo_tools import Cosmo
c=Cosmo()

the defualt comoslogy is

c

FlatLambdaCDM(name="Planck13", H0=67.77 km / (Mpc s), Om0=0.30712, Tcmb0=2.7255 K, Neff=3.046, m_nu=[0.   0.   0.06] eV, Ob0=0.048252)

c.get_DL_cm(z=1)

2.095e+28

you can choose a different cosmology model from astropy eg

from astropy.cosmology import WMAP9 as cw
c=Cosmo(astropy_cosmo=cw)

c

FlatLambdaCDM(name="WMAP9", H0=69.32 km / (Mpc s), Om0=0.2865, Tcmb0=2.725 K, Neff=3.04, m_nu=[0. 0. 0.] eV, Ob0=0.04628)

c.get_DL_cm(z=1)

2.075e+28

or you can set directly a luminosity distance in cm

c=Cosmo(DL_cm=1e28)

c

cosmology is not defined, the luminosity distance has been set to 1e+28 cm

Changing cosmology model in data#

from jetset.data_loader import Data
from jetset.data_loader import ObsData
from jetset.cosmo_tools import Cosmo
from jetset.test_data_helper import  test_SEDs
c=Cosmo()
data_table=Data.from_file(test_SEDs[1])
sed_data=ObsData(data_table=data_table,cosmo=c)

Changing cosmology model in Jet models#

at object instantiation

from jetset.jet_model import Jet
my_jet=Jet(cosmo=c)

===> setting C threads to 12

or for already built models

my_jet.cosmo=c
Contents