This notebook will let me look at the format of the DS14a stuff and prep to cache it, for makign my covmats.
In [1]:
import numpy as np
import h5py
In [15]:
from matplotlib import pyplot as plt
%matplotlib inline
import seaborn as sns
sns.set()
In [2]:
fname = '/scratch/users/swmclau2/Darksky/ds14_a_halos_1.0000.hdf5'
In [3]:
f = h5py.File(fname, 'r')
In [6]:
f['halos'].keys()
Out[6]:
[u'subbox_000',
u'subbox_001',
u'subbox_002',
u'subbox_003',
u'subbox_004',
u'subbox_005',
u'subbox_006',
u'subbox_007',
u'subbox_010',
u'subbox_011',
u'subbox_012',
u'subbox_013',
u'subbox_014',
u'subbox_015',
u'subbox_016',
u'subbox_017',
u'subbox_020',
u'subbox_021',
u'subbox_022',
u'subbox_023',
u'subbox_024',
u'subbox_025',
u'subbox_026',
u'subbox_027',
u'subbox_030',
u'subbox_031',
u'subbox_032',
u'subbox_033',
u'subbox_034',
u'subbox_035',
u'subbox_036',
u'subbox_037',
u'subbox_040',
u'subbox_041',
u'subbox_042',
u'subbox_043',
u'subbox_044',
u'subbox_045',
u'subbox_046',
u'subbox_047',
u'subbox_050',
u'subbox_051',
u'subbox_052',
u'subbox_053',
u'subbox_054',
u'subbox_055',
u'subbox_056',
u'subbox_057',
u'subbox_060',
u'subbox_061',
u'subbox_062',
u'subbox_063',
u'subbox_064',
u'subbox_065',
u'subbox_066',
u'subbox_067',
u'subbox_070',
u'subbox_071',
u'subbox_072',
u'subbox_073',
u'subbox_074',
u'subbox_075',
u'subbox_076',
u'subbox_077',
u'subbox_100',
u'subbox_101',
u'subbox_102',
u'subbox_103',
u'subbox_104',
u'subbox_105',
u'subbox_106',
u'subbox_107',
u'subbox_110',
u'subbox_111',
u'subbox_112',
u'subbox_113',
u'subbox_114',
u'subbox_115',
u'subbox_116',
u'subbox_117',
u'subbox_120',
u'subbox_121',
u'subbox_122',
u'subbox_123',
u'subbox_124',
u'subbox_125',
u'subbox_126',
u'subbox_127',
u'subbox_130',
u'subbox_131',
u'subbox_132',
u'subbox_133',
u'subbox_134',
u'subbox_135',
u'subbox_136',
u'subbox_137',
u'subbox_140',
u'subbox_141',
u'subbox_142',
u'subbox_143',
u'subbox_144',
u'subbox_145',
u'subbox_146',
u'subbox_147',
u'subbox_150',
u'subbox_151',
u'subbox_152',
u'subbox_153',
u'subbox_154',
u'subbox_155',
u'subbox_156',
u'subbox_157',
u'subbox_160',
u'subbox_161',
u'subbox_162',
u'subbox_163',
u'subbox_164',
u'subbox_165',
u'subbox_166',
u'subbox_167',
u'subbox_170',
u'subbox_171',
u'subbox_172',
u'subbox_173',
u'subbox_174',
u'subbox_175',
u'subbox_176',
u'subbox_177',
u'subbox_200',
u'subbox_201',
u'subbox_202',
u'subbox_203',
u'subbox_204',
u'subbox_205',
u'subbox_206',
u'subbox_207',
u'subbox_210',
u'subbox_211',
u'subbox_212',
u'subbox_213',
u'subbox_214',
u'subbox_215',
u'subbox_216',
u'subbox_217',
u'subbox_220',
u'subbox_221',
u'subbox_222',
u'subbox_223',
u'subbox_224',
u'subbox_225',
u'subbox_226',
u'subbox_227',
u'subbox_230',
u'subbox_231',
u'subbox_232',
u'subbox_233',
u'subbox_234',
u'subbox_235',
u'subbox_236',
u'subbox_237',
u'subbox_240',
u'subbox_241',
u'subbox_242',
u'subbox_243',
u'subbox_244',
u'subbox_245',
u'subbox_246',
u'subbox_247',
u'subbox_250',
u'subbox_251',
u'subbox_252',
u'subbox_253',
u'subbox_254',
u'subbox_255',
u'subbox_256',
u'subbox_257',
u'subbox_260',
u'subbox_261',
u'subbox_262',
u'subbox_263',
u'subbox_264',
u'subbox_265',
u'subbox_266',
u'subbox_267',
u'subbox_270',
u'subbox_271',
u'subbox_272',
u'subbox_273',
u'subbox_274',
u'subbox_275',
u'subbox_276',
u'subbox_277',
u'subbox_300',
u'subbox_301',
u'subbox_302',
u'subbox_303',
u'subbox_304',
u'subbox_305',
u'subbox_306',
u'subbox_307',
u'subbox_310',
u'subbox_311',
u'subbox_312',
u'subbox_313',
u'subbox_314',
u'subbox_315',
u'subbox_316',
u'subbox_317',
u'subbox_320',
u'subbox_321',
u'subbox_322',
u'subbox_323',
u'subbox_324',
u'subbox_325',
u'subbox_326',
u'subbox_327',
u'subbox_330',
u'subbox_331',
u'subbox_332',
u'subbox_333',
u'subbox_334',
u'subbox_335',
u'subbox_336',
u'subbox_337',
u'subbox_340',
u'subbox_341',
u'subbox_342',
u'subbox_343',
u'subbox_344',
u'subbox_345',
u'subbox_346',
u'subbox_347',
u'subbox_350',
u'subbox_351',
u'subbox_352',
u'subbox_353',
u'subbox_354',
u'subbox_355',
u'subbox_356',
u'subbox_357',
u'subbox_360',
u'subbox_361',
u'subbox_362',
u'subbox_363',
u'subbox_364',
u'subbox_365',
u'subbox_366',
u'subbox_367',
u'subbox_370',
u'subbox_371',
u'subbox_372',
u'subbox_373',
u'subbox_374',
u'subbox_375',
u'subbox_376',
u'subbox_377',
u'subbox_400',
u'subbox_401',
u'subbox_402',
u'subbox_403',
u'subbox_404',
u'subbox_405',
u'subbox_406',
u'subbox_407',
u'subbox_410',
u'subbox_411',
u'subbox_412',
u'subbox_413',
u'subbox_414',
u'subbox_415',
u'subbox_416',
u'subbox_417',
u'subbox_420',
u'subbox_421',
u'subbox_422',
u'subbox_423',
u'subbox_424',
u'subbox_425',
u'subbox_426',
u'subbox_427',
u'subbox_430',
u'subbox_431',
u'subbox_432',
u'subbox_433',
u'subbox_434',
u'subbox_435',
u'subbox_436',
u'subbox_437',
u'subbox_440',
u'subbox_441',
u'subbox_442',
u'subbox_443',
u'subbox_444',
u'subbox_445',
u'subbox_446',
u'subbox_447',
u'subbox_450',
u'subbox_451',
u'subbox_452',
u'subbox_453',
u'subbox_454',
u'subbox_455',
u'subbox_456',
u'subbox_457',
u'subbox_460',
u'subbox_461',
u'subbox_462',
u'subbox_463',
u'subbox_464',
u'subbox_465',
u'subbox_466',
u'subbox_467',
u'subbox_470',
u'subbox_471',
u'subbox_472',
u'subbox_473',
u'subbox_474',
u'subbox_475',
u'subbox_476',
u'subbox_477',
u'subbox_500',
u'subbox_501',
u'subbox_502',
u'subbox_503',
u'subbox_504',
u'subbox_505',
u'subbox_506',
u'subbox_507',
u'subbox_510',
u'subbox_511',
u'subbox_512',
u'subbox_513',
u'subbox_514',
u'subbox_515',
u'subbox_516',
u'subbox_517',
u'subbox_520',
u'subbox_521',
u'subbox_522',
u'subbox_523',
u'subbox_524',
u'subbox_525',
u'subbox_526',
u'subbox_527',
u'subbox_530',
u'subbox_531',
u'subbox_532',
u'subbox_533',
u'subbox_534',
u'subbox_535',
u'subbox_536',
u'subbox_537',
u'subbox_540',
u'subbox_541',
u'subbox_542',
u'subbox_543',
u'subbox_544',
u'subbox_545',
u'subbox_546',
u'subbox_547',
u'subbox_550',
u'subbox_551',
u'subbox_552',
u'subbox_553',
u'subbox_554',
u'subbox_555',
u'subbox_556',
u'subbox_557',
u'subbox_560',
u'subbox_561',
u'subbox_562',
u'subbox_563',
u'subbox_564',
u'subbox_565',
u'subbox_566',
u'subbox_567',
u'subbox_570',
u'subbox_571',
u'subbox_572',
u'subbox_573',
u'subbox_574',
u'subbox_575',
u'subbox_576',
u'subbox_577',
u'subbox_600',
u'subbox_601',
u'subbox_602',
u'subbox_603',
u'subbox_604',
u'subbox_605',
u'subbox_606',
u'subbox_607',
u'subbox_610',
u'subbox_611',
u'subbox_612',
u'subbox_613',
u'subbox_614',
u'subbox_615',
u'subbox_616',
u'subbox_617',
u'subbox_620',
u'subbox_621',
u'subbox_622',
u'subbox_623',
u'subbox_624',
u'subbox_625',
u'subbox_626',
u'subbox_627',
u'subbox_630',
u'subbox_631',
u'subbox_632',
u'subbox_633',
u'subbox_634',
u'subbox_635',
u'subbox_636',
u'subbox_637',
u'subbox_640',
u'subbox_641',
u'subbox_642',
u'subbox_643',
u'subbox_644',
u'subbox_645',
u'subbox_646',
u'subbox_647',
u'subbox_650',
u'subbox_651',
u'subbox_652',
u'subbox_653',
u'subbox_654',
u'subbox_655',
u'subbox_656',
u'subbox_657',
u'subbox_660',
u'subbox_661',
u'subbox_662',
u'subbox_663',
u'subbox_664',
u'subbox_665',
u'subbox_666',
u'subbox_667',
u'subbox_670',
u'subbox_671',
u'subbox_672',
u'subbox_673',
u'subbox_674',
u'subbox_675',
u'subbox_676',
u'subbox_677',
u'subbox_700',
u'subbox_701',
u'subbox_702',
u'subbox_703',
u'subbox_704',
u'subbox_705',
u'subbox_706',
u'subbox_707',
u'subbox_710',
u'subbox_711',
u'subbox_712',
u'subbox_713',
u'subbox_714',
u'subbox_715',
u'subbox_716',
u'subbox_717',
u'subbox_720',
u'subbox_721',
u'subbox_722',
u'subbox_723',
u'subbox_724',
u'subbox_725',
u'subbox_726',
u'subbox_727',
u'subbox_730',
u'subbox_731',
u'subbox_732',
u'subbox_733',
u'subbox_734',
u'subbox_735',
u'subbox_736',
u'subbox_737',
u'subbox_740',
u'subbox_741',
u'subbox_742',
u'subbox_743',
u'subbox_744',
u'subbox_745',
u'subbox_746',
u'subbox_747',
u'subbox_750',
u'subbox_751',
u'subbox_752',
u'subbox_753',
u'subbox_754',
u'subbox_755',
u'subbox_756',
u'subbox_757',
u'subbox_758',
u'subbox_760',
u'subbox_761',
u'subbox_762',
u'subbox_763',
u'subbox_764',
u'subbox_765',
u'subbox_766',
u'subbox_767',
u'subbox_768',
u'subbox_770',
u'subbox_771',
u'subbox_772',
u'subbox_773',
u'subbox_774',
u'subbox_775',
u'subbox_776',
u'subbox_777',
u'subbox_778',
u'subbox_780',
u'subbox_781',
u'subbox_782',
u'subbox_783',
u'subbox_784',
u'subbox_785',
u'subbox_786',
u'subbox_787',
u'subbox_850',
u'subbox_851',
u'subbox_852',
u'subbox_853',
u'subbox_854',
u'subbox_855',
u'subbox_856',
u'subbox_857',
u'subbox_858',
u'subbox_860',
u'subbox_861',
u'subbox_862',
u'subbox_863',
u'subbox_864',
u'subbox_865',
u'subbox_866',
u'subbox_867',
u'subbox_868',
u'subbox_870',
u'subbox_871',
u'subbox_872',
u'subbox_873',
u'subbox_874',
u'subbox_875',
u'subbox_876',
u'subbox_877',
u'subbox_878',
u'subbox_880',
u'subbox_881',
u'subbox_882',
u'subbox_883',
u'subbox_884',
u'subbox_885',
u'subbox_886',
u'subbox_887']
In [27]:
halo_table = f['halos']['subbox_010'].value
In [28]:
halo_table[0, :]
Out[28]:
array([ 7.72566376e+01, 1.00098010e+03, 2.50001511e+02,
3.31900597e+12, 3.64358783e-01, 2.41735611e+02,
0.00000000e+00], dtype=float32)
In [29]:
#x y z m200b r2oob vmax id
In [30]:
mass_cut = halo_table[:,3] > 1e15
x_pos, y_pos = halo_table[mass_cut,0], halo_table[mass_cut,1]
In [31]:
plt.scatter(x_pos, y_pos, alpha = 0.1)
Out[31]:
<matplotlib.collections.PathCollection at 0x7f3cfc51c9d0>
In [ ]:
Content source: mclaughlin6464/pearce
Similar notebooks: