from glob import glob
import numpy as np
from gen_plot import plotz

sph_files = glob("../1mil_hgg/volfrac*")

sph_count = []

for f in sph_files:
    sph_count.append(np.loadtxt(f))


plotz(masslist, deltarrs, style="hist_quant", axdim=[1, 3], xlog=True,
      xbins=len(sph_files), ybins=10, ymin=0.0, ymax=1.0
      name="/home/kilian/public_html/tinker/quantile_mass-delta.png")

import h5py as h5
from gen_plot import plotz


"""
This script plots the points in the randoms mask and in the veto mask
together in RA and Dec for the survey and cap given in the call
arguments.
"""


if len(sys.argv) != 2:

    print "usage: python plot_masks.py <[CMASS, LOWZ]/[NGC, SGC]>"

dirname = "../dat/" + sys.argv[1]

fgood = h5.File(dirname + "/srch_pts.hdf5")
fbad = h5.File(dirname + "/veto.hdf5")

good_pts = fgood["good_pts"]
bad_pts = fbad["bad_pts"]

plotz([good_pts[0, :], bad_pts[0, :]], [bad_pts[1, :], bad_pts[1, :]],
      style='scatter', overplot=True, colours=["g", "r"], alpha=0.1,
      name="./test_figs/{0}/masks.png".format(dirname))


fgood.close()
fbad.close()