Changeset 639 for BAORadio/AmasNancay/trunk/imgdrift.pic

Timestamp:

Jan 17, 2012, 2:15:02 PM (12 years ago)

Author:

campagne

Message:

improve interactive drift analysis (jec)

File:

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

BAORadio/AmasNancay/trunk/imgdrift.pic

@@ -10 +10 @@
 set path "/sps/baoradio/AmasNancay/JEC/${srcMaj}"
 
+#first script executed
 clearscript image
+
+#mergeimg should be executed before timeevol, fwhnInt
+clearscript mergeimg
+clearscript timeevol
+clearscript fwhmInt
+
+#script independant from mergeimg
+clearscript timeevol2
 #######################
 defscript image
@@ -37 +46 @@
 #######################
 
-for ic ${fcycle}:${maxcycle}
-echo "process $ic"
-image $ic
-end
-
+
+#######################
+defscript mergeimg
+#
 del ich0Tot
 del ich1Tot
@@ -66 +74 @@
 newwin 
 disp ich1Tot 'lut=lin,-0.1,3.5 colbr128 showcmap=top'
-
-#intensite en ft du temps
+endscript
+#######################
+defscript timeevol
+#
+#intensite en ft du temps avec image cumulee
 del rg0Tot
 del rg1Tot
@@ -96 +107 @@
 plot2d rg0Tot (n*170)/273 val 1 'blue cpts notit nsta'
 plot2d rg1Tot (n*170)/273 val 1 'same red cpts notit nsta'
-settitle "${srcMaj} Drift Scan [1405,1415]MHz 2cycles Ch 0 (blue) Ch 1 (red)"
+settitle "${srcMaj} Drift Scan [1405,1415]MHz  all cycles Ch 0 (blue) Ch 1 (red)"
 setaxelabels 't(sec)' 'I (a.u)'
-
+endscript
+
+#######################
+defscript fwhmInt
+#
 #spectre de frequence canal en temps [120,146] FWHM.
 #
@@ -124 +139 @@
 plot2d s0 (n/8192)*250+1250 val n>0 'blue cpts notit nsta'
 plot2d s1 (n/8192)*250+1250 val 1 'same red cpts notit nsta'
-settitle "${srcMaj} Drift Scan 2cycles Ch 0 (blue) Ch 1 (red)"
+settitle "${srcMaj} Drift Scan all cycles Ch 0 (blue) Ch 1 (red)"
 setaxelabels 'Freq. (MHz)' 'I (a.u)'
-
+endscript
+#######################
+#######################
+defscript timeevol2
+
+#append all images 
+del ch0evol
+del ch1evol
+c++exec sa_size_t nrows=(${nframes}+1)*(${lcycle}-${fcycle}+1); \
+        sa_size_t ncols=ich${fcycle}0.NCols(); \
+        TMatrix<r_4>ch0evol(nrows,ncols); TMatrix<r_4>ch1evol(nrows,ncols); \
+        KeepObj(ch0evol); KeepObj(ch1evol);
+
+for ic ${fcycle}:${maxcycle}
+  c++exec sa_size_t curRow=(${ic}-${fcycle})*(${nframes}+1); sa_size_t nextm1CurRow=curRow+(${nframes}+1)-1; \
+          cout << "ic, curRow, next : "<<${ic}<<" "<<curRow<<" "<<nextm1CurRow<<endl; \
+          ch0evol(Range(curRow,nextm1CurRow),Range::all())=ich${ic}0(Range(0,${nframes}),Range::all()); \
+          ch1evol(Range(curRow,nextm1CurRow),Range::all())=ich${ic}1(Range(0,${nframes}),Range::all());
+end
+
+echo "Coucou......"
+
+#intensite en ft du temps avec image cumulee
+del i0evol
+del i1evol
+#1405,1415 MHz => range [5080,5405]
+c++exec TMatrix<r_4>tmp0(ch0evol(Range::all(),Range(5080,5405)),false); \
+        TVector<r_4>i0evol(ch0evol.NRows()); \
+        for(sa_size_t i=0;i<tmp0.NRows();++i){ \
+          double mean,sigma; \
+          MeanSigma(tmp0(Range(i),Range::all()),mean,sigma); \
+          i0evol(i) = mean; \
+        } \
+        KeepObj(i0evol); \
+        TMatrix<r_4>tmp1(ch1evol(Range::all(),Range(5080,5405)),false); \
+        TVector<r_4>i1evol(ch1evol.NRows()); \
+        for(sa_size_t i=0;i<tmp1.NRows();++i){ \
+          double mean,sigma; \
+          MeanSigma(tmp1(Range(i),Range::all()),mean,sigma); \
+          i1evol(i) = mean; \
+        } \
+        KeepObj(i1evol);
+
+
+
+newwin
+setaxesatt 'grid'
+graphicatt '' 
+plot2d i0evol n  val 1 'blue cpts notit nsta'
+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)'
+        
+endscript
+###################
+#for ic ${fcycle}:${maxcycle}
+#echo "process $ic"
+#image $ic
+#end
+