Changeset 630 for BAORadio/AmasNancay/trunk/dabchrom_plot.pic
- Timestamp:
- Jan 13, 2012, 11:46:16 AM (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
BAORadio/AmasNancay/trunk/dabchrom_plot.pic
r611 r630 48 48 49 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.; 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 54 81 55 82 ## Save Mean(delta(first,second DAB)) per channel into a ppf file … … 140 167 addtext 1440 3 "Ch0*18 (A.U.)" "bigfont" 141 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"
Note: See TracChangeset
for help on using the changeset viewer.