[611] | 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" |
---|