notebook.community

Edit and run 





In [1]:



    


import numpy as np
from astropy.cosmology import LambdaCDM
from astropy.io import fits



    











In [2]:



    


theta_bins = np.logspace(np.log10(2.5), np.log10(250), 21)/60
print theta_bins



    


















    






[ 0.04166667  0.05245523  0.06603722  0.08313593  0.10466193  0.13176157
  0.16587799  0.20882801  0.26289889  0.3309701   0.41666667  0.52455225
  0.66037216  0.8313593   1.04661935  1.31761569  1.65877988  2.08828014
  2.62898894  3.30970098  4.16666667]

















In [3]:



    


fname = '/u/ki/jderose/public_html/bcc/measurement/y3/3x2pt/buzzard/flock/buzzard-2/tpt_Y3_v0.fits'
hdulist = fits.open(fname)



    











In [4]:



    


hdulist.info()



    


















    






Filename: /u/ki/jderose/public_html/bcc/measurement/y3/3x2pt/buzzard/flock/buzzard-2/tpt_Y3_v0.fits
No.    Name         Type      Cards   Dimensions   Format
  0  PRIMARY     PrimaryHDU       4   ()      
  1  xip         BinTableHDU     32   200R x 5C   [K, K, K, D, D]   
  2  xim         BinTableHDU     32   200R x 5C   [K, K, K, D, D]   
  3  gammat      BinTableHDU     32   400R x 5C   [K, K, K, D, D]   
  4  wtheta      BinTableHDU     32   300R x 5C   [K, K, K, D, D]   
  5  nz_shear_bpz  BinTableHDU     32   250R x 7C   [D, D, D, D, D, D, D]   
  6  nz_shear_true  BinTableHDU     32   250R x 7C   [D, D, D, D, D, D, D]   
  7  nz_pos_zrm  BinTableHDU     31   400R x 8C   [D, D, D, D, D, D, D, D]   
  8  nz_pos_zspec  BinTableHDU     31   400R x 8C   [D, D, D, D, D, D, D, D]   

















In [5]:



    


z_bins = np.array([0.15, 0.3, 0.45, 0.6, 0.75, 0.9])



    











In [6]:



    


wtheta = hdulist[4]
# includes cross, make sure to only take autos
zbin = 1
y = []
for row in wtheta.data:
    print row
    if row[0] == zbin and row[1] == zbin:
        y.append(row[3])
        
y = np.array(y)
print y



    


















    






(1, 1, 0, 0.45695746584854985, 2.823621287743753)
(1, 1, 1, 0.36215490680577156, 3.557146025674375)
(1, 1, 2, 0.29041139011879313, 4.4823793795798936)
(1, 1, 3, 0.23118317670343497, 5.6400730550718423)
(1, 1, 4, 0.19100328415347961, 7.1046610696008692)
(1, 1, 5, 0.15100299741973047, 8.9441097105123593)
(1, 1, 6, 0.13108268339458876, 11.262638241945165)
(1, 1, 7, 0.10749353219459415, 14.177008926475571)
(1, 1, 8, 0.086800552571491038, 17.847698728733025)
(1, 1, 9, 0.075283022559461496, 22.471916045640072)
(1, 1, 10, 0.062911948852758512, 28.288981962942596)
(1, 1, 11, 0.051891487004727226, 35.614102833654002)
(1, 1, 12, 0.041154875101075204, 44.834108999916808)
(1, 1, 13, 0.030233398732561853, 56.447535720056472)
(1, 1, 14, 0.022234977712605405, 71.064011754584612)
(1, 1, 15, 0.015794000030595636, 89.46250579903645)
(1, 1, 16, 0.010559505277520233, 112.62320328217676)
(1, 1, 17, 0.0067023981889799003, 141.77573475161819)
(1, 1, 18, 0.004538967516791963, 178.48268026087032)
(1, 1, 19, 0.002665683634473909, 224.6601331679756)
(1, 2, 0, 0.062681787362790106, 2.8290407198474585)
(1, 2, 1, 0.046305558861309042, 3.5604431184209635)
(1, 2, 2, 0.036822684075634771, 4.4838863126467245)
(1, 2, 3, 0.029015685199586774, 5.6442448048837868)
(1, 2, 4, 0.024375867090192151, 7.1069292315432291)
(1, 2, 5, 0.021807222083970824, 8.9478438450845026)
(1, 2, 6, 0.014181525053317892, 11.263210470419624)
(1, 2, 7, 0.014233689264610072, 14.181790241527827)
(1, 2, 8, 0.0097318054632305247, 17.853725380478405)
(1, 2, 9, 0.0081705627571717222, 22.474547151422257)
(1, 2, 10, 0.0054650184335542154, 28.295934787910923)
(1, 2, 11, 0.004527557634235616, 35.621058746931077)
(1, 2, 12, 0.0033038218183915535, 44.843690061555591)
(1, 2, 13, 0.002230519091482897, 56.455675565218719)
(1, 2, 14, 0.00061106074651157832, 71.072736461945823)
(1, 2, 15, 0.00071091367798590491, 89.470486654555302)
(1, 2, 16, 0.0008436409504731778, 112.63293992907248)
(1, 2, 17, 0.00056859521754434232, 141.78966809679642)
(1, 2, 18, 0.00072991636780382336, 178.48698928995756)
(1, 2, 19, 0.00011970968828811287, 224.67503262595082)
(1, 3, 0, 0.0056013299813476806, 2.828897790106192)
(1, 3, 1, 0.00011015175784118856, 3.562267032088589)
(1, 3, 2, 0.00049940577241262927, 4.4838353155391388)
(1, 3, 3, 0.0022984864418030488, 5.6454226177098255)
(1, 3, 4, -0.0029135659259559737, 7.1081675524022288)
(1, 3, 5, -0.0021014798949743831, 8.9479704338620678)
(1, 3, 6, -0.0017298028105887911, 11.265056748423463)
(1, 3, 7, -0.0036759444885173542, 14.181315405257227)
(1, 3, 8, -0.0019821865162525589, 17.854417064368942)
(1, 3, 9, -0.0013061276923152876, 22.476185293434689)
(1, 3, 10, -0.0020336386191268552, 28.295092662537705)
(1, 3, 11, -0.0019037109713265709, 35.62268927813551)
(1, 3, 12, -0.00074422442238226657, 44.84607221274856)
(1, 3, 13, -0.0010562344992945613, 56.45720860460699)
(1, 3, 14, -0.00062130794547296129, 71.072656712554902)
(1, 3, 15, -1.0611418074154596e-05, 89.471438669353944)
(1, 3, 16, -0.00068679485048630748, 112.63172781923566)
(1, 3, 17, -0.0009666718188312476, 141.78441207643775)
(1, 3, 18, -0.00090881992380154595, 178.4887529822806)
(1, 3, 19, -0.00040354441223792672, 224.67670807740737)
(1, 4, 0, 0.0028716697280496179, 2.8300531350393618)
(1, 4, 1, -0.0022893815710856732, 3.5640544491068726)
(1, 4, 2, 0.0024457229598506057, 4.4843225343952904)
(1, 4, 3, -0.0046156116255739602, 5.6449597212109168)
(1, 4, 4, -0.0031916813815122858, 7.1076982102597901)
(1, 4, 5, -0.0031295038107924267, 8.9489265134878391)
(1, 4, 6, 0.0014292317457993788, 11.264893055220758)
(1, 4, 7, -0.00049778229828275331, 14.180463075134131)
(1, 4, 8, 0.00023443716293191087, 17.854214756894468)
(1, 4, 9, -8.6665085375825426e-05, 22.476584709717248)
(1, 4, 10, -0.001391852517103401, 28.296612001719364)
(1, 4, 11, -0.0020407995976366423, 35.62139268067687)
(1, 4, 12, -0.0023438420497415425, 44.844369503441825)
(1, 4, 13, -0.0014969752277719233, 56.458694336280537)
(1, 4, 14, -0.00075655579735580476, 71.072173421652522)
(1, 4, 15, -0.00073620254791828717, 89.47265984087133)
(1, 4, 16, -0.00027654673073776003, 112.63107580655631)
(1, 4, 17, -0.00022797119927273383, 141.78860830572472)
(1, 4, 18, -0.0012226187206394437, 178.48521435761225)
(1, 4, 19, -0.00082181194728008227, 224.67895791836435)
(1, 5, 0, 0.0012888120447552187, 2.8298426100587517)
(1, 5, 1, -0.0026916288290046489, 3.5641559684798114)
(1, 5, 2, -0.0078655534507456517, 4.484002317258601)
(1, 5, 3, -0.0049898315733780318, 5.6438297553706995)
(1, 5, 4, -0.00302549964186947, 7.1072287165639842)
(1, 5, 5, -0.0050666705755585033, 8.9456028393181111)
(1, 5, 6, -0.00085334260006416111, 11.267455142278603)
(1, 5, 7, -0.0056272573886391166, 14.180899584306497)
(1, 5, 8, -0.0043446939595392304, 17.853160989861983)
(1, 5, 9, -0.00046219555125933933, 22.475483225812358)
(1, 5, 10, -0.00025301207675365762, 28.295836297557496)
(1, 5, 11, -0.00073796063529523272, 35.621242861719928)
(1, 5, 12, -0.0017150032158056018, 44.842408458616042)
(1, 5, 13, 0.0003890360784539979, 56.458240650268195)
(1, 5, 14, -0.00064679561166196503, 71.075322075048561)
(1, 5, 15, 0.00019092876802030497, 89.474296972833628)
(1, 5, 16, -0.00046089029395191385, 112.63407434160955)
(1, 5, 17, -0.00093506153591567853, 141.79002478443118)
(1, 5, 18, -0.0013414562557877683, 178.47882863347388)
(1, 5, 19, -0.0010894683605736776, 224.67903345516638)
(2, 2, 0, 0.44291303299153767, 2.8251488513797227)
(2, 2, 1, 0.33615296427854585, 3.5575379288401567)
(2, 2, 2, 0.26826407014019177, 4.4795015222905743)
(2, 2, 3, 0.21152713290350522, 5.6428720953865161)
(2, 2, 4, 0.17787163093311181, 7.1042534631650831)
(2, 2, 5, 0.14418750264523444, 8.9435169899261133)
(2, 2, 6, 0.11888569066209813, 11.260317504487126)
(2, 2, 7, 0.097832428659035561, 14.176019162515644)
(2, 2, 8, 0.079691889666470214, 17.847348437676882)
(2, 2, 9, 0.064197443143085953, 22.471073940519229)
(2, 2, 10, 0.050909286513896229, 28.289741783251866)
(2, 2, 11, 0.039995038954464823, 35.61559330924122)
(2, 2, 12, 0.030602524260995811, 44.836712342195028)
(2, 2, 13, 0.022261850128123262, 56.448297946199176)
(2, 2, 14, 0.015059935709551824, 71.063846249515009)
(2, 2, 15, 0.010123216857683965, 89.465090239778817)
(2, 2, 16, 0.0060508888585274019, 112.62420831485296)
(2, 2, 17, 0.0036079869872942326, 141.77999656891794)
(2, 2, 18, 0.0015060783190812073, 178.48312460900101)
(2, 2, 19, 0.00082333720333242072, 224.67897319641474)
(2, 3, 0, 0.044512507215608975, 2.8292056321838963)
(2, 3, 1, 0.032521856339183149, 3.5615663808062292)
(2, 3, 2, 0.025053863844049222, 4.4845746267883833)
(2, 3, 3, 0.022425064569129984, 5.6452680437241352)
(2, 3, 4, 0.018159327532949968, 7.1066729996371745)
(2, 3, 5, 0.01560390793829536, 8.9482687330301207)
(2, 3, 6, 0.01220551077069159, 11.264068714279066)
(2, 3, 7, 0.010184371070229776, 14.181146493206578)
(2, 3, 8, 0.0087313787768825788, 17.85271345117966)
(2, 3, 9, 0.0069844911697630203, 22.47553830300485)
(2, 3, 10, 0.0056216725822514895, 28.295021508134855)
(2, 3, 11, 0.004577824697447319, 35.621571741159244)
(2, 3, 12, 0.0032350209995712039, 44.843513702175223)
(2, 3, 13, 0.0021704166664792894, 56.45609525483183)
(2, 3, 14, 0.0012385990147339078, 71.071178296665508)
(2, 3, 15, 0.0010005381996411908, 89.470154458032241)
(2, 3, 16, 0.00083863262760255236, 112.63450257030657)
(2, 3, 17, 0.00082243645111614006, 141.78381169750978)
(2, 3, 18, -0.00029547644072196747, 178.48324005931525)
(2, 3, 19, 2.3143413351094135e-05, 224.68075419266566)
(2, 4, 0, 0.0020106236704139089, 2.8299808488721632)
(2, 4, 1, -0.0015929297003773725, 3.561535745827582)
(2, 4, 2, 0.0029781516567348736, 4.4835314551182242)
(2, 4, 3, 0.0015865451580800701, 5.6460884987382434)
(2, 4, 4, 0.0010670978454909547, 7.1078547367378686)
(2, 4, 5, 0.0026861699170324174, 8.9484846466428589)
(2, 4, 6, 0.0020812928629945684, 11.26496724207699)
(2, 4, 7, 0.0011881996668334414, 14.181546485207658)
(2, 4, 8, 0.0017821418994638431, 17.854111440848261)
(2, 4, 9, 0.0009860789734961262, 22.477584120764288)
(2, 4, 10, 0.0018780060257800137, 28.295803709308814)
(2, 4, 11, 0.0018238682717880265, 35.62220088620429)
(2, 4, 12, 0.0017271273399191759, 44.844802615188698)
(2, 4, 13, 0.0012930919260915634, 56.457299831585047)
(2, 4, 14, 0.00059701895800699963, 71.069173261645446)
(2, 4, 15, 0.0002674470608941488, 89.470533546899205)
(2, 4, 16, 0.00057032726330596658, 112.63402695552755)
(2, 4, 17, 0.00021219551445496209, 141.78652581153003)
(2, 4, 18, -0.00017694366447875329, 178.48715934374601)
(2, 4, 19, -0.00014939671559419151, 224.68296051169736)
(2, 5, 0, 0.012432366817991721, 2.828157664626985)
(2, 5, 1, 0.015758381613352333, 3.564876706201257)
(2, 5, 2, 0.0059401807083342808, 4.4825153554253472)
(2, 5, 3, 0.0089566565712027486, 5.6468846460934792)
(2, 5, 4, 0.0056524968796522624, 7.1052773362970649)
(2, 5, 5, 0.0088271320564853144, 8.9486232183812451)
(2, 5, 6, 0.0053970161197038841, 11.265606214780792)
(2, 5, 7, 0.0020630245645845125, 14.178926769002729)
(2, 5, 8, 0.0022030622111344446, 17.852222208038604)
(2, 5, 9, 0.0023802518024306496, 22.477533580698399)
(2, 5, 10, 0.0022802232466613915, 28.295760404040838)
(2, 5, 11, 0.0014536199540971224, 35.62162867619265)
(2, 5, 12, 0.00098546364232447725, 44.84527214280785)
(2, 5, 13, 0.0015189234619416754, 56.45956259007054)
(2, 5, 14, 0.00099012331025344455, 71.069703774237951)
(2, 5, 15, 0.00013536486665052892, 89.472178900315967)
(2, 5, 16, -0.00015807832250034913, 112.63302742357389)
(2, 5, 17, -0.00010291963507590373, 141.78332527752318)
(2, 5, 18, 5.3351718442570922e-05, 178.48564047838306)
(2, 5, 19, -0.00034124686360398267, 224.67530007392523)
(3, 3, 0, 0.33440796444369741, 2.8263192221786815)
(3, 3, 1, 0.26551515760519578, 3.5591669737910978)
(3, 3, 2, 0.21719484855920707, 4.48167658532918)
(3, 3, 3, 0.1759839874757172, 5.6426126297821018)
(3, 3, 4, 0.14895634882254413, 7.1043198779289281)
(3, 3, 5, 0.12267721207560284, 8.9449236262009624)
(3, 3, 6, 0.10246548414525539, 11.260378716815099)
(3, 3, 7, 0.084659513700350228, 14.177500964528258)
(3, 3, 8, 0.069568588662311664, 17.848773458050051)
(3, 3, 9, 0.056076892075017766, 22.470881566932583)
(3, 3, 10, 0.044433620063752781, 28.289616648157651)
(3, 3, 11, 0.034678430853648523, 35.615845576650642)
(3, 3, 12, 0.025034874063516788, 44.837261007221706)
(3, 3, 13, 0.017192189587774055, 56.448432558843827)
(3, 3, 14, 0.011393862418196196, 71.06689910945417)
(3, 3, 15, 0.0071092797633219816, 89.466539327148084)
(3, 3, 16, 0.00391830705835305, 112.62748057848884)
(3, 3, 17, 0.0022793010465048591, 141.78231748399696)
(3, 3, 18, 0.00125068012089006, 178.4814188749128)
(3, 3, 19, 0.0010765775455649077, 224.67326602439317)
(3, 4, 0, 0.056458524498937776, 2.8289169621643504)
(3, 4, 1, 0.045304256112249038, 3.5619931920510286)
(3, 4, 2, 0.03839858780457385, 4.4843756115063895)
(3, 4, 3, 0.033048464578671125, 5.6456564900322066)
(3, 4, 4, 0.027537529888630615, 7.107152017313199)
(3, 4, 5, 0.021335503548855907, 8.9469771182934572)
(3, 4, 6, 0.017816741322095979, 11.264109200382897)
(3, 4, 7, 0.015539607139098839, 14.180915560088378)
(3, 4, 8, 0.012710807240550032, 17.851866490055759)
(3, 4, 9, 0.0097550510160876545, 22.474750783075081)
(3, 4, 10, 0.007207857917120768, 28.29509815963242)
(3, 4, 11, 0.0051736278875759286, 35.619810362957544)
(3, 4, 12, 0.0038747977422150894, 44.843081736608887)
(3, 4, 13, 0.0023147623964514465, 56.456172216258352)
(3, 4, 14, 0.0015166463150183276, 71.070491496092984)
(3, 4, 15, 0.0009932275740215026, 89.471874149455445)
(3, 4, 16, 0.00096634537160991349, 112.63146663337781)
(3, 4, 17, 0.00096238896969107028, 141.78439737665724)
(3, 4, 18, 0.00030105088701182643, 178.48491599669507)
(3, 4, 19, 0.00021999674000596265, 224.67564715642703)
(3, 5, 0, 0.0049090151357742929, 2.8298305989271815)
(3, 5, 1, 0.0011535736151568302, 3.5618042733490829)
(3, 5, 2, 0.0026685773082431798, 4.4837958014153765)
(3, 5, 3, 0.00099534673116967193, 5.6457735472808706)
(3, 5, 4, -0.00051450675754333066, 7.1073140582202239)
(3, 5, 5, 0.0026955965344334378, 8.9483861353007175)
(3, 5, 6, 0.00068645615772977803, 11.263062802327864)
(3, 5, 7, 0.002105088410889236, 14.181274924178155)
(3, 5, 8, 0.00056288712410386833, 17.853176920070933)
(3, 5, 9, 0.00039456598460191019, 22.473654872253892)
(3, 5, 10, 0.00044328222607363814, 28.296258113202562)
(3, 5, 11, -0.00015086515594202457, 35.622917224640524)
(3, 5, 12, 0.00021235751853919505, 44.845335308688995)
(3, 5, 13, -0.00063488310844303671, 56.455281400912583)
(3, 5, 14, -0.00016012746799158961, 71.071994262327962)
(3, 5, 15, -1.7119701655608174e-05, 89.471042685048829)
(3, 5, 16, 0.0002998088129945503, 112.63363474205515)
(3, 5, 17, -0.00012287504010779907, 141.78281039483767)
(3, 5, 18, 0.00020278637935601551, 178.48306987248662)
(3, 5, 19, 2.9988158449982913e-05, 224.67513503843836)
(4, 4, 0, 0.4977577202543142, 2.8253563952893201)
(4, 4, 1, 0.38686650327976574, 3.5586385362986075)
(4, 4, 2, 0.3146338515944534, 4.4814757002664267)
(4, 4, 3, 0.26243109599924991, 5.6400758580894825)
(4, 4, 4, 0.21364762964600742, 7.1023438716826792)
(4, 4, 5, 0.17483390560442674, 8.9411794822614485)
(4, 4, 6, 0.14916900082099729, 11.257939649691856)
(4, 4, 7, 0.11882015474795091, 14.174236986442191)
(4, 4, 8, 0.092226480069477798, 17.845632763592679)
(4, 4, 9, 0.073330114587813244, 22.469009358428966)
(4, 4, 10, 0.057160508981274727, 28.287419850654732)
(4, 4, 11, 0.044176912207443859, 35.614297436075738)
(4, 4, 12, 0.032000819135871557, 44.835405954882368)
(4, 4, 13, 0.019189757900171821, 56.446261753609349)
(4, 4, 14, 0.011406198953262951, 71.064499696279839)
(4, 4, 15, 0.0067631469773993614, 89.468320922360391)
(4, 4, 16, 0.0033608553061440978, 112.62571649979526)
(4, 4, 17, 0.00093822191151461968, 141.78668436600051)
(4, 4, 18, 0.00062074717807965806, 178.49062627628948)
(4, 4, 19, 7.6409194132451969e-05, 224.67401255972001)
(4, 5, 0, 0.02283217764535557, 2.8285452603478398)
(4, 5, 1, 0.027143838640308272, 3.5619017821749299)
(4, 5, 2, 0.018387630495069491, 4.4851765197834679)
(4, 5, 3, 0.02159128222451609, 5.644020869597604)
(4, 5, 4, 0.021017377400627923, 7.1068701182331786)
(4, 5, 5, 0.009376205078539538, 8.950406337944111)
(4, 5, 6, 0.012756989313495386, 11.263437019739532)
(4, 5, 7, 0.0080085891061448637, 14.178886618292546)
(4, 5, 8, 0.004974304982678839, 17.854663316574698)
(4, 5, 9, 0.0057847285280170305, 22.473770669749332)
(4, 5, 10, 0.0040973703756086253, 28.292746240223153)
(4, 5, 11, 0.0056567484389360902, 35.621423299637279)
(4, 5, 12, 0.0038003356881214983, 44.843195490948439)
(4, 5, 13, 0.0017221442691083275, 56.451678056046802)
(4, 5, 14, 0.0012566804751472302, 71.071150410008499)
(4, 5, 15, 0.00051253492338500998, 89.469962288494003)
(4, 5, 16, 0.00049177902382400562, 112.6313291792479)
(4, 5, 17, -8.1911068198835247e-05, 141.78235674330736)
(4, 5, 18, -0.0001310947899328966, 178.48471254612707)
(4, 5, 19, -0.00013983307587455429, 224.68041785477422)
(5, 5, 0, 0.2931836225473517, 2.8292919298921726)
(5, 5, 1, 0.31276971520235042, 3.5610325719402027)
(5, 5, 2, 0.22125858742370993, 4.4823911146670552)
(5, 5, 3, 0.22558600334517268, 5.6437938455310235)
(5, 5, 4, 0.17879430286481829, 7.104666248052407)
(5, 5, 5, 0.15007387991959195, 8.9462448717272913)
(5, 5, 6, 0.1198475414172141, 11.256639676708279)
(5, 5, 7, 0.092632477839272309, 14.17089431008857)
(5, 5, 8, 0.074423451276596581, 17.846801005157428)
(5, 5, 9, 0.063349562008801388, 22.469167238787875)
(5, 5, 10, 0.045179037099225224, 28.293220720317631)
(5, 5, 11, 0.033846345339259266, 35.616195479448244)
(5, 5, 12, 0.023904136471663152, 44.835232624377085)
(5, 5, 13, 0.01620018253437381, 56.442361953352027)
(5, 5, 14, 0.0093013319703546954, 71.063864692850544)
(5, 5, 15, 0.005390211563153998, 89.472011493360071)
(5, 5, 16, -0.00054851205580428125, 112.62485839271983)
(5, 5, 17, 0.00054618728605118822, 141.78431939427128)
(5, 5, 18, 0.00072401109893683441, 178.47309233145495)
(5, 5, 19, -0.001100387317727151, 224.67519354803539)
[ 0.45695747  0.36215491  0.29041139  0.23118318  0.19100328  0.151003
  0.13108268  0.10749353  0.08680055  0.07528302  0.06291195  0.05189149
  0.04115488  0.0302334   0.02223498  0.015794    0.01055951  0.0067024
  0.00453897  0.00266568]

















