In [1]:
import pandas as pd
from lsst.sims.photUtils import Bandpass
import sncosmo
from astropy.cosmology import FlatLambdaCDM
In [2]:
dc2_df = pd.read_csv('/Users/Bryce/Documents/DESC/Twinkles/data/dc2_sne_cat.csv')
In [4]:
dc2_df['MB'].head()
Out[4]:
0 -19.534964
1 -19.534964
2 -19.570532
3 -19.570532
4 -19.602254
Name: MB, dtype: float64
In [5]:
#DC2 Cosmology
#FlatLambdaCDM(H0=71 km / (Mpc s), Om0=0.265, Tcmb0=0 K, Neff=3.04, m_nu=None, Ob0=0.045)
In [6]:
source = sncosmo.get_source('salt2-extended')
In [7]:
model = sncosmo.Model(source=source)
In [8]:
model.param_names
Out[8]:
['z', 't0', 'x0', 'x1', 'c']
In [9]:
model.parameters
Out[9]:
array([ 0., 0., 1., 0., 0.])
In [10]:
# Danny's x1
model.set(x1 = 0.)
In [12]:
cosmo = FlatLambdaCDM(H0=71, Om0=0.265, Tcmb0=0, Neff=3.04, m_nu=None, Ob0=0.045)
In [13]:
cosmo
Out[13]:
FlatLambdaCDM(H0=71 km / (Mpc s), Om0=0.265, Tcmb0=0 K, Neff=3.04, m_nu=None, Ob0=0.045)
In [14]:
x0 = []
for sn_row in dc2_df.iterrows():
z = sn_row[1]['zs']
MB = sn_row[1]['MB']
model.set(z=z)
model.set_source_peakabsmag(MB, 'bessellb', 'ab', cosmo=cosmo)
x0.append(model.get('x0'))
In [15]:
x0
Out[15]:
[5.022016731155119e-06,
5.0220167311551029e-06,
9.6968507978990683e-06,
9.6968507978990683e-06,
3.9881024553384965e-06,
3.9881024553384965e-06,
2.9517452686003477e-06,
2.9517452686003477e-06,
8.4912574036901685e-06,
8.4912574036901414e-06,
1.3060588891085317e-05,
1.3060588891085275e-05,
4.2019579443778416e-06,
4.201957944377856e-06,
2.4870983940948158e-06,
2.4870983940948158e-06,
5.5856155427907436e-06,
5.5856155427907258e-06,
4.4988821081651216e-06,
4.4988821081651216e-06,
4.6018726367863331e-06,
4.6018726367863483e-06,
4.3424516256507524e-05,
4.3424516256507666e-05,
5.9546210533973505e-06,
5.954621053397331e-06,
6.6542725979596752e-06,
6.6542725979596752e-06,
3.5383320893094676e-06,
3.5383320893094795e-06,
4.181990160993982e-06,
4.1819901609939684e-06,
6.2031401092717371e-06,
6.2031401092717371e-06,
6.2031401092717371e-06,
2.0386548344483967e-06,
2.0386548344484035e-06,
7.4587333440154906e-06,
7.4587333440155152e-06,
2.4085723149986547e-06,
2.408572314998647e-06,
7.8569833279442584e-06,
7.8569833279442839e-06,
3.4308635687124276e-06,
3.4308635687124276e-06,
3.4308635687124276e-06,
3.4308635687124276e-06,
4.3494407241230156e-06,
4.34944072412303e-06,
4.9572756120309853e-06,
4.9572756120310014e-06,
5.2457038471454818e-06,
5.2457038471454818e-06,
3.4661636017808191e-06,
3.4661636017808191e-06,
3.4661636017808191e-06,
3.4661636017808191e-06,
3.9479015197410726e-06,
3.9479015197410726e-06,
1.2177390426442126e-05,
1.2177390426442126e-05,
5.5262872420783548e-06,
5.5262872420783548e-06,
2.6629464703828036e-06,
2.6629464703828036e-06,
3.6824815791927196e-06,
3.6824815791927196e-06,
5.8881822872196878e-06,
5.8881822872196878e-06,
1.3583144953085215e-05,
1.3583144953085215e-05,
1.3583144953085215e-05,
1.3583144953085215e-05,
7.4563057952679826e-06,
7.4563057952680072e-06,
6.7347818959626931e-06,
6.7347818959626931e-06,
8.4086264498113991e-06,
8.4086264498114279e-06,
3.5521934652567865e-06,
3.5521934652567983e-06,
2.0647263587776478e-06,
2.0647263587776478e-06,
5.2607387125669116e-06,
5.2607387125669116e-06,
5.2607387125669116e-06,
5.2607387125669116e-06,
7.4844370783101359e-06,
7.4844370783101359e-06,
3.3147158353634455e-06,
3.3147158353634455e-06,
3.0537487721375022e-06,
3.0537487721375124e-06,
3.0537487721375026e-06,
3.0537487721375128e-06,
2.7403156784576762e-06,
2.7403156784576762e-06,
2.7403156784576762e-06,
1.2202547544997149e-05,
1.2202547544997149e-05,
5.52892837808651e-06,
5.52892837808651e-06,
7.8444920277072496e-06,
7.8444920277072496e-06,
3.8882466899964485e-06,
3.8882466899964612e-06,
2.5890686562445894e-06,
2.5890686562445894e-06,
3.7325395324735301e-06,
3.7325395324735301e-06,
6.3748602624039295e-06,
6.3748602624039295e-06,
2.4513526049577837e-06,
2.4513526049577837e-06,
3.4177215975699314e-06,
3.4177215975699314e-06,
2.6457237135368426e-06,
2.6457237135368341e-06,
3.1050403498946809e-05,
3.1050403498946809e-05,
4.3117154669124262e-06,
4.3117154669124262e-06,
4.3117154669124262e-06,
4.3117154669124262e-06,
1.2210672315521821e-05,
1.2210672315521862e-05,
4.045179739106224e-06,
4.0451797391062375e-06,
4.0451797391062375e-06,
4.0451797391062375e-06,
4.0091071012569166e-06,
4.0091071012569039e-06,
3.026627419555696e-06,
3.0266274195556862e-06,
6.706630899956195e-06,
6.706630899956195e-06,
1.4695660609260097e-05,
1.4695660609260097e-05,
4.2018945378489357e-06,
4.2018945378489357e-06,
1.0717860750788253e-05,
1.0717860750788288e-05,
1.0717860750788254e-05,
1.071786075078829e-05,
2.0657626823004806e-05,
2.0657626823004806e-05,
8.4374269759417456e-06,
8.4374269759417744e-06,
2.0376266894729319e-06,
2.0376266894729319e-06,
2.4844119032282461e-06,
2.4844119032282546e-06,
6.588074262646423e-06,
6.5880742626464019e-06,
5.5525638486153679e-06,
5.5525638486153501e-06,
3.065356244076975e-06,
3.0653562440769652e-06,
5.0158708818155609e-06,
5.0158708818155609e-06,
3.1311671596279417e-06,
3.1311671596279315e-06,
2.5422322012335912e-06,
2.5422322012335996e-06,
7.9828837568911161e-06,
7.9828837568911161e-06,
2.4182504022461275e-06,
2.4182504022461275e-06,
1.0316925660120619e-05,
1.0316925660120619e-05,
4.2738383624932329e-06,
4.2738383624932329e-06,
1.0741152128907432e-05,
1.0741152128907432e-05,
6.4503559662619676e-06,
6.4503559662619676e-06,
5.9551130270719985e-06,
5.955113027071979e-06,
4.3043617888631991e-06,
4.3043617888631991e-06,
4.4626618074752097e-06,
4.4626618074752097e-06,
3.1431005291473502e-05,
3.1431005291473502e-05,
6.7037390722301754e-06,
6.7037390722301754e-06,
4.0365303262909257e-06,
4.0365303262909257e-06,
1.2045003784519826e-05,
1.2045003784519826e-05,
2.7587270040142531e-06,
2.7587270040142531e-06,
2.7587270040142531e-06,
2.7587270040142531e-06,
1.8338597109632443e-05,
1.8338597109632443e-05,
8.4065012378056275e-06,
8.4065012378056275e-06,
8.9040249794710044e-06,
8.9040249794709756e-06,
2.0813026290196211e-06,
2.0813026290196279e-06,
4.7938064768214855e-06,
4.7938064768214855e-06,
2.1619096382639186e-06,
2.1619096382639186e-06,
4.5693244500640831e-06,
4.5693244500640831e-06,
2.0400682816319715e-06,
2.0400682816319783e-06,
4.1086332615455279e-05,
4.1086332615455143e-05,
6.8302446279457543e-06,
6.8302446279457322e-06,
1.1811232191803063e-05,
1.1811232191803063e-05,
5.187339181245427e-06,
5.18733918124541e-06,
6.2951246666651178e-06,
6.2951246666651178e-06,
2.397602773282556e-05,
2.397602773282556e-05,
2.397602773282556e-05,
2.397602773282556e-05,
6.1507700683411881e-06,
6.1507700683411881e-06,
2.4740036859294057e-06,
2.4740036859294057e-06,
8.561235509276238e-06,
8.5612355092762668e-06,
6.1045772959770469e-06,
6.1045772959770274e-06,
5.418372186598501e-06,
5.4183721865985188e-06,
5.812528292350473e-06,
5.8125282923504543e-06,
3.8818130774934176e-06,
3.8818130774934176e-06,
1.1361819279613686e-05,
1.1361819279613686e-05,
4.5516912515974357e-06,
4.5516912515974357e-06,
1.4380348444416505e-05,
1.4380348444416552e-05,
1.4380348444416507e-05,
1.4380348444416554e-05,
8.3731504269620118e-06,
8.3731504269620118e-06,
5.8417265247125974e-06,
5.8417265247126168e-06,
5.6003687981536411e-06,
5.6003687981536411e-06,
4.5849560049982565e-06,
4.5849560049982421e-06,
1.6332834138640842e-05,
1.6332834138640788e-05,
3.5017396094838331e-06,
3.5017396094838217e-06,
3.0857333550263858e-06,
3.0857333550263858e-06,
1.6781192887075853e-05,
1.6781192887075853e-05,
1.263077105297761e-05,
1.2630771052977652e-05,
4.0775953454341672e-06,
4.0775953454341545e-06,
5.524707451574162e-06,
5.5247074515741442e-06,
8.187265389004685e-06,
8.187265389004685e-06,
3.6881461594376561e-06,
3.6881461594376561e-06,
6.131421743972388e-06,
6.1314217439724083e-06,
5.7011461339351988e-06,
5.7011461339351802e-06,
2.437959165404278e-05,
2.437959165404278e-05,
2.6904532381807461e-06,
2.6904532381807461e-06,
2.6904532381807461e-06,
2.6904532381807461e-06,
4.5240019547775452e-06,
4.5240019547775452e-06,
4.5987974721611396e-06,
4.5987974721611252e-06,
4.5786592267916988e-06,
4.5786592267916988e-06,
1.3691550871675002e-05,
1.3691550871675002e-05,
3.6528675188206681e-06,
3.6528675188206681e-06,
3.6528675188206681e-06,
3.6528675188206681e-06,
1.1442313214094053e-05,
1.1442313214094016e-05,
2.1084268272840856e-06,
2.1084268272840856e-06,
4.604667572993871e-06,
4.604667572993871e-06,
3.541878754302482e-06,
3.541878754302482e-06,
2.4712983295995457e-06,
2.4712983295995457e-06,
3.2160752442867211e-06,
3.2160752442867211e-06,
4.7712372173264702e-06,
4.7712372173264702e-06,
3.0498638056598993e-06,
3.0498638056598993e-06,
4.9239585984817311e-06,
4.9239585984817311e-06,
4.8675331336438833e-06,
4.8675331336438833e-06,
2.2842116895015913e-05,
2.2842116895015988e-05,
5.2028385462794186e-06,
5.2028385462794356e-06,
2.4484305175616957e-06,
2.4484305175616957e-06,
2.4484305175616957e-06,
2.4484305175616957e-06,
3.805122688590902e-06,
3.805122688590902e-06,
6.2039194083724685e-06,
6.2039194083724685e-06,
7.859474638509711e-06,
7.859474638509711e-06,
2.3563173015085562e-06,
2.3563173015085562e-06,
4.6487826874205356e-06,
4.6487826874205509e-06,
4.920835540026979e-06,
4.920835540026979e-06,
3.1002188074331109e-06,
3.1002188074331109e-06,
1.1571611705965168e-05,
1.1571611705965168e-05,
1.1571611705965168e-05,
1.1571611705965168e-05,
5.7539918181985429e-06,
5.7539918181985429e-06,
2.6581806760784891e-06,
2.6581806760784891e-06,
1.0760945706238462e-05,
1.0760945706238462e-05,
5.2998031115939831e-06,
5.2998031115939831e-06,
3.0887665719483977e-06,
3.0887665719484079e-06,
4.2606858369615226e-06,
4.2606858369615226e-06,
4.2606858369615226e-06,
4.2606858369615226e-06,
4.1430951043065601e-06,
4.1430951043065601e-06,
5.7543495133236009e-06,
5.7543495133236009e-06,
5.6331602528101404e-06,
5.6331602528101404e-06,
1.7342408858244875e-05,
1.7342408858244875e-05,
5.1908985267619736e-06,
5.1908985267619905e-06,
4.3488034570528234e-06,
4.3488034570528234e-06,
1.8651746973416981e-06,
1.8651746973416981e-06,
3.672253715936682e-06,
3.672253715936682e-06,
2.5208032544039278e-05,
2.5208032544039363e-05,
5.2005783797670513e-06,
5.2005783797670343e-06,
3.3884363328508701e-06,
3.3884363328508701e-06,
5.6431056135758699e-06,
5.6431056135758885e-06,
6.3388688819680414e-06,
6.3388688819680626e-06,
6.3388688819680626e-06,
6.3388688819680626e-06,
3.3982376804483831e-06,
3.398237680448372e-06,
7.1668914309445455e-06,
7.1668914309445692e-06,
2.7042349795473491e-06,
2.7042349795473491e-06,
3.0890260207922095e-06,
3.0890260207922095e-06,
3.0890260207922095e-06,
3.0890260207922095e-06,
3.2936847229500861e-06,
3.2936847229500971e-06,
2.6936349250843682e-06,
2.6936349250843682e-06,
6.4948310002373927e-06,
6.4948310002373715e-06,
6.4948310002373935e-06,
6.4948310002373724e-06,
2.1785188972625627e-06,
2.1785188972625555e-06,
2.3435448098010919e-06,
2.3435448098010919e-06,
2.3435448098010919e-06,
2.3435448098010919e-06,
1.1094041836462614e-05,
1.1094041836462652e-05,
3.1073196832777517e-06,
3.1073196832777618e-06,
4.7122161511763364e-06,
4.7122161511763364e-06,
1.1912997409969216e-05,
1.1912997409969177e-05,
2.9888604561338583e-06,
2.9888604561338685e-06,
2.9888604561338685e-06,
2.9888604561338685e-06,
3.8954333164825489e-06,
3.8954333164825489e-06,
6.980575127899503e-06,
6.980575127899503e-06,
2.5201174493050669e-06,
2.5201174493050588e-06,
2.9468316076886584e-06,
2.9468316076886584e-06,
3.9796943332087164e-06,
3.9796943332087164e-06,
4.0641106554163235e-06,
4.064110655416337e-06,
1.017201263382257e-05,
1.017201263382257e-05,
4.4757880792995693e-06,
4.4757880792995693e-06,
7.0217952804828419e-06,
7.0217952804828419e-06,
2.7585035064061759e-05,
2.7585035064061759e-05,
2.7468068786076194e-06,
2.7468068786076194e-06,
3.0624736015771189e-06,
3.0624736015771189e-06,
3.0624736015771189e-06,
3.0624736015771189e-06,
4.2877762273073198e-06,
4.2877762273073198e-06,
1.3966685844165795e-05,
1.3966685844165795e-05,
5.8995765176332076e-06,
5.8995765176332076e-06,
6.1454353516261332e-06,
6.1454353516261332e-06,
5.5201703708033082e-06,
5.5201703708032905e-06,
1.2743649024067882e-05,
1.2743649024067882e-05,
2.7622382343385722e-06,
2.7622382343385722e-06,
3.7533931974289334e-06,
3.7533931974289334e-06,
1.5306241622610695e-05,
1.5306241622610644e-05,
8.5259377857258447e-06,
8.5259377857258447e-06,
2.0374111084166058e-06,
2.0374111084166058e-06,
3.1182886073897005e-06,
3.1182886073897111e-06,
4.9656023261123705e-06,
4.9656023261123544e-06,
4.1402436342911069e-06,
4.1402436342911069e-06,
4.1402436342911069e-06,
4.1402436342911069e-06,
1.5376840189918021e-06,
1.5376840189918021e-06,
1.5376840189918021e-06,
1.5376840189918021e-06,
2.9966735091331917e-06,
2.9966735091331917e-06,
8.6328484778784868e-06,
8.6328484778784868e-06,
3.0604192526118529e-06,
3.0604192526118529e-06,
7.1709565788631362e-06,
7.1709565788631599e-06,
3.6922185085676639e-06,
3.692218508567652e-06,
1.1020334819007432e-05,
1.1020334819007432e-05,
2.0016042421086667e-06,
2.0016042421086667e-06,
4.4232867552549383e-06,
4.4232867552549239e-06,
5.2836140006335989e-06,
5.2836140006336167e-06,
1.1853004749578995e-05,
1.1853004749578995e-05,
6.5732047561319609e-06,
6.5732047561319609e-06,
5.9155012575787651e-06,
5.9155012575787651e-06,
2.9694489456560951e-06,
2.9694489456560853e-06,
1.0972525285532315e-05,
1.0972525285532352e-05,
1.0972525285532352e-05,
1.0972525285532352e-05,
3.1713547023865003e-06,
3.1713547023865003e-06,
2.7165951386654565e-06,
2.7165951386654565e-06,
2.7165951386654565e-06,
2.7165951386654565e-06,
3.1486481583619144e-05,
3.1486481583619144e-05,
1.7219339116384607e-06,
1.7219339116384664e-06,
3.404882358910482e-06,
3.404882358910471e-06,
1.2231314885653497e-05,
1.2231314885653537e-05,
2.3892970864016528e-06,
2.3892970864016452e-06,
6.75020701313782e-06,
6.7502070131377759e-06,
4.6040553643151166e-06,
4.6040553643151013e-06,
3.5936997015482976e-06,
3.5936997015482976e-06,
1.1461488540107763e-05,
1.1461488540107763e-05,
3.8306185861074729e-06,
3.8306185861074729e-06,
3.8927875333266242e-06,
3.8927875333266242e-06,
3.8927875333266242e-06,
3.8927875333266242e-06,
5.0780340226894953e-06,
5.0780340226894953e-06,
2.3436157357573897e-06,
2.3436157357573897e-06,
5.5189823069558284e-06,
5.5189823069558284e-06,
3.9642608121409904e-06,
3.9642608121409904e-06,
3.9642608121409904e-06,
3.9642608121409904e-06,
4.2029751428517103e-06,
4.2029751428517103e-06,
5.4932017153671445e-06,
5.4932017153671445e-06,
4.5048651102772048e-06,
4.5048651102772048e-06,
2.7362652721320154e-06,
2.7362652721320154e-06,
9.16437206720232e-06,
9.16437206720232e-06,
3.5220217714750667e-06,
3.5220217714750667e-06,
8.7073156271316794e-06,
8.7073156271316794e-06,
1.1506668151455907e-05,
1.1506668151455907e-05,
3.0250852858638723e-06,
3.0250852858638723e-06,
3.0250852858638723e-06,
3.0250852858638723e-06,
6.6894371345811248e-06,
6.6894371345811248e-06,
6.7119881029194506e-06,
6.7119881029194506e-06,
5.7209256591358748e-06,
5.7209256591358748e-06,
2.4892855745019102e-06,
2.4892855745019022e-06,
2.4892855745019107e-06,
2.4892855745019026e-06,
6.9754905679119374e-06,
6.9754905679119374e-06,
2.5302436146944877e-05,
2.5302436146944877e-05,
3.8942395990105829e-06,
3.8942395990105956e-06,
8.5599332354909339e-06,
8.5599332354909339e-06,
4.4653967897492783e-06,
4.4653967897492783e-06,
1.6598336446707787e-05,
1.6598336446707787e-05,
1.2258201320069363e-05,
1.2258201320069363e-05,
8.7651263548684881e-06,
8.7651263548684881e-06,
3.6852140899516712e-06,
3.6852140899516712e-06,
7.3000852456940727e-06,
7.3000852456940727e-06,
9.8315737095137571e-06,
9.8315737095137571e-06,
1.2020887935184066e-05,
1.2020887935184066e-05,
2.8154095027525707e-06,
2.8154095027525707e-06,
2.8154095027525707e-06,
2.8154095027525707e-06,
2.7485326600004095e-06,
2.7485326600004095e-06,
5.1090388995422056e-06,
5.1090388995422056e-06,
3.0377746808520834e-06,
3.0377746808520936e-06,
2.1044313312859879e-05,
2.1044313312859879e-05,
3.6888481250880026e-06,
3.6888481250879908e-06,
5.864322301642213e-06,
5.864322301642213e-06,
5.864322301642213e-06,
5.864322301642213e-06,
1.1946884097260775e-05,
1.1946884097260775e-05,
9.6075484670890398e-06,
9.6075484670890398e-06,
2.1158364013845745e-06,
2.1158364013845817e-06,
1.2157855905691063e-05,
1.2157855905691063e-05,
1.9471146941894962e-06,
1.9471146941895026e-06,
2.3674609937496361e-06,
2.3674609937496361e-06,
4.5528347202921819e-06,
4.5528347202921972e-06,
1.5980206430783044e-05,
1.5980206430783044e-05,
1.2426360215850715e-05,
1.2426360215850715e-05,
2.5676103494579796e-06,
2.5676103494579796e-06,
4.5049515732173432e-06,
4.5049515732173432e-06,
7.4667243059855885e-06,
7.4667243059855885e-06,
8.0598674285556254e-06,
8.0598674285556254e-06,
1.5900485201261083e-05,
1.5900485201261137e-05,
6.5184646921073976e-06,
6.5184646921073976e-06,
5.6707570733873469e-06,
5.6707570733873469e-06,
4.3775744894671009e-06,
4.3775744894671009e-06,
2.6823846479925018e-06,
2.6823846479925018e-06,
1.0345585913329086e-05,
1.0345585913329086e-05,
3.8045096432799733e-06,
3.8045096432799733e-06,
2.9799208600695708e-05,
2.9799208600695708e-05,
2.1466757698641991e-06,
2.1466757698642063e-06,
2.1466757698642063e-06,
2.1466757698642063e-06,
4.6241383325928062e-06,
4.624138332592791e-06,
6.8984594971020387e-06,
6.8984594971020387e-06,
7.700329639439457e-06,
7.7003296394394316e-06,
1.8585215458396523e-06,
1.8585215458396585e-06,
2.5497109076416124e-06,
2.5497109076416209e-06,
2.3096920836577815e-06,
2.3096920836577739e-06,
1.5833261533448465e-05,
1.5833261533448414e-05,
2.4246332485205648e-06,
2.4246332485205648e-06,
7.4906367060767402e-06,
7.4906367060767402e-06,
2.536754320695234e-06,
2.536754320695234e-06,
2.536754320695234e-06,
2.536754320695234e-06,
2.1241784696035182e-06,
2.1241784696035182e-06,
7.3958120172860696e-06,
7.3958120172860696e-06,
3.0851667540266744e-06,
3.0851667540266744e-06,
2.7868036021547507e-06,
2.7868036021547507e-06,
2.3779342566921821e-06,
2.3779342566921821e-06,
2.5637965256368412e-06,
2.5637965256368412e-06,
2.2826683326323479e-06,
2.2826683326323555e-06,
7.0461396717574539e-06,
7.0461396717574539e-06,
6.4079945660734448e-06,
6.407994566073466e-06,
1.5639480215120487e-05,
1.5639480215120436e-05,
9.2977112546480898e-06,
9.2977112546480593e-06,
4.4522358619702566e-06,
4.4522358619702566e-06,
2.2759295085051992e-06,
2.2759295085051992e-06,
3.017050257177756e-06,
3.017050257177756e-06,
1.8058102616136639e-05,
1.8058102616136639e-05,
2.9538702405644664e-06,
2.9538702405644664e-06,
2.8768558824385744e-06,
2.8768558824385744e-06,
6.8678395607919293e-06,
6.8678395607919521e-06,
4.5993146899210522e-06,
4.5993146899210522e-06,
1.2453186455452114e-05,
1.2453186455452114e-05,
5.622102141719828e-06,
5.622102141719828e-06,
6.5065881243658253e-06,
6.5065881243658253e-06,
3.0307215367530768e-05,
3.030721536753067e-05,
8.4544121081975283e-06,
8.4544121081975571e-06,
2.6871242326102952e-06,
2.6871242326102952e-06,
1.0049903032756367e-05,
1.0049903032756367e-05,
3.4578609058343109e-06,
3.4578609058343223e-06,
4.9596829187364606e-06,
4.9596829187364606e-06,
7.5365133563458287e-06,
7.5365133563458287e-06,
9.1366016798976709e-06,
9.1366016798976709e-06,
4.3508491464610029e-06,
4.3508491464610029e-06,
5.4876688783356249e-06,
5.4876688783356249e-06,
2.9047079522329073e-06,
2.9047079522329073e-06,
9.5558843161517347e-06,
9.5558843161517347e-06,
1.5507357008991288e-05,
1.5507357008991237e-05,
3.6458271996395158e-06,
3.6458271996395158e-06,
4.8784645187157337e-06,
4.8784645187157337e-06,
1.0844864362655698e-05,
1.0844864362655698e-05,
1.0844864362655698e-05,
1.0844864362655698e-05,
5.3913228008568611e-06,
5.3913228008568611e-06,
4.7490126650263642e-06,
4.7490126650263803e-06,
4.749012665026365e-06,
4.7490126650263811e-06,
3.7838145671962729e-06,
3.7838145671962729e-06,
2.6034740778919736e-06,
2.6034740778919651e-06,
1.3182880671821328e-05,
1.3182880671821328e-05,
5.5225946716667717e-06,
5.5225946716667717e-06,
1.0105407807381765e-05,
1.0105407807381765e-05,
6.2179031629582023e-06,
6.2179031629582023e-06,
9.5704384439842113e-06,
9.5704384439842113e-06,
9.5704384439842113e-06,
9.5704384439842113e-06,
2.5920111783957609e-06,
2.5920111783957609e-06,
4.6860820939350253e-06,
4.6860820939350101e-06,
4.7335837297230066e-06,
4.7335837297230066e-06,
4.7335837297230066e-06,
4.7335837297230066e-06,
8.3415482711847023e-06,
8.3415482711847294e-06,
2.2021274553523634e-06,
2.2021274553523706e-06,
1.1748486480027262e-05,
1.1748486480027225e-05,
4.498113145177278e-06,
4.4981131451772932e-06,
7.093390455774828e-06,
7.093390455774828e-06,
2.6342598817234034e-06,
2.6342598817234034e-06,
7.4440288622649579e-06,
7.4440288622649579e-06,
5.9974206126465701e-06,
5.9974206126465701e-06,
5.9974206126465701e-06,
5.9974206126465701e-06,
2.2322195192882222e-06,
2.2322195192882222e-06,
2.8432885369060113e-06,
2.8432885369060206e-06,
4.3069633970501432e-06,
4.3069633970501296e-06,
5.6841549633625055e-06,
5.6841549633625055e-06,
7.6877512211032536e-06,
7.6877512211032536e-06,
5.0994147950563788e-06,
5.0994147950563788e-06,
5.6700027332742849e-06,
5.6700027332742849e-06,
1.33535638190829e-05,
1.33535638190829e-05,
4.2842563958222573e-06,
4.2842563958222438e-06,
5.3455063494841826e-06,
5.3455063494842004e-06,
2.4801996554816396e-06,
2.4801996554816396e-06,
7.8967850400284471e-06,
7.8967850400284216e-06,
5.4002331032772981e-06,
5.4002331032772981e-06,
1.0087011711616552e-05,
1.0087011711616586e-05,
2.7480426209865213e-06,
2.7480426209865213e-06,
2.7480426209865213e-06,
2.7480426209865213e-06,
9.5906211762135208e-06,
9.5906211762135208e-06,
2.5812248933947347e-06,
2.5812248933947431e-06,
1.308696103595273e-05,
1.308696103595273e-05,
2.8488966394244263e-06,
2.8488966394244263e-06,
2.8488966394244263e-06,
2.8488966394244263e-06,
6.5845310217880837e-06,
6.5845310217880837e-06,
2.469407704938089e-06,
2.469407704938089e-06,
2.469407704938089e-06,
2.469407704938089e-06,
4.0839489167156115e-06,
4.083948916715598e-06,
4.0839489167156115e-06,
4.083948916715598e-06,
5.807070244217728e-06,
5.807070244217728e-06,
3.1735768886781905e-06,
3.1735768886781905e-06,
3.1735768886781905e-06,
3.1735768886781905e-06,
1.1083179304577234e-05,
1.1083179304577234e-05,
4.2332906016316662e-06,
4.2332906016316527e-06,
4.2916770489382128e-05,
4.2916770489382271e-05,
3.5979828734211148e-06,
3.5979828734211267e-06,
1.1060898186244606e-05,
1.1060898186244606e-05,
9.6449896632619698e-06,
9.6449896632619393e-06,
2.6200011555559538e-06,
2.6200011555559454e-06,
6.6429922894492435e-06,
6.6429922894492435e-06,
7.6005622336497618e-06,
7.6005622336497618e-06,
2.8703814943532416e-06,
2.8703814943532416e-06,
2.470743937932078e-06,
2.470743937932078e-06,
1.692592860641768e-05,
1.6925928606417626e-05,
1.2411010758459374e-05,
1.2411010758459414e-05,
4.7846138517818049e-06,
4.784613851781821e-06,
3.3943403657698852e-06,
3.3943403657698852e-06,
9.6405934130404105e-06,
9.64059341304038e-06,
3.7089589517299131e-06,
3.7089589517299131e-06,
4.5764784666867309e-06,
4.5764784666867461e-06,
4.146440477196742e-06,
4.1464404771967556e-06,
3.2081321448154949e-06,
3.2081321448155055e-06,
4.0897582324914742e-06,
4.0897582324914607e-06,
1.4743593520155599e-05,
1.4743593520155599e-05,
7.6538738305332816e-06,
7.6538738305332562e-06,
1.0326292648866057e-05,
1.0326292648866057e-05,
5.037593953469442e-06,
5.037593953469442e-06,
6.5252955692634284e-06,
6.5252955692634284e-06,
8.8298615289439642e-06,
8.8298615289439354e-06,
3.9493307481852874e-06,
3.9493307481853009e-06,
2.2686492704910944e-06,
2.2686492704910944e-06,
5.8627620584156257e-06,
5.8627620584156257e-06,
7.0598552055305331e-06,
7.0598552055305331e-06,
2.7830466705828831e-06,
2.7830466705828743e-06,
1.8682596451835844e-06,
1.8682596451835785e-06,
5.2642185634922693e-06,
5.2642185634922693e-06,
2.7015945047483668e-06,
2.7015945047483579e-06,
4.8976988879809357e-06,
4.8976988879809518e-06,
7.195276346841727e-06,
7.195276346841727e-06,
2.6976574830346886e-06,
2.6976574830346886e-06,
2.348769907445799e-06,
2.3487699074458066e-06,
4.3719615765360043e-06,
4.3719615765360043e-06,
4.3719615765360043e-06,
4.3719615765360043e-06,
2.7715276068037098e-06,
2.7715276068037098e-06,
1.8174035695663269e-05,
1.8174035695663211e-05,
1.8174035695663272e-05,
1.8174035695663214e-05,
7.4227653848903402e-06,
7.4227653848903164e-06,
2.630202726683042e-06,
2.6302027266830336e-06,
3.3075373119923601e-06,
3.3075373119923601e-06,
4.2406130175371966e-06,
4.240613017537211e-06,
4.240613017537211e-06,
4.240613017537211e-06,
5.7782581134882681e-06,
5.7782581134882495e-06,
1.2675911323158365e-05,
1.2675911323158365e-05,
1.559266276804977e-06,
1.5592662768049823e-06,
4.405966505861971e-06,
4.4059665058619854e-06,
4.9469667248780752e-06,
4.9469667248780591e-06,
8.1134401379502934e-06,
8.1134401379502934e-06,
3.2864324786096616e-06,
3.2864324786096616e-06,
5.4905796770422926e-06,
5.4905796770422749e-06,
3.3842474984026035e-06,
3.3842474984026149e-06,
6.2075611651340541e-06,
6.2075611651340541e-06,
6.2075611651340541e-06,
6.2075611651340541e-06,
2.5078353063154506e-06,
2.5078353063154426e-06,
7.9368690254363151e-06,
7.9368690254362897e-06,
7.9368690254362897e-06,
7.9368690254362897e-06,
3.3070065395908575e-06,
3.3070065395908685e-06,
...]
In [99]:
dc2_df['x0'] = x0
In [101]:
dc2_df['x1'] = 0.
dc2_df['c'] = 0.
In [103]:
dc2_df.to_csv('/Users/Bryce/Documents/DESC/Twinkles/data/dc2_sne_cat.csv', index=False)
In [75]:
imSimBand = Bandpass()
print(imSimBand)
imSimBand.imsimBandpass()
<lsst.sims.photUtils.Bandpass.Bandpass instance at 0x11588dfc8>
In [79]:
imSimBand.wavelen[1990:2010]
Out[79]:
array([ 499. , 499.1, 499.2, 499.3, 499.4, 499.5, 499.6, 499.7,
499.8, 499.9, 500. , 500.1, 500.2, 500.3, 500.4, 500.5,
500.6, 500.7, 500.8, 500.9])
In [81]:
imSimBand.sb[1990:2010]
Out[81]:
array([ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 1., 0., 0.,
0., 0., 0., 0., 0., 0., 0.])
In [84]:
dc2_df.iloc[1]
Out[84]:
t0 14.945
x -0.125
y 0.1134
sigma 171.981
gamma 0.0738287
e 0.31621
theta_e 79.1602
theta_gamma 134.088
zs 0.734945
zl 0.399013
snx 0.102703
sny -0.149628
MB -19.535
mu -0.9887
sysno 105
imno 1
mag 23.6654
x0 5.00048e-06
x1 1
c 0
twinkles_sysno 1000
bulge_magnorm 20.8513
r_eff 1.05134
lens_sed Burst.62E09.04Z.spec
t_start 60367.3
Name: 1, dtype: object
In [ ]:
Content source: DarkEnergyScienceCollaboration/Twinkles
Similar notebooks: