if __name__=="__main__": nstart = args.input nend = args.output dataDDnf = correlation.openDataFront('data/galaxies_DR9_CMASS_North.fits', nstart, nend) dataRRnf = correlation.openDataFront('data/randoms_DR9_CMASS_North.fits', nstart, nend) dataDDnb = correlation.openDataBack('data/galaxies_DR9_CMASS_North.fits', nstart, nend) dataRRnb = correlation.openDataBack('data/randoms_DR9_CMASS_North.fits', nstart, nend) # dataDDs = correlation.openData('../Downloads/galaxies_DR9_CMASS_South.fits', npoints) # dataRRs = correlation.openData('../Downloads/randoms_DR9_CMASS_South.fits', npoints) print "Concat Now" DpRnf = dataDDnf + dataRRnf DpRnb = dataDDnb + dataRRnb # DpRs = dataDDs + dataRRs print "First hist" histDDnf = correlation.histBB(dataDDnf) histDpRnf = correlation.histBB(DpRnf) histRRnf = correlation.histBB(dataRRnf) histDDnb = correlation.histBB(dataDDnb) histDpRnb = correlation.histBB(DpRnb) histRRnb = correlation.histBB(dataRRnb) histDDn = histDDnf.Clone("histDDn") histDDn.Reset() histDDn.Add(histDDnf) histDDn.Add(histDDnb) histDpRn = histDDnf.Clone("histDDn") histDpRn.Reset() histDpRn.Add(histDpRnf) histDpRn.Add(histDpRnb) histRRn = histDDnf.Clone("histDDn")
def renameData(table): return [correlation.cartesian(row) for row in table] if __name__=="__main__": td1 = cuts(resize(openTable('data/galaxies_DR9_CMASS_North.fits')), 205, 207, 0, 2) tr1 = cuts(openTable('data/randoms_DR9_CMASS_North.fits'), 205, 207, 0, 2) tdpr1 = vstack([td1, tr1]) td = renameData(td1) tdpr = renameData(tdpr1) tr = renameData(tr1) print "Concat" print len(td) print len(tr) print len(tdpr) histDD = correlation.histBB(td) histDpR = correlation.histBB(tdpr) histRR = correlation.histBB(tr) print "correlation done" histD = histDD.Clone("histD") histD.Reset() histD.Add(histDD) histDR = histDD.Clone("histD") histDR.Reset() histDR.Add(histDpR) histR = histDD.Clone("histD") histR.Reset() histR.Add(histRR) corrhist = correlationHistogram(histD, histDR, histR)