From ef04f91a93690e5e353b4a6c9e133806e02cbcef Mon Sep 17 00:00:00 2001 From: elvinpoole Date: Wed, 9 Oct 2024 15:50:56 -0400 Subject: [PATCH] updated DES maglim flag cut and documented --- txpipe/lens_selector.py | 6 +++++- 1 file changed, 5 insertions(+), 1 deletion(-) diff --git a/txpipe/lens_selector.py b/txpipe/lens_selector.py index f5373b27..99b54fff 100755 --- a/txpipe/lens_selector.py +++ b/txpipe/lens_selector.py @@ -295,7 +295,11 @@ def select_lens_DESmaglim(self, phot_data): cut1 = (mag_i > bright_limit) #cut very bright gaalxies cut2 = (mag_i < a*z + b) #the z-dependant mag cut cut3 = (sg == 3) #star/galaxy separator - cut4 = (np.bitwise_and(flags,120) == 0) #flags + + # flags, selects objects flagged by bits 2^1, 2^2, 2^3, 2^4, 2^5, 2^6 but not 2^7 + # from Rodriguez-Monroy et al + cut4 = (np.bitwise_and(flags,126) == 0) + s = (cut1 & cut2 & cut3 & cut4).astype(np.int8) return s