In [ ]:



    


np.savetxt('/u/ki/swmclau2/Git/pearce/bin/mcmc/buzzard2_wt_%d%d.npy'%(zbin,zbin), y)



    











In [7]:



    


hdulist.info()



    


















    






Filename: /u/ki/jderose/public_html/bcc/measurement/y3/3x2pt/buzzard/flock/buzzard-2/tpt_Y3_v0.fits
No.    Name         Type      Cards   Dimensions   Format
  0  PRIMARY     PrimaryHDU       4   ()      
  1  xip         BinTableHDU     32   200R x 5C   [K, K, K, D, D]   
  2  xim         BinTableHDU     32   200R x 5C   [K, K, K, D, D]   
  3  gammat      BinTableHDU     32   400R x 5C   [K, K, K, D, D]   
  4  wtheta      BinTableHDU     32   300R x 5C   ['K', 'K', 'K', 'D', 'D']   
  5  nz_shear_bpz  BinTableHDU     32   250R x 7C   [D, D, D, D, D, D, D]   
  6  nz_shear_true  BinTableHDU     32   250R x 7C   [D, D, D, D, D, D, D]   
  7  nz_pos_zrm  BinTableHDU     31   400R x 8C   [D, D, D, D, D, D, D, D]   
  8  nz_pos_zspec  BinTableHDU     31   400R x 8C   [D, D, D, D, D, D, D, D]   

















