| | 83 | |
| | 84 | ''' WSdiff processing ''' by JRF |
| | 85 | |
| | 86 | Added the OSS.WS.nightlyscience label back into PS2 processing early in the morning, after everything else had finished. The majority of diffs cleared no problem. There was one that I thought was stuck: |
| | 87 | {{{ |
| | 88 | 0 ippx021 BUSY 9195.68 0.0.1.8bd5 0 diff_skycell.pl --threads @MAX_THREADS@ --diff_id 50495 --skycell_id skycell.1370.020 --diff_skyfile_id 2720978 --outroot neb://ipp137.1/gpc2/OSS.WS.nightlyscience/2019/05/20/RINGS.V3/skycell.1370.020/RINGS.V3.skycell.1370.020.WS.dif.50495 --redirect-output --run-state new --reduction SWEETSPOT_WS --dbname gpc2 --verbose |
| | 89 | }}} |
| | 90 | It's been there for a good 2 and half hours. |
| | 91 | |
| | 92 | I found some info about WSdiffs being stuck in an older log: |
| | 93 | {{{ |
| | 94 | past log : http://svn.pan-starrs.ifa.hawaii.edu/trac/ipp/wiki/PS1_IPP_Czarlog_20180618 |
| | 95 | -- turn diff.revert.off, kill ppSub job, wait to fault and manually set quality 42 |
| | 96 | difftool -dbname gpc2 -updatediffskyfile -fault 0 -set_quality 42 -diff_id 50495 -skycell_id skycell.1370.020 |
| | 97 | }}} |
| | 98 | The job actually finished while I was looking this up... so nevermind. |
| | 99 | |
| | 100 | To ensure that we don't have to manually intervene during the night, by removing the label again, I've changed the {{{~ippps2/stdscience/input}}} file so that all WS labels are not automatically added on restart. This means that the label can be re-addedto pantasks at our leisure in the morning. |
| | 101 | |
| | 102 | |
| | 103 | ''' plots ''' by JRF |
| | 104 | |
| | 105 | czarplot.pl can be run manually to generate some plots, for example: |
| | 106 | {{{ |
| | 107 | cd ~ippitc/src/ippMonitor/czartool ; ./czarplot.pl -rt -l OSS.WS.nightlyscience -p '12 hour' |
| | 108 | }}} |
| | 109 | This one gives you a **r**ate plot and a cumulative **t**ime of the past 12 hours for the specific label. This is useful as at present automatic generation of these plots is disabled (they were barely ever used and greatly slowed down IPPmonitor). |
| | 110 | |
| | 111 | |
| | 112 | ''' Checking out diffs ''' |
| | 113 | |
| | 114 | To view the quality of the resulting diffs in ds9 try the following query (from MEH): |
| | 115 | {{{ |
| | 116 | mysql -hscidbm -uippuser -pippuser gpc2 -B -e ' select path_base,diffSkyfile.skycell_id,sources from diffRun join diffSkyfile using (diff_id) join diffInputSkyfile using (diff_id) join warpRun on (diffInputSkyfile.warp1=warpRun.warp_id) join fakeRun using (fake_id) join camRun using (cam_id) join chipRun using (chip_id) join rawExp using (exp_id) where exp_name="o8507h0371o" and quality=0 and diff_mode=1 group by diffSkyfile.skycell_id limit 5; ' | grep -v path_base | awk '{print "neb-ls -p "$1".fits"}' | bash | xargs ~/bin/ds9 -zscale -prefs nancolor orange -geometry 900x800 -frame last -zoom 0.3 -match frame -tile grid gap 3 |
| | 117 | }}} |
| | 118 | Make sure to change the exposure name to the one you desire, and the limit to how many you want to view at once (also note the ds9 path) |
| | 119 | |
| | 120 | For last night, one of the first visits from the first chunk: |
| | 121 | o8624h0111o |
| | 122 | one of the second vists of the second chunk: |
| | 123 | o8624h0181o |
| | 124 | |
| | 125 | The first is a very poor diff where there is a huge overlap in the diffs -> suggesting that the pointing is off. This is also evident when looking at the number of sources detected by the diffs, where we see >10000 deteections; normally less than 1000 are detected. |
| | 126 | |
| | 127 | An example query that will show you how many sources were detected and the fraction of good pixels in the warp and diff stage. Looking at a percentage difference in the good_frac gives an indication on the quality: larger -ve frac -> a worse diff |
| | 128 | {{{ |
| | 129 | SELECT |
| | 130 | diffSkyfile.skycell_id, |
| | 131 | warpSkyfile.skycell_id, |
| | 132 | sources, |
| | 133 | diffSkyfile.good_frac difffrac, |
| | 134 | warpSkyfile.good_frac warpfrac, |
| | 135 | diffSkyfile.good_frac - warpSkyfile.good_frac / warpSkyfile.good_frac as good_stuff |
| | 136 | FROM |
| | 137 | diffRun |
| | 138 | JOIN diffSkyfile USING (diff_id) |
| | 139 | JOIN diffInputSkyfile USING (diff_id) |
| | 140 | JOIN warpRun on (diffInputSkyfile.warp1 = warpRun.warp_id) |
| | 141 | JOIN fakeRun USING (fake_id) |
| | 142 | JOIN camRun USING (cam_id) |
| | 143 | JOIN chipRun USING (chip_id) |
| | 144 | JOIN rawExp USING (exp_id) |
| | 145 | JOIN warpSkyfile USING (warp_id) |
| | 146 | WHERE |
| | 147 | exp_name="o8507h0371o" and diff_mode=1 |
| | 148 | AND |
| | 149 | warpSkyfile.skycell_id = diffSkyfile.skycell_id |
| | 150 | AND |
| | 151 | warpSkyfile.good_frac > 0.5 |
| | 152 | GROUP BY |
| | 153 | diffSkyfile.skycell_id |
| | 154 | ORDER BY |
| | 155 | good_stuff DESC |
| | 156 | LIMIT 20; |
| | 157 | }}} |
| | 158 | |
| | 159 | As a on liner: |
| | 160 | {{{ |
| | 161 | SELECT diffSkyfile.skycell_id, warpSkyfile.skycell_id, sources, diffSkyfile.good_frac difffrac, warpSkyfile.good_frac warpfrac, diffSkyfile.good_frac - warpSkyfile.good_frac / warpSkyfile.good_frac as good_stuff FROM diffRun JOIN diffSkyfile USING (diff_id) JOIN diffInputSkyfile USING (diff_id) JOIN warpRun on (diffInputSkyfile.warp1 = warpRun.warp_id) JOIN fakeRun USING (fake_id) JOIN camRun USING (cam_id) JOIN chipRun USING (chip_id) JOIN rawExp USING (exp_id) JOIN warpSkyfile USING (warp_id) WHERE exp_name="o8507h0371o" and diff_mode=1 AND warpSkyfile.skycell_id = diffSkyfile.skycell_id AND warpSkyfile.good_frac > 0.5 GROUP BY diffSkyfile.skycell_id ORDER BY good_stuff DESC LIMIT 20; |
| | 162 | }}} |