LSSTDESC / SprintWeek2022

Meeting repository for the LSST DESC 2022 Sprint Week
Creative Commons Zero v1.0 Universal
1 stars 0 forks source link

Interface to DSPS in RAIL #12

Closed torluca closed 2 years ago

torluca commented 2 years ago

Interface to DSPS in RAIL

The aim of this sprint is to start the implementation of an interface in RAIL to the stellar population synthesis code "DSPS" (Hearin+21).

Contacts: Luca Tortorelli, Andrew Hearin Day/Time: 10/17, 18/10, 19/10, from 10 am CET to 9 pm CET Main communication channel: Slack GitHub repo: https://github.com/LSSTDESC/RAIL/tree/issue/268/dspsinterface Zoom room (if applicable): https://lmu-munich.zoom.us/j/95536091814?pwd=VlBJUmY2Tmg4VFd6T3FMcTg1dTFZZz09

Goals and deliverable

The goal for the week is to accomplish the milestone: "Design and implementation of the interface in rail.creation that allows the generation of intrinsic SEDs with DSPS." detailed in the GER-LMU workplan, together with the deliverable: "Documented contributions in the GitHub repository of RAIL."

Resources and skills needed

Knowledge of Python, RAIL and how stellar population synthesis codes work is needed.

Detailed description

We will use the base class SedGenerator (see rail_fsps) that takes as input data in the form of a file passed by tables_io and provide as output another tables_io readable file. SedGenerator is in turns a derived class of RailStage. The class DSPSSedGenerator will therefore be a derived class of SedGenerator that needs an input table of galaxy properties, as well as configuration options for which stellar population component the user might want to add.

torluca commented 2 years ago

Zoom Room #15 https://stanford.zoom.us/j/94385118949?pwd=TWhGT3Nsd3ZBU0l5RFJ1NUJsNENjdz09 on Monday morning

torluca commented 2 years ago

First day report: Following a discussion with Andrew Hearin, it has been implemented in RAIL.creation.galaxy_modelling two derived classes, DSPSSingleSedModeler and DSPSPopulationSedModeler, that interfaces with DSPS to genereate a single or populations of rest frame SEDs based on provided grid of ages, metallicities and templates SSPs.

torluca commented 2 years ago

Final report: We completed the implementation of the interface to DSPS by creating a new class, DSPSPhotometryCreator, in rail.creation.galaxy_modelling.dsps_photometry_creator.py, that generates synthetic photometric data from provided rest-frame seds. The photometry depends on the provided SEDs, filters, galaxy redshift and cosmological parameters. The code can be found in the folder src/rail/creation/galaxy_modelling. Only unit tests and demo notebook are missing.