# 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: