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(dab01.NElts()); KeepObj(meandab0d); \ |
---|
51 | TVector<r_4> meandab1d(dab01.NElts()); KeepObj(meandab1d); \ |
---|
52 | meandab0d=dab01d+dab02d; meandab0d/=2.; \ |
---|
53 | meandab1d=dab11d+dab12d; meandab1d/=2.; |
---|
54 | |
---|
55 | ## Save Mean(delta(first,second DAB)) per channel into a ppf file |
---|
56 | #if ( ppf == "1" ) then |
---|
57 | # del daboff |
---|
58 | # c++exec TMatrix<r_4> daboff(2,dab01.NElts()); KeepObj(daboff); \ |
---|
59 | # daboff.Row(0)=meandab0d; daboff.Row(1)=meandab1d; |
---|
60 | # saveppf daboff "${toppath}/${source}/daboff_${date}_${srclower}_calibcycle8.fits" |
---|
61 | #endif |
---|
62 | |
---|
63 | # Plot results |
---|
64 | #-------------- |
---|
65 | graphicatt "xylimits=1250,1500,0,3" |
---|
66 | |
---|
67 | ## Baseline |
---|
68 | newwin 1 1 |
---|
69 | plot2d b0 (n/8192)*250+1250 val n>0 "blue cpts nstat notit" |
---|
70 | plot2d b1 (n/8192)*250+1250 val n>0 "same red cpts nstat notit" |
---|
71 | settitle "Baseline raw $source $date $mode Ch0 blue Ch1 red" |
---|
72 | |
---|
73 | newwin 1 1 |
---|
74 | plot2d b0f (n/8192)*250+1250 val n>0 "blue cpts nstat notit" |
---|
75 | plot2d b1f (n/8192)*250+1250 val n>0 "same red cpts nstat notit" |
---|
76 | settitle "Baseline filt $source $date $mode Ch0 blue Ch1 red" |
---|
77 | |
---|
78 | ## First DAB |
---|
79 | #newwin 1 1 |
---|
80 | #plot2d dab01 (n/8192)*250+1250 val n>0 "blue cpts nstat notit" |
---|
81 | #plot2d dab11 (n/8192)*250+1250 val n>0 "same red cpts nstat notit" |
---|
82 | #settitle "First DAB raw $source $date $mode Ch0 blue Ch1 red" |
---|
83 | |
---|
84 | newwin 1 1 |
---|
85 | plot2d dab01f (n/8192)*250+1250 val n>0 "blue cpts nstat notit" |
---|
86 | plot2d dab11f (n/8192)*250+1250 val n>0 "same red cpts nstat notit" |
---|
87 | settitle "First DAB filt $source $date $mode Ch0 blue Ch1 red" |
---|
88 | |
---|
89 | ## Second DAB |
---|
90 | #newwin 1 1 |
---|
91 | #plot2d dab02 (n/8192)*250+1250 val n>0 "blue cpts nstat notit" |
---|
92 | #plot2d dab12 (n/8192)*250+1250 val n>0 "same red cpts nstat notit" |
---|
93 | #settitle "Second DAB raw $source $date $mode Ch0 blue Ch1 red" |
---|
94 | |
---|
95 | newwin 1 1 |
---|
96 | plot2d dab02f (n/8192)*250+1250 val n>0 "blue cpts nstat notit" |
---|
97 | plot2d dab12f (n/8192)*250+1250 val n>0 "same red cpts nstat notit" |
---|
98 | settitle "Second DAB filt $source $date $mode Ch0 blue Ch1 red" |
---|
99 | |
---|
100 | ## Delta=DAB-baseline + Calibrated DAB in Jansky |
---|
101 | graphicatt "xylimits=1100,1600,0,30" |
---|
102 | newwin 1 1 |
---|
103 | n/pl nt_dabjy.dabjy%dabfreq 1 ! "marker=star,4 turquoise notit nstat" |
---|
104 | plot2d dab01d (n/8192)*250+1250 val*18 n>0 "same blue cpts nstat notit" |
---|
105 | plot2d dab11d (n/8192)*250+1250 val*22 n>0 "same red cpts nstat notit" |
---|
106 | settitle "First DAB delta $source $date $mode" |
---|
107 | setaxelabels "Freq (MHz)" "Intensity" "bigboldfont" |
---|
108 | addline 1420 4 1430 4 "turquoise line=solid,3" |
---|
109 | addline 1420 3 1430 3 "blue line=solid,3" |
---|
110 | addline 1420 2 1430 2 "red line=solid,3" |
---|
111 | addtext 1440 4 "DAB Voie E (Jy)" "bigfont" |
---|
112 | addtext 1440 3 "Ch0*18 (A.U.)" "bigfont" |
---|
113 | addtext 1440 2 "Ch1*22 (A.U.)" "bigfont" |
---|
114 | |
---|
115 | newwin 1 1 |
---|
116 | n/pl nt_dabjy.dabjy%dabfreq 1 ! "marker=star,4 turquoise notit nstat" |
---|
117 | plot2d dab02d (n/8192)*250+1250 val*18 n>0 "same blue cpts nstat notit" |
---|
118 | plot2d dab12d (n/8192)*250+1250 val*22 n>0 "same red cpts nstat notit" |
---|
119 | settitle "Second DAB delta $source $date $mode" |
---|
120 | setaxelabels "Freq (MHz)" "Intensity" "bigboldfont" |
---|
121 | addline 1420 4 1430 4 "turquoise line=solid,3" |
---|
122 | addline 1420 3 1430 3 "blue line=solid,3" |
---|
123 | addline 1420 2 1430 2 "red line=solid,3" |
---|
124 | addtext 1440 4 "DAB Voie E (Jy)" "bigfont" |
---|
125 | addtext 1440 3 "Ch0*18 (A.U.)" "bigfont" |
---|
126 | addtext 1440 2 "Ch1*22 (A.U.)" "bigfont" |
---|
127 | |
---|
128 | ## Delta=Mean(DAB-baseline) + Calibrated DAB in Jansky |
---|
129 | graphicatt "xylimits=1100,1600,0,30" |
---|
130 | newwin 1 1 |
---|
131 | n/pl nt_dabjy.dabjy%dabfreq 1 ! "marker=star,4 turquoise notit nstat" |
---|
132 | plot2d meandab0d (n/8192)*250+1250 val*18 n>0 "same blue cpts nstat notit" |
---|
133 | plot2d meandab1d (n/8192)*250+1250 val*23 n>0 "same red cpts nstat notit" |
---|
134 | settitle "Mean DAB delta $source $date $mode" |
---|
135 | setaxelabels "Freq (MHz)" "Intensity" "bigboldfont" |
---|
136 | addline 1420 4 1430 4 "turquoise line=solid,3" |
---|
137 | addline 1420 3 1430 3 "blue line=solid,3" |
---|
138 | addline 1420 2 1430 2 "red line=solid,3" |
---|
139 | addtext 1440 4 "DAB Voie E (Jy)" "bigfont" |
---|
140 | addtext 1440 3 "Ch0*18 (A.U.)" "bigfont" |
---|
141 | addtext 1440 2 "Ch1*23 (A.U.)" "bigfont" |
---|