Changeset 588 for BAORadio/AmasNancay/trunk/etude_mergeAna.pic
- Timestamp:
- Oct 24, 2011, 5:10:52 PM (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
BAORadio/AmasNancay/trunk/etude_mergeAna.pic
r582 r588 27 27 clearscript sigmaRaw1420Side 28 28 ############################################################################ 29 30 29 defscript evolAtCalibFreq 31 30 … … 147 146 defscript sigmaRaw1420Side 148 147 149 150 newh1d hSigRaw0All -0.01 0.01 25 151 newh1d hSigRaw1All -0.01 0.01 25 152 153 n/proj hSigRaw0All onoffevol.onoffRaw01420side 154 n/proj hSigRaw1All onoffevol.onoffRaw11420side 155 156 echo "Resultats sur ${ncycles} cycles (Ch0, Ch1): mean (${hSigRaw0All.mean}, ${hSigRaw1All.mean}) sigma (${hSigRaw0All.sigma}, ${hSigRaw1All.sigma})" 157 158 newwin 1 1 159 disp hSigRaw0All "blue cpts notit" 160 disp hSigRaw1All "red same cpts notit" 161 settitle "Mean and Sigma for 1 cycle" 162 163 #Ts les 10cycles 164 n10cycles = floor(${ncycles}/10) 165 x10min = ${hSigRaw0All.xmin} 166 x10max = ${hSigRaw0All.xmax} 167 n10bin = ${hSigRaw0All.nbin} 168 148 c++compile rebining 149 c++link rebining.so dorebin 150 151 152 153 154 #col 13 = onoffRaw01420side 155 #col 14 = onoffRaw11420side 156 157 ntcol2var onoffevol 13 linCh0All 158 line2vec vecCh0All $linCh0All 159 160 ntcol2var onoffevol 14 linCh1All 161 line2vec vecCh1All $linCh1All 162 163 164 set packVal "1 10 25 50 100" 165 set sigmaCh0 "" 166 set errsigCh0 "" 167 set sigmaCh1 "" 168 set errsigCh1 "" 169 170 171 foreach ipack ( $packVal ) 172 173 #Packing per ipack 174 del invec 175 del outvec 176 cp vecCh0All invec 177 call dorebin $ipack 178 mean0 = ${outvec.sum}/${outvec.size} 179 sigma0 = sqrt(${outvec.sumsq}/${outvec.size}-${mean0}*${mean0}) 180 errsig0 = ${sigma0}/sqrt(2*${outvec.size}) 181 set sigmaCh0 "${sigmaCh0} ${sigma0}" 182 set errsigCh0 "${errsigCh0} ${errsig0}" 183 # 184 del invec 185 del outvec 186 cp vecCh1All invec 187 call dorebin $ipack 188 mean1 = ${outvec.sum}/${outvec.size} 189 sigma1 = sqrt(${outvec.sumsq}/${outvec.size}-${mean1}*${mean1}) 190 errsig1 = ${sigma0}/sqrt(2*${outvec.size}) 191 set sigmaCh1 "${sigmaCh1} ${sigma1}" 192 set errsigCh1 "${errsigCh1} ${errsig1}" 193 end 194 195 196 newnt nt0 x y ey 197 newnt nt1 x y ey 198 199 set x0 ( $packVal ) 200 set y0 ( $sigmaCh0 ) 201 set ey0 ( $errsigCh0 ) 202 203 set x1 ( $packVal ) 204 set y1 ( $sigmaCh1 ) 205 set ey1 ( $errsigCh1 ) 169 206 170 for i 0:${n10cycles} 171 ic = ${i}*10+1 172 jc = ${ic}+9 173 newh1d hSigRaw010c${i} ${x10min} ${x10max} ${n10bin} 174 newh1d hSigRaw110c${i} ${x10min} ${x10max} ${n10bin} 175 176 n/proj hSigRaw010c${i} onoffevol.onoffRaw01420side (cycle>=${ic})&&(cycle<=${jc}) 177 cp hSigRaw010c${i} h0 178 mSig010c[i] = ${h0.mean} 179 180 n/proj hSigRaw110c${i} onoffevol.onoffRaw11420side (cycle>=${ic})&&(cycle<=${jc}) 181 cp hSigRaw110c${i} h1 182 mSig110c[i] = ${h1.mean} 207 for i 0:$#x0 208 line2nt nt0 $x0[i] $y0[i] $ey0[i] 209 line2nt nt1 $x1[i] $y1[i] $ey1[i] 183 210 end 184 211 185 212 186 line2vec hSigRaw010c $mSig010c 187 line2vec hSigRaw110c $mSig110c 188 189 190 191 192 newwin 1 1 193 disp hSigRaw010c "blue cpts notit" 194 disp hSigRaw110c "red same cpts notit" 195 settitle "Mean and Sigma for 10 cycles" 213 newwin 1 1 214 plot2de nt0 x y 0 ey 1 "blue marker=fcircle,9 notit nsta" 215 plot2de nt1 x y 0 ey 1 "same red marker=fcircle,9 notit nsta" 216 n = ${#x0}-1 217 func $y0[0]/sqrt(x) $x0[0] $x0[n] 100 "same" 218 settitle "Sigma [1418,1419]U[1422,1423]MHz ${source} Ch 0 (blue) Ch 1 (red)" ' ' $defatt 219 setaxelabels "num of cycles" "Sigma (a.u)" $axedefatt 220 196 221 197 222 … … 389 414 390 415 graphicatt "" 391 setaxesatt " "416 setaxesatt "$axedefatt" 392 417 393 418
Note: See TracChangeset
for help on using the changeset viewer.