Identifying Optical GW Counterparts with Image Processing Noemi Glaeser, University of South Carolina SIST Final Talk 6 August 2018
Overview • Background – GW – DES • The DES-GW pipeline – SE – DiffImg • Next steps 2 06/08/18 N. Glaeser | SIST Final Talk
Gravitational Waves LIGO/T. Pyle 2016 “Perpetuating changes in the gravitational field caused by accelerating masses” 3 06/08/18 N. Glaeser | SIST Final Talk
Gravitational Waves LIGO/T. Pyle 2016 Binary Black Hole (BBH) Black Hole-Neutron Star (BH-NS) Binary Neutron Star (BNS) 4 06/08/18 N. Glaeser | SIST Final Talk
Gravitational Waves LIGO/Caltech/MIT/LSC 2017 Wave shape → event type Amplitude → distance 5 06/08/18 N. Glaeser | SIST Final Talk
Gravitational Waves LIGO/Virgo/NASA/Leo Singer 2017 Dark Energy Survey Collaboration (Milky Way image: Axel Mellinger) 6 06/08/18 N. Glaeser | SIST Final Talk
Dark Energy Survey (DES) Reidar Hahn (Fermilab) DES Collaboration DES Collaboration DECam, assembled at Fermilab, now mounted on the Blanco telescope in Chile 7 06/08/18 N. Glaeser | SIST Final Talk
Dark Energy Survey (DES) DES Collaboration The DES footprint spans ~5000 square degrees (⅛ of the sky) 8 06/08/18 N. Glaeser | SIST Final Talk
Dark Energy Survey (DES): Why? Independent measure of the Hubble parameter and other cosmological parameters DES Collaboration 9 06/08/18 N. Glaeser | SIST Final Talk
Dark Energy Survey (DES): GW170817 GW170817 Optical Counterpart composite detection images (Soares-Santos et al. 2017) • The first optical counterpart of a GW event was observed by DES in August 2017! – One point on the distance-redshift plot • In O3 we expect ~10x as many events 10 06/08/18 N. Glaeser | SIST Final Talk
Dark Energy Survey (DES): Why? Independent measure of the Hubble parameter and other cosmological parameters DES Collaboration 11 06/08/18 N. Glaeser | SIST Final Talk
DES-GW Image Processing Pipeline 12 06/08/18 N. Glaeser | SIST Final Talk
DES-GW Pipeline Full Pipeline 13 06/08/18 N. Glaeser | SIST Final Talk
DES-GW Pipeline • Single Epoch (SE) processing • Intermediate steps • DiffImg Full Pipeline 14 06/08/18 N. Glaeser | SIST Final Talk
DES-GW Pipeline • Single Epoch (SE) processing • Intermediate steps • DiffImg Full Pipeline 15 06/08/18 N. Glaeser | SIST Final Talk
DES-GW Pipeline • Single Epoch (SE) processing • Intermediate steps • DiffImg Full Pipeline 16 06/08/18 N. Glaeser | SIST Final Talk
DES-GW Pipeline • Single Epoch (SE) processing • Intermediate steps • DiffImg Runtime: ~5-8 hr Full Pipeline 17 06/08/18 N. Glaeser | SIST Final Talk
DES-GW Pipeline • Single Epoch (SE) processing • Intermediate steps • DiffImg Runtime: ~5-8 hr Goal: <1 hr Full Pipeline 18 06/08/18 N. Glaeser | SIST Final Talk
DES-GW Pipeline • Single Epoch (SE) processing • Intermediate steps • DiffImg Runtime: ~5-8 hr Goal: <1 hr Full Pipeline 19 06/08/18 N. Glaeser | SIST Final Talk
Image Processing Pipeline: SE processing • Image Correction – Raw Images → “science-ready” • Astrometric calibration: 2MASS • Object Cataloguing • Photometric Calibration: 2MASS Single Epoch Processing 20 06/08/18 N. Glaeser | SIST Final Talk
Image Processing Pipeline: SE processing • Image Correction – Raw Images → “science-ready” • Astrometric calibration: 2MASS • Object Cataloguing • Photometric Calibration: 2MASS Single Epoch Processing 21 06/08/18 N. Glaeser | SIST Final Talk
Image Processing Pipeline: SE processing • Image Correction – Raw Images → “science-ready” • Astrometric calibration: 2MASS • Object Cataloguing • Photometric Calibration: 2MASS Single Epoch Processing 22 06/08/18 N. Glaeser | SIST Final Talk
Image Processing Pipeline: SE processing • Image Correction – Raw Images → “science-ready” • Astrometric calibration: 2MASS • Object Cataloguing • Photometric Calibration: 2MASS Single Epoch Processing 23 06/08/18 N. Glaeser | SIST Final Talk
Image Processing Pipeline: SE processing • Image Correction – Raw Images → “science-ready” • Astrometric calibration: 2MASS • Object Cataloguing • Photometric Calibration: 2MASS Single Epoch Processing 24 06/08/18 N. Glaeser | SIST Final Talk
Image Processing Pipeline: SE processing • Image Correction – Raw Images → “science-ready” • Astrometric calibration: 2MASS • Object Cataloguing • Photometric Calibration: 2MASS Single Epoch Processing 25 06/08/18 N. Glaeser | SIST Final Talk
Image Processing Pipeline: SE processing • Image Correction – Raw Images → “science-ready” • Astrometric calibration: 2MASS – CCD-by-CCD • Object Cataloguing • Photometric Calibration: 2MASS – CCD-by-CCD Single Epoch Processing 26 06/08/18 N. Glaeser | SIST Final Talk
Image Processing Pipeline: SE processing • Image Correction – Raw Images → “science-ready” • Astrometric calibration: 2MASS – CCD-by-CCD • Object Cataloguing • Photometric Calibration: 2MASS – CCD-by-CCD Single Epoch Processing 27 06/08/18 N. Glaeser | SIST Final Talk
Image Processing Pipeline: SE processing • Image Correction – Raw Images → “science-ready” • Astrometric calibration: GAIA-DR2 – CCD-by-CCD • Object Cataloguing • Photometric Calibration: GAIA-DR2 – CCD-by-CCD Single Epoch Processing 28 06/08/18 N. Glaeser | SIST Final Talk
Image Processing Pipeline: SE processing • Image Correction – Raw Images → “science-ready” • Astrometric calibration: GAIA-DR2 – CCD-by-CCD • Object Cataloguing • Photometric Calibration: GAIA-DR2 – CCD-by-CCD Single Epoch Processing 29 06/08/18 N. Glaeser | SIST Final Talk
Image Processing Pipeline: SE Modifications Raw image • “Fast” CCDs don’t have to wait for slower ones CCD62 ... CCD2 CCD1 Photometric calibration 30 06/08/18 N. Glaeser | SIST Final Talk
Image Processing Pipeline: SE Modifications Raw image • “Fast” CCDs don’t have to wait for slower ones CCD62 ... CCD2 Photometric calibration CCD1 31 06/08/18 N. Glaeser | SIST Final Talk
Image Processing Pipeline: SE Modifications Raw image • “Fast” CCDs don’t have to wait for slower ones CCD62 ... Photometric calibration CCD2 CCD1 32 06/08/18 N. Glaeser | SIST Final Talk
Image Processing Pipeline: SE Modifications • “Fast” CCDs don’t have to wait for slower ones Photometric calibration CCD62 ... CCD2 CCD1 33 06/08/18 N. Glaeser | SIST Final Talk
Image Processing Pipeline: SE Modifications • “Fast” CCDs don’t have to wait for slower ones Photometric calibration Science-ready image 34 06/08/18 N. Glaeser | SIST Final Talk
Image Processing Pipeline: SE Modifications • “Fast” CCDs don’t have to wait for slower ones Photometric calibration Science-ready to DiffImg image 35 06/08/18 N. Glaeser | SIST Final Talk
Image Processing Pipeline: SE Modifications Raw image • “Fast” CCDs don’t have to wait for slower ones CCD1 CCD2 ... CCD62 Photometric calibration 36 06/08/18 N. Glaeser | SIST Final Talk
Image Processing Pipeline: SE Modifications Raw image • “Fast” CCDs don’t have to wait for slower ones CCD1 ... CCD62 Photometric calibration CCD2 37 06/08/18 N. Glaeser | SIST Final Talk
Image Processing Pipeline: SE Modifications Raw image • “Fast” CCDs don’t have to wait for slower ones CCD1 ... CCD62 Photometric calibration CCD2 to DiffImg 38 06/08/18 N. Glaeser | SIST Final Talk
Image Processing Pipeline: SE Modifications • “Fast” CCDs don’t have to wait for slower ones Single Epoch Parallelization Speedup ( series : 0.9625 h; parallel : µ=0.15, median=0.13, σ=0.05) 39 06/08/18 N. Glaeser | SIST Final Talk
Image Processing Pipeline: SE Modifications • “Fast” CCDs don’t have to wait for slower ones • But… Single Epoch Parallelization Speedup ( series : 0.9625 h; parallel : µ=0.15, median=0.13, σ=0.05) 40 06/08/18 N. Glaeser | SIST Final Talk
Image Processing Pipeline: SE Modifications • “Fast” CCDs don’t have to wait for slower ones • But… – do the results still make sense? – did we break something down the line? – Test on GW170817 Single Epoch Parallelization Speedup ( series : 0.9625 h; parallel : µ=0.15, median=0.13, σ=0.05) 41 06/08/18 N. Glaeser | SIST Final Talk
Image Processing Pipeline: SE Modifications • Run newly SE-processed images through DiffImg Make sure we still – identify the counterpart 42 06/08/18 N. Glaeser | SIST Final Talk
Image Processing Pipeline: SE Modifications • Run newly SE-processed images through DiffImg Make sure we still – identify the counterpart • 135 unmatched objects (45%)! Non-matching candidates with old SE (all 2MASS) vs. new SE (all GAIA) 43 06/08/18 N. Glaeser | SIST Final Talk
Recommend
More recommend
