Project

General

Profile

Wiki » History » Version 81

Shantanu Desai, 10/23/2013 05:30 PM

1 1 Shantanu Desai
h1. Wiki
2 1 Shantanu Desai
3 2 Shantanu Desai
This can be divided into 3 steps
4 2 Shantanu Desai
5 3 Shantanu Desai
h2. Step 1 run modelfitting photometry
6 2 Shantanu Desai
7 2 Shantanu Desai
you can find the scripts at /big/users/shantanu/VLTmfp/catalogs/runmfp_b.sh (for b band and so on)
8 2 Shantanu Desai
9 2 Shantanu Desai
note that since initial catalogs will be uncalibrated , use -mag_zeropoint for b,I and z bands.
10 1 Shantanu Desai
The zeropoints for HST bands is :
11 3 Shantanu Desai
F606W 26.493mag, F814W 25.947mag.
12 3 Shantanu Desai
the sextractor param file used in the scripts is 
13 3 Shantanu Desai
/big/des_home/stable/etc/sex.param_sersic
14 3 Shantanu Desai
you can edit to remove params you don't need.
15 3 Shantanu Desai
16 4 Shantanu Desai
h2. Step 2 Produce a matched catalog across all 4 bands
17 3 Shantanu Desai
18 3 Shantanu Desai
This is the most tedious step.  Unfortunately you cannot run sextractor in dual image mode
19 3 Shantanu Desai
since coadds are of different size.
20 3 Shantanu Desai
The way I did it was to run a modified version of a code 
21 3 Shantanu Desai
used in DESDM called coadd_grabmatches. I first ran it to match R vs I band and then used
22 3 Shantanu Desai
that match catalog to match against z band and so on.
23 3 Shantanu Desai
A code to do this is at /home/moon/shantanu/desdm/devel/ImageProc/trunk/src/coadd_grabmatches_singleband.c
24 3 Shantanu Desai
There is one line in the code (below # 266) which needs to be changed depending upon number of coluns
25 3 Shantanu Desai
in input catalog.
26 3 Shantanu Desai
Once this is done you need to make a file with the matched magnitudes in all four bands
27 3 Shantanu Desai
(its possible to automate this script with some effort so that it does not need to be changed)
28 3 Shantanu Desai
29 3 Shantanu Desai
30 4 Shantanu Desai
h3. Step 3 Doing the stellar locus
31 3 Shantanu Desai
32 3 Shantanu Desai
One you have the catalogs all you need to do is to do a fit by eye or quantitative to SED fits provided by Joerg.
33 3 Shantanu Desai
Here is my 2nd order polynomial  curve fit done to these relations
34 3 Shantanu Desai
35 3 Shantanu Desai
1) b-R vs R-I
36 3 Shantanu Desai
b -R = -0.2904*x^2+1.319x +0.2941
37 3 Shantanu Desai
(where x=R-I)
38 3 Shantanu Desai
39 3 Shantanu Desai
2) R-I vs I -z
40 3 Shantanu Desai
R-I = -1.016*x^2+ 3.532*x +0.1294
41 3 Shantanu Desai
where x=I-z
42 3 Shantanu Desai
43 3 Shantanu Desai
3) b-R vs I-z
44 3 Shantanu Desai
R-I = -2.20378*x^2+ 3.96368*x +0.386878
45 1 Shantanu Desai
46 1 Shantanu Desai
47 4 Shantanu Desai
You can either use (1) & (2) or (2)  & (3) to get the zeropoints.
48 4 Shantanu Desai
49 1 Shantanu Desai
I did the plotting in ROOT, but probably it can be done in any software.
50 4 Shantanu Desai
51 4 Shantanu Desai
52 4 Shantanu Desai
h4. Step 4 Rerun sextractor with the new zeropoints
53 4 Shantanu Desai
Same as (1) but use the new zeropoints now.
54 5 Shantanu Desai
55 5 Shantanu Desai
h3. Direct calibration of I-band
56 6 Shantanu Desai
57 5 Shantanu Desai
Joerg pointed out that I band can be directly calibrated using
58 5 Shantanu Desai
59 5 Shantanu Desai
I_BESS = a0 + a1 * F814W + a2 * (F606W - F814W)
60 5 Shantanu Desai
61 5 Shantanu Desai
with
62 5 Shantanu Desai
63 5 Shantanu Desai
a0 = -1.009
64 5 Shantanu Desai
a1 =  1.067
65 5 Shantanu Desai
a2 = -0.001
66 7 Shantanu Desai
67 13 Shantanu Desai
!{width: 10%}Izeropoint.png!
68 8 Shantanu Desai
69 8 Shantanu Desai
This allows us to get direct calibration of I-zeropoint= 28.37+/- 0.02 (where error = sigma of fit/sqrt(no of points))
70 9 Shantanu Desai
71 9 Shantanu Desai
We then recalibrated (b-R) vs (R-I)
72 14 Shantanu Desai
73 15 Shantanu Desai
!{width: 10%}plotslrfit.png!
74 16 Shantanu Desai
75 16 Shantanu Desai
Also a histogram wrt stellar locus
76 16 Shantanu Desai
77 16 Shantanu Desai
!{width: 10%}slrhisto.png!
78 17 Shantanu Desai
79 40 Shantanu Desai
the zeropoint for b band = 26.88 +/- 0.02
80 40 Shantanu Desai
z band 28.02 +/- 0.012
81 18 Shantanu Desai
82 18 Shantanu Desai
83 18 Shantanu Desai
h1. Results for other clusters
84 19 Shantanu Desai
85 19 Shantanu Desai
h2. SPT-CLJ0102-4915
86 21 Shantanu Desai
87 19 Shantanu Desai
I band zero-point = 28.52 +/- 0.0137
88 20 Shantanu Desai
b band zp = 26.87 +/- 0.016
89 32 Shantanu Desai
z band zp = 28.21 +/- 0.015 (I-0.31)
90 24 Shantanu Desai
91 24 Shantanu Desai
h2. SPT-CLJ0533-5005
92 26 Shantanu Desai
93 27 Shantanu Desai
I band  = 28.52 +/- 0.0123
94 1 Shantanu Desai
b band= 26.69 +/- 0.014
95 30 Shantanu Desai
z band = (I-z+0.4) = 28.12 +/- 0.006
96 28 Shantanu Desai
97 28 Shantanu Desai
h3. SPT-CLJ0559-5249
98 28 Shantanu Desai
99 1 Shantanu Desai
I-band 28.32 +/- 0.005
100 29 Shantanu Desai
b band 27.2 +/- 0.023
101 29 Shantanu Desai
z band (I-z+0.44) 27.88 +/- 0.0147
102 33 Shantanu Desai
103 33 Shantanu Desai
h3. SPT-CLJ2331-5051
104 33 Shantanu Desai
105 33 Shantanu Desai
I-band  28.32 +/- 0.012
106 33 Shantanu Desai
b-band 27.07 +/- 0.017
107 33 Shantanu Desai
z-band  I-z+0.42   27.9 +/- 0.013
108 34 Shantanu Desai
109 35 Shantanu Desai
h4. SPT-CLJ2337-5942
110 34 Shantanu Desai
111 34 Shantanu Desai
I-band 28.38 +/- 0.0123
112 36 Shantanu Desai
b band 27.05 +/- 0.0198
113 36 Shantanu Desai
z (I-z+0.37)  28.01 +/- 0.0076
114 37 Shantanu Desai
115 37 Shantanu Desai
h4. SPT-CLJ2341-5119
116 38 Shantanu Desai
117 38 Shantanu Desai
I-band 28.38 +/- 0.0096
118 38 Shantanu Desai
 b-band           27.09 +/- 0.021
