From b416fba6d7191424fd4114f92f43ab0cf91257c7 Mon Sep 17 00:00:00 2001 From: michaelJwilson Date: Fri, 17 Dec 2021 10:13:27 -0800 Subject: [PATCH] Add desi grab. --- desi.py | 13 ++ docs/nb/desi_qa.ipynb | 335 ++++++++++++++++++++++++++++++++++++++++++ 2 files changed, 348 insertions(+) create mode 100644 docs/nb/desi_qa.ipynb diff --git a/desi.py b/desi.py index b091148a..9e334a64 100644 --- a/desi.py +++ b/desi.py @@ -49,6 +49,8 @@ def tile2rosette(tile): # e.g. 280/28027/redrock-sv3-bright-28027.fits tabs = [] +print('Gathering DESI zs.') + for x in fpaths: zbest = Table.read(x, hdu='REDSHIFTS') fmap = Table.read(x, hdu='FIBERMAP') @@ -78,6 +80,8 @@ def tile2rosette(tile): # DESI c = SkyCoord(ra=desi_zs['TARGET_RA']*u.degree, dec=desi_zs['TARGET_DEC']*u.degree) +print('Matching DESI to GAMA Gold.') + # GAMA catalog = SkyCoord(ra=gold['RA'], dec=gold['DEC']) idx, d2d, d3d = c.match_to_catalog_3d(catalog) @@ -88,5 +92,14 @@ def tile2rosette(tile): desi_zs = hstack([desi_zs, gold_match]) desi_zs['GAMA_SEP'] = d2d +max_sep = 1.0 * u.arcsec + +print(np.mean(desi_zs['GAMA_SEP'] < max_sep)) + +opath = fpath.replace('gama', 'desi') + +print('Writing {}'.format(opath)) +desi_zs.write(opath, format='fits', overwrite=True) +print('Done.') diff --git a/docs/nb/desi_qa.ipynb b/docs/nb/desi_qa.ipynb new file mode 100644 index 00000000..79948573 --- /dev/null +++ b/docs/nb/desi_qa.ipynb @@ -0,0 +1,335 @@ +{ + "cells": [ + { + "cell_type": "code", + "execution_count": 61, + "id": "c87ba29f-4042-476a-a5e8-910b14965f28", + "metadata": {}, + "outputs": [], + "source": [ + "import numpy as np\n", + "import pylab as pl\n", + "import astropy.units as u\n", + "import matplotlib.pyplot as plt\n", + "\n", + "from astropy.table import Table" + ] + }, + { + "cell_type": "code", + "execution_count": 23, + "id": "6cbbb969-22ce-4bdd-8d89-5077558e53c5", + "metadata": {}, + "outputs": [ + { + "data": { + "text/html": [ + "Table length=113687\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "
CATAIDOBJIDRADECR_PETROZGAMANQSPECIDSURVEY_CLASSUMAG_DRED_SDSSGMAG_DRED_SDSSRMAG_DRED_SDSSIMAG_DRED_SDSSZMAG_DRED_SDSSLUMDISTDISTMODFIELD
degdegmagmagmagmagmagmag
int32int64float64float64float32float32int16bytes33int16float32float32float32float32float32float64float64bytes4
144805587724650868310204181.58795236-1.3518315219.082380.148634G12_Y3_005_153620.32484219.3540618.9679918.76835418.587437496.994443234932738.48175766513651G12
266731587726016159678736218.917384392.744148918.4931680.230334G15_Y4_212_374621.83661719.81377618.4314817.94876317.627686810.834522940065339.54466115695725G15
273412587726031176335475183.53808651.4036218618.442690.0453342dF300493619.4954318.55638118.31552118.07634717.974678140.833014875150135.74352238260108G12
543632588848898852979250212.09167127-1.038158219.3363930.185474G15_Y1_AX1_260621.22465320.07184419.19798318.88099718.871613635.161206003633139.01441982325185G15
419350587728879260270976139.758867462.825441419.2664220.199244G09_Y2_008_225622.87487820.39433719.29180718.77583718.345142688.232674973577639.18867643808146G9
575631588848899895460171140.71190199-0.0233768519.3572580.195014G09_Y2_016_056621.16621820.22690219.36333718.92238818.59357671.848339977607539.136356242171374G9
3887131588010930834112659136.68691275-1.6218226218.0364070.054744G09_Y4_210_012619.5872218.38993518.02149617.77270317.665688171.2854394976590436.168602231798175G9
320105587726032264757360216.994582511.7522162218.8480510.184454G15_Y1_EN1_247620.49638419.40386218.78446218.40174518.247435631.260532933904839.001043188163564G15
550502588848899356492084135.93356841-0.590812719.489370.250574G09_Y6_060_154522.1317920.6926619.45462818.88904618.616196892.671049796986539.75345725201822G9
...................................................
287944587726031709930027176.029674361.7682520819.5806920.224064G12_Y2_032_241522.34701520.68018219.61660219.09822318.889364785.803511847651439.47656982696438G12
288729587726031711633601179.900065571.7440209219.6505320.161764G12_Y6_058_241522.42066620.63369819.54914719.11395318.780981545.591830627066138.68433929412785G12
386186587727944562114974130.142603972.130772518.0444260.050314G09_Y1_AX2_054619.38177118.30162617.96287517.79958717.704718156.897865906031335.97808518219416G9
750405588848901003804894219.651524510.8365096919.261030.166634G15_Y3_031_089421.59509820.01795419.23169518.9960418.707644563.80003444980938.755625490361155G15
343400587726032802939285220.014128942.1401005119.5717470.137444G15_Y4_224_033521.86443520.40737719.57758519.06966418.692848456.151636802880238.29554618787875G15
136793587724650329079970176.12962846-1.8271458417.726950.105084G12_Y3_039_168620.34066818.49436617.56527117.20823116.879065341.0706804435168437.66422193815578G12
106558587722984440463597216.579216590.9072643517.9613740.06484G15_Y6_099_158619.57778518.37752517.8985217.64722817.474823204.2949332232951536.55128797852565G15
545982588848898857370290222.15670281-0.8957220519.021470.23334G15_Y3_008_097621.2782620.14663719.0059318.49265118.198502822.744540276546639.57632504501355G15
549037588848899353936322130.13443705-0.4522776117.6750130.135954G09_Y1_AS1_026620.48630318.59862917.55425817.11374916.741615450.753454751582538.269695320053714G9
561643588848899377725582184.38997161-0.4580199818.8365170.122734G12_Y1_HN1_067621.96440519.64311218.71840118.3211317.998386403.275640979747938.02800995075748G12
" + ], + "text/plain": [ + "\n", + " CATAID OBJID RA ... DISTMOD FIELD \n", + " deg ... \n", + " int32 int64 float64 ... float64 bytes4\n", + "------- ------------------ ------------ ... ------------------ ------\n", + " 144805 587724650868310204 181.58795236 ... 38.48175766513651 G12\n", + " 266731 587726016159678736 218.91738439 ... 39.54466115695725 G15\n", + " 273412 587726031176335475 183.5380865 ... 35.74352238260108 G12\n", + " 543632 588848898852979250 212.09167127 ... 39.01441982325185 G15\n", + " 419350 587728879260270976 139.75886746 ... 39.18867643808146 G9\n", + " 575631 588848899895460171 140.71190199 ... 39.136356242171374 G9\n", + "3887131 588010930834112659 136.68691275 ... 36.168602231798175 G9\n", + " 320105 587726032264757360 216.99458251 ... 39.001043188163564 G15\n", + " 550502 588848899356492084 135.93356841 ... 39.75345725201822 G9\n", + " ... ... ... ... ... ...\n", + " 287944 587726031709930027 176.02967436 ... 39.47656982696438 G12\n", + " 288729 587726031711633601 179.90006557 ... 38.68433929412785 G12\n", + " 386186 587727944562114974 130.14260397 ... 35.97808518219416 G9\n", + " 750405 588848901003804894 219.65152451 ... 38.755625490361155 G15\n", + " 343400 587726032802939285 220.01412894 ... 38.29554618787875 G15\n", + " 136793 587724650329079970 176.12962846 ... 37.66422193815578 G12\n", + " 106558 587722984440463597 216.57921659 ... 36.55128797852565 G15\n", + " 545982 588848898857370290 222.15670281 ... 39.57632504501355 G15\n", + " 549037 588848899353936322 130.13443705 ... 38.269695320053714 G9\n", + " 561643 588848899377725582 184.38997161 ... 38.02800995075748 G12" + ] + }, + "execution_count": 23, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "gama = Table.read('/global/cscratch1/sd/mjwilson/norberg//GAMA4/gama_gold.fits')\n", + "gama" + ] + }, + { + "cell_type": "code", + "execution_count": 38, + "id": "ced22df1-cf2e-4669-8661-cabe3bb79222", + "metadata": {}, + "outputs": [ + { + "data": { + "text/html": [ + "Table length=176662\n", + "
\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "
TARGETIDCHI2COEFF [10]ZZERRZWARNNPIXELSSPECTYPESUBTYPENCOEFFDELTACHI2COADD_FIBERSTATUSTARGET_RATARGET_DECPMRAPMDECREF_EPOCHFA_TARGETFA_TYPEOBJTYPESUBPRIORITYOBSCONDITIONSRELEASEBRICKIDBRICK_OBJIDMORPHTYPEFLUX_GFLUX_RFLUX_ZFLUX_IVAR_GFLUX_IVAR_RFLUX_IVAR_ZMASKBITSREF_IDREF_CATGAIA_PHOT_G_MEAN_MAGGAIA_PHOT_BP_MEAN_MAGGAIA_PHOT_RP_MEAN_MAGPARALLAXBRICKNAMEEBVFLUX_W1FLUX_W2FLUX_IVAR_W1FLUX_IVAR_W2FIBERFLUX_GFIBERFLUX_RFIBERFLUX_ZFIBERTOTFLUX_GFIBERTOTFLUX_RFIBERTOTFLUX_ZSERSICSHAPE_RSHAPE_E1SHAPE_E2PHOTSYSPRIORITY_INITNUMOBS_INITSV3_DESI_TARGETSV3_BGS_TARGETSV3_MWS_TARGETSV3_SCND_TARGETDESI_TARGETBGS_TARGETMWS_TARGETPLATE_RAPLATE_DECCOADD_NUMEXPCOADD_EXPTIMECOADD_NUMNIGHTCOADD_NUMTILEMEAN_DELTA_XRMS_DELTA_XMEAN_DELTA_YRMS_DELTA_YMEAN_FIBER_RASTD_FIBER_RAMEAN_FIBER_DECSTD_FIBER_DECMEAN_PSF_TO_FIBER_SPECFLUXCATAIDOBJIDRADECR_PETROZGAMANQSPECIDSURVEY_CLASSUMAG_DRED_SDSSGMAG_DRED_SDSSRMAG_DRED_SDSSIMAG_DRED_SDSSZMAG_DRED_SDSSLUMDISTDISTMODFIELDGAMA_SEP
degdegmagmagmagmagmagmagdeg
int64float64float64float64float64int64int64bytes6bytes20int64float64int32float64float64float32float32float32int64uint8bytes3float64int32int16int32int32bytes4float32float32float32float32float32float32int16int64bytes2float32float32float32float32bytes8float32float32float32float32float32float32float32float32float32float32float32float32float32float32float32bytes1int64int64int64int64int64int64int64int64int64float64float64int16float32int16int16float32float32float32float32float64float32float64float32float32int32int64float64float64float32float32int16bytes33int16float32float32float32float32float32float64float64bytes4float64
396277514890167088978.298949219286649.6114366625287 .. -21.7947638075317620.288896333377219874.945844833026809e-0507900GALAXY10311.12725171819330178.98596856089202-1.558652906414740.00.02015.511529215046068469761TGT0.840784064550535151690103224434996SER5.326890521.65908644.19788217.58514131.4912427.18373503602245146248298496G220.66958220.27453418.632970.01788m0150.0265700755.93694737.8923761.87226140.46116392.36342339.60965719.6096212.36342339.60965719.6096213.07258060.8982461-0.35175312-0.3856709S1021009115292150460684697613107400000178.98596856089202-1.558652906414741306.176211-0.0030.0030.0040.004178.985957191010130.0-1.55866734777889570.00.74884486185097587725041702076632178.9933246-1.5781292218.0023940.1482742dF299026619.59061618.43720417.90775517.5733317.354822495.6722064135246338.47597283802337G120.020818202744412463
396277514932061968463.270686775446221.26270114472453 .. -4.182073916657330.147565157524701481.0347126423482268e-0507914GALAXY10610.50013980269430179.00684662808823-1.58543590490194770.00.02015.511529215046068469761TGT0.399348976996212865169010322444180SER9.51099117.95502729.486227201.2548599.925620.759195000.00.00.00.01791m0150.0262674829.40115217.2494451.87085020.41616662.07454543.91636566.4315612.07529143.917676.4325131.35764962.47018380.00839696-0.72231585S1021009115292150460684697613107400000179.00684662808823-1.58543590490194771306.176211-0.0010.0010.0010.001179.0068428057580.0-1.58543947983935720.00.72444165185109587725041702076895179.00689353-1.5854494919.3921590.147484G12_Y2_044_188620.96919.94911219.28565619.0199618.83842492.7725415244080438.46323249946821G124.881250127241795e-05
396277575288097238585.26397619769604.5944099773179 .. 13.4169005604530120.275557428476754467.047210453836566e-0507909GALAXY10278.79910986125470178.76341769464946-1.3600046436645550.00.02015.511529215046068469761TGT0.94820613269253765169010323883251DEV6.070861324.14773449.132526185.2272353.75711414.891148000.00.00.00.01788m0120.02560180660.17595737.541111.61797760.432454562.12305288.44475117.1822342.12305288.44475217.1822344.01.2388442-0.087256990.32913506S1021009115292150460684697613107400000178.76341769464946-1.3600046436645551306.1762110.0010.0010.00.0178.76342159542490.0-1.36000476131918170.00.74281824144058587724650867130498178.7684619-1.3558212518.5548760.167874G12_Y1_DS1_060620.20260619.1657818.49138518.10161817.781433568.451934129680938.77346874243836G120.006552137142048083
396277514890144558650.939685828984416.96115551275443 .. 4.8697328762844970.166818942060010052.4017364386746455e-0507900GALAXY10174.267326049506660178.8741785599472-1.48742107795557590.00.02015.511529215046068469761TGT0.809587948205768551690103224432743EXP4.52767379.89430816.32525191.19542.28577426.197313204800.00.00.00.01788m0150.02722166919.929779.8082371.7115090.456015651.92880884.21501836.95462851.92880884.21501836.95462851.00.78911877-0.153922450.15664464S102100911529215046068469766554500000178.8741785599472-1.48742107795557591306.1762110.0030.003-0.0030.003178.87418994899860.0-1.48741027132866030.00.7542349185068587725041702011155178.85126445-1.5003737219.4562230.162884G12_Y3_006_161522.24576820.04062319.3464418.94675618.842228549.770598777447638.70090755230861G120.02631483428506429
396277514890136089649.093964904547904.0885701427658 .. 45.8993108400710.142905282424063531.1196772547171301e-0507890GALAXY102951.06369724869730178.8316242558849-1.48077409064269450.00.02015.511529215046068469761TGT0.808233022436806851690103224431896SER96.03242193.96255347.4955733.1665816.3601193.062924000.00.00.00.01788m0150.026511215279.6854222.626940.479031740.1530563810.2538520.71032737.10379810.25437220.71084637.1040953.65610625.2337050.096955430.26879966S1021009115292150460684697613107400000178.8316242558849-1.48077409064269451306.1762110.00.0-0.0030.003178.831623905671880.0-1.48076295496748060.00.7251477185046587725041702011010178.83164024-1.4807622317.0726910.142955371704239182866432618.81292317.49883716.81205416.3827216.084978476.196379333347338.38893044503849G121.9899652496315195e-05
396277514890145718585.085429362953293.12257447550013 .. -3.998581494788910.448271394714513050.0001176558145159310607870GALAXY10190.137391012161970178.88223322608434-1.45928741141925470.00.02015.511529215046090096731TGT0.3951255915969361751790103224432859DEV3.343165218.03135744.755554110.6660120.12917712.197061000.00.00.00.01788m0150.0273530480.8757343.2218251.41571370.384066761.00882285.44108413.5052931.00882285.44108413.5052934.01.2382280.05217518-0.030236574S1021009115292150460900967313107400000178.88223322608434-1.45928741141925471306.1762110.0010.001-0.0020.002178.882236926884220.0-1.45928006212723750.00.7411678185061587725041702011083178.86422165-1.4491656317.530560.143295371704788938680320620.03337118.35415317.55133417.11694116.805912477.4374812915807538.39458255176949G120.020655710777290157
396277514890135658576.984808493406732.433116304454 .. 31.0914751487186260.164997370399199483.415068247628022e-0507887GALAXY10182.276091990293930178.82967123169004-1.45954892076117650.00.02015.511529215046068469761TGT0.867977265195502751690103224431853SER8.23886420.36367838.18573230.07874123.50455526.720825000.00.00.00.01788m0150.02611202734.27368523.8896522.10338280.506451962.3134795.71813510.722582.3134795.71813510.722581.2575241.20307370.15278563-0.06945849S1021009115292150460684697613107400000178.82967123169004-1.45954892076117651306.1762110.0030.003-0.0020.002178.829682702401560.0-1.45954182078132530.00.7406577185078587725041702011347178.82967047-1.4595469919.2716030.165114G12_Y2_026_137621.61876120.24319.28707718.86235218.59563558.106415577495738.73358507436336G122.0754840503260076e-06
396277514890133258832.4502959484231048.708750054115 .. 16.0291555923882070.14371158419496956.772776033222978e-0507874GALAXY10429.301304315216840178.8172435272043-1.52269474696273050.00.02015.511529215046068469761TGT0.695620475092322651690103224431613SER29.00128473.189445136.0306551.66032417.8651455.380392000.00.00.00.01788m0150.027112927113.4580470.16831.04272170.297476535.92303214.94773127.7820075.92303614.94774327.782034.45315172.3590236-0.00598993150.09250433S1021009115292150460684697613107400000178.8172435272043-1.52269474696273051306.1762220.62550.88388380.8511.2027888178.819777878365939.112456-1.525923159202719111.6085670.36190474185045587725041702011001178.81725062-1.522689417.9970070.143694G12_Y1_DS2_091620.50700418.90831417.91440817.49366817.167404478.898237520215538.401216193600455G128.880441323715524e-06
396277514890132908921.767743099481134.053182832482 .. 31.8881993586457250.151047225128789763.465762556696129e-0507892GALAXY10682.54269363172350178.81570364028272-1.49808755873221160.00.02015.511529215046068469761TGT0.984465852753200851690103224431578SER24.27738661.2458116.911736171.0524380.0796220.01898203602293769574513536G220.05745719.30178617.8093870.01788m0150.026537275123.9150892.205431.3404750.39580047.30572218.43051735.1819047.30573918.43054835.1819462.60188061.5548697-0.209723030.38728243S1021009115292150460684697613107400000178.81570364028272-1.49808755873221161306.176211-0.0020.002-0.0020.002178.81569562033510.0-1.49807991500712760.00.73144543185044587725041702010999178.81571751-1.4980914218.0807290.151064G12_Y1_DS1_077620.70589318.89657217.96280517.50604617.22875505.9338778058283338.52046880573549G121.4392601038633543e-05
......................................................................................................................................................................................................................................................................................................
396278179225933618614.806522250175391.53594722569454 .. -9.248169858633920.321033549096336251.3793656245724903e-0507897GALAXY10501.735346555709840178.530808526838171.2896388928054710.00.02015.511529215046068469761TGT0.34884421718820695169010338282593EXP3.35719977.421169811.123369391.018142.0446233.122307000.00.00.00.01786p0120.01863030910.3224327.78412062.57260730.57876231.45258813.21098074.8128431.45258813.21098074.8128431.00.79827740.070216680.25059018S102000911529215046068469766553700000178.530808526838171.2896388928054711453.512110.0050.005-0.0150.015178.530829954632340.01.2896956041652420.00.73778564272264587726031174172835178.52623381.2946930519.3732470.195794G12_Y1_CX1_060621.42875120.39089219.29406418.84725218.578283674.86417690104939.14608187780051G120.006816303101388629
396278179225958089728.754658986312797.7507271138024 .. 0.00.000376614802404819837.352109143944327e-0607902STARK51983.76110015148330178.673996273031781.2708709561777538-3.2828703-1.85674832015.523058430092136939521TGT0.0775558732040940751690103382823040PSF28.65992446.63581556.854122270.054893.8079652.74489603892236699509601920G218.44221318.77518717.85983-0.0108683831786p0120.01913285415.2684777.0526782.63846180.617566622.2766136.24880244.19122722.2766136.24880244.1912270.00.00.00.0S10150092305843009213693952012800000178.673996273031781.27087095617775381453.512110.0050.005-0.010.01178.67401695649750.01.27090916443016820.00.789272297587726031174238330178.663393191.2744872518.448080.154154G12_Y1_CX2_059621.30567219.47096818.38210517.94776317.684658517.337204310988238.56886855912413G120.011200334371606792
396278179225957968480.874703511596379.58385447872996 .. 12.9025489367311560.441178790973802477.84398085487234e-0507919GALAXY10313.279093079268930178.673371080723461.30046777509584270.00.02015.511529215046090096731TGT0.1017307432727907551790103382823028DEV3.897287619.38752645.74866158.6428541.25037410.792119000.00.00.00.01786p0120.01852879195.71306659.247451.39054510.40333671.16075785.77432913.625661.16075785.77432913.625664.01.47456880.205497180.10994947S1021009115292150460900967313107400000178.673371080723461.30046777509584271453.512110.0030.003-0.0130.013178.673384407100540.01.30051681323278450.00.7196229272314587726031174238456178.684850171.3030849119.1870330.154644G12_Y2_042_385621.20609519.9037719.30225418.87267918.707264519.149195390211438.57646092626833G120.011770763633231547
396278179225964939258.8609282311521902.43554440192 .. 0.0-2.1808198749777563e-055.42258301679921e-0607874STARK52343.8705102119130178.715377365040271.3740217150071953.2078822-8.3909562015.523058430092136939521TGT0.892784764884696651690103382823725PSF15.0301734.53500751.269516417.80438116.4015555.07328803892239998044498560G218.82343719.31318718.0522210.92681461786p0120.01745743325.87915614.2322322.45575170.5830080511.67137726.81746739.81231311.67138926.817539.812380.00.00.00.0S101400923058430092136939520330000178.715377365040271.3740217150071951453.512110.0050.005-0.0060.006178.715397519507970.01.3740446186489260.00.789272298587726031174238355178.727347651.3727050216.1388150.07735579861576178231296619.1879117.24188816.18227215.64559515.254873245.9570329575875336.95429622699561G120.012039065494049993
396278179267876098756.046881139278248.54919308133788 .. -2.21024283843309280.196522838117804551.1618982387471274e-0507894GALAXY10664.75760614871980178.780813456325151.36742420951332890.00.02015.557646075230342348801TGT0.90149074530693075759010338283537SER4.4369687.99191612.144446485.03915153.908735.813503000.00.00.00.01788p0120.01864770810.57269911.9329042.4814850.559604942.00728063.61553625.49413732.0073033.61559995.49422650.50.7542698-0.205266090.06280421S10200095764607523034234880655370131072000178.780813456325151.36742420951332891453.512110.0050.005-0.0050.005178.780833521592880.01.36744339770242850.00.7536605272327587726031174303863178.78784321.3694555418.6748470.152144G12_Y1_CX1_014620.5844319.42326418.65760218.29789218.119825509.914980852482938.53748885671906G120.007315423137352873
396278239665862808786.541831851006180.80091945264329 .. -4.7271203829866780.272525763197189135.6429948866294346e-0507921GALAXY10151.202917851507660178.827317294593261.39672035761762810.00.02015.511529215046068469761TGT0.698575476597628451690103397231448SER3.209938812.00995524.542437420.80658101.1330931.163475000.00.00.00.01788p0150.0189985729.55284323.9223062.271050.54011841.72335456.44791413.1763631.72806296.45981913.1973552.12828370.63613750.1195642950.5805457S102000911529215046068469766553700000178.827317294593261.39672035761762811453.512110.0060.006-0.0040.004178.827341150625240.01.39673586753463260.00.7593031272335587726031174303986178.826678711.3941389219.0874390.152364G12_Y1_CX2_020621.93816820.3094319.18215418.74475918.41411510.7265169902798338.5409420353477G120.002659204508869675
3962782396658502711275.056098937988863.4450481066694 .. 15.3157219097009880.234411974542486771.3818521613110174e-0607896GALAXY1027430.285377502440178.76070200732021.38757269093643250.00.02015.511529215046068469761TGT0.087899380988755835169010339723195SER4.83330447.9370714.836026315.94775104.0369930.404774000.00.00.00.01788p0150.0180448833.85077265.18481542.673840.562917052.48922284.0877082.49062472.49140174.09114172.49472986.00.625308160.032758776-0.42403778S102000911529215046068469766553700000178.76070200732021.38757269093643251453.512110.0030.003-0.0090.009178.760714778402440.01.3876064146551850.00.75779647272322587726031174238759178.761774961.4053068419.3454150.157144G12_Y1_CX2_003621.67354420.3412919.34869418.88506918.473942528.409691591287238.61485387217287G120.017766558100515917
396278300105795558869.9024840657662197.099799673631 .. 0.0-0.000115948974147631956.253736084395648e-0607862STARK51735.27052121711310179.12345978280371.655529475770051-2.1955686-6.25915432015.523058430092136939521TGT0.523536075131863151690103411642659PSF18.06181736.9983352.811565373.03452109.5655454.0363303892279958420113280G218.71702619.24661818.02450.248473421791p0170.02174891919.7632966.64837272.3685210.585514814.02745528.7342341.01535814.02745528.7342341.0153580.00.00.00.0S101500923058430092136939520102400000179.12345978280371.6555294757700511491.121311-0.010.01-0.0060.006179.123421176028730.01.65555049601852970.00.789288542587726031711305740179.095287241.6562040419.4490010.227314G12_Y1_CX1_385522.11621320.67433419.39242218.87351818.491077798.758982378053639.51207877508292G120.02816885646974523
396278239707797467802.604053474963296.60473963667107 .. -4.6385407929139450.31153694104272044.0410130104760964e-0507899GALAXY10229.371178586035970179.027850292945261.59097989887162640.00.02015.511529215046068469761TGT0.25163167247824065169010339724610REX3.1120969.89797420.614836343.08878103.2465115.841923000.00.00.00.01791p0150.02184936431.285973-2.07654552.1629640.50492561.40866034.48022279.3311061.41612984.5003769.3677261.00.59886970.00.0S102000911529215046068469766553700000179.027850292945261.59097989887162641396.408511-0.0010.001-0.0010.001179.02784646492940.01.59098346125806930.00.76247287219772587726013994762400179.038878791.5771829618.6250020.179994G12_Y1_CX1_390620.5585219.31672518.51490618.07055517.807049614.254304806564338.9417410438899G120.01766041683263398
396278300105783297754.733414856717527.5785436129403 .. -7.7482408116727730.213979353661882785.8908448892770265e-0507890GALAXY10192.81308423727750179.064467520175131.6427545670253540.00.02015.511529215046068469761TGT0.731175760217351551690103411641433DEV4.55417414.53792928.321898197.991745.99648320.6546000.00.00.00.01791p0170.02169811728.21156720.6439572.10500140.50590792.52180228.05014915.6828062.52180228.05014915.6828064.00.50519973-0.03364423-0.4733275S102000911529215046068469766553700000179.064467520175131.6427545670253541427.103811-0.010.010.0050.005179.06442794022250.01.64273543510576480.00.7630033219762587726013994762281179.06448221.6427632719.7437460.213584G12_Y2_014_372522.12273820.90774519.69423719.23598718.883875744.308290212011339.3587642812638G121.706053700392514e-05
" + ], + "text/plain": [ + "\n", + " TARGETID CHI2 ... FIELD GAMA_SEP \n", + " ... deg \n", + " int64 float64 ... bytes4 float64 \n", + "----------------- ------------------ ... ------ ----------------------\n", + "39627751489016708 8978.298949219286 ... G12 0.020818202744412463\n", + "39627751493206196 8463.270686775446 ... G12 4.881250127241795e-05\n", + "39627757528809723 8585.26397619769 ... G12 0.006552137142048083\n", + "39627751489014455 8650.939685828984 ... G12 0.02631483428506429\n", + "39627751489013608 9649.093964904547 ... G12 1.9899652496315195e-05\n", + "39627751489014571 8585.085429362953 ... G12 0.020655710777290157\n", + "39627751489013565 8576.984808493406 ... G12 2.0754840503260076e-06\n", + "39627751489013325 8832.450295948423 ... G12 8.880441323715524e-06\n", + "39627751489013290 8921.76774309948 ... G12 1.4392601038633543e-05\n", + " ... ... ... ... ...\n", + "39627817922593361 8614.806522250175 ... G12 0.006816303101388629\n", + "39627817922595808 9728.75465898631 ... G12 0.011200334371606792\n", + "39627817922595796 8480.874703511596 ... G12 0.011770763633231547\n", + "39627817922596493 9258.860928231152 ... G12 0.012039065494049993\n", + "39627817926787609 8756.046881139278 ... G12 0.007315423137352873\n", + "39627823966586280 8786.541831851006 ... G12 0.002659204508869675\n", + "39627823966585027 11275.056098937988 ... G12 0.017766558100515917\n", + "39627830010579555 8869.902484065766 ... G12 0.02816885646974523\n", + "39627823970779746 7802.604053474963 ... G12 0.01766041683263398\n", + "39627830010578329 7754.733414856717 ... G12 1.706053700392514e-05" + ] + }, + "execution_count": 38, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "desi = Table.read('/global/cscratch1/sd/mjwilson/norberg//GAMA4/desi_gold.fits')\n", + "desi" + ] + }, + { + "cell_type": "code", + "execution_count": 54, + "id": "ba7a8b21-48a3-4cb0-9ce1-21f8859b77d8", + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "0.1910088191008819" + ] + }, + "execution_count": 54, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "good_match = desi['GAMA_SEP'] < 1. * u.arcsec\n", + "np.mean(good_match)" + ] + }, + { + "cell_type": "code", + "execution_count": 55, + "id": "d4f81b3b-72cf-4f9c-9187-16a79567130c", + "metadata": {}, + "outputs": [ + { + "data": { + "text/html": [ + "Table length=33744\n", + "
\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "\n", + "
TARGETIDCHI2COEFF [10]ZZERRZWARNNPIXELSSPECTYPESUBTYPENCOEFFDELTACHI2COADD_FIBERSTATUSTARGET_RATARGET_DECPMRAPMDECREF_EPOCHFA_TARGETFA_TYPEOBJTYPESUBPRIORITYOBSCONDITIONSRELEASEBRICKIDBRICK_OBJIDMORPHTYPEFLUX_GFLUX_RFLUX_ZFLUX_IVAR_GFLUX_IVAR_RFLUX_IVAR_ZMASKBITSREF_IDREF_CATGAIA_PHOT_G_MEAN_MAGGAIA_PHOT_BP_MEAN_MAGGAIA_PHOT_RP_MEAN_MAGPARALLAXBRICKNAMEEBVFLUX_W1FLUX_W2FLUX_IVAR_W1FLUX_IVAR_W2FIBERFLUX_GFIBERFLUX_RFIBERFLUX_ZFIBERTOTFLUX_GFIBERTOTFLUX_RFIBERTOTFLUX_ZSERSICSHAPE_RSHAPE_E1SHAPE_E2PHOTSYSPRIORITY_INITNUMOBS_INITSV3_DESI_TARGETSV3_BGS_TARGETSV3_MWS_TARGETSV3_SCND_TARGETDESI_TARGETBGS_TARGETMWS_TARGETPLATE_RAPLATE_DECCOADD_NUMEXPCOADD_EXPTIMECOADD_NUMNIGHTCOADD_NUMTILEMEAN_DELTA_XRMS_DELTA_XMEAN_DELTA_YRMS_DELTA_YMEAN_FIBER_RASTD_FIBER_RAMEAN_FIBER_DECSTD_FIBER_DECMEAN_PSF_TO_FIBER_SPECFLUXCATAIDOBJIDRADECR_PETROZGAMANQSPECIDSURVEY_CLASSUMAG_DRED_SDSSGMAG_DRED_SDSSRMAG_DRED_SDSSIMAG_DRED_SDSSZMAG_DRED_SDSSLUMDISTDISTMODFIELDGAMA_SEP
degdegmagmagmagmagmagmagdeg
int64float64float64float64float64int64int64bytes6bytes20int64float64int32float64float64float32float32float32int64uint8bytes3float64int32int16int32int32bytes4float32float32float32float32float32float32int16int64bytes2float32float32float32float32bytes8float32float32float32float32float32float32float32float32float32float32float32float32float32float32float32bytes1int64int64int64int64int64int64int64int64int64float64float64int16float32int16int16float32float32float32float32float64float32float64float32float32int32int64float64float64float32float32int16bytes33int16float32float32float32float32float32float64float64bytes4float64
396277514932061968463.270686775446221.26270114472453 .. -4.182073916657330.147565157524701481.0347126423482268e-0507914GALAXY10610.50013980269430179.00684662808823-1.58543590490194770.00.02015.511529215046068469761TGT0.399348976996212865169010322444180SER9.51099117.95502729.486227201.2548599.925620.759195000.00.00.00.01791m0150.0262674829.40115217.2494451.87085020.41616662.07454543.91636566.4315612.07529143.917676.4325131.35764962.47018380.00839696-0.72231585S1021009115292150460684697613107400000179.00684662808823-1.58543590490194771306.176211-0.0010.0010.0010.001179.0068428057580.0-1.58543947983935720.00.72444165185109587725041702076895179.00689353-1.5854494919.3921590.147484G12_Y2_044_188620.96919.94911219.28565619.0199618.83842492.7725415244080438.46323249946821G124.881250127241795e-05
396277514890136089649.093964904547904.0885701427658 .. 45.8993108400710.142905282424063531.1196772547171301e-0507890GALAXY102951.06369724869730178.8316242558849-1.48077409064269450.00.02015.511529215046068469761TGT0.808233022436806851690103224431896SER96.03242193.96255347.4955733.1665816.3601193.062924000.00.00.00.01788m0150.026511215279.6854222.626940.479031740.1530563810.2538520.71032737.10379810.25437220.71084637.1040953.65610625.2337050.096955430.26879966S1021009115292150460684697613107400000178.8316242558849-1.48077409064269451306.1762110.00.0-0.0030.003178.831623905671880.0-1.48076295496748060.00.7251477185046587725041702011010178.83164024-1.4807622317.0726910.142955371704239182866432618.81292317.49883716.81205416.3827216.084978476.196379333347338.38893044503849G121.9899652496315195e-05
396277514890135658576.984808493406732.433116304454 .. 31.0914751487186260.164997370399199483.415068247628022e-0507887GALAXY10182.276091990293930178.82967123169004-1.45954892076117650.00.02015.511529215046068469761TGT0.867977265195502751690103224431853SER8.23886420.36367838.18573230.07874123.50455526.720825000.00.00.00.01788m0150.02611202734.27368523.8896522.10338280.506451962.3134795.71813510.722582.3134795.71813510.722581.2575241.20307370.15278563-0.06945849S1021009115292150460684697613107400000178.82967123169004-1.45954892076117651306.1762110.0030.003-0.0020.002178.829682702401560.0-1.45954182078132530.00.7406577185078587725041702011347178.82967047-1.4595469919.2716030.165114G12_Y2_026_137621.61876120.24319.28707718.86235218.59563558.106415577495738.73358507436336G122.0754840503260076e-06
396277514890133258832.4502959484231048.708750054115 .. 16.0291555923882070.14371158419496956.772776033222978e-0507874GALAXY10429.301304315216840178.8172435272043-1.52269474696273050.00.02015.511529215046068469761TGT0.695620475092322651690103224431613SER29.00128473.189445136.0306551.66032417.8651455.380392000.00.00.00.01788m0150.027112927113.4580470.16831.04272170.297476535.92303214.94773127.7820075.92303614.94774327.782034.45315172.3590236-0.00598993150.09250433S1021009115292150460684697613107400000178.8172435272043-1.52269474696273051306.1762220.62550.88388380.8511.2027888178.819777878365939.112456-1.525923159202719111.6085670.36190474185045587725041702011001178.81725062-1.522689417.9970070.143694G12_Y1_DS2_091620.50700418.90831417.91440817.49366817.167404478.898237520215538.401216193600455G128.880441323715524e-06
396277514890132908921.767743099481134.053182832482 .. 31.8881993586457250.151047225128789763.465762556696129e-0507892GALAXY10682.54269363172350178.81570364028272-1.49808755873221160.00.02015.511529215046068469761TGT0.984465852753200851690103224431578SER24.27738661.2458116.911736171.0524380.0796220.01898203602293769574513536G220.05745719.30178617.8093870.01788m0150.026537275123.9150892.205431.3404750.39580047.30572218.43051735.1819047.30573918.43054835.1819462.60188061.5548697-0.209723030.38728243S1021009115292150460684697613107400000178.81570364028272-1.49808755873221161306.176211-0.0020.002-0.0020.002178.81569562033510.0-1.49807991500712760.00.73144543185044587725041702010999178.81571751-1.4980914218.0807290.151064G12_Y1_DS1_077620.70589318.89657217.96280517.50604617.22875505.9338778058283338.52046880573549G121.4392601038633543e-05
396277514890154738628.5112837255-5.698035012981374 .. -3.64457971351108160.174607231274479252.4759095888423346e-0507880GALAXY10291.83911931514740178.92539961412555-1.49665805449036230.00.02015.511529215046068469761TGT0.428897564600674251690103224433761SER6.507001415.07444730.169344313.85968127.3754127.39917000.00.00.00.01788m0150.02759210448.49127646.9349251.73484230.421616081.48499123.44020496.88507751.4850183.44022976.88509941.49784052.20550920.60992080.07544761S1021009115292150460684697613107400000178.92539961412555-1.49665805449036231306.176211-0.0010.001-0.0030.003178.925395392076870.0-1.49664681830335630.00.7261786185106587725041702076744178.9253884-1.4966925419.481590.17464G12_Y2_026_129521.2927820.36214619.5471419.06635918.589636593.810510604586938.86823940228173G123.6261842145141034e-05
396277514848222078712.22295193933620.3478713134267 .. 25.8921680811750720.166625378547711787.202557314946325e-0507916GALAXY10153.80181655427440178.73361188944136-1.41946752645383230.00.02015.511529215046068469761TGT0.96771737558553551690103224424799DEV10.22210327.57901450.78948174.1618780.20431517.930693204800.00.00.00.01786m0150.02417336452.52982331.9682871.44090710.37398852.66019327.177143613.217422.6606017.17831913.2197134.01.7598299-0.18761422-0.08575299S1021009115292150460684697613107400000178.73361188944136-1.41946752645383231306.176211-0.0040.0040.00.0178.733596365018490.0-1.41946701524665770.00.7286044144033587724650867065134178.73362171-1.4194531919.2737430.166564G12_Y1_DS1_084622.16289920.25073819.15458318.7194418.522907563.537553225975638.754614310121795G121.7375790599941835e-05
3962775148901361818748.8085707128050.008060513384451986 .. 0.00.143752999078508641.0088499081094984e-0507865QSO414428.1442729532720178.83191473124245-1.4088208304854346-0.078048610.85639912015.511529215046068469761TGT0.4208009796340498351690103224431906SER44.94024786.10391131.91595130.6225371.13521.23931703602297751007783296G219.1793618.59502417.4785440.0101444371788m0150.026070394428.45615784.347530.662926550.1580370816.24059731.11640547.6720516.24088731.11830347.677394.2822940.93999150.037090283-0.048969235S1021009115292150460684697613107400000178.83191473124245-1.40882083048543461453.512220.0820.118105890.3040.4292144178.83226913267931.3180335-1.40996944901834594.122830.37455797144068587724650867130551178.83194042-1.4088287217.7088850.144642dF297953619.6241518.33401917.67460317.28059817.148766482.2240026014756638.41624411749395G122.6865550455542155e-05
396277514890124319577.251358687878747.9734566799025 .. 5.0303211658376290.226259029925227341.06018178675957e-0507880GALAXY102021.1548483967780178.77832203354774-1.47079318201055640.00.02015.511529215046068469761TGT0.98580442590659945169010322443719SER30.94845468.86559123.4790269.7047331.472427.1526175000.00.00.00.01788m0150.025220707140.34793120.4697950.95843860.273347654.727732710.52001218.862844.72871610.522017518.8657072.34105852.6005730.084363590.17305252S1021009115292150460684697613107400000178.77832203354774-1.47079318201055641306.1762110.0010.001-0.0030.003178.77832554295350.0-1.47078217081514470.00.72435594185052587725041702011060178.77831555-1.4707916118.0349160.2283142dF299062619.91620618.9234318.07956517.63412517.336027802.753607342740639.52291132918415G126.669326376680688e-06
......................................................................................................................................................................................................................................................................................................
396278239665866238774.540522575378819.0996430863098 .. 20.3637844946790.226629369968461114.7377131645285844e-0607905GALAXY103094.3512668609620178.845008286977761.45862548002423640.00.02015.511529215046068469761TGT0.395588105524020851690103397231791REX8.79516914.35095119.554956352.88324110.8759127.716251000.00.00.00.01788p0150.01860162415.38016712.3989822.46291660.56004953.41908655.5788747.601913.41912675.5789457.60200071.00.72115140.00.0S102000911529215046068469766553700000178.845008286977761.45862548002423641396.408511-0.0090.0090.0040.004178.84497270680520.01.45860963219874180.00.7529406219744587726013994696994178.845014851.4585916419.5914150.22674G12_Y2_033_213520.75381720.1673619.56801619.34571319.146425796.324167341597639.50544948082788G123.447016960091287e-05
396278239665859278633.59043803811407.57460436300875 .. -1.01887395855194040.15294329711091131.1822916203372724e-0507928GALAXY10644.54144069552420178.810600861924141.43806728067917520.00.02015.511529215046068469761TGT0.953081538352244151690103397231095SER28.14318753.9167988.06263127.8162235.486519.82751000.00.00.00.01788p0150.01817997475.3514555.927581.52118220.415324183.94113837.550442712.3321853.94118627.550509512.3322580.62743431.74220010.065064530.079854846S1021009115292150460684697613107400000178.810600861924141.43806728067917521315.6995110.0080.008-0.020.02178.810634519077640.01.438142944755040.00.68552214219732587726013994696921178.810612971.4380612118.2475740.152894G12_Y6_060_335619.89840518.83554518.14762717.78190617.56995512.682442416577638.54924222312786G121.3541281718852934e-05
396278179225948838570.0262544751171007.7877356615194 .. 38.479630673168280.07066472820578813.2095733086307396e-0607850GALAXY103118.0185405611990178.620250093182341.2489956268887840.00.02015.511529215046068469761TGT0.716970199483690151690103382822115REX19.63025329.33417938.1727177.2663653.2545315.078373000.00.00.00.01786p0120.01800961620.67508324.6434232.1472760.511603364.48960546.70897538.730424.50343946.72913558.7499141.01.1973390.00.0S1021009115292150460684697613107400000178.620250093182341.2489956268887841315.6995110.0080.008-0.0260.026178.620284482892430.01.2490954084930050.00.68244445272259587726031174172810178.620300731.249039518.8851870.070734G12_Y1_CX2_042620.1784719.20386918.86583318.58326718.461788223.9705876754210336.75095494772087G126.699043876701257e-05
396278239665862648402.992509886855316.1233490956164 .. -7.3370024041488520.152594985266718940.0001667453632356858507883GALAXY1098.068192416569220178.826658115052421.39412912011138260.00.02015.511529215046068469761TGT0.662801099905710951690103397231432SER8.3637521.73786241.122578145.6683534.63345311.029137000.00.00.00.01788p0150.01902444329.22277318.1213021.57427510.382850972.44656976.35877312.0292022.4498696.36415512.0339956.01.65813480.08906218-0.06584132S1021009115292150460684697613107400000178.826658115052421.39412912011138261427.103811-0.0060.0060.0130.013178.82663319613970.01.39407950263908640.00.71608126272335587726031174303986178.826678711.3941389219.0874390.152364G12_Y1_CX2_020621.93816820.3094319.18215418.74475918.41411510.7265169902798338.5409420353477G122.2802162403838536e-05
396278239665869909520.32962005958505.7526165386138 .. 5.5890337180031240.226212586484777935.2774781728126445e-0507908GALAXY10161.030135642737150178.865367073631.4958206749360440.00.02015.511529215046068469761TGT0.2555860476802378551690103397232158SER6.21133519.39540137.90628251.8898664.80450416.27137000.00.00.00.01788p0150.01947448634.7609720.24481.40943940.396826452.02951966.33734112.3856692.03024366.339036512.3915496.01.40791640.092984820.117882766S1021009115292150460684697613107400000178.865367073631.4958206749360441427.103811-0.0060.0060.0060.006178.865342983443640.01.49579769115265520.00.73033994219754587726013994697154178.865374481.4958012719.5254480.226473G12_Y6_060_337522.1572820.61704419.52377319.01620518.820854795.406539056982139.50294578505199G122.076941215168854e-05
396278239665859698679.489034175873869.131382491431 .. 4.1706380610367660.1435550765739276.7357141023634935e-0607928GALAXY101283.41035270690920178.81225530865241.41448400573949360.00.02015.511529215046068469761TGT0.380749825023218151690103397231137REX9.51823916.8200326.069818368.20688114.38473530.660776000.00.00.00.01788p0150.01852253528.02073727.5574552.2534430.5370294.1096857.262374411.256154.11100867.269427311.27465251.00.63595290.00.0S1021009115292150460684697613107400000178.81225530865241.41448400573949361427.103811-0.0040.0040.0080.008178.81223876940330.01.41445351380690480.00.75924104272334587726031174303976178.81224541.4144522419.4785480.143594G12_Y3_016_399521.0395220.05551719.47109819.12314618.951698478.532972335706738.39955933658423G123.327437111473794e-05
396278179267875779066.6542344316841684.088389664666 .. 61.2464990143421840.152321470436103174.3124557370809455e-0507895GALAXY10897.63689505704680178.778761502887051.35793519536994720.00.02015.511529215046068469761TGT0.65302926070741825169010338283505SER23.86339861.066452110.03711148.5420436.6615113.28902803892244808407867136G219.97027419.44065517.933780.01788p0120.018761592107.23865569.753731.22924920.35084627.433583319.02254534.2771847.433583319.02254534.2771846.01.4389416-0.065962880.046224613S1021009115292150460684697613107400000178.778761502887051.35793519536994721427.103811-0.0050.0050.0110.011178.778740706669540.01.35789305670389050.00.72617835272299587726031174238368178.778774021.3579348318.2187180.152314G12_Y1_CX1_018621.14812519.31712318.24160817.860417.571278510.5420495081566538.54015758732808G121.2518930471032511e-05
396278239665855768618.5778959281741099.9542151303162 .. -13.0642142969362160.18031044959984415.011526325310797e-0507896GALAXY10614.13190765399490178.79257399689821.5050951837930340.00.02015.511529215046068469761TGT0.97899693624197715169010339723744DEV9.83287427.4246452.593864407.42358103.85631633.0994503892248454836491136G220.26997819.81333218.3616180.01788p0150.01841319748.88733.587352.0174820.52709995.234691614.59995627.9992035.234699714.5999928.4252434.00.470173270.0985848460.19844648S1021009115292150460684697613107400000178.79257399689821.5050951837930341427.103811-0.0060.0060.0090.009178.792549575499660.01.50506108199878640.00.766291713629587726013994631360178.792573741.5050575119.0323770.180174G12_Y2_042_368422.35725420.11621319.0074718.61402918.249283614.939105238334738.944160557996156G123.767466831602611e-05
396278239665864818604.0819264054395.82173609472824 .. -17.9564088604483520.14966775380682499.765666841871356e-0607888GALAXY10553.70037919282910178.837819452373741.5156753600167670.00.02015.511529215046068469761TGT0.734795739970341251690103397231649SER16.26198836.4023370.93679165.8285245.2681611.10314000.00.00.00.01788p0150.01939417100.9450978.485071.34554590.375973022.89369946.477522412.6226722.89369946.477522412.6226721.10357052.84445450.6107672-0.3437365S1021009115292150460684697613107400000178.837819452373741.5156753600167671427.103811-0.0020.0020.0150.015178.83780986261240.01.51561957869376340.00.7031302219743587726013994696990178.837838351.5156683218.709240.149894G12_Y1_CX2_007620.7751319.55502718.69822918.25606217.913076501.6266008308940638.50190279349681G122.01601653025351e-05
396278300105783297754.733414856717527.5785436129403 .. -7.7482408116727730.213979353661882785.8908448892770265e-0507890GALAXY10192.81308423727750179.064467520175131.6427545670253540.00.02015.511529215046068469761TGT0.731175760217351551690103411641433DEV4.55417414.53792928.321898197.991745.99648320.6546000.00.00.00.01791p0170.02169811728.21156720.6439572.10500140.50590792.52180228.05014915.6828062.52180228.05014915.6828064.00.50519973-0.03364423-0.4733275S102000911529215046068469766553700000179.064467520175131.6427545670253541427.103811-0.010.010.0050.005179.06442794022250.01.64273543510576480.00.7630033219762587726013994762281179.06448221.6427632719.7437460.213584G12_Y2_014_372522.12273820.90774519.69423719.23598718.883875744.308290212011339.3587642812638G121.706053700392514e-05
" + ], + "text/plain": [ + "\n", + " TARGETID CHI2 ... FIELD GAMA_SEP \n", + " ... deg \n", + " int64 float64 ... bytes4 float64 \n", + "----------------- ------------------ ... ------ ----------------------\n", + "39627751493206196 8463.270686775446 ... G12 4.881250127241795e-05\n", + "39627751489013608 9649.093964904547 ... G12 1.9899652496315195e-05\n", + "39627751489013565 8576.984808493406 ... G12 2.0754840503260076e-06\n", + "39627751489013325 8832.450295948423 ... G12 8.880441323715524e-06\n", + "39627751489013290 8921.76774309948 ... G12 1.4392601038633543e-05\n", + "39627751489015473 8628.5112837255 ... G12 3.6261842145141034e-05\n", + "39627751484822207 8712.22295193933 ... G12 1.7375790599941835e-05\n", + "39627751489013618 18748.808570712805 ... G12 2.6865550455542155e-05\n", + "39627751489012431 9577.251358687878 ... G12 6.669326376680688e-06\n", + " ... ... ... ... ...\n", + "39627823966586623 8774.540522575378 ... G12 3.447016960091287e-05\n", + "39627823966585927 8633.59043803811 ... G12 1.3541281718852934e-05\n", + "39627817922594883 8570.026254475117 ... G12 6.699043876701257e-05\n", + "39627823966586264 8402.992509886855 ... G12 2.2802162403838536e-05\n", + "39627823966586990 9520.32962005958 ... G12 2.076941215168854e-05\n", + "39627823966585969 8679.489034175873 ... G12 3.327437111473794e-05\n", + "39627817926787577 9066.654234431684 ... G12 1.2518930471032511e-05\n", + "39627823966585576 8618.577895928174 ... G12 3.767466831602611e-05\n", + "39627823966586481 8604.08192640543 ... G12 2.01601653025351e-05\n", + "39627830010578329 7754.733414856717 ... G12 1.706053700392514e-05" + ] + }, + "execution_count": 55, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "match = desi[good_match]\n", + "match" + ] + }, + { + "cell_type": "code", + "execution_count": 69, + "id": "10213f47-7c4b-4f45-a195-59c9be52cd0a", + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "Text(0, 0.5, 'Declination [deg.]')" + ] + }, + "execution_count": 69, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAmAAAAE9CAYAAACsk95kAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+j8jraAAAgAElEQVR4nO3deZxkVXk38N/T1d3V66w9G8wm+6AwjTSDsrjwilFUwAUXYhRcEKN5JcSFRN8Y1CQYE6PvG01CRNBEXFAQFIlAVFZh6IGeYRkQGGaYAWaf6Zleq7vqef+oWzN3bt/l3Fv31r1d9ft+Pv3p6rrnnnu6q6vqqbM8R1QVRERERFQ7TWk3gIiIiKjRMAAjIiIiqjEGYEREREQ1xgCMiIiIqMYYgBERERHVGAMwIiIiohprTrsBYfT09Ojy5cvTbgYRERFRoDVr1uxU1Xlux6ZVALZ8+XL09/en3QwiIiKiQCKyyesYhyCJiIiIaowBGBEREVGNMQAjIiIiqjEGYEREREQ1xgCMiIiIqMYYgBERERHVGAMwIiIiohpjAEZERERUY9MqEWuWrdm0Bz97eAsEwDteuRgnL5uddpOIiIgooxiAxWDNpj1417/eD7V+/tFDz+MnHzuNQRgRERG54hBkldZs2oML/u1g8AUAxRLw0e89lFqbiIiIKNsYgFVhzaY9eOe/3o+STj22e2QC5//LvbVvFBEREWVe6gGYiORE5BER+WXabQnr3+561vf4wJZBrNm0p0atISIiouki9QAMwKcArE+7EVH87qntgWW+etu0/NWIiIgoQakGYCKyGMBbAHwnzXZENVF0GXt0WPM8e8CIiIjoUGn3gH0DwGcBlFJuR2KKdfubERERUVSpBWAi8lYA21V1TUC5S0SkX0T6d+zYUaPWxeuyHz2SdhOIiIgoQ9LsATsdwLkishHAjwCcJSL/5Sykqlerap+q9s2bN6/WbXS1ZtMefP6mR43L3zzwYoKtISIioukmtQBMVf9SVRer6nIA7wXwG1V9f1rtMVXJ+/WDB583PkcBXPUrTsYnIiKisrTngE07n/zBGte8X0GuvmdD/I0hIiKiaSkTAZiq/k5V35p2O4Ks2bQHL+0bj3RuScGcYERERAQgIwHYdFFtTq8bH94SU0uIiIhoOmMAFsIjm/dWdf7D7AEjIiIiMAALxSTxqp/nd4/E1BIiIiKazhiA1dDYRDHtJhAREVEGMAAzFMcE+io70IiIiKhOMAAz9IUQiVf9XB8ifxgRERHVJwZghp7ZMRRLPd+9l/nAiIiIGh0DMEPVTsCv2DkULY8YERER1Q8GYDU2ND6ZdhOIiIgoZQzAamyylHYLiIiIKG0MwIiIiIhqjAEYERERUY0xADMQ9ybaV/2quj0liYiIaHpjAGbggQ27Yq3vJ2u4KTcREVEjYwBmYP/oRKz1FZUz8YmIiBoZAzADP3s43h6rlib+2YmIiBoZIwEDO4cKsdYnsdZGRERE0w0DMANx76FdKHIIkoiIqJExAEvBZEzbGhEREdH0xAAsBRPsASMiImpoDMACXP/g87HXWWAPGBERUUNjABbgW797JpF6407uSkRERNMHA7AAe4bjXQFZ8dXbmA2fiIioUTEAC6CazHDhM9uHEqmXiIiIso8BWICxiWQmzHe3tSRSLxEREWUfA7AASU2Xb8kxHSsREVGjSi0AE5E2EVktImtF5HERuTKttqShJcfYl4iIqFE1p3jtcQBnqeqQiLQAuFdEblPVByLXuHk1sPEeYPmZwJJVsTUUAFa034ueznXYOXwi1o+eYVy+ONmJXPMwOkYW4MjxZjxQWoGH9RjsG5+MtX1ERETTXtT38QTf/5OqP7UATMuz2ysz0Vusr+gjfptXA987F5gcBSDAzCXA0Fag6FjFKDlg9jKgVAT2bwWK40AuD8xeDrR2ATueBHKtwILjgXnH4mfNv8Nv5+7G92d1YYsI0LMJ3XoLOkollEQwAcGyQhPaJ1vxYtsoVIDOomJLaw6HbuG9CVsmJnDljptx+GgeM4ZHgS9bh0rF8q+eawO6eoDiJDCyq9xWlAAtlf8yTU2AarnNFU0tQGmi/Dt3LQC65gHFifLvPbQNkGag72Lg7IbqYCQioiT0Xwf87u+AoR3l96SZS4EZC4HJ8XJw8sQtwODzQEuntfFxE1CaPPheWyoCkyOH1tnUCsxcDEyMlN+3pxCgyQpXShMGjRQg1wIs7gNOeA8wugvY/iSw6V5gdLB8HWkCtHiwPASAbc63NJd/v2VnAJ09wKM/OXgs11p+Xy45OlI+fEeo4EySWuVndHGRHIA1AI4C8C1V/ZxLmUsAXAIAS5cuPXnTpk3uld3zT8D/fCn2Nt7Q1Ykv9cypNKYcAJmqlHec958vbsPK8QKkltPATr+MQRgREUXXfx3wy0+l3YpscwRhIrJGVfvciqY6EUlVi6raC2AxgFUi8gqXMlerap+q9s2bN8+7sufuTaSNd3R2lG/YgynnV+W482e38wD87wU9tQ2+gEOjdyIiorDu/XraLci++75pXDQTM8FVdS+A3wF4U+RKNkefOuZHKj1Xfj1fzmjK/nMl+LLZk8vhhq7OmFpoqGNuba9HRET1ZWxv2i3Ivl3mu+ekuQpynojMsm63A3gDgCcjV9icj6llBw3kW3F/R3v5B3vvlimX4Kvy843dXTG0MITCcG2vR0RE9WXW0rRbkH0hOjvS7AFbBOC3IrIOwEMA7lDVX0aubcbiuNp1QH9bW/lGUODl7B2r/Oxz3vxijVdB7t9W2+sREVF9YQAWrH22cdE0V0GuA3BSbBU2t8ZWVUXf2BiaMBMle4DlFlT5DUHa2eq5eHB/DC0MoXtBba9HRET1pYvvI4G6fOaqO2RiDlgslp8Ze5W94wW8bsRaLhtmCDJCjxgREVGmjXMP40ArLzQuWj8B2PhgItW+0OzRSajqPTHfma7CHnhZt2+p9ST8vZtrez0iIqovL/Sn3YLs2/aEcdH6CcAQf+/SQL4VT7V6DG1WesT8gjDA8/g6r3qTwt43IiKqxuGu6azI7sFvGxetnwBs4crYqzzQS+U3pysg0PI6d28uV2XrQioWyrsFEBERRbHvhbRbkH2FUeOi9ROAbR2IvcrAPiOXoUUA7sGYI5/YkYXC1DJJ23hP7a9JRET1Yc9zabcg+5rM1zbWTwCWwBDk24as3Flhth9yy/1lv8/6vr/WPWCQRBYqEBFRgzjh3Wm3IPvG9xkXrZ8ALIEhyN7xApZO2Db+9AvEKpPy/YIvmy3NNQ7Alr06mR3iiYioMZx9JSC17jyYZuYdY1y0fgKwBIYgAeCiSr4uj0AqkMc5e5pyGMjXciI+J+ETERElKkRS+PoJwIZ2xF6lAvhxZcugoODLb1WkM2WFCCC2TPu1MLyzdtciIqL6s3k1oMW0W5FtT/zcuGhqmfBjFyL7rCkB8FxrS8iTDDLlW8FY39hYtIZFkcBOAURE1EC4kCtY0XyBXf30gC3sjb3KgXwrDvlT+iVfDSuuekzt3VLb6xERUX1pN99oumE1mXd21E8A9uiPY6+yv63t0C2IvLYj8gqmvAI2q56aDkGO7WEeMCIiim50V9otyL6FJxgXrZ8AbOvjsVc5s2iNdUftrXIGbI5cYAfqrxV2HxMRUVTsAQu28ynjovUTgLXE35u0vrJK0WQCfpDKKkrbasrBWucCYx4wIiKKij1gwQrDxkXrJwBbHP8eVbtytj+PSS+YXxl78GUNTdZ0Ev7MpcwDRkRE0Y2ZJxltXOYjZvUTgJ1+WexVzi2WDv4QZu6XQbnDJovoHa/hdkRD22p3LSIiqj9/uC3tFmRf92HGResnAEvAuUPD5fSl9gDKPrHemffLayjSPonfKv9HwyMh4uQYFCeCyxAREXmZqOGozXT17u8ZF62fAGzt9bFX2TtewPzJyUPvtE+sN8mObw/YbOXXteVrnJu+FFyEiIjIy0zzLO+NKdy7ev0EYAmEMzd0dWJbsyNXrbMHrHKfZ7McaSws61pbatsDBgD919X6ikREVC/mHZt2CzJOgbU/NC5dPwHYyvfFvknoTW7bELnlAouwR+SkSO13Z1x/c62vSERE9SI/I+0WZN/QduOi9ROALVkFnPZnsVY56Mzh5fxycuT5mnK/Ta7WmfABYOGJtb8mERHVh63r0m5B9u3dZFy0fgIwAHjillirG2+y/XkqPV/2L2cQ5TLU6GWyqQkD+Rrvz9jGTy9ERBRRCv0G087ezcZF6ycAu+OLwJ4NsVZ5ZMFKE+FcBVlhmprCPlfMNhH/F12dMbXUELMYExFRVHs3pt2C7AvR0VE/AdijP4m9yq2VCfjOOWB+Qhyv+YcJZjEmIqKoVpybdguyb85RxkXrJwBrjn8rojG3fRyr4Vg1ee6Q+ZYFseBWREREFNXZVwLSHFyuke180rhoagGYiCwRkd+KyHoReVxEPlVVhYMvxNSyg44r2DLVm+T7CqGlVKptJnyAWxEREVF0m1cDOhlcrpF1mE/1SbMHbBLAX6jqCgCvAvAJETk+cm3F+IOZi/ftL98wCa5CpqJoqfUqyKaW2l6PiIjqy33fSLsF2bf4FOOiqfUlqupLAF6ybu8XkfUADgfwRKQKe44Bdj4VXwNRzoS/YGIS21oi/JkCsuTnq2hXJDHnSKMGs3k1cOcXgZfWHfywI01AczvQvQA49eNA30XlchvvKQ93//bvgefvB9AETI4CLR1AvhsY3QuUJsrPkVxL+cNBvhNo7QSGdwEtbeXbQ9vLW5/kcsCs5cCe58pbanUvKm/3se2Jcm67hSeWJ75ufxJ44ufl9rV0AEe+/uAesRvvKR/feC8wZzmweFV5Sf2K84AFxwO3Xg7s2QgcfgowZxkAARauBJ65A9h0HzA5DsxcUm734X3AvhfKrzc9xwJv+Bv2LlNj2L817RZkX4hcaaJp5KNyNkJkOYC7AbxCVT23W+/r69P+/n73g5tXA9ecHWu7BvKtuHjRAkzpcI2QePUQqsipYmDTlurqCevDd/CNgsLbvBr47psALabdkowS4MO387lF9e+OL7IXLMispcBljx74UUTWqGqfW9HUJ+GLSBeAnwG4zC34EpFLRKRfRPp37NjhXVGI9P+m+tvaysFXUP4vN36JWgEURWqfB2zjPbW9HtWHjfcw+PKlwH3fTLsRRMnb9UzaLci+wohx0VQDMBFpQTn4+oGq3uhWRlWvVtU+Ve2bN2+ed2Uh0v+b6hsbqzTi0AMm+b/8esmsY/1t8a/c9MVVkBQF/2+C7Xo67RYQJW8n/88DnfR+46JproIUANcAWK+qX0+rHX56xwuYXfT45G/fjihgvtcBjjQUBwK8mshxiISiWbIKaO1OuxXZ1tGTdguIktdzdNotyL6zrzQummYP2OkA/gTAWSIyYH2dE726ZOayeU5d90vOajBE2VMs1jYNxZxltbsW1R9JfbYCEaXtqHjnWdelzauNi6a5CvJeAFXOZrfpWhBbVXa+M1+8er287rf1lM316llLyt7N5X8M9oJRWJtXA+ODabci24Z3pt0CouQ9+uO0W5B9a683fp/1/FgrInMMvmbF1uhqLVyZSLWnjXjMAzNhH6YEDgnMduVqnBaiNMFJ+BQNJ5gHK02k3QKi5G2LliWqsZj3K/n1gL1offnVlgOw1PhqSXrmjmSqbbUSmIZNPREwLyyVNWWcTE1RcOJtsMEap5QhSkOu5hksp5+V7zMu6jexY72qHqGqL/P6ApCd3Z33v5RItZuiJGEF/IchAZy/v8b7QKKJw48UDSfeBiuyB4wawBLzLO8NK0QvoV8A9mqD803K1MacIxOpdmapFHud+WIRl++t9ZyaUqjJgUQHxDDx9oauTrx+8WHoXbYYfUsOw0UL51edB28g34pPzZ+L9y1agBu6OqtuY1W4SIEawenVbdncENbfbFzUs3tHVQNzJJiUqZmRZCbBdpYirq70GYIcz+UwkG+t/WbcISYHEh1Q5cTbK+bOwa3dBwOkYi6HNbkcLlq0ANe9tC3S88BZ52P5Obizox3/vj2lyfAa/wc1osxZsqq8rR0TM3sLETJE+tgmIr+Mcl6iFp6YSLUtUdNbBKyQvHZmCnmVNtxV+2vS9LdnY+RTb+jqPBgoOXaUKCJaMuKPze+ZWieA+zvaccXcOZHbWp34FnQTZdbm1dEWpDWSFx82Lhq13/yjEc9LTpv5BphhTCT0wrojl0IGkMHNtb8mTX/Lzoh86n/NsD5oeOTKe6Y53PPghq5O3N/RPrVO6/at3Z213+ILADrSCvyIamjtDwGwt9dXiPmgkQIwVU1mxns1Elrht3zS4I8Z5hNBJQv+aAqjt01cwUIR5KPPr9qW83iJiRgw3dTddcj5bnXekvZ8MKJ6lcCWf3VnzhHGRQM/forIo5g6qjkIoB/AV1Q1GyshE8pPcvHgfvymo8M/rUSYFBXWZt7DXm9MSTquio0GqHEN7Yh86kST7f/c+Ryyngu3dHUazwMzKfVgPoUPGt2Lan9Noprj8GOgCfMMByb9/7ehnLbqeuvn96I84WEQwHUA3haudQl55PuJVNs7XkBLqXToG0kMUvk3nn9cGlel6a5rXuRTO0slFCpJh62Ay7kn6s+6unDu0LBREDbRFPxh58XmGic5BhJLg0OUKQntOFNXBl8wLmoSVZyuqn+pqo9aX58H8FpV/SqA5RGbGL/m8JN5TZ1YsA1DVjsB0Tr/3KFa5wEDE7FSNPmZkU+ddM6htAdf1gT6kpgvSmkxeP5NNjXVfjJ+QquwiTJl5fvABScB2s03CDIJwLpE5NTKDyKyCoA1EQOT4VqWoLF9iVV92Z695RuVNw2/NwGDN4jZtd6Iu4IpKCiKCFtYVZ4FJa++Xsew/TqDYcOBfCs22yftuz3XbHPLvj4reuAYHvOAUQNYsoo574K87q+Mi5r8JT8C4Dsi8pyIPAfgOwA+KiKdAP4+WgsTMJLcVLTe8YLnno5T+B2z6jiwv2Qt5dprf02qD90LQxVXAJfO78HJyxZj2HB+5E4rN56XG7o6cdHC+RixTwXwWFlZuf9XtZyMn08hrQxRGpjzzlsuD/RdZFw8MABT1YdU9QQAvQBOUtUTVXW1qg6r6k+itzRmpeQSw93Q1Rl+L0g3Vh0dacwAE06epIhOvyxU8Uvn9+D+jnYURIBKwOTXM2w9L74yZ7br4YF8K77SMwdFW84vz94v2/35BHax8DQxUrtrEaVl82pwIr6P4jjQf51x8cAATEQWiMg1AH6kqntF5HgR+XAVTUzG8LbEqr6js6N8I2j40VAqI+gm6TSI3IRcYdzfbs3HdCReDfJUa4trL9gtXZ3lzEOVOvzqs90/WMuVxnF8QCPKugjTERrOg/9qXNTkFeo6AL8GcJj18x8AhPtIPM2dPWx9uvVLRWFCFU0A3pbGBPxcCivDqD6EXGHcWXT0PKkGf3Cxgiq3zPgbWnwWa/vUO9jkP6wZqwR74Ikygwu5go3sMS5qEoD1WEONJQBQ1UmU01JkS4ITA4+e8Ok9Mu0Rs8p9cO++dCbgc/sIiirkCuMx54cUkx6wSoLisUPnRw7kW/Gw33ZFfj1hUsOkrNwbj4gAoMU8D6FJ1DIsInNhDfyKyKtQzgGWLZLc1j4HPpW7zT8x7RGzys1IKxAqcQiSIpp3jHHRgXwrRmMc+rt2RnflhSfS+RtCbnVERD44BBls5mLjoiavlJcDuAXAkSJyH4DvA/izaC1LUIJTMA58KnessjogKKiyhmByqulsQURUjYW9xkUP9DiZBkyO59S1Mw5dTfhgW8is9o7n4p40krIS1av2uWm3IPvmmSc8D/x4qKoPi8hrARyLcpjzlKpmsDsluSHI3vECOopFjHjNowp6s7Em7/eOj6O3kMLwI8DcLRTdqFmKl4F8K35W2asxaL6kLRGr3caWlkN+Hq5yB4qWkkLB1JFEsTB8LWho40PGRT0DMBF5h8ehY0QEqnpj2HYlKuHcJD2lEp53BmAmk/JtZZ5xvLnUFAMwispw4m1/W9uhqxX92MvYniMmme49uZy7tTnH4IsoLuwBC/ZCv3FRvx6wyh6P8wGcBuA31s+vB/A7ANkKwErJ9ixN2VLFdEXkIcviy6uyOAmf6lHf2BgEM6FhVwvbyj5tpaLoHS+Ez7/nUnZfzHu4empmomNqAPd9M+0WZN+Kc42Ler46qerFqnoxypPvj1fVd6rqOwG8vPoW1oGwk4Kt8m7L7GuDARhFZDjxtne8gEUTtt3J3IJ+rw8CVtBWEsEvrHlk35k5o3ysitQvKlKbLYlWvC24DNF0t29L2i3IvqduMy5q8vFwuaq+ZPt5GwDzZVHT2EC+FZ+aPxcXLlqA7mqyatvyIO0TwXdmzogtP9FAvhVfnjsbX5o7279OfkKnqELk/mkOWiFskEB1p7WKcm/Y3itncGfV+dPKvLQk5Wu47RFRWma/LO0WZN/OPxgXNVmj/TsR+TWAH6LcjfJeAL+N1rLpYyDfiosXLsCk8/0iSjJWW/qKa2eVP9ULZuL1IyO4eHB/5CHJG7o68aWeOQd+/ml3F77/0jaP+tgDRhEtWVXe46w47lvs67Nm4vlWa55jFUmLB635ih1QTNngx69ej/uLNfnf50wzIkKozg6TvSA/CeDfAKxEeT/Iq1U1e2koYnbtjO5y8BW0/5wpR8+AAvhNRwc+vHB+pN6wG7o68aW5sw/UV6nzT+f3uJ/Q0hH6GkQHLDs9sMiv7Ckoqhg2fMYK4nIlj/0eQ+p2qyduuzcmfw2itI1lLwVo5sxfYVzUqI9fVW9S1T+3vm6K3DAHEfmuiGwXkcfiqjMON3R14jeV/R+jJF11cqvDepMqeGy/4mcg31ru+XJ5o9ufy+FjbkFYaXLqfUQmNq8GNvwmsNisYohs8M7tiWy3h6yhx+EoCV0r9djq29Zcgy2JmKCSGkH7rLRbkH27zIcgPV/hROSXQSeblAlwHYA3VVlHWS5kwkYfh2y+7bYBsLMnzLRnzCOB6y86w83PmpKZ31H//R3tU99wxvaG2qWd6ADD4GJ7mKSnPr1kzVaP1YxKQBem59n+fLVd5xdJb0nEnSaoEXQfFlym0XXOMy7qNwfsDBG5xee4ADje+EouVPVuEVleTR0HtM8ChrbFUtXZwyP4fXube+b7ajbjdrIStG5obcUNXZ24wHCT7pn2Nya3IEwV35g9E9dt3XHosfU3A30XVd9uaixj+wKLDORbsSfMhu+V/13XhKzl5922Sn1hsup7lN0Z4/ZIRA1r55NptyD7QgxB+gVg5xmcn3hCKxG5BMAlALB06VLvgsV4mqIALhgaLs+v8vqU7nyhryYos96EburuMg7ABg3emNa1ugy5dHjMDyPys3VdYJEDWwjZgyov9uMu5cZzOXx91kwUw66C9L1muKqIyEVzWmmUppFt642LegZgqnpXLI2pkqpeDeBqAOjr6/N+Ga0mTYQLQQyv2SFWgm0MMXyzzz4U6lH/hMjUpK+7nzW+BtEBK84DnvWfA/akW8DvxS9Is479uDvknpKB14ynGu/62cNGDaB9dtotyL49zxkXrZ9Xjdk+vWMhCMqT8DWOF/4QdezP5cqZvw081O4xB8x+XbfJ/fz0QglpQ8jFKgEri0dNe7/85ofZjs0tJrtVGeYenWz9RFlw0gfSbsE0YN51Uz8B2KxlVVdxQ1cnLly0AF+u5NaqNgjzWel1COs6/zQ7OGP3DV2d2BDUW2ZL+nqIeQ2RP5fitv7mwCILKxnwTSfMB5Qzfgnze45WetoAnGs4vB/ZzMOTrZ8oCxZUNe27MYRYEJhqACYiPwTwewDHisgWEflw5Mr2bqqqLV+fNRNf6pmDR/Ot5Rd/k022gzjnkAXUORzwqf9j83vwpZ45GAma7Gxd51pnxv2Fvf7nEblZeGJgkUfz1ouOac48rzlg1nmhdp5wftBx1AUAP+pKOBu+wd+IaNpjupVgIVJ1BAZgInK6iNwhIn8QkQ0i8pyIbKiqgRZVfZ+qLlLVFlVdrKrXRK5sl/m4q9NAvhXXOvedC9izzv7p2rNchVfg5bJ1ynmHLXAtesXcObi/o92/PkddEODztkz5GN0VfB6R03hw8sUh57+k6f+8k3VeqCkAfolfrft/1d2ZbC6wthnJ1U2UFe1z025B9oWY6mPSA3YNgK8DOAPAKQD6rO/ZItGnzN/S1Vme/GXyphHmjSHMvnjW7Q2trVPeKAbyrbg14qTk51taDtYXYk8/ooP8/+cG8q1Qe++tX5Bl0rOsitG45mBaH5gUCJ3wOBQ+t6gRbF2bdguyb+EJxkVNArBBVb1NVber6q7KV/TWJSQfPH/Ky66cy5uHvYfL/mbgkml7ioDhEE/Wdb7hmAt2S5cj+ArZ83btzG7g9MvKe/oRhbXyfYB4D3s/5JUY2JRLT/AMkyFI0+eUVa5vbCxC4wx0H8bnFjUI5nMJNPco46ImAdhvReRrIvJqEXll5St66xIyMWXbXmM9zhVSQcvk7d/duJUJ8ea0oaXlkJ93OZNIhtxrb3uuGdj1jHF5oinUe5uhUyqBjdsHFhMuz5PRA9etojetVva/WN6uiajerbwQyed0meZCzJPzS8Racar1vc92nwI4K0STkldFioXjxgtANw7Nzm3CJOlkBK2O67+QM3mYbG2qsOpZVpgA9r8UR9OoEa293vdw73gBUixCK4tDvJ4TIZ4rY0HJhsM876znaX9b26F58eK09nr2glH9W7IKaO0CCvvTbkl2dS8yLhr4zq6qr6+qMbXS0hH51CmZ5cNsfWL/HiZ48zHiWA0ZmHbCi9We9flWYM6RVbeLGtTQjsAiTU1NONBH5vXBxPQDi8kClwhzMRMbggTAXgFqGFWMNjWE0z9lXNRkFeRMEfm6iPRbX/8kItEnXCVlJNq0NIVtM+ywwZPbnCyvT/4hOP+9Q2/J4vBcawt+ueHhquog8nJDVyeKHhvDB/KbLxmwsjFM3aKaXO8XAOzemFzdRFSXTN7ZvwtgP4B3W1/7AFybZKMiaY8WE/60qxMbKtuouE22N9BTKpX/kPZ8RM55LSEm7Rebmg5kxR/It6JkmlvJrV5rBdj/5BPftpPqlv//3U3dVo4tt7PYBDMAACAASURBVP9Tt9umOcKi9CZ7fBhSa2uuxLz0SHJ1E2UJd1Txd983jYuaBGBHquoXVXWD9XUlgCMiNy4pc8xXHtj910zbJsJ2IXqyOoolfO+lbXj3/iGc7JyQ7FefG6vcjdab2i/sKyCrmNzc3zkZ7lyiii733HQVg00+zx3nbeeHnKCgze1nPz4LZxJNQ7HopOTqJsqK/uuAiYR3lZjuXlpnXNQkABsVkTMqP4jI6bAtUsqMF9ZEOm1fmOE9jx6y/bkm9I4X8H927cGOoMnDdj5DLRPW951VDj9WDMZUDzWghSt9D4+H3Yjab3WxPSgL+4HDq4fNa2uuOL3sjOAyRNOdwbZkDS/EAgWTV86PA/iWiGwUkU0A/gXApRGblhyfZfL+54Up6z68Ye9b2tIcccWiw4R1rCfMliw+mmJYHEANauuA7+H5kyH3gawISk7sd5+TW+Dm+IBzYBP7uOXyTMRKjYFbbgUL8ToYGICp6oCqrgRwIoATVPUkVc1eOtxF/p/SPYX5UOzxRrA/lzswZyuOVZAAMGwN67ytsolwlfU2KbBm055qm0UNyftJMpBvxWNtjn0g02AQuM0vJjQMX5wILkNUD7jlVrBZS4yLegZgIvJ+6/vlInI5gI8A+Ijt52w54d3Rzqs2XrJe5L9j7SUZeaDPsaHwHmso82lHUlbPjYf96gUwkWvCP9/wjaito0a28n2eh/rb2sw2r0+DY0jy4r1J5S4qBeZKI6oL3Asy2KxlxkX94gWrSwfdLl9dUduWmAgbTQ/kW7HLLcdW2G2GAOyvdo6VY4JyuzX0eGdnx8HjLuVc2+dcjWmVF/RX10Yihz63RSd+G9mb3J/QcPlvKpvZJ8EgVxrRtBfhfZa8eU5YUtV/t27eqar32Y9ZE/GzZfuToU+5pavz4Kd3k+SO9kz5jrIdWsJAvhWTlUAsTLJIl9VgQ7kcBvKtmD1ZDFefz+rNXSMcv6cIfLbW6B0voKVUwoT9A0jY/F1uKyPjCMLsz1VV/Ly7E5fvHay+Xjd7NiVTL1GWcK5jsO1PGBc16bb5f4b3pWvTvdWd77YKy6uMy9L5vU1NB5e5h0kZ4ZExvIjyJtq3djs24o5ISoonR7hSiyIIeNGdskzENN+Xm4i5+EzqSzQTXsRE0ERUZ4bNe8M9e8BE5NUATgMwzzHnawaAiHvjJCjCa/UK5x6QQLhAx1Z2vKkJD1USPYbdp87j5yfdEsRG1ASXN0qiKg3kW6fu1BCmB8vkA09MJpOcpzbvuOTqJsqKEBtNN6zDTzEu6tcD1oryXK9mHDr/ax+Ad1XRvGRE2J9qfd4jwIm4nH5d0GqwMG9KYc8PyDhebBKcNut7ZtcnsvN50T2k19fOrxe42udBxLLjTU3JZcMPMexANG1xEn6wOcuNi/rNAbsLwF0icp2qZn+CQzH8AMOunC3+jNIL5piv0mK6xYrzem7lVNFc0kPL+rXLLfu4/bsqJmb8wb99RG58hiDLk/BnVtfrG7Wc8zlr8Jz6RVdnMntCju6Ov06irOEkfAPmHxxNsoaOiMjXALwcwIFMhqp6VviGZUvV2eGdL/5RzvWxsxIghnkj8joOoGXfMaatIzpoySrPQ73jBXQUixjJhZyV4PWhIkogZ/IBxRLXzhIujUmoXqIMGduXdguyb2GvcVGTV6MfAHgSwMsAXAlgI4CHorQrUce9NfQph2yhYhLkhB06cSsfIlCbCPNGFFRvqYT7937Q+NpEB2xe7XloIN9qHnx57CRxCPv/sulzxdmz7HPeYI5bchFFxjlgwQJ2DrEzeTWaq6rXAJhQ1btU9UMAXhW1bYmZcVjoU96xf6h8w+SF3v4J2yNv0YTb0J9bPUGirCJzBmEem4EzGz6FtvaHnoc854DZ+ex56vo/HnUvyMo5Ph9I9oTtqTO/cEL1EmVI98K0WzANmL8WmARglX02XhKRt4jISQAWR2lWoh75r1DFFcDRExM4bMJjexK/3qvKC7yjjOeeeGF7zqw3nommkG8WbvPAHP79rmfD1UnkM6dhZtGWp+5Acffg37v6mBOv+gxjLp9IaNug4lgy9RJlyemXpd2C7Ftovi2iSQD2FRGZCeAvAHwawHcA/Hm0liWoEG4V5Np8Ky5eOB8vtljT4JxvIEG9V/YeMev7C5WNuH1SS/iaUr/Zab5sQd2K9nvxyPPsAaOQVl7oeWiw0qMUNcVDmJx5Yep0sp4HZ4wkFCgJhzapAWzjat9AIRYqBE7CV9VfWjcHAbw+YpOSNxkuALt2Znc5L5Db6qmwbwjWi/u4yQRf00nGUXoFPJK6Vo71dK7Dhr2vCV8vNbYlq4Bcq+tKY9dVkCZJjSvHK6J+aDFlPc/X51uBygb3cZq1PP46ibJm/c1ptyD7QqTqCAzARGQegI8CWG4vb80Fm7a252y/uulWRG6sN5GcKorOeWJR54BFLR+Qd2nn8IkoTDIdK4XUf51nmpfe8QJmF4vuc6uC/nejrO6N+jy1rEsqD1i+M7gM0XS38ETg2d+k3YpsC9EDZtJvfjOAmQDuBHCr7atqIvImEXlKRJ4RkSuqqy3cEMCygjUXxCTpqV9vlPVmUHRO0A+zpUqEzb8DOc7rLBaxfvSMKBsGUKN75Pu+h48ouMyritqDa7/ttuglYu90xd6w8ypN7d6QTL1EWdI2I+0WZF+I/TJN8oB1qOrnorfGnYjkAHwLwNkAtgB4SERuUdVog8wiobYjerDdsXrLbwgkaBil2l4vtwStbr1yJkM6Hr/PSGL5j6juBSQ5frS1ZeqdUXquqpkG4MYliDsnieFHAFwFSQ2Bm3HHyuRd+Zcick4C114F4BlV3aCqBQA/AnBe5Nq0GKr4vqYQ2e7t352cn9KrXdHllerCfsyNIyv/lMOVKqprHTWinPew3Q1dnSiYpHYwyVXn7AHzS18R0Vmjo7HVdYiu+cnUS0TTS4hcaSYB2KdQDsJGRWSfiOwXkTjS4R4OYLPt5y3WfTWRj5IqwmAo0mj/u7iX3Ye4PgMwCu2kD3geurOzo3zD74NDRdg5YVF2mPCqz7p97Yzu6PX5GdqRTL1EWbL2+rRbkH3P3Wtc1GQVZEKvWK6xwJRXWxG5BMAlALB06dLYLj7qN+QYZo85+zle95t+gne7bpWTjg+w6uAUfAptwfEoP12n/o/PnrTlAatVSokgPs+ZjS0uw6VxKOxPpl6iLHn2rrRbkH2bHzAu6hmAichxqvqkiLzS7biqPhyhaXZbACyx/bwYwIsu17kawNUA0NfX5/1xWHKhhiFdX5693kSqeVMJk+bCbS5ZmGvHFawR2W28B14TLJ+vzP8ynZvoJ+r/vZPPubNL4aYqEJHNvhfSbkH25c0XKvj1gP0Fyukn/snlmAKodjPuhwAcLSIvA/ACgPcC8M74GGTG4cDg80ZFB/KtmIg7caJfj1WtgiKf3rqeyRL4GZ0i8clrM29yEghK7eC3QtitXIKO9Nr5ompc5EINoH0WMLQt7VZk2+v+0rioZwCmqh+1vieSfFVVJ0XkkwB+DSAH4Luq+njkCvdtMS7a39ZW7gKL8wU/zgnDYXsLDM5tVvaMUUQ+eW3OGB3Dbzo7wq9ydOPc6ivO4Uer7rcltQpyWfa2xyWK3cr3Afd9I+1WZFvfRcZF/YYg3+F3oqreaN4izzp+BeBX1dZTZv5i7Zq9G4hnqDEu9lWVYYYxPd6AdjczAKOIlp8J5PJAcXzKofWV3i+voCfs/MdqP8gEDIX2jvun1Ihs8apk6iXKkl1Pp92C7Nu8urx7iAG/fvO3+Xy9tcomxq+5zbjoyvECciXHdPRqgqewbxom14o6jOmWfwlAoUlw2qzvmdVBZLdkFfDmf3A9tCvn8xLitjIyKJ1LpUy1q5I9rjGQVCb8reuSqZcoS3Y+k3YLsu/OvzEu6jcEeXEcbamZCfO9IAXlXzzW6bhhgrA4thcKe74qJmb8AdgLXP/g87jw1PhWlFID2DrgevegSYLfML3McT83HImNr53ZjW9uN98qxNjCE+OvkyhrOucCO9NuRMbteMq4aOCrp4j8nYjMsv08W0S+ErFpCQr3aXhu0dEDVu3wo0n+ryBuZas933Z/y75jAAC3PfaSeZ1Em1cD/de6HtoTZmsf0//vOIfyHc/Ljc0JpaHYxZ4BagDzjk27Bdk3Y5FxUZOlO29W1b2VH1R1D4AkMuPX1Elj1nyWoAz2psMhQLTtV/zOr9wX9fqOtty/94MAgLmdCQ3DUH3yyey8fNJlH0gvpotGElwNmVgaip2cG0MNYGX0RAUNY5b56JJJAJYTkXzlBxFpB5D3KT8tPNJu/QpeeYfsQ4qmQx1Ocb2RxPyG9NgLg7HWR3XOJw3FxYNWcpO4d3dIYrcIADOdPd9x6Tk6mXqJsmTJKqD7sLRbkXHm79cmm3H/F4D/EZFrUR7n+xCAaT+be1fQ3JWwQU+UzbNTMjgaoteCyCcNRe94wT/1Q9RUEgnpcS6+icvpn0qmXqKsKfH9w1fXPOOigT1gqvoPAL4CYAWAlwP4snXftJZzvsYnuT9jWud7/E6jBWYDpxCWn+l7+IRKagfn/1s1qSTi/uCSdB6wJ29Npl6irOlemHYLsu1F8xXRpumb1wP4b1X9CwD3iEhS+0NWIdwLdrd6TMKPa/PfMAImz0dma0+TraqRCQZgFEJATpsDKyGTSGwco5ZiMbk8YAM/SKZeoqyZtSztFmTb1rXGRU1WQX4UwE8B/Lt11+EAfh6pYUlqnx2qeMFtzle1Q4Ym+Y7cVLs1i8G1mm1lEhzhoXq0ebXnoRu6Og/uB2milv98jmtNmKTMiIp7TFKj+MPtabcg25rbjYuavCJ9AsDpAPYBgKo+DWB+pIYlaebhoYq3ug2XeE3ENxVXVv0wbTAMGidtRZqyNy2Nsmzt9Z6H7ujsKN8wfe6YpmuJY2eJWs6/7DSf90E0bd3xRaCUUC9yvTjitcZFTQKwcVU98BcXkWaETbpVC2P7QhXfH/dm3EmIY389S5NtiDZ7Dx5l2tAOz0NnD1sJkKud/+WVwmW6dNe+6k/TbgFR8h79SdotyL6CeVJ4kyjkLhH5KwDtInI2gBsA/CJi05IzOXWfOj8TJu8PcW4tFFW1W7JYx5YWDj7U7ACjUHxW9VwwNGz+/x/leZLBVcSuQmzASzRtdXinpCHLc3cZFzUJwK4AsAPAowA+hvLm2V+I1LAkhXhxVwBqEoaYBj+VT+pR54B5MV0N5lfmQJ4zw/JETgt7fQ8b58IPSBQ8rd3xxbRbQJS8xX1ptyD71Hw+qEkaihLKk+7/VFXfpar/oZrBV8sQQYUAUJOBuDDBT9TNs+Ng8HCMyOSB2+3N02D4lbIjYFVPLJm10nz+xGEdh2aoATATfjAx357N851Yyv5GRHYCeBLAUyKyQ0T+OoYmxq8QLr9PV9gYMmjPuqAtjaKIYwsiy46Wgw91e2uI/fuIAj6stIRNblrtsHpEklQSVoBDM9QYtj2RdguyL2e+UZBfV8hlKK9+PEVV56rqHACnAjhdRP68uhYmIOQcsJFqJwk776sMF6b5yd3nTYuL5CmylRfCb+agxPFcchNzECZJpqFYwqEZagAPfDvtFmTf5KhxUb9XpA8AeJ+qPle5Q1U3AHi/dSxb8p3GRQfyrZioNuWEmzAbdyfBuQ2SjT0PWCuHICmMJauAEy7wPFycJsOFiT4zQ2S/Jpq2pslzPV3mrzR+78QtqrpzStWqOwCEyLxYI/NXGBftb2ub2lsVxz+WSQ9YHAlaI5xrD8Cak+wJoPq0+1nvY0l96Ig5DUVzkh+OtjIAowZw6sfTbkH2xTEHDIBftrUMZmIzD1aeabb2II+7x8qkrjBBVZi2BfwuY7kcVrTfCwDYP87NVCmEzauBFx7xPHzk5OTUO6uZE2k/J8bhylySAdh0yCtIVK0Fx6fdguzr7DEu6veqsVJE9rl87QdwQtWNjNuuZ4yL9re3lW/UYs5W2CDKLkzbKr+L2/Wsenpm3Q0AOHJel3m9RGuvh99axy/s3FO+4RY4RXl+RdmX1eA640n2/OYzuD0uUdw23pN2C7IvxLZkzV4HVHV6LZUb3WNcdPHEBLY11+jXs7+ZuN22/2wyfOm3lYvhm93JS8Ptm0mNLqV5HzF/OEp0DtjCE5OsnSgb+q9LuwXZFyIWqZ9+8xBR52V7B8s3ajFhvnIN53wzv2EWr+GbqElZrXp27n0NAOD3G3YF10NUsfJ9vocfsvcoxyGJlC4A8kk+38cHk6ubKCsGN6fdguxrimcO2PTSbJ57o3e8gAUTtnkrcb0wu61C9HpTMspeH3II0jATf56rICmMgNw/p4yNlW8kEeDEtToZwIX7hqqvy8uI+adeoumLqyADLT7FuGj9vBM3hVuY+fIJ2zqCJOaBxVFn2EUCXqs6K3PAOssrtfaPuUyaJvKy/mbfwyvHC8kEX9XO0XR8ILl8b4K9VKPsVaYG4LMvLFkaMgBrNc8DBgCDIboJY2OwafYht2Me0tk5XJ6nsmWveaI4Iqw4z/dw6P9S0z1Tq92twvH8Gci3hqsvlPp5KSXy1Dk/7RZk32M3Ghetn1eNYrjMGE+2eK4/mCqurVNMhh2D7o/a06CK9aNnAABmtIX43Yn6LgJa/Vf5tUf9v0w6b55tgUt/W1u4thHRoYa3p92C7AsxHzSVAExELhCRx0WkJCLx7OExY5Fx0YF8K4ZzLj1gXp/Mq5z8HlqUQC2oSts+eExDQaEFDD0cXwiRW87vf9i5mrfaDyCVuZGq6KvMVUvCBHuVqQGMJziPsl4seLlx0bR6wB4D8A4Ad8dWY4jX5WtnWp/mnS/u1WbGryaPl1s91SaxtN9ty4H0/O6RKK2iRrV5NbB7Q0CheBMaN9luHyLic6oF5cU3iSkxuTE1gDw/vAc64T3GRVMJwFR1vao+FWulO82r25FLaAguruHCaoO/gOv2LpkVvn5qXAbJF/fEOKdSVNFWtNLKxNSzfMREwgGScrt7agBzj0y7Bdm3da1x0czPARORS0SkX0T6d+zY4V0wRB6wvtGAZfNRe57czqtmdVjYNBSV7y7nNduGII9ewKzdFMLyMwOLtMTYA9ZUUuT8pvaHmTdm3f7CroTTROSSnOBPlBVMQxFoh3lnUGIBmIjcKSKPuXz5L6lyUNWrVbVPVfvmzfOZhzL3KOM6Z/itMgyRUf4QboFPNcvoq9mn0uW8SdsQ5P5RDpdQCEtWBU7Cb4lxBDIvwLv2W3NNojwHHDnxcqrJDj8CwOJ4prISZdrwzrRbkH0htkVMLABT1Teo6itcvvyTCkUVYmy6zy9xpFsyVdM3AbdhwDSGId3OFTmwGTcz4VNoAcP27wgKmNyeDx5lc1bOrtnFKob1bM+BtpL3PpaxmbE4+WsQpa3n6LRbkH2FYeOimR+CNLbvJeOivaaJI+PYUDjMuUFtCjoekMqikoiVmfAptIkqVhDae5UNFroMW72186sJwGxm1mLLsU33Jn8NorQddXbaLci+nHlS+LTSULxdRLYAeDWAW0Xk11VXOhKuVycXV26vOAUFa1WuqixOlpPVcg4YhVby3z3hzs6O8g2P3tcwWqz/180xLZapyfT4HHOMUQPgjg+xSmsV5E2qulhV86q6QFX/qOpKQ06CDVyzFXUuWBRegZ5zf8eorN8j1/4CAOAdr+RwCYXl/3/4hmErtUk12eut2yePjZd/junpt6epBi9zHbOTvwZR2gwW5DS8SfP5pvUzFhXy03LBtLep2gDIeX6YLYf8gsAq2nXyMr5ZUEhN/h9wLhgaRkeYIUPn88BlmLJnMp6+q8hZ+sPgGxMRAaF25amfAGx8f3XnV5MF36+upJK9RphbtnPva8JfjwiASaLVfJjqKr27lZXCtt7eSm9aXHO3DvSoJSnEyieiaeu+b6bdguwLmK5hVz8BWIg8YADQ6lwZ5RUYRVnNGDawiiPjfUAOslypdGAvSKLQDF5UjgizHREQ+Dx5u3NlZdhce6oQABfvq/LDmYlN9yd/DaK07Xw67RZkX8dc46L1E4A1h5sE22n/oZqNtE0CobDX8ivvtbrSLQu/7b4ZtRiGoYZ22V5rE9oqU6/c2F1OKXPB0DCOKBQO1mm6Yb3NMYVC8jnAAGAywX0mibKi0zy4aFhLX2VctH4CsBVvC1X8tJGAbPgVca9MDLtKLGjj4qDzrDLLBhcaNI7IQ6f/ZtyAld7FL+eWYXLh+cWDvW3DlS2OIqaEeaqlBQP5GmSpP+6tyV+DKG3zjku7BdnXFfxaWVE/AdhEuA2mj5q0XuSryVYfVdheMy8h2l1s4QbcVIW8WeqSwHlbfit7VQEFLh48OGSYV8Mkql5JlUXQ38YUEUSxyM9IuwXZt7DXuGj9BGAhx6ZnVlZsxTU0Z715HPKzVzmTwClo3kvQeRXWtZ7tqME8GKpfhmleFvqtXDToxWotFQ8ZMpw0/ZDhs2n3gZ0vkvTUbclfgyhtG+9JuwXZFyJXWv0EYJ09oYoP5hxDG35bpJgGQPb3gDDzVYLml0XZlNuxeGAk13RgKyKiUDavBrY9blR0ospVwyc4JvIXoiYDs+o8YWy8NnPAmkOtASWanroXpd2C7GtvxEn4844NVbxvbKz80u63H6QzADKZcxVWUpPjncv8Ud6KqImb2VNYG+8BYDYUuHwiwkbvtufOkZOHrrac4ZxTZvxhqIptxKLgJHxqBHOPSrsF2bd1wLho/QRguzeGKt47XsBcr8SRVeTaSp1b710lD9jwiTiyp9PlJCIfy88EJHDvCAC2lA9h0kXYjr1t6NCNbN/vrC9kQDWvaJ6Tpyohsl8TTVsPXp12C7Jvx1PGResnAKt2M1xnHjC/njEvtcgZFlS/V+9d5dRaLzig6W/JKqClw6ho73gBJ1QSn3qlUHGyyv31zt1ThgsvGBpGm9sHJcN0L/YJ/Yni04oawSQXcwUabsQ5YCHS/1dMGSypMgO+VIZL3HJy+fE4nqvUZ/pG5izj3Npl1t3Z6amj6SXEKuOzRq3huJA9yBc4er8qPrvHkV/MbSGLy/91vlSqzfwvAMi11OY6RGnKca5joJL5NIz6CcAimPLLV9k71By1N8vjvOaY57FM5MbwoTOOiKUuajCGQ5CAbdVhNat3bS4YGsaCCdtQouHzYW7RMIVFHBaZLz0nmrYW96Xdguwb32dctH4CsKZwm3EDwPn7rU/c1aZ6qDTBq3xQolWP+lzjaNN2ugyHthTbcOGpS83OJ7Irmfck9Y4XsGDSY+6Vx//vYQHB0j/u3OV7vtv2YR8ZNH8hrNriU2p3LaK0hAguGlaP+YLA+gnAWs0SRdpdvncQSyurtqKkeqiwXvAjf+L2qK/JbWJJlMz7qmgCN+OmiDavDn3KP+6wBUx+Q/uGwVLveAFnjZgPg/YUi55Dmol4/KbaXYsoLfteSrsF2feGvzEuWj8BmGGmbqfBpir+BLY3lmMLBfzR8MiU+yPVZ71J5dyygJv01rkMa544lMNLeH20dlFji5B8sXe8cHBz7oAPDceOF4yCpTP8tg9zPC9OHB83bmss9r1Q2+sRpWEG84AF2vaEcdH6CcDaZoY+5Yq5cw4mZA2zgtFlheHhE5MYyjUden8VgZgAmOuzvYpv21yGY94ytgvvXcRPLxRBiMSCdrOdqxc9Ehx/Yfceo/rWV/Z09EpmbPX2imrtVj9WcHELNYJZy9JuQfY98n3jovUTgDWH33D3rs728g2XFYO+XMr0lEpY59z0t4qNuo8ZL+Aje61hmTBz1JzXtH7elxP8yaLN4dpDBITaWsPOmVTV7flwbKGAlYYrFQOfTdZz5KLB/bVb/VjBBJXUCEbNPiw1tBC7YtRPALb8zNCnHO41UdjJb2si6/YwBHubzFeK+dYPYKJJcMHQ8MFeBJdeLaM6re+vHC3g8N43RmsfNbYIzy2gnFS1CfBMDtwE4Au79hin0DqQpNUn2fDFe/fh8r2DkdpblXxX7a9JVGuD/BAfqH22cdH6CcCe/Z/Qp3xhlxXNBwU0Hr1K9tv97W04x+sNIoKWUrmO0/zmvXhxmQP2os4tJ9QkCmvJqkirjHvHC/jCzt0Ht/yyBUpLJybwvRe3heqp6h0v4LSR0fIP9l5h63/93fuH0gm+AE5OpsbAHR+Cje41Llo/AdiOP4Q+pXe8gFbnXnMRzSoVcfneQbQGbW/kxmVoZsLatHFrc25qmaChTbcAsa3N/xwiP53zI512wdAwvv/SNhzhWG389qER9BbCv5j/+/adeMv+YbSWbD28qmhWnbKNUU0VajznjCgN3HQ+2K5njIvWTwAWIlGk3WRM+z1W5rGc4LUZcSXfl2FPVmVT4w2tVWbYtq63pMg3CKpCFaufescLuHLnbrSpIqeKNtWDyVojuGrXbqx5fgv+euduvMJKT3Ht1u21n/dld8yb07s2Ua1EnWbTSEoenTAuwo8rZFXEbPHtpRKGc4b/VJUgymU+1sH5KSHa6LOlSmUV1xGFCaxpr+Kf3mrzY51FoP86oO+i6HVR45pzJPDCmsin944X8B9bt+OhtjacMjYWS7B0wdBwbXN9ETW6MSZiDXTkWcZF66cHLGLX6OJJ82gVgPtkeFuvViFMXjGPoHHpxMSBN6jL9jr2wYtoe3NzqOWxRIfY/WzVVfSOF/DRwX3p9lQl5Ymb024BUfKK5vscNqz5xxkXrZ8ALOIqpPDJK3Do5F/rqzLHqm/UcNK8VzJJlJfRV/SOF9Bin1fmVW/A/e/YPwR0L/RvE5EXzv3w55Y0majeTI6m3YJsa2oOtWo8lQBMRL4mIk+KyDoRuUlEZlVdaYRErADw9v1D5RsuqSU8OXvB9OAGxDNckrQG1uG4b70jn9ix9l46t/N8gjkAQKmEdw6NAKdf5t8mIk/xbAhftxigUiPo4od4X6/+ZKhsA2n1gN0B4BWqeiKAPwD4y6prHNwS6bQLhoYhpNi/+gAAEvpJREFUzpWQIeeT5bV0YFilb2zsYO6joOz6HkHfrtyhD8s7nEGisz637Pi2Tb57SiU05WcwDQVFF/H51TCYCZ8awYT5fqwNaf0toYqnEoCp6u2qWsmC+gCAxVVXWog+GXeRPQAL0xPmone8gOX2PfA8Ju0fOO522+Fo58pKW3Dly6pzf1MOKAwFlyfyMhl91WJD4NAMNQKugvS3ewOwebVx8SzMAfsQgNuqrqU1eibq4+z5iEyDG5u5xUN70JbbAya/JK4Vjus56zuQw8stYPNrq3Xs5LExbpVC0fVfBwxtS7sV2Sb1s6CcyBPfR4Ktvd64aGIBmIjcKSKPuXydZyvzeQCTAH7gU88lItIvIv07duzwvuCy0yK39eLB/eX0EX5zqXyOfWTw0KW5F+/bj1zluPPLrW7bzwLgXMfS+gM5k+xtCarL0l0s4tvbdgIr3+N6nCjQeq7wC6QhV1MTTUfzjkm7BdOA+RSmxAIwVX2Dqr7C5etmABCRDwJ4K4A/VvXuxlHVq1W1T1X75s2b533B0z/lkYw1+I+xcryAiwdtG1/bvzvvcwRTp42M4l2OgGnleAHXvbgNF+wbwlnDIzhirHgw474zGLN9FwX+z87dUzYn7h0v4OK9jva5tdXRviZVfHvbTpSkJfJ+fkRYcV5wmUbXNiPtFhAlb+WFkZOexyofsOiu51jguLcCnQtQjgGC4oAm7zKHnexzmjNRehOw8n0B1zoolX5zEXkTgM8BeK2qxjOrb8kq4EP/Xe7+G9oBdM0r/7MAwMZ7gO1PAk/fDkyMAd3zgTlHATueBDrmYm9+IU7bsBmdE3txT2czZGw+UOrAzuETMTv/LCZm/AFtY7OBXAHDLfuxu3UCJQHO3zeMj+0Zwa6jzkfPy04sJ6nbug6r287ATY+8gDePrsZtxVVYW/pfAIAV7feiZ+6v8HxbASURLBnpxuxCO+YVC9iRa8WGoVdjzfg4LljxHLDwRGB8H/DUfwPDO3DRzgnMnRjCbzrbMbHvBOiM9Xi6rYT5kyV0TubRNjYbw227sCM/iVmTwDHjglfua8djY2ei7bw/xXGcgE9R9V0E9H8H2Ppo2i3Jrld+MO0WECWv8j573zeBnU8Dza3A2GB5BGlipPyetewM4Pjzyj3nO58BBp8vn5vLlz+o9P4xcNxbynXsf6kcTO18EmhuK29kvfxMYNfTwK5ngY655Tq3ry9nmO9eBLz7e+V2bF59sB09RwMtHcAL/cCKc4GzrzzY5s2ryzFA+1xgdFf5+zN3ABvvA4rj5UBt1UcPlnnmdmDrY0BLO3Dqx63Xv+vKv48CePHhcidH38Xl6/zso8BTtwFzlgNv+XqoxW7i0/mUGBF5BkAewC7rrgdU9dKg8/r6+rS/vz+RNi2/4tZI573x+AW4+gN9U+5/2RW3+ibF99KVz+GxK9805f4PXPMg7n56Z+j6BMBzV70lQkuILJtXA9ecnXYrsq3vQ8Bb/zntVhBRxojIGlWdGiQgvVWQR6nqElXttb4Cg6+kvebonkjnfey1R7rf/5ojItX3V+cc73r/qUfMjVTfysXR8qMRHbDxnrRbkH1bHkq7BUQ0zWRhFWQmfP/Dp2JWe7gR2SYAJy+b7XrsinNWhK4PAC48danr/a86Ym6k7S5//skzwp9EZMf5g8F2PZd2C4hommEAZjPwxT8KVb5vuXvwVfHZN60IVV9nq/fkxpOXzcbfnn9CqPpyTF5OVBsH0hoSEZlhAFaFz73ZP8Dy6s3y8ievWhZY37I5Hcb1efXOEYUSIq9Nw2ppT7sFRDTNMACLqKerNdYAp6s1hyvOCe4xm9XhXPbqLShAJDLDrtRA+e60W0BE0wwDMAfTeVuXn32sUbl8s9mf+K/e4j753uk9p5j1ql36miPYA0bxCJHXpmG1cbELEYXDAMzBZN7Wpa85wnh48eLTlgeWOb/3MOP6Ljx1KZqb/HskVizsNupNIzKyZFVw4sNGNzqYdguIaJphAOZw4alLcX7vYZ7HexfPDBXcXHHOCvT6pIJY2J3HN957Uqg2fuSMl/ke/8rbw03WJ/K1eTUwzgDDV6v53EwiIoABmKtvvPck/J1LENPVmouU1uHnnzwDl77mCNeZNN96v882Bx6uOGeFZ5DIoUeKHfOABev2/tBGROSGAZiHC09dio1XvQW9i2eiuUnQu3gmHvvS1Az1pq44ZwWeu+otOL/3MHTlc1ixsBs/+/hpkYOlb7z3JPzs46dhxcJutOYE87pa8XdvP4FDjxS/9mhJgBvKjifTbgERTTOp7AU5ncSdyDTscKOfk5fNxm2XvSa2+ohcje4KLtPoioW0W0BE0wx7wIjI39i+tFuQfYWhtFtARNMMAzAi8sc5YMGazPPzEREBDMCIKEj3wrRbkH2lYtotIKJphgEYEfk76o1ptyD7imNpt4CIphkGYETkj5Pwg2naDSCi6YYBGBH5W35m2i3IPuF+mUQUDgMwIvK3ZBWQa027FdkmfCklonD4qkFEwRZweytfpYm0W0BE0wwDMCIKtu+FtFtARFRXGIARUbDRPWm3gIiorjAAI6Jg845JuwVERHWFARgRBev7SNotICKqKwzAiMjf5tXAbZ9JuxVERHWFARgR+Vt7PVAspN2KjGMeMCIKhwEYEQVgcBGMfyMiCocBGBH5W7gy7RZMA6W0G0BE00wqAZiIfFlE1onIgIjcLiKHpdEOIjKwdSDtFhAR1Z20esC+pqonqmovgF8C+OuU2kFEQbb0p90CIqK6k0oApqr7bD92AtA02kFEBva9lHYLiIjqTnNaFxaRvwXwAQCDAF7vU+4SAJcAwNKlS2vTOCI6aN4xwKadabeCiKiuJNYDJiJ3ishjLl/nAYCqfl5VlwD4AYBPetWjqlerap+q9s2bNy+p5hKRlzdcmXYLsu+Is9JuARFNM4n1gKnqGwyLXg/gVgBfTKotRFSFJauAD98B3Ho5sPUxHJgx0GLNHmhqAQpDQFMOWHBC+fbkGFCcAIZ3lO+XXPm+Rb3AKz8A3P4FoLDfdpEmoLWzXKY0Ub4r1wp09ABaAiZHgYlxoDhWPiY5oLkNOGwlkGsDnn8AKI6Xy0oTkJ9R/nliZOrvk8sDpUlAtVxWJ6eWaWop119pY/dhwJzlwK5ngZFdQKlYPleagOVnAh+4KZ6/NRE1jFSGIEXkaFV92vrxXABPptEOIjK0ZBVw6b3x1dd3UXx1ERFNQ2nNAbtKRI5FOXnOJgCXptQOIiIioppLJQBT1XemcV0iIiKiLGAmfCIiIqIaYwBGREREVGMMwIiIiIhqjAEYERERUY0xACMiIiKqMQZgRERERDXGAIyIiIioxkRV026DMRHZgXLi1lrqAcCdiOsDH8v6wceyfvCxrB98LKdapqquG1lPqwAsDSLSr6p9abeDqsfHsn7wsawffCzrBx/LcDgESURERFRjDMCIiIiIaowBWLCr024AxYaPZf3gY1k/+FjWDz6WIXAOGBEREVGNsQeMiIiIqMYaOgATkSUi8lsRWS8ij4vIp6z7vyYiT4rIOhG5SURm2c75SxF5RkSeEpE/Sq/1ZOfzWM4RkTtE5Gnr+2zbOXwsM0hEvisi20XkMdt9vSLygIgMiEi/iKyyHePjmFFuj6V1/59Zj9fjIvIPtvv5WGaU12NpHfu0iKiI9Nju42MZRFUb9gvAIgCvtG53A/gDgOMBvBFAs3X/VwF81bp9PIC1APIAXgbgWQC5tH8Pfvk+lv8A4Arr/iv4WGb/C8BrALwSwGO2+24H8Gbr9jkAfsfHMftfHo/l6wHcCSBv/Tyfj2X2v9weS+v+JQB+jXKOzh4+luZfDd0DpqovqerD1u39ANYDOFxVb1fVSavYAwAWW7fPA/AjVR1X1ecAPANglbNeqj2vxxLlx+x7VrHvATjfus3HMqNU9W4Au513A5hh3Z4J4EXrNh/HDPN4LD8O4CpVHbfKbLfu52OZYR6PJQD8M4DPovwcreBjaaChAzA7EVkO4CQADzoOfQjAbdbtwwFsth3bYt1HGeJ4LBeo6ktAOUgDMN8qxsdyerkMwNdEZDOAfwTwl9b9fBynn2MAnCkiD4rIXSJyinU/H8tpRkTOBfCCqq51HOJjaYABGAAR6QLwMwCXqeo+2/2fBzAJ4AeVu1xO5zLSDPF6LN2KutzHxzK7Pg7gz1V1CYA/B3CNdT8fx+mnGcBsAK8C8BkAPxERAR/LaUVEOgB8HsBfux12uY+PpUPDB2Ai0oLyG/YPVPVG2/0fBPBWAH+s1qA2ylH8Etvpi3FwKIRS5vFYbhORRdbxRQAqwx18LKeXDwKoPKY34OBwBh/H6WcLgBu1bDWAEsp7CPKxnF6ORHl+11oR2Yjy4/WwiCwEH0sjDR2AWZ+6rgGwXlW/brv/TQA+B+BcVR2xnXILgPeKSF5EXgbgaACra9lmcuf1WKL8mH3Quv1BADfb7udjOX28COC11u2zADxt3ebjOP38HOXHECJyDIBWlDdw5mM5jajqo6o6X1WXq+pylIOuV6rqVvCxNNKcdgNSdjqAPwHwqIgMWPf9FYD/i/LqjTvK7+t4QFUvVdXHReQnAJ5AeWjyE6paTKHdNJXXY3kVykMcHwbwPIALAICPZXaJyA8BvA5Aj4hsAfBFAB8F8E0RaQYwBuASgI9j1nk8lt8F8F0rnUEBwAetUQY+lhnm9liq6jVuZfm8NMNM+EREREQ11tBDkERERERpYABGREREVGMMwIiIiIhqjAEYERERUY0xACMiIiKqMQZgRGRERIoiMiAij4nIL0RklnX/YSLyU4PzhzzuP19Ejo+7vXETkV9Vfucq67lORJ4TkUs9jrv+nSJc5z0i8oyI/DKO+ogoXgzAiMjUqKr2quorUN6U9xMAoKovquq7qqj3fACZD8BU9RxV3RtTdZ9R1X+LqS5XqvpjAB9J8hpEFB0DMCKK4vewNtcVkeVWUk2ISIeI/ERE1onIj60Nl/sqJ4nI34rIWhF5QEQWiMhpAM5FeaPtARE50n4REXmbVccjInKniCyw7n+tVX7AOtZt3f9ZEXnUusZV1n1Hish/i8gaEblHRI6z7r9ORP6viNwvIhtE5F3W/YtE5G5bb9+Z1v0bRaTHun25dewxEbnM9ndYLyL/ISKPi8jtItIe9IcUkZeJyO9F5CER+bLj2Ges+9eJyJW2+/+PiDwpIneIyA9F5NOhHj0iSh0DMCIKRURyAP4XytuNOP0pgD2qeiKALwM42XasE+VdJVYCuBvAR1X1fquez1i9a8866rsXwKtU9SQAPwLwWev+T6OcXbsXwJkARkXkzSj3pp1qXeMfrLJXA/gzVT3ZOu/btvoXATgD5X1fr7LuuxDAr626VwIYsJWHiJwM4GIAp6K8ofRHReQk6/DRAL6lqi8HsBfAO13+Rk7fBPCvqnoKgK2267zRqm8VgF4AJ4vIa6yA9p0ATgLwDgB9U6skoqxr9K2IiMhcu7XN03IAawDc4VLmDJQDCqjqYyKyznasAKAyH2kNgLMNrrkYwI+lvJF6K4DnrPvvA/B1EfkByhs7bxGRNwC4trJ/q6ruFpEuAKcBuMHaVgwobzNW8XNVLQF4otK7BuAhlLfKabGOHxKAWb/jTao6DAAiciPKQeAtAJ6zlV+D8t8qyOk4GKj9J4CvWrffaH09Yv3chXJA1g3gZlUdta7/C4NrEFHGsAeMiEyNWr1Cy1AOhj7hUkZc7quY0IN7nxVh9gHw/wH4F1U9AcDHALQBgKpehfL8pnYAD1jDigLAubdaE4C9Vu9a5WuF7fi4s+2qejeA1wB4AcB/isgHQvyO9vpMf0e4tLtynb+3tfsoa+89v+sT0TTBAIyIQlHVQQD/G8CnrV4iu3sBvBsArJWNJxhUuR/lXh03M1EOhADgg5U7ReRIVX1UVb8KoB/AcQBuB/AhEemwysxR1X0AnhORC6z7RERW+jVGRJYB2K6q/wHgGgCvdBS5G8D51ny3TgBvB3CPwe/p5T4A77Vu/7Ht/l9bv0+X1a7DRWQ+yn/jt4lIm3XsLVVcm4hSwgCMiEJT1UcArMXBwKHi2wDmWUOPnwOwDsBgQHU/AvAZazL9kY5jf4Py8OE9AHba7r/MmgC/FsAogNtU9b9RHgbst4ZKKxPT/xjAh62yjwM4L6A9rwMwICKPoDw0+E37QVV9GMB1AFYDeBDAd6y/R1SfAvAJEXkI5YCzcp3bAVwP4Pci8iiAnwLoVtWHrN9zLYAbUQ5ABwFARC71Sm9BRNkiB0cEiIiqY03Qb1HVMSuY+h8Ax6hqIeWmZYaIXAfgl6oamDvNp44uVR2yevvuBnCJFRg6y70OwKdV9a1Rr0VEyeAkfCKKUweA31pDkwLg4wy+phgE8GUR6akiF9jV1hBvG4DveQRf7wHwRZQXAxBRxrAHjIiIiKjGOAeMiIiIqMYYgBERERHVGAMwIiIiohpjAEZERERUYwzAiIiIiGqMARgRERFRjf1/+nzzcLi3JpcAAAAASUVORK5CYII=\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "fig, ax = plt.subplots(1,1, figsize=(10,5))\n", + "\n", + "pl.plot(desi['TARGET_RA'], desi['TARGET_DEC'], marker='.', lw=0.0)\n", + "\n", + "pl.plot(gama['RA'], gama['DEC'], marker='.', lw=0.0)\n", + "pl.plot(match['RA'], match['DEC'], marker='.', lw=0.0)\n", + "\n", + "pl.xlim(225., 128.)\n", + "\n", + "pl.xlabel('Right ascension [deg.]')\n", + "pl.ylabel('Declination [deg.]')" + ] + }, + { + "cell_type": "markdown", + "id": "e3f53dac-e4ef-48f9-b7aa-aa2a3f842c12", + "metadata": {}, + "source": [ + "# Done." + ] + } + ], + "metadata": { + "kernelspec": { + "display_name": "DESI master", + "language": "python", + "name": "desi-master" + }, + "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.8.3" + } + }, + "nbformat": 4, + "nbformat_minor": 5 +}