{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "import pylab as plt\n",
    "import pymoc\n",
    "import xidplus\n",
    "import numpy as np\n",
    "%matplotlib inline\n",
    "from astropy.table import Table, join"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "import seaborn as sns"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This notebook uses all the raw data from the CIGALE predictions and photoz catalogue, maps, PSF and relevant MOCs to create XID+ prior object and relevant tiling scheme"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Read in MOCs\n",
    "The selection functions required are the main MOC associated with the masterlist. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "Sel_func=pymoc.MOC()\n",
    "Sel_func.read('../../dmu4/dmu4_sm_GAMA-12/data/holes_GAMA-12_ukidss_k_O16_20180419_MOC.fits')\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Read in CIGALE predictions catalogue"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 39,
   "metadata": {},
   "outputs": [],
   "source": [
    "cigale=Table.read('../../dmu28/dmu28_GAMA-12/data/GAMA-12_Ldust_prediction_results.fits')\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 40,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "cigale['id'].name = 'help_id'\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 41,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<Table length=1427798>\n",
       "<table id=\"table4618292696\" class=\"table-striped table-bordered table-condensed\">\n",
       "<thead><tr><th>help_id</th><th>bayes.dust.luminosity</th><th>bayes.dust.luminosity_err</th><th>best.omegacam_u</th><th>best.omegacam_g</th><th>best.suprime_g</th><th>best.decam_g</th><th>best.gpc1_g</th><th>best.suprime_r</th><th>best.gpc1_r</th><th>best.omegacam_r</th><th>best.decam_r</th><th>best.gpc1_i</th><th>best.omegacam_i</th><th>best.suprime_i</th><th>best.gpc1_z</th><th>best.vista_z</th><th>best.suprime_z</th><th>best.decam_z</th><th>best.gpc1_y</th><th>best.suprime_y</th><th>best.vista_y</th><th>best.ukidss_y</th><th>best.ukidss_j</th><th>best.vista_j</th><th>best.ukidss_h</th><th>best.vista_h</th><th>best.vista_ks</th><th>best.ukidss_k</th></tr></thead>\n",
       "<thead><tr><th></th><th></th><th></th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th></tr></thead>\n",
       "<thead><tr><th>str28</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th></tr></thead>\n",
       "<tr><td>HELP_J112916.236+004455.105</td><td>9.65287793768e+37</td><td>1.17666784092e+38</td><td>0.000980295023014</td><td>0.00152945692592</td><td>0.00155522672771</td><td>0.00156031374351</td><td>0.0015811572678</td><td>0.00319442648463</td><td>0.00323499658544</td><td>0.00344452392723</td><td>0.00365353497765</td><td>0.00506550119885</td><td>0.00507429427253</td><td>0.00516587092685</td><td>0.00571054914042</td><td>0.005799151538</td><td>0.00589362525503</td><td>0.00606803848146</td><td>0.00634927742142</td><td>0.00641461877416</td><td>0.00659292795451</td><td>0.00665987676509</td><td>0.00809632604272</td><td>0.00811067438847</td><td>0.0106524893003</td><td>0.0106792419785</td><td>0.0122238062342</td><td>0.0124163888182</td></tr>\n",
       "<tr><td>HELP_J112916.735+004525.955</td><td>4.99604151924e+38</td><td>6.01405625316e+38</td><td>0.00228961737982</td><td>0.00416181870533</td><td>0.00425796740958</td><td>0.00427614671171</td><td>0.00434788511322</td><td>0.00959801620012</td><td>0.00965743796182</td><td>0.0103856772571</td><td>0.0110859039015</td><td>0.0168555632412</td><td>0.0170860313911</td><td>0.0175084982515</td><td>0.0223159250776</td><td>0.0228707747602</td><td>0.0235845969178</td><td>0.0248448303069</td><td>0.027143746692</td><td>0.027827712213</td><td>0.0299828410087</td><td>0.0306283344946</td><td>0.0437204990637</td><td>0.0439311275299</td><td>0.0693150570334</td><td>0.0696774039528</td><td>0.0906166212251</td><td>0.0936494380324</td></tr>\n",
       "<tr><td>HELP_J112916.868+004538.834</td><td>6.74395178932e+35</td><td>2.11732503197e+35</td><td>0.00797875041146</td><td>0.0298080393801</td><td>0.0305873658148</td><td>0.0307507600584</td><td>0.0316232736597</td><td>0.0574327371112</td><td>0.0577872312662</td><td>0.0598883825059</td><td>0.0617535017724</td><td>0.0808511378306</td><td>0.081672435631</td><td>0.0833186989953</td><td>0.101873791141</td><td>0.104195668445</td><td>0.106982286305</td><td>0.114021556098</td><td>0.125815868503</td><td>0.12929421768</td><td>0.137854161055</td><td>0.14010272913</td><td>0.163531125537</td><td>0.163858604646</td><td>0.207051791153</td><td>0.208521374093</td><td>0.195515497984</td><td>0.19135190687</td></tr>\n",
       "<tr><td>HELP_J112916.993+004528.665</td><td>1.06941218163e+37</td><td>1.18786921122e+37</td><td>0.000958880860804</td><td>0.00204094231687</td><td>0.00212978078874</td><td>0.00215053957163</td><td>0.00221212661902</td><td>0.00438028344746</td><td>0.00439633094758</td><td>0.00453091385133</td><td>0.00465871777647</td><td>0.00602217800415</td><td>0.00607213083905</td><td>0.0061811261165</td><td>0.00718156074182</td><td>0.00729981335821</td><td>0.00744050683883</td><td>0.00781017044537</td><td>0.00848844401419</td><td>0.00862261078026</td><td>0.00902709450123</td><td>0.00914284994064</td><td>0.012212649293</td><td>0.0122324762736</td><td>0.0145243616455</td><td>0.0145698906336</td><td>0.0179420683456</td><td>0.0178088588752</td></tr>\n",
       "<tr><td>HELP_J112917.687+004454.217</td><td>8.27142268505e+37</td><td>5.40913585177e+37</td><td>0.00137585337232</td><td>0.00538198184203</td><td>0.0058148696356</td><td>0.00590330344897</td><td>0.00614237388744</td><td>0.0182970759978</td><td>0.0183700985544</td><td>0.0194974307305</td><td>0.0205248944548</td><td>0.0319632835879</td><td>0.0324309746237</td><td>0.0334546451974</td><td>0.0437975562519</td><td>0.0447575259366</td><td>0.0461367954266</td><td>0.0498327120644</td><td>0.0556160601284</td><td>0.0580288181608</td><td>0.0628588008841</td><td>0.0640299734127</td><td>0.0980593652049</td><td>0.0984818456491</td><td>0.137654757191</td><td>0.138522057534</td><td>0.196027629972</td><td>0.198087233649</td></tr>\n",
       "<tr><td>HELP_J112917.704+004533.215</td><td>2.62197498446e+36</td><td>1.3701873132e+36</td><td>0.016534572991</td><td>0.0840618775116</td><td>0.0873091714541</td><td>0.0879454585475</td><td>0.0909989571154</td><td>0.195372872856</td><td>0.196313152927</td><td>0.205333356537</td><td>0.214011374115</td><td>0.30729033966</td><td>0.313252209594</td><td>0.321932819252</td><td>0.423012204598</td><td>0.435123206332</td><td>0.446370918813</td><td>0.483403907101</td><td>0.546884818768</td><td>0.566406555256</td><td>0.614831497502</td><td>0.626988793227</td><td>0.783608565337</td><td>0.785897174058</td><td>1.1039665533</td><td>1.1128944361</td><td>1.08867447396</td><td>1.06836625457</td></tr>\n",
       "<tr><td>HELP_J112918.031+004449.789</td><td>2.47229097952e+37</td><td>1.38352426361e+37</td><td>0.012063642223</td><td>0.0116270626371</td><td>0.0117979015648</td><td>0.0118222894812</td><td>0.0118366856847</td><td>0.0142852318447</td><td>0.0142829495391</td><td>0.014203194741</td><td>0.0141328553595</td><td>0.0132481827668</td><td>0.0133045894053</td><td>0.0132304404813</td><td>0.0127280861781</td><td>0.0126761603511</td><td>0.0125914341498</td><td>0.0126227125551</td><td>0.0125912927721</td><td>0.0125654024821</td><td>0.0126517230037</td><td>0.0126882500024</td><td>0.0134582586489</td><td>0.0134731307007</td><td>0.0135681982561</td><td>0.0135586563671</td><td>0.0141270715672</td><td>0.0143151654038</td></tr>\n",
       "<tr><td>HELP_J112918.084+004525.273</td><td>2.53947597807e+38</td><td>2.73391542285e+38</td><td>0.00174774603812</td><td>0.00301834946916</td><td>0.00308757079135</td><td>0.0031019382768</td><td>0.0031546914562</td><td>0.00653613466426</td><td>0.00658659521449</td><td>0.00703650922845</td><td>0.00748522603927</td><td>0.0119562680755</td><td>0.0120329268754</td><td>0.0123283305898</td><td>0.0149195883756</td><td>0.0152551995612</td><td>0.0157617182251</td><td>0.0165621229933</td><td>0.0179115010131</td><td>0.0182560765003</td><td>0.0191542214764</td><td>0.0194137045489</td><td>0.0255225608339</td><td>0.0256080599769</td><td>0.0369537359343</td><td>0.0371265850063</td><td>0.0449156944897</td><td>0.0457657667788</td></tr>\n",
       "<tr><td>HELP_J112918.563+004439.998</td><td>8.53276321612e+37</td><td>8.3862544751e+37</td><td>0.00522196151719</td><td>0.00804356017641</td><td>0.00825797029335</td><td>0.00829458625391</td><td>0.0084252846109</td><td>0.0128446362232</td><td>0.0128626795311</td><td>0.0130471347724</td><td>0.0132327659064</td><td>0.0147622416521</td><td>0.0148708568175</td><td>0.0149796962934</td><td>0.0161558633422</td><td>0.0162565301368</td><td>0.0163955971252</td><td>0.0167788897438</td><td>0.0174330543329</td><td>0.0176522400693</td><td>0.0184232429783</td><td>0.0186294326799</td><td>0.0228008123269</td><td>0.0228548116973</td><td>0.0371404379802</td><td>0.037448806764</td><td>0.0456766979859</td><td>0.0470656329732</td></tr>\n",
       "<tr><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td></tr>\n",
       "<tr><td>HELP_J122825.517-013639.700</td><td>2.41696074215e+41</td><td>2.16798644453e+41</td><td>0.00584736446799</td><td>0.025446839426</td><td>0.0269831207279</td><td>0.0272854328242</td><td>0.0283459368974</td><td>0.0850294114319</td><td>0.0843397434377</td><td>0.0926537479955</td><td>0.101908844007</td><td>0.225221925813</td><td>0.230352977847</td><td>0.242045554482</td><td>0.357928639928</td><td>0.368759134527</td><td>0.381708053067</td><td>0.417048662004</td><td>0.47759138436</td><td>0.501988640341</td><td>0.563892818398</td><td>0.580464808623</td><td>0.926010173912</td><td>0.933478251972</td><td>1.80074297894</td><td>1.82394091446</td><td>2.97207465549</td><td>3.11702279515</td></tr>\n",
       "<tr><td>HELP_J122825.688-013656.756</td><td>1.50956081754e+38</td><td>1.97889595044e+38</td><td>0.000779256277935</td><td>0.0016412152759</td><td>0.00169421321418</td><td>0.0017037741803</td><td>0.00174393813242</td><td>0.00518403739423</td><td>0.00522190886792</td><td>0.00570746347014</td><td>0.00615347395043</td><td>0.00976218890319</td><td>0.00991515631907</td><td>0.0101475602301</td><td>0.0131839973841</td><td>0.0135200605425</td><td>0.0139506747469</td><td>0.0145696204963</td><td>0.0156940139329</td><td>0.0160600078884</td><td>0.0172507443855</td><td>0.0174911329778</td><td>0.0236920817891</td><td>0.0238421357289</td><td>0.0351594379785</td><td>0.0352995234062</td><td>0.0430604888042</td><td>0.0443455638067</td></tr>\n",
       "<tr><td>HELP_J122825.798-013712.032</td><td>2.14805714074e+36</td><td>1.47680621342e+36</td><td>0.00409576465093</td><td>0.00845498801097</td><td>0.00869953108459</td><td>0.00875489990811</td><td>0.00897530312772</td><td>0.0119930133259</td><td>0.0119546556664</td><td>0.0121248049756</td><td>0.0123041597412</td><td>0.0136530209198</td><td>0.0137328097856</td><td>0.013797267028</td><td>0.0150262671606</td><td>0.0152080813341</td><td>0.0154178298916</td><td>0.015690385073</td><td>0.0161516369506</td><td>0.0163623503855</td><td>0.0171793608417</td><td>0.0173776624253</td><td>0.0192548141121</td><td>0.0192945000987</td><td>0.0211320408895</td><td>0.0211853374912</td><td>0.0220530326294</td><td>0.0215201768733</td></tr>\n",
       "<tr><td>HELP_J122825.799-013701.692</td><td>4.28322184014e+37</td><td>4.9210598117e+37</td><td>0.00273673184984</td><td>0.00369623563625</td><td>0.00373527427781</td><td>0.00374430934721</td><td>0.0037699887006</td><td>0.0073081729775</td><td>0.00740499092791</td><td>0.0078043236369</td><td>0.00807701526021</td><td>0.00907830633842</td><td>0.00915677128173</td><td>0.0091669397766</td><td>0.00981306403673</td><td>0.00987883719133</td><td>0.00996093144176</td><td>0.0100834181341</td><td>0.0102303620919</td><td>0.0103267828018</td><td>0.0105917509741</td><td>0.0106230068723</td><td>0.0121449195962</td><td>0.012204707927</td><td>0.0145782589472</td><td>0.0145867872266</td><td>0.0158722440665</td><td>0.0161073868445</td></tr>\n",
       "<tr><td>HELP_J122825.877-013755.825</td><td>4.43386678634e+38</td><td>5.09131852008e+38</td><td>0.000132283593125</td><td>0.000441118274717</td><td>0.000459182922406</td><td>0.000462519193721</td><td>0.000476413400313</td><td>0.00105571031941</td><td>0.00106010027923</td><td>0.00111447094317</td><td>0.00116748768923</td><td>0.00213236917757</td><td>0.00217436469534</td><td>0.00228141005749</td><td>0.00332512058435</td><td>0.00341618227865</td><td>0.00356574548659</td><td>0.00375611084075</td><td>0.00410227882471</td><td>0.00420900122049</td><td>0.00456678092155</td><td>0.00466049134526</td><td>0.0065452172812</td><td>0.00657780105983</td><td>0.0105144653162</td><td>0.0106233456705</td><td>0.0159286774673</td><td>0.016404330287</td></tr>\n",
       "<tr><td>HELP_J122826.138-013818.211</td><td>2.50458512808e+36</td><td>8.92150969757e+35</td><td>0.019536877316</td><td>0.0313734708931</td><td>0.0318879108337</td><td>0.0320085446911</td><td>0.0325418790243</td><td>0.0361820125396</td><td>0.0360251680581</td><td>0.036239196461</td><td>0.0363591466653</td><td>0.0372503869887</td><td>0.037510170013</td><td>0.0374997520303</td><td>0.0393017039618</td><td>0.0395854239418</td><td>0.0398142152834</td><td>0.0402904819483</td><td>0.0409471206129</td><td>0.0414053373402</td><td>0.043203387577</td><td>0.0436166760491</td><td>0.0459442420835</td><td>0.0460183793708</td><td>0.0483148508369</td><td>0.0484033109054</td><td>0.0479365115214</td><td>0.0466315409021</td></tr>\n",
       "<tr><td>HELP_J122826.531-013729.496</td><td>4.95938707418e+38</td><td>3.70590375119e+38</td><td>0.000700024693841</td><td>0.00155607983451</td><td>0.00162342258177</td><td>0.00163690540196</td><td>0.00167950363534</td><td>0.00386859751955</td><td>0.00380038540541</td><td>0.00418479610097</td><td>0.00463704791405</td><td>0.0107797907471</td><td>0.0109466937146</td><td>0.0114809534548</td><td>0.0154796026675</td><td>0.015737475243</td><td>0.0159947263059</td><td>0.0171309653971</td><td>0.0189352873112</td><td>0.0197415940529</td><td>0.0216373799162</td><td>0.0221286668903</td><td>0.0307841418842</td><td>0.0309767633824</td><td>0.0494875999151</td><td>0.0499538010009</td><td>0.0679912747101</td><td>0.0701424743446</td></tr>\n",
       "<tr><td>HELP_J122827.358-013618.662</td><td>1.20499668671e+39</td><td>1.6571915369e+39</td><td>0.000147082824979</td><td>0.000428127839656</td><td>0.000451150022557</td><td>0.000455727664497</td><td>0.000470779894341</td><td>0.00128264091654</td><td>0.00125683445288</td><td>0.00139770550129</td><td>0.00156568590424</td><td>0.00412455696042</td><td>0.004218154283</td><td>0.00446607270148</td><td>0.00660992769147</td><td>0.00677106864889</td><td>0.00694866607941</td><td>0.00752220133806</td><td>0.00847610695007</td><td>0.00890457285108</td><td>0.010007536247</td><td>0.0102860602568</td><td>0.015676571298</td><td>0.0157941921292</td><td>0.0281231825466</td><td>0.0284438210156</td><td>0.0426205295113</td><td>0.0443098305765</td></tr>\n",
       "<tr><td>HELP_J122829.089-013711.804</td><td>2.65636176092e+39</td><td>1.4112316383e+39</td><td>0.00716952899517</td><td>0.0218542495742</td><td>0.0228917105666</td><td>0.0230969601511</td><td>0.0238636222247</td><td>0.0948014718452</td><td>0.0957259882033</td><td>0.104986538139</td><td>0.113585357926</td><td>0.192024942628</td><td>0.19630936081</td><td>0.202838335013</td><td>0.287043773003</td><td>0.295300427953</td><td>0.306009196585</td><td>0.326607777237</td><td>0.36135082405</td><td>0.375843291657</td><td>0.409367279679</td><td>0.417742870715</td><td>0.640979866031</td><td>0.645353079634</td><td>1.04780865522</td><td>1.05413321538</td><td>1.4415402699</td><td>1.49841442715</td></tr>\n",
       "<tr><td>HELP_J122829.372-013659.895</td><td>6.84338249493e+37</td><td>8.08634710532e+37</td><td>0.00129959778414</td><td>0.00163818066884</td><td>0.001662871293</td><td>0.00166745277485</td><td>0.00168427989875</td><td>0.00320539268524</td><td>0.00322347666598</td><td>0.00345029898004</td><td>0.00367862099597</td><td>0.00537169699745</td><td>0.00539896452352</td><td>0.00551351465593</td><td>0.00647767980332</td><td>0.0066026518953</td><td>0.006744097808</td><td>0.00696991499573</td><td>0.00736377341283</td><td>0.00746329440193</td><td>0.00776632335094</td><td>0.00786772542183</td><td>0.00973783930807</td><td>0.00975055343527</td><td>0.0130954367986</td><td>0.0131267966265</td><td>0.0146719576511</td><td>0.0149235576397</td></tr>\n",
       "</table>"
      ],
      "text/plain": [
       "<Table length=1427798>\n",
       "          help_id            bayes.dust.luminosity ...  best.ukidss_k \n",
       "                                                   ...       mJy      \n",
       "           str28                    float64        ...     float64    \n",
       "---------------------------- --------------------- ... ---------------\n",
       "HELP_J112916.236+004455.105      9.65287793768e+37 ... 0.0124163888182\n",
       "HELP_J112916.735+004525.955      4.99604151924e+38 ... 0.0936494380324\n",
       "HELP_J112916.868+004538.834      6.74395178932e+35 ...   0.19135190687\n",
       "HELP_J112916.993+004528.665      1.06941218163e+37 ... 0.0178088588752\n",
       "HELP_J112917.687+004454.217      8.27142268505e+37 ...  0.198087233649\n",
       "HELP_J112917.704+004533.215      2.62197498446e+36 ...   1.06836625457\n",
       "HELP_J112918.031+004449.789      2.47229097952e+37 ... 0.0143151654038\n",
       "HELP_J112918.084+004525.273      2.53947597807e+38 ... 0.0457657667788\n",
       "HELP_J112918.563+004439.998      8.53276321612e+37 ... 0.0470656329732\n",
       "                         ...                   ... ...             ...\n",
       "HELP_J122825.517-013639.700      2.41696074215e+41 ...   3.11702279515\n",
       "HELP_J122825.688-013656.756      1.50956081754e+38 ... 0.0443455638067\n",
       "HELP_J122825.798-013712.032      2.14805714074e+36 ... 0.0215201768733\n",
       "HELP_J122825.799-013701.692      4.28322184014e+37 ... 0.0161073868445\n",
       "HELP_J122825.877-013755.825      4.43386678634e+38 ...  0.016404330287\n",
       "HELP_J122826.138-013818.211      2.50458512808e+36 ... 0.0466315409021\n",
       "HELP_J122826.531-013729.496      4.95938707418e+38 ... 0.0701424743446\n",
       "HELP_J122827.358-013618.662      1.20499668671e+39 ... 0.0443098305765\n",
       "HELP_J122829.089-013711.804      2.65636176092e+39 ...   1.49841442715\n",
       "HELP_J122829.372-013659.895      6.84338249493e+37 ... 0.0149235576397"
      ]
     },
     "execution_count": 41,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "cigale"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 48,
   "metadata": {},
   "outputs": [],
   "source": [
    "for i in range(0,len(cigale['help_id'])):\n",
    "    cigale['help_id'][i]=cigale['help_id'][i].strip()\n",
    "               "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Read in photoz"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 49,
   "metadata": {},
   "outputs": [],
   "source": [
    "photoz=Table.read('../../dmu24/dmu24_GAMA-12/data/master_catalogue_gama-12_20171210_photoz_20180410_r_optimised.fits')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 50,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "&lt;Table length=8591676&gt;\n",
       "<table id=\"table4617346576\" class=\"table-striped table-bordered table-condensed\">\n",
       "<thead><tr><th>help_id</th><th>RA</th><th>DEC</th><th>id</th><th>z1_median</th><th>z1_min</th><th>z1_max</th><th>z1_area</th><th>z2_median</th><th>z2_min</th><th>z2_max</th><th>z2_area</th><th>za_hb</th><th>chi_r_eazy</th><th>chi_r_atlas</th><th>chi_r_cosmos</th><th>chi_r_stellar</th><th>stellar_type</th></tr></thead>\n",
       "<thead><tr><th>str28</th><th>float64</th><th>float64</th><th>int64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>str6</th></tr></thead>\n",
       "<tr><td>HELP_J115809.625-010119.203</td><td>179.540104115</td><td>-1.02200083733</td><td>415493</td><td>0.6959</td><td>0.2271</td><td>1.1167</td><td>0.776</td><td>1.7302</td><td>1.6579</td><td>1.8051</td><td>0.023</td><td>0.716334851089</td><td>0.002686655</td><td>0.0150886325</td><td>0.0201354525</td><td>0.3918725</td><td>m3iii</td></tr>\n",
       "<tr><td>HELP_J115850.554-004148.083</td><td>179.710640849</td><td>-0.6966898527</td><td>415494</td><td>0.8684</td><td>0.4865</td><td>1.214</td><td>0.752</td><td>1.7465</td><td>1.6739</td><td>1.822</td><td>0.043</td><td>0.872100152911</td><td>0.04407825</td><td>0.02533725</td><td>0.043383625</td><td>3.34674</td><td>m4iii</td></tr>\n",
       "<tr><td>HELP_J114948.465+001147.728</td><td>177.451936426</td><td>0.196591224527</td><td>415495</td><td>0.887</td><td>0.4999</td><td>1.2677</td><td>0.724</td><td>1.6913</td><td>1.595</td><td>1.7884</td><td>0.07</td><td>0.87771645337</td><td>0.35133575</td><td>0.08987975</td><td>0.029299375</td><td>4.303075</td><td>m4v</td></tr>\n",
       "<tr><td>HELP_J120125.808-000108.647</td><td>180.357532881</td><td>-0.0190684827784</td><td>415496</td><td>0.9072</td><td>0.583</td><td>1.2542</td><td>0.75</td><td>1.6694</td><td>1.6106</td><td>1.7305</td><td>0.043</td><td>0.894666650493</td><td>0.0040718675</td><td>0.0216934375</td><td>0.028596175</td><td>1.58375425</td><td>m4iii</td></tr>\n",
       "<tr><td>HELP_J120150.572+005703.469</td><td>180.460717674</td><td>0.950963732611</td><td>415497</td><td>1.6554</td><td>1.4222</td><td>1.822</td><td>0.799</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>1.65789758247</td><td>1.186299</td><td>4.031635</td><td>2.62545</td><td>4.5122975</td><td>m7iii</td></tr>\n",
       "<tr><td>HELP_J120422.567-005913.736</td><td>181.09402893</td><td>-0.987148769669</td><td>415498</td><td>0.8732</td><td>0.5548</td><td>1.2474</td><td>0.702</td><td>1.5987</td><td>1.4735</td><td>1.7224</td><td>0.091</td><td>0.888999651538</td><td>0.00577969</td><td>0.25491375</td><td>0.28173975</td><td>0.54081825</td><td>m4iii</td></tr>\n",
       "<tr><td>HELP_J115939.276+003916.795</td><td>179.913648297</td><td>0.654665172133</td><td>415499</td><td>0.7295</td><td>0.1838</td><td>1.3019</td><td>0.615</td><td>1.9048</td><td>1.595</td><td>2.2389</td><td>0.183</td><td>0.716334851089</td><td>0.00098887325</td><td>0.0336144</td><td>0.04352285</td><td>1.25826925</td><td>rk3iii</td></tr>\n",
       "<tr><td>HELP_J115154.051+010414.748</td><td>177.975212106</td><td>1.0707633185</td><td>415500</td><td>0.9787</td><td>0.6508</td><td>1.3088</td><td>0.655</td><td>1.6778</td><td>1.5335</td><td>1.822</td><td>0.137</td><td>0.993648875848</td><td>0.199398525</td><td>0.020380665</td><td>0.048394</td><td>4.22107</td><td>m5iii</td></tr>\n",
       "<tr><td>HELP_J120457.625+010258.007</td><td>181.240105818</td><td>1.04944633237</td><td>415501</td><td>0.5947</td><td>0.1284</td><td>1.2542</td><td>0.709</td><td>1.8528</td><td>1.6819</td><td>2.0323</td><td>0.087</td><td>0.395846438326</td><td>0.006772955</td><td>0.0390784</td><td>0.2812145</td><td>0.90032575</td><td>rk0iii</td></tr>\n",
       "<tr><td>HELP_J120309.306-002609.591</td><td>180.788776986</td><td>-0.435997558437</td><td>415502</td><td>1.0283</td><td>0.7112</td><td>1.4735</td><td>0.659</td><td>1.7953</td><td>1.6262</td><td>1.9694</td><td>0.134</td><td>0.844269625561</td><td>0.0212020025</td><td>0.0434579</td><td>0.72502975</td><td>7.4572625</td><td>m4iii</td></tr>\n",
       "<tr><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td></tr>\n",
       "<tr><td>HELP_J120423.125-021109.933</td><td>181.096354437</td><td>-2.18609262636</td><td>11298533</td><td>0.6941</td><td>0.2053</td><td>1.2273</td><td>0.798</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>0.621378066675</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td></td></tr>\n",
       "<tr><td>HELP_J113919.112-000250.618</td><td>174.829631382</td><td>-0.0473937819863</td><td>11298534</td><td>1.0065</td><td>0.4513</td><td>1.6499</td><td>0.797</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>0.795213132021</td><td>0.000552244</td><td>0.0026483475</td><td>0.0208553275</td><td>9.59029</td><td>m4iii</td></tr>\n",
       "<tr><td>HELP_J122117.451+012256.663</td><td>185.322711149</td><td>1.38240630631</td><td>11298535</td><td>0.0591</td><td>0.001</td><td>0.1217</td><td>0.785</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>0.0596240452813</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td></td></tr>\n",
       "<tr><td>HELP_J115724.509+002904.083</td><td>179.352122416</td><td>0.484467577981</td><td>11298536</td><td>1.1821</td><td>0.9581</td><td>1.4005</td><td>0.789</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>1.16807747171</td><td>7.79985428571</td><td>8.00293</td><td>6.00342857143</td><td>13.0244171429</td><td>m7iii</td></tr>\n",
       "<tr><td>HELP_J120943.941-010034.499</td><td>182.433086068</td><td>-1.00958303725</td><td>11298537</td><td>0.8944</td><td>0.2569</td><td>1.595</td><td>0.799</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>0.784490130825</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td></td></tr>\n",
       "<tr><td>HELP_J115115.399-024125.001</td><td>177.814162851</td><td>-2.6902779989</td><td>11298538</td><td>0.4614</td><td>0.2345</td><td>0.6908</td><td>0.795</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>0.412672122158</td><td>0.0533738</td><td>0.19438265</td><td>0.34900625</td><td>1.772971</td><td>k3v</td></tr>\n",
       "<tr><td>HELP_J120158.333-022805.652</td><td>180.493052433</td><td>-2.46823661983</td><td>11298540</td><td>0.6168</td><td>0.3186</td><td>0.889</td><td>0.791</td><td>0.2607</td><td>0.2569</td><td>0.2645</td><td>0.001</td><td>0.660702188287</td><td>0.174539625</td><td>0.471560375</td><td>0.435440375</td><td>8.78018125</td><td>k3iii</td></tr>\n",
       "<tr><td>HELP_J120724.943+005318.682</td><td>181.853930565</td><td>0.888522915897</td><td>11298541</td><td>1.6154</td><td>1.4005</td><td>1.839</td><td>0.675</td><td>2.1784</td><td>2.0689</td><td>2.2878</td><td>0.082</td><td>1.63411931628</td><td>11.33364</td><td>12.47032</td><td>12.38136</td><td>30.7325</td><td>f0iii</td></tr>\n",
       "<tr><td>HELP_J115811.275+002726.998</td><td>179.546980197</td><td>0.457499520518</td><td>11298542</td><td>0.1416</td><td>0.0756</td><td>0.2017</td><td>0.509</td><td>0.048</td><td>0.0161</td><td>0.0724</td><td>0.279</td><td>0.0564546812376</td><td>0.573593454545</td><td>0.231829636364</td><td>0.157118272727</td><td>0.609553363636</td><td>f8v</td></tr>\n",
       "<tr><td>HELP_J114712.920-012059.469</td><td>176.803832175</td><td>-1.34985256706</td><td>11298544</td><td>1.3586</td><td>0.9061</td><td>1.822</td><td>0.795</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>1.32271677357</td><td>0.0260693428571</td><td>0.179097285714</td><td>0.398644571429</td><td>16.4203142857</td><td>m5iii</td></tr>\n",
       "</table>"
      ],
      "text/plain": [
       "<Table length=8591676>\n",
       "          help_id                 RA      ... chi_r_stellar  stellar_type\n",
       "           str28               float64    ...    float64         str6    \n",
       "--------------------------- ------------- ... -------------- ------------\n",
       "HELP_J115809.625-010119.203 179.540104115 ...      0.3918725        m3iii\n",
       "HELP_J115850.554-004148.083 179.710640849 ...        3.34674        m4iii\n",
       "HELP_J114948.465+001147.728 177.451936426 ...       4.303075          m4v\n",
       "HELP_J120125.808-000108.647 180.357532881 ...     1.58375425        m4iii\n",
       "HELP_J120150.572+005703.469 180.460717674 ...      4.5122975        m7iii\n",
       "HELP_J120422.567-005913.736  181.09402893 ...     0.54081825        m4iii\n",
       "HELP_J115939.276+003916.795 179.913648297 ...     1.25826925       rk3iii\n",
       "HELP_J115154.051+010414.748 177.975212106 ...        4.22107        m5iii\n",
       "HELP_J120457.625+010258.007 181.240105818 ...     0.90032575       rk0iii\n",
       "HELP_J120309.306-002609.591 180.788776986 ...      7.4572625        m4iii\n",
       "                        ...           ... ...            ...          ...\n",
       "HELP_J120423.125-021109.933 181.096354437 ...          -99.0             \n",
       "HELP_J113919.112-000250.618 174.829631382 ...        9.59029        m4iii\n",
       "HELP_J122117.451+012256.663 185.322711149 ...          -99.0             \n",
       "HELP_J115724.509+002904.083 179.352122416 ...  13.0244171429        m7iii\n",
       "HELP_J120943.941-010034.499 182.433086068 ...          -99.0             \n",
       "HELP_J115115.399-024125.001 177.814162851 ...       1.772971          k3v\n",
       "HELP_J120158.333-022805.652 180.493052433 ...     8.78018125        k3iii\n",
       "HELP_J120724.943+005318.682 181.853930565 ...        30.7325        f0iii\n",
       "HELP_J115811.275+002726.998 179.546980197 ... 0.609553363636          f8v\n",
       "HELP_J114712.920-012059.469 176.803832175 ...  16.4203142857        m5iii"
      ]
     },
     "execution_count": 50,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "photoz"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Join CIGALE and photoz tables"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 51,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "prior=join(cigale,photoz,keys='help_id')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 52,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "&lt;Table length=1427798&gt;\n",
       "<table id=\"table4618309248\" class=\"table-striped table-bordered table-condensed\">\n",
       "<thead><tr><th>help_id</th><th>bayes.dust.luminosity</th><th>bayes.dust.luminosity_err</th><th>best.omegacam_u</th><th>best.omegacam_g</th><th>best.suprime_g</th><th>best.decam_g</th><th>best.gpc1_g</th><th>best.suprime_r</th><th>best.gpc1_r</th><th>best.omegacam_r</th><th>best.decam_r</th><th>best.gpc1_i</th><th>best.omegacam_i</th><th>best.suprime_i</th><th>best.gpc1_z</th><th>best.vista_z</th><th>best.suprime_z</th><th>best.decam_z</th><th>best.gpc1_y</th><th>best.suprime_y</th><th>best.vista_y</th><th>best.ukidss_y</th><th>best.ukidss_j</th><th>best.vista_j</th><th>best.ukidss_h</th><th>best.vista_h</th><th>best.vista_ks</th><th>best.ukidss_k</th><th>RA</th><th>DEC</th><th>id</th><th>z1_median</th><th>z1_min</th><th>z1_max</th><th>z1_area</th><th>z2_median</th><th>z2_min</th><th>z2_max</th><th>z2_area</th><th>za_hb</th><th>chi_r_eazy</th><th>chi_r_atlas</th><th>chi_r_cosmos</th><th>chi_r_stellar</th><th>stellar_type</th></tr></thead>\n",
       "<thead><tr><th></th><th></th><th></th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th>mJy</th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th></tr></thead>\n",
       "<thead><tr><th>str28</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>int64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>str6</th></tr></thead>\n",
       "<tr><td>HELP_J112916.236+004455.105</td><td>9.65287793768e+37</td><td>1.17666784092e+38</td><td>0.000980295023014</td><td>0.00152945692592</td><td>0.00155522672771</td><td>0.00156031374351</td><td>0.0015811572678</td><td>0.00319442648463</td><td>0.00323499658544</td><td>0.00344452392723</td><td>0.00365353497765</td><td>0.00506550119885</td><td>0.00507429427253</td><td>0.00516587092685</td><td>0.00571054914042</td><td>0.005799151538</td><td>0.00589362525503</td><td>0.00606803848146</td><td>0.00634927742142</td><td>0.00641461877416</td><td>0.00659292795451</td><td>0.00665987676509</td><td>0.00809632604272</td><td>0.00811067438847</td><td>0.0106524893003</td><td>0.0106792419785</td><td>0.0122238062342</td><td>0.0124163888182</td><td>172.317651842</td><td>0.748640313787</td><td>9401913</td><td>0.5866</td><td>0.1732</td><td>1.0359</td><td>0.8</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>0.527093516722</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td></td></tr>\n",
       "<tr><td>HELP_J112916.735+004525.955</td><td>4.99604151924e+38</td><td>6.01405625316e+38</td><td>0.00228961737982</td><td>0.00416181870533</td><td>0.00425796740958</td><td>0.00427614671171</td><td>0.00434788511322</td><td>0.00959801620012</td><td>0.00965743796182</td><td>0.0103856772571</td><td>0.0110859039015</td><td>0.0168555632412</td><td>0.0170860313911</td><td>0.0175084982515</td><td>0.0223159250776</td><td>0.0228707747602</td><td>0.0235845969178</td><td>0.0248448303069</td><td>0.027143746692</td><td>0.027827712213</td><td>0.0299828410087</td><td>0.0306283344946</td><td>0.0437204990637</td><td>0.0439311275299</td><td>0.0693150570334</td><td>0.0696774039528</td><td>0.0906166212251</td><td>0.0936494380324</td><td>172.319729599</td><td>0.757209664906</td><td>10715740</td><td>0.5535</td><td>0.2162</td><td>0.9175</td><td>0.795</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>0.442606570742</td><td>0.0850949</td><td>0.0989438666667</td><td>0.0923512666667</td><td>9.717435</td><td>m3ii</td></tr>\n",
       "<tr><td>HELP_J112916.868+004538.834</td><td>6.74395178932e+35</td><td>2.11732503197e+35</td><td>0.00797875041146</td><td>0.0298080393801</td><td>0.0305873658148</td><td>0.0307507600584</td><td>0.0316232736597</td><td>0.0574327371112</td><td>0.0577872312662</td><td>0.0598883825059</td><td>0.0617535017724</td><td>0.0808511378306</td><td>0.081672435631</td><td>0.0833186989953</td><td>0.101873791141</td><td>0.104195668445</td><td>0.106982286305</td><td>0.114021556098</td><td>0.125815868503</td><td>0.12929421768</td><td>0.137854161055</td><td>0.14010272913</td><td>0.163531125537</td><td>0.163858604646</td><td>0.207051791153</td><td>0.208521374093</td><td>0.195515497984</td><td>0.19135190687</td><td>172.320283365</td><td>0.760787329227</td><td>9485122</td><td>0.0486</td><td>0.001</td><td>0.0984</td><td>0.788</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>0.0501443637558</td><td>0.215585285714</td><td>0.608331428571</td><td>2.41235428571</td><td>1.37624242857</td><td>rk1iii</td></tr>\n",
       "<tr><td>HELP_J112916.993+004528.665</td><td>1.06941218163e+37</td><td>1.18786921122e+37</td><td>0.000958880860804</td><td>0.00204094231687</td><td>0.00212978078874</td><td>0.00215053957163</td><td>0.00221212661902</td><td>0.00438028344746</td><td>0.00439633094758</td><td>0.00453091385133</td><td>0.00465871777647</td><td>0.00602217800415</td><td>0.00607213083905</td><td>0.0061811261165</td><td>0.00718156074182</td><td>0.00729981335821</td><td>0.00744050683883</td><td>0.00781017044537</td><td>0.00848844401419</td><td>0.00862261078026</td><td>0.00902709450123</td><td>0.00914284994064</td><td>0.012212649293</td><td>0.0122324762736</td><td>0.0145243616455</td><td>0.0145698906336</td><td>0.0179420683456</td><td>0.0178088588752</td><td>172.320805757</td><td>0.757962388233</td><td>10222287</td><td>0.2516</td><td>0.0628</td><td>0.482</td><td>0.805</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>0.198091697744</td><td>0.0029134775</td><td>0.01379993</td><td>0.13480295</td><td>0.14214465</td><td>m3v</td></tr>\n",
       "<tr><td>HELP_J112917.687+004454.217</td><td>8.27142268505e+37</td><td>5.40913585177e+37</td><td>0.00137585337232</td><td>0.00538198184203</td><td>0.0058148696356</td><td>0.00590330344897</td><td>0.00614237388744</td><td>0.0182970759978</td><td>0.0183700985544</td><td>0.0194974307305</td><td>0.0205248944548</td><td>0.0319632835879</td><td>0.0324309746237</td><td>0.0334546451974</td><td>0.0437975562519</td><td>0.0447575259366</td><td>0.0461367954266</td><td>0.0498327120644</td><td>0.0556160601284</td><td>0.0580288181608</td><td>0.0628588008841</td><td>0.0640299734127</td><td>0.0980593652049</td><td>0.0984818456491</td><td>0.137654757191</td><td>0.138522057534</td><td>0.196027629972</td><td>0.198087233649</td><td>172.323696237</td><td>0.748393544076</td><td>10385325</td><td>0.2803</td><td>0.1318</td><td>0.4469</td><td>0.757</td><td>0.4747</td><td>0.46</td><td>0.4909</td><td>0.029</td><td>0.291259067401</td><td>0.547619142857</td><td>0.566588571429</td><td>0.602499571429</td><td>12.1982628571</td><td>m3ii</td></tr>\n",
       "<tr><td>HELP_J112917.704+004533.215</td><td>2.62197498446e+36</td><td>1.3701873132e+36</td><td>0.016534572991</td><td>0.0840618775116</td><td>0.0873091714541</td><td>0.0879454585475</td><td>0.0909989571154</td><td>0.195372872856</td><td>0.196313152927</td><td>0.205333356537</td><td>0.214011374115</td><td>0.30729033966</td><td>0.313252209594</td><td>0.321932819252</td><td>0.423012204598</td><td>0.435123206332</td><td>0.446370918813</td><td>0.483403907101</td><td>0.546884818768</td><td>0.566406555256</td><td>0.614831497502</td><td>0.626988793227</td><td>0.783608565337</td><td>0.785897174058</td><td>1.1039665533</td><td>1.1128944361</td><td>1.08867447396</td><td>1.06836625457</td><td>172.323767131</td><td>0.759226479068</td><td>9998042</td><td>0.027</td><td>0.001</td><td>0.0565</td><td>0.757</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>0.0345329813244</td><td>1.06796071429</td><td>0.793394857143</td><td>1.14390671429</td><td>0.332926142857</td><td>k4i</td></tr>\n",
       "<tr><td>HELP_J112918.031+004449.789</td><td>2.47229097952e+37</td><td>1.38352426361e+37</td><td>0.012063642223</td><td>0.0116270626371</td><td>0.0117979015648</td><td>0.0118222894812</td><td>0.0118366856847</td><td>0.0142852318447</td><td>0.0142829495391</td><td>0.014203194741</td><td>0.0141328553595</td><td>0.0132481827668</td><td>0.0133045894053</td><td>0.0132304404813</td><td>0.0127280861781</td><td>0.0126761603511</td><td>0.0125914341498</td><td>0.0126227125551</td><td>0.0125912927721</td><td>0.0125654024821</td><td>0.0126517230037</td><td>0.0126882500024</td><td>0.0134582586489</td><td>0.0134731307007</td><td>0.0135681982561</td><td>0.0135586563671</td><td>0.0141270715672</td><td>0.0143151654038</td><td>172.325127324</td><td>0.747163665978</td><td>9960668</td><td>0.383</td><td>0.1116</td><td>0.6757</td><td>0.797</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>0.350601972086</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td></td></tr>\n",
       "<tr><td>HELP_J112918.084+004525.273</td><td>2.53947597807e+38</td><td>2.73391542285e+38</td><td>0.00174774603812</td><td>0.00301834946916</td><td>0.00308757079135</td><td>0.0031019382768</td><td>0.0031546914562</td><td>0.00653613466426</td><td>0.00658659521449</td><td>0.00703650922845</td><td>0.00748522603927</td><td>0.0119562680755</td><td>0.0120329268754</td><td>0.0123283305898</td><td>0.0149195883756</td><td>0.0152551995612</td><td>0.0157617182251</td><td>0.0165621229933</td><td>0.0179115010131</td><td>0.0182560765003</td><td>0.0191542214764</td><td>0.0194137045489</td><td>0.0255225608339</td><td>0.0256080599769</td><td>0.0369537359343</td><td>0.0371265850063</td><td>0.0449156944897</td><td>0.0457657667788</td><td>172.325348384</td><td>0.757020306313</td><td>10308562</td><td>0.6138</td><td>0.2162</td><td>1.0604</td><td>0.799</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>0.47317532873</td><td>0.1921243</td><td>0.2061254</td><td>0.2078884</td><td>8.228578</td><td>m3ii</td></tr>\n",
       "<tr><td>HELP_J112918.563+004439.998</td><td>8.53276321612e+37</td><td>8.3862544751e+37</td><td>0.00522196151719</td><td>0.00804356017641</td><td>0.00825797029335</td><td>0.00829458625391</td><td>0.0084252846109</td><td>0.0128446362232</td><td>0.0128626795311</td><td>0.0130471347724</td><td>0.0132327659064</td><td>0.0147622416521</td><td>0.0148708568175</td><td>0.0149796962934</td><td>0.0161558633422</td><td>0.0162565301368</td><td>0.0163955971252</td><td>0.0167788897438</td><td>0.0174330543329</td><td>0.0176522400693</td><td>0.0184232429783</td><td>0.0186294326799</td><td>0.0228008123269</td><td>0.0228548116973</td><td>0.0371404379802</td><td>0.037448806764</td><td>0.0456766979859</td><td>0.0470656329732</td><td>172.327347312</td><td>0.744443831151</td><td>9879505</td><td>0.3719</td><td>0.1083</td><td>0.6557</td><td>0.796</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>0.338519123815</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td></td></tr>\n",
       "<tr><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td></tr>\n",
       "<tr><td>HELP_J122825.517-013639.700</td><td>2.41696074215e+41</td><td>2.16798644453e+41</td><td>0.00584736446799</td><td>0.025446839426</td><td>0.0269831207279</td><td>0.0272854328242</td><td>0.0283459368974</td><td>0.0850294114319</td><td>0.0843397434377</td><td>0.0926537479955</td><td>0.101908844007</td><td>0.225221925813</td><td>0.230352977847</td><td>0.242045554482</td><td>0.357928639928</td><td>0.368759134527</td><td>0.381708053067</td><td>0.417048662004</td><td>0.47759138436</td><td>0.501988640341</td><td>0.563892818398</td><td>0.580464808623</td><td>0.926010173912</td><td>0.933478251972</td><td>1.80074297894</td><td>1.82394091446</td><td>2.97207465549</td><td>3.11702279515</td><td>187.106322308</td><td>-1.61102767494</td><td>9563672</td><td>0.7806</td><td>0.001</td><td>2.5969</td><td>0.8</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>0.587734064661</td><td>23.7258090909</td><td>20.9273818182</td><td>16.2084363636</td><td>6.80967090909</td><td>m4v</td></tr>\n",
       "<tr><td>HELP_J122825.688-013656.756</td><td>1.50956081754e+38</td><td>1.97889595044e+38</td><td>0.000779256277935</td><td>0.0016412152759</td><td>0.00169421321418</td><td>0.0017037741803</td><td>0.00174393813242</td><td>0.00518403739423</td><td>0.00522190886792</td><td>0.00570746347014</td><td>0.00615347395043</td><td>0.00976218890319</td><td>0.00991515631907</td><td>0.0101475602301</td><td>0.0131839973841</td><td>0.0135200605425</td><td>0.0139506747469</td><td>0.0145696204963</td><td>0.0156940139329</td><td>0.0160600078884</td><td>0.0172507443855</td><td>0.0174911329778</td><td>0.0236920817891</td><td>0.0238421357289</td><td>0.0351594379785</td><td>0.0352995234062</td><td>0.0430604888042</td><td>0.0443455638067</td><td>187.107034377</td><td>-1.6157654852</td><td>11147197</td><td>0.5226</td><td>0.3107</td><td>0.7112</td><td>0.794</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>0.573529768951</td><td>0.418236125</td><td>0.476232</td><td>0.769594125</td><td>8.9608125</td><td>m2iii</td></tr>\n",
       "<tr><td>HELP_J122825.798-013712.032</td><td>2.14805714074e+36</td><td>1.47680621342e+36</td><td>0.00409576465093</td><td>0.00845498801097</td><td>0.00869953108459</td><td>0.00875489990811</td><td>0.00897530312772</td><td>0.0119930133259</td><td>0.0119546556664</td><td>0.0121248049756</td><td>0.0123041597412</td><td>0.0136530209198</td><td>0.0137328097856</td><td>0.013797267028</td><td>0.0150262671606</td><td>0.0152080813341</td><td>0.0154178298916</td><td>0.015690385073</td><td>0.0161516369506</td><td>0.0163623503855</td><td>0.0171793608417</td><td>0.0173776624253</td><td>0.0192548141121</td><td>0.0192945000987</td><td>0.0211320408895</td><td>0.0211853374912</td><td>0.0220530326294</td><td>0.0215201768733</td><td>187.107491817</td><td>-1.62000889478</td><td>9387777</td><td>0.1547</td><td>0.115</td><td>0.2053</td><td>0.786</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>0.142023735815</td><td>0.89923625</td><td>1.65035375</td><td>2.54497</td><td>1.32692</td><td>rf8v</td></tr>\n",
       "<tr><td>HELP_J122825.799-013701.692</td><td>4.28322184014e+37</td><td>4.9210598117e+37</td><td>0.00273673184984</td><td>0.00369623563625</td><td>0.00373527427781</td><td>0.00374430934721</td><td>0.0037699887006</td><td>0.0073081729775</td><td>0.00740499092791</td><td>0.0078043236369</td><td>0.00807701526021</td><td>0.00907830633842</td><td>0.00915677128173</td><td>0.0091669397766</td><td>0.00981306403673</td><td>0.00987883719133</td><td>0.00996093144176</td><td>0.0100834181341</td><td>0.0102303620919</td><td>0.0103267828018</td><td>0.0105917509741</td><td>0.0106230068723</td><td>0.0121449195962</td><td>0.012204707927</td><td>0.0145782589472</td><td>0.0145867872266</td><td>0.0158722440665</td><td>0.0161073868445</td><td>187.107494141</td><td>-1.61713662341</td><td>11012312</td><td>0.5038</td><td>0.371</td><td>0.636</td><td>0.793</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>0.527093516722</td><td>0.4023016</td><td>0.814329</td><td>1.930058</td><td>6.207836</td><td>k3v</td></tr>\n",
       "<tr><td>HELP_J122825.877-013755.825</td><td>4.43386678634e+38</td><td>5.09131852008e+38</td><td>0.000132283593125</td><td>0.000441118274717</td><td>0.000459182922406</td><td>0.000462519193721</td><td>0.000476413400313</td><td>0.00105571031941</td><td>0.00106010027923</td><td>0.00111447094317</td><td>0.00116748768923</td><td>0.00213236917757</td><td>0.00217436469534</td><td>0.00228141005749</td><td>0.00332512058435</td><td>0.00341618227865</td><td>0.00356574548659</td><td>0.00375611084075</td><td>0.00410227882471</td><td>0.00420900122049</td><td>0.00456678092155</td><td>0.00466049134526</td><td>0.0065452172812</td><td>0.00657780105983</td><td>0.0105144653162</td><td>0.0106233456705</td><td>0.0159286774673</td><td>0.016404330287</td><td>187.10782164</td><td>-1.63217370607</td><td>10505331</td><td>0.9357</td><td>0.6707</td><td>1.234</td><td>0.555</td><td>0.5272</td><td>0.4254</td><td>0.6311</td><td>0.136</td><td>0.789843601218</td><td>0.77898</td><td>1.21890016667</td><td>0.828495333333</td><td>4.04519166667</td><td>m3ii</td></tr>\n",
       "<tr><td>HELP_J122826.138-013818.211</td><td>2.50458512808e+36</td><td>8.92150969757e+35</td><td>0.019536877316</td><td>0.0313734708931</td><td>0.0318879108337</td><td>0.0320085446911</td><td>0.0325418790243</td><td>0.0361820125396</td><td>0.0360251680581</td><td>0.036239196461</td><td>0.0363591466653</td><td>0.0372503869887</td><td>0.037510170013</td><td>0.0374997520303</td><td>0.0393017039618</td><td>0.0395854239418</td><td>0.0398142152834</td><td>0.0402904819483</td><td>0.0409471206129</td><td>0.0414053373402</td><td>0.043203387577</td><td>0.0436166760491</td><td>0.0459442420835</td><td>0.0460183793708</td><td>0.0483148508369</td><td>0.0484033109054</td><td>0.0479365115214</td><td>0.0466315409021</td><td>187.108910013</td><td>-1.63839194657</td><td>9332438</td><td>0.1356</td><td>0.1116</td><td>0.1662</td><td>0.712</td><td>0.0482</td><td>0.0407</td><td>0.0565</td><td>0.058</td><td>0.13180688451</td><td>4.282342</td><td>3.06846</td><td>10.38794</td><td>0.6893214</td><td>wf8v</td></tr>\n",
       "<tr><td>HELP_J122826.531-013729.496</td><td>4.95938707418e+38</td><td>3.70590375119e+38</td><td>0.000700024693841</td><td>0.00155607983451</td><td>0.00162342258177</td><td>0.00163690540196</td><td>0.00167950363534</td><td>0.00386859751955</td><td>0.00380038540541</td><td>0.00418479610097</td><td>0.00463704791405</td><td>0.0107797907471</td><td>0.0109466937146</td><td>0.0114809534548</td><td>0.0154796026675</td><td>0.015737475243</td><td>0.0159947263059</td><td>0.0171309653971</td><td>0.0189352873112</td><td>0.0197415940529</td><td>0.0216373799162</td><td>0.0221286668903</td><td>0.0307841418842</td><td>0.0309767633824</td><td>0.0494875999151</td><td>0.0499538010009</td><td>0.0679912747101</td><td>0.0701424743446</td><td>187.110546588</td><td>-1.62485998468</td><td>10358854</td><td>0.7838</td><td>0.5595</td><td>0.9758</td><td>0.798</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>0.844269625561</td><td>9.95825454545</td><td>7.87271636364</td><td>4.86007636364</td><td>1.66697272727</td><td>m4v</td></tr>\n",
       "<tr><td>HELP_J122827.358-013618.662</td><td>1.20499668671e+39</td><td>1.6571915369e+39</td><td>0.000147082824979</td><td>0.000428127839656</td><td>0.000451150022557</td><td>0.000455727664497</td><td>0.000470779894341</td><td>0.00128264091654</td><td>0.00125683445288</td><td>0.00139770550129</td><td>0.00156568590424</td><td>0.00412455696042</td><td>0.004218154283</td><td>0.00446607270148</td><td>0.00660992769147</td><td>0.00677106864889</td><td>0.00694866607941</td><td>0.00752220133806</td><td>0.00847610695007</td><td>0.00890457285108</td><td>0.010007536247</td><td>0.0102860602568</td><td>0.015676571298</td><td>0.0157941921292</td><td>0.0281231825466</td><td>0.0284438210156</td><td>0.0426205295113</td><td>0.0443098305765</td><td>187.113992233</td><td>-1.60518393681</td><td>10199242</td><td>0.8035</td><td>0.4688</td><td>1.1231</td><td>0.793</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>0.784490130825</td><td>0.503333833333</td><td>0.6441775</td><td>0.2474265</td><td>8.20423</td><td>m5iii</td></tr>\n",
       "<tr><td>HELP_J122829.089-013711.804</td><td>2.65636176092e+39</td><td>1.4112316383e+39</td><td>0.00716952899517</td><td>0.0218542495742</td><td>0.0228917105666</td><td>0.0230969601511</td><td>0.0238636222247</td><td>0.0948014718452</td><td>0.0957259882033</td><td>0.104986538139</td><td>0.113585357926</td><td>0.192024942628</td><td>0.19630936081</td><td>0.202838335013</td><td>0.287043773003</td><td>0.295300427953</td><td>0.306009196585</td><td>0.326607777237</td><td>0.36135082405</td><td>0.375843291657</td><td>0.409367279679</td><td>0.417742870715</td><td>0.640979866031</td><td>0.645353079634</td><td>1.04780865522</td><td>1.05413321538</td><td>1.4415402699</td><td>1.49841442715</td><td>187.121203187</td><td>-1.61994543488</td><td>9386359</td><td>0.4813</td><td>0.001</td><td>1.5108</td><td>0.799</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>-99.0</td><td>0.001</td><td>22.48</td><td>20.0994363636</td><td>12.0607181818</td><td>2.73683727273</td><td>m2p5v</td></tr>\n",
       "<tr><td>HELP_J122829.372-013659.895</td><td>6.84338249493e+37</td><td>8.08634710532e+37</td><td>0.00129959778414</td><td>0.00163818066884</td><td>0.001662871293</td><td>0.00166745277485</td><td>0.00168427989875</td><td>0.00320539268524</td><td>0.00322347666598</td><td>0.00345029898004</td><td>0.00367862099597</td><td>0.00537169699745</td><td>0.00539896452352</td><td>0.00551351465593</td><td>0.00647767980332</td><td>0.0066026518953</td><td>0.006744097808</td><td>0.00696991499573</td><td>0.00736377341283</td><td>0.00746329440193</td><td>0.00776632335094</td><td>0.00786772542183</td><td>0.00973783930807</td><td>0.00975055343527</td><td>0.0130954367986</td><td>0.0131267966265</td><td>0.0146719576511</td><td>0.0149235576397</td><td>187.122383773</td><td>-1.61663739853</td><td>10338374</td><td>0.56</td><td>0.3875</td><td>0.7685</td><td>0.772</td><td>0.2359</td><td>0.2162</td><td>0.2569</td><td>0.021</td><td>0.442606570742</td><td>0.0914002333333</td><td>0.326045333333</td><td>0.526406333333</td><td>4.244995</td><td>wk3iii</td></tr>\n",
       "</table>"
      ],
      "text/plain": [
       "<Table length=1427798>\n",
       "          help_id           bayes.dust.luminosity ... stellar_type\n",
       "                                                  ...             \n",
       "           str28                   float64        ...     str6    \n",
       "--------------------------- --------------------- ... ------------\n",
       "HELP_J112916.236+004455.105     9.65287793768e+37 ...             \n",
       "HELP_J112916.735+004525.955     4.99604151924e+38 ...         m3ii\n",
       "HELP_J112916.868+004538.834     6.74395178932e+35 ...       rk1iii\n",
       "HELP_J112916.993+004528.665     1.06941218163e+37 ...          m3v\n",
       "HELP_J112917.687+004454.217     8.27142268505e+37 ...         m3ii\n",
       "HELP_J112917.704+004533.215     2.62197498446e+36 ...          k4i\n",
       "HELP_J112918.031+004449.789     2.47229097952e+37 ...             \n",
       "HELP_J112918.084+004525.273     2.53947597807e+38 ...         m3ii\n",
       "HELP_J112918.563+004439.998     8.53276321612e+37 ...             \n",
       "                        ...                   ... ...          ...\n",
       "HELP_J122825.517-013639.700     2.41696074215e+41 ...          m4v\n",
       "HELP_J122825.688-013656.756     1.50956081754e+38 ...        m2iii\n",
       "HELP_J122825.798-013712.032     2.14805714074e+36 ...         rf8v\n",
       "HELP_J122825.799-013701.692     4.28322184014e+37 ...          k3v\n",
       "HELP_J122825.877-013755.825     4.43386678634e+38 ...         m3ii\n",
       "HELP_J122826.138-013818.211     2.50458512808e+36 ...         wf8v\n",
       "HELP_J122826.531-013729.496     4.95938707418e+38 ...          m4v\n",
       "HELP_J122827.358-013618.662     1.20499668671e+39 ...        m5iii\n",
       "HELP_J122829.089-013711.804     2.65636176092e+39 ...        m2p5v\n",
       "HELP_J122829.372-013659.895     6.84338249493e+37 ...       wk3iii"
      ]
     },
     "execution_count": 52,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "prior"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 53,
   "metadata": {},
   "outputs": [],
   "source": [
    "from astropy.cosmology import Planck15 as cosmo\n",
    "from astropy import units as u\n",
    "f_pred=prior['bayes.dust.luminosity']/(4*np.pi*cosmo.luminosity_distance(prior['z1_median']).to(u.cm))\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 54,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "prior=prior[np.isfinite(f_pred.value)][np.log10(f_pred.value[np.isfinite(f_pred.value)])>8.5]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 55,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "prior['DEC'].name='Dec'"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Read in Maps"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 56,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "\n",
    "pswfits='../../dmu19/dmu19_HELP-SPIRE-maps/data/GAMA-12_SPIRE250_v0.9.fits'#SPIRE 250 map\n",
    "pmwfits='../../dmu19/dmu19_HELP-SPIRE-maps/data/GAMA-12_SPIRE350_v0.9.fits'#SPIRE 350 map\n",
    "plwfits='../../dmu19/dmu19_HELP-SPIRE-maps/data/GAMA-12_SPIRE500_v0.9.fits'#SPIRE 500 map\n",
    "\n",
    "#output folder\n",
    "output_folder='./'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 57,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "from astropy.io import fits\n",
    "from astropy import wcs\n",
    "\n",
    "#-----250-------------\n",
    "hdulist = fits.open(pswfits)\n",
    "im250phdu=hdulist[0].header\n",
    "im250hdu=hdulist[1].header\n",
    "\n",
    "im250=hdulist[1].data*1.0E3 #convert to mJy\n",
    "nim250=hdulist[3].data*1.0E3 #convert to mJy\n",
    "w_250 = wcs.WCS(hdulist[1].header)\n",
    "pixsize250=3600.0*w_250.wcs.cdelt #pixel size (in arcseconds)\n",
    "hdulist.close()\n",
    "#-----350-------------\n",
    "hdulist = fits.open(pmwfits)\n",
    "im350phdu=hdulist[0].header\n",
    "im350hdu=hdulist[1].header\n",
    "\n",
    "im350=hdulist[1].data*1.0E3 #convert to mJy\n",
    "nim350=hdulist[3].data*1.0E3 #convert to mJy\n",
    "w_350 = wcs.WCS(hdulist[1].header)\n",
    "pixsize350=3600.0*w_350.wcs.cdelt #pixel size (in arcseconds)\n",
    "hdulist.close()\n",
    "#-----500-------------\n",
    "hdulist = fits.open(plwfits)\n",
    "im500phdu=hdulist[0].header\n",
    "im500hdu=hdulist[1].header\n",
    "im500=hdulist[1].data*1.0E3 #convert to mJy\n",
    "nim500=hdulist[3].data*1.0E3 #convert to mJy\n",
    "w_500 = wcs.WCS(hdulist[1].header)\n",
    "pixsize500=3600.0*w_500.wcs.cdelt #pixel size (in arcseconds)\n",
    "hdulist.close()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 58,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "## Set XID+ prior class"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 59,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "#---prior250--------\n",
    "prior250=xidplus.prior(im250,nim250,im250phdu,im250hdu, moc=Sel_func)#Initialise with map, uncertianty map, wcs info and primary header\n",
    "prior250.prior_cat(prior['RA'] ,prior['Dec'] ,'GAMA12_Ldust_prediction_results.fits',ID=prior['help_id'] )#Set input catalogue\n",
    "prior250.prior_bkg(-5.0,5)#Set prior on background (assumes Gaussian pdf with mu and sigma)\n",
    "#---prior350--------\n",
    "prior350=xidplus.prior(im350,nim350,im350phdu,im350hdu, moc=Sel_func)\n",
    "prior350.prior_cat(prior['RA'] ,prior['Dec'] ,'GAMA12_Ldust_prediction_results.fits',ID=prior['help_id'] )\n",
    "prior350.prior_bkg(-5.0,5)\n",
    "\n",
    "#---prior500--------\n",
    "prior500=xidplus.prior(im500,nim500,im500phdu,im500hdu, moc=Sel_func)\n",
    "prior500.prior_cat(prior['RA'] ,prior['Dec'] ,'GAMA12_Ldust_prediction_results.fits',ID=prior['help_id'] )\n",
    "prior500.prior_bkg(-5.0,5)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 60,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "#pixsize array (size of pixels in arcseconds)\n",
    "pixsize=np.array([pixsize250,pixsize350,pixsize500])\n",
    "#point response function for the three bands\n",
    "prfsize=np.array([18.15,25.15,36.3])\n",
    "#use Gaussian2DKernel to create prf (requires stddev rather than fwhm hence pfwhm/2.355)\n",
    "from astropy.convolution import Gaussian2DKernel\n",
    "\n",
    "##---------fit using Gaussian beam-----------------------\n",
    "prf250=Gaussian2DKernel(prfsize[0]/2.355,x_size=101,y_size=101)\n",
    "prf250.normalize(mode='peak')\n",
    "prf350=Gaussian2DKernel(prfsize[1]/2.355,x_size=101,y_size=101)\n",
    "prf350.normalize(mode='peak')\n",
    "prf500=Gaussian2DKernel(prfsize[2]/2.355,x_size=101,y_size=101)\n",
    "prf500.normalize(mode='peak')\n",
    "\n",
    "pind250=np.arange(0,101,1)*1.0/pixsize[0,1] #get 250 scale in terms of pixel scale of map\n",
    "pind350=np.arange(0,101,1)*1.0/pixsize[1,1] #get 350 scale in terms of pixel scale of map\n",
    "pind500=np.arange(0,101,1)*1.0/pixsize[2,1] #get 500 scale in terms of pixel scale of map\n",
    "\n",
    "prior250.set_prf(prf250.array,pind250,pind250)#requires psf as 2d grid, and x and y bins for grid (in pixel scale)\n",
    "prior350.set_prf(prf350.array,pind350,pind350)\n",
    "prior500.set_prf(prf500.array,pind500,pind500)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 61,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "----- There are 5013 tiles required for input catalogue and 100 large tiles\n",
      "writing total_bytes=1493072309...\n",
      "writing bytes [0, 1073741824)... done.\n",
      "writing bytes [1073741824, 1493072309)... done.\n"
     ]
    },
    {
     "ename": "SystemExit",
     "evalue": "",
     "output_type": "error",
     "traceback": [
      "An exception has occurred, use %tb to see the full traceback.\n",
      "\u001b[0;31mSystemExit\u001b[0m\n"
     ]
    }
   ],
   "source": [
    "import pickle\n",
    "#from moc, get healpix pixels at a given order\n",
    "from xidplus import moc_routines\n",
    "order=9\n",
    "tiles=moc_routines.get_HEALPix_pixels(order,prior250.sra,prior250.sdec,unique=True)\n",
    "order_large=6\n",
    "tiles_large=moc_routines.get_HEALPix_pixels(order_large,prior250.sra,prior250.sdec,unique=True)\n",
    "print('----- There are '+str(len(tiles))+' tiles required for input catalogue and '+str(len(tiles_large))+' large tiles')\n",
    "output_folder='./'\n",
    "outfile=output_folder+'Master_prior.pkl'\n",
    "xidplus.io.pickle_dump({'priors':[prior250,prior350,prior500],'tiles':tiles,'order':order,'version':xidplus.io.git_version()},outfile)\n",
    "outfile=output_folder+'Tiles.pkl'\n",
    "with open(outfile, 'wb') as f:\n",
    "    pickle.dump({'tiles':tiles,'order':order,'tiles_large':tiles_large,'order_large':order_large,'version':xidplus.io.git_version()},f)\n",
    "raise SystemExit()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "prior250.nsrc"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "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.6.8"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}