From 85e0f36c99b7ecee42758e62412478b7060dbedf Mon Sep 17 00:00:00 2001 From: Guilhem Lavaux Date: Tue, 16 Feb 2021 12:56:00 +0200 Subject: [PATCH] Allow mask of any resolution --- c_tools/mock/generateFromCatalog.cpp | 22 +++++++++------------- 1 file changed, 9 insertions(+), 13 deletions(-) diff --git a/c_tools/mock/generateFromCatalog.cpp b/c_tools/mock/generateFromCatalog.cpp index b98a18c..4d29926 100644 --- a/c_tools/mock/generateFromCatalog.cpp +++ b/c_tools/mock/generateFromCatalog.cpp @@ -377,10 +377,10 @@ void generateSurfaceMask(generateFromCatalog_info& args , output_data.uniqueID.push_back(-1); insertion++; - fprintf(fp, "%e %e %e\n", - (p.xyz[0]), - (p.xyz[1]), - (p.xyz[2])); + // fprintf(fp, "%e %e %e\n", + // (p.xyz[0]), + // (p.xyz[1]), + // (p.xyz[2])); } } } @@ -445,10 +445,10 @@ void generateSurfaceMask(generateFromCatalog_info& args , output_data.redshift.push_back(-1); output_data.uniqueID.push_back(-1); insertion++; - fprintf(fp, "%e %e %e\n", - (p.xyz[0]), - (p.xyz[1]), - (p.xyz[2])); +// fprintf(fp, "%e %e %e\n", +// (p.xyz[0]), +// (p.xyz[1]), +// (p.xyz[2])); } fclose(fp); // END PMS @@ -606,7 +606,6 @@ int main(int argc, char **argv) cout << "Loading data " << args_info.catalog_arg << "..." << endl; vector data; - Healpix_Map o_mask; vector pixel_list; vector full_mask_list; ParticleData output_data; @@ -615,12 +614,9 @@ int main(int argc, char **argv) cout << "Loading mask..." << endl; - read_Healpix_map_from_fits(args_info.mask_arg, o_mask); - Healpix_Map mask; - mask.SetNside(128, RING); - mask.Import(o_mask); + read_Healpix_map_from_fits(args_info.mask_arg, mask); computeContourPixels(mask,pixel_list); computeMaskPixels(mask,full_mask_list);