In [1]:

    

import pandas as pd

guo = pd.read_table('/Users/user/Code/Pangloss/data/GGL_los_8_0_0_0_0_N_4096_ang_4_Guo_galaxies_on_plane_27_to_63.images.txt')


    

In [2]:

    

guo.columns


    

    
Out[2]:





Index([u'GalID', u'HaloID', u'SubhaloID', u'Type', u'PlaneNumber', u'z_spec',
       u'pos_0[rad]', u'pos_1[rad]', u'Dc_los[Mpc/h]', u'M_Halo[M_sol/h]',
       u'M_Subhalo[M_sol/h]', u'M_Stellar[M_sol/h]', u'mag_SDSS_u',
       u'mag_SDSS_g', u'mag_SDSS_r', u'mag_SDSS_i', u'mag_SDSS_z', u'mag_J',
       u'mag_H', u'mag_K'],
      dtype='object')

In [3]:

    

len(guo['GalID'].unique())


    

    
Out[3]:





287818

In [4]:

    

len(guo['SubhaloID'].unique())


    

    
Out[4]:





262598

In [5]:

    

len(guo['HaloID'].unique())


    

    
Out[5]:





262598

In [6]:

    

guo[guo['SubhaloID']==32002300000469]


    

    
Out[6]:






  

    

      

      
GalID
      
HaloID
      
SubhaloID
      
Type
      
PlaneNumber
      
z_spec
      
pos_0[rad]
      
pos_1[rad]
      
Dc_los[Mpc/h]
      
M_Halo[M_sol/h]
      
M_Subhalo[M_sol/h]
      
M_Stellar[M_sol/h]
      
mag_SDSS_u
      
mag_SDSS_g
      
mag_SDSS_r
      
mag_SDSS_i
      
mag_SDSS_z
      
mag_J
      
mag_H
      
mag_K
    
  
  

    

      
7
      
23003339000771
      
23003339001608
      
32002300000469
      
2
      
32
      
1.99478
      
-0.030579
      
-0.030002
      
3783.62
      
1.155350e+13
      
0.000000e+00
      
8.061260e+09
      
26.2996
      
25.5214
      
25.5081
      
25.3058
      
24.7935
      
23.8448
      
23.4391
      
23.0691
    
    

      
9
      
23003339001406
      
23003339001608
      
32002300000469
      
2
      
32
      
1.99143
      
-0.030633
      
-0.029939
      
3783.72
      
1.155350e+13
      
0.000000e+00
      
6.010610e+09
      
26.1708
      
25.5023
      
25.4905
      
25.3534
      
24.9863
      
24.0450
      
23.6678
      
23.3826
    
    

      
10
      
23003339001615
      
23003339001608
      
32002300000469
      
2
      
32
      
1.99222
      
-0.030625
      
-0.029962
      
3783.82
      
1.155350e+13
      
0.000000e+00
      
2.007910e+10
      
26.4384
      
25.4217
      
25.4436
      
25.1767
      
24.5282
      
23.4385
      
22.9097
      
22.4723
    
    

      
12
      
23003339001583
      
23003339001608
      
32002300000469
      
2
      
32
      
1.98974
      
-0.030629
      
-0.030004
      
3783.91
      
1.155350e+13
      
0.000000e+00
      
2.421100e+10
      
25.9685
      
24.9399
      
24.8657
      
24.6463
      
24.0885
      
22.7949
      
22.3892
      
22.0492
    
    

      
46
      
23003339001659
      
23003339001608
      
32002300000469
      
2
      
32
      
1.98626
      
-0.030607
      
-0.030033
      
3783.75
      
1.155350e+13
      
0.000000e+00
      
1.273490e+10
      
27.2130
      
26.1418
      
26.0659
      
25.7819
      
25.0998
      
23.8364
      
23.3265
      
22.8836
    
    

      
56
      
23003339001242
      
23003339001608
      
32002300000469
      
0
      
32
      
1.99093
      
-0.030600
      
-0.030025
      
3783.79
      
1.155350e+13
      
1.015140e+13
      
5.455610e+10
      
26.1306
      
24.9598
      
24.9767
      
24.6975
      
23.9679
      
22.9020
      
22.3320
      
21.7575
    
    

      
210
      
23003339005111
      
23003339001608
      
32002300000469
      
2
      
32
      
1.98261
      
-0.030610
      
-0.030027
      
3783.85
      
1.155350e+13
      
0.000000e+00
      
3.053020e+09
      
26.1573
      
25.8840
      
25.8328
      
25.8362
      
25.6672
      
24.7054
      
24.3242
      
24.0957
    
  

In [9]:

    

guo['pos_0[rad]'].describe()


    

    
Out[9]:





count    287822.000000
mean         -0.026113
std           0.005050
min          -0.034898
25%          -0.030429
50%          -0.026116
75%          -0.021708
max          -0.017453
Name: pos_0[rad], dtype: float64

In [10]:

    

guo['pos_1[rad]'].describe()


    

    
Out[10]:





count    287822.000000
mean         -0.026166
std           0.005034
min          -0.034898
25%          -0.030497
50%          -0.026196
75%          -0.021799
max          -0.017453
Name: pos_1[rad], dtype: float64

In [14]:

    

guo['M_Subhalo[M_sol/h]'].idxmax(axis=1)


    

    
Out[14]:





96995

In [17]:

    

guo[96995:96996]


    

    
Out[17]:






  

    

      

      
GalID
      
HaloID
      
SubhaloID
      
Type
      
PlaneNumber
      
z_spec
      
pos_0[rad]
      
pos_1[rad]
      
Dc_los[Mpc/h]
      
M_Halo[M_sol/h]
      
M_Subhalo[M_sol/h]
      
M_Stellar[M_sol/h]
      
mag_SDSS_u
      
mag_SDSS_g
      
mag_SDSS_r
      
mag_SDSS_i
      
mag_SDSS_z
      
mag_J
      
mag_H
      
mag_K
    
  
  

    

      
96995
      
76001684000066
      
76001684000066
      
46007200000197
      
0
      
46
      
0.622979
      
-0.021112
      
-0.024584
      
1626.2
      
1.825550e+14
      
2.126830e+14
      
4.750840e+11
      
23.6679
      
21.5568
      
19.8954
      
18.7935
      
18.3071
      
17.671
      
17.298
      
17.1647
    
  

In [ ]: