Project

General

Profile

Galaxy depth and completeness using Mangle Mask » History » Version 18

Diego Capozzi, 11/27/2013 07:54 AM

1 1 Diego Capozzi
h1. Galaxy depth and completeness using Mangle Mask
2 1 Diego Capozzi
 
3 10 Diego Capozzi
*PAGE UNDER CONSTRUCTION!*
4 10 Diego Capozzi
5 2 Diego Capozzi
Here I (Diego) describe the results of a study aiming at identifying homogeneous depth regions using the information contained in the Mangle mask, using SVA1 data. These regions will be constructed in the molygon space, so they won't necessarily be constituted by contiguous molygons.  In addition, relations between depths estimates measured with various types of magnitude (the one used for constructing the Mangle mask is a 2"-aperture magnitude, MAG_APER4) is explored so to be useful for data selection purposes, which can be different according to the science carried out. A first attempt to study these issues in detail has been carried out and described here: [[Some tests on depth with Mangle Mask]]. Here I implement the old tests with the selection cuts that have been discussed during the SVA1 telecons and in particular among Eli, Nacho and myself. These cuts also include the latest star/galaxy separation criteria described here: [[A Modest Proposal for Preliminary Star/Galaxy Separation]].
6 1 Diego Capozzi
7 2 Diego Capozzi
Galaxy completeness cannot be studied at the moment, but I use the surface brightness information to identify what could possibly be the magnitude at which the survey is complete. However, without having a deeper reference for a standard 10=tilings field, at which completeness level this magnitude value corresponds to cannot be inferred.
8 1 Diego Capozzi
9 1 Diego Capozzi
10 2 Diego Capozzi
I point out that the study described below is carried out only on the following fields: SPT (E & W); El Gordo; RXJ2248. As of today, there is no SVA1 data-based Mangle Mask for tew Bullet Cluster field. SN fields (as of now, a mangle mask is available only for SN_S) have been excluded on purpose, as my aim is to characterise a region with homogeneous depth for the standard 10-tilings Survey (from now on referred to as Wide Survey). In addition, I only considered tiles where the zero-point offsets are minimised according to Huan Lin's "galaxy-locus criteria":.
11 1 Diego Capozzi
12 2 Diego Capozzi
Finally, molygons with area>3 str and with i-band 2"-aperture mag=0 (the latter corresponding to regions that are masked by bright stars) have been excluded (see [[Some tests on depth with Mangle Mask]] for more details).
13 1 Diego Capozzi
14 1 Diego Capozzi
h2. 1. Studying galaxy depth via galaxy counts, using Mangle mask: comparing MAG_AUTO, MAG_DETMODEL & MAG_APER4 depths against molygon ones. 
15 9 Diego Capozzi
16 1 Diego Capozzi
17 1 Diego Capozzi
* Data Selection
18 1 Diego Capozzi
** Galaxy selection step 1 (from SVA1_COADD_GRIZY table):  
19 1 Diego Capozzi
### Contained in SPT (E & W), El Gordo & RXJ2248 fields
20 1 Diego Capozzi
### Contained in good tiles according to Huan Lin's "galaxy locus" zeropoint offsets criteria 
21 2 Diego Capozzi
### MAGERR_APER4_I<0.11
22 2 Diego Capozzi
### -1<MAG_DETMODEL_G-MAG_DETMODEL_R<4
23 2 Diego Capozzi
### -1<MAG_DETMODEL_I-MAG_DETMODEL_Z<4
24 2 Diego Capozzi
### NITER_MODEL_GRIZ>0
25 2 Diego Capozzi
### DEC>-61 (avoid LMC)
26 2 Diego Capozzi
### FLAGS_I<4
27 2 Diego Capozzi
### (((SPREAD_MODEL_I+(3*SPREADERR_MODEL_I))>0.003) and ((MAG_AUTO_I>12 and MAG_AUTO_I<18 and CLASS_STAR_I<0.3) or (MAG_AUTO_I>18 and MAG_AUTO_I<25))) and MAG_PSF<30
28 1 Diego Capozzi
** Molygon selection step 1 (from MOLYGON table):
29 1 Diego Capozzi
### I-band molygons contained in SPT (E & W), El Gordo & RXJ2248 fields
30 1 Diego Capozzi
### I-band molygons in tiles according to Huan Lin's "galaxy locus" zeropoint offsets criteria 
31 1 Diego Capozzi
### MAG_LIMIT (i-band)>0 (Avoiding regions masked by bright stars)
32 1 Diego Capozzi
### AREA_STR<3 [Avoiding molygons with strangely large area (these are few outliers and reducing the cut to 1 or less won't change the extracted molygons)]
33 1 Diego Capozzi
** Galaxy selection step 2 (from COADD_OBJECTS_MOLYGON table): 
34 1 Diego Capozzi
### Contained in molygons selected in previous point (step necessary to extract information linking galaxies to the molygons they belong to)
35 1 Diego Capozzi
** Galaxy selection step 3:
36 1 Diego Capozzi
### Selecting properties from SVA1_COADD_GRIZY of galaxies identified in step 1 which are contained in the identified molygons using the information obtained in "Galaxy selection step 2"
37 1 Diego Capozzi
** Molygon selection step 2:
38 2 Diego Capozzi
### Checking the actual molygons used to identify the unique molygons ID values associated with the galaxy sample identified in "Galaxy selection step 3"   
39 1 Diego Capozzi
40 2 Diego Capozzi
* For the identified galaxy sample, values of magnitude (APER_MAG4 excluded, for which corrections are not available) and surface brightness have been corrected for zeropoint offset according to Huan Lin's tables.
41 1 Diego Capozzi
42 1 Diego Capozzi
43 3 Diego Capozzi
* The importance of what surface brightness measurement is used.
44 3 Diego Capozzi
In the previous study, the surface brightness measure that was used was model surface brightness [either effective (MU_EFF_MODEL) or that at the brightness peak (MU_MAX_MODEL)]. This measurement was always found to show a gausian-like distribution rather than the common steeply increasing distribution with a sharp drop at the faint end. This was proven to be independent on the depth inhomogeneity over the sky (see GEWG report on SVA1 for details). For this study I used MU_MAX, which is a model-independent measure of surface brightness. In the plot in Figure 1 below, I compare the distributions of MU_EFF_MODEL (whose shape is equivalent to the one of MU_MAX_MODEL) and MU_MAX for the galaxy catalogue identified above. One can see the great difference between the two distributions, and that the MU_MAX distribution is definitely closer to the expected distribution, despite its drop still lacking the expected sharpness. 
45 3 Diego Capozzi
46 3 Diego Capozzi
47 1 Diego Capozzi
!{width:900px}surf_bright_hist_comparison.gif!
48 15 Diego Capozzi
*Figure 1*
49 3 Diego Capozzi
50 3 Diego Capozzi
51 3 Diego Capozzi
52 5 Diego Capozzi
This comparison indicates that there is something going on with model-based measures of surface brightness. I point out that DC6 data did not show such a gaussian-like shape of model-dependent surface brightness measures, as shown in Figure 2 below. 
53 1 Diego Capozzi
54 1 Diego Capozzi
55 4 Diego Capozzi
!{width:600px}Fig4b_dc6_rep.jpg!
56 15 Diego Capozzi
*Figure2:* Distribution of model-dependent surface brightness measure. See black full line.
57 2 Diego Capozzi
58 1 Diego Capozzi
* I-band depth comparison. What has been done:
59 1 Diego Capozzi
** Binning molygons in mag_limit (2-arcsec aperture-magnitude depth associated with molygons as output of the Mangle masking process) bins of 0.2 mag. 
60 2 Diego Capozzi
** In each bin, distributions of MAG_APER4, MU_MAX, MAG_AUTO and MAG_DETMODEL are analysed. Bin size used for all these distributions: 0.05 mag
61 1 Diego Capozzi
** For each distribution the depth is found as the peak in the number counts [if the peak is not unique, but more than 1 with the same value are seen, then the first peak (corresponding to the brightest magnituded) is taken]. The presence of more peaks should only happen with smaller galaxy samples, i.e. for the brighter molygon-depth bins. In general, the peak should be one and at magnitude fainter than this peak, galaxy number counts should decrease.
62 1 Diego Capozzi
** Galaxy depth is measured. In general, the process for defining galaxy depth for galaxies in molygons within a given mag_limit bin, works as follows:
63 2 Diego Capozzi
### Find peak in the magnitude distribution. This peak is defined as the galaxy magnitude depth. 
64 2 Diego Capozzi
65 2 Diego Capozzi
** Approximate value for galaxy magnitude limit at which the sample is complete at an unknown completeness level is measured. In general, the process for defining this value for galaxies in molygons within a given mag_limit bin, works as follows:
66 2 Diego Capozzi
### Find peak in the surface brigthness (in this case MU_MAX) distribution. NOTE: the shape of this distribution does not show the typical sharp drop after the peak. The situation is definitely better than when using MU_MAX_MODEL (see below for more details).
67 1 Diego Capozzi
### Select galaxies with: reasonable lower limit<surface brightness values<surface brightness peak. In this study,  reasonable lower limit=16.
68 2 Diego Capozzi
69 2 Diego Capozzi
70 7 Diego Capozzi
Note that, as expected, the surface brightness strongly influences the galaxy completeness limit.
71 2 Diego Capozzi
72 8 Diego Capozzi
The following plots show the results obtained for molygons with 22.6<MAG_LIMIT<24.2.
73 2 Diego Capozzi
74 1 Diego Capozzi
75 8 Diego Capozzi
*PLOT CAPTION*
76 8 Diego Capozzi
upper left panel: 2-arcsec aperture-magnitude distribution for the identified galaxy sample within the selected molygons
77 8 Diego Capozzi
upper centre panel: MU_MAX distributions for the identified galaxies in: a) all the molygons selected in the "Data Selection" section; b) molygons in the the MAG_LIMIT bin analysed
78 8 Diego Capozzi
upper right panel: MAG_AUTO distributions for the identified galaxies in the molygons contained in the analysed MAG_LIMIT bin: with and without MU_MAX selection
79 8 Diego Capozzi
lower left panel: MAG_DETMODEL distributions for the identified galaxies in the molygons contained in the analysed MAG_LIMIT bin: with and without MU_MAX selection
80 8 Diego Capozzi
lower centre panel: MAG_APER4 distributions for the identified galaxies in the molygons contained in the analysed MAG_LIMIT bin: with and without MU_MAX selection
81 1 Diego Capozzi
lower right panel: comparison between MAG_AUTO, MAG_DETMODEL & MAG_APER4 distribution for galaxies selected also in MU_MAX
82 9 Diego Capozzi
83 9 Diego Capozzi
TO BE UPDATED!
84 9 Diego Capozzi
85 9 Diego Capozzi
+22.6<MAG_LIMIT<22.8+
86 11 Diego Capozzi
!{width:900px}Depth_comparison_magerr_aper4_selection_bin3.jpg!
87 15 Diego Capozzi
*Figure 3*
88 9 Diego Capozzi
89 1 Diego Capozzi
+22.8<MAG_LIMIT<23.0+
90 1 Diego Capozzi
!{width:900px}Depth_comparison_magerr_aper4_selection_bin4.jpg!
91 15 Diego Capozzi
*Figure 4*
92 1 Diego Capozzi
93 1 Diego Capozzi
+23.0<MAG_LIMIT<23.2+
94 11 Diego Capozzi
!{width:900px}Depth_comparison_magerr_aper4_selection_bin5.jpg!
95 15 Diego Capozzi
*Figure 5*
96 1 Diego Capozzi
97 9 Diego Capozzi
+23.2<MAG_LIMIT<23.4+
98 11 Diego Capozzi
!{width:900px}Depth_comparison_magerr_aper4_selection_bin6.jpg!
99 15 Diego Capozzi
*Figure 6*
100 9 Diego Capozzi
101 9 Diego Capozzi
+23.4<MAG_LIMIT<23.6+
102 11 Diego Capozzi
!{width:900px}Depth_comparison_magerr_aper4_selection_bin7.jpg!
103 15 Diego Capozzi
*Figure 7*
104 1 Diego Capozzi
105 1 Diego Capozzi
+23.6<MAG_LIMIT<23.8+
106 11 Diego Capozzi
!{width:900px}Depth_comparison_magerr_aper4_selection_bin8.jpg!
107 15 Diego Capozzi
*Figure 8* 
108 1 Diego Capozzi
109 9 Diego Capozzi
+23.8<MAG_LIMIT<24.0+ 
110 1 Diego Capozzi
!{width:900px}Depth_comparison_magerr_aper4_selection_bin9.jpg!
111 15 Diego Capozzi
*Figure 9*
112 1 Diego Capozzi
113 9 Diego Capozzi
+24.0<MAG_LIMIT<24.2+
114 1 Diego Capozzi
!{width:900px}Depth_comparison_magerr_aper4_selection_bin10.jpg!
115 15 Diego Capozzi
*Figure 10*
116 1 Diego Capozzi
117 1 Diego Capozzi
*Galaxy depth comparison:*
118 15 Diego Capozzi
The galaxy depths measured via number counts in each MAG_LIMIT bin are compared with those given as the output of the mangle masking process in Figure 11 (left panel). The same is done for galaxy completeness limits (i.e., magnitude values at which the sample is complete at an unknown completeness level, calculated using also surface brightness cuts) in the right panel of Figure 11. I point out that the "Reliable" region inFigure 11 plots represents MAG_LIMIT bins where the galaxy sample is large enough not to show multiple peaks in the magnitude distributions.
119 1 Diego Capozzi
120 1 Diego Capozzi
Note that, as expected, the galaxy completeness values are brighter (~0.5 mag brighter) than depth ones, a results that also shows how surface brightness cuts influences the completeness of a catalogue. Another result to note is that, differently form the results of the previous study, the number-counts derived depth follows very closely the one given by the Mangle mask (based on finding the 10-sigma cut a  magnitude vs magnitude error plot). The main reason for this agreement (which should be anyway expected and is actually reassuring, given that two different methods for estimating depth are consistent with each other, as they should) is the use of magerr_aper4<0.11 in the catalogue selection. This is a natural consequence of the fact that the way the Mangle mask depth is calculated makes use of these mag_aper4 errors. 
121 1 Diego Capozzi
122 1 Diego Capozzi
!{width:400px}Depths_magerr_aper4_selection.jpg! !{width:400px}Completeness_magerr_aper4_selection.jpg!
123 1 Diego Capozzi
*Figure 11:* Case of MAGERR_APER4<0.11 used in sample selection. _Left panel_: Mangle Mask depth vs. number-counts depth; _right panel_: Mangle Mask depth vs. number-count vs. number-counts completeness
124 1 Diego Capozzi
125 1 Diego Capozzi
126 15 Diego Capozzi
The improved star/galaxy separation does not make a significant difference. In fact, when repeating the analysis using a catalogue selected as described in the _Data Selection_ Section but using magerr_auto<0.11 instead of magerr_aper4<0.11, one gets the same results as in the old analysis, i.e. number-counts depth values for APER_MAG4 are significantly different than those given by the Mangle mask, as shown in the left panel of Figure 12 (the distribution plots for the magerr_auto<0.11 are shown in the _Additional Plots_ Section). The improved star/galaxy separation has probably more effect for lower detection level (higher magnitude errors). [Eli, did you cut your sample to a 10-sigma level when doing your latest tests on star/galaxy separation?] 
127 1 Diego Capozzi
128 1 Diego Capozzi
!{width:400px}Depths_magerr_auto_selection.jpg! !{width:400px}Completeness_magerr_auto_selection.jpg!
129 1 Diego Capozzi
*Figure 12:* Case of MAGERR_AUTO<0.11 used in sample selection. _Left panel_: Mangle Mask depth vs. number-counts depth; _right panel_: Mangle Mask depth vs. number-count vs. number-counts completeness
130 15 Diego Capozzi
131 15 Diego Capozzi
132 15 Diego Capozzi
As expected, in all cases, the integrated-magnitude depth and completeness, as measured via MAG_DETMODEL or MAG_AUTO, are always ~0.5 mag shallower than the aperture-magnitude ones. It is also good news that MAG_DETMODEL and MAG_AUTO depths are always close to each other. 
133 15 Diego Capozzi
134 15 Diego Capozzi
It is important to notice that one can use the Mangle mask depth only when selecting a galaxy sample similarly to what done for the Mangle mask process (e.g., using MAGERR_APER4). If for any scientific purposes the galaxy sample selection is carried out differently, the Mangle mask cannot be used or can only be used by applying corrections (e.g., magnitude scaling corrections), even when using the same magnitudes (e.g., MAG_APER4), as Figure 12 demonstrates. For instance, if for studying the galaxy luminosity function one wants to select a galaxy catalogue out to the 10-sigma level depth for MAG_AUTO (or MAG_DETMODEL), one cannot use mangle mask depth values only corrected for flux differences between aperture and integrated magnitude, to trace the 10-sigma level depth for MAG_AUTO (or MAG_DETMODEL). A correction taking also into account that the 10-sigma level for the integrated magnitude is different from the one corresponding to the aperture magnitude. This fact is clear if one compares left panels of Figures 11 and 12. For instance, the sample selected as described in the _Data Selection_ Section (so using MAGERR_APER4<0.11), contains about 2 millions galaxies more than the one selected in the same way but using the MAGERR_AUTO<0.11 criterion. This difference in galaxy selection causes the difference in depth and completeness values (see Figures 11 and 12). 
135 15 Diego Capozzi
136 15 Diego Capozzi
We can probably model these magnitude differences so to calibrate them in the "Reliable" region identified in the plots in Figures 11 and 12.
137 15 Diego Capozzi
138 15 Diego Capozzi
h3. Additional plots (the MAGERR_AUTO<0.11 case)
139 15 Diego Capozzi
140 15 Diego Capozzi
+22.6<MAG_LIMIT<22.8+
141 15 Diego Capozzi
!{width:900px}Depth_comparison_magerr_auto_selection_bin3.jpg!
142 15 Diego Capozzi
*Figure 13*
143 15 Diego Capozzi
144 15 Diego Capozzi
+22.8<MAG_LIMIT<23.0+
145 15 Diego Capozzi
!{width:900px}Depth_comparison_magerr_aauto_selection_bin4.jpg!
146 15 Diego Capozzi
*Figure 14*
147 15 Diego Capozzi
148 15 Diego Capozzi
+23.0<MAG_LIMIT<23.2+
149 15 Diego Capozzi
!{width:900px}Depth_comparison_magerr_auto_selection_bin5.jpg!
150 15 Diego Capozzi
*Figure 15*
151 15 Diego Capozzi
152 15 Diego Capozzi
+23.2<MAG_LIMIT<23.4+
153 15 Diego Capozzi
!{width:900px}Depth_comparison_magerr_auto_selection_bin6.jpg!
154 15 Diego Capozzi
*Figure 16*
155 15 Diego Capozzi
156 15 Diego Capozzi
+23.4<MAG_LIMIT<23.6+
157 15 Diego Capozzi
!{width:900px}Depth_comparison_magerr_auto_selection_bin7.jpg!
158 15 Diego Capozzi
*Figure 17*
159 15 Diego Capozzi
160 15 Diego Capozzi
+23.6<MAG_LIMIT<23.8+
161 15 Diego Capozzi
!{width:900px}Depth_comparison_magerr_auto_selection_bin8.jpg!
162 15 Diego Capozzi
*Figure 18* 
163 15 Diego Capozzi
164 15 Diego Capozzi
+23.8<MAG_LIMIT<24.0+ 
165 15 Diego Capozzi
!{width:900px}Depth_comparison_magerr_auto_selection_bin9.jpg!
166 15 Diego Capozzi
*Figure 19*
167 15 Diego Capozzi
168 15 Diego Capozzi
+24.0<MAG_LIMIT<24.2+
169 15 Diego Capozzi
!{width:900px}Depth_comparison_magerr_auto_selection_bin10.jpg!
170 15 Diego Capozzi
*Figure 20*
171 15 Diego Capozzi
172 15 Diego Capozzi
h2. Final Remarks
173 15 Diego Capozzi
174 15 Diego Capozzi
On average, molygons correctly tell us how the aperture magnitude depth varies over the surveyed area. However there are still suspicious molygons with too low a total exposure time for the depth value they are associated with. Below I report the tables obtained in the old analysis reporting the expected stellar magnitude in the i-band at several S/N levels (what we want to consider though is S/N~10 level) for two values of total exposure time. These are calculated by using the DECam Exposure Time Calculator for stellar photometry available on the "CTIO DECam webpage":http://www.ctio.noao.edu/noao/content/Dark-Energy-Camera-DECam. I have considered the case of best observing conditions, i.e. when the sky brightness is the lowest possible (so to speak, with new moon). This choice is to show that even in the best observing scenario, the depth reached in some molygons doesn't seem realisitic. One has to take into account that the Exposure Calculator is made for stellar photometry, so one should roughly subtract 0.5/1 mag from the stellar magnitude value it provides, to get a realistic value for galaxies.
175 15 Diego Capozzi
176 15 Diego Capozzi
h3. Case of Total exposure time fixed at 540 s (i.e. 6 x pointings/90 s exposures)
177 15 Diego Capozzi
178 15 Diego Capozzi
|_.Stellar i-band magnitude limit (mag)      |_.Galaxy  i-band magnitude limit (mag)                  |_.S/N|
179 15 Diego Capozzi
|26                                                            |~25/25.5                                                                |1.9|
180 15 Diego Capozzi
|25.5                                                         |~24.5/25                                                                |3.1|
181 15 Diego Capozzi
|25                                                            |~24/24.5                                                                |4.9|
182 15 Diego Capozzi
|24.5                                                         |~23.5/24                                                                |7.7| 
183 15 Diego Capozzi
|24                                                            |~23/23.5                                                                |12.1| 
184 1 Diego Capozzi
|23.5                                                         |~22.5/23                                                                |19.1|
185 15 Diego Capozzi
|23                                                            |~22/22.5                                                                |30.0|
186 17 Diego Capozzi
187 16 Diego Capozzi
*Table 1*
188 15 Diego Capozzi
189 15 Diego Capozzi
h3. Case of Total exposure time fixed at 450 s (i.e. 5 x pointings/90 s exposures)
190 15 Diego Capozzi
191 15 Diego Capozzi
|_.Stellar i-band magnitude limit (mag)      |_.Galaxy  i-band magnitude limit (mag)                  |_.S/N|
192 15 Diego Capozzi
|26                                                            |~25/25.5                                                                |1.8|
193 15 Diego Capozzi
|25.5                                                         |~24.5/25                                                                |2.8|
194 15 Diego Capozzi
|25                                                            |~24/24.5                                                                |4.4|
195 15 Diego Capozzi
|24.5                                                         |~23.5/24                                                                |7.0| 
196 15 Diego Capozzi
|24                                                            |~23/23.5                                                                |11.1| 
197 1 Diego Capozzi
|23.5                                                         |~22.5/23                                                                |17.4|
198 15 Diego Capozzi
|23                                                            |~22/22.5                                                                |27.4|
199 17 Diego Capozzi
200 16 Diego Capozzi
*Table 2*
201 15 Diego Capozzi
202 15 Diego Capozzi
I also report some plots (Figures 21, 22 and 23) summarising the properties of the molygons identified in the _Data Selection_ Section (MAGERR_APER4<0.11 case), which can be inspected having the results of Tables 1 and 2 in mind.
203 15 Diego Capozzi
204 15 Diego Capozzi
h2. Next steps (to do list)
205 15 Diego Capozzi
206 15 Diego Capozzi
* Identify 3-4 different regions in the used footprint constituted by molygons with depth deeper than a minimum depth value and totalling, for each region, a reasonable area (at least ~100 sq. deg)
207 15 Diego Capozzi
* For each region, select galaxies out to the identified minimum depth value, so to homogenise depth over the identified region
208 15 Diego Capozzi
* Study Luminosity and Stellar Mass Functions (LF & MF) for the samples identified in the 3-4 identified regions
209 15 Diego Capozzi
* Study completeness as soon as this becomes possible and implement this into the study of LF and MF
210 15 Diego Capozzi
* Rescale MAG_DETMODEL values to i-band MAG_AUTO
211 15 Diego Capozzi
* Repeat all the analysis (including depth determination)