In [1]:

    

import h5py
import numpy as np


    

In [2]:

    

fname = '/home/users/swmclau2/scratch/DES_emu.hdf5'


    

In [3]:

    

f = h5py.File(fname, 'r')


    

In [9]:

    

ids = f['Box_000']['z=0.550']['halos'][:,0]


    

In [20]:

    

mass = f['Box_000']['z=0.550']['halos'][:,1]


    

In [25]:

    

x = f['Box_000']['z=0.550']['halos'][:,2]


    

In [11]:

    

from collections import Counter


    

In [12]:

    

c = Counter(ids)


    

In [17]:

    

len(c), len(c)*16


    

    
Out[17]:





(47477, 759632)

In [16]:

    

len(ids)


    

    
Out[16]:





729494

In [36]:

    

selected_halo = ids == 2


    

In [37]:

    

mass[selected_halo]


    

    
Out[37]:





array([  1.50586800e+13,   9.46545600e+12,   8.60496000e+12,
         7.74446400e+12,   7.14211680e+12,   5.59322400e+12,
         5.24902560e+12,   4.90482720e+12,   3.87223200e+12,
         3.78618240e+12,   3.61408320e+12,   3.52803360e+12,
         2.92568640e+12,   2.83963680e+12,   2.58148800e+12,
         2.58148800e+12])

In [38]:

    

x[selected_halo]


    

    
Out[38]:





array([ 754.42022705,  751.40588379,  252.80305481,  250.41143799,
          2.90892649,  500.62332153,    5.09253693,    2.3887012 ,
        252.60365295,  750.7487793 ,  501.49395752,    1.44789207,
        253.6940918 ,  751.0380249 ,  503.23608398,  500.82861328])

In [45]:

    

mass.max()/mass.min()


    

    
Out[45]:





122.63333333333334

In [ ]: