Changeset 591 for BAORadio/AmasNancay/trunk/etude_mergeAna.pic

Timestamp:

Oct 26, 2011, 11:36:56 AM (13 years ago)

Author:

campagne

Message:

continue improve analysis (jec)

File:

• BAORadio/AmasNancay/trunk/etude_mergeAna.pic (modified) (7 diffs)

BAORadio/AmasNancay/trunk/etude_mergeAna.pic

@@ -18 +18 @@
 
 echo $nbins1420
-
 
 
@@ -26 +25 @@
 clearscript IntAt1420
 clearscript sigmaRaw1420Side
+clearscript sigmaRaw1400a1420
 ############################################################################
 defscript evolAtCalibFreq
@@ -150 +150 @@
 
 
-
-
-#col 13 = onoffRaw01420side
-#col 14 = onoffRaw11420side
+#col 13 = onoffRaw01420side   Ch 0
+#col 14 = onoffRaw11420side   Ch 1
 
 ntcol2var onoffevol 13 linCh0All
@@ -161 +159 @@
 line2vec vecCh1All $linCh1All
 
-
+#packing values only valid for 500 cycles...
 set packVal "1 10 25 50 100"
 set sigmaCh0 ""
@@ -188 +186 @@
   mean1 = ${outvec.sum}/${outvec.size}
   sigma1 = sqrt(${outvec.sumsq}/${outvec.size}-${mean1}*${mean1})
-  errsig1 = ${sigma0}/sqrt(2*${outvec.size})    
+  errsig1 = ${sigma1}/sqrt(2*${outvec.size})    
   set sigmaCh1 "${sigmaCh1} ${sigma1}"  
   set errsigCh1 "${errsigCh1} ${errsig1}"
@@ -222 +220 @@
 
 endscript
+##################################################
+defscript sigmaRaw1400a1420
+
+c++compile rebining
+c++link rebining.so dorebin
+
+
+#col 19 = onoffRaw0f14001420   Ch 0
+#col 20 = onoffRaw1f14001420  Ch 1
+
+ntcol2var onoffevol 19 linCh0All
+line2vec vecCh0All $linCh0All
+
+ntcol2var onoffevol 20 linCh1All
+line2vec vecCh1All $linCh1All
+
+#packing values only valid for 500 cycles...
+set packVal "1 10 25 50 100"
+set sigmaCh0 ""
+set errsigCh0 ""
+set sigmaCh1 ""
+set errsigCh1 ""
+
+
+foreach ipack ( $packVal )
+
+#Packing per ipack
+  del invec 
+  del outvec
+  cp vecCh0All invec
+  call dorebin $ipack
+  mean0 = ${outvec.sum}/${outvec.size}
+  sigma0 = sqrt(${outvec.sumsq}/${outvec.size}-${mean0}*${mean0})
+  errsig0 = ${sigma0}/sqrt(2*${outvec.size})    
+  set sigmaCh0 "${sigmaCh0} ${sigma0}"  
+  set errsigCh0 "${errsigCh0} ${errsig0}"       
+#
+  del invec
+  del outvec
+  cp vecCh1All invec
+  call dorebin $ipack
+  mean1 = ${outvec.sum}/${outvec.size}
+  sigma1 = sqrt(${outvec.sumsq}/${outvec.size}-${mean1}*${mean1})
+  errsig1 = ${sigma1}/sqrt(2*${outvec.size})    
+  set sigmaCh1 "${sigmaCh1} ${sigma1}"  
+  set errsigCh1 "${errsigCh1} ${errsig1}"
+end
+
+
+newnt nt0 x y ey
+newnt nt1 x y ey
+
+set x0 ( $packVal )
+set y0 ( $sigmaCh0 )
+set ey0 ( $errsigCh0 )
+
+set x1 ( $packVal )
+set y1 ( $sigmaCh1 )
+set ey1 ( $errsigCh1 )
+        
+for i 0:$#x0
+ line2nt nt0 $x0[i] $y0[i] $ey0[i]
+ line2nt nt1 $x1[i] $y1[i] $ey1[i]
+end
+
+
+newwin 1 1
+plot2de nt0 x y 0 ey 1 "blue marker=fcircle,9 notit nsta"     
+plot2de nt1 x y 0 ey 1 "same red  marker=fcircle,9 notit nsta"
+n = ${#x0}-1
+func $y0[0]/sqrt(x) $x0[0] $x0[n]  100 "same"
+settitle "Sigma [1400,1420]MHz ${source} Ch 0 (blue) Ch 1 (red)" ' ' $defatt
+setaxelabels "num of cycles" "Sigma (a.u)" $axedefatt
+
+
+
+endscript
 ###########################################
+
 defscript evolAt1420SideFreq
 
@@ -361 +437 @@
 # Calibration coeff. per cycles
 #
+
+graphicatt "xylimits=1250,1500,-0.01,0.01"
+del sonoffovoff0 sonoffovoff1
+objaoper meanOvOffNoCalib row 0 sonoffovoff0
+objaoper meanOvOffNoCalib row 1 sonoffovoff1
+newwin 1 1
+plot2d sonoffovoff0 (n/8192)*250+1250 val n>0 "blue cpts notit nsta"
+plot2d sonoffovoff1 (n/8192)*250+1250 val n>0 "red same cpts notit nsta"
+settitle "Raw (ON-OFF)/OFF  ${source} Ch 0 (blue) Ch 1 (red) ${ncycles}cycles" ' ' $defatt
+setaxelabels "Freq. (MHz)" "I (a.u)" $axedefatt