ComCam.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "%matplotlib ipympl\n",
    "\n",
    "import lsst.daf.persistence as dafPersist\n",
    "import lsst.afw.cameraGeom.utils as cgUtils"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "import lsst.afw.display as afwDisplay\n",
    "\n",
    "afwDisplay.setDefaultBackend(\"matplotlib\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "butler = dafPersist.Butler(\"/datasets/comCam/repo\")\n",
    "camera = butler.get(\"camera\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "if False:\n",
    "    cgUtils.plotFocalPlane(camera)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "['DARK', 'FE55', 'FLAT', 'TRAP', 'LAMBDA', 'SFLAT_500']"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "butler.queryMetadata(\"raw\", [\"testType\"], run='4625')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "ename": "NoResults",
     "evalue": "No locations for get: datasetType:raw_visitInfo dataId:DataId(initialdata={'visit': 272461572, 'ccd': 'S11'}, tag=set())",
     "output_type": "error",
     "traceback": [
      "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[0;31mNoResults\u001b[0m                                 Traceback (most recent call last)",
      "\u001b[0;32m<ipython-input-8-e25635731aa5>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m\u001b[0m\n\u001b[1;32m      6\u001b[0m     \u001b[0;32mif\u001b[0m \u001b[0mdataId\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;34m\"ccd\"\u001b[0m\u001b[0;34m]\u001b[0m \u001b[0;34m!=\u001b[0m \u001b[0;34m'S11'\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m      7\u001b[0m         \u001b[0;32mcontinue\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m----> 8\u001b[0;31m     \u001b[0mexpTime\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mbutler\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mget\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m'raw_visitInfo'\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mdataId\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mgetExposureTime\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m      9\u001b[0m     \u001b[0;32mif\u001b[0m \u001b[0mexpTime\u001b[0m \u001b[0;34m>\u001b[0m \u001b[0;36m10\u001b[0m \u001b[0;32mand\u001b[0m \u001b[0mexpTime\u001b[0m \u001b[0;34m<\u001b[0m \u001b[0;36m50\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     10\u001b[0m         \u001b[0;32mif\u001b[0m \u001b[0mbutler\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mget\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m'raw_md'\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mdataId\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mget\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m\"MONOCH-WAVELENG\"\u001b[0m\u001b[0;34m)\u001b[0m \u001b[0;34m>\u001b[0m \u001b[0;36m950\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m/opt/lsst/software/stack/stack/miniconda3-4.5.4-fcd27eb/Linux64/daf_persistence/16.0-17-g7c01f5c+3/python/lsst/daf/persistence/butler.py\u001b[0m in \u001b[0;36mget\u001b[0;34m(self, datasetType, dataId, immediate, **rest)\u001b[0m\n\u001b[1;32m   1378\u001b[0m         \u001b[0mlocation\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_locate\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mdatasetType\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mdataId\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mwrite\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;32mFalse\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1379\u001b[0m         \u001b[0;32mif\u001b[0m \u001b[0mlocation\u001b[0m \u001b[0;32mis\u001b[0m \u001b[0;32mNone\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m-> 1380\u001b[0;31m             \u001b[0;32mraise\u001b[0m \u001b[0mNoResults\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m\"No locations for get:\"\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mdatasetType\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mdataId\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m   1381\u001b[0m         \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mlog\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mdebug\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m\"Get type=%s keys=%s from %s\"\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mdatasetType\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mdataId\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mstr\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mlocation\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1382\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;31mNoResults\u001b[0m: No locations for get: datasetType:raw_visitInfo dataId:DataId(initialdata={'visit': 272461572, 'ccd': 'S11'}, tag=set())"
     ]
    }
   ],
   "source": [
    "keys = (\"visit\", \"ccd\",)\n",
    "#dataIds = [did for did in [dict(zip(keys, _)) for _ in butler.queryMetadata(\"raw\", keys, run='4625', testType='LAMBDA')] if butler.datasetExists('raw', did)]\n",
    "dataIds = [dict(zip(keys, _)) for _ in butler.queryMetadata(\"raw\", keys, run='4625', testType='LAMBDA')]\n",
    "\n",
    "for dataId in dataIds:\n",
    "    if dataId[\"ccd\"] != 'S11':\n",
    "        continue\n",
    "    expTime = butler.get('raw_visitInfo', dataId).getExposureTime()\n",
    "    if expTime > 10 and expTime < 50:\n",
    "        if butler.get('raw_md', dataId).get(\"MONOCH-WAVELENG\") > 950:\n",
    "            break"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "{'visit': 272461572, 'ccd': 'S00'}\n",
      "['/datasets/comCam/repo/raw/4625/S00/S00-272461572.fits[0]']\n"
     ]
    }
   ],
   "source": [
    "for did in [dict(zip(keys, _)) for _ in butler.queryMetadata(\"raw\", keys, run='4625', testType='LAMBDA')]:\n",
    "    print(did)\n",
    "    break\n",
    "    \n",
    "print(butler.get('raw_filename', did))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "application/vnd.jupyter.widget-view+json": {
       "model_id": "787fdf85107245e6aced2f6f73a25daa",
       "version_major": 2,
       "version_minor": 0
      },
      "text/plain": [
       "FigureCanvasNbAgg()"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "ename": "NoResults",
     "evalue": "No locations for get: datasetType:raw_amp dataId:DataId(initialdata={'visit': 272461572, 'ccd': 'S22', 'run': '4625', 'channel': 1}, tag=set())",
     "output_type": "error",
     "traceback": [
      "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[0;31mNoResults\u001b[0m                                 Traceback (most recent call last)",
      "\u001b[0;32m<ipython-input-8-1a340e6cd9a8>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m\u001b[0m\n\u001b[1;32m      1\u001b[0m \u001b[0mdisp2\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mafwDisplay\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mDisplay\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;36m2\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mreopenPlot\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;32mTrue\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m      2\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m----> 3\u001b[0;31m \u001b[0mraw\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mbutler\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mget\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m'raw'\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mdataId\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mccd\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;34m'S22'\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m      4\u001b[0m \u001b[0mdisp2\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mmtv\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mraw\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m~/LSST/daf/persistence/python/lsst/daf/persistence/butler.py\u001b[0m in \u001b[0;36mget\u001b[0;34m(self, datasetType, dataId, immediate, **rest)\u001b[0m\n\u001b[1;32m   1394\u001b[0m                 \u001b[0;32mreturn\u001b[0m \u001b[0mlocation\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mmapper\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mstandardize\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mlocation\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mdatasetType\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0minnerCallback\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mdataId\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1395\u001b[0m         \u001b[0;32mif\u001b[0m \u001b[0mimmediate\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m-> 1396\u001b[0;31m             \u001b[0;32mreturn\u001b[0m \u001b[0mcallback\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m   1397\u001b[0m         \u001b[0;32mreturn\u001b[0m \u001b[0mReadProxy\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mcallback\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1398\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m~/LSST/daf/persistence/python/lsst/daf/persistence/butler.py\u001b[0m in \u001b[0;36mcallback\u001b[0;34m()\u001b[0m\n\u001b[1;32m   1392\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1393\u001b[0m             \u001b[0;32mdef\u001b[0m \u001b[0mcallback\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m-> 1394\u001b[0;31m                 \u001b[0;32mreturn\u001b[0m \u001b[0mlocation\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mmapper\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mstandardize\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mlocation\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mdatasetType\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0minnerCallback\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mdataId\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m   1395\u001b[0m         \u001b[0;32mif\u001b[0m \u001b[0mimmediate\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1396\u001b[0m             \u001b[0;32mreturn\u001b[0m \u001b[0mcallback\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m~/LSST/daf/persistence/python/lsst/daf/persistence/butler.py\u001b[0m in \u001b[0;36mcallback\u001b[0;34m()\u001b[0m\n\u001b[1;32m   1387\u001b[0m         \u001b[0;32melse\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1388\u001b[0m             \u001b[0;32mdef\u001b[0m \u001b[0mcallback\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m-> 1389\u001b[0;31m                 \u001b[0;32mreturn\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_read\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mlocation\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m   1390\u001b[0m         \u001b[0;32mif\u001b[0m \u001b[0mlocation\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mmapper\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mcanStandardize\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mlocation\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mdatasetType\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1391\u001b[0m             \u001b[0minnerCallback\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mcallback\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m~/LSST/daf/persistence/python/lsst/daf/persistence/butler.py\u001b[0m in \u001b[0;36m_read\u001b[0;34m(self, location)\u001b[0m\n\u001b[1;32m   1576\u001b[0m                 \u001b[0;32mif\u001b[0m \u001b[0mcomponentInfo\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0msubset\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1577\u001b[0m                     \u001b[0msubset\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0msubset\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mdatasetType\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mcomponentInfo\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mdatasetType\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mdataId\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mlocation\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mdataId\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m-> 1578\u001b[0;31m                     \u001b[0mcomponentInfo\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mobj\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0;34m[\u001b[0m\u001b[0mobj\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mget\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m \u001b[0;32mfor\u001b[0m \u001b[0mobj\u001b[0m \u001b[0;32min\u001b[0m \u001b[0msubset\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m   1579\u001b[0m                 \u001b[0;32melse\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1580\u001b[0m                     \u001b[0mobj\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mget\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mcomponentInfo\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mdatasetType\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mlocation\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mdataId\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mimmediate\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;32mTrue\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m~/LSST/daf/persistence/python/lsst/daf/persistence/butler.py\u001b[0m in \u001b[0;36m<listcomp>\u001b[0;34m(.0)\u001b[0m\n\u001b[1;32m   1576\u001b[0m                 \u001b[0;32mif\u001b[0m \u001b[0mcomponentInfo\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0msubset\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1577\u001b[0m                     \u001b[0msubset\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0msubset\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mdatasetType\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mcomponentInfo\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mdatasetType\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mdataId\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mlocation\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mdataId\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m-> 1578\u001b[0;31m                     \u001b[0mcomponentInfo\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mobj\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0;34m[\u001b[0m\u001b[0mobj\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mget\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m \u001b[0;32mfor\u001b[0m \u001b[0mobj\u001b[0m \u001b[0;32min\u001b[0m \u001b[0msubset\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m   1579\u001b[0m                 \u001b[0;32melse\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1580\u001b[0m                     \u001b[0mobj\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mget\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mcomponentInfo\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mdatasetType\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mlocation\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mdataId\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mimmediate\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;32mTrue\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m~/LSST/daf/persistence/python/lsst/daf/persistence/butlerSubset.py\u001b[0m in \u001b[0;36mget\u001b[0;34m(self, datasetType, **rest)\u001b[0m\n\u001b[1;32m    196\u001b[0m         \u001b[0;32mif\u001b[0m \u001b[0mdatasetType\u001b[0m \u001b[0;32mis\u001b[0m \u001b[0;32mNone\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    197\u001b[0m             \u001b[0mdatasetType\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mbutlerSubset\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mdatasetType\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 198\u001b[0;31m         \u001b[0;32mreturn\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mbutlerSubset\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mbutler\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mget\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mdatasetType\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mdataId\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m**\u001b[0m\u001b[0mrest\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    199\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    200\u001b[0m     \u001b[0;32mdef\u001b[0m \u001b[0mput\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mobj\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mdatasetType\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;32mNone\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mdoBackup\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;32mFalse\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m**\u001b[0m\u001b[0mrest\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m~/LSST/daf/persistence/python/lsst/daf/persistence/butler.py\u001b[0m in \u001b[0;36mget\u001b[0;34m(self, datasetType, dataId, immediate, **rest)\u001b[0m\n\u001b[1;32m   1378\u001b[0m         \u001b[0mlocation\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_locate\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mdatasetType\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mdataId\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mwrite\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;32mFalse\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1379\u001b[0m         \u001b[0;32mif\u001b[0m \u001b[0mlocation\u001b[0m \u001b[0;32mis\u001b[0m \u001b[0;32mNone\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m-> 1380\u001b[0;31m             \u001b[0;32mraise\u001b[0m \u001b[0mNoResults\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m\"No locations for get:\"\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mdatasetType\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mdataId\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m   1381\u001b[0m         \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mlog\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mdebug\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m\"Get type=%s keys=%s from %s\"\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mdatasetType\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mdataId\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mstr\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mlocation\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1382\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;31mNoResults\u001b[0m: No locations for get: datasetType:raw_amp dataId:DataId(initialdata={'visit': 272461572, 'ccd': 'S22', 'run': '4625', 'channel': 1}, tag=set())"
     ]
    }
   ],
   "source": [
    "disp2 = afwDisplay.Display(2, reopenPlot=True)\n",
    "\n",
    "raw = butler.get('raw', dataId, ccd='S22')\n",
    "disp2.mtv(raw)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "collapsed": false,
    "scrolled": true
   },
   "outputs": [],
   "source": [
    "from lsst.ip.isr import AssembleCcdTask\n",
    "\n",
    "config = AssembleCcdTask.ConfigClass()\n",
    "config.doTrim = False\n",
    "\n",
    "assembleTask = AssembleCcdTask(config=config)\n",
    "\n",
    "exposure = assembleTask.assembleCcd(raw)\n",
    "if not True:\n",
    "    disp2.mtv(exposure)\n",
    "else:\n",
    "    disp2.erase()\n",
    "cgUtils.overlayCcdBoxes(exposure.getDetector(), display=disp2, isTrimmed=config.doTrim) # , ignoreBBoxes=[\"ccd\", \"raw\"])\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {
    "collapsed": false,
    "scrolled": false
   },
   "outputs": [],
   "source": [
    "imageSource = cgUtils.ButlerImage(butler, type='raw', verbose=True, \n",
    "                                  callback=lambda im, ccd, imageSource : \n",
    "                                           cgUtils.rawCallback(im, ccd, imageSource, subtractBias=True),\n",
    "                                  isTrimmed=True, dataId=dataId)\n",
    "\n",
    "cgUtils.showCamera(camera, imageSource=imageSource, display=disp1, showWcs=False, binSize=1); None"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "if True:\n",
    "    mos = cgUtils.showCamera(camera, overlay=True, display=disp1)\n",
    "else:\n",
    "    mos = cgUtils.showCcd(camera[1], overlay=True, display=disp1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAmoAAAJQCAYAAAA+M0i0AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3X20ZXV95/nPV0Cgg2lQ8KmgLLUhEQQZLG0fOmlFI7ZP\nZGZNjOl0QhKXEKWJOk4SCZnuSfe4ltMmqGslGBCNrI49hqgJxNEoYjCdRYs8qJQUooDIw0J5GIwB\npRT4zh/3VLxFqm7dgjp3/6739Vqr1j1nn33P+bLdFm/2Ofvs6u4AADCeR0w9AAAA2yfUAAAGJdQA\nAAYl1AAABiXUAAAGJdQAAAYl1AAABiXUAAAGJdQAAAa159QD7A4HHnhgb9iwYeoxAAB26vLLL7+j\nuw9azro/EqG2YcOGXHbZZVOPAQCwU1X1jeWu661PAIBBCTUAgEEJNQCAQQk1AIBBCTUAgEEJNQCA\nQQk1AIBBCTUAgEEJNQCAQQk1AIBBCTUAgEEJNQCAQQk1AIBBCTUAgEEJNQCAQQk1AIBBCTUAgEEJ\nNQCAQQk1AIBBTRpqVbV/VX24qr5SVVdX1XOr6tFVdUFVfW3284ApZwQAmMrUR9TeneSvu/snkzwj\nydVJ3prkwu4+NMmFs/sAAGvOZKFWVf88yU8neV+SdPf3u/vbSY5Pcs5stXOS/Ow0EwIATGvKI2pP\nTnJ7kj+pqi9U1dlV9WNJHtfdt87W+WaSx002IQDAhPac+LWPSXJKd19SVe/Og97m7O6uqt7eL1fV\niUlOTJL169fPe9YkyUknrcjLwC4788yFn/ZRRmY/ZTXYup+OYsojajcnubm7L5nd/3AWwu1bVfWE\nJJn9vG17v9zdZ3X3xu7eeNBBB63IwAAAK2myUOvubya5qap+YrboRUk2Jzk/yQmzZSckOW+C8QAA\nJjflW59JckqSD1bVI5Ncn+RXsxCP51bVa5N8I8mrJ5wPAGAyk4Zad38xycbtPPSilZ4FAGA0U3+P\nGgAAOyDUAAAGJdQAAAYl1AAABiXUAAAGJdQAAAYl1AAABiXUAAAGJdQAAAYl1AAABiXUAAAGJdQA\nAAYl1AAABiXUAAAGJdQAAAYl1AAABiXUAAAGJdQAAAYl1AAABiXUAAAGJdQAAAYl1AAABiXUAAAG\nJdQAAAYl1AAABiXUAAAGJdQAAAYl1AAABiXUAAAGJdQAAAYl1AAABiXUAAAGJdQAAAYl1AAABiXU\nAAAGJdQAAAYl1AAABiXUAAAGJdQAAAYl1AAABiXUAAAGJdQAAAYl1AAABiXUAAAGJdQAAAYl1AAA\nBiXUAAAGJdQAAAYl1AAABiXUAAAGJdQAAAYl1AAABiXUAAAGJdQAAAYl1AAABiXUAAAGJdQAAAYl\n1AAABiXUAAAGJdQAAAYl1AAABiXUAAAGJdQAAAYl1AAABiXUAAAGJdQAAAYl1AAABiXUAAAGJdQA\nAAYl1AAABiXUAAAGtefUA/Cj67rrzs1dd23OXnvtl333PSj33HNL7r33jjz3uafnyitPT9UjklSO\nPPKNU4/KGrbUfnrddedm8+Yz8spXXjT1mKxxS+2nV111RrZsuTNJ5ZhjfnfqUdnNhBpzc889N2Wf\nfQ7MunXH5oADDk+SXHTRryVJ7r77pjzvee/MxRe/ecoRYcn99KlPfXW+9a2LpxwPkiy9nx5xxBty\n//3fz8UX/8aUIzInk4daVe2R5LIkt3T3K6rq0Un+LMmGJDckeXV33zXdhDxURx31ltx994255poP\n5DGPOTp33HF5nv70U7ZZp6ommg4WLGc/hakttZ/ed9+9ufTS03L00adOPCXzMMJn1N6Y5OpF99+a\n5MLuPjTJhbP7rEI33HB+rr/+w7n33jty9dVn5t5778itt342Dzxwf/bbb302bXpXHvWoDVOPyRq3\n1H56yy2fyZ13fiFXX/3eqcdkjVtqP73wwp/PHnvsnZtv/tTUYzIHkx5Rq6qDk7w8yduS/G+zxccn\necHs9jlJLkry2ys9Gw/fhg2v2uFjRx3lLU/GsNR+um7dsVm37tgVnAa2b6n99LjjzlvBSVhpUx9R\ne1eS30rywKJlj+vuW2e3v5nkcSs+FQDAACYLtap6RZLbuvvyHa3T3Z2kd/D7J1bVZVV12e233z6v\nMQEAJjPlEbXnJ3lVVd2Q5ENJjq2qP03yrap6QpLMft62vV/u7rO6e2N3bzzooINWamYAgBUzWah1\n96ndfXB3b0jymiSf6e5/l+T8JCfMVjshiTffAYA1aerPqG3P25P8TFV9LcmLZ/cBANacyb9HLUm6\n+6IsnN2Z7r4zyYumnAcAYAQjHlEDACBCDQBgWEINAGBQQg0AYFBCDQBgUEINAGBQQg0AYFBCDQBg\nUEINAGBQQg0AYFBCDQBgUEINAGBQQg0AYFBCDQBgUEINAGBQQg0AYFBCDQBgUEINAGBQQg0AYFBC\nDQBgUEINAGBQQg0AYFBCDQBgUEINAGBQQg0AYFBCDQBgUEINAGBQQg0AYFBCDQBgUEINAGBQQg0A\nYFBCDQBgUEINAGBQQg0AYFBCDQBgUEINAGBQQg0AYFBCDQBgUEINAGBQQg0AYFBCDQBgUEINAGBQ\nQg0AYFBCDQBgUEINAGBQQg0AYFBCDQBgUEINAGBQQg0AYFBCDQBgUEINAGBQQg0AYFBCDQBgUEIN\nAGBQQg0AYFBCDQBgUEINAGBQQg0AYFBCDQBgUEINAGBQQg0AYFBCDQBgUEINAGBQQg0AYFBCDQBg\nUEINAGBQQg0AYFBCDQBgUEINAGBQQg0AYFBCDQBgUEINAGBQQg0AYFCThVpVHVJVf1NVm6vqqqp6\n42z5o6vqgqr62uznAVPNCAAwpSmPqN2X5C3dfXiS5yQ5uaoOT/LWJBd296FJLpzdBwBYcyYLte6+\ntbuvmN3+hyRXJ1mX5Pgk58xWOyfJz04zIQDAtKq7p54hVbUhyd8meXqSG7t7/9nySnLX1vsP+p0T\nk5yYJOvXr3/mN77xjRWbFwDgoaqqy7t743LWnfxkgqraL8lHkrypu7+z+LFeqMjtlmR3n9XdG7t7\n40EHHbQCkwIArKw9p3zxqtorC5H2we7+6Gzxt6rqCd19a1U9Iclt0024rZNOmnoC2L4zz1z4aR9l\nZPZTVoOt++kopjzrs5K8L8nV3X36oofOT3LC7PYJSc5b6dkAAEYw5RG15yf5pSSbquqLs2W/k+Tt\nSc6tqtcm+UaSV080HwDApCYLte7+uyS1g4dftJKzAACMaPKTCQAA2D6hBgAwKKEGADAooQYAMCih\nBgAwKKEGADAooQYAMCihBgAwKKEGADAooQYAMCihBgAwKKEGADAooQYAMCihBgAwKKEGADAooQYA\nMCihBgAwKKEGADAooQYAMCihBgAwKKEGADAooQYAMCihBgAwKKEGADAooQYAMCihBgAwKKEGADAo\noQYAMCihBgAwKKEGADAooQYAMCihBgAwqD2Xs1JVPSLJM5I8Mcn3kny5u2+b52AAAGvdkqFWVU9N\n8ttJXpzka0luT7JPksOq6rtJzkxyTnc/MO9BAQDWmp0dUfu/krwnyUnd3YsfqKrHJvm3SX4pyTnz\nGQ8AYO1aMtS6+xeWeOy2JO/a7RMBAJBk+Z9R2yPJy5NsWPw73X36fMYCAGBZoZbkr5Lcm2RTEp9H\nAwBYAcsNtYO7+6i5TgIAwDaW+z1qn6iql8x1EgAAtrHcI2qfS/IXs+9T+0GSStLd/eNzmwwAYI1b\nbqidnuS5STY9+Gs6AACYj+W+9XlTFq5GINIAAFbIco+oXZ/koqr6RJItWxf6eg4AgPlZbqh9ffbn\nkbM/AADM2bJCrbt/b96DAACwreVemWBjktOSPCnbXpnAd6sBAMzJct/6/GCS34wrEwAArJjlhtrt\n3X3+XCcBAGAbyw21/1hVZye5MNue9fnRuUwFAMCyQ+1Xk/xkkr3yw7c+O4lQAwCYk+WG2rO6+yfm\nOgkAANtY7pUJLq6qw+c6CQAA21juEbXnJPliVX09C59R23pRdl/PAQAwJ8sNtZfOdQoAAP6JJUOt\nqvbr7ru7+xs7W2f3jwYAsLbt7DNq51XVH1TVT1fVj21dWFVPqarXVtUn42gbAMBcLHlErbtfVFUv\nS3JSkudX1QFJ7ktyTZL/N8kJ3f3N+Y8JALD27PQzat398SQfX4FZAABYZLlfzwEAwAoTagAAgxJq\nAACDWjLUqurjVbVhZUYBAGCxnR1R+5Mkn6qq06pqr5UYCACABTv7eo4/r6pPJPk/klxWVf81yQOL\nHj99zvMBAKxZy7mE1PeT3JNk7ySPyqJQAwBgfnZ2CamXJjk9yflJjunu767IVAAA7PSI2mlJfq67\nr1qJYQAA+KGdfUbtp1ZqEAAAtuV71AAABiXUAAAGtexQq6onVdWLZ7f3rapHzW8sAACW8/UcqarX\nJTkxyaOTPDXJwUn+OMmL5jXY7IzTdyfZI8nZ3f32eb0W83Hddefmrrs2Z6+99st++63P5s1n5CUv\n+cvsvff+ufLK01P1iCSVI49849SjsoYttZ9ed9252bz5jLzylRdNPSZr3FL76VVXnZEtW+5MUjnm\nmN+delR2s2WFWpKTkzw7ySVJ0t1fq6rHzmuoqtojyR8l+ZkkNye5tKrO7+7N83pNdr977rkp++xz\nYNatOzYHHHB47rrrh//z3X33TXne896Ziy9+84QTwtL76VOf+up861sXTzgdLFhqPz3iiDfk/vu/\nn4sv/o0JJ2RelhtqW7r7+1WVJKmqPZP03KZaiMJru/v62et9KMnxSYTaKnLUUW/J3XffmGuu+UAe\n85ijt7vO1n0KprKc/RSmttR+et999+bSS0/L0UefOtF0zNNyP6P22ar6nST7VtXPJPnzJH81v7Gy\nLslNi+7fPFvGKnLDDefn+us/nHvvvSP/7J89Ibfd9rls3vyePPDA/dlvv/XZtOldedSjNkw9Jmvc\nUvvpLbd8Jnfe+YVcffV7px6TNW6p/fTCC38+e+yxd26++VNTj8kcLPeI2luTvDbJpiQnJfl4krPn\nNdRyVNWJWfjcXNavXz/lKOzAhg2v2ub+y1721/94+6ijvOXJGJbaT9etOzbr1h270iPBP7HUfnrc\nceet9DisoGWFWnc/kOS9sz8r4ZYkhyy6f/Bs2eKZzkpyVpJs3Lhxnm/DAgBMYllvfVbV86vqgqr6\nalVdX1Vfr6rr5zjXpUkOraonV9Ujk7wmC9cbBQBYM5b71uf7krw5yeVJ7p/fOAu6+76q+vdJPpmF\nr+d4v+uNAgBrzXJD7e+7+xNzneRBuvvjWfgsHADAmrTcUPubqnpHko8m2bJ1YXdfMZepAABYdqj9\ny9nPjYuWdRKnQwEAzMlyz/p84bwHAQBgW8s9opaqenmSI5Lss3VZd/+neQwFAMDyv57jj5P8fJJT\nklSSn0vypDnOBQCw5i33ElLP6+5fTnJXd/9ekucmOWx+YwEAsNxQ+97s53er6olJfpDkCfMZCQCA\nZPmfUftYVe2f5B1JrsjCGZ+TXusTAOBH3XLP+vzPs5sfqaqPJdmnu/9+fmMBALBkqFXV/7LEY+nu\nj+7+kQAASHZ+RO2VSzzWWbhSAQAAc7BkqHX3r67UIAAAbGvJsz6r6l9W1Zeq6u6q+h9V9bSVGgwA\nYK3b2ddz/FGS/z3JY5KcnuRdc58IAIAkOw+1R3T3Bd29pbv/PMlBKzEUAAA7P5lg/wed+bnNfWd9\nAgDMz85C7bPZ9szPxfed9QkAMEfO+gQAGNRyr/UJAMAKE2oAAIMSagAAg3rI1/pMnPUJADBPrvUJ\nADAoZ30CAAxqZ0fU/lFVvTzJEUn22bqsu//TPIYCAGCZJxNU1R8n+fkkpySpJD+X5ElznAsAYM1b\n7lmfz+vuX05yV3f/XpLnJjlsfmMBALDcUPve7Od3q+qJSX6Q5AnzGQkAgGT5n1H7WFXtn+QdSa7I\nwhmfZ89tKgAAlhdq3f2fZzc/UlUfS7JPd//9/MYCAGC5JxOcPDuilu7ekuQRVfWGuU4GALDGLfcz\naq/r7m9vvdPddyV53XxGAgAgWX6o7VFVtfVOVe2R5JHzGQkAgGT5JxP8dZI/q6ozZ/dPmi0DAGBO\nlhtqv52FOHv97P4FcdYnAMBcLfeszweq6n1J/i4LX81xTXffP9fJAADWuGWFWlW9IMk5SW7IwiWk\nDqmqE7r7b+c3GgDA2rbctz7/IMlLuvuaJKmqw5L8P0meOa/BAADWuuWe9bnX1khLku7+apK95jMS\nAADJ8o+oXVZVZyf509n9X0xy2XxGAgAgWX6ovT7JyUl+Y3b/vyc5Yy4TAQCQZCehVlXru/vG2WWj\nTp/9AQBgBezsM2p/ufVGVX1kzrMAALDIzkKtFt1+yjwHAQBgWzsLtd7BbQAA5mxnJxM8o6q+k4Uj\na/vObmd2v7v7x+c6HQDAGrZkqHX3His1CAAA21ruF94CALDChBoAwKCEGgDAoIQaAMCghBoAwKCE\nGgDAoIQaAMCghBoAwKCEGgDAoIQaAMCghBoAwKCEGgDAoIQaAMCghBoAwKCEGgDAoIQaAMCghBoA\nwKCEGgDAoIQaAMCghBoAwKCEGgDAoIQaAMCghBoAwKCEGgDAoIQaAMCghBoAwKAmCbWqekdVfaWq\nrqyqv6iq/Rc9dmpVXVtV11TVcVPMBwAwgqmOqF2Q5OndfVSSryY5NUmq6vAkr0lyRJKXJjmjqvaY\naEYAgElNEmrd/anuvm9293NJDp7dPj7Jh7p7S3d/Pcm1SZ49xYwAAFMb4TNqv5bkE7Pb65LctOix\nm2fLAADWnD3n9cRV9ekkj9/OQ6d193mzdU5Lcl+SDz6E5z8xyYlJsn79+ocxKQDAmOYWat394qUe\nr6pfSfKKJC/q7p4tviXJIYtWO3i2bHvPf1aSs5Jk48aNvb11AABWs/phI63gi1a9NMnpSf51d9++\naPkRSf5bFj6X9sQkFyY5tLvvX+r5Nm7c2JdddtkcJwYA2D2q6vLu3ricded2RG0n/jDJ3kkuqKok\n+Vx3/3p3X1VV5ybZnIW3RE/eWaStpJNOmnoC2L4zz1z4aR9lZPZTVoOt++koJgm17v4XSzz2tiRv\nW8FxAACGNMJZnwAAbIdQAwAYlFADABiUUAMAGJRQAwAYlFADABiUUAMAGJRQAwAYlFADABiUUAMA\nGJRQAwAYlFADABiUUAMAGJRQAwAYlFADABiUUAMAGJRQAwAYlFADABiUUAMAGJRQAwAYlFADABiU\nUAMAGJRQAwAYlFADABiUUAMAGJRQAwAYlFADABiUUAMAGJRQAwAYlFADABiUUAMAGJRQAwAYlFAD\nABiUUAMAGJRQAwAYlFADABiUUAMAGJRQAwAYlFADABiUUAMAGJRQAwAYlFADABiUUAMAGJRQAwAY\nlFADABiUUAMAGJRQAwAYlFADABiUUAMAGJRQAwAYlFADABiUUAMAGJRQAwAYlFADABiUUAMAGJRQ\nAwAYlFADABiUUAMAGJRQAwAYlFADABiUUAMAGJRQAwAYlFADABiUUAMAGJRQAwAYlFADABiUUAMA\nGJRQAwAYlFADABiUUAMAGJRQAwAYlFADABiUUAMAGNSkoVZVb6mqrqoDFy07taquraprquq4KecD\nAJjSnlO9cFUdkuQlSW5ctOzwJK9JckSSJyb5dFUd1t33TzMlAMB0Jgu1JO9M8ltJzlu07PgkH+ru\nLUm+XlXXJnl2kv8xwXw8TNddd27uumtz9tprv+y770G599478oMf3JNnPvM/5MorT0/VI5JUjjzy\njVOPyhq21H563XXnZvPmM/LKV1409ZiscUvtp1dddUa2bLkzSeWYY3536lHZzSZ567Oqjk9yS3d/\n6UEPrUty06L7N8+Wbe85Tqyqy6rqsttvv31Ok/Jw3HPPTdlnnwOzfv3LcscdX8xRR70lSbJly7dz\n99035cgj35R/+Icbph2SNW+p/fSpT311HvOYoyeeEJbeT4844g15xjN+O/fcc/PEUzIPczuiVlWf\nTvL47Tx0WpLfycLbng9Zd5+V5Kwk2bhxYz+c52I+jjrqLbn77htzzTUfyK23XpQkqapt1nnwfVhp\ny9lPYWpL7af33XdvLr30tBx99KkTTsi8zC3UuvvF21teVUcmeXKSL812soOTXFFVz05yS5JDFq1+\n8GwZq9ANN5yf73zn2tx77x3ZsOFnc+WVf5Ak2Xvv/bPffuuzadO78qhHbZh2SNa8pfbTW275TO68\n8wu5+ur35mlPe93Ek7KWLbWffvKTx+eAA47IzTd/yn76I2jFP6PW3ZuSPHbr/aq6IcnG7r6jqs5P\n8t+q6vQsnExwaJLPr/SM7B4bNrxqh48dddSbV3AS2LGl9tN1647NunXHruA0sH1L7afHHXfeDh9j\n9ZvyZIJ/oruvqqpzk2xOcl+Sk53xCQCsVZOHWndveND9tyV52zTTAACMw5UJAAAGJdQAAAYl1AAA\nBiXUAAAGJdQAAAYl1AAABiXUAAAGJdQAAAYl1AAABiXUAAAGJdQAAAYl1AAABiXUAAAGJdQAAAYl\n1AAABiXUAAAGJdQAAAYl1AAABiXUAAAGJdQAAAYl1AAABiXUAAAGJdQAAAYl1AAABiXUAAAGJdQA\nAAYl1AAABiXUAAAGJdQAAAYl1AAABiXUAAAGJdQAAAYl1AAABiXUAAAGJdQAAAYl1AAABiXUAAAG\nJdQAAAYl1AAABiXUAAAGJdQAAAYl1AAABiXUAAAGJdQAAAYl1AAABiXUAAAGJdQAAAYl1AAABiXU\nAAAGJdQAAAYl1AAABiXUAAAGJdQAAAYl1AAABiXUAAAGJdQAAAYl1AAABiXUAAAGJdQAAAYl1AAA\nBiXUAAAGJdQAAAYl1AAABiXUAAAGJdQAAAYl1AAABiXUAAAGJdQAAAYl1AAABiXUAAAGJdQAAAYl\n1AAABjVZqFXVKVX1laq6qqr+y6Llp1bVtVV1TVUdN9V8AABT23OKF62qFyY5PskzuntLVT12tvzw\nJK9JckSSJyb5dFUd1t33TzEnAMCUpjqi9vokb+/uLUnS3bfNlh+f5EPdvaW7v57k2iTPnmhGAIBJ\nTRVqhyX5qaq6pKo+W1XPmi1fl+SmRevdPFsGALDmzO2tz6r6dJLHb+eh02av++gkz0nyrCTnVtVT\ndvH5T0xyYpKsX7/+4Q0LADCguYVad794R49V1euTfLS7O8nnq+qBJAcmuSXJIYtWPXi2bHvPf1aS\ns5Jk48aNvbvmBgAYxVRvff5lkhcmSVUdluSRSe5Icn6S11TV3lX15CSHJvn8RDMCAExqkrM+k7w/\nyfur6stJvp/khNnRtauq6twkm5Pcl+RkZ3wCAGtVLfTR6rZx48a+7LLLph4DAGCnqury7t64nHVd\nmQAAYFBCDQBgUEINAGBQQg0AYFBCDQBgUEINAGBQQg0AYFBCDQBgUEINAGBQQg0AYFBCDQBgUEIN\nAGBQQg0AYFBCDQBgUEINAGBQQg0AYFBCDQBgUEINAGBQQg0AYFBCDQBgUEINAGBQQg0AYFBCDQBg\nUEINAGBQ1d1Tz/CwVdXtSb4x4QgHJrljwtf/UWE7Pny24e5hO+4etuPuYTvuHiNtxyd190HLWfFH\nItSmVlWXdffGqedY7WzHh8823D1sx93Ddtw9bMfdY7VuR299AgAMSqgBAAxKqO0eZ009wI8I2/Hh\nsw13D9tx97Addw/bcfdYldvRZ9QAAAbliBoAwKCE2sNQVadU1Veq6qqq+i+Llp9aVddW1TVVddyU\nM64WVfWWquqqOnDRMttxmarqHbN98cqq+ouq2n/RY7bjLqiql8621bVV9dap51kNquqQqvqbqto8\n+/vwjbPlj66qC6rqa7OfB0w962pQVXtU1Req6mOz+7bjLqqq/avqw7O/F6+uqueu1u0o1B6iqnph\nkuOTPKO7j0jy+7Plhyd5TZIjkrw0yRlVtcdkg64CVXVIkpckuXHRMttx11yQ5OndfVSSryY5NbEd\nd9Vs2/xRkn+T5PAkvzDbhiztviRv6e7Dkzwnycmz7fbWJBd296FJLpzdZ+femOTqRfdtx1337iR/\n3d0/meQZWdieq3I7CrWH7vVJ3t7dW5Kku2+bLT8+yYe6e0t3fz3JtUmePdGMq8U7k/xWksUfmLQd\nd0F3f6q775vd/VySg2e3bcdd8+wk13b39d39/SQfysI2ZAndfWt3XzG7/Q9Z+Jfiuixsu3Nmq52T\n5GenmXD1qKqDk7w8ydmLFtuOu6Cq/nmSn07yviTp7u9397ezSrejUHvoDkvyU1V1SVV9tqqeNVu+\nLslNi9a7ebaM7aiq45Pc0t1fetBDtuND92tJPjG7bTvuGtvrYaqqDUn+pySXJHlcd986e+ibSR43\n0Virybuy8B+uDyxaZjvumicnuT3Jn8zeQj67qn4sq3Q77jn1ACOrqk8nefx2HjotC9vu0Vk4zP+s\nJOdW1VNWcLxVYyfb8Xey8LYnO7HUduzu82brnJaFt6E+uJKzQZJU1X5JPpLkTd39nar6x8e6u6vK\n1wwsoapekeS27r68ql6wvXVsx2XZM8kxSU7p7kuq6t150Nucq2k7CrUldPeLd/RYVb0+yUd74ftN\nPl9VD2ThOmK3JDlk0aoHz5atWTvajlV1ZBb+y+dLs7/QD05yRVU9O7bjP7HU/pgkVfUrSV6R5EX9\nw+/dsR13je31EFXVXlmItA9290dni79VVU/o7lur6glJbtvxM5Dk+UleVVUvS7JPkh+vqj+N7bir\nbk5yc3dfMrv/4SyE2qrcjt76fOj+MskLk6SqDkvyyCxc7PX8JK+pqr2r6slJDk3y+cmmHFh3b+ru\nx3b3hu7ekIX/cx3T3d+M7bhLquqlWXi75FXd/d1FD9mOu+bSJIdW1ZOr6pFZOBHj/IlnGl4t/JfW\n+5Jc3d2nL3ro/CQnzG6fkOS8lZ5tNenuU7v74Nnfh69J8pnu/nexHXfJ7N8hN1XVT8wWvSjJ5qzS\n7eiI2kP3/iTvr6ovJ/l+khNmRzGuqqpzs7BT3Jfk5O6+f8I5V6Xuth13zR8m2TvJBbOjk5/r7l+3\nHXdNd9/NGnDAAAAEO0lEQVRXVf8+ySeT7JHk/d191cRjrQbPT/JLSTZV1Rdny34nyduz8LGQ1yb5\nRpJXTzTfamc77rpTknxw9h9c1yf51SwcnFp129GVCQAABuWtTwCAQQk1AIBBCTUAgEEJNQCAQQk1\nAIBBCTUAgEEJNWCuqur+qvrioj8bduNzf6Cq/tcdLP/67PWuqKrnLrX+bpznkNnrPnp2/4DZ/Q3b\nWXff2XWC95jDHJ+uqgN29/MCK0+oAfP2ve4+etGfG1bodX+zu4/OwqVjzlyJF+zum5K8JwtfUJrZ\nz7N28M/8a1m4DN08voD4vyZ5wxyeF1hhQg1YcVW1T1X9SVVtqqovVNXWy7HtUVW/X1Vfrqorq+qU\n2fL/UFWXzpafVYuv9r1zf5vkX2xnhu0+Z1VdVFX/d1V9vqq+WlU/tWi2d8x+58qqOmkHr/fOJM+p\nqjcl+VdJfn8H6/1iZpewqaoXzI6unVdV11fV26vqF2czbKqqp87W+0BVvaeqPjdb7wVV9f6qurqq\nPrDouc9P8gu7sI2AQQk1YN72XfS251/Mlp2cpLv7yCwExTlVtU+SE5NsSHJ0dx+V5IOz9f+wu5/V\n3U9Psm8WLj6/XK9Msmk7y5d6zj27+9lJ3pTkP86WvTbJ33f3s5I8K8nrZtdP3UZ3/yDJb2Yh2N40\nu7+N2WVtnvKgI23PSPLrSZ6WhcsxHTab4ewsXA5nqwOSPDfJm7MQZO9MckSSI6vq6NkMdyXZu6oe\ns8OtAqwKQg2Yt8Vvff7Ps2X/KsmfJkl3fyUL1907LMmLk5zZ3ffNHvv/Zuu/sKouqapNSY7NQpjs\nzDtm1508MQuR9WBLPedHZz8vz0I4JslLkvzy7DkvSfKYLFzkfnv+TZJbkzx9B48fmOTbD1p2aXff\n2t1bklyX5FOz5ZsWzZAkfzW7rvCmJN/q7k3d/UCSqx603m1JnriD1wdWCRdlB4Y2O9J2RpKN3X1T\nVf2fSfZZxq/+Znd/+CE+55bZz/vzw78nK8kp3f3Jncx7dJKfSfKcJH9XVR/q7lsftNr3tvPPsGXR\n7QcW3X8g2/5dvWU762xvvX1mrwOsYo6oAVP471n4jFaq6rAk65Nck+SCJCdV1Z6zxx6dHwbNHVW1\nX5LdcdbmQ3nOTyZ5fVXttXXuqvqxxSvMPuf2niy85XljkndkO59Rm701uccsGHe72RyPT3LDPJ4f\nWDlCDZjCGUkeMXvb8c+S/MrsLb+zk9yY5Mqq+lKSf9vd307y3iRfzkIsXfpwX/whPufZSTYnuaKq\nvpyFM0kf/K7E65Lc2N0XzO6fkeRpVfWvt/N8n8rCW8Dz8Mwkn9v6FjKwetXCRx0AWElVdUySN3f3\nL83hud+d5PzuvnB3PzewshxRA5hAd1+R5G/m8YW3Sb4s0uBHgyNqAACDckQNAGBQQg0AYFBCDQBg\nUEINAGBQQg0AYFD/P31BTqqP8Dj6AAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7fe4eb286b90>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "cgUtils.plotFocalPlane(comCam)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Playground"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "S00 E2V-CCD250-220-Dev\n",
      "S01 E2V-CCD250-239-Dev\n",
      "S02 E2V-CCD250-154-Dev\n",
      "S10 E2V-CCD250-165-Dev\n",
      "S11 E2V-CCD250-130-Dev\n",
      "S12 E2V-CCD250-153-Dev\n",
      "S20 E2V-CCD250-163-Dev\n",
      "S21 E2V-CCD250-216-Dev\n",
      "S22 E2V-CCD250-252-Dev\n"
     ]
    }
   ],
   "source": [
    "import lsst.afw.image as afwImage\n",
    "for ccd in ['S%d%d' % (i, j) for i in range(3) for j in range(3)]:\n",
    "    print ccd, afwImage.readMetadata(butler.get('raw_filename', dataId, ccd=ccd)[0], 0).toOrderedDict()[\"LSST_NUM\"]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "OrderedDict([(u'SIMPLE', True),\n",
       "             (u'BITPIX', 8L),\n",
       "             (u'NAXIS', 0L),\n",
       "             (u'EXTEND', True),\n",
       "             (u'DATE', u'2017-05-06T19:50:55.862'),\n",
       "             (u'DATE-OBS', u'2017-05-06T19:50:37.613'),\n",
       "             (u'MJD', 57879.827035439666),\n",
       "             (u'MJD-OBS', 57879.826824224554),\n",
       "             (u'RUNNUM', u'4947D'),\n",
       "             (u'IMAGETAG', u'1494100237613'),\n",
       "             (u'CONTNUM', u'000018910d7f'),\n",
       "             (u'SEQFILE', u'/home/ts8prod/sequencers/seq-e2v-shorterp-v3.seq'),\n",
       "             (u'CCD_MANU', u'E2V'),\n",
       "             (u'CCD_TYPE', u'CCD250'),\n",
       "             (u'CCD_SERN', u'16013-12-04'),\n",
       "             (u'LSST_NUM', u'E2V-CCD250-252-Dev'),\n",
       "             (u'RAFTNAME', u'LCA-11021_RTM-005-Dev'),\n",
       "             (u'REBNAME', u'R00.Reb2'),\n",
       "             (u'TESTTYPE', u'LAMBDA'),\n",
       "             (u'IMGTYPE', u'FLAT'),\n",
       "             (u'SEQNUM', 2L),\n",
       "             (u'DETSIZE', u'[1:4096,1:4004]'),\n",
       "             (u'TEMP_SET', -100.0),\n",
       "             (u'CCDTEMP', -99.512757),\n",
       "             (u'MONOWL', 960.01),\n",
       "             (u'FILTER', u'550CutOn'),\n",
       "             (u'EXPTIME', 47.393),\n",
       "             (u'SHUT_DEL', 0.0),\n",
       "             (u'FILENAME',\n",
       "              u'/home/ts8prod/sequencers/seq-e2v-shorterp-v3.seq'),\n",
       "             (u'BINX', 1L),\n",
       "             (u'BINY', 1L),\n",
       "             (u'HEADVER', 1L),\n",
       "             (u'CCDGAIN', 1.0),\n",
       "             (u'CCDNOISE', 0.0),\n",
       "             (u'MONDIODE', -11.374909785220137),\n",
       "             (u'CCDBSS', 50.2),\n",
       "             (u'FILTPOS', 3L),\n",
       "             (u'MONOCH-WAVELENG', 960.01),\n",
       "             (u'MONOCH-SLIT_A', 240.0),\n",
       "             (u'MONOCH-SLIT_B', 240.0),\n",
       "             (u'MONOCH-SLIT_C', 0.0),\n",
       "             (u'MONOCH-BANDPASS', 0.0),\n",
       "             (u'MONOCH-FILT_1', 3.0),\n",
       "             (u'MONOCH-MSTEPS', 1200.0),\n",
       "             (u'MONOCH-GRATING', 1.0),\n",
       "             (u'AMP0-IDN', u'KEITHLEY-INSTRUMENTS-INC.,MODEL-6487'),\n",
       "             (u'AMP0-AZERO', u'F'),\n",
       "             (u'AMP0-COUNT', 2048L),\n",
       "             (u'AMP0-MEAS_NPLC', 1L),\n",
       "             (u'AMP2-IDN', u'KEITHLEY-INSTRUMENTS-INC.,MODEL-6487'),\n",
       "             (u'AMP2-CURRENT', 0.0),\n",
       "             (u'AMP2-ZERO_CHECK', u'off'),\n",
       "             (u'CHECKSUM', u'7jMbAhKZ7hKaAhKY'),\n",
       "             (u'DATASUM', u'0')])"
      ]
     },
     "execution_count": 17,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "afwImage.readMetadata(butler.get('raw_filename', dataId, ccd=ccd)[0], 0).toOrderedDict()"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "LSST",
   "language": "python",
   "name": "lsst"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}