Skies and Universes Johan Comparat (UAM CSIC IFT) SelGIFS workshop, - - PowerPoint PPT Presentation

skies and universes
SMART_READER_LITE
LIVE PREVIEW

Skies and Universes Johan Comparat (UAM CSIC IFT) SelGIFS workshop, - - PowerPoint PPT Presentation

Skies and Universes Johan Comparat (UAM CSIC IFT) SelGIFS workshop, Universidad autonoma de madrid november 2016 Skies and universes, rationale To overcome the cultural barrier and communicate or co-work efficiently, one needs a solid structure


slide-1
SLIDE 1

Skies and Universes

Johan Comparat (UAM CSIC IFT) SelGIFS workshop, Universidad autonoma de madrid november 2016

slide-2
SLIDE 2

Skies and universes, rationale

To overcome the cultural barrier and communicate or co-work efficiently, one needs a solid structure that in the case of astrophysics or cosmology can be reduced to a cluster of computers where data, code and resources are in open access and documented. It aims to be a very special place to connect ideas to data and theory in the most versatile, fluid manner. Currently exist many high-quality databases related to telescopes and virtual observatories: www.sdss.org, http://vizier.u-strasbg.fr/, etc. But these are not directly connected with scientific code repositories, theoretical predictions nor computing facilities. This multiplication of interfaces between the scientist and these three entities, leads to important time losses. Such a co-working space also enables better reproducibility and easier teamwork.

11/21/2016

  • J. Comparat

2

slide-3
SLIDE 3

11/21/2016

  • J. Comparat

3

code data simulations http://projects.ift.uam-csic.es/skies-universes/

Collaborators:

  • J. Comparat
  • F. Prada
  • G. Favole
  • S. Rodriguez-Torres
  • P. Zarrouk
  • E. Burtin
  • C. Maraston
  • V. Gonzalez-Perez
slide-4
SLIDE 4

Software & Database

Classes Libraries

DATA RELATED Class GalaxySurvey{SurveyName}

  • Overall calibration methods
  • Converts flux => luminosities

Class GalaxySpectrum{SurveyName}

  • Opens spectra
  • Spectra calibration
  • plotting

SCIENCE RELATED Class LineLuminosityFunction

  • derives LF from a catalog: completeness, weighting, uncertainty

Class ModelLuminosityFunction

  • fitModel: Schechter or Saunders

Libraries

  • filterList
  • Decam, megacam (cfht), sdss broad band filter curves
  • LineFittingLibrary
  • Integrates spectrum, fits emission line (single or doublet) models (Gaussian, Lorentzian, pseudoVoigt)
  • lib_plot
  • Produce standard plots from the catalogs
  • lineListAir / Vac
  • Line list to consider. Depends on pyneb

11/21/2016

  • J. Comparat

4

slide-5
SLIDE 5

Emission-line-luminosity-selected galaxy samples

Based on Comparat et al. 2016. MNRAS. Arxiv 1605.02875

11/21/2016

  • J. Comparat

5

slide-6
SLIDE 6

DATA

  • Use major spectroscopic magnitude limited surveys that cover more

than a square degree

  • Fold data into a common data base (skies and universes) to enhance

control over flux limits determination, weighting and completeness

11/21/2016

  • J. Comparat

6

slide-7
SLIDE 7

DATA

11/21/2016

  • J. Comparat

7

slide-8
SLIDE 8

Workflow DATA SCIENCE

  • flux calibration for DEEP2 spectra
  • fit emission line to spectra (VVDS, DEEP2, VIPERS)
  • create emission line catalogs
  • convert observed flux into luminosity
  • Compute TSR and SSR
  • Determine magnitude limit and flux limit

11/21/2016

  • J. Comparat

8

  • Estimate line luminosity function
  • Convert references to Planck cosmology
  • Fit Saunders and Schechter models
slide-9
SLIDE 9

Fitting Example

Signal to noise is low therefore convergence is tricky. It depends on * resolution * redshift * line flux * estimation of the continuum * estimation of the uncertainty

  • n the flux

* …

11/21/2016

  • J. Comparat

9

slide-10
SLIDE 10

Results, [OII] LF

11/21/2016

  • J. Comparat

10

slide-11
SLIDE 11

[OIII] 5007 LF

11/21/2016

  • J. Comparat

11

slide-12
SLIDE 12

H beta LF: new !

11/21/2016

  • J. Comparat

12

slide-13
SLIDE 13

Fit results

11/21/2016

  • J. Comparat

13

slide-14
SLIDE 14

Results

  • Skies and universes started a year

ago

  • It hosts a few projects:
  • 2 published (1507.04356,

1605.02875), the latter is 100% end- to-end reproducible.

  • 2 in the referee process (1611.05457)
  • 4 in preparation
  • And many more to come :-)
  • New collaborators are welcome !
  • Emission line LFs
  • The models plugged-in to an analysis

allow comparison to any observed field: “where do I stand with respect to the mean density of line emitters.”

  • Great tool to plan observations,

understand and forward model completeness

  • In preparation
  • relation between line LF and stellar

populations models

  • Conditional emission line LFs

11/21/2016

  • J. Comparat

14

slide-15
SLIDE 15

Exercise

11/22/2016

  • J. Comparat

15

Arxiv: 1507.07979

slide-16
SLIDE 16

Fitting emission line models

  • In the folder “ELG-

composite” taken from Zhu, Comparat et al. arxiv:1507.07979. It contains :

  • High SNR ELG stack:

average properties of

  • bserved ELGs
  • Two scripts
  • Three data files
  • Open and read

“pysu_tutorial.py”

  • 1. With the composite spectrum:
  • 1.a. Fit emission line models to the the [OIII] 5007 line existing

in the composite spectrum provided here :

  • 1.b. plot the best emission line models found along with the

data, discuss the chi2 / degrees of freedom found.

  • 1.c. Artificially increase / decrease the uncertainty on the flux

(errors) and look at how the line model parameters vary. Discuss the importance of the error.

  • 1.d. Fit models to other lines: [OII] doublet, H beta. Discuss

the challenges corresponding to each lines.

  • 2. With the stellar population model (SPM):
  • 2.a. compare the SPM spectrum provided and the composite
  • spectra. Is it a reasonable model ?
  • 2.b subtract the SPM from the composite and fit emission line

models to the residual spectrum. Be careful with the estimation of the errors on the residual spectrum

11/21/2016

  • J. Comparat

16

slide-17
SLIDE 17

11/22/2016

  • J. Comparat

17