# sex_caller('lhn1n1_2010apr_r_stack_fc_fix.fits', 'lhn1n1_2010apr_r_stack_fc_fix')
# sex_caller('lhn1n1_2010dec_z_stack_fc_fix.fits', 'lhn1n1_2010dec_z_stack_fc_fix')
# Read catalogs
catag_r = CatalogReader("lhn1n1_2010apr_r_stack_fc_fix.cat")
catag_r.read()
catag_z = CatalogReader("lhn1n1_2010dec_z_stack_fc_fix.cat")
catag_z.read()
# Give value to the cross-matching radius
r = 3
# Create object for cross-matching
crossmatching = CrossMatching(catag_r.fcat, catag_z.fcat)
crossmatching.kdtree(n=r * 1e-06)
crossmatching.catalog_writter("lhn1n1_crossmatching_1to2", compare="1to2")
print "\n"
crossmatching.catalog_writter("lhn1n1_crossmatching_2to1", compare="2to1")
if crossmatching.cont1to2 < crossmatching.cont2to1:
catag_final_1 = CatalogReader("lhn1n1_crossmatching_1to2.fcat")
catag_final_1.read()
catag_final_2 = CatalogReader("lhn1n1_crossmatching_2to1.fcat")
catag_final_2.read()
crossmatching_final = CrossMatching(catag_final_1.fcat, catag_final_2.fcat)
crossmatching_final.kdtree(n=r * 1e-06)
crossmatching_final.catalog_writter("lhn1n1_crossmatching_final", compare="2to1")
if crossmatching.cont1to2 > crossmatching.cont2to1:
