Changeset 630 for BAORadio

Timestamp:

Jan 13, 2012, 11:46:16 AM (12 years ago)

Author:

torrento

Message:

Added calibration coefficient spectrum

File:

• BAORadio/AmasNancay/trunk/dabchrom_plot.pic (modified) (2 diffs)

BAORadio/AmasNancay/trunk/dabchrom_plot.pic

@@ -48 +48 @@
 
 ## Mean(delta(first,second DAB)) per channel
-c++exec TVector<r_4> meandab0d(dab01.NElts()); KeepObj(meandab0d); \
-        TVector<r_4> meandab1d(dab01.NElts()); KeepObj(meandab1d); \
-        meandab0d=dab01d+dab02d; meandab0d/=2.; \
-        meandab1d=dab11d+dab12d; meandab1d/=2.;
+c++exec TVector<r_4> meandab0d(dab01d,false); KeepObj(meandab0d); \
+        TVector<r_4> meandab1d(dab11d,false); KeepObj(meandab1d); \
+        meandab0d+=dab02d; meandab0d/=2.; \
+        meandab1d+=dab12d; meandab1d/=2.;
+
+## Rebin BAO and calibrated data in 0.5Hz/bin to perform division
+exptovec vf nt_dabjy dabfreq
+exptovec vjy nt_dabjy dabjy
+c++exec TVector<r_4> dumm0(vf), dumm1(vjy);\
+        TVector<r_4> vcalfreq=dumm0.SubVector(Range(320,854)).TransposeSelf();\
+        KeepObj(vcalfreq);\
+        TVector<r_4> vcaljy=dumm1.SubVector(Range(320,854)).TransposeSelf();\
+        KeepObj(vcaljy);
+
+del invec0 invec1 invec2 invec3 val0 val1
+cp meandab0d invec0
+cp meandab1d invec1
+cp vcaljy invec2
+cp vcalfreq invec3
+ls
+
+c++compile interpol
+c++link interpol.so dointerpol
+call dointerpol
+
+c++exec TVector<r_8> calcoeff0f1400a1420=calcoeff0.SubVector(Range(300,340)); \
+        KeepObj(calcoeff0f1400a1420);\
+        TVector<r_8> calcoeff1f1400a1420=calcoeff1.SubVector(Range(300,340)); \
+        KeepObj(calcoeff1f1400a1420);
+projh1d h0 calcoeff0f1400a1420 val
+projh1d h1 calcoeff1f1400a1420 val
 
 ## Save Mean(delta(first,second DAB)) per channel into a ppf file
@@ -140 +167 @@
 addtext 1440 3 "Ch0*18 (A.U.)" "bigfont"
 addtext 1440 2 "Ch1*23 (A.U.)" "bigfont"
+
+graphicatt ""
+newwin 1 1
+plot2d meandab0d (n/8192)*250+1250 val n>0 "blue cpts nstat notit"
+plot2d meandab1d (n/8192)*250+1250 val n>0 "same red cpts nstat notit"
+settitle "Raw MeanDAB Ch0(blue) Ch1(red)"
+
+## Rebinned DAB spectra
+graphicatt ""
+newwin 1 1
+plot2d val0 (n/500)*250+1250 val 1 "cpts blue notit nstat"
+plot2d val1 (n/500)*250+1250 val 1 "cpts same red notit nstat"
+settitle "BAO DAB spectra Ch0 (blue) Ch1(red)"
+setaxelabels "Freq. (MHz)" "I (a.u.)" "bigboldfont"
+
+newwin 1 1
+plot2d val2 (n/500)*250+1250 val 1 "cpts blue notit nstat"
+settitle "Calibrated DAB spectra"
+setaxelabels "Freq. (MHz)" "I (Jy)" "bigboldfont"
+
+## Calibration coefficients Jy/a.u. in the whole bandwidth
+newwin 1 1
+plot2d calcoeff1 (n/500)*250+1250 val 1 "cpts red notit nstat"
+plot2d calcoeff0 (n/500)*250+1250 val 1 "cpts same blue notit nstat"
+settitle "Calibration coefficients spectrum Ch0(blue) Ch1(red)"
+setaxelabels "Freq. (MHz)" "Coeff Jy/a.u."  "bigboldfont"
+
+## Calibration coefficients Jy/a.u. in [1400,1420]MHz
+newwin 1 1
+disp h0 "notit"
+settitle "CalCoeff (Jy/a.u.) Ch0"
+newwin 1 1
+disp h1 "notit"
+settitle "CalCoeff (Jy/a.u.) Ch1"