Context Navigation

etude_onoff.pic@ 639

Last change on this file since 639 was 556, checked in by campagne, 14 years ago
generalise the process to download the files for final analysis (jec)
File size: 10.0 KB

Line
1
2	set cols ( blue red black green purple orange cyan violet yellow )
3
4	set toppath "/sps/baoradio/AmasNancay/JEC"
5	set source "Abell${1}"
6	set srclower "abell${1}"
7	set date "${2}"
8
9
10	#####################################
11	#Plots ON-OFF with all calibrations
12	#####################################
13	openppf ${toppath}/${source}/diffOnOff_${date}_${srclower}.ppf
14
15
16	#specONOFF<i> Matrix 2 x 8192 : i=0...
17	#specONOFFMean Matrix 2 x 8192
18	#specONOFF2ChanMean Vector 8192 x 1
19
20
21	#####################################
22	#Plots with different Means: per channel and all
23	#####################################
24	del sallcyclev0
25	objaoper specONOFFMean row 0 sallcyclev0
26	del sallcyclev1
27	objaoper specONOFFMean row 1 sallcyclev1
28
29	newwin 1 1
30	graphicatt "xylimits=1250,1500,-0.5,0.5"
31	plot2d specONOFF2ChanMean (n/8192)*250+1250 val n>0 "black cpts nsta notit"
32	plot2d sallcyclev0 (n/8192)*250+1250 val n>0 "same blue cpts nsta notit"
33	plot2d sallcyclev1 (n/8192)*250+1250 val n>0 "same red cpts nsta notit"
34	settitle "ON-OFF ${source}/${date} Ch 0 (blue) Ch 1 (red) Mean (black) All cycles"
35	setaxelabels "Freq. (MHz)" "I (Jy)"
36
37	#####################################
38	#Examine fluctuations sigma for the combination of all cycles
39	#####################################
40	cp sallcyclev0 specv0
41	cp sallcyclev1 specv1
42	c++compile reducing
43	c++link reducing.so doreduc
44	call doreduc
45
46	newwin 1 1
47	graphicatt "xylimits=1250,1500,-0.01,0.01"
48	del redMeanv0
49	del redMeanv1
50	objaoper reducedMeanDiffOnOff row 0 redMeanv0
51	objaoper reducedMeanDiffOnOff row 1 redMeanv1
52	plot2d redMeanv0 (n/32)*250+1250 val n>0 "blue marker=fcircle,7 nsta notit"
53	plot2d redMeanv1 (n/32)*250+1250 val n>0 "same red marker=fcircle,7 nsta notit"
54	settitle "Red Mean ON-OFF ${source} Ch 0 (blue) Ch 1 (red) ${tmpvec.sum} cycles"
55	setaxelabels "Freq. (MHz)" "I (Jy)"
56
57	newwin 1 1
58	graphicatt "xylimits=1250,1500,0,0.001"
59	del redSigmav0
60	del redSigmav1
61	objaoper reducedSigmaDiffOnOff row 0 redSigmav0
62	objaoper reducedSigmaDiffOnOff row 1 redSigmav1
63	plot2d redSigmav0 (n/32)*250+1250 val n>0 "blue marker=fcircle,7 nsta notit"
64	plot2d redSigmav1 (n/32)*250+1250 val n>0 "same red marker=fcircle,7 nsta notit"
65	settitle "ON-OFF ${source} Ch 0 (blue) Ch 1 (red) ${tmpvec.sum} cycles"
66	setaxelabels "Freq. (MHz)" "I (Jy)"
67
68
69	#####################################
70	#Plots ON-OFF with all calibrations: chan 0 & 1 separaed and cycle per cycle
71	#####################################
72	newwin 1 1
73	graphicatt "xylimits=1250,1500,-0.5,0.5"
74	del sc1v0
75	objaoper specONOFF1 row 0 sc1v0
76	plot2d sc1v0 (n/8192)*250+1250 val n>0 "$cols[0] cpts nsta notit"
77	for ic 2:13
78	k = ${ic}-floor(${ic}/${#cols})*${#cols}
79	del sc${ic}v0
80	objaoper specONOFF${ic} row 0 sc${ic}v0
81	plot2d sc${ic}v0 (n/8192)*250+1250 val n>0 "same $cols[k] cpts nsta notit"
82	end
83	settitle "ON-OFF ${source}/${date} Ch 0 cycle per cycle"
84	setaxelabels "Freq. (MHz)" "I (Jy)"
85
86	newwin 1 1
87	graphicatt "xylimits=1250,1500,-0.5,0.5"
88	del sc1v1
89	objaoper specONOFF1 row 1 sc1v1
90	plot2d sc1v1 (n/8192)*250+1250 val n>0 "$cols[0] cpts nsta notit"
91	for ic 2:13
92	k = ${ic}-floor(${ic}/${#cols})*${#cols}
93	del sc${ic}v1
94	objaoper specONOFF${ic} row 1 sc${ic}v1
95	plot2d sc${ic}v1 (n/8192)*250+1250 val n>0 "same $cols[k] cpts nsta notit"
96	end
97	settitle "ON-OFF ${source}/${date} Ch 1 cycle per cycle"
98	setaxelabels "Freq. (MHz)" "I (Jy)"
99
100
101	#####################################
102	# Plots RAW no calibration at all
103	#####################################
104	openppf ${toppath}/${source}/dataRaw_${date}_${srclower}.ppf
105	newwin 2 2
106	graphicatt "xylimits=1250,1500,0.95,1.05"
107
108	### #
109	### On Ch 0 # On Ch 1
110	######################
111	### Off Ch 0# Off Ch 1
112	### #
113
114	set mode On
115	set first 1
116	for ic 1:13
117	k = ${ic}-floor(${ic}/${#cols})*${#cols}
118	del sm${mode}c${ic}v0
119	objaoper specRaw${mode}${ic} row 0 sm${mode}c${ic}v0
120	if ( ${first} == 1 ) then
121	plot2d sm${mode}c${ic}v0 (n/8192)*250+1250 val n>0 "$cols[k] cpts nsta notit"
122	else
123	plot2d sm${mode}c${ic}v0 (n/8192)*250+1250 val n>0 "same $cols[k] cpts nsta notit"
124	endif
125	set first 0
126	end
127	settitle "${source}/${date} Raw: On Ch 0"
128
129	set first 1
130	for ic 1:13
131	k = ${ic}-floor(${ic}/${#cols})*${#cols}
132	del sm${mode}c${ic}v1
133	objaoper specRaw${mode}${ic} row 1 sm${mode}c${ic}v1
134	if ( ${first} == 1 ) then
135	plot2d sm${mode}c${ic}v1 (n/8192)*250+1250 val n>0 "$cols[k] cpts nsta notit"
136	else
137	plot2d sm${mode}c${ic}v1 (n/8192)*250+1250 val n>0 "same $cols[k] cpts nsta notit"
138	endif
139	set first 0
140	end
141	settitle "${source}/${date} Raw: On Ch 1"
142
143	set mode Off
144	set first 1
145	for ic 1:13
146	k = ${ic}-floor(${ic}/${#cols})*${#cols}
147	del sm${mode}c${ic}v0
148	objaoper specRaw${mode}${ic} row 0 sm${mode}c${ic}v0
149	if ( ${first} == 1 ) then
150	plot2d sm${mode}c${ic}v0 (n/8192)*250+1250 val n>0 "$cols[k] cpts nsta notit"
151	else
152	plot2d sm${mode}c${ic}v0 (n/8192)*250+1250 val n>0 "same $cols[k] cpts nsta notit"
153	endif
154	set first 0
155	end
156	settitle "${source}/${date} Raw: Off Ch 0"
157
158	set first 1
159	for ic 1:13
160	k = ${ic}-floor(${ic}/${#cols})*${#cols}
161	del sm${mode}c${ic}v1
162	objaoper specRaw${mode}${ic} row 1 sm${mode}c${ic}v1
163	if ( ${first} == 1 ) then
164	plot2d sm${mode}c${ic}v1 (n/8192)*250+1250 val n>0 "$cols[k] cpts nsta notit"
165	else
166	plot2d sm${mode}c${ic}v1 (n/8192)*250+1250 val n>0 "same $cols[k] cpts nsta notit"
167	endif
168	set first 0
169	end
170	settitle "${source}/${date} Raw: Off Ch 1"
171
172	#####################################
173	# Plots BAO calibration included
174	#####################################
175	openppf ${toppath}/${source}/dataBAOCalib_${date}_${srclower}.ppf
176	newwin 2 2
177	#graphicatt "xylimits=1250,1500,0.95,1.05"
178	graphicatt ""
179
180
181	### #
182	### On Ch 0 # On Ch 1
183	######################
184	### Off Ch 0# Off Ch 1
185	### #
186
187	set mode On
188	set first 1
189	for ic 1:13
190	k = ${ic}-floor(${ic}/${#cols})*${#cols}
191	del sm${mode}c${ic}v0
192	objaoper specBAOCalib${mode}${ic} row 0 sm${mode}c${ic}v0
193	if ( ${first} == 1 ) then
194	plot2d sm${mode}c${ic}v0 (n/8192)*250+1250 val n>0 "$cols[k] cpts nsta notit"
195	else
196	plot2d sm${mode}c${ic}v0 (n/8192)*250+1250 val n>0 "same $cols[k] cpts nsta notit"
197	endif
198	set first 0
199	end
200	settitle "${source}/${date} BAOCalib: On Ch 0"
201
202	set first 1
203	for ic 1:13
204	k = ${ic}-floor(${ic}/${#cols})*${#cols}
205	del sm${mode}c${ic}v1
206	objaoper specBAOCalib${mode}${ic} row 1 sm${mode}c${ic}v1
207	if ( ${first} == 1 ) then
208	plot2d sm${mode}c${ic}v1 (n/8192)*250+1250 val n>0 "$cols[k] cpts nsta notit"
209	else
210	plot2d sm${mode}c${ic}v1 (n/8192)*250+1250 val n>0 "same $cols[k] cpts nsta notit"
211	endif
212	set first 0
213	end
214	settitle "${source}/${date} BAOCalib: On Ch 1"
215
216	set mode Off
217	set first 1
218	for ic 1:13
219	k = ${ic}-floor(${ic}/${#cols})*${#cols}
220	del sm${mode}c${ic}v0
221	objaoper specBAOCalib${mode}${ic} row 0 sm${mode}c${ic}v0
222	if ( ${first} == 1 ) then
223	plot2d sm${mode}c${ic}v0 (n/8192)*250+1250 val n>0 "$cols[k] cpts nsta notit"
224	else
225	plot2d sm${mode}c${ic}v0 (n/8192)*250+1250 val n>0 "same $cols[k] cpts nsta notit"
226	endif
227	set first 0
228	end
229	settitle "${source}/${date} BAOCalib: Off Ch 0"
230
231	set first 1
232	for ic 1:13
233	k = ${ic}-floor(${ic}/${#cols})*${#cols}
234	del sm${mode}c${ic}v1
235	objaoper specBAOCalib${mode}${ic} row 1 sm${mode}c${ic}v1
236	if ( ${first} == 1 ) then
237	plot2d sm${mode}c${ic}v1 (n/8192)*250+1250 val n>0 "$cols[k] cpts nsta notit"
238	else
239	plot2d sm${mode}c${ic}v1 (n/8192)*250+1250 val n>0 "same $cols[k] cpts nsta notit"
240	endif
241	set first 0
242	end
243	settitle "${source}/${date} BAOCalib: Off Ch 1"
244
245	#####################################
246	# Plots BAO+RT calibration included
247	#####################################
248	openppf ${toppath}/${source}/dataBAORTCalib_${date}_${srclower}.ppf
249	newwin 2 2
250	#graphicatt "xylimits=1250,1500,0.95,1.05"
251	graphicatt ""
252
253
254	### #
255	### On Ch 0 # On Ch 1
256	######################
257	### Off Ch 0# Off Ch 1
258	### #
259
260	set mode On
261	set first 1
262	for ic 1:13
263	k = ${ic}-floor(${ic}/${#cols})*${#cols}
264	del sm${mode}c${ic}v0
265	objaoper specBAORTCalib${mode}${ic} row 0 sm${mode}c${ic}v0
266	if ( ${first} == 1 ) then
267	plot2d sm${mode}c${ic}v0 (n/8192)*250+1250 val n>0 "$cols[k] cpts nsta notit"
268	else
269	plot2d sm${mode}c${ic}v0 (n/8192)*250+1250 val n>0 "same $cols[k] cpts nsta notit"
270	endif
271	set first 0
272	end
273	settitle "${source}/${date} BAORTCalib: On Ch 0"
274
275	set first 1
276	for ic 1:13
277	k = ${ic}-floor(${ic}/${#cols})*${#cols}
278	del sm${mode}c${ic}v1
279	objaoper specBAORTCalib${mode}${ic} row 1 sm${mode}c${ic}v1
280	if ( ${first} == 1 ) then
281	plot2d sm${mode}c${ic}v1 (n/8192)*250+1250 val n>0 "$cols[k] cpts nsta notit"
282	else
283	plot2d sm${mode}c${ic}v1 (n/8192)*250+1250 val n>0 "same $cols[k] cpts nsta notit"
284	endif
285	set first 0
286	end
287	settitle "${source}/${date} BAORTCalib: On Ch 1"
288
289	set mode Off
290	set first 1
291	for ic 1:13
292	k = ${ic}-floor(${ic}/${#cols})*${#cols}
293	del sm${mode}c${ic}v0
294	objaoper specBAORTCalib${mode}${ic} row 0 sm${mode}c${ic}v0
295	if ( ${first} == 1 ) then
296	plot2d sm${mode}c${ic}v0 (n/8192)*250+1250 val n>0 "$cols[k] cpts nsta notit"
297	else
298	plot2d sm${mode}c${ic}v0 (n/8192)*250+1250 val n>0 "same $cols[k] cpts nsta notit"
299	endif
300	set first 0
301	end
302	settitle "${source}/${date} BAORTCalib: Off Ch 0"
303
304	set first 1
305	for ic 1:13
306	k = ${ic}-floor(${ic}/${#cols})*${#cols}
307	del sm${mode}c${ic}v1
308	objaoper specBAORTCalib${mode}${ic} row 1 sm${mode}c${ic}v1
309	if ( ${first} == 1 ) then
310	plot2d sm${mode}c${ic}v1 (n/8192)*250+1250 val n>0 "$cols[k] cpts nsta notit"
311	else
312	plot2d sm${mode}c${ic}v1 (n/8192)*250+1250 val n>0 "same $cols[k] cpts nsta notit"
313	endif
314	set first 0
315	end
316	settitle "${source}/${date} BAORTCalib: Off Ch 1"
317

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

Original Format