modular_test_07062015.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "from glob import glob\n",
    "\n",
    "import des_utils\n",
    "import des_io\n",
    "\n",
    "import matplotlib.pyplot as plt\n",
    "%matplotlib inline\n",
    "%config InlineBackEnd.figure_format='retina'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "sim_files = glob(\"../ZD_KNDES_TEST/*DAT\")\n",
    "thelist,theheaders = des_utils.get_all_obs(sim_files)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "shallow_list,deep_list = des_utils.get_depth_lists(thelist)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "zbandinfo = des_utils.get_band_info(shallow_list,'z')\n",
    "ibandinfo = des_utils.get_band_info(shallow_list,'i')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "ztrig = des_utils.trigger_selector(zbandinfo)\n",
    "itrig = des_utils.trigger_selector(ibandinfo)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "zsellist,isellist = des_utils.common_trignite_selector(zbandinfo,ibandinfo)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "zSNR_sel = des_utils.get_SNR_selector(zbandinfo)\n",
    "iSNR_sel = des_utils.get_SNR_selector(ibandinfo)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "679\n"
     ]
    }
   ],
   "source": [
    "trig_flags_list,anytrigs = des_utils.get_trig_flags_list(ztrig,itrig,zsellist,isellist,zSNR_sel,iSNR_sel)\n",
    "print np.sum(anytrigs)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "MJDtriglist = des_utils.get_trig_MJD_list(zbandinfo,trig_flags_list,zsellist)\n",
    "detection_flags_list = des_utils.get_detection_flags_list(MJDtriglist,zbandinfo,ibandinfo)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "sim_triggers, sim_colors, sim_ifluxes,simdet,SNIDset,deltaT = des_utils.extract_colors(sim_files)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[] 0\n",
      "[] 0\n"
     ]
    }
   ],
   "source": [
    "SNIDset2 = set()\n",
    "for i in range(0,len(detection_flags_list)):\n",
    "    if detection_flags_list[i] and anytrigs[i]:\n",
    "        SNIDset2.add(theheaders[i]['SNID'])\n",
    "        if not (theheaders[i]['SNID'] in SNIDset):\n",
    "            print i\n",
    "            print theheaders[i]['SNID']\n",
    "            print \"zbandinfo: \", zbandinfo[i]\n",
    "            print \"ibandinfo: \", ibandinfo[i]\n",
    "            print \"MJDtriglist: \", MJDtriglist[i]\n",
    "            \n",
    "print list(set(SNIDset)-set(SNIDset2)), len(list(set(SNIDset)-set(SNIDset2)))\n",
    "print list(set(SNIDset2)-set(SNIDset)), len(list(set(SNIDset2)-set(SNIDset)))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "9012\n"
     ]
    }
   ],
   "source": [
    "print theheaders[900]['SNID']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "True\n"
     ]
    }
   ],
   "source": [
    "print detection_flags_list[900]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[False False False False False False False False False False False False\n",
      " False False False]\n"
     ]
    }
   ],
   "source": [
    "print zSNR_sel[900]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": [
       "iVBORw0KGgoAAAANSUhEUgAAAXMAAAEACAYAAABBDJb9AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\n",
       "AAALEgAACxIB0t1+/AAAEMxJREFUeJzt3X/oXfddx/HXK2nKlgUWayWJbUb6h9NOxMQfdWzVfe0c\n",
       "hCFZBalWZGUMKaJrqSjtitCv7o/9gM2hwhCWjq/bqIyVldR1LrHm0oo0Y5q0adMaC620mnyrtov9\n",
       "gdjZt3/c821ubu695/c5937u8wFfcu69557zzud7vq+cfO77nK8jQgCAxbap7wIAAPUR5gCQAMIc\n",
       "ABJAmANAAghzAEgAYQ4ACZgZ5rbfYvuY7RO2T9n+ZPb8ZbaP2D5t+7Dt7d2UCwCYxHl95ra3RsRr\n",
       "ti+R9PeSfl/SAUn/GRGfsX27pB+IiDvaLxcAMEnuNEtEvJYtXipps6SXNAzztez5NUnXt1IdAKCQ\n",
       "3DC3vcn2CUnrko5GxBOSdkTEerbKuqQdLdYIAMhxSd4KEfGGpL223y7p27Z/cez1sM09AQCgR7lh\n",
       "viEiztn+pqSflrRue2dEnLW9S9IL4+sT8ABQTUS4ypumfkm6XNL2bPmtkh6S9H5Jn5F0e/b8HZI+\n",
       "NeG9MWvbfX1JWu27BmqipmWsi5oK1xRV3pd3Zr5L0prtTRrOr385Ih60fVzS12x/VNKzkm4o/a8I\n",
       "AKAxM8M8Ik5K+qkJz78o6ZfaKgoAUM4yXgE66LuACQZ9FzDBoO8CJhj0XcAEg74LmGLQdwETDPou\n",
       "YIJB3wU0JfeiocobtiOqTOIDwBKrmp3LeGYOAMkhzAEgAYQ5ACSAMAeABBDmAJAAwhwAEkCYA0AC\n",
       "CHMASEDhuyYC82r8Dp1crIZlxJk5EhHZF7CcCHMASABhDgAJIMwBIAGEOQAkgDAHgAQQ5gCQAPrM\n",
       "sXDG+8oBEOZYWKN5zjVCANMsAJAAwhwAEkCYA0ACCHMASABhDgAJIMwBIAGEOQAkgDAHgATMDHPb\n",
       "u20ftf2E7cdt35I9v2r7edvHs6/93ZQLAJjEEdOvjLa9U9LOiDhhe5ukf5R0vaQbJL0cEZ+b8d7g\n",
       "13ehDcPL+cevAI03lznusMiqZufMy/kj4qyks9nyK7aflHTFxj5LVwkAaEXhOXPbeyTtk/RI9tTH\n",
       "bD9q+6Dt7S3UBgAoqFCYZ1MsX5d0a0S8IukLkq6StFfSGUmfba1CAECu3Lsm2t4i6V5JX4mI+yQp\n",
       "Il4Yef2Lku6f8t7VkYeDiBjUKRYAUmN7RdJK7e3kfABqSWuS/isibht5fldEnMmWb5P0sxHxG2Pv\n",
       "5QNQtIIPQJGyqtmZF+bXSnpI0mM6/9Nyp6QbNZxiCUnPSLo5ItabKAjIQ5gjZa2EeR2EOdpCmCNl\n",
       "VbOTK0ABIAGEOQAkgDAHgATwC52BEcP5+POYf8ei4MwcuEjowg9YgflHmANAAghzAEgAYQ4ACSDM\n",
       "ASABhDkAJIAwB4AE0GeOZNEzjmXCmTkSR884lgNhDgAJIMwBIAGEOQAkgDAHgAQQ5gCQAFoTMbdo\n",
       "LQSK48wcc47WQqAIwhwAEkCYA0ACCHMASABhDgAJIMwBIAGEOQAkgD5zzJXx3vJlQ289quLMHHNo\n",
       "2XvLl/3vjyoIcwBIwMwwt73b9lHbT9h+3PYt2fOX2T5i+7Ttw7a3d1MuAGCSvDPz1yXdFhE/Lund\n",
       "kn7H9tWS7pB0JCLeKenB7DEAoCczwzwizkbEiWz5FUlPSrpC0gFJa9lqa5Kub7NIAMBshefMbe+R\n",
       "tE/SMUk7ImI9e2ld0o7GKwMAFFaoNdH2Nkn3Sro1Il62z3dLRURMayezvTrycBARg+qlAu1Y9nZI\n",
       "9Mv2iqSV2tuJmH0c294i6a8lfSsiPp8995SklYg4a3uXpKMR8WNj7wt6ZFHWMFg3jklrdHnjeLpw\n",
       "naLreWrPdpF9dqVozUhX1ezM62axpIOSTm0EeeaQpJuy5Zsk3Vd2xwCA5sw8M7d9raSHJD2m86cL\n",
       "H5f0HUlfk/QOSc9KuiEivjf2Xs7MURpn5pyZL7uq2Zk7zVIVYY4qCHPCfNm1Ms0CAFgMhDkAJIC7\n",
       "JmJhLEILIXc9RF84M8cCWZS7CS5KnUgJYQ4ACSDMASABhDkAJIAwB4AEEOYAkADCHAASQJ85kjPj\n",
       "lsxvPt/PZfrn0X+OpnFmjgRN6/Puu/+77/0jZYQ5ACSAMAeABBDmAJAAwhwAEkCYA0ACaE0ECmj6\n",
       "9rtF2yT7bKfEYiHMgcIm/2q5i3+FXdltNbEelh3TLACQAMIcABJAmANAAghzAEgAYQ4ACSDMASAB\n",
       "tCZiKTXdNw70jTDHkqrSGw7ML6ZZACABuWFu+27b67ZPjjy3avt528ezr/3tlgkAmKXImfmXJI2H\n",
       "dUj6XETsy77+pvnSAABF5YZ5RDws6aUJLzHRCABzos6c+cdsP2r7oO3tjVUEACitajfLFyT9cbb8\n",
       "CUmflfTR8ZVsr448HETEoOL+gF7UbWGkBRJ5bK9IWqm9nYj8Y832Hkn3R8RPFH3NdnD/ZZQ1DL+8\n",
       "W81Ouu1smfc0+f762xr9OZn19+fnaTlUzc5K0yy2d408/BVJJ6etCwBoX+40i+17JL1P0uW2n5N0\n",
       "l6QV23s1PG14RtLNrVYJAJip0DRLpQ0zzYIKmGZhmmXZdTrNAgCYL4Q5ACSAMAeABHDXRPSOXmyg\n",
       "Ps7MMSdCF36QCKAMwhwAEkCYA0ACCHMASABhDgAJIMwBIAG0JqIXtCOWNzpmky73njSm3AJgeRDm\n",
       "6NHoPUiQr8h4jd8nBsuCaRYASABhDgAJIMwBIAGEOQAkgDAHgAQQ5gCQAFoT0Zrxvmd6noH2cGaO\n",
       "lnFrW6ALhDkAJIAwB4AEEOYAkADCHAASQJgDQAJoTURttCB2L+92uFg+nJmjIbQgdovxxoUIcwBI\n",
       "QG6Y277b9rrtkyPPXWb7iO3Ttg/b3t5umQCAWYqcmX9J0v6x5+6QdCQi3inpwewxAKAnuWEeEQ9L\n",
       "emns6QOS1rLlNUnXN1wXAKCEqnPmOyJiPVtel7SjoXoAABXUbk2MiJj2m9Ztr448HETEoO7+MB+m\n",
       "fc8xX2hhnH+2VySt1N1O1TBft70zIs7a3iXphUkrRcRq5cqwAIr8tnj0i+/RvMtOcgcbj23fVWU7\n",
       "VadZDkm6KVu+SdJ9FbcDAGhAkdbEeyT9g6Qftf2c7Y9I+pSkD9g+Lem67DEAoCeOaGfq03YwR5eu\n",
       "4Vzs6H/hzy9vfN/H1xk9Hma9v9xy3++vv6264zJ5vGfvB/OranZyBSgAJIAwB4AEEOYAkABugYvO\n",
       "bmFLb/pkjAuaQJgj00U/8viHcxhiXFAf0ywAkADCHAASQJgDQAIIcwBIAGEOAAmgmwVYMrNaIbnk\n",
       "f3FxZg4spdD5lsjRZSwqwhwAEkCYA0ACCHMASABhDgAJIMwBIAGEOQAkgD7zJVX2tqtl1ueWrkD3\n",
       "CPOlVva2t0Vv1drF7XQBjGKaBQASQJgDQAIIcwBIAGEOAAkgzAEgAXSzYCpaDNNS5Ps5vg63xF0c\n",
       "hDly0GaYjiLfy6Ltp5g3TLMAQAJqnZnbflbSf0v6P0mvR8Q1TRQFACin7jRLSFqJiBebKAYAUE0T\n",
       "0yxMrAFAz+qGeUj6W9vftf1bTRQEACiv7jTLeyPijO0fknTE9lMR8XAThQEAiqsV5hFxJvvzP2x/\n",
       "Q9I1kt4Mc9urI6sPImJQZ38ALjYv1wPQo16N7RVJK7W3E1HtOLC9VdLmiHjZ9tskHZb0RxFxOHs9\n",
       "+GbOr+EP3mjf8fnlje/b9HUufk+55brvp5auainzMzx+vPDzX03V7KxzZr5D0jdsb2znqxtBDgDo\n",
       "VuUwj4hnJO1tsBYAQEVcAQoACSDMASABhDkAJIC7JjZsUptY1U/1m9xW3f0CmG+EeSuavI1oH7ck\n",
       "5ba3wKJhmgUAEkCYA0ACCHMASABhDgAJIMwBIAF0s/SkjzvM0XKIsoocM3WO3Wk/B3215S4ywrxX\n",
       "fbQA0naIMorcabHJfUx6vqn9pI1pFgBIAGEOAAkgzAEgAYQ5ACSAMAeABBDmAJAAWhM7MNozm9cr\n",
       "W6UXvEyvLrDoZh3XTfWi5/3stNFbXxdh3omyPblVenjp1cUyaavnPW8fTe2n+ZqZZgGABBDmAJAA\n",
       "whwAEkCYA0ACev8A1PZmSVeNPf2vEfF6H/UAwCLqPcwlvV3SaWnnq8OH61ul2CPpuaIbaLKNaNat\n",
       "N9tsiSraRlikzZGWRHRp2vHWRJtt1e2VWbeLW/h2cfveeQhzSVv/Vzqzbbj8g69KL1bYxnjrUJ3W\n",
       "n1ltSF3fBrTKetzmFl0q2hZb5Gdn1rFbtcW3zP6rqPJ3aR5z5gCQgMphbnu/7ads/4vt25ssCgBQ\n",
       "TqUwzz60/HNJ+yW9S9KNtq9usrC22F7pu4bFMOi7gAkGfRewQAZ9FzDBoO8Cklb1zPwaSU9HxLNZ\n",
       "18lfSfpQc2W1aqXvAhbDoO8CJhj0XcACGfRdwASDvgtIWtUwv0IXdps8nz0HAOhB1W6Whlvf/meL\n",
       "dN254fK5tzW7bQBInyPK57Ltd0tajYj92eOPS3ojIj49sg69zgBQQZW+9Kphfomkf5b0fkn/Luk7\n",
       "km6MiCdLbwwAUFulaZaI+L7t35X0bUmbJR0kyAGgP5XOzAEA86X2FaBFLh6y/afZ64/a3ld3n3Vr\n",
       "sr1i+5zt49nXH3ZQ0922122fnLFO1+M0s6aexmm37aO2n7D9uO1bpqzX2VgVqanrsbL9FtvHbJ+w\n",
       "fcr2J6es1/UxlVtXH8dVtt/N2f7un/J6p2OVV1PpcYqIyl8aTrE8LWmPpC2STki6emydD0p6IFv+\n",
       "OUmP1NlnQzWtSDrUZh0T6vp5SfsknZzyeqfjVLCmPsZpp6S92fI2DT+b6fuYKlJTH2O1NfvzEkmP\n",
       "SLq272OqYF2dj1W239+T9NVJ++5xrGbVVGqc6p6ZF7l46ICkNUmKiGOSttveUXO/dWuSOr4LVUQ8\n",
       "LOmlGat0PU5FapK6H6ezEXEiW35F0pOSfnhstU7HqmBNUvdj9Vq2eKmGJzHjd6jr/JgqWJfU8VjZ\n",
       "vlLDwP7ilH13PlYFatKM5y9SN8yLXDw0aZ0ra+63bk0h6T3Zf6cesP2uFuspqutxKqLXcbK9R8P/\n",
       "ORwbe6m3sZpRU+djZXuT7ROS1iUdjYhTY6v0Mk4F6urjuPoTSX8g6Y0pr/cxVnk1lRqnumFe9NPT\n",
       "WffGbFqRbf+TpN0R8ZOS/kzSfS3WU0aX41REb+Nke5ukr0u6NTsbvmiVscetj1VOTZ2PVUS8ERF7\n",
       "NQydX5hy36HOx6lAXZ2Ole1flvRCRBzX7DPdzsaqYE2lxqlumP+bpN0jj3dr+C/arHWuzJ5rS25N\n",
       "EfHyxn8FI+JbkrbYvqzFmoroepxy9TVOtrdIulfSVyJi0gHc+Vjl1dTnMRUR5yR9U9LPjL3U6zE1\n",
       "ra4exuo9kg7YfkbSPZKus/2XY+t0PVa5NZUdp7ph/l1JP2J7j+1LJf2apENj6xyS9GHpzStHvxcR\n",
       "6zX3W6sm2ztsb/xGkGs0bNGs8hsxmtT1OOXqY5yy/R2UdCoiPj9ltU7HqkhNXY+V7cttb8+W3yrp\n",
       "A5KOj63W+TFVpK6uxyoi7oyI3RFxlaRfl/R3EfHhsdU6HasiNZUdp1q/aSimXDxk++bs9b+IiAds\n",
       "f9D205JelfSROvtsoiZJvyrpt21/X9JrGg5mq2zfI+l9ki63/ZykuzTstullnIrUpB7GSdJ7Jf2m\n",
       "pMdsb4TAnZLesVFXD2OVW5O6H6tdktZsb9LwpOzLEfFgnz97RetSP8fVqJCkORirmTWp5Dhx0RAA\n",
       "JIBfGwcACSDMASABhDkAJIAwB4AEEOYAkADCHAASQJgDQAIIcwBIwP8DahCdRyzkih8AAAAASUVO\n",
       "RK5CYII=\n"
      ],
      "text/plain": [
       "<matplotlib.figure.Figure at 0x106bea890>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "colors = []\n",
    "for i in range(0,len(trig_flags_list)):\n",
    "    if anytrigs[i]:\n",
    "        imag = -2.5*np.log(ibandinfo[i][2][itrig[i] & isellist[i]][0])\n",
    "        zmag = -2.5*np.log(zbandinfo[i][2][ztrig[i] & zsellist[i]][0])\n",
    "        colors.append(imag-zmag)\n",
    "plt.hist(colors,bins=100);"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "1762\n"
     ]
    },
    {
     "data": {
      "image/png": [
       "iVBORw0KGgoAAAANSUhEUgAAAXQAAAEACAYAAACj0I2EAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\n",
       "AAALEgAACxIB0t1+/AAAEA9JREFUeJzt3X2sZHddx/H3h25reVA2Dcm2lJJWw4YH5aEx0CgJF6Sk\n",
       "IdqWGIskmIrgP2JAgujWJrqJEShGwGj4Q3nIqmm1gmmoQOjN2htITIpIWwrbdcGk0kL2Vp5qgajF\n",
       "fv1jztLZy717Z87MvWfmd9+v5KbnzJwz8+29M5/9zu+c85tUFZKk5fe4oQuQJM2HgS5JjTDQJakR\n",
       "BrokNcJAl6RGGOiS1IgzBnqSDyZZT3LP2G3nJVlNciLJbUn2j913XZIvJTme5BU7Wbgk6XTbdegf\n",
       "Aq7YcNshYLWqDgJHu3WSPBt4NfDsbp/3JfETgCTtkjMGblV9GvjWhpuvBI50y0eAq7vlq4CbquqR\n",
       "qroP+DLwwvmVKkk6kz4d9IGqWu+W14ED3fJTgQfGtnsAuHCG2iRJU5hpSKRG8wacae4A5xWQpF2y\n",
       "r8c+60nOr6qTSS4AHuxu/ypw0dh2T+tuO00SQ16SeqiqnOn+PoH+UeBa4Ibuv7eM3X5jknczGmp5\n",
       "BvCZPkUtsiSHq+rw0HX0Zf3Dsv7hLHPtMFkzfMZAT3IT8BLgKUnuB34feCdwc5LXA/cB1wBU1bEk\n",
       "NwPHgO8Dv1FO5ShJu+aMgV5Vr9nirpdvsf3bgbfPWpQkaXqeJz69taELmNHa0AXMaG3oAma0NnQB\n",
       "M1obuoAZrA1dwE7Lbo+KJKllHkOXpCFMkp19DopKmtJmB7RsbDRvBrq0a8Yz3SzX/DmGLkmNMNAl\n",
       "qREGuiQ1wkCXpEYY6JLUCM9ykRbIxtMbPbVR07BDlxbOdrNSS5sz0CWpEQa6JDXCQJekRhjoktQI\n",
       "A12SGmGgS1IjDHRJaoQXFkk7ZJIv9ZXmyQ5d2lFeJKTdY6BLUiMMdElqhIEuSY0w0CWpEQa6JDXC\n",
       "QJekRhjoktQIA12SGmGgS1IjDHRJaoSBLkmNMNAlqREGuiQ1wkCXpEYY6JLUCANdkhrRO9CTXJfk\n",
       "i0nuSXJjkh9Jcl6S1SQnktyWZP88i5Ukba1XoCe5GPh14NKq+ingLOCXgUPAalUdBI5265KkXdC3\n",
       "Q/8v4BHgCUn2AU8AvgZcCRzptjkCXD1zhdICSlLjP0PXI0HPQK+qbwJ/AnyFUZB/u6pWgQNVtd5t\n",
       "tg4cmEuV0kLy+0K1WPb12SnJTwC/BVwMPAT8fZLXjm9TVVt2LkkOj62uVdVanzokqVVJVoCVqfap\n",
       "mr7DSPJq4PKqekO3/ivAZcDLgJdW1ckkFwC3V9UzN+xbVZWpn1RaIKNm5dR7J2z2mt64zend/GT7\n",
       "+F7RKZNkZ98x9OPAZUkenyTAy4FjwK3Atd021wK39Hx8SdKUeg25VNXdSf4K+CzwKPA54C+AHwVu\n",
       "TvJ64D7gmjnVKUnaRq8hl5me0CEXNcAhF+22nRxykSQtGANdkhphoEtSIwx0SWqEgS5JjTDQJakR\n",
       "vc5Dl3Q6J+jSIrBDl+ZifKIuJ+3SMAx0SWqEgS5JjTDQJakRBrokNcJAl6RGeNqiNLAznfI4fp8z\n",
       "L2o7dujSQtjqVEdPgdTkDHRJaoSBLkmNMNAlqREeFJUGMu38Lxu39yCpNjLQpcGMf9/oNNtPs4/2\n",
       "EgNdGrNVF+xsiloGBrr0Q7bqnO2Qtdg8KCpJjTDQJakRBrokNcJAl6RGGOiS1AgDXZIaYaBLUiMM\n",
       "dElqhIEuSY0w0CWpEQa6JDXCQJekRhjoktSI3oGeZH+SDye5N8mxJC9Kcl6S1SQnktyWZP88i5X0\n",
       "w5LU+M/Q9Wg4s3Tofwp8vKqeBTwXOA4cAlar6iBwtFuXtOOK06f31V6UqulfBEmeDNxZVT++4fbj\n",
       "wEuqaj3J+cBaVT1zwzblV2dpUY063MfmQz/9Cy42zodePZfns//mtcWvpmvUJNnZt0O/BPjPJB9K\n",
       "8rkkf5nkicCBqlrvtlkHDvR8fGkhOIyhZdI30PcBlwLvq6pLge+yYXilRq2/bwQtOV/GWh59v4Lu\n",
       "AeCBqvqXbv3DwHXAySTnV9XJJBcAD262c5LDY6trVbXWsw5JalKSFWBlqn36jKF3T/Yp4A1VdaIL\n",
       "6Cd0d32jqm5IcgjYX1WHNuznGLoW1sbx6J0cA3cMXdOYJDtnCfTnAe8HzgH+HXgdcBZwM/B04D7g\n",
       "mqr69rRFSUMx0LWodjTQ+zLQtcgMdC2qSbKz7xi6pIF59o028tJ/aWl5Bo5OZ6BLUiMMdElqhIEu\n",
       "SY0w0CWpEQa6JDXCQJekRhjoktQIA12SGmGgS1IjDHRJaoRzuUiN2jjXi5N2tc8OXWqa873sJQa6\n",
       "JDXCQJekRhjoktQIA12SGmGgS1IjDHRJaoSBLkmNMNAlqREGuiQ1wkCXpEYY6JLUCCfn0p63cRIr\n",
       "aVnZoUuAk1ipBXboUmP8xLF32aFLzfHTxl5loEtSIwx0SWqEgS5JjTDQJakRBrokNcJAl6RGGOiS\n",
       "1IiZAj3JWUnuTHJrt35ektUkJ5LclmT/fMqUJG1n1g79zcAxHruK4RCwWlUHgaPduiRpF/QO9CRP\n",
       "A14JvB9Id/OVwJFu+Qhw9UzVSZImNkuH/h7gbcCjY7cdqKr1bnkdODDD40uSptAr0JP8PPBgVd3J\n",
       "Y935aarKCSUkaRf1nW3xZ4Ark7wSOBf4sSR/DawnOb+qTia5AHhws52THB5bXauqtZ51SFKTkqwA\n",
       "K1PtM2qkZ3rSlwC/XVW/kORdwDeq6oYkh4D9VXVow/ZVVZt29dIQRtPNnnofhO2XJ91usfb3fbfc\n",
       "JsnOeZ2HfupV807g8iQngJd165KkXTBzhz71E9qha8HYoWsZTJKdfmORtEds/CYjA749Brq0Z2zs\n",
       "8NUa53KRpEYY6JLUCANdkhphoEtSIwx0SWqEZ7loT9p4Cp/UAjt07WHOH6e2GOiS1AgDXZIa4Ri6\n",
       "9gzHzU83/vvYbBoApwpYPnbo2mMcN3/MJL8Lf1/LxECXpEYY6JLUCANdkhrhQVE1y4Og0/N3ttwM\n",
       "dDXOOcCnN/6NR1omDrlIUiMMdElqhEMukhw7b4QduiS8gKgNBrokNcJAl6RGGOiS1AgPiqopHtzT\n",
       "XmaHrgZ5gE97k4EuSY0w0CWpEQa6JDXCQJekRhjoktQIA12SGmGgS1IjDHRJaoSBLkmN6BXoSS5K\n",
       "cnuSLyb5QpI3dbefl2Q1yYkktyXZP99yJUlb6duhPwK8paqeA1wGvDHJs4BDwGpVHQSOduuSpF3Q\n",
       "K9Cr6mRV3dUtfwe4F7gQuBI40m12BLh6HkVKkrY38xh6kouBFwB3AAeqar27ax04MOvjS5ImM9P0\n",
       "uUmeBHwEeHNVPZzkB/dVVW01lWmSw2Ora1W1NksdktSaJCvAylT7VPWbZjTJ2cA/Ap+oqvd2tx0H\n",
       "VqrqZJILgNur6pkb9quqyg8/ojS7URNx6jUdTp9Gd3x9luW9ub/v22FNkp19z3IJ8AHg2Kkw73wU\n",
       "uLZbvha4pc/jS5Km16tDT/Ji4FPA53nsn/DrgM8ANwNPB+4Drqmqb2/Y1w5dO8YO3Q69VZNkZ+8h\n",
       "l74MdO0kA91Ab9Uk2el3imrp+T2i0oiX/qsRfo+oZIeupbGxE3cIYHf5+198BrqWzPj4rnbXxjF4\n",
       "LRqHXCSpEQa6JDXCQJekRhjoktQIA12SGmGgS1IjDHRJaoSBLkmNMNAlqRFeKaql5aRc0uns0LXE\n",
       "nJBLGmegS1IjDHRJaoRj6FooTtG6fPr8zfw77ww7dC0gx8aXT5+/mX/neTPQJakRDrlI6mW700Yd\n",
       "Vtl9duiSeppkyMRhld1kh66F5sVD0uQMdC04v8dSmpSBrkE4virNn2PoGpDjq9I82aFrS5uNX9tJ\n",
       "S4vLQNc2HMOWloVDLpLUCDt0LQRPT2xHn7/l+D4O6/Vnh64F4QHSdjivy1Ds0HWaaS7nPtVJTXoK\n",
       "ol343javv78H67dmoGsTp94vm71HtrrvTPv02U7tmecBdg/Wb8ZA19w7ZztxaRgGujrz7JzHH8tO\n",
       "Stotcz8omuSKJMeTfCnJ78778ZdZkhr/2enH2riNnbOWyTSvXV/jI3MN9CRnAX8OXAE8G3hNkmfN\n",
       "8zmGlmRltkeY59H8aaYvHW03e/1DWxu6gBmtDV3AjNZ2+fkmfb9sv93yv/a3N+8O/YXAl6vqvqp6\n",
       "BPhb4KpJd07y5OTsh5Nzvzf6OfvhJBfPucZZrQxdwIxWhi5gNmtDFzCjtaELmNHa0AXMYmXoAnba\n",
       "vMfQLwTuH1t/AHjRFPsH9p0DXztntHrwu/D1H/yjM8npcbOe0jTp/jv10W4eQzG78TzSTjnTEOK0\n",
       "+89yOuNWedM3IyY9nXeWmucd6HMIif/ZB7/40Gj5oSdu/RRn+n+e9UDcpPtvdvBvswOB0x5wnPWx\n",
       "PIVQy2yRTo3drpbtnmd3a07V/Bq1JJcBh6vqim79OuDRqrphbBs7Q0nqYbvufd6Bvg/4N+DngK8B\n",
       "nwFeU1X3zu1JJEmbmuuQS1V9P8lvAp8EzgI+YJhL0u6Ya4cuSRrOoLMtJnlrkkeTnDdkHdNK8odJ\n",
       "7k5yV5KjSS4auqZpJPnjJPd2/w//kOTJQ9c0jSS/lOSLSf4vyaVD1zOJZb7gLskHk6wnuWfoWvpI\n",
       "clGS27vXzBeSvGnomqaR5Nwkd3R5cyzJO7badrBA70LwcuA/hqphBu+qqudV1fOBW4A/GLqgKd0G\n",
       "PKeqngecAK4buJ5p3QO8CvjU0IVMooEL7j7EqPZl9Qjwlqp6DnAZ8MZl+v1X1X8DL+3y5rnAS5O8\n",
       "eLNth+zQ3w38zoDP31tVPTy2+iTg60PV0kdVrVbVo93qHcDThqxnWlV1vKpODF3HFGa64G5oVfVp\n",
       "4FtD19FXVZ2sqru65e8A9wJPHbaq6VTV97rFcxgdn/zmZtsNEuhJrgIeqKrPD/H885Dkj5J8BbgW\n",
       "eOfQ9czg14CPD11E4za74O7CgWrZ07orz1/AqJFZGkkel+QuYB24vaqObbbdjs22mGQVOH+Tu65n\n",
       "9BH/FeOb71QdfZ2h/t+rqlur6nrg+iSHgPcAr9vVArexXf3dNtcD/1tVN+5qcROYpP4l4pkHCyDJ\n",
       "k4APA2/uOvWl0X2ifn53vOuTSVaqam3jdjsW6FV1+Wa3J/lJ4BLg7iQw+rj/r0leWFUP7lQ909qq\n",
       "/k3cyAJ2uNvVn+RXgVcyumZg4Uzx+18GXwXGD5xfxKhL1y5JcjbwEeBvquqWoevpq6oeSvIx4KfZ\n",
       "ZGKdXR9yqaovVNWBqrqkqi5h9MK+dJHCfDtJnjG2ehVw51C19JHkCuBtwFXdAZdltnCf7jbxWeAZ\n",
       "SS5Ocg7wauCjA9e0Z2TUOX4AOFZV7x26nmkleUqS/d3y4xmdTLJp5izCl0Qv48fRdyS5pxvTWgHe\n",
       "OnA90/ozRgdzV5PcmeR9Qxc0jSSvSnI/ozMWPpbkE0PXdCZV9X3g1AV3x4C/W6YL7pLcBPwzcDDJ\n",
       "/UkWanhxAj8LvJbR2SF3dj/LdNbOBcA/dXlzB3BrVR3dbEMvLJKkRixChy5JmgMDXZIaYaBLUiMM\n",
       "dElqhIEuSY0w0CWpEQa6JDXCQJekRvw/uNev75ze9HYAAAAASUVORK5CYII=\n"
      ],
      "text/plain": [
       "<matplotlib.figure.Figure at 0x103ea0910>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "data_files = glob(\"../DESY1_numepoch1_iz_real_text/*.dat\")\n",
    "datalist, dataheaders = des_utils.get_all_obs(data_files)\n",
    "shallow_list_data,deep_list_data = des_utils.get_depth_lists(datalist)\n",
    "zbandinfo_data = des_utils.get_band_info(shallow_list_data,'z')\n",
    "ibandinfo_data = des_utils.get_band_info(shallow_list_data,'i')\n",
    "ztrig_data = des_utils.trigger_selector(zbandinfo_data)\n",
    "itrig_data = des_utils.trigger_selector(ibandinfo_data)\n",
    "zsellist_data,isellist_data = des_utils.common_trignite_selector(zbandinfo_data,ibandinfo_data)\n",
    "trig_flags_list_data,anytrigs_data = des_utils.get_trig_flags_list(ztrig_data,itrig_data,zsellist_data,isellist_data)\n",
    "print np.sum(anytrigs_data)\n",
    "MJDtriglist_data = des_utils.get_trig_MJD_list(zbandinfo_data,trig_flags_list_data,zsellist_data)\n",
    "detection_flags_list_data = des_utils.get_detection_flags_list(MJDtriglist_data,zbandinfo_data,ibandinfo_data)\n",
    "colors_data = []\n",
    "for i in range(0,len(trig_flags_list_data)):\n",
    "    if anytrigs_data[i]:\n",
    "        imag_data = -2.5*np.log(ibandinfo_data[i][2][itrig_data[i] & isellist_data[i]][0])\n",
    "        zmag_data = -2.5*np.log(zbandinfo_data[i][2][ztrig_data[i] & zsellist_data[i]][0])\n",
    "        colors_data.append(imag_data-zmag_data)\n",
    "plt.hist(colors_data,bins=100);"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.text.Text at 0x107b80390>"
      ]
     },
     "execution_count": 17,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": [
       "iVBORw0KGgoAAAANSUhEUgAAAYMAAAEZCAYAAAB1mUk3AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\n",
       "AAALEgAACxIB0t1+/AAAHhBJREFUeJzt3X+cXXV95/HXm8BkSIZkEn4kgGECKfQRU/mlYtRYg4Bm\n",
       "LYLWgtBFQdlut+su1lUrPx5K6o8W6S9kbR91V0RAUEEqj2ItJVIHte6CKAEkYbOykgQ1E4EhMhEw\n",
       "k3z2j/M9mTOTe2fuJHPPuXfu+/l4zCP3/Lrnc09m7ud8fx5FBGZm1tn2qzoAMzOrnpOBmZk5GZiZ\n",
       "mZOBmZnhZGBmZjgZmJkZTgbWoiStlnTTFLzPjyT99lTENIlz9ku6uM62xZJ2SSr9b0/SFyR9vOzz\n",
       "WntwMrBWNSUDYCLityLi21PxXpM5LVMU/xRrOK7xEppNT04G1nSS9q86BttNDe7XisnMmsjJwJpC\n",
       "0hOS/kTSw8BzkvaTtFzS9yQNSlor6fWF/Y+WdK+kX0q6GzhkEuc6RNLX0/s+LenbhW1PSHpDer1a\n",
       "0m2SbkrneVjSsZIukzQgaaOkMwrHXiTp8bTv/5P0+5O4BL8h6T5J2yTdIWneONfptMLyqOqx8a5Z\n",
       "A9flJEk/TPF/GegubJuXrtlWSc9IulPSkWnbJ4HXAZ+R9Jyka9P6T0valD7TA5JWTOJ6WItzMrBm\n",
       "Og/4d0AvcDjwdeBjETEP+CBwu6SD0763AN8HDgY+DlxI4e5U0kOSzqtzng8Am8kSyGHAZYVtY+9w\n",
       "zwRuBOYBDwJr0voj0nk/m843G/g0sCoi5gCvBtY2+LkFvAt4d/rcw8C1dfYdW3VT/MxHMs41k3Sp\n",
       "pDtrBiB1AXcAN6TPehvw9sL7C7gOOCr9PA98BiAirgC+A7w3Ig6KiEvSMfcDJ6T3uwW4LZ3HpgEn\n",
       "A2uWAK6NiJ9GxIvABcA3IuIugIj4JvAA8DuSjgJeAXwkInZExHeAOylUaUTECRHx5Trn+jXZl+7i\n",
       "iNgZEf82Tlzfjog1EbET+CpZ8rkqLX8FWCxpTtp3F/AySQdGxEBErJvEZ78xItZFxK+AjwDnSmqk\n",
       "iqa4T91rlpavioi31Hmf5cD+EfHpdE1uJ0u2pGOfiYivRcQLETEE/BkwttQxKt6IuDkiBiNiV0T8\n",
       "NTAT+M0GPpO1AScDa6bNhdd9wDmpumNQ0iDwWmAh2V35YEQ8X9h/4yTO8xfAj4G7U7XOh8fZd2vh\n",
       "9fPAUzEyW2N+/p6I2A68A/hPwM9SlcpkvviKn30TcACTqPpKxrtmEzkC+OmYdRtJX/CSZkn6bKqm\n",
       "2gbcC8wdk7BGlaokfVDSOknPpljm7sVnshblZGDNVPwy2QTcFBHzCj8HRcTVwM+BeZJmFfbvo8FG\n",
       "zIgYiogPRsQS4Czgv0k6dZ+Dj7g7It5I9uX7GPA/J3H4UWNe7wCeqrHfdmB2YXkhI597vGs2kZ8D\n",
       "R45ZV7ymHwCOA06JiLlkpQIxUhoYmwheB3wIOCcielO11TYab5C2FudkYGX5IvAWSW+UNENSt6SV\n",
       "ko6MiI1k1R9/KumA1DB5ZqNvLOl3JP1Guqv9JbCTrIpnr0k6TNLZqe1gB9mX9s60LR8rcFS9w4EL\n",
       "JC1NCe5jwG2FEkjRWuA8SftLegVZvX6u7jVr4CN8DxiWdEm6pr8LvLKwvYesJLRN0nzgyjHHDwBL\n",
       "CssHkbV9PCWpS9JHgTnYtOFkYKWIiCeBs4HLyapqNpHdnea/g78PvAp4BvgoWcPnbsoGj51f5+2P\n",
       "JWsIfo7sS/BvI+LeWmGwZ2mj3vJ+wPvJqlqeJutd80dp2yLgCfashim+x43AF8ju0LuAS8Zsz32E\n",
       "7Et3EFgN3Lx7p/rXLK/quVzSN2oGELED+F3gohT/ucDthV2uAQ4kK618D/jnMXF9Gvi91NPoGuCu\n",
       "9LMhffbnUzw2TahZD7eR9Hmyhq6tEfGytG4+WSNdH9kv1LkR8WzadhnwHrK7r0si4u6mBGa2jyRd\n",
       "QfZ7PZlqI7OW1sxk8DpgiKxXRZ4MriZrsLs6NfLNi4hLJb2UrKvaK8nqOb8JHBcR+1TUNzOzxjSt\n",
       "mih1Dxwcs/osRor/NwBvTa/PBr6UuhU+QdYz5JRmxWZmZqOV3WawICIG0usBYEF6fQTwZGG/J9mz\n",
       "J4SZmTVJZQ3IqWfFeHVUnhvFzKwkZU8gNiBpYURskXQ4IwOAfkrWQyP3Emr01JDkBGFmthciYvwx\n",
       "IRHRtB9gMfBIYflq4MPp9aVk0wAAvJSsv3UXcDTwOKlxe8z7RTPj3YfPubrqGByTY+rEuBxTwzHF\n",
       "RPs0rWQg6UtkoxoPkbSZrO/4VcCtyuZJf4Ks7zMRsU7SrcA6soEt/znSJzAzs+ZrWjKIiHoDhE6v\n",
       "s/+fkU2WZWZmJfMI5KnRX3UANfRXHUAN/VUHUEN/1QHU0F91AHX0Vx1ADf1VB1BDf9UB7I2mDTpr\n",
       "BkkREzWCmJnZKI18d/pxhNaxpKUroC/NGLpxe8T671YbkVl1nAysg/XNhrvStNKrPC+/dTS3GZiZ\n",
       "mZOBmZk5GZiZGU4GZmaGk4GZmeHeRGZ1jXQ9dbdTm/5cMjCrK+96mo9FMJu+nAzMzMzJwMzMnAzM\n",
       "zAwnAzMzw8nAzMxwMjAzM5wMzMwMDzqzDjQymGxoGXBv1fGYtQKXDKwD5YPJerqrjsSsVTgZmJmZ\n",
       "k4GZmTkZmJkZTgZmZoaTgZmZ4WRgZmY4GZiZGU4GZmaGk4GZmeFkYGZmOBmYmRlOBmZmhpOBmZnh\n",
       "ZGBmZjgZmJkZTgZmZkZFyUDSZZIelfSIpFskzZQ0X9IaSRsk3S2pt4rYzMw6UenJQNJi4A+AkyPi\n",
       "ZcAM4DzgUmBNRBwH3JOWzcysBFU8A/mXwA5glqSdwCzgZ8BlwOvTPjcA/Tgh2F4YecYxwMbtEeu/\n",
       "W21EZq2v9JJBRDwD/BWwiSwJPBsRa4AFETGQdhsAFpQdm00X+TOO73pqJCmY2XhKLxlIWgL8MbAY\n",
       "2AbcJumC4j4REZKizvGrC4v9EdHfnEjNzNqTpJXAyskcU0U10SuA70XE0wCS/gF4NbBF0sKI2CLp\n",
       "cGBrrYMjYnVpkZqZtaF0k9yfL0u6cqJjquhN9BiwXNKBkgScDqwD7gQuTPtcCNxRQWxmZh2p9JJB\n",
       "RDwk6UbgAWAX8EPgfwAHAbdKuhh4Aji37NjMzDpVFdVERMTVwNVjVj9DVkowM7OSeQSymZk5GZiZ\n",
       "mZOBmZnhZGBmZjgZmJkZTgZmZkZFXUvNyjO4TFoFbFwCfY9n64aWAfdWGZVZq3HJwKa5ed1pwrre\n",
       "kcnrerqrjsqs1TgZmJmZk4GZmTkZmJkZTgZmZoaTgZmZ4a6lZqOMfn5y3gU1754KfqayTVdOBmaj\n",
       "5M9PBliVuqDm3VMBVh1STVxmzeVqIjMzczIwMzMnAzMzw20GZkneSDzRvEW75zpyQ7JNKy4ZmAEj\n",
       "jcQTzVu0e66j2eXEZVYOlwysLY10AR25Qx9Z51lJzSbLJQNrU3kX0OIder7Os5KaTZaTgZmZORmY\n",
       "mZmTgZmZ4WRgZmY4GZiZGU4GZmaGk4GZmeFkYGZmOBmYmRlOBmZmhpOBmZnhZGBmZjgZmJkZFSUD\n",
       "Sb2SvippvaR1kl4lab6kNZI2SLpbUm8VsZlNlrR0hbTqTdkU2mbtqaqSwaeBb0TEUuB44DHgUmBN\n",
       "RBwH3JOWzdpArem0zdpL6clA0lzgdRHxeYCIGI6IbcBZwA1ptxuAt5Ydm5lZp6riSWdHA7+QdD1w\n",
       "AvAD4I+BBRExkPYZABZUEJu1nfyZxOAnnJntvSqqifYHTgb+LiJOBrYzpkooIgKICmKztpM/k9hP\n",
       "ODPbF1WUDJ4EnoyI76flrwKXAVskLYyILZIOB7bWOljS6sJif0T0NzNYM7N2I2klsHIyx5SeDNKX\n",
       "/WZJx0XEBuB04NH0cyHwqfTvHXWOX11WrGZm7SjdJPfny5KunOiYKkoGAP8VuFlSF/A48G5gBnCr\n",
       "pIuBJ4BzK4rNzKzjVJIMIuIh4JU1Np1edixmZlZdycCszbkXk00vTgZmeyXvxQSwyr2YrO15biIz\n",
       "M3MyMDMzJwMzM8PJwMzMcDIwMzMaSAaS7mlknZmZta+6XUslHQjMAg6VNL+waQ5wZLMDMzOz8ow3\n",
       "zuAPgfcBR5BNM517DvhMM4MyM7Ny1U0GEXENcI2kSyLi2hJjMjOzkk04AjkirpX0GmBxcf+IuLGJ\n",
       "cZmZWYkmTAaSvggcA6wFdhY2ORmY1SEtXZE9E3nj9oj13606HrOJNDI30cuBl6anj5lZQ/pmZ3MX\n",
       "rTqk6kjMGtHIOIMfAYc3OxAzM6tOIyWDQ4F1ku4HXkzrIiLOal5YZmZWpkaSwepmB2FmZtVqpDdR\n",
       "fwlxmJlZhRrpTTQE5I3HXcABwFBEzGlmYGZmVp5GSgY9+WtJ+wFnAcubGZSZmZVrUrOWRsSuiLgD\n",
       "WDXhzmZm1jYaqSZ6e2FxP7JxB883LSIzMytdI72J3sJIm8Ew8ARwdrMCMjOz8jXSZnBRCXGYmVmF\n",
       "Gnm4zSJJX5P0i/Rzu6SXlBGcmZmVo5FqouuBm4Fz0/K/T+vOaFZQZrWMTP4GMLQMuLfSgMymkUZ6\n",
       "Ex0aEddHxI708wXgsCbHZVZDPvnbXU9BT3fV0ZhNJ42UDJ6W9E7gFkDAecBTTY3KrC0NLpPyXtcu\n",
       "uVh7aaRk8G6yKqItwM+Bc9I6MxtlXrdLLtauGikZfAx4V0QMAkiaD/wl8J5mBmZmZuVppGRwQp4I\n",
       "ACLiGeDk5oVkZmZlayQZKJUG8oX5wIzmhWRmZmVrpJror4D/JelWsgbkc4BPNjUqMzMrVSMjkG+U\n",
       "9APgDWTTUrwtItY1PTIzMytNIyUDIuJR4NEmx2JmZhWZ1BTWZmY2PVWWDCTNkPSgpDvT8nxJayRt\n",
       "kHS3pN6qYjMz6zRVlgzeB6xjZHrsS4E1EXEccE9aNjOzElSSDNKsp28GPkfWQwmyx2nekF7fALy1\n",
       "gtDMzDpSVSWDvwE+BOwqrFsQEQPp9QCwoPSozMw6VOnJQNKZwNaIeJCRUsEoERGMVB+ZmVmTNdS1\n",
       "dIq9BjhL0puBbmCOpJuAAUkLI2KLpMOBrbUOlrS6sNgfEf3NDtjMrJ1IWgmsnMwxpSeDiLgcuBxA\n",
       "0uuBD0bEOyVdDVwIfCr9e0ed41eXFKqZWVtKN8n9+bKkKyc6phXGGeTVQVcBZ0jaQDba+arqQjIz\n",
       "6yxVVBPtFhH3kh4AkmZDPb3KeMzMOlWlycBs+is+/Wzj9oj13600HLM6nAzMmip/+hnAqkOqjcWs\n",
       "vlZoMzAzs4o5GZiZmZOBmZk5GZiZGU4GZmaGexNZG5CWroC+2TC0jDQuxcymlksG1gb6ZmfdM3u6\n",
       "q47EbLpyMjAzMycDMzNzm4G1qJF2Apg+bQX51BR7Tksx8nk9ZYVVwyUDa1F5O8F0aivIp6bIk1xR\n",
       "/nlrbTNrPicDMzNzMjAzMycDMzPDDcjWYjptgFmnfV5rXS4ZWIvptAFmnfZ5rVU5GZiZmZOBmZm5\n",
       "zcCsAsXnIrutwFqDk4FZ6UY9F9ltBdYSnAzMbFKWSiv6YPZG2L4+wlNnTBNuMzCzSemD2XfBU33g\n",
       "qTOmEScDMzNzNZFVb3rOUDo9LZVWHAz+P5qGnAysBeQDr8ANqq0tVQ35/2gacjWRmZk5GZiZmZOB\n",
       "mZnhZGBmZjgZmFmTLJVWLJVWVB2HNca9icysKTworb24ZGBm+8QlgOnBJQMz2ycuAUwPLhmYmVn5\n",
       "yUDSIknfkvSopB9JuiStny9pjaQNku6W1Ft2bGa2J1cDdYYqSgY7gPdHxDJgOfBeSUuBS4E1EXEc\n",
       "cE9aNrOK9cFsVwVNf6Ung4jYEhFr0+shYD1wJHAWcEPa7QbgrWXHZmYjlkorVklvGsompqvk/C6R\n",
       "lKfSNgNJi4GTgPuABRExkDYNAAsqCsvMGHluQU9FE9O5RFKuynoTSeoBbgfeFxHPSdq9LSJCUtQ5\n",
       "bnVhsT8i+psZp5lNTv4ktCFY1gOPVh1PJ5K0Elg5mWMqSQaSDiBLBDdFxB1p9YCkhRGxRdLhwNZa\n",
       "x0bE6pLCNLO9kJcoVnmq68qkm+T+fFnSlRMdU0VvIgHXAesi4prCpn8ELkyvLwTuGHusmVVjEJZN\n",
       "dfuB2wRaSxVtBq8FLgBOlfRg+lkFXAWcIWkD8Ia0bGYtYB50T3X7gdsEWkvp1UQR8V3qJ6HTy4zF\n",
       "zMwyno7CKjPy7GM/97gdpaojhmDZDmCVxEbY3rcX75U3Om+E7euzG0YrmaejsArlzz7ucUNjGypW\n",
       "HeWv97baJ290drVRdVwysFKMlAI2bo9Y7zu/ugaXZU1oMJ2uVeHOfwlAHzxeq+vpICxbKu0uHeSl\n",
       "D5cYms8lAytJXgro853fuOZ1Z9dpel2rwp1/bx/01muMngfdxdLBvpY4rHFOBmY2ZQZhWVXTV9i+\n",
       "cTWRmU2ZeR5o1racDMysaYo9jqjRY2yi7VYeVxOZWdNMNFitGYPZbO84GZiZmauJzGy04qyjuOqm\n",
       "Y7hkYGajVP0cA6uGk4GZ7bZUWtGKXUPzWVOLs5x61tOp5WoiM9stDe5quRJB3tC8Cg7J13kg2tRy\n",
       "MrCSFadb8AR1rSC/u/Z0D53NycBKlk+3ALCq5e5AO5HvsA2cDMw6TjtPFz12IjubOm5ANusw7Txd\n",
       "9NiJ7GzquGRgU87TVU+tkesJ413TRvZbKq04eBqNH/AU11PHJQNrAk9XPbXy6znRNZ14vz6YPZ3G\n",
       "D3iK66njZGBmZq4mMmtdeTfcPbvgjlcltJCLT5zHugWTaWgtzh7aM3UfYK+0UiydxCUDs5aVd8Ot\n",
       "9Yzo+lVChzHUfTbd2yZTddJKs4e2UiydxCUDayIPMLNypQF0x/fB425UnhwnA2siDzCzcqXSUO/Y\n",
       "qStsYk4Gtk8a7fZo5VjIxSeKDUuAx4vri9NS5/Xwrpu3IrcZ2D5qtNujleEwhrpnMXzA2PW1pqV2\n",
       "3bwVORl0MGnpCmnVm7K7e5uuWnVa6qmUT3Fd/Jxp6gr/bjfIyaCjeXBYJ5huA81qqVXK8dQVk+Nk\n",
       "YGZmbkC2qVR/kJRNlcau8WZ2LtnEornS0u1vKi+4UkxVw7ef4zCaSwY2hcYbJGVTo7FrPExP1wG8\n",
       "ett0rAKcqobvPpjtaqQRLhl0mNFdQfO7y+LgsKx76ERdRke2uxSwtxZy8Ykz+cmxL3LxwJZx9pvL\n",
       "5r4DufjELVy3tvYexdLC3BfztTvZcMxSnp47BC/Uu4P+BAefOMjwsfvxwovwYp29rBO4ZNBxil1B\n",
       "87vL/G6z2Jg8UZfRfLtLAXvrMIa6l9O9/TCGxr2GvQx3jb9P7dLCHIa7zqZ723h30NuZ3z2DRUO7\n",
       "6J452fhtenHJoEPs2518fue5cQn0pcFMLhE0w1w29/XxjuVb6XlhvNLCWHkpQ8wabFpwbajwvIMl\n",
       "AI1MU1F8Elw6pi2fCjdZLhl0jH25k8/vPPt69yxV2FTqZbjrQo58dqLSwlh5KaPWgLNOVnjeQW9f\n",
       "mqZionaC4pPg2vmpcJPVUiUDZbef1wAzgM9FxKcqDmnK7ctTwMY71tNCtK+dbDhmLr/u2TZmfV5K\n",
       "+BXPHTb2mGI7wlw294mB7ROdZ7z2gRfRoZdx7PJBtrr9oI5CiWFJoxPhFY8BHm60RFJFKaRlSgaS\n",
       "ZgCfAVYBLwXOl7S02qgaI2ll43vvy0Cv8Y4dXcc/uZjKctHLq45gT9XHNIfhrl6Gu/LlNdy3FEZK\n",
       "CbPYtcdNW7EdoZfhrkZKBOO1DwQzDziUM58dr/3gIqj8Wo21CY4u61yFUsK4JYzi317xmMmUSKY4\n",
       "9Ia0TDIATgF+HBFPRMQO4MvA2Y0cKB39MulNK7OfQ0r75ShYWcE5J7Ky6gD29P2W+zKpFdNBfGLx\n",
       "Mt718lkM7HFHXoYfs7klb4K+34LJ4Gk4puoYalhZdQB7o5WqiY4ENheWnwRe1dihi2fD1TtgUzd8\n",
       "aPcd1lRXqzT2fsVG1vx1Yw2vY2Jq6NjxGoZrdyO1icxjcOZp7JjxZYa7flVnn7ls7uvllN+by9FP\n",
       "bmRz39gqHhjddRRgJj85Fk45CLYtzNfVs5MNxxxQqB7ayYZjjucdFM+Vr8uqkQ4eqvU+m9m5ZJC5\n",
       "87LXDMLoKqF8XVG+fRubjppsddE2Zh61i+6Zn2D/AXgaGKmeKq4r7n8ZRy0fZOu42/cmlnqKg9Z2\n",
       "AGMbmNP8RvcOwrJU3Br1uvheY6qB6IPHZ8JR452/ONitWDXUNyWfbu+1Uskg9v7QoWG49SD4djfs\n",
       "KLxP49UqjZ2nkfcrNrLmrxtteB0VU4PHjtcwXKsbqU1kB/vFT3ihaweq+zvZy3DXETzfcyFHPlus\n",
       "4ikqdh3NXx/B8z2NdCedw3BXsXpoTqoyKp5rzjjVSLlherpmsGhoBouGhjmwC0ZXCeXrivLte9Pd\n",
       "dBfdM2ewaGg783d/vrx6qriuuH8ey3jbp7Lra3HQWq0G5rwr7jzorvW6aEw1UO9d8FQ3jFtlVxzs\n",
       "VnXVUJEi9uE7eApJWg6sjohVafkyYFexEVmq/8dpZmb1RYTG295KyWB/4P8ApwE/A+4Hzo+I9ZUG\n",
       "ZmbWAVqmzSAihiX9F+BfyLqWXudEYGZWjpYpGZiZWXVaqQF5UiR9QNIuSfNbIJaPS3pI0lpJ90ha\n",
       "VHVMAJL+QtL6FNs/SJrbAjGdI+lRSTslnVxxLKskPSbp/0r6cJWxpHg+L2lA0iNVx5KTtEjSt9L/\n",
       "2Y8kXdICMXVLui/9va2T9OdVx5STNEPSg5LurDqWnKQnJD2c4rq/3n5tmQzSl+0ZwMaqY0mujogT\n",
       "IuJE4A7gyqoDSu4GlkXECcAG4LKK4wF4BHgb8O0qg2jRQY7Xp3hayQ7g/RGxDFgOvLfq6xQRLwCn\n",
       "pr+344FT1TqPt3wfsI596h055QJYGREnRcQp9XZqy2QA/DXwJ1UHkYuI5wqLPcBTVcVSFBFrImJX\n",
       "WrwPeEmV8QBExGMRsaHqONiHQY7NEhHfoUa//ypFxJaIWJteDwHrgSOqjQoiIh8C0kXWxvhMheEA\n",
       "IOklwJuBzwHj9typwITxtF0ykHQ28GREPFx1LEWSPilpE3AhcFXV8dTwHuAbVQfRQmoNcjyyolja\n",
       "gqTFwElkNxaVkrSfpLXAAPCtiFhXdUzA3wAfAnZNtGPJAvimpAck/UG9nVqmN1GRpDXAwhqbriCr\n",
       "6nhjcfeKY7o8Iu6MiCuAKyRdSvZL8e5WiCvtcwXw64i4pVViagGtVIxveZJ6gK8C70slhEqlEu+J\n",
       "qR3sXyStjIj+quKRdCawNSIebMF5wV4bET+XdCiwRtJjqRQ6Sksmg4g4o9Z6Sb9FNjHVQ5Igq/b4\n",
       "gaRTImJrFTHVcAsl3oFPFJeki8iKrqeVEhCTulZV+ilQbOhfRFY6sDEkHQDcDnwxIu6oOp6iiNgm\n",
       "6Z+AVwD9FYbyGuAsSW8mG6k8R9KNEfGuCmMCICJ+nv79haSvkVWR7pEM2qqaKCJ+FBELIuLoiDia\n",
       "7I/35GYngolIOraweDbwYFWxFKUpwT8EnJ0a3VpNlfWqDwDHSlosqQt4B/CPFcbTkpTddV0HrIuI\n",
       "a6qOB0DSIZJ60+sDyTqTVPo3FxGXR8Si9L10HvCvrZAIJM2SdFB6PZusVqVmb7W2SgY1tEpR/88l\n",
       "PZLqMFcCH6g4ntx/J2vQXpO6lf1d1QFJepukzWQ9U/5J0j9XEUdEDAP5IMd1wFeqHuQo6UvA94Dj\n",
       "JG2WVEpV4wReC1xA1mPnwfRTdY+nw4F/TX9v9wF3RsQ9Fcc0Vqt8Ny0AvlO4Vl+PiLtr7ehBZ2Zm\n",
       "1vYlAzMzmwJOBmZm5mRgZmZOBmZmhpOBmZnhZGBmZjgZWIeT9G8lnadf0svLOJfZ3nAysI4WEa8t\n",
       "61RMYiCSJP9tWqn8C2cdTVLNSdcKo20flPQrSa8bs32GpL9MI88fSo9sRdJpkn6YHiZyXZrqYux7\n",
       "n5+2PyLpqsL6ofSea8lGaJuVxsnAOl3Nu/X0IJCTgI8C3yebJqLoPwJHASekhwfdLKmb7AE150bE\n",
       "8WQTQf5R8SBJR5BNcX4qcCLwyjQtO8As4H9HxIkRMfZ8Zk3lZGBWR5qA8GqyL/edYzafBnw2f3hQ\n",
       "RAwCvwn8JCJ+nPa5Afjt4lsCrwT6I+Lp9J43F/bZSTY7qFnpnAzM2P1wogcl/TAt9wBfAf5DRAzU\n",
       "O2zM8thSRq1ZWWvtk697ITxZmFXEycAMiIgrUtXQyWnV54HrI6Jeb6M1wB+mZykjaR7Zc6YXS1qS\n",
       "9nkno+fYD+B+4PWSDk7HngfcO7WfxmzynAys0+1xJy6pD3g78J5CI/LJY3b7HLAJeDg1+J6fnhnx\n",
       "buA2SQ8Dw8DfjzpZxBbgUuBbwFrggcLT31wqsMp4CmszM3PJwMzMnAzMzAwnAzMzw8nAzMxwMjAz\n",
       "M5wMzMwMJwMzM8PJwMzMgP8P0hZJEPVztrYAAAAASUVORK5CYII=\n"
      ],
      "text/plain": [
       "<matplotlib.figure.Figure at 0x107b5a050>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.hist(colors_data,bins=100,alpha=0.2);\n",
    "plt.hist(colors,bins=100,alpha=0.2,color='r');\n",
    "plt.xlabel('i-z color')\n",
    "plt.ylabel('count')\n",
    "plt.title('red: sims, blue: data')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "946\n",
      "198\n"
     ]
    }
   ],
   "source": [
    "print np.sum(detection_flags_list_data)\n",
    "print np.sum(detection_flags_list)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "reload(des_utils)\n",
    "zbandinfo_deep = des_utils.get_band_info(deep_list,'z')\n",
    "ibandinfo_deep = des_utils.get_band_info(deep_list,'i')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "z_cntrigs_deep, i_cntrigs_deep= des_utils.common_trignite_selector(zbandinfo_deep,ibandinfo_deep,1,0)\n",
    "deeptriggers = 0\n",
    "for h in range(0,len(z_cntrigs_deep)):\n",
    "    if np.any(z_cntrigs_deep[h]):\n",
    "        deeptriggers = deeptriggers + 1\n",
    "        "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "303\n"
     ]
    }
   ],
   "source": [
    "print deeptriggers"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": [
       "iVBORw0KGgoAAAANSUhEUgAAAXQAAAEACAYAAACj0I2EAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\n",
       "AAALEgAACxIB0t1+/AAAEA9JREFUeJzt3X2sZHddx/H3h25reVA2Dcm2lJJWw4YH5aEx0CgJF6Sk\n",
       "IdqWGIskmIrgP2JAgujWJrqJEShGwGj4Q3nIqmm1gmmoQOjN2htITIpIWwrbdcGk0kL2Vp5qgajF\n",
       "fv1jztLZy717Z87MvWfmd9+v5KbnzJwz8+29M5/9zu+c85tUFZKk5fe4oQuQJM2HgS5JjTDQJakR\n",
       "BrokNcJAl6RGGOiS1IgzBnqSDyZZT3LP2G3nJVlNciLJbUn2j913XZIvJTme5BU7Wbgk6XTbdegf\n",
       "Aq7YcNshYLWqDgJHu3WSPBt4NfDsbp/3JfETgCTtkjMGblV9GvjWhpuvBI50y0eAq7vlq4CbquqR\n",
       "qroP+DLwwvmVKkk6kz4d9IGqWu+W14ED3fJTgQfGtnsAuHCG2iRJU5hpSKRG8wacae4A5xWQpF2y\n",
       "r8c+60nOr6qTSS4AHuxu/ypw0dh2T+tuO00SQ16SeqiqnOn+PoH+UeBa4Ibuv7eM3X5jknczGmp5\n",
       "BvCZPkUtsiSHq+rw0HX0Zf3Dsv7hLHPtMFkzfMZAT3IT8BLgKUnuB34feCdwc5LXA/cB1wBU1bEk\n",
       "NwPHgO8Dv1FO5ShJu+aMgV5Vr9nirpdvsf3bgbfPWpQkaXqeJz69taELmNHa0AXMaG3oAma0NnQB\n",
       "M1obuoAZrA1dwE7Lbo+KJKllHkOXpCFMkp19DopKmtJmB7RsbDRvBrq0a8Yz3SzX/DmGLkmNMNAl\n",
       "qREGuiQ1wkCXpEYY6JLUCM9ykRbIxtMbPbVR07BDlxbOdrNSS5sz0CWpEQa6JDXCQJekRhjoktQI\n",
       "A12SGmGgS1IjDHRJaoQXFkk7ZJIv9ZXmyQ5d2lFeJKTdY6BLUiMMdElqhIEuSY0w0CWpEQa6JDXC\n",
       "QJekRhjoktQIA12SGmGgS1IjDHRJaoSBLkmNMNAlqREGuiQ1wkCXpEYY6JLUCANdkhrRO9CTXJfk\n",
       "i0nuSXJjkh9Jcl6S1SQnktyWZP88i5Ukba1XoCe5GPh14NKq+ingLOCXgUPAalUdBI5265KkXdC3\n",
       "Q/8v4BHgCUn2AU8AvgZcCRzptjkCXD1zhdICSlLjP0PXI0HPQK+qbwJ/AnyFUZB/u6pWgQNVtd5t\n",
       "tg4cmEuV0kLy+0K1WPb12SnJTwC/BVwMPAT8fZLXjm9TVVt2LkkOj62uVdVanzokqVVJVoCVqfap\n",
       "mr7DSPJq4PKqekO3/ivAZcDLgJdW1ckkFwC3V9UzN+xbVZWpn1RaIKNm5dR7J2z2mt64zend/GT7\n",
       "+F7RKZNkZ98x9OPAZUkenyTAy4FjwK3Atd021wK39Hx8SdKUeg25VNXdSf4K+CzwKPA54C+AHwVu\n",
       "TvJ64D7gmjnVKUnaRq8hl5me0CEXNcAhF+22nRxykSQtGANdkhphoEtSIwx0SWqEgS5JjTDQJakR\n",
       "vc5Dl3Q6J+jSIrBDl+ZifKIuJ+3SMAx0SWqEgS5JjTDQJakRBrokNcJAl6RGeNqiNLAznfI4fp8z\n",
       "L2o7dujSQtjqVEdPgdTkDHRJaoSBLkmNMNAlqREeFJUGMu38Lxu39yCpNjLQpcGMf9/oNNtPs4/2\n",
       "EgNdGrNVF+xsiloGBrr0Q7bqnO2Qtdg8KCpJjTDQJakRBrokNcJAl6RGGOiS1AgDXZIaYaBLUiMM\n",
       "dElqhIEuSY0w0CWpEQa6JDXCQJekRhjoktSI3oGeZH+SDye5N8mxJC9Kcl6S1SQnktyWZP88i5X0\n",
       "w5LU+M/Q9Wg4s3Tofwp8vKqeBTwXOA4cAlar6iBwtFuXtOOK06f31V6UqulfBEmeDNxZVT++4fbj\n",
       "wEuqaj3J+cBaVT1zwzblV2dpUY063MfmQz/9Cy42zodePZfns//mtcWvpmvUJNnZt0O/BPjPJB9K\n",
       "8rkkf5nkicCBqlrvtlkHDvR8fGkhOIyhZdI30PcBlwLvq6pLge+yYXilRq2/bwQtOV/GWh59v4Lu\n",
       "AeCBqvqXbv3DwHXAySTnV9XJJBcAD262c5LDY6trVbXWsw5JalKSFWBlqn36jKF3T/Yp4A1VdaIL\n",
       "6Cd0d32jqm5IcgjYX1WHNuznGLoW1sbx6J0cA3cMXdOYJDtnCfTnAe8HzgH+HXgdcBZwM/B04D7g\n",
       "mqr69rRFSUMx0LWodjTQ+zLQtcgMdC2qSbKz7xi6pIF59o028tJ/aWl5Bo5OZ6BLUiMMdElqhIEu\n",
       "SY0w0CWpEQa6JDXCQJekRhjoktQIA12SGmGgS1IjDHRJaoRzuUiN2jjXi5N2tc8OXWqa873sJQa6\n",
       "JDXCQJekRhjoktQIA12SGmGgS1IjDHRJaoSBLkmNMNAlqREGuiQ1wkCXpEYY6JLUCCfn0p63cRIr\n",
       "aVnZoUuAk1ipBXboUmP8xLF32aFLzfHTxl5loEtSIwx0SWqEgS5JjTDQJakRBrokNcJAl6RGGOiS\n",
       "1IiZAj3JWUnuTHJrt35ektUkJ5LclmT/fMqUJG1n1g79zcAxHruK4RCwWlUHgaPduiRpF/QO9CRP\n",
       "A14JvB9Id/OVwJFu+Qhw9UzVSZImNkuH/h7gbcCjY7cdqKr1bnkdODDD40uSptAr0JP8PPBgVd3J\n",
       "Y935aarKCSUkaRf1nW3xZ4Ark7wSOBf4sSR/DawnOb+qTia5AHhws52THB5bXauqtZ51SFKTkqwA\n",
       "K1PtM2qkZ3rSlwC/XVW/kORdwDeq6oYkh4D9VXVow/ZVVZt29dIQRtPNnnofhO2XJ91usfb3fbfc\n",
       "JsnOeZ2HfupV807g8iQngJd165KkXTBzhz71E9qha8HYoWsZTJKdfmORtEds/CYjA749Brq0Z2zs\n",
       "8NUa53KRpEYY6JLUCANdkhphoEtSIwx0SWqEZ7loT9p4Cp/UAjt07WHOH6e2GOiS1AgDXZIa4Ri6\n",
       "9gzHzU83/vvYbBoApwpYPnbo2mMcN3/MJL8Lf1/LxECXpEYY6JLUCANdkhrhQVE1y4Og0/N3ttwM\n",
       "dDXOOcCnN/6NR1omDrlIUiMMdElqhEMukhw7b4QduiS8gKgNBrokNcJAl6RGGOiS1AgPiqopHtzT\n",
       "XmaHrgZ5gE97k4EuSY0w0CWpEQa6JDXCQJekRhjoktQIA12SGmGgS1IjDHRJaoSBLkmN6BXoSS5K\n",
       "cnuSLyb5QpI3dbefl2Q1yYkktyXZP99yJUlb6duhPwK8paqeA1wGvDHJs4BDwGpVHQSOduuSpF3Q\n",
       "K9Cr6mRV3dUtfwe4F7gQuBI40m12BLh6HkVKkrY38xh6kouBFwB3AAeqar27ax04MOvjS5ImM9P0\n",
       "uUmeBHwEeHNVPZzkB/dVVW01lWmSw2Ora1W1NksdktSaJCvAylT7VPWbZjTJ2cA/Ap+oqvd2tx0H\n",
       "VqrqZJILgNur6pkb9quqyg8/ojS7URNx6jUdTp9Gd3x9luW9ub/v22FNkp19z3IJ8AHg2Kkw73wU\n",
       "uLZbvha4pc/jS5Km16tDT/Ji4FPA53nsn/DrgM8ANwNPB+4Drqmqb2/Y1w5dO8YO3Q69VZNkZ+8h\n",
       "l74MdO0kA91Ab9Uk2el3imrp+T2i0oiX/qsRfo+oZIeupbGxE3cIYHf5+198BrqWzPj4rnbXxjF4\n",
       "LRqHXCSpEQa6JDXCQJekRhjoktQIA12SGmGgS1IjDHRJaoSBLkmNMNAlqRFeKaql5aRc0uns0LXE\n",
       "nJBLGmegS1IjDHRJaoRj6FooTtG6fPr8zfw77ww7dC0gx8aXT5+/mX/neTPQJakRDrlI6mW700Yd\n",
       "Vtl9duiSeppkyMRhld1kh66F5sVD0uQMdC04v8dSmpSBrkE4virNn2PoGpDjq9I82aFrS5uNX9tJ\n",
       "S4vLQNc2HMOWloVDLpLUCDt0LQRPT2xHn7/l+D4O6/Vnh64F4QHSdjivy1Ds0HWaaS7nPtVJTXoK\n",
       "ol343javv78H67dmoGsTp94vm71HtrrvTPv02U7tmecBdg/Wb8ZA19w7ZztxaRgGujrz7JzHH8tO\n",
       "Stotcz8omuSKJMeTfCnJ78778ZdZkhr/2enH2riNnbOWyTSvXV/jI3MN9CRnAX8OXAE8G3hNkmfN\n",
       "8zmGlmRltkeY59H8aaYvHW03e/1DWxu6gBmtDV3AjNZ2+fkmfb9sv93yv/a3N+8O/YXAl6vqvqp6\n",
       "BPhb4KpJd07y5OTsh5Nzvzf6OfvhJBfPucZZrQxdwIxWhi5gNmtDFzCjtaELmNHa0AXMYmXoAnba\n",
       "vMfQLwTuH1t/AHjRFPsH9p0DXztntHrwu/D1H/yjM8npcbOe0jTp/jv10W4eQzG78TzSTjnTEOK0\n",
       "+89yOuNWedM3IyY9nXeWmucd6HMIif/ZB7/40Gj5oSdu/RRn+n+e9UDcpPtvdvBvswOB0x5wnPWx\n",
       "PIVQy2yRTo3drpbtnmd3a07V/Bq1JJcBh6vqim79OuDRqrphbBs7Q0nqYbvufd6Bvg/4N+DngK8B\n",
       "nwFeU1X3zu1JJEmbmuuQS1V9P8lvAp8EzgI+YJhL0u6Ya4cuSRrOoLMtJnlrkkeTnDdkHdNK8odJ\n",
       "7k5yV5KjSS4auqZpJPnjJPd2/w//kOTJQ9c0jSS/lOSLSf4vyaVD1zOJZb7gLskHk6wnuWfoWvpI\n",
       "clGS27vXzBeSvGnomqaR5Nwkd3R5cyzJO7badrBA70LwcuA/hqphBu+qqudV1fOBW4A/GLqgKd0G\n",
       "PKeqngecAK4buJ5p3QO8CvjU0IVMooEL7j7EqPZl9Qjwlqp6DnAZ8MZl+v1X1X8DL+3y5rnAS5O8\n",
       "eLNth+zQ3w38zoDP31tVPTy2+iTg60PV0kdVrVbVo93qHcDThqxnWlV1vKpODF3HFGa64G5oVfVp\n",
       "4FtD19FXVZ2sqru65e8A9wJPHbaq6VTV97rFcxgdn/zmZtsNEuhJrgIeqKrPD/H885Dkj5J8BbgW\n",
       "eOfQ9czg14CPD11E4za74O7CgWrZ07orz1/AqJFZGkkel+QuYB24vaqObbbdjs22mGQVOH+Tu65n\n",
       "9BH/FeOb71QdfZ2h/t+rqlur6nrg+iSHgPcAr9vVArexXf3dNtcD/1tVN+5qcROYpP4l4pkHCyDJ\n",
       "k4APA2/uOvWl0X2ifn53vOuTSVaqam3jdjsW6FV1+Wa3J/lJ4BLg7iQw+rj/r0leWFUP7lQ909qq\n",
       "/k3cyAJ2uNvVn+RXgVcyumZg4Uzx+18GXwXGD5xfxKhL1y5JcjbwEeBvquqWoevpq6oeSvIx4KfZ\n",
       "ZGKdXR9yqaovVNWBqrqkqi5h9MK+dJHCfDtJnjG2ehVw51C19JHkCuBtwFXdAZdltnCf7jbxWeAZ\n",
       "SS5Ocg7wauCjA9e0Z2TUOX4AOFZV7x26nmkleUqS/d3y4xmdTLJp5izCl0Qv48fRdyS5pxvTWgHe\n",
       "OnA90/ozRgdzV5PcmeR9Qxc0jSSvSnI/ozMWPpbkE0PXdCZV9X3g1AV3x4C/W6YL7pLcBPwzcDDJ\n",
       "/UkWanhxAj8LvJbR2SF3dj/LdNbOBcA/dXlzB3BrVR3dbEMvLJKkRixChy5JmgMDXZIaYaBLUiMM\n",
       "dElqhIEuSY0w0CWpEQa6JDXCQJekRvw/uNev75ze9HYAAAAASUVORK5CYII=\n"
      ],
      "text/plain": [
       "<matplotlib.figure.Figure at 0x10f25fbd0>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "reload(des_utils)\n",
    "zbandinfo_data_deep = des_utils.get_band_info(deep_list_data,'z')\n",
    "ibandinfo_data_deep = des_utils.get_band_info(deep_list_data,'i')\n",
    "ztrig_data_deep = des_utils.trigger_selector(zbandinfo_data_deep)\n",
    "itrig_data_deep = des_utils.trigger_selector(ibandinfo_data_deep)\n",
    "zsellist_data_deep, isellist_data_deep= des_utils.common_trignite_selector(zbandinfo_data_deep,ibandinfo_data_deep,1,1)\n",
    "trig_flags_list_data_deep,anytrigs_data_deep = des_utils.get_trig_flags_list(ztrig_data_deep,itrig_data_deep,zsellist_data_deep,isellist_data_deep)\n",
    "MJDtriglist_data_deep = des_utils.get_trig_MJD_list(zbandinfo_data_deep,trig_flags_list_data_deep,zsellist_data_deep)\n",
    "detection_flags_list_data_deep = des_utils.get_detection_flags_list(MJDtriglist_data_deep,zbandinfo_data_deep,ibandinfo_data_deep)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "553"
      ]
     },
     "execution_count": 23,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "colors_data_deep1 = []\n",
    "for i in range(0,len(zsellist_data_data_deep)):\n",
    "    if anytrigs_data_deep[i]:\n",
    "        imag_data_deep = -2.5*np.log(ibandinfo_data_deep[i][2][itrig_data_deep[i] & isellist_data_deep[i]][0])\n",
    "        zmag_data_deep = -2.5*np.log(zbandinfo_data_deep[i][2][ztrig_data_deep[i] & zsellist_data_deep[i]][0])\n",
    "        colors_data_deep1.append(imag_data_deep-zmag_data_deep)\n",
    "plt.hist(colors_data,bins=100);"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0"
      ]
     },
     "execution_count": 24,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "np.sum(detection_flags_list_data_deep)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 2",
   "language": "python",
   "name": "python2"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.9"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}