diff --git a/mytools/generateMock.cpp b/mytools/generateMock.cpp index 2d36ed8..5c56e45 100644 --- a/mytools/generateMock.cpp +++ b/mytools/generateMock.cpp @@ -100,7 +100,7 @@ Interpolate make_cosmological_redshift(double OM, double OL, double z0, double z return buildFromVector(pairs); } -void metricTransform(SimuData *data, int axis, bool reshift) +void metricTransform(SimuData *data, int axis, bool reshift, bool pecvel) { int x0, x1, x2; @@ -144,7 +144,8 @@ void metricTransform(SimuData *data, int axis, bool reshift) // Distorted redshift z = (z_vs_D.compute(reduced_red)-z_base)*LIGHT_SPEED/100.; // Add peculiar velocity - z += v/100; + if (pecvel) + z += v/100; } } @@ -321,7 +322,7 @@ int main(int argc, char **argv) cout << "Omega_M = " << simu->Omega_M << endl; cout << "Omega_Lambda = " << simu->Omega_Lambda << endl; - metricTransform(simu, args_info.axis_arg, args_info.preReShift_flag); + metricTransform(simu, args_info.axis_arg, args_info.preReShift_flag, args_info.peculiarVelocities_flag); makeBox(simu, simuOut, args_info); delete simu; diff --git a/mytools/generateMock.ggo b/mytools/generateMock.ggo index a581ef2..d71114a 100644 --- a/mytools/generateMock.ggo +++ b/mytools/generateMock.ggo @@ -19,4 +19,5 @@ option "rangeY_max" - "Maximum range in Y for making the box" double require option "rangeZ_min" - "Minimum range in Z for making the box (after distortion)" double required option "rangeZ_max" - "Maximum range in Z for making the box (after distortion)" double required -option "preReShift" - "Reshift the zero of the Z axis" flag off \ No newline at end of file +option "preReShift" - "Reshift the zero of the Z axis" flag off +option "peculiarVelocities" - "Added peculiar velocities distortion" flag off \ No newline at end of file