From c088a8780979a3d8b91fd3d22ad9ef116c2c4e8a Mon Sep 17 00:00:00 2001
From: "P.M. Sutter" <sutter@iap.fr>
Date: Thu, 25 Apr 2013 11:05:12 -0500
Subject: [PATCH] allow multiple stacks/slices for AP analysis within a
 redshift bin

---
 python_tools/pipeline_source/prepareCatalogs.in.py | 12 +++++++++---
 1 file changed, 9 insertions(+), 3 deletions(-)

diff --git a/python_tools/pipeline_source/prepareCatalogs.in.py b/python_tools/pipeline_source/prepareCatalogs.in.py
index fbb0c36..5bcaa7e 100644
--- a/python_tools/pipeline_source/prepareCatalogs.in.py
+++ b/python_tools/pipeline_source/prepareCatalogs.in.py
@@ -295,9 +295,15 @@ newSample.addStack({zMin}, {zMax}, 2*{minRadius}+18, 2*{minRadius}+24, True, Fal
                                          useLightCone=useLightCone,
                                          subsample=str(subsample).strip('[]')))
 
-          scriptFile.write(stackInfo.format(zMin=sliceMin,
-                                           zMax=sliceMax,
-                                           minRadius=minRadius))
+          for iAPSlice in xrange(numAPSlices):
+            sliceWidth = float(sliceMax) - float(sliceMin)
+            sliceAPMin = float(sliceMin) + iAPSlice*sliceWidth/numAPSlices
+            sliceAPMax = float(sliceMin) + (iAPSlice+1)*sliceWidth/numAPSlices
+            sliceAPMin = "%0.2f" % sliceAPMin
+            sliceAPMax = "%0.2f" % sliceAPMax
+            scriptFile.write(stackInfo.format(zMin=sliceAPMin,
+                                              zMax=sliceAPMax,
+                                              minRadius=minRadius))
 
 
   scriptFile.close()