Changeset 639 for BAORadio/AmasNancay/trunk
- Timestamp:
- Jan 17, 2012, 2:15:02 PM (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
BAORadio/AmasNancay/trunk/imgdrift.pic
r638 r639 10 10 set path "/sps/baoradio/AmasNancay/JEC/${srcMaj}" 11 11 12 #first script executed 12 13 clearscript image 14 15 #mergeimg should be executed before timeevol, fwhnInt 16 clearscript mergeimg 17 clearscript timeevol 18 clearscript fwhmInt 19 20 #script independant from mergeimg 21 clearscript timeevol2 13 22 ####################### 14 23 defscript image … … 37 46 ####################### 38 47 39 for ic ${fcycle}:${maxcycle} 40 echo "process $ic" 41 image $ic 42 end 43 48 49 ####################### 50 defscript mergeimg 51 # 44 52 del ich0Tot 45 53 del ich1Tot … … 66 74 newwin 67 75 disp ich1Tot 'lut=lin,-0.1,3.5 colbr128 showcmap=top' 68 69 #intensite en ft du temps 76 endscript 77 ####################### 78 defscript timeevol 79 # 80 #intensite en ft du temps avec image cumulee 70 81 del rg0Tot 71 82 del rg1Tot … … 96 107 plot2d rg0Tot (n*170)/273 val 1 'blue cpts notit nsta' 97 108 plot2d rg1Tot (n*170)/273 val 1 'same red cpts notit nsta' 98 settitle "${srcMaj} Drift Scan [1405,1415]MHz 2cycles Ch 0 (blue) Ch 1 (red)"109 settitle "${srcMaj} Drift Scan [1405,1415]MHz all cycles Ch 0 (blue) Ch 1 (red)" 99 110 setaxelabels 't(sec)' 'I (a.u)' 100 111 endscript 112 113 ####################### 114 defscript fwhmInt 115 # 101 116 #spectre de frequence canal en temps [120,146] FWHM. 102 117 # … … 124 139 plot2d s0 (n/8192)*250+1250 val n>0 'blue cpts notit nsta' 125 140 plot2d s1 (n/8192)*250+1250 val 1 'same red cpts notit nsta' 126 settitle "${srcMaj} Drift Scan 2cycles Ch 0 (blue) Ch 1 (red)"141 settitle "${srcMaj} Drift Scan all cycles Ch 0 (blue) Ch 1 (red)" 127 142 setaxelabels 'Freq. (MHz)' 'I (a.u)' 128 143 endscript 144 ####################### 145 ####################### 146 defscript timeevol2 147 148 #append all images 149 del ch0evol 150 del ch1evol 151 c++exec sa_size_t nrows=(${nframes}+1)*(${lcycle}-${fcycle}+1); \ 152 sa_size_t ncols=ich${fcycle}0.NCols(); \ 153 TMatrix<r_4>ch0evol(nrows,ncols); TMatrix<r_4>ch1evol(nrows,ncols); \ 154 KeepObj(ch0evol); KeepObj(ch1evol); 155 156 for ic ${fcycle}:${maxcycle} 157 c++exec sa_size_t curRow=(${ic}-${fcycle})*(${nframes}+1); sa_size_t nextm1CurRow=curRow+(${nframes}+1)-1; \ 158 cout << "ic, curRow, next : "<<${ic}<<" "<<curRow<<" "<<nextm1CurRow<<endl; \ 159 ch0evol(Range(curRow,nextm1CurRow),Range::all())=ich${ic}0(Range(0,${nframes}),Range::all()); \ 160 ch1evol(Range(curRow,nextm1CurRow),Range::all())=ich${ic}1(Range(0,${nframes}),Range::all()); 161 end 162 163 echo "Coucou......" 164 165 #intensite en ft du temps avec image cumulee 166 del i0evol 167 del i1evol 168 #1405,1415 MHz => range [5080,5405] 169 c++exec TMatrix<r_4>tmp0(ch0evol(Range::all(),Range(5080,5405)),false); \ 170 TVector<r_4>i0evol(ch0evol.NRows()); \ 171 for(sa_size_t i=0;i<tmp0.NRows();++i){ \ 172 double mean,sigma; \ 173 MeanSigma(tmp0(Range(i),Range::all()),mean,sigma); \ 174 i0evol(i) = mean; \ 175 } \ 176 KeepObj(i0evol); \ 177 TMatrix<r_4>tmp1(ch1evol(Range::all(),Range(5080,5405)),false); \ 178 TVector<r_4>i1evol(ch1evol.NRows()); \ 179 for(sa_size_t i=0;i<tmp1.NRows();++i){ \ 180 double mean,sigma; \ 181 MeanSigma(tmp1(Range(i),Range::all()),mean,sigma); \ 182 i1evol(i) = mean; \ 183 } \ 184 KeepObj(i1evol); 185 186 187 188 newwin 189 setaxesatt 'grid' 190 graphicatt '' 191 plot2d i0evol n val 1 'blue cpts notit nsta' 192 plot2d i1evol n val 1 'same red cpts notit nsta' 193 settitle "${srcMaj} Drift Scan [1405,1415]MHz all cycles Ch 0 (blue) Ch 1 (red)" 194 setaxelabels 't(a.u)' 'I (a.u)' 195 196 endscript 197 ################### 198 #for ic ${fcycle}:${maxcycle} 199 #echo "process $ic" 200 #image $ic 201 #end 202
Note: See TracChangeset
for help on using the changeset viewer.