In [8]:



    


nz_zspec = hdulist[8]
N = 0#np.zeros((5,))
for row in nz_zspec.data:
    N+=row[2+zbin]
    #for idx, n in enumerate(row[3:]):
    #    N[idx]+=n
print N



    


















    






216485.0

















In [9]:



    


cat_filename = '/afs/slac/u/ki/jderose/public_html/bcc/catalog/redmagic/y3/buzzard/flock/buzzard-0/a/buzzard-0_1.6_y3_run_redmapper_v6.4.20_redmagic_highdens_0.5-10.fit'



    











In [10]:



    


hdulist = fits.open(cat_filename)
catalog = hdulist[1].data



    











In [11]:



    


N = np.sum(np.logical_and(z_bins[zbin-1] < catalog['ZREDMAGIC'], catalog['ZREDMAGIC'] < z_bins[zbin] ))



    











In [12]:



    


N



    


















    
Out[12]:








226994

















In [13]:



    


area = 5063 #sq degrees
full_sky = 41253 #sq degrees



    











In [14]:



    


buzzard = LambdaCDM(H0=70, Om0=0.286, Ode0=0.714, Tcmb0=2.725, Neff=3.04)
#volIn, volOut = buzzard.comoving_volume(z_bins[zbin-1]), buzzard.comoving_volume(z_bins[zbin])
volIn, volOut = buzzard.comoving_volume(z_bins[zbin-1]), buzzard.comoving_volume(z_bins[zbin])

fullsky_volume = volOut-volIn
survey_volume = fullsky_volume*area/full_sky
nd = N/survey_volume
print nd



    


















    






0.000298352370804 1 / Mpc3
np.savetxt('/u/ki/swmclau2/Git/pearce/bin/mcmc/buzzard2_nd_%d%d.npy'%(zbin, zbin),np.array(nd.value).reshape(1, 1)) 




In [16]:



    


mbins = np.logspace(np.log10(5e12), 15, 21)
mbc = mbins[1:]+mbins[:-1]



    











In [17]:



    


hod_measurement_fname = '/u/ki/jderose/notebooks/des/redmagic_hod_buzzard1_v1.3_mean_5e12.npy'
hod = np.load(hod_measurement_fname)



    











In [18]:



    


hod.shape



    


















    
Out[18]:








(20, 5, 3)

















In [19]:



    


mbins.shape



    


















    
Out[19]:








(21,)

















In [20]:



    


cen_hod = hod[:,0,1]/hod[:,0,0]
sat_hod = hod[:,0,2]/hod[:,0,0]

plt.plot(mbc, cen_hod)#+sat_hod)
plt.xscale('log')
plt.xlim(5e12, 1e15)



    


















    






---------------------------------------------------------------------------
NameError                                 Traceback (most recent call last)
<ipython-input-20-f79abad96f61> in <module>()
      2 sat_hod = hod[:,0,2]/hod[:,0,0]
      3 
----> 4 plt.plot(mbc, cen_hod)#+sat_hod)
      5 plt.xscale('log')
      6 plt.xlim(5e12, 1e15)

NameError: name 'plt' is not defined
















In [ ]:

Content source: mclaughlin6464/pearce

Similar notebooks:

notebook.community | gallery | about