119 38 Shantanu Desai
z band (I-z+0.39) 27.99 +/- 0.013
120 38 Shantanu Desai
121 38 Shantanu Desai
h4. SPT-CLJ2342-5411
122 38 Shantanu Desai
123 38 Shantanu Desai
I-band 28.3 +/- 0.017
124 38 Shantanu Desai
b band 27.07  +/- 0.021
125 38 Shantanu Desai
z band  27.91 +/- 0.007 
126 39 Shantanu Desai
127 39 Shantanu Desai
128 39 Shantanu Desai
h4. SPT-CLJ2359-5009
129 39 Shantanu Desai
130 39 Shantanu Desai
I-band  28.41 +/- 0.0088
131 39 Shantanu Desai
b-band 26.92 +/- 0.023
132 39 Shantanu Desai
z-band  I-z+0.39     28.02 +/- 0.01
133 41 Shantanu Desai
134 41 Shantanu Desai
Jan 17 update:
135 42 Shantanu Desai
Joerg provided a new relation I_BESS - F814W = -0.052 + 0.0095 * (F606W - F814W) only valid for
136 42 Shantanu Desai
only valid for (F606W - F814W) < 1.7.
137 43 Shantanu Desai
138 43 Shantanu Desai
139 74 Shantanu Desai
h2. Jan 19
140 45 Shantanu Desai
141 45 Shantanu Desai
h2. SPT-CLJ000-5748
142 46 Shantanu Desai
143 43 Shantanu Desai
I :  27.86 +/- 0.0129
144 73 Shantanu Desai
b (b-r+2.37)  27.37 +- 0.045
145 73 Shantanu Desai
z (I-z+0.56) 27.3 +/- 0.017
146 45 Shantanu Desai
147 45 Shantanu Desai
h2. SPT-CLJ0102-4915
148 1 Shantanu Desai
149 46 Shantanu Desai
I : 27.89 +/- 0.014
150 73 Shantanu Desai
b : 27.47 +/- 0.0317
151 73 Shantanu Desai
z : (I-z+0.51)  27.38 +/- 0.022
152 47 Shantanu Desai
153 55 Shantanu Desai
h2. SPT-CLJ0205-5829
154 55 Shantanu Desai
155 56 Shantanu Desai
I 27.81 +/- 0.0155
156 73 Shantanu Desai
b 27.52 +/- 0.032
157 73 Shantanu Desai
z (I -z + 0.55)    27.26 +/ - 0.02
158 55 Shantanu Desai
159 55 Shantanu Desai
160 55 Shantanu Desai
161 47 Shantanu Desai
h2. SPT-CLJ0533-5005
162 47 Shantanu Desai
163 47 Shantanu Desai
I : 27.8 +/- 0.0196
164 73 Shantanu Desai
b  27.44  +/- 0.024 
165 73 Shantanu Desai
z  (I-z + 0.62)  27.18 +/- 0.02
166 48 Shantanu Desai
167 48 Shantanu Desai
168 50 Shantanu Desai
h2. SPT-CLJ0559-5009 
169 1 Shantanu Desai
170 50 Shantanu Desai
(same as 0533-5005)
171 50 Shantanu Desai
172 50 Shantanu Desai
I:   27.8 +/- 0.0174
173 73 Shantanu Desai
b  27.52 +/- 0.0276
174 73 Shantanu Desai
z  (I-z+0.62) 27.18 +/- 0.0223
175 51 Shantanu Desai
176 54 Shantanu Desai
h2. SPT-CLJ2331-5051
177 52 Shantanu Desai
178 52 Shantanu Desai
I : 27.83 +/- 0.0085
179 73 Shantanu Desai
b  27.54  +/- 0.0194
180 73 Shantanu Desai
z     (I-z+0.61)    27.22      +/- 0.0174
181 58 Shantanu Desai
182 77 Shantanu Desai
Inon reports some problems with the calibrated catalog. (seems fine after fixing this)
183 75 Shantanu Desai
SD: errors seem a bit larger
184 75 Shantanu Desai
185 58 Shantanu Desai
h2. SPT-CLJ2337-5942
186 58 Shantanu Desai
187 59 Shantanu Desai
I   27.81 +/- 0.008 
188 59 Shantanu Desai
b  27.565    +/- 0.0139
189 59 Shantanu Desai
z   (I-z+0.55)     27.26          +/- 0.0105
190 60 Shantanu Desai
191 60 Shantanu Desai
192 60 Shantanu Desai
h2. SPT-CLJ2341-5119
193 60 Shantanu Desai
194 60 Shantanu Desai
I  27.81+/- 0.0109
195 73 Shantanu Desai
b 27.6 +/- 0.0241
196 73 Shantanu Desai
z  (I-z + 0.565) 27.245 +/- 0.0154
197 64 Shantanu Desai
198 64 Shantanu Desai
199 64 Shantanu Desai
h2. SPT-CLJ2040-5726
200 64 Shantanu Desai
201 64 Shantanu Desai
I  27.8  +/- 0.0162
202 73 Shantanu Desai
b 27.46         +/- 0.0215
203 73 Shantanu Desai
z  (I-z + 0.5)   27.3 +/- 0.0297 (large asymmetric tail in SLR histogram) 
204 65 Shantanu Desai
205 65 Shantanu Desai
206 65 Shantanu Desai
h2. SPT-CLJ0615-5746
207 65 Shantanu Desai
208 65 Shantanu Desai
I : 27.7 +/- 0.0155
209 73 Shantanu Desai
b  27.54        +/-  0.0234 (a slight asymmetric tail)
210 73 Shantanu Desai
z (I-z+0.65) 27.05  +/- 0.01669
211 66 Shantanu Desai
212 66 Shantanu Desai
213 66 Shantanu Desai
h2. SPT-CLJ0546-5345
214 66 Shantanu Desai
215 66 Shantanu Desai
I  27.79 +/- 0.0155
216 73 Shantanu Desai
b 27.41  +/- 0.0283
217 73 Shantanu Desai
z (I-z+0.62)  27.17  +/- 0.0197
218 67 Shantanu Desai
219 67 Shantanu Desai
220 67 Shantanu Desai
h2.  SPT-CLJ2342-5411
221 67 Shantanu Desai
222 69 Shantanu Desai
I 27.79 +/- 0.011
223 73 Shantanu Desai
b  27.49  +/- 0.021
224 73 Shantanu Desai
z  (I-z+0.57)  27.22 +/- 0.0147
225 69 Shantanu Desai
226 69 Shantanu Desai
h2. SPT-CLJ2106-5844
227 70 Shantanu Desai
228 70 Shantanu Desai
I 27.74 +/- 0.0179
229 73 Shantanu Desai
b  27.49   +/- 0.023
230 73 Shantanu Desai
z  (I-z+0.49)  27.25 +/- 0.0218
231 71 Shantanu Desai
232 76 Shantanu Desai
March 2013
233 76 Shantanu Desai
zeropoints re-calculated after cataloging in dual image mode with euclid style reductions
234 76 Shantanu Desai
235 80 Shantanu Desai
I 27.74 +/- 0.007
236 80 Shantanu Desai
b 27.49 +/- 0.011
237 81 Shantanu Desai
z (I-z+1.1) 26.64 +/- 0.022 -> WRONG (see below)
238 76 Shantanu Desai
239 76 Shantanu Desai
240 81 Shantanu Desai
October 2013 Update : the z-band zeropoints from March were wrong. The new zeropoints are same as in January
241 79 Shantanu Desai
242 80 Shantanu Desai
z  (I-z+0.49)  27.25 +/- 0.0074
243 79 Shantanu Desai
244 79 Shantanu Desai
245 72 Shantanu Desai
h2. SPT-CLJ2359-5009
246 71 Shantanu Desai
247 71 Shantanu Desai
I 27.79 +/- 0.0056
248 73 Shantanu Desai
b 27.53 +/- 0.0272
249 73 Shantanu Desai
z (I-z+0.57) 27.22 +/- 0.0146
Redmine Appliance - Powered by TurnKey Linux