DES SN 5YR Tutorial Part 3: Load LCPlot file¶
Dark Energy Survey Supernova Program
Tutorial Part 3 Load LC Plot files Contact author: Bruno Sánchez
Import our utils module to load the LCPLOT files in the Release
[1]:
import warnings
warnings.filterwarnings("ignore")
[2]:
import os
[3]:
from dessndr import utils, data
[4]:
lcplot_data = utils.read_lcplot(
os.path.join(
data.DES5YRDR_DATA_ROOT,
'DES-SN5YR/0_DATA/DES5YR_SALT3_LCFIT.LCPLOT.gz'
)
)
Inspect around. These files contain the lightcurves twice. Once for the DATA and once for the MODEL fitted fluxes.
[5]:
lcplot_data.iloc[0]
[5]:
SNID 1562098
MJD 57701.07
PHASE -19.64
FLUXCAL 12.53
FLUXCALERR 8.9893
DATA_MODEL 1
BAND r
CHI2 1.3
IGNOREME 1
Name: 0, dtype: object
[6]:
len(lcplot_data.SNID.unique())
[6]:
1635
[7]:
lcplot_data.SNID.unique()[:10]
[7]:
array([1562098, 1702083, 1435094, 1309227, 1324542, 1320474, 1290816,
1298246, 1311632, 1302738])
[8]:
lcplot_data.head()
[8]:
| SNID | MJD | PHASE | FLUXCAL | FLUXCALERR | DATA_MODEL | BAND | CHI2 | IGNOREME | |
|---|---|---|---|---|---|---|---|---|---|
| 0 | 1562098 | 57701.070 | -19.64 | 12.530 | 8.9893 | 1 | r | 1.30 | 1 |
| 1 | 1562098 | 57708.094 | -12.62 | 102.850 | 41.4340 | 1 | r | 1.32 | 1 |
| 2 | 1562098 | 57709.035 | -11.68 | 67.120 | 5.4472 | 1 | r | 1.21 | 1 |
| 3 | 1562098 | 57714.188 | -6.53 | 74.187 | 4.4849 | 1 | r | 1.33 | 1 |
| 4 | 1562098 | 57718.070 | -2.64 | 80.931 | 5.5057 | 1 | r | 0.26 | 1 |
Use the loaded data to look at light-curves with the SALT3 fitted models!
[9]:
from dessndr import plot
[10]:
ax = plot.plot_lcplot(lcplot_data, cid=1562098)
[11]:
ax = plot.plot_lcplot(lcplot_data, cid=1702083)
[ ]: