Summary of MD04.pv2_stacktest.20130903

(back to MediumDeepFields summary)

This is a test sample for PV2 update processing, refstack preparation with extra test versions and deep stack tests using MD04 while it is an inactive field.

PV2: update type reprocessing

Deep stacks:

Reference stack -- this should be useful for making the MD04 refstack for return to field early November 2013. May move section for refstack to its own page for consistency.

Other stacks:

  • in and out of season stacks -- haven't been done yet, need sample for Taiwan meeting
  • fewer N input refstack

Status:

  • prior 9/3 many issues being worked out for PV2 and MD (MD10.refstack.20130715)
  • 9/3 PV2 updates started for MD04 after M31 reprocessing finished
  • 9/6 PV2 updates finished
  • 10/3 test all-in deep stack skycell run
  • 10/10 started full all-in deep stack running
  • 11/03 deep stack test with cuts started -- deep stack tests having problems
  • 11/22 (great amount of updating to do)
  • 11/26 all stack and staticsky should be finished and available on datastore now -- all stack types should have label and data_group to indicate what they are

PV2 updates

MD04

PV1 camera label -- data_group: mix of processing labels/data_groups, but all of PV1 type (cam_id>=434615)

label data_group
MD04.pv1.20130322 MD04.pv1.20130322
MD04.refstack.20121125 MD04.refstack.20121125
MD04.nightlyscience MD04.<obsdate 20121129--20130525>

There have been some test reprocessing that had to be avoided -- czw.MDVDtest%, czw.MDVD.MD04.%

PV2 camera (as well as chip, warp) label MD04.pv2.20130903, data_group left to what originally set as in order to be able to track backwards to the processing group.

With ipp047 data disk down, have found ~3% of exposures missing a related imfile (~half are the burn.tbl files that can be regenerated) and indicates serious need of fixing the replication of raw data. Also found some replicated but on a host (ipp064X) that no longer exists.. it is also possible some replicated products may be on the same host if push to the stsci nodes with 3 partitions..

Number totals through Sept. 2013: (will be more by end of survey)

  • not all exposures pass camera stage, FWHM>12 pixels or other quality issues
    | filter  | count(cam_id) | count(warp_id) |
    +---------+---------------+----------------+
    | g.00000 |           596 |            577 | 
    | r.00000 |           588 |            571 | 
    | i.00000 |           719 |            713 | 
    | z.00000 |           711 |            710 | 
    | y.00000 |           471 |            465 | 
    

PV2 warps will also be kept available (not cleaned up) for making deep stacks but also will be a accessible via the PSS/PSPS interface more quickly (release name MD04.PV2).

3PI overlap

Unclear what want to do here still, with LAP running, and how useful it will be in the stack -- like MD10.refstack.20130715#MD3PIstack, will make a sample refstack with 3PI addition


Deep Stacks

Like with MD09, different versions of stacks will be produced for testing (not all will be full-field?)

  • all --
  • moderate input cuts --
  • best cuts >2010 and moderate input cuts --

All

Basically all warps available with ~FWHM <13 pixels to avoid any badly measured PSF warps.

  • stack -- label, data_group = MD04.alldeeptest.20131010
  • staticsky -- label, data_group = MD04.alldeeptest.20131010.Nx20131122, MD04.alldeeptest.20131010.<#>x20131122 -- where <#> is the number of filters available
  • skycal -- label, data_group = MD04.alldeeptest.20131010.Nx20131122.cal20131122, MD04.alldeeptest.20131010.Nx20131122.cal20131122

Updated stack photometry (staticsky+skycal stage) run for next MD loading into the PSPS/sandbox --

  • staticsky -- label, data_group = MD04.alldeeptest.20131010.Nx20131216, MD04.alldeeptest.20131010.<#>x20131216 -- where <#> is the number of filters available
  • skycal -- label, data_group = MD04.alldeeptest.20131010.Nx20131216.cal20131216, MD04.alldeeptest.20131010.Nx20131216.cal20131216

Cut

All warps FWHM <10 pixels and input warp cuts

  • ZPT > ZPTmean-0.33
    g: 24.11
    r: 24.32
    i: 24.07
    z: 23.81
    y: 22.79
    
  • elongation <1.5
  • background:
    g: <400
    r: <500
    i: <1800
    z: <2500
    y: <5000
    
  • stack -- label, data_group = MD04.cutdeeptest.20131103
  • staticsky -- label, data_group = MD04.cutdeeptest.20131103.Nx20131122, MD04.cutdeeptest.20131103.<#>x20131122 -- where <#> is the number of filters available
  • skycal -- label, data_group = MD04.cutdeeptest.20131103.Nx20131122.cal20131126, MD04.cutdeeptest.20131103.Nx20131122.cal20131126

Best Cut

All warps FWHM <10 pixels and the same input warp cuts, only season years 2-4 (excludes 2009/2010 season)

  • stack -- label, data_group = MD04.cut2009deeptest.20131103
  • staticsky -- label, data_group = MD04.cut2009deeptest.20131103.Nx20131122, MD04.cut2009deeptest.20131103.<#>x20131122 -- where <#> is the number of filters available
  • skycal -- label, data_group = MD04.cut2009deeptest.20131103.Nx20131122.cal20131122, MD04.cut2009deeptest.20131103.Nx20131122.cal20131122

Summary

Sample script runtime (dtime_script) for skycell.045

| stack_id | label                     | filter  | skycell_id  | state | dtime_script | reject_images | count(warp_id) |
+----------+---------------------------+---------+-------------+-------+--------------+---------------+----------------+
|  2830328 | MD04.alldeeptest.20131010 | g.00000 | skycell.045 | full  |        55222 |             0 |            577 | 
|  2854300 | MD04.alldeeptest.20131010 | r.00000 | skycell.045 | full  |        52963 |             0 |            571 | 
|  2854220 | MD04.alldeeptest.20131010 | i.00000 | skycell.045 | full  |        73856 |             0 |            713 | 
|  2830487 | MD04.alldeeptest.20131010 | z.00000 | skycell.045 | full  |        67972 |             0 |            710 | 
|  2830408 | MD04.alldeeptest.20131010 | y.00000 | skycell.045 | full  |        42284 |             0 |            465 | 


| stack_id | label                     | filter  | skycell_id  | state | dtime_script | reject_images | count(warp_id) |
+----------+---------------------------+---------+-------------+-------+--------------+---------------+----------------+
|  2900916 | MD04.cutdeeptest.20131103 | g.00000 | skycell.045 | full  |        52056 |            27 |            534 | 
|  2900995 | MD04.cutdeeptest.20131103 | r.00000 | skycell.045 | full  |        52519 |             9 |            515 | 
|  2901073 | MD04.cutdeeptest.20131103 | i.00000 | skycell.045 | full  |        65598 |            12 |            648 | 
|  2901152 | MD04.cutdeeptest.20131103 | z.00000 | skycell.045 | full  |        59987 |             7 |            601 | 
|  2901231 | MD04.cutdeeptest.20131103 | y.00000 | skycell.045 | full  |        37040 |             5 |            391 | 


| stack_id | label                         | filter  | skycell_id  | state | dtime_script | reject_images | count(warp_id) |
+----------+-------------------------------+---------+-------------+-------+--------------+---------------+----------------+
|  2885595 | MD04.cut2009deeptest.20131103 | g.00000 | skycell.045 | full  |        34472 |            22 |            374 | 
|  2885770 | MD04.cut2009deeptest.20131103 | r.00000 | skycell.045 | full  |        34300 |             6 |            364 | 
|  2885830 | MD04.cut2009deeptest.20131103 | i.00000 | skycell.045 | full  |        36800 |             4 |            376 | 
|  2885890 | MD04.cut2009deeptest.20131103 | z.00000 | skycell.045 | full  |        42155 |             7 |            413 | 
|  2885950 | MD04.cut2009deeptest.20131103 | y.00000 | skycell.045 | full  |        21221 |             5 |            227 | 

  • all deep does have other skycells with rejected inputs (i.e., s078,082 g-band reject 10 inputs as FWHM>14 pixels)
  • in cut and cut2009, the reject_images come from the <10 pixel input cut in the configuration, and s045 is an above average number

Other comparison details:

http://ps1sc.ifa.hawaii.edu/PS1wiki/index.php/Kp12tw_stacking_analysis good comparison of various stacks including these.


Reference Stacks

In addition to the normal reference stack for a season, test samples of planned final reference stack types also being made.

| exp_name    | filter  | dateobs             | comment                         |
+-------------+---------+---------------------+---------------------------------+
| o4777g0224o | g.00000 | 2008-11-07 14:56:59 | MD04 g MJD54777 MD04 center     | 
..
| o4803g0354o | r.00000 | 2008-12-03 15:25:16 | MD04 r MJD54803 MD04 center     | 
---- pre-survey

| o5161g0222o | g.00000 | 2009-11-26 13:55:38 | MD04 g,r MJD55161 g MD04 center | 
..
| o5321g0014o | i.00000 | 2010-05-05 06:20:20 | MD04 i N5321 MD04 center        | 
---- yr1 -- 11/26/2009 -- 5/5/2010

| o5545g0436o | g.00000 | 2010-12-15 12:12:35 | MD04 g,r N5545 g MD04 center    | 
..
| o5681g0089o | i.00000 | 2011-04-30 07:34:50 | MD04 i N5681 MD04 center        | 
---- yr2 -- 12/15/2010 -- 4/30/2011

| o5892g0155o | g.00000 | 2011-11-27 14:41:57 | MD04 g,r N5892 g MD04 center    | 
..
| o6056g0114o | z.00000 | 2012-05-09 07:38:31 | MD04 z N6056 MD04 center        | 
---- yr3 -- 11/27/2011 -- 5/09/2012

| o6260g0448o | y.00000 | 2012-11-29 15:14:41 | MD04 y N6260 MD04 center        | 
..
| o6437g0061o | y.00000 | 2013-05-25 05:50:54 | MD04 y N6437 MD04 center        | 
---- yr4 -- 11/29/2012 -- 5/25 2013

---- yr5 -- actively being observed since 2013-11-08

(more details to follow)

Basic Reference Stack

Default reference stacks for nightly diffims

  • reduction MD_REF_STACK_1DG, FWHM <6 pixel and 1D convolution
  • best ~60 input warps excluding 2009--2010 year/season
  • stack -- label, data_group = MD04.refstack.20131018, MD04.refstack.20131018_1dg -- _1dg to indicate 1D convolution method, but is now default for all stacks
  • staticsky -- label, data_group = MD04.refstack.20131018.Nx20131122, MD04.refstack.20131018.<F>x20131122 -- where <F> is the number of filters available for the skycell
  • skycal -- label, data_group = MD04.refstack.20131018.Nx20131122.cal20131122, MD04.refstack.20131018.Nx20131122.cal20131122

Summary info --

| stack_id | filter  | skycell_id  | data_group                 | dtime_stack | reject_images | count(warp_id) |
+----------+---------+-------------+----------------------------+-------------+---------------+----------------+
|  2870506 | z.00000 | skycell.055 | MD04.refstack.20131018_1dg |     6912.08 |            10 |             75 | 
|  2870566 | i.00000 | skycell.055 | MD04.refstack.20131018_1dg |     6530.24 |            13 |             75 | 
|  2870746 | y.00000 | skycell.055 | MD04.refstack.20131018_1dg |     6004.42 |            10 |             68 | 
|  2870806 | r.00000 | skycell.055 | MD04.refstack.20131018_1dg |     6299.23 |            14 |             76 | 
|  2870866 | g.00000 | skycell.055 | MD04.refstack.20131018_1dg |     5184.99 |            16 |             74 |


| stack_id | filter  | skycell_id  | data_group             | dtime_stack | reject_images | count(warp_id) |
+----------+---------+-------------+------------------------+-------------+---------------+----------------+
|  1878094 | y.00000 | skycell.055 | MD04.refstack.20121125 |     13235.3 |            25 |             92 | 
|  1878155 | z.00000 | skycell.055 | MD04.refstack.20121125 |     14771.3 |            18 |             89 | 
|  1878215 | i.00000 | skycell.055 | MD04.refstack.20121125 |     13255.8 |            20 |             92 | 
|  1878275 | r.00000 | skycell.055 | MD04.refstack.20121125 |     14253.2 |            27 |             90 | 
|  1878335 | g.00000 | skycell.055 | MD04.refstack.20121125 |     12648.1 |            21 |             89 | 

reftest sample summary --

  • label: MD04.reftest.20131018
  • data_group:
    • MD04.reftest.20131018.c#f<FWHM cut>_1dg where #=comparison run, <FWHM cut>=cut <fwhm_major from camProcessedExp x10>
    • MD04.reftest.20131018.all.<FWHM cut>_1dg where all means 2009--2012 data possible, <FWHM cut>=cut <fwhm_major from camProcessedExp x10>
    • MD04.reftest.20131018.2012.<FWHM cut>_1dg where 2012 means 2012 data only, <FWHM cut>=cut <fwhm_major from camProcessedExp x10>
      | stack_id | filter  | skycell_id  | data_group                        | dtime_stack | reject_images | count(warp_id) |
      +----------+---------+-------------+-----------------------------------+-------------+---------------+----------------+
      |  2858194 | y.00000 | skycell.055 | MD04.reftest.20131018.all.40_1dg  |     5362.59 |             6 |             64 | 
      |  2858195 | y.00000 | skycell.055 | MD04.reftest.20131018.2012.40_1dg |     2466.19 |             3 |             30 | 
      |  2858196 | z.00000 | skycell.055 | MD04.reftest.20131018.all.39_1dg  |     5625.43 |            12 |             70 | 
      |  2858197 | z.00000 | skycell.055 | MD04.reftest.20131018.2012.39_1dg |     775.345 |             2 |             13 | 
      |  2858198 | i.00000 | skycell.055 | MD04.reftest.20131018.all.40_1dg  |     7143.17 |            14 |             86 | 
      |  2858199 | i.00000 | skycell.055 | MD04.reftest.20131018.2012.40_1dg |     824.543 |             4 |             16 | 
      |  2858200 | r.00000 | skycell.055 | MD04.reftest.20131018.all.47_1dg  |     8400.79 |            18 |            105 | 
      |  2858201 | r.00000 | skycell.055 | MD04.reftest.20131018.2012.47_1dg |     697.182 |             1 |             11 | 
      |  2858202 | g.00000 | skycell.055 | MD04.reftest.20131018.all.54_1dg  |       10185 |            21 |            132 | 
      |  2858203 | g.00000 | skycell.055 | MD04.reftest.20131018.2012.54_1dg |     513.141 |             4 |             12 | 
      
      |  2858204 | y.00000 | skycell.055 | MD04.reftest.20131018.c1f38_1dg   |     3047.56 |            10 |             42 | 
      |  2858205 | y.00000 | skycell.055 | MD04.reftest.20131018.c1f40_1dg   |     4244.98 |            10 |             57 | 
      |  2858206 | y.00000 | skycell.055 | MD04.reftest.20131018.c1f42_1dg   |      6810.5 |            12 |             85 | 
      |  2858207 | y.00000 | skycell.055 | MD04.reftest.20131018.c1f45_1dg   |     8526.39 |            22 |            120 | 
      
      |  2858208 | z.00000 | skycell.055 | MD04.reftest.20131018.c1f36_1dg   |     1924.18 |             7 |             29 | 
      |  2858209 | z.00000 | skycell.055 | MD04.reftest.20131018.c1f39_1dg   |     4697.62 |            12 |             60 | 
      |  2858210 | z.00000 | skycell.055 | MD04.reftest.20131018.c1f40_1dg   |     6659.67 |            10 |             76 | 
      |  2858211 | z.00000 | skycell.055 | MD04.reftest.20131018.c1f42_1dg   |     11243.9 |            21 |            128 | 
      
      |  2858212 | i.00000 | skycell.055 | MD04.reftest.20131018.c1f38_1dg   |     3108.67 |             6 |             36 | 
      |  2858213 | i.00000 | skycell.055 | MD04.reftest.20131018.c1f40_1dg   |        4866 |            11 |             65 | 
      |  2858214 | i.00000 | skycell.055 | MD04.reftest.20131018.c1f41_1dg   |     7960.16 |            15 |             94 | 
      |  2858215 | i.00000 | skycell.055 | MD04.reftest.20131018.c1f42_1dg   |     9922.86 |            20 |            116 | 
      
      |  2858216 | r.00000 | skycell.055 | MD04.reftest.20131018.c1f41_1dg   |     2058.32 |             3 |             27 | 
      |  2858217 | r.00000 | skycell.055 | MD04.reftest.20131018.c1f44_1dg   |     4757.83 |             8 |             61 | 
      |  2858218 | r.00000 | skycell.055 | MD04.reftest.20131018.c1f46_1dg   |     7017.18 |            18 |             90 | 
      |  2858219 | r.00000 | skycell.055 | MD04.reftest.20131018.c1f48_1dg   |     8713.09 |            24 |            120 | 
      
      |  2858220 | g.00000 | skycell.055 | MD04.reftest.20131018.c1f46_1dg   |     1985.38 |             6 |             29 | 
      |  2858221 | g.00000 | skycell.055 | MD04.reftest.20131018.c1f50_1dg   |     4995.98 |            13 |             67 | 
      |  2858222 | g.00000 | skycell.055 | MD04.reftest.20131018.c1f51_1dg   |     6537.38 |            16 |             85 | 
      |  2858223 | g.00000 | skycell.055 | MD04.reftest.20131018.c1f53_1dg   |     10531.3 |            25 |            132 |
      
      

IQ Stack

Basically a best reference stack possible sample -- 10-15 inputs, uses MD_DEEP_STACK_1DG so no extra cuts on FWHM made in order to look for bad warps need to be excluded in future runs

  • stack -- label, data_group = MD04.iqtest.20131118
  • staticsky -- label, data_group = MD04.iqtest.20131118.Nx20131122, MD04.iqtest.20131118.<#>x20131122 -- where <#> is the number of filters available
  • skycal -- label, data_group = MD04.iqtest.20131118.Nx20131122.cal20131122, MD04.iqtest.20131118.Nx20131122.cal20131122

Summary --

| stack_id | filter  | skycell_id  | data_group           | dtime_stack | reject_images | count(warp_id) |
+----------+---------+-------------+----------------------+-------------+---------------+----------------+
|  2922630 | g.00000 | skycell.055 | MD04.iqtest.20131118 |     1275.24 |             0 |             16 | 
|  2922690 | r.00000 | skycell.055 | MD04.iqtest.20131118 |     1278.36 |             0 |             17 | 
|  2922750 | i.00000 | skycell.055 | MD04.iqtest.20131118 |     1330.77 |             0 |             16 | 
|  2922810 | z.00000 | skycell.055 | MD04.iqtest.20131118 |     1383.91 |             0 |             17 | 
|  2922870 | y.00000 | skycell.055 | MD04.iqtest.20131118 |     1164.04 |             0 |             16 | 

Reference Yearly/Seasonal Stack

Basically a best reference stack possible using the specific year/season data -- 30-35 warp inputs, used MD_DEEP_STACK_1DG so no extra cuts on FWHM made in order to look for bad warps need to be excluded in future runs

  • stack -- label, data_group = MD04.yr<#>.reftest.20131118, where <#> = 1-4
  • staticsky -- label, data_group = MD04.yr<Y>.reftest.20131118.Nx20131122, MD04.yr<YY>.reftest.20131118.<F>x20131122 -- where <Y> is the year/season, <F> is the number of filters available
  • skycal -- label, data_group = MD04.yr<Y>.reftest.20131118.Nx20131122.cal20131126, MD04.yr<Y>.reftest.20131118.Nx20131122.cal20131126 -- where <Y> is the year/season

Comparison of yr1 and yr2-4 may help show if 2009/2010 season should be excluded from other stack types due to baffle improvements.

Will need to look for any target events lost due to fewer inputs (field edge events) and compare to the MD04.iqtest.20131118.

Summary --

| stack_id | filter  | skycell_id  | data_group                | dtime_stack | reject_images | count(warp_id) |
+----------+---------+-------------+---------------------------+-------------+---------------+----------------+
|  2922939 | g.00000 | skycell.055 | MD04.yr1.reftest.20131118 |     2975.37 |             0 |             36 | 
|  2923010 | r.00000 | skycell.055 | MD04.yr1.reftest.20131118 |     2736.69 |             0 |             35 | 
|  2923083 | i.00000 | skycell.055 | MD04.yr1.reftest.20131118 |     3076.07 |             0 |             36 | 
|  2923153 | z.00000 | skycell.055 | MD04.yr1.reftest.20131118 |     2916.97 |             0 |             38 | 
|  2923224 | y.00000 | skycell.055 | MD04.yr1.reftest.20131118 |     2453.73 |             0 |             34 | 

| stack_id | filter  | skycell_id  | data_group                | dtime_stack | reject_images | count(warp_id) |
+----------+---------+-------------+---------------------------+-------------+---------------+----------------+
|  2923285 | g.00000 | skycell.055 | MD04.yr2.reftest.20131118 |     3151.29 |             0 |             34 | 
|  2923345 | r.00000 | skycell.055 | MD04.yr2.reftest.20131118 |     2884.53 |             0 |             32 | 
|  2923405 | i.00000 | skycell.055 | MD04.yr2.reftest.20131118 |     2641.78 |             0 |             33 | 
|  2923465 | z.00000 | skycell.055 | MD04.yr2.reftest.20131118 |     2575.58 |             0 |             31 | 
|  2923525 | y.00000 | skycell.055 | MD04.yr2.reftest.20131118 |      1164.6 |             0 |             16 | 

| stack_id | filter  | skycell_id  | data_group                | dtime_stack | reject_images | count(warp_id) |
+----------+---------+-------------+---------------------------+-------------+---------------+----------------+
|  2923585 | g.00000 | skycell.055 | MD04.yr3.reftest.20131118 |     2455.01 |             0 |             30 | 
|  2923645 | r.00000 | skycell.055 | MD04.yr3.reftest.20131118 |     2369.33 |             0 |             31 | 
|  2923705 | i.00000 | skycell.055 | MD04.yr3.reftest.20131118 |     3532.97 |             0 |             36 | 
|  2923765 | z.00000 | skycell.055 | MD04.yr3.reftest.20131118 |     3856.02 |             0 |             33 | 
|  2923825 | y.00000 | skycell.055 | MD04.yr3.reftest.20131118 |     2659.61 |             0 |             34 | 

| stack_id | filter  | skycell_id  | data_group                | dtime_stack | reject_images | count(warp_id) |
+----------+---------+-------------+---------------------------+-------------+---------------+----------------+
|  2923885 | g.00000 | skycell.055 | MD04.yr4.reftest.20131118 |     2684.92 |             0 |             32 | 
|  2923945 | r.00000 | skycell.055 | MD04.yr4.reftest.20131118 |     2965.94 |             0 |             32 | 
|  2924005 | i.00000 | skycell.055 | MD04.yr4.reftest.20131118 |      2862.7 |             0 |             36 | 
|  2924065 | z.00000 | skycell.055 | MD04.yr4.reftest.20131118 |     2646.95 |             0 |             32 | 
|  2924125 | y.00000 | skycell.055 | MD04.yr4.reftest.20131118 |     2798.91 |             0 |             34 | 

Reference Out-of-Season Stack

Two types possible, all exposures before season (default during active observing) and all exposures before/after season (or N seasons)

Not currently possible with current tools unless add warp_id manually for each skycell -- this method may even be necessary for all final stacks to exclude problematic exposure/chips in a warp skycell

Sample test for not season 2 and 3, only including warps from season 1 (2009/2010) and 4 (2012/2013)

  • stack -- label, data_group = MD04.notyr2_3.reftest.20131124
  • staticsky -- label, data_group = MD04.notyr2_3.reftest.20131124.Nx20131124, MD04.notyr2_3.reftest.20131124.<F>x20131124 -- where <F> is the number of filters available
  • skycal -- label, data_group = MD04.notyr2_3.reftest.20131124.Nx20131124.cal20131126, MD04.notyr2_3.reftest.20131124.Nx20131124.cal20131126

Summary --

| stack_id | filter  | skycell_id  | data_group                     | dtime_stack | reject_images | count(warp_id) |
+----------+---------+-------------+--------------------------------+-------------+---------------+----------------+
|  2937260 | g.00000 | skycell.055 | MD04.notyr2_3.reftest.20131124 |      2367.6 |             0 |             32 | 
|  2937325 | r.00000 | skycell.055 | MD04.notyr2_3.reftest.20131124 |     2475.84 |             0 |             32 | 
|  2937385 | i.00000 | skycell.055 | MD04.notyr2_3.reftest.20131124 |     2495.37 |             0 |             31 | 
|  2937448 | z.00000 | skycell.055 | MD04.notyr2_3.reftest.20131124 |     3058.89 |             0 |             37 | 
|  2937511 | y.00000 | skycell.055 | MD04.notyr2_3.reftest.20131124 |     2247.48 |             0 |             30 |