guilhem_lavaux/vide_public
0
0
Fork 0
mirror of https://bitbucket.org/cosmicvoids/vide_public.git synced 2025-07-04 15:21:11 +00:00
Code Issues Projects Releases Packages Wiki Activity Actions

corrected bug in reporting of RA and Dec

Browse source
This commit is contained in:
P.M. Sutter 2015-09-06 15:26:32 -05:00
parent 3bb09067cd
commit c9cb9ea1cc
1 changed files with 15 additions and 5 deletions
Download patch file Download diff file Expand all files Collapse all files

20
c_tools/stacking/pruneVoids.cpp
View file

@ -921,6 +921,7 @@ int main(int argc, char **argv) {
} }
voids.resize(iGood); voids.resize(iGood);
/*
iGood = 0; iGood = 0;
for (iVoid = 0; iVoid < voids.size(); iVoid++) { for (iVoid = 0; iVoid < voids.size(); iVoid++) {
// just in case // just in case
@ -932,6 +933,8 @@ int main(int argc, char **argv) {
} }
} }
voids.resize(iGood); voids.resize(iGood);
*/
printf(" 4th filter: rejected %d outside redshift boundaries\n", numNearZ); printf(" 4th filter: rejected %d outside redshift boundaries\n", numNearZ);
// take only top-level voids // take only top-level voids
@ -1167,12 +1170,19 @@ void outputVoids(string outputDir, string sampleName, string prefix,
outVoid.level, outVoid.level,
outVoid.numChildren, outVoid.numChildren,
outVoid.centralDen); outVoid.centralDen);
double phi = atan2(outVoid.macrocenter[1]-boxLen[1]/2.,
outVoid.macrocenter[0]-boxLen[0]/2.);
if (phi < 0) phi += 2.*M_PI;
double RA = phi * 180./M_PI;
double theta = acos((outVoid.macrocenter[2]-boxLen[2]/2.) /
outVoid.redshiftInMpc);
double dec = (M_PI/2. - theta) * 180./M_PI;
fprintf(fpSkyPositions, "%.2f %.2f %.5f %.2f %d\n", fprintf(fpSkyPositions, "%.2f %.2f %.5f %.2f %d\n",
atan((outVoid.macrocenter[1]-boxLen[1]/2.) / RA,
(outVoid.macrocenter[0]-boxLen[0]/2.)) * 180/M_PI + 180, dec,
asin((outVoid.macrocenter[2]-boxLen[2]/2.) /
outVoid.redshiftInMpc) * 180/M_PI,
outVoid.redshift, outVoid.redshift,
outVoid.radius, outVoid.radius,
outVoid.voidID); outVoid.voidID);