#!/usr/bin/python

import h5py

#---------------------------------
# Read GLAM-LRG catalogue
#--------------------------------

hf = h5py.File('z0.60/DESI_GLAM_LRG_z0.60_HOD_0.hdf5', 'r') # See elements hf.keys()

# Example: read postions and velocities

pos = hf.get('pos').value
vel = hf.get('vel').value

# Example: read correlation functions

s = hf.get('Correlation_function').get('s').value
xi_r =  hf.get('Correlation_function').get('xi_r').value
xi_0 = hf.get('Correlation_function').get('xi_0').value
xi_2 = hf.get('Correlation_function').get('xi_2').value