Initial/ongoing comparison of IPP to SuperMACHO/ESSENCE Photpipe reductions using GPC1 images from the 2009.01 release.

[[PageOutline]]

*** '''WARNING''' work in progress with irregular documentation - please highlight areas lacking ***
(don't forget to note the IPP svn-head revision # being used if different from a default, say 23243 post-pswarp fix?
also test runs with all default set options, need to address this somehow as well.)
*** note - many skycell.02877.31 typo'd as 02879.31 including data runs, need to fix

== Action Items ==

* Come up with explicit test and verification of pswarp vs. SWarp, ppStack vs SWarp:stack, and apply to idealized(simtest) and real(release.200901) images 
  - Does the output match what we think it should?
  - Verify that the actual estimation of the "best" lowest-common-denominator PSF is correct.
  - Verify that the generated convolved stacked image matches what it theoretically should.
  - Verify that the covariance map is correct (Paul has already done this).

* Learn what ppSub expects the stacked image to have so we can pass information from SWarp-generated co-add as if it were an IPP-generated coadd
  - ppSub needs SWarp:stack put into skycell dimensions (use swarp_single to setup each img first)
  - ppSub needs SWarp:stack x flipped 
  - ppSub needs wt and mk from SWarp:stack 

* Develop non-convolving version of ppStack. [DONE]
  - Median and un/weighted sums/average combines that keeps track of the image, variance, mask, and PSF information appropriately.

  * Paul has added the non-convolving feature to ppStack on trunk, r23602. The behaviour is controlled by the PPSTACK recipe:
    * Convolution can be turned off by setting CONVOLVE = FALSE.
    * All rejection can be turned off by setting ITER = 0.
    * ie in recipes/ppStack.config: CONVOLVE BOOL FALSE
    * note: gpc1/ppStack.config is not set to be used in gpc1/camera.config by default, will need to add if making modifications there.


* Run psphot on ppSub images from non-convolved stacks.  Compare results with release.200901 diffim.  Search for transients in new subtractions.
  - Prototype with Rest/Huber hacks on 10 skycells
  - Full run once re-developed ppStack exists

== IPP revisions and configurations ==
Need to be sure to note version used in FITS header

|| r23243 || post-pswarp fix, much of the groucho data originally run with this ||
|| r23602 || post-ppStack convolution option fix ||
|| ~ r23685 || pswarp improvement to avoid arbitrary cut of X% pixels ||
|| r23739 || keep eye out for problems after inverse option added (A-B,B-A) ||
|| r23769 || psphot improvements, issue with r23790 as PP adding substamp outputs in ppSub ||
|| || ||

* turn off compression on images: modify recipes/filerules-split.mdc
 
 change the *OUTPUT* FITS.TYPE from COMP_ to NONE?

* turn off convolution in ppStack: add to recipes/ppStack.config

{{{
 CONVOLVE        BOOL    FALSE
 ITER            S32     0
}}}

* to save intermediate images in ppStack, add to the ppStack call: -recipe PPSTACK DEBUG

* get coords for stamps (all/used/rejected) and lesson on -trace -- summarized ticket #1218
{{{
build the IPP without optimisation (no '-optimize' flag to psbuild)

To get all trace messages regarding PSF-matching, try adding "-trace psModules.imcombine 10" (LOT of extra information)
You can direct this to a file using "-tracedest /path/to/file". For more fine-grained control, try:

#1 rejected stamps: "-trace psModules.imcombine.pmSubtractionRejectStamps 3" (see pmSubtraction.c:866)
#2 used stamps: "-trace psModules.imcombine.pmSubtractionCalculateDeviations 5" (see pmSubtractionEquation.c:1111)
#3 all possible stamps: "-trace psModules.imcombine.pmSubtractionStampsSet 9" (see pmSubtractionStamps.c:438)
}}}

== Basic Data Locations ==
All original data (particularly the chip level used to remake the warps) is from the release.200901 set, processed with ipp-magic-v0.

* Odyessy:
 original data (36475/chip,warp,diff etc): /n/data1/panstarrs/GPC1/MD01_20081002_1/

 workspace (swarps, ipp-rerun collections): /n/data1/panstarrs/GPC1/MD01_20081002_1/workspace/

 IPP source: /odyssey/apps/src/ipp/ipp-SVN_top

* groucho: 

 photpipe runs (pswarp hotpants etc): http://groucho.pha.jhu.edu/PS/smsn_workspace/

 local ipp re-runs (uploaded to odyssey): http://groucho.pha.jhu.edu/PS/ippsvn/ipp_vs_photpipe/

* huevon:

 /data/pipe/v10.0/PS1/GPC1/workspace/

 http://huevon.physics.harvard.edu/arest/smsn_GPC1/
 
 http://huevon.physics.harvard.edu/arest/ipp_GPC1/

== Chip/OTA positions pre/post-refub ==

Consider the OTAs not being replaced/changed during refurbishment. From JT - This was the original configuration (with the OTA ID numbers).  Note that this is OTA0y -> OTA7y, i.e. looking down on the FP and reversed from standard camera coords. 

{{{
 OTA07   2-16b2  1-18b1  2-10b2  1-21a1  2-12a1  2-16a2  OTA77	
 2-09b1  1-09b1  2-07b2  1-02b1  2-07a2  1-10a2  2-14a1  2-05a1
 2-16b1  2-13b1  1-01b1  1-01b2  1-21a2  2-12a2  2-17a2  2-16a1
 1-09a1  1-09a2  2-04a2  2-04a1  1-05a2  2-09a1  1-07a1  1-18a1
 1-13a1  1-14a2  1-02a1  1-02a2  2-11b2  1-05b2  1-02b2  1-25b1
 2-18a1  1-17a2  2-09a2  2-01a2  1-07b1  1-04b1  2-23b2  1-14b1
 2-15a1  2-13a2  1-04a1  2-01a1  1-03b2  2-09b2  1-03b1  2-04b2
 OTA00   2-22a1  1-19a1  1-12a1  2-18b2  1-17b1  2-22b2  OTA70
}}}

This is the present configuration.  As you can see a lot of devices have moved as well as being replaced.

{{{
 OTA07   2-07b2 [1-18b1] 2-23b2  3-14a1  3-11a2  1-09a1  OTA77
[2-09b1] 1-02b1  3-04b1  3-14b1  3-02a1  3-02a2  3-10a1 [2-05a1]
 1-09b1  3-08b1  3-05b2 [1-01b2] 3-04a2  3-03a2  2-04a1  1-21a1
 1-01b1  3-05b1  3-10b1  3-03b1  2-09a1  3-05a2 [1-07a1] 1-05a2
 2-04a2  3-06a1 [1-02a1][1-02a2][2-11b2][1-05b2][1-02b2] 3-11b2
[2-18a1] 3-11a1 [2-09a2][2-01a2][1-07b1][1-04b1] 3-03b2 [1-14b1]
[2-15a1] 2-01a1 [1-04a1] 3-05a1 [1-03b2][2-09b2][1-03b1] 2-13b1
 OTA00  [2-22a1][1-19a1] 3-13a2  3-02b1 [1-17b1] 2-16b1  OTA70
}}}

* The ones not being changed are in bold (would like them to be... for now []).  A sample that may be good to work with are XY
 * 33,43 - center tent
 * 32,42,52,53,41 - ring of better focus
 * 50,61,72 - outer edge


== Skycell Layout ==

Same as default used for 2009.01 release and Basic_GPC1_Processing_for_End_Users
{{{
skycells -level 5 -scale 0.2 -nx 8 -ny 8 -fix-ns -padding 0.05 -D CATDIR default 
}}}

 `odyssey:/n/data1/panstarrs/GPC1/MD01_20081002_1/workspace/su_pswarp/default`


== Test chip/skycell Sample ==

Hard to clearly tell where chip boundaries are, so make skycell set with just single chip. (these chosen because they are already mostly made and provide a center vs better focus ring sample). May need another sample skycell near chip XY22 with other chip contribution to test multi-chip psf matching for ppStack/ppSub.

* '''skycell.02877.31''': focus on chip XY33 (~50% of field, center of camera and tent region)
{{{
 RA,Dec center = 02:23:37.250 -04:20:03.53 J2000  35.90521  -4.33431
 X,Y center (size) = 2754, 2003 (5508,4006)
 platescale = 0.2
}}}

 i.e. `odyssey:/n/data1/panstarrs/GPC1/MD01_20081002_1/36475/chip/o4741g0234o.36475.ch.11948.XY33.ch.fits`

 [[Image(chip33_all_skycell02877.31.ds9.jpg,120px)]]
 [[Image(chip33_skycell02877.31.ds9.jpg,120px)]]
 [[Image(zoom_chip33_skycell02877.31.ds9.jpg,120px)]]

 || MD01_20081002_1/36475/chip/o4741g0234o.36475.ch.11948.XY33.ch.fits
 || MD01_20081002_1/36476/chip/o4741g0235o.36476.ch.11949.XY33.ch.fits
 || MD01_20081002_1/36477/chip/o4741g0236o.36477.ch.11950.XY33.ch.fits
 || MD01_20081002_1/36478/chip/o4741g0237o.36478.ch.11951.XY33.ch.fits
 || MD01_20081002_1/36479/chip/o4741g0238o.36479.ch.11952.XY33.ch.fits
 || MD01_20081002_1/36480/chip/o4741g0240o.36480.ch.11953.XY33.ch.fits
 || MD01_20081002_1/36481/chip/o4741g0241o.36481.ch.11954.XY33.ch.fits
 || MD01_20081002_1/36625/chip/o4741g0239o.36625.ch.11955.XY33.ch.fits



* '''skycell.02879.00''':  focus on chip XY22 (~90% of field, in better focus region around center. Note, large masked region in o4741g0237o.36478.ch.11951.XY22.ch.fits)
{{{
 RA,Dec center = 02:24:32.502 -05:00:06.29 J2000  36.09376  -5.00175
 X,Y center (size) = 2754, 2003 (5508,4006)
 platescale = 0.2
}}}

 i.e. `odyssey:/n/data1/panstarrs/GPC1/MD01_20081002_1/36475/chip/o4741g0234o.36475.ch.11948.XY22.ch.fits`

 [[Image(chip22_skycell02879.00.ds9.jpg,120px)]]
 [[Image(zoom_chip22_skycell02879.00.ds9.jpg,120px)]]

 || MD01_20081002_1/36475/chip/o4741g0234o.36475.ch.11948.XY22.ch.fits
 || MD01_20081002_1/36476/chip/o4741g0235o.36476.ch.11949.XY22.ch.fits
 || MD01_20081002_1/36477/chip/o4741g0236o.36477.ch.11950.XY22.ch.fits
 || MD01_20081002_1/36478/chip/o4741g0237o.36478.ch.11951.XY22.ch.fits
 || MD01_20081002_1/36479/chip/o4741g0238o.36479.ch.11952.XY22.ch.fits
 || MD01_20081002_1/36480/chip/o4741g0240o.36480.ch.11953.XY22.ch.fits
 || MD01_20081002_1/36481/chip/o4741g0241o.36481.ch.11954.XY22.ch.fits
 || MD01_20081002_1/36625/chip/o4741g0239o.36625.ch.11955.XY22.ch.fits



== Warping (PSwarp and SWarp) ==

Check SWarp and PSWarp on the same chip to same skycell level using the following 

* pswarps - uncompressed, single FITS files in 

 `odyssey:/n/data1/panstarrs/GPC1/MD01_20081002_1/workspace/su_pswarp/`

 http://groucho.pha.jhu.edu/PS/ippsvn/ipp_vs_photpipe/su_pswarps/

{{{
fully filled skycell: suo4741g0234o.36475.wrp.11948.all.skycell.02877.31.fits
                      suo4741g0234o.36475.wrp.11948.all.skycell.02879.00.fits
single chip skycell: suo4741g0234o.36475.wrp.11948.XY33.skycell.02877.31.fits
                     suo4741g0234o.36475.wrp.11948.XY22.skycell.02879.00.fits
}}}

* SWarps: Armin notes some correlated noise still apparent in pswarps compared to SWarps Lanczos3 vs 4?

* need to look at Lanczos w/ masking and boundaries - mask growth still different w/ photpipe
 * photpipe has good/ok/poor set, pswarp?


== Differencing (ppSub and hotpants) ==

Look at before stacking to since ppStack uses the same routines/convolution as ppSub. Try ppSub on both pswarps and swarps to look at differences.

'''ppSub''' [pswarp-pswarp, swarp-swarp, pswarp-ppStack, pswarp-swarp:stack]

 `odyssey:/n/data1/panstarrs/GPC1/MD01_20081002_1/workspace/meh_ppsub_pswarp/`
 
 http://groucho.pha.jhu.edu/PS/ippsvn/ipp_vs_photpipe/ppsub_pswarps/


* ppSub [pswarp X-pswarp Y]: 

 `odyssey: subisis.36478_36476.wrp.XY33.skycell.02877.31.fits`
 
 http://groucho.pha.jhu.edu/PS/ippsvn/ipp_vs_photpipe/ppsub_pswarps/subisis.36478_36476.wrp.XY33.skycell.02877.31.fits

 `odyssey: subisis.36477_36475.wrp.XY22.skycell.02879.00.fits`
 
 http://groucho.pha.jhu.edu/PS/ippsvn/ipp_vs_photpipe/ppsub_pswarps/subisis.36477_36475.wrp.XY22.skycell.02879.00.fits

 || xy33/02877.31, default settings, 36478-36476 || xy22/02879.00, default settings, 36477-36475 ||
 || [[Image(subisis_default_n78_r76_xy33_skycell02877.31.ds9.jpg,120px)]] ||   [[Image(subisis_default_n77_r75_xy22_skycell02879.00.ds9.jpg,120px)]] ||

 Several other combos for skycell.02877.31, XY33 - 

 || subisis.36475_36479.wrp.XY33.skycell.02877.31.fits ||
 || subisis.36477_36478.wrp.XY33.skycell.02877.31.fits ||
 || subisis.36477_36479.wrp.XY33.skycell.02877.31.fits || tetris pattern ||
 || subisis.36479_36477.wrp.XY33.skycell.02877.31.fits || reversed no pattern ||
 || subisis.36477_36625.wrp.XY33.skycell.02877.31.fits ||
 || subisis.36478_36476.wrp.XY33.skycell.02877.31.fits ||
 || subisis.36478_36478.wrp.XY33.skycell.02877.31.fits ||
 || subisis.36480_36479.wrp.XY33.skycell.02877.31.fits ||
 || subisis.36481_36479.wrp.XY33.skycell.02877.31.fits ||
 || subisis.36625_36477.wrp.XY33.skycell.02877.31.fits ||
 || subisis.36625_36479.wrp.XY33.skycell.02877.31.fits ||
 || 36481-36475 || can't make stamps ||
 || 36477-36475 || can't make stamps ||


* ppSub [pswarp X-pswarp X]: silly test subtracting a warp from itself, 36478 for skycell.02877.31, XY33

 `odyssey: subisis.36478_36478.wrp.XY33.skycell.02877.31.fits`
 
 http://groucho.pha.jhu.edu/PS/ippsvn/ipp_vs_photpipe/ppsub_pswarps/subisis.36478_36478.wrp.XY33.skycell.02877.31.fits

 || 36478-36478 for skycell.02877.31 XY33 ||
 || [[Image(subisis_default_n78_r78_xy33_skycell02877.31.ds9.jpg,120px)]] ||

* ppSub [pswarp X-pswarp Y]: fine-tune using 36625-36476

 PSWarps:

 `odyssey:/n/data1/panstarrs/GPC1/MD01_20081002_1/workspace/`
 `su_pswarp/suo4741g0239o.36625.wrp.11955.XY33.skycell.02877.31.fits`
 `su_pswarp/suo4741g0235o.36476.wrp.11949.XY33.skycell.02877.31.fits`

 ppSub: initial run

 `odyssey:/n/data1/panstarrs/GPC1/MD01_20081002_1/workspace/`
 `meh_ppsub_pswarps/subisis.36625_36476.wrp.XY33.skycell.02877.31.fits`

 (sub)Stamp locations (from hardcoded log output line in  pmSubtractionStamps.c : pmSubtractionStampsExtract):
{{{
2009-04-04 18:22:10Z|groucho             |I|psModules.imcombine
    Found 9 stamps
2009-04-04 18:22:10Z|groucho             |I|psModules.imcombine
    Stamp 70 (x y size): 2890 469 35
2009-04-04 18:22:10Z|groucho             |I|psModules.imcombine
    Stamp 94 (x y size): 2071 757 35
2009-04-04 18:22:10Z|groucho             |I|psModules.imcombine
    Stamp 99 (x y size): 3033 711 35
2009-04-04 18:22:10Z|groucho             |I|psModules.imcombine
    Stamp 127 (x y size): 2960 770 35
2009-04-04 18:22:10Z|groucho             |I|psModules.imcombine
    Stamp 137 (x y size): 5054 842 35
2009-04-04 18:22:10Z|groucho             |I|psModules.imcombine
    Stamp 146 (x y size): 1357 1040 35
2009-04-04 18:22:10Z|groucho             |I|psModules.imcombine
    Stamp 152 (x y size): 2419 1060 35
2009-04-04 18:22:10Z|groucho             |I|psModules.imcombine
    Stamp 174 (x y size): 1356 1211 35
2009-04-04 18:22:10Z|groucho             |I|psModules.imcombine
    Stamp 375 (x y size): 2238 2668 35
}}}

 Using the -trace options, the used stamps are 
{{{
     Deviation for stamp 70 (2890,469): 11.264941
     Deviation for stamp 94 (2071,757): 0.347919
     Deviation for stamp 99 (3033,711): 0.488299
     Deviation for stamp 127 (2960,770): 1.450665
     Deviation for stamp 137 (5054,842): 0.512477
     Deviation for stamp 146 (1357,1040): 3.164435
     Deviation for stamp 152 (2419,1060): 0.342433
     Deviation for stamp 174 (1356,1211): 0.387268
     Deviation for stamp 375 (2238,2668): 2.065892
     Deviation for stamp 94 (2071,757): 0.433812
     Deviation for stamp 99 (3033,711): 0.527412
     Deviation for stamp 127 (2960,770): 1.187763
     Deviation for stamp 137 (5054,842): 0.543499
     Deviation for stamp 146 (1357,1040): 2.416274
     Deviation for stamp 152 (2419,1060): 0.405773
     Deviation for stamp 174 (1356,1211): 0.495316
     Deviation for stamp 375 (2238,2668): 1.806831

 A stamp is rejected in a later cut(?) to reduce the RMD deviation from 3.10 to 0.97
   Rejecting stamp 70 (2890,469)

 There are some 234 tested, so 225 rejected initially because of bad mask.

}}}


* ppSub [SWarp-SWarp:stack]: need SWarp'd images and setup for what ppSub needs added to header, weight/mask files

* ppSub [pswarp-SWarp:stack]: need to learn what ppSub needs added to header, how to make weight and mask files

* ppSub [pswarp-ppStack]: need to setup ppStack w+w/o convolution 

 http://groucho.pha.jhu.edu/PS/ippsvn/ipp_vs_photpipe/ppsub_pswarp_ppstack

 * w/o convolution - quick run example of epoch 36479-stk4a

 http://groucho.pha.jhu.edu/PS/ippsvn/ipp_vs_photpipe/ppsub_pswarp_ppstack_r23602/

 http://groucho.pha.jhu.edu/PS/ippsvn/ipp_vs_photpipe/ppsub_pswarp_ppstack_r23602/subisis.36479_stk4a.skycell.02877.31.XY33.fits

 * look at psf used for convolution - FITS table w/ coefficients and desc. of image residuals
 http://groucho.pha.jhu.edu/PS/ippsvn/ipp_vs_photpipe/ppsub_pswarp_ppstack/skycell.02877.31.XY33.stk4a.target.psf


'''Hotpants''' [pswarp-pswarp, pswarp-ppStack, swarp-ppStack, pswarp-SWarp:stack]

* hotpants [pswarp-pswarp]: 
 * pswarps converted to photpipe format with convert_ippwarp2smsn.pl
  * fields are skycell.02877.31==f02877_31, 02879.00==f02879.00

 * test sample, fully filled skycells - will make amp 0 soon

 36480-36475: http://groucho.pha.jhu.edu/PS/smsn_workspace/ut081002meh80_ut081002meh75/1/

 36478-36475: http://groucho.pha.jhu.edu/PS/smsn_workspace/ut081002meh78_ut081002meh75/1/

 * for skycells with single chip only inclusion, amp=chip#

 36477-36475: http://groucho.pha.jhu.edu/PS/smsn_workspace/ut081002meh77_ut081002meh75/22

 || 36475_XY22 || 36477_XY22 ||
 || ppSub || hotpants ||
 [[Image(ppsub_hotpants_default_n77_r75_xy22_skycell02879.00.ds9.jpg,120px)]] 

 36625-26479: http://groucho.pha.jhu.edu/PS/smsn_workspace/ut081002meh625_ut081002meh79/33

 || 36479_XY33 || 36625_XY33 ||
 || hotpants || ppSub ||
 [[Image(ppsub_hotpants_default_n625_r79_xy33_skycell02877.31.ds9.jpg,120px)]]

 Initial sample test runs --

 http://groucho.pha.jhu.edu/PS/smsn_workspace/

 || 36476-36475 || FSIG=8.752>2.000 || 2/2 stamps ||
 || 36477-36475 || FSIG=66.432>2.000 || 3/3 ||
 || 36477-36476 || FSIG=6.761>2.000 || 2/2 || 
 || 36478-36477 || FSCAT=1.1428>0.5000 || 7/7 ||
 || || ||
 || 36479-36475 || FSIG=8.802>2.000 || 7/8 || 
 || 36479-36476 || FSIG=10.883>2.000 || 7/8 ||
 || 36479-36477 || FSIG=4.941>2.000 || 7/8 ||
 || 36479-36478 || OK? some, parts of diffim look poor || 7/8 || 
 || || ||
 || 36625-36475 || FSIG=118.939>2.000 || 4/4 ||
 || 36625-36476 || FSCAT=0.9864>0.5000 || 7/8 ||
 || 36625-36477 || FSIG=50.972>2.000 || 6/6  ||
 || 36625-36478 || FSCAT=2.0320>0.5000 || 6/6 ||
 || 36625-36479 || OK || 7/7 ||
 || || ||
 || 36480-36475 || FSIG=33.954>2.000 || 3/3 ||
 || 36480-36476 || FSIG=51.571>2.000 || 3/3 ||
 || 36480-36479 || FSIG=28.267>2.000 || 5/5 ||
 || 36481-36479 || OK || 7/8 ||
 || 36481-36480 || OK || 7/8 ||


* hotpants [pswarp X-pswarp Y]: fine-tune using 36625-36476

 http://groucho.pha.jhu.edu/PS/smsn_workspace/ut081002meh625_ut081002meh76/33/


* hotpants [pswarp-ppStack]: convert the ppStack file like a pswarp using convert_ippwarp2smsn.pl

 http://groucho.pha.jhu.edu/PS/smsn_workspace/ut081002meh80_ut081002mehstk4a/1/


* hotpants [swarp-ppStack]: need to setup ppStack w+w/o convolution


* hotpants [pswarp-SWarp:stack]: need photpipe setup for SWarp:stack, weight/mask files?


== Stacking (ppStack and SWarp) ==

'''ppStack''' [pswarps,pswarps:diffims, SWarps,SWarp:diffims]

 `odyssey:/n/data1/panstarrs/GPC1/MD01_20081002_1/workspace/meh_ppstack`
 
 http://groucho.pha.jhu.edu/PS/ippsvn/ipp_vs_photpipe/ppstack_pswarps/

* ppStack [pswarps WXYZ]: test combination 1st 4 imgs(a), 2nd 4(b), all 8

 || skycell.02879.31.stk4a.fits ||
 || skycell.02879.31.stk4b.fits ||
 || skycell.02879.31.stk8.fits ||

[[Image(stktest.ds9.jpg,120px)]]


* ppStack [pswarps XXXX]: silly test combining 4 copies of the same image

 `odyssey: skycell.02879.31.stk4a78.fits`

 http://groucho.pha.jhu.edu/PS/ippsvn/ipp_vs_photpipe/ppstack_pswarps/skycell.02879.31.stk4a78.fits

 [[Image(skycell.02879.31.stk4a78.ds9.jpg,120px)]]

* ppStack [pswarps, r23602 test]: dupe of previous set test with convolution turned off

 http://groucho.pha.jhu.edu/PS/ippsvn/ipp_vs_photpipe/ppsub_pswarp_ppstack_r23602/

 http://groucho.pha.jhu.edu/PS/ippsvn/ipp_vs_photpipe/ppsub_pswarp_ppstack_r23602/skycell.02877.31.XY33.stk4a.fits

 http://groucho.pha.jhu.edu/PS/ippsvn/ipp_vs_photpipe/ppsub_pswarp_ppstack_r23602/skycell.02877.31.XY33.stk4a78.fits
  

  * silly test: suo4741g0237o.36478.wrp.11951.XY33.skycell.02877.31.fits - (skycell.02877.31.XY33.stk4a78.fits/4.0)

 http://groucho.pha.jhu.edu/PS/ippsvn/ipp_vs_photpipe/ppsub_pswarp_ppstack_r23602/test.36478-skycell.02877.31.XY33.stk4a78.fits
  
   Similar to the pswarp-pswarp case.


* ppStack [SWarps]: SWarp skycells need to be built

* ppStack [ppSubs:pswarps]: stack of ppSub:pswarp diffims (and for SWarps)

* ppStack [hotpants:pswarps]: stack of hotpants:pswarp diffims (and for SWarps)


'''SWarp''' [pswarps, pswarp:diffims, SWarps,SWarp:diffims]

* SWarp [pswarps]: more info needed from Armin and needs to be resampled into skycell boundaries.

 `odyssey:/n/data1/panstarrs/GPC1/MD01_20081002_1/workspace/MYDIR`

 http://groucho.pha.jhu.edu/PS/ippsvn/ipp_vs_photpipe/swarpstk_pswarps

 Note that a set of 4 okay, doing set of has an odd result (in what way Armin?)

* SWarp [SWarps]: SWarp skycells need to be built

* SWarp [ppSub:pswarps]: SWarp stack of ppSub:pswarp diffims (and for SWarps)



Initial stack test stats --
{{{
PSWarp
 imstats suo4741g0234o.36475.wrp.11948.all.skycell.02877.31.fits -box 4948 2331 80 80
 STATS: bkg:1.608557 sigma_bkg:24.464390 Nused:6382 Nclipped:18 lowest_val:-85.380753 highest_val:88.475899

SWarp:stack
 imstats 36475-36478.fits -box 4948 2331 80 80
 STATS: bkg:-0.985735 sigma_bkg:12.967378 Nused:6382 Nclipped:18 lowest_val:-52.730080 highest_val:55.483871

ppStack first 4
 imstats skycell.02879.31.stk4a.fits -box 4948 2331 80 80
 STATS: bkg:-0.324937 sigma_bkg:2.308740 Nused:6225 Nclipped:175 lowest_val:-5.852595 highest_val:10.131362

ppStack last 4
 imstats skycell.02879.31.stk4b.fits -box 4948 2331 80 80
 STATS: bkg:-7.943039 sigma_bkg:6.104924 Nused:6304 Nclipped:6 lowest_val:-25.602289 highest_val:12.588408

ppStack all 8
 imstats skycell.02879.31.stk8.fits -box 4948 2331 80 80
 STATS: bkg:-2.741246 sigma_bkg:5.151842 Nused:6216 Nclipped:184 lowest_val:-17.833969 highest_val:21.132132

 -----

PSWarp, new position 
 imstats suo4741g0234o.36475.wrp.11948.all.skycell.02877.31.fits -box 2650 2180 80 80
 STATS: bkg:4.214430 sigma_bkg:28.889675 Nused:3614 Nclipped:149 lowest_val:-4511.794922 highest_val:4742.093750

ppStack all 8, new position
 imstats skycell.02879.31.stk8.fits -box 2650 2180 80 80
 STATS: bkg:0.223009 sigma_bkg:5.304905 Nused:6396 Nclipped:4 lowest_val:-14.943192 highest_val:17.027874
}}}


== Photometry (psphot and dophot/pixie) ==

* psphot: -visual option.  MEH after ~17/18 continues fails with seg fault, right after plotting peak as mag vs cr_n_sig and ext_n_sig