1 | set user "AST" |
---|
2 | set source "Abell1205" |
---|
3 | set srclower "abell1205" |
---|
4 | set date "20110415" |
---|
5 | set cycle "8" |
---|
6 | set mode "Off" |
---|
7 | set toppath "/sps/baoradio/AmasNancay/${user}" |
---|
8 | set path "${toppath}/${source}/dab" |
---|
9 | set ppf "0" |
---|
10 | |
---|
11 | # Read DAB calibration in Jansky (Voie E / PKS1127-14) |
---|
12 | #----------------------------------------------------- |
---|
13 | n/read nt_dabjy ${toppath}/DAB_Jy.dat " " dabfreq,0 dabjy,1 |
---|
14 | newwin 1 1 |
---|
15 | n/pl nt_dabjy.dabjy%dabfreq 1 ! "marker=star,4 magenta notit" |
---|
16 | settitle "DAB Voie E / PKS1127-14. Juin 2011" |
---|
17 | |
---|
18 | |
---|
19 | # Read DAB calibration in A.U. (from BAORadio data) |
---|
20 | #-------------------------------------------------- |
---|
21 | openppf ${path}/dabchrom_${date}_${srclower}_${mode}.ppf |
---|
22 | ls |
---|
23 | |
---|
24 | ## Baseline |
---|
25 | del b0 b1 b0f b1f |
---|
26 | objaoper /home/base row 0 b0 |
---|
27 | objaoper /home/base row 1 b1 |
---|
28 | objaoper /home/basefilt row 0 b0f |
---|
29 | objaoper /home/basefilt row 1 b1f |
---|
30 | |
---|
31 | ## First DAB |
---|
32 | del dab01 dab11 dab01f dab11f dab01n dab11n |
---|
33 | objaoper /home/firstdab row 0 dab01 |
---|
34 | objaoper /home/firstdab row 1 dab11 |
---|
35 | objaoper /home/firstdabfilt row 0 dab01f |
---|
36 | objaoper /home/firstdabfilt row 1 dab11f |
---|
37 | objaoper /home/firstdelta row 0 dab01d |
---|
38 | objaoper /home/firstdelta row 1 dab11d |
---|
39 | |
---|
40 | ## Second DAB |
---|
41 | del dab02 dab12 dab02f dab12f dab02n dab12n |
---|
42 | objaoper /home/seconddab row 0 dab02 |
---|
43 | objaoper /home/seconddab row 1 dab12 |
---|
44 | objaoper /home/seconddabfilt row 0 dab02f |
---|
45 | objaoper /home/seconddabfilt row 1 dab12f |
---|
46 | objaoper /home/seconddelta row 0 dab02d |
---|
47 | objaoper /home/seconddelta row 1 dab12d |
---|
48 | |
---|
49 | ## Mean(delta(first,second DAB)) per channel |
---|
50 | c++exec TVector<r_4> meandab0d(dab01d,false); KeepObj(meandab0d); \ |
---|
51 | TVector<r_4> meandab1d(dab11d,false); KeepObj(meandab1d); \ |
---|
52 | meandab0d+=dab02d; meandab0d/=2.; \ |
---|
53 | meandab1d+=dab12d; meandab1d/=2.; |
---|
54 | |
---|
55 | ## Rebin BAO and calibrated data in 0.5Hz/bin to perform division |
---|
56 | exptovec vf nt_dabjy dabfreq |
---|
57 | exptovec vjy nt_dabjy dabjy |
---|
58 | c++exec TVector<r_4> dumm0(vf), dumm1(vjy);\ |
---|
59 | TVector<r_4> vcalfreq=dumm0.SubVector(Range(320,854)).TransposeSelf();\ |
---|
60 | KeepObj(vcalfreq);\ |
---|
61 | TVector<r_4> vcaljy=dumm1.SubVector(Range(320,854)).TransposeSelf();\ |
---|
62 | KeepObj(vcaljy); |
---|
63 | |
---|
64 | del invec0 invec1 invec2 invec3 val0 val1 |
---|
65 | cp meandab0d invec0 |
---|
66 | cp meandab1d invec1 |
---|
67 | cp vcaljy invec2 |
---|
68 | cp vcalfreq invec3 |
---|
69 | ls |
---|
70 | |
---|
71 | c++compile interpol |
---|
72 | c++link interpol.so dointerpol |
---|
73 | call dointerpol |
---|
74 | |
---|
75 | c++exec TVector<r_8> calcoeff0f1400a1420=calcoeff0.SubVector(Range(300,340)); \ |
---|
76 | KeepObj(calcoeff0f1400a1420);\ |
---|
77 | TVector<r_8> calcoeff1f1400a1420=calcoeff1.SubVector(Range(300,340)); \ |
---|
78 | KeepObj(calcoeff1f1400a1420); |
---|
79 | projh1d h0 calcoeff0f1400a1420 val |
---|
80 | projh1d h1 calcoeff1f1400a1420 val |
---|
81 | |
---|
82 | ## Save Mean(delta(first,second DAB)) per channel into a ppf file |
---|
83 | #if ( ppf == "1" ) then |
---|
84 | # del daboff |
---|
85 | # c++exec TMatrix<r_4> daboff(2,dab01.NElts()); KeepObj(daboff); \ |
---|
86 | # daboff.Row(0)=meandab0d; daboff.Row(1)=meandab1d; |
---|
87 | # saveppf daboff "${toppath}/${source}/daboff_${date}_${srclower}_calibcycle8.fits" |
---|
88 | #endif |
---|
89 | |
---|
90 | # Plot results |
---|
91 | #-------------- |
---|
92 | graphicatt "xylimits=1250,1500,0,3" |
---|
93 | |
---|
94 | ## Baseline |
---|
95 | newwin 1 1 |
---|
96 | plot2d b0 (n/8192)*250+1250 val n>0 "blue cpts nstat notit" |
---|
97 | plot2d b1 (n/8192)*250+1250 val n>0 "same red cpts nstat notit" |
---|
98 | settitle "Baseline raw $source $date $mode Ch0 blue Ch1 red" |
---|
99 | |
---|
100 | newwin 1 1 |
---|
101 | plot2d b0f (n/8192)*250+1250 val n>0 "blue cpts nstat notit" |
---|
102 | plot2d b1f (n/8192)*250+1250 val n>0 "same red cpts nstat notit" |
---|
103 | settitle "Baseline filt $source $date $mode Ch0 blue Ch1 red" |
---|
104 | |
---|
105 | ## First DAB |
---|
106 | #newwin 1 1 |
---|
107 | #plot2d dab01 (n/8192)*250+1250 val n>0 "blue cpts nstat notit" |
---|
108 | #plot2d dab11 (n/8192)*250+1250 val n>0 "same red cpts nstat notit" |
---|
109 | #settitle "First DAB raw $source $date $mode Ch0 blue Ch1 red" |
---|
110 | |
---|
111 | newwin 1 1 |
---|
112 | plot2d dab01f (n/8192)*250+1250 val n>0 "blue cpts nstat notit" |
---|
113 | plot2d dab11f (n/8192)*250+1250 val n>0 "same red cpts nstat notit" |
---|
114 | settitle "First DAB filt $source $date $mode Ch0 blue Ch1 red" |
---|
115 | |
---|
116 | ## Second DAB |
---|
117 | #newwin 1 1 |
---|
118 | #plot2d dab02 (n/8192)*250+1250 val n>0 "blue cpts nstat notit" |
---|
119 | #plot2d dab12 (n/8192)*250+1250 val n>0 "same red cpts nstat notit" |
---|
120 | #settitle "Second DAB raw $source $date $mode Ch0 blue Ch1 red" |
---|
121 | |
---|
122 | newwin 1 1 |
---|
123 | plot2d dab02f (n/8192)*250+1250 val n>0 "blue cpts nstat notit" |
---|
124 | plot2d dab12f (n/8192)*250+1250 val n>0 "same red cpts nstat notit" |
---|
125 | settitle "Second DAB filt $source $date $mode Ch0 blue Ch1 red" |
---|
126 | |
---|
127 | ## Delta=DAB-baseline + Calibrated DAB in Jansky |
---|
128 | graphicatt "xylimits=1100,1600,0,30" |
---|
129 | newwin 1 1 |
---|
130 | n/pl nt_dabjy.dabjy%dabfreq 1 ! "marker=star,4 turquoise notit nstat" |
---|
131 | plot2d dab01d (n/8192)*250+1250 val*18 n>0 "same blue cpts nstat notit" |
---|
132 | plot2d dab11d (n/8192)*250+1250 val*22 n>0 "same red cpts nstat notit" |
---|
133 | settitle "First DAB delta $source $date $mode" |
---|
134 | setaxelabels "Freq (MHz)" "Intensity" "bigboldfont" |
---|
135 | addline 1420 4 1430 4 "turquoise line=solid,3" |
---|
136 | addline 1420 3 1430 3 "blue line=solid,3" |
---|
137 | addline 1420 2 1430 2 "red line=solid,3" |
---|
138 | addtext 1440 4 "DAB Voie E (Jy)" "bigfont" |
---|
139 | addtext 1440 3 "Ch0*18 (A.U.)" "bigfont" |
---|
140 | addtext 1440 2 "Ch1*22 (A.U.)" "bigfont" |
---|
141 | |
---|
142 | newwin 1 1 |
---|
143 | n/pl nt_dabjy.dabjy%dabfreq 1 ! "marker=star,4 turquoise notit nstat" |
---|
144 | plot2d dab02d (n/8192)*250+1250 val*18 n>0 "same blue cpts nstat notit" |
---|
145 | plot2d dab12d (n/8192)*250+1250 val*22 n>0 "same red cpts nstat notit" |
---|
146 | settitle "Second DAB delta $source $date $mode" |
---|
147 | setaxelabels "Freq (MHz)" "Intensity" "bigboldfont" |
---|
148 | addline 1420 4 1430 4 "turquoise line=solid,3" |
---|
149 | addline 1420 3 1430 3 "blue line=solid,3" |
---|
150 | addline 1420 2 1430 2 "red line=solid,3" |
---|
151 | addtext 1440 4 "DAB Voie E (Jy)" "bigfont" |
---|
152 | addtext 1440 3 "Ch0*18 (A.U.)" "bigfont" |
---|
153 | addtext 1440 2 "Ch1*22 (A.U.)" "bigfont" |
---|
154 | |
---|
155 | ## Delta=Mean(DAB-baseline) + Calibrated DAB in Jansky |
---|
156 | graphicatt "xylimits=1100,1600,0,30" |
---|
157 | newwin 1 1 |
---|
158 | n/pl nt_dabjy.dabjy%dabfreq 1 ! "marker=star,4 turquoise notit nstat" |
---|
159 | plot2d meandab0d (n/8192)*250+1250 val*18 n>0 "same blue cpts nstat notit" |
---|
160 | plot2d meandab1d (n/8192)*250+1250 val*23 n>0 "same red cpts nstat notit" |
---|
161 | settitle "Mean DAB delta $source $date $mode" |
---|
162 | setaxelabels "Freq (MHz)" "Intensity" "bigboldfont" |
---|
163 | addline 1420 4 1430 4 "turquoise line=solid,3" |
---|
164 | addline 1420 3 1430 3 "blue line=solid,3" |
---|
165 | addline 1420 2 1430 2 "red line=solid,3" |
---|
166 | addtext 1440 4 "DAB Voie E (Jy)" "bigfont" |
---|
167 | addtext 1440 3 "Ch0*18 (A.U.)" "bigfont" |
---|
168 | addtext 1440 2 "Ch1*23 (A.U.)" "bigfont" |
---|
169 | |
---|
170 | graphicatt "" |
---|
171 | newwin 1 1 |
---|
172 | plot2d meandab0d (n/8192)*250+1250 val n>0 "blue cpts nstat notit" |
---|
173 | plot2d meandab1d (n/8192)*250+1250 val n>0 "same red cpts nstat notit" |
---|
174 | settitle "Raw MeanDAB Ch0(blue) Ch1(red)" |
---|
175 | |
---|
176 | ## Rebinned DAB spectra |
---|
177 | graphicatt "" |
---|
178 | newwin 1 1 |
---|
179 | plot2d val0 (n/500)*250+1250 val 1 "cpts blue notit nstat" |
---|
180 | plot2d val1 (n/500)*250+1250 val 1 "cpts same red notit nstat" |
---|
181 | settitle "BAO DAB spectra Ch0 (blue) Ch1(red)" |
---|
182 | setaxelabels "Freq. (MHz)" "I (a.u.)" "bigboldfont" |
---|
183 | |
---|
184 | newwin 1 1 |
---|
185 | plot2d val2 (n/500)*250+1250 val 1 "cpts blue notit nstat" |
---|
186 | settitle "Calibrated DAB spectra" |
---|
187 | setaxelabels "Freq. (MHz)" "I (Jy)" "bigboldfont" |
---|
188 | |
---|
189 | ## Calibration coefficients Jy/a.u. in the whole bandwidth |
---|
190 | newwin 1 1 |
---|
191 | plot2d calcoeff1 (n/500)*250+1250 val 1 "cpts red notit nstat" |
---|
192 | plot2d calcoeff0 (n/500)*250+1250 val 1 "cpts same blue notit nstat" |
---|
193 | settitle "Calibration coefficients spectrum Ch0(blue) Ch1(red)" |
---|
194 | setaxelabels "Freq. (MHz)" "Coeff Jy/a.u." "bigboldfont" |
---|
195 | |
---|
196 | ## Calibration coefficients Jy/a.u. in [1400,1420]MHz |
---|
197 | newwin 1 1 |
---|
198 | disp h0 "notit" |
---|
199 | settitle "CalCoeff (Jy/a.u.) Ch0" |
---|
200 | newwin 1 1 |
---|
201 | disp h1 "notit" |
---|
202 | settitle "CalCoeff (Jy/a.u.) Ch1" |
---|