Changeset 640 for BAORadio

Timestamp:

Jan 18, 2012, 3:14:26 PM (13 years ago)

Author:

campagne

Message:

add timeevol2 for Tsys survey

File:

• BAORadio/AmasNancay/trunk/imgdrift.pic (modified) (4 diffs)

BAORadio/AmasNancay/trunk/imgdrift.pic

@@ -12 +12 @@
 #first script executed
 clearscript image
+clearscript doImage
 
 #mergeimg should be executed before timeevol, fwhnInt
@@ -23 +24 @@
 defscript image
 #
+#compute (ON-OFF)/OFFn with OFF computed from the first 30s of each image
 set cycle  $1
 openppf ${path}/img_${date}_${source}_cycle${cycle}.ppf 
 del ich0
+del zCh0
 objaoper img${cycle} slicexz 0 ich0
 c++exec TMatrix<r_4>zCh0(ich0(Range(0),Range::all()),false); \
 for (sa_size_t i=1;i<50;i++){zCh0+=ich0(Range(i),Range::all());} \
-zCh0/=50.; \
-for (sa_size_t i=0;i<ich0.NRows();i++){ich0(Range(i),Range::all())-=zCh0;}
+zCh0/=50.; KeepObj(zCh0);\
+for (sa_size_t i=0;i<ich0.NRows();i++){ich0(Range(i),Range::all())-=zCh0;ich0(Range(i),Range::all()).Div(zCh0);}
 #
 rename ich0 ich${cycle}0
+rename zCh0 tsys${cycle}0
 
 del ich1
+del zCh1
 objaoper img${cycle} slicexz 1 ich1
 c++exec TMatrix<r_4>zCh1(ich1(Range(0),Range::all()),false); \
 for (sa_size_t i=1;i<50;i++){zCh1+=ich1(Range(i),Range::all());} \
-zCh1/=50.; \
-for (sa_size_t i=0;i<ich1.NRows();i++){ich1(Range(i),Range::all())-=zCh1;}
+zCh1/=50.; KeepObj(zCh1);\
+for (sa_size_t i=0;i<ich1.NRows();i++){ich1(Range(i),Range::all())-=zCh1;ich1(Range(i),Range::all()).Div(zCh1);}
 #
 rename ich1 ich${cycle}1
-
-endscript
-#######################
+rename zCh1 tsys${cycle}1
+
+endscript
+#######################
+###################
+defscript doImage
+
+for ic ${fcycle}:${maxcycle}
+echo "process ${ic}"
+image $ic
+end
+endscript
 
 
@@ -161 +175 @@
 end
 
-echo "Coucou......"
-
 #intensite en ft du temps avec image cumulee
 del i0evol
@@ -192 +204 @@
 plot2d i1evol n  val 1 'same red cpts notit nsta'
 settitle "${srcMaj} Drift Scan [1405,1415]MHz all cycles Ch 0 (blue) Ch 1 (red)"
-setaxelabels 't(a.u)' 'I (a.u)'
+setaxelabels 't(a.u)' '(I-OFF)/OFF (a.u)'
+
+
+#Tsys en ft du temps
+c++exec sa_size_t nrows=(${lcycle}-${fcycle}+1); \
+        sa_size_t ncols=ich${fcycle}0.NCols(); \
+        TMatrix<r_4>Tsys0evol(nrows,ncols); TMatrix<r_4>Tsys1evol(nrows,ncols); \
+        KeepObj(Tsys0evol); KeepObj(Tsys1evol);
+
+for ic ${fcycle}:${maxcycle}
+  c++exec sa_size_t curRow=(${ic}-${fcycle}); \
+          cout << "ic, curRow, next : "<<${ic}<<" "<<curRow<<endl; \
+          Tsys0evol(Range(curRow),Range::all())=tsys${ic}0(Range(0),Range::all()); \
+          Tsys1evol(Range(curRow),Range::all())=tsys${ic}1(Range(0),Range::all());
+end
+
+del T0evol
+del T1evol
+#1405,1415 MHz => range [5080,5405]
+c++exec TMatrix<r_4>tmp0(Tsys0evol(Range::all(),Range(5080,5405)),false); \
+        TVector<r_4>T0evol(Tsys0evol.NRows()); \
+        for(sa_size_t i=0;i<tmp0.NRows();++i){ \
+          double mean,sigma; \
+          MeanSigma(tmp0(Range(i),Range::all()),mean,sigma); \
+          T0evol(i) = mean; \
+        } \
+        KeepObj(T0evol); \
+        TMatrix<r_4>tmp1(Tsys1evol(Range::all(),Range(5080,5405)),false); \
+        TVector<r_4>T1evol(Tsys1evol.NRows()); \
+        for(sa_size_t i=0;i<tmp1.NRows();++i){ \
+          double mean,sigma; \
+          MeanSigma(tmp1(Range(i),Range::all()),mean,sigma); \
+          T1evol(i) = mean; \
+        } \
+        KeepObj(T1evol);
+
+newwin
+setaxesatt 'grid'
+graphicatt '' 
+plot2d T0evol n+${fcycle}  val 1 'blue cpts notit nsta'
+plot2d T1evol n+${fcycle}  val 1 'same red cpts notit nsta'
+settitle "${srcMaj} Drift Scan [1405,1415]MHz all cycles Ch 0 (blue) Ch 1 (red)"
+setaxelabels 'cycle num.' 'OFF/Gain (a.u)'
+
         
 endscript
-###################
-#for ic ${fcycle}:${maxcycle}
-#echo "process $ic"
-#image $ic
-#end
